PWGLF/STRANGENESS/LambdaK0PbPb/MainStreamAnalysis/AnalysisStrange.C

   1 const char *anatype = "ESD";
   2
   3 void AnalysisStrange()
   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 base root libraries
  20    gSystem->Load("libTree");
  21    gSystem->Load("libGeom");
  22    gSystem->Load("libVMC");
  23    gSystem->Load("libPhysics");
  24
  25    gSystem->Load("libMinuit");
  26
  27 // Load analysis framework libraries
  28    gSystem->Load("libSTEERBase");
  29    gSystem->Load("libESD");
  30    gSystem->Load("libAOD");
  31    gSystem->Load("libANALYSIS");
  32    gSystem->Load("libOADB");
  33    gSystem->Load("libANALYSISalice");
  34    gSystem->Load("libCORRFW");
  35
  36 // include path
  37    TString intPath = gInterpreter->GetIncludePath();
  38    TObjArray *listpaths = intPath.Tokenize(" ");
  39    TIter nextpath(listpaths);
  40    TObjString *pname;
  41    while ((pname=(TObjString*)nextpath())) {
  42       TString current = pname->GetName();
  43       if (current.Contains("AliRoot") || current.Contains("ALICE_ROOT")) continue;
  44       gSystem->AddIncludePath(current);
  45    }
  46    if (listpaths) delete listpaths;
  47    gROOT->ProcessLine(".include $ALICE_ROOT/include");
  48    printf("Include path: %s\n", gSystem->GetIncludePath());
  49
  50 // Add aditional AliRoot libraries
  51
  52 // analysis source to be compiled at runtime (if any)
  53    gROOT->ProcessLine(".L AliAnalysisCentralitySelector.cxx+g");
  54    gROOT->ProcessLine(".L AliAnalysisTaskPerformanceStrange.cxx+g");
  55
  56 // connect to AliEn and make the chain
  57    if (!TGrid::Connect("alien://")) return;
  58 // read the analysis manager from file
  59    AliAnalysisManager *mgr = AliAnalysisAlien::LoadAnalysisManager("analysis.root");
  60    if (!mgr) return;
  61    mgr->PrintStatus();
  62    AliLog::SetGlobalLogLevel(AliLog::kError);
  63    TChain *chain = CreateChain("wn.xml", anatype);
  64
  65    mgr->StartAnalysis("localfile", chain);
  66    timer.Stop();
  67    timer.Print();
  68 }
  69
  70 //________________________________________________________________________________
  71 TChain* CreateChain(const char *xmlfile, const char *type="ESD")
  72 {
  73 // Create a chain using url's from xml file
  74    TString filename;
  75    Int_t run = 0;
  76    TString treename = type;
  77    treename.ToLower();
  78    treename += "Tree";
  79    printf("***************************************\n");
  80    printf("    Getting chain of trees %s\n", treename.Data());
  81    printf("***************************************\n");
  82    TAlienCollection *coll = TAlienCollection::Open(xmlfile);
  83    if (!coll) {
  84       ::Error("CreateChain", "Cannot create an AliEn collection from %s", xmlfile);
  85       return NULL;
  86    }
  87    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
  88    TChain *chain = new TChain(treename);
  89    coll->Reset();
  90    while (coll->Next()) {
  91       filename = coll->GetTURL();
  92       if (mgr) {
  93          Int_t nrun = AliAnalysisManager::GetRunFromAlienPath(filename);
  94          if (nrun && nrun != run) {
  95             printf("### Run number detected from chain: %d\n", nrun);
  96             mgr->SetRunFromPath(nrun);
  97             run = nrun;
  98          }
  99       }
 100       chain->Add(filename);
 101    }
 102    if (!chain->GetNtrees()) {
 103       ::Error("CreateChain", "No tree found from collection %s", xmlfile);
 104       return NULL;
 105    }
 106    return chain;
 107 }
 108