ANALYSIS/CalibMacros/runCalibTrain.C

   1 /*
   2   Template of calibration/filtering  macro using ESD
   3
   4   Example:
   5   .L $ALICE_ROOT/ANALYSIS/macros/runCalibTrain.C
   6   runCalibTrain(105160);
   7
   8 */
   9
  10 void runCalibTrain(TString runNumberString, const char *inFileName = "AliESDs.root")
  11 {
  12   gROOT->Macro("LoadLibraries.C");
  13   gROOT->LoadMacro("ConfigCalibTrain.C");
  14   gROOT->LoadMacro("AddTaskFilterFriend.C");
  15   gROOT->LoadMacro("AddTaskFilterFriendSecond.C");
  16   gROOT->LoadMacro("AddTaskAddObject.C");
  17
  18   // detector tasks
  19
  20   gROOT->LoadMacro("AddTaskTPCCalib.C");
  21
  22   AliLog::SetClassDebugLevel("AliESDEvent",19);
  23   TChain *chain = new TChain("esdTree");
  24
  25   // Steering input chain
  26
  27   chain->Add(inFileName);
  28   Int_t runNumber = runNumberString.Atoi();
  29   printf("runNumber from runCalibTrain = %d\n",runNumber);
  30   ConfigCalibTrain(runNumber, "raw://");
  31
  32   AliAnalysisManager *mgr  = new AliAnalysisManager("ESD to ESD", "Analysis Manager");
  33   // mgr->SetDebugLevel(3);
  34
  35   // Input
  36
  37   AliESDInputHandler* inpHandler = new AliESDInputHandler();
  38   inpHandler->SetActiveBranches("ESDfriend*");
  39   mgr->SetInputEventHandler  (inpHandler);
  40
  41   // Output
  42
  43   AliESDHandler* esdHandler   = new AliESDHandler();
  44   mgr->SetOutputEventHandler(esdHandler);
  45   esdHandler->SetOutputFileName("AliESDfriends_v1.root");
  46   // Steering Tasks
  47
  48   AliAnalysisTaskFilterFriend* filter = AddTaskFilterFriend();
  49   AliAnalysisTaskFilterFriendSecond* filter2 = AddTaskFilterFriendSecond();
  50   AliAnalysisTaskAddObject* add = AddTaskAddObject();
  51
  52   // Detector Tasks
  53
  54   AliAnalysisTask* tTPC = AddTaskTPCCalib(runNumber);
  55
  56   // Run the analysis
  57
  58   if (!mgr->InitAnalysis()) {
  59     printf("Analysis cannot be started, returning\n");
  60     return;
  61   }
  62
  63   mgr->PrintStatus();
  64   mgr->StartAnalysis("local", chain);
  65
  66   return;
  67 }