2 Template of calibration/filtering macro using ESD
4 .L $ALICE_ROOT/ANALYSIS/macros/runCalibTrain.C
9 #if !defined(__CINT__) || defined(__MAKECINT__)
12 #include "AliAnalysisManager.h"
13 #include "AliESDInputHandler.h"
14 #include "AliAnalysisTaskESDfilter.h"
15 #include "AliAnalysisDataContainer.h"
18 void runCalibTrain(TString runNumberString, const char *inFileName = "AliESDs.root",
19 const char *outFileName = "AliESDs_v1.root")
22 gSystem->Load("libTree");
23 gSystem->Load("libGeom");
24 gSystem->Load("libPhysics");
25 gSystem->Load("libVMC");
26 gSystem->Load("libSTEERBase");
27 gSystem->Load("libESD");
28 gSystem->Load("libAOD");
30 gSystem->Load("libANALYSIS");
31 gSystem->Load("libANALYSISalice");
32 gSystem->Load("libANALYSIScalib");
33 gSystem->Load("libCORRFW");
34 gSystem->Load("libPWG3muon");
35 gROOT->LoadMacro("ConfigCalibTrain.C");
36 gROOT->LoadMacro("AddTaskCopyESD.C");
37 gROOT->LoadMacro("AddTaskFilterFriend.C");
38 gROOT->LoadMacro("AddTaskFilterFriendSecond.C");
39 gROOT->LoadMacro("AddTaskAddObject.C");
43 gROOT->LoadMacro("AddTaskTPCCalib.C");
45 AliLog::SetClassDebugLevel("AliESDEvent",19);
46 TChain *chain = new TChain("esdTree");
48 // Steering input chain
50 chain->Add(inFileName);
51 Int_t runNumber = runNumberString.Atoi();
52 printf("runNumber from runCalibTrain = %d\n",runNumber);
53 ConfigCalibTrain(runNumber, "raw://");
55 AliAnalysisManager *mgr = new AliAnalysisManager("ESD to ESD", "Analysis Manager");
56 // mgr->SetDebugLevel(3);
60 AliESDInputHandler* inpHandler = new AliESDInputHandler();
61 inpHandler->SetActiveBranches("ESDfriend*");
62 mgr->SetInputEventHandler (inpHandler);
66 AliESDHandler* esdHandler = new AliESDHandler();
67 esdHandler->SetOutputFileName(outFileName);
68 mgr->SetOutputEventHandler(esdHandler);
72 AliAnalysisTaskCopyESD *copy = AddTaskCopyESD();
73 AliAnalysisTaskFilterFriend* filter = AddTaskFilterFriend();
74 AliAnalysisTaskFilterFriendSecond* filter2 = AddTaskFilterFriendSecond();
75 AliAnalysisTaskAddObject* add = AddTaskAddObject();
79 AliAnalysisTask* tTPC = AddTaskTPCCalib(runNumber);
83 if (!mgr->InitAnalysis()) {
84 printf("Analysis cannot be started, returning\n");
89 mgr->StartAnalysis("local", chain);