PWGLF/SPECTRA/IdentifiedHighPt/grid/macro_HighPtDeDx_lhc10b_Data_ESDs.C

   1 const char *anatype = "ESD";
   2
   3 void macro_HighPtDeDx_lhc10b_Data_ESDs()
   4 {
   5 // Analysis using ESD data
   6 // Automatically generated analysis steering macro executed in grid subjobs
   7
   8    TStopwatch timer;
   9    timer.Start();
  10
  11 // Set temporary merging directory to current one
  12    gSystem->Setenv("TMPDIR", gSystem->pwd());
  13
  14 // Set temporary compilation directory to current one
  15    gSystem->SetBuildDir(gSystem->pwd(), kTRUE);
  16
  17 // Reset existing include path and add current directory first in the search
  18    gSystem->SetIncludePath("-I.");
  19 // Load analysis framework libraries
  20    gSystem->Load("libANALYSIS");
  21    gSystem->Load("libOADB");
  22    gSystem->Load("libANALYSISalice");
  23    gSystem->Load("libCORRFW");
  24
  25 // include path
  26    TString intPath = gInterpreter->GetIncludePath();
  27    TObjArray *listpaths = intPath.Tokenize(" ");
  28    TIter nextpath(listpaths);
  29    TObjString *pname;
  30    while ((pname=(TObjString*)nextpath())) {
  31       TString current = pname->GetName();
  32       if (current.Contains("AliRoot") || current.Contains("ALICE_ROOT")) continue;
  33       gSystem->AddIncludePath(current);
  34    }
  35    if (listpaths) delete listpaths;
  36    gSystem->AddIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include ");
  37    gROOT->ProcessLine(".include $ALICE_ROOT/include");
  38    printf("Include path: %s\n", gSystem->GetIncludePath());
  39
  40 // Add aditional AliRoot libraries
  41
  42 // analysis source to be compiled at runtime (if any)
  43    gROOT->ProcessLine(".L DebugClasses.C+g");
  44    gROOT->ProcessLine(".L AliAnalysisTaskHighPtDeDx.cxx+g");
  45    gROOT->ProcessLine(".L AliAnalysisTaskHighPtDeDxV0.cxx+g");
  46
  47 // fast xrootd reading enabled
  48    printf("!!! You requested FastRead option. Using xrootd flags to reduce timeouts. Note that this may skip some files that could be accessed !!!");
  49    gEnv->SetValue("XNet.ConnectTimeout",50);
  50    gEnv->SetValue("XNet.RequestTimeout",50);
  51    gEnv->SetValue("XNet.MaxRedirectCount",2);
  52    gEnv->SetValue("XNet.ReconnectTimeout",50);
  53    gEnv->SetValue("XNet.FirstConnectMaxCnt",1);
  54
  55 // connect to AliEn and make the chain
  56    if (!TGrid::Connect("alien://")) return;
  57 // read the analysis manager from file
  58    AliAnalysisManager *mgr = AliAnalysisAlien::LoadAnalysisManager("HighPtDeDx_lhc10b_Data_ESDs.root");
  59    if (!mgr) return;
  60    mgr->PrintStatus();
  61    AliLog::SetGlobalLogLevel(AliLog::kError);
  62    TChain *chain = CreateChain("wn.xml", anatype);
  63
  64    mgr->StartAnalysis("localfile", chain);
  65    timer.Stop();
  66    timer.Print();
  67 }
  68
  69 //________________________________________________________________________________
  70 TChain* CreateChain(const char *xmlfile, const char *type="ESD")
  71 {
  72 // Create a chain using url's from xml file
  73    TString filename;
  74    Int_t run = 0;
  75    TString treename = type;
  76    treename.ToLower();
  77    treename += "Tree";
  78    printf("***************************************\n");
  79    printf("    Getting chain of trees %s\n", treename.Data());
  80    printf("***************************************\n");
  81    TAlienCollection *coll = TAlienCollection::Open(xmlfile);
  82    if (!coll) {
  83       ::Error("CreateChain", "Cannot create an AliEn collection from %s", xmlfile);
  84       return NULL;
  85    }
  86    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  87    TChain *chain = new TChain(treename);
  88    coll->Reset();
  89    while (coll->Next()) {
  90       filename = coll->GetTURL();
  91       if (mgr) {
  92          Int_t nrun = AliAnalysisManager::GetRunFromAlienPath(filename);
  93          if (nrun && nrun != run) {
  94             printf("### Run number detected from chain: %d\n", nrun);
  95             mgr->SetRunFromPath(nrun);
  96             run = nrun;
  97          }
  98       }
  99       chain->Add(filename);
 100    }
 101    if (!chain->GetNtrees()) {
 102       ::Error("CreateChain", "No tree found from collection %s", xmlfile);
 103       return NULL;
 104    }
 105    return chain;
 106 }
 107