Float_t etaMax[numberOfSyst] = {0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,1.0,0.6,0.4}; // eta cuts\r
\r
Bool_t kUsePID = kFALSE;\r
+Bool_t kUseHBTCut = kFALSE;\r
+Bool_t kUseConversionCut = kFALSE;\r
\r
//______________________________________________________________________________\r
void runBalanceFunctionPsi(\r
gSystem->Load("libAOD.so");\r
gSystem->Load("libANALYSIS.so");\r
gSystem->Load("libANALYSISalice.so");\r
+ gSystem->Load("libEventMixing.so");\r
+ gSystem->Load("libCORRFW.so");\r
+ gSystem->Load("libPWGTools.so");\r
gSystem->Load("libPWGCFebye.so");\r
\r
// additional\r
TChain* chain = 0x0;\r
if((!bAOD)&&(!bMCtruth)) {\r
chain = new TChain("esdTree");\r
- for(Int_t i = 0; i < 8; i++) {\r
+ for(Int_t i = 0; i < 4; i++) {\r
filename = "/data/alice2/pchrist/HeavyIons/Data/2011/Set";\r
filename += i; filename += "/AliESDs.root";\r
chain->Add(filename.Data());\r
}\r
else if((bAOD)&&(!bMCtruth)) {\r
chain = new TChain("aodTree");\r
- for(Int_t i = 0; i < 8; i++) {\r
- filename = "/data/alice2/pchrist/HeavyIons/Data/2011/Set";\r
+ for(Int_t i = 0; i < 4; i++) {\r
+ filename = "/project/alice/pchrist/Data/2011/Set";\r
filename += i; filename += "/AliAOD.root";\r
chain->Add(filename.Data());\r
}\r
\r
//Add the BF task (all centralities)\r
gROOT->LoadMacro("AddTaskBalancePsiCentralityTrain.C"); \r
- AliAnalysisTaskBFPsi *task = AddTaskBalancePsiCentralityTrain(0,100,1,"V0M",vZ[0],DCAxy[0],DCAz[0],ptMin[0],ptMax[0],etaMin[0],etaMax[0],-1,-1,kUsePID);\r
+ AliAnalysisTaskBFPsi *task = AddTaskBalancePsiCentralityTrain(0,100,1,"V0M",vZ[0],DCAxy[0],DCAz[0],ptMin[0],ptMax[0],etaMin[0],etaMax[0],-1,-1,kUseHBTCut,kUseConversionCut,kUsePID);\r
\r
// enable debug printouts\r
//mgr->SetDebugLevel(2);\r