PWG3/vertexingHF/RunAnalysisAODVertexingHF.C

   1 void RunAnalysisAODVertexingHF()\r
   2 {\r
   3   //\r
   4   // Test macro for AliAnalysisTaskSE's for heavy-flavour candidates\r
   5   // It has the structure of a Analysis Train:\r
   6   // - in this macro, change things related to running mode\r
   7   //   and input preparation \r
   8   // - add your task using a AddTaskXXX macro \r
   9   //\r
  10   // A.Dainese, andrea.dainese@lnl.infn.it\r
  11   // "grid" mode added by R.Bala, bala@to.infn.it\r
  12   //\r
  13 \r
  14   TString analysisMode = "local"; // "local" or "grid"\r
  15   TString inputMode    = "list"; // "list" or "xml"\r
  16 \r
  17 \r
  18   Bool_t useParFiles=kFALSE;\r
  19 \r
  20 \r
  21   gROOT->LoadMacro("$ALICE_ROOT/PWG3/vertexingHF/LoadLibraries.C");\r
  22   gROOT->LoadMacro("$ALICE_ROOT/PWG3/vertexingHF/MakeAODInputChain.C");\r
  23   LoadLibraries(useParFiles);\r
  24 \r
  25   if(analysisMode=="grid") TGrid::Connect("alien:",0,0,"t");\r
  26 \r
  27 \r
  28 \r
  29 \r
  30   //-------------------------------------------------------------------\r
  31   // Prepare input chain\r
  32   TChain *chainAOD = 0;\r
  33   \r
  34   if(inputMode=="list") {\r
  35     // Local files\r
  36     chainAOD = MakeAODInputChain();// with this it reads ./AliAOD.root and ./AliAOD.VertexingHF.root\r
  37     //chainAOD = MakeAODInputChain("alien:///alice/cern.ch/user/r/rbala/analysis/out_lhcw/290001/",2,2);\r
  38   } else if(inputMode=="xml") {\r
  39     // xml\r
  40     chainAOD = MakeAODInputChain("collection_aod.xml","collection_aodHF.xml");\r
  41   }\r
  42 \r
  43   // Create the analysis manager\r
  44   AliAnalysisManager *mgr  = new AliAnalysisManager("My Manager","My Manager");\r
  45   mgr->SetDebugLevel(10);\r
  46   \r
  47 \r
  48   // Input\r
  49   AliAODInputHandler *inputHandler = new AliAODInputHandler();\r
  50   mgr->SetInputEventHandler(inputHandler);\r
  51   //-------------------------------------------------------------------\r
  52 \r
  53   \r
  54   //-------------------------------------------------------------------\r
  55   // Analysis tasks    \r
  56   //\r
  57   gROOT->LoadMacro("$ALICE_ROOT/PWG3/vertexingHF/AddTaskCompareHF.C");\r
  58   AliAnalysisTaskSECompareHF *cmpTask = AddTaskCompareHF();\r
  59   //gROOT->LoadMacro("$ALICE_ROOT/PWG3/vertexingHF/AddTaskSelectHF.C");\r
  60   //AliAnalysisTaskSESelectHF *seleTask = AddTaskSelectHF();\r
  61   //gROOT->LoadMacro("$ALICE_ROOT/PWG3/vertexingHF/AddTaskBkgLikeSign.C");\r
  62   //AliAnalysisTaskSEBkgLikeSignJPSI *lsTask = AddTaskBkgLikeSign();\r
  63   //gROOT->LoadMacro("$ALICE_ROOT/PWG3/vertexingHF/AddTaskBtoJPSItoEle.C");\r
  64   //AliAnalysisTaskSEBtoJPSItoEle *jpsiTask = AddTaskBtoJPSItoEle();\r
  65   gROOT->LoadMacro("$ALICE_ROOT/PWG3/vertexingHF/AddTaskCF.C");\r
  66   AliCFHeavyFlavourTask *cfTask = AddTaskCF();\r
  67   //gROOT->LoadMacro("$ALICE_ROOT/PWG3/vertexingHF/AddTaskCFMultiVar.C");\r
  68   //AliCFHeavyFlavourTaskMultiVar *cfmvTask = AddTaskCFMultiVar();\r
  69 \r
  70 \r
  71   //-------------------------------------------------------------------\r
  72 \r
  73   //\r
  74   // Run the analysis\r
  75   //    \r
  76   printf("CHAIN HAS %d ENTRIES\n",(Int_t)chainAOD->GetEntries());\r
  77   if(mgr->InitAnalysis()) {\r
  78     mgr->PrintStatus();\r
  79     mgr->StartAnalysis("local",chainAOD);\r
  80   }\r
  81 \r
  82   return;\r
  83 }\r