TENDER becomes Tender, removing .so
[u/mrichter/AliRoot.git] / PWGGA / GammaConv / macros / AddTask_GammaConvDalitzV1_pPb.C
1 void AddTask_GammaConvDalitzV1_pPb(    Int_t trainConfig = 1,
2                                        Bool_t isMC       = kFALSE, //run MC 
3                                        Bool_t enableQAMesonTask = kTRUE, //enable QA in AliAnalysisTaskGammaConvDalitzV1
4                                        Bool_t enableDoMesonChic = kFALSE, // enable additional Chic analysis
5                                        TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
6                                        Bool_t doWeighting = kFALSE,  //enable Weighting
7                                        TString generatorName = "DPMJET",                                
8                                        TString cutnumberAODBranch = "0000000060084001001500000"
9                                   ) {
10
11
12    
13    cout<<"*********Parameters*******"<<endl;
14    cout<<"trainConfig: "<<trainConfig<<endl;
15    cout<<"isMC: "<<isMC<<endl;
16    cout<<"enableQAMesonTask: "<<enableQAMesonTask<<endl;
17    cout<<"enableDoMesonChic: "<<enableDoMesonChic<<endl;
18    cout<<"fileNameInputForWeighting: "<<fileNameInputForWeighting.Data()<<endl;
19    cout<<"doWeighting: "<<doWeighting<<endl;
20    cout<<"generatorName: "<<generatorName.Data()<<endl;
21    cout<<"cutnumberAODBranch: "<<cutnumberAODBranch.Data()<<endl;
22
23   // ================= Load Librariers =================================
24    gSystem->Load("libCore");  
25    gSystem->Load("libTree");
26    gSystem->Load("libGeom");
27    gSystem->Load("libVMC");
28    gSystem->Load("libPhysics");
29    gSystem->Load("libMinuit");
30    gSystem->Load("libSTEERBase");
31    gSystem->Load("libESD");
32    gSystem->Load("libAOD");
33    gSystem->Load("libANALYSIS");
34    gSystem->Load("libANALYSISalice");  
35    gSystem->Load("libPWGGAGammaConv");
36    gSystem->Load("libCDB");
37    gSystem->Load("libSTEER");
38    gSystem->Load("libSTEERBase");
39    gSystem->Load("libTender");
40    gSystem->Load("libTenderSupplies");
41
42
43    cout<<"Entro 0"<<endl;
44    
45    Int_t isHeavyIon = 2;
46
47    // ================== GetAnalysisManager ===============================
48    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
49    if (!mgr) {
50       Error(Form("AddTask_GammaConvDalitzV1_pPb_%i",trainConfig), "No analysis manager found.");
51       return ;
52    }
53
54    // ================== GetInputEventHandler =============================
55    AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
56    
57    //========= Add PID Reponse to ANALYSIS manager ====
58    if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
59       gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
60       AddTaskPIDResponse(isMC);
61    }
62    
63    //=========  Set Cutnumber for V0Reader ================================
64    TString cutnumberEvent = "8000000";
65    
66    TString cutnumberPhoton="060000084001001500000000";   //Online  V0 finder
67    
68    TString ElecCuts      = "90005400000002000000";            //Electron Cuts
69    
70    Bool_t doEtaShift = kFALSE;
71   
72
73
74    AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
75    
76    //========= Add V0 Reader to  ANALYSIS manager if not yet existent =====
77    if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
78                 AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
79                 
80                 fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
81                 fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
82                 fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
83                 
84                 if (!mgr) {
85                         Error("AddTask_V0ReaderV1", "No analysis manager found.");
86                         return;
87                 }
88         
89                 AliConvEventCuts *fEventCuts=NULL;
90                 if(cutnumberEvent!=""){
91                         fEventCuts= new AliConvEventCuts(cutnumberEvent.Data(),cutnumberEvent.Data());
92                         fEventCuts->SetPreSelectionCutFlag(kTRUE);
93                         if(fEventCuts->InitializeCutsFromCutString(cutnumberEvent.Data())){
94                                 fEventCuts->DoEtaShift(doEtaShift);
95                                 fV0ReaderV1->SetEventCuts(fEventCuts);
96                                 fEventCuts->SetFillCutHistograms("",kTRUE);
97                         }
98                 }
99
100                 // Set AnalysisCut Number
101                 AliConversionPhotonCuts *fCuts=NULL;
102                 if(cutnumberPhoton!=""){
103                         fCuts= new AliConversionPhotonCuts(cutnumberPhoton.Data(),cutnumberPhoton.Data());
104                         fCuts->SetPreSelectionCutFlag(kTRUE);
105                         fCuts->SetIsHeavyIon(isHeavyIon);
106                         if(fCuts->InitializeCutsFromCutString(cutnumberPhoton.Data())){
107                                 fV0ReaderV1->SetConversionCuts(fCuts);
108                                 fCuts->SetFillCutHistograms("",kTRUE);
109                         }
110                 }
111                 if(inputHandler->IsA()==AliAODInputHandler::Class()){
112                 // AOD mode
113                         cout << "AOD handler: adding " << cutnumberAODBranch.Data() << " as conversion branch" << endl;
114                         fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
115                 }
116                 fV0ReaderV1->Init();
117         
118                 AliLog::SetGlobalLogLevel(AliLog::kInfo);
119         
120                 //connect input V0Reader
121                 mgr->AddTask(fV0ReaderV1);
122                 mgr->ConnectInput(fV0ReaderV1,0,cinput);
123         
124     }
125
126    //================================================
127    //========= Add Electron Selector ================
128
129
130    if( !(AliDalitzElectronSelector*)mgr->GetTask("ElectronSelector") ){
131
132    AliDalitzElectronSelector *fElectronSelector = new AliDalitzElectronSelector("ElectronSelector");
133
134    // Set AnalysisCut Number
135
136    AliDalitzElectronCuts *fElecCuts=0;
137
138    
139
140     if( ElecCuts!=""){
141
142        fElecCuts= new AliDalitzElectronCuts(ElecCuts.Data(),ElecCuts.Data());
143
144             if(fElecCuts->InitializeCutsFromCutString(ElecCuts.Data())){
145
146                 fElectronSelector->SetDalitzElectronCuts(fElecCuts);
147
148                 fElecCuts->SetFillCutHistograms("",kTRUE);
149
150             }
151
152     }
153
154     fElectronSelector->Init();
155     mgr->AddTask(fElectronSelector);
156     
157     AliAnalysisDataContainer *cinput1  = mgr->GetCommonInputContainer();
158
159     //connect input V0Reader
160
161     mgr->ConnectInput (fElectronSelector,0,cinput1);
162
163     }
164
165
166
167     cout<<"Entro"<<endl;
168    //================================================
169    //========= Add task to the ANALYSIS manager =====
170    //================================================
171    //            find input container
172    
173   
174  
175    AliAnalysisTaskGammaConvDalitzV1 *task=NULL;
176
177    task= new AliAnalysisTaskGammaConvDalitzV1(Form("GammaConvDalitzV1_%i",trainConfig));
178
179    task->SetIsHeavyIon(isHeavyIon);
180    task->SetIsMC(isMC);
181
182
183
184    // Cut Numbers to use in Analysis
185    Int_t numberOfCuts = 4;
186
187    
188    
189    TString *eventCutArray   = new TString[numberOfCuts];
190    TString *photonCutArray  = new TString[numberOfCuts];
191    TString *ElecCutarray    = new TString[numberOfCuts];
192    TString *MesonCutarray   = new TString[numberOfCuts];
193
194    Bool_t doEtaShiftIndCuts = kFALSE;
195    TString stringShift = "";
196
197    // Shifting in pPb direction
198
199    doEtaShiftIndCuts = kFALSE;
200    stringShift = "pPb";
201
202
203    
204    
205 if( trainConfig == 1 ) {  // No eta shift |Y| < 0.8
206         
207         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011
208         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007900000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Chi2 < 15
209         eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007800000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Chi2 < 20
210         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007100000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Chi2 < 50
211         
212 }  else if( trainConfig == 2 ) {  // No eta shift |Y| < 0.8
213                                                          
214         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603002200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Qt < 0.7
215         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603003200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Qt < 0.5
216         eventCutArray[2]="8000011"; photonCutArray[2] = "002000093653007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec sec 0.3 GeV Low and 3.5 High momentum
217         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093601007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec sec 0.5 GeV Low and 5.0 High momentum
218         
219 }  else if( trainConfig == 3 ) {  // No eta shift |Y| < 0.8
220
221         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093803007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec  sec  2.0sigmas Low and  1 High momentum
222         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90435400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec primary 2.0sigmas Low and 0 High momentum
223         eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90477400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec primary 0.3 GeV Low and 3.5 High momentum
224         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475200233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec primary 0.5 GeV Low and 5.0 High momentum
225         
226 }  else if( trainConfig == 4 ) {  // No eta shift  |Y| < 0.8
227   
228         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90425400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec  primary  2.0sigmas Low and -1 High momentum 
229         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400133102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + SPD first layer
230         eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233302623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + PsiPair cut 0.52
231         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031005009000"; //standard cut Pi0 pPb 00-100 Standard cut + Alpha cut < 0.7 
232         
233 }  else if( trainConfig == 5 ) { // No eta shift  |Y| < 0.8     
234   
235         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90375400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary   electron -5,5
236         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90575400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary   electron -3,5
237         eventCutArray[2]="8000011"; photonCutArray[2] = "002000091603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -5,5
238         eventCutArray[3]="8000011"; photonCutArray[3] = "002000092603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -3,5
239         
240 } else if ( trainConfig == 6 ) { //No eta shift   |Y| < 0.8
241         
242         eventCutArray[0]="8000011"; photonCutArray[0] = "042000093603007200000000"; ElecCutarray[0] = "90475400235102623710"; MesonCutarray[0] = "01032035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Y < 0.70  and prim and sec e |eta| < 0.75 //NOTE revisar
243         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233102633710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Single prim Pt cut > 0.150
244         eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102631710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + Single prim Pt cut > 0.100
245         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102622710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + DCAxy < 1 cm
246         
247         
248 } else if ( trainConfig == 7 ) {  // No eta shift |Y| < 0.8
249
250         eventCutArray[0]="8000011"; photonCutArray[0] = "002000493603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Single sec  Pt cut > 0.075
251         eventCutArray[1]="8000011"; photonCutArray[1] = "002000193603007200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Single sec  Pt cut > 0.100
252         eventCutArray[2]="8000011"; photonCutArray[2] = "002000083603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Findable Cls sec  > 0.35
253         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400273102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Findable Cls prim > 0.60
254         
255         
256 } else if ( trainConfig == 8 ) {  //No eta shift |Y| < 0.8
257         
258         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623810"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.015 < InvMass(e+,e-) < 0.050
259         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233102623910"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.025 < InvMass(e+,e-) < 0.035
260         eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603001200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + qT < 0.1
261         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011
262         
263 } else if ( trainConfig == 9 ) {  //No eta shift |Y| < 0.8
264
265   
266         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102723710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Stardad cut +100 events background
267         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233101623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Background method V0 multiplicity
268         eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035000000"; //standard cut Pi0 PbPb 00-100 + No extra smearing
269         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400253102621710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 + Old Standard
270
271 } else if( trainConfig == 10 ) {  // No eta shift |Y| < 0.8 + AddedSignals
272         
273         eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011
274         eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007900000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Chi2 < 15
275         eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007800000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Chi2 < 20
276         eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007100000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Chi2 < 50
277         
278 }  else if( trainConfig == 11 ) {  // No eta shift |Y| < 0.8 + AddedSignals
279
280         eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603002200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Qt < 0.7
281         eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603003200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Qt < 0.5
282         eventCutArray[2]="8000012"; photonCutArray[2] = "002000093653007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec sec 0.3 GeV Low and 3.5 High momentum
283         eventCutArray[3]="8000012"; photonCutArray[3] = "002000093601007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec sec 0.5 GeV Low and 5.0 High momentum
284         
285 }  else if( trainConfig == 12 ) {  // No eta shift |Y| < 0.8 + AddedSignals
286
287         eventCutArray[0]="8000012"; photonCutArray[0] = "002000093803007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec  sec  2.0sigmas Low and  1 High momentum
288         eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90435400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec primary 2.0sigmas Low and 0 High momentum
289         eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90477400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec primary 0.3 GeV Low and 3.5 High momentum
290         eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475200233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + dEdx pion rejec primary 0.5 GeV Low and 5.0 High momentum
291         
292 }  else if( trainConfig == 13 ) {  // No eta shift  |Y| < 0.8 + AddedSignals
293   
294         eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90425400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx pion rejec  primary  2.0sigmas Low and -1 High momentum 
295         eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400133102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + SPD first layer
296         eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233302623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + PsiPair cut 0.52
297         eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031005009000"; //standard cut Pi0 pPb 00-100 Standard cut + Alpha cut < 0.7 
298         
299 }  else if( trainConfig == 14 ) { // No eta shift  |Y| < 0.8 + AddedSignals
300   
301         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90375400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary   electron -5,5
302         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90575400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx primary   electron -3,5
303         eventCutArray[2]="8000011"; photonCutArray[2] = "002000091603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -5,5
304         eventCutArray[3]="8000011"; photonCutArray[3] = "002000092603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 New Standard cut + dEdx secondary electron -3,5
305         
306 } else if ( trainConfig == 15 ) { //No eta shift   |Y| < 0.8 + AddedSignals
307         
308         eventCutArray[0]="8000012"; photonCutArray[0] = "042000093603007200000000"; ElecCutarray[0] = "90475400235102623710"; MesonCutarray[0] = "01032035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Y < 0.70  and prim and sec e |eta| < 0.75 //NOTE revisar
309         eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233102633710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Single prim Pt cut > 0.150
310         eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102631710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + Single prim Pt cut > 0.100
311         eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102622710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + DCAxy < 1 cm
312         
313         
314 } else if ( trainConfig == 16 ) {  // No eta shift |Y| < 0.8 
315
316         eventCutArray[0]="8000011"; photonCutArray[0] = "002000493603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Single sec  Pt cut > 0.075
317         eventCutArray[1]="8000011"; photonCutArray[1] = "002000193603007200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Single sec  Pt cut > 0.100
318         eventCutArray[2]="8000011"; photonCutArray[2] = "002000083603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Findable Cls sec  > 0.35
319         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400273102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Findable Cls prim > 0.60
320         
321         
322 } else if ( trainConfig == 17 ) {  //No eta shift |Y| < 0.8 + AddedSignals
323         
324         eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623810"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.015 < InvMass(e+,e-) < 0.050
325         eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233102623910"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + 0.025 < InvMass(e+,e-) < 0.035
326         eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603001200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100 Standard cut + qT < 0.1
327         eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011
328         
329 } else if ( trainConfig == 18 ) {  //No eta shift |Y| < 0.8 + AddedSignals
330   
331         eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102723710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Stardad cut +100 events background
332         eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233101623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  New Standard cut + Background method V0 multiplicity
333         eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102623710"; MesonCutarray[2] = "01031035000000"; //standard cut Pi0 PbPb 00-100 + No extra smearing
334         eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400253102621710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100 + Old Standard
335
336 } else if ( trainConfig == 19 ) {
337         
338         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011
339         eventCutArray[1]="8000011"; photonCutArray[1] = "032000093603007200000000"; ElecCutarray[1] = "90475400239102623710"; MesonCutarray[1] = "01033035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.6 and |Gamma_eta| < 0.65 and |e+_eta| < 0.65 and |e-_eta| < 0.65 
340         eventCutArray[2]="8000011"; photonCutArray[2] = "042000093603007200000000"; ElecCutarray[2] = "90475400235102623710"; MesonCutarray[2] = "01032035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.7 and |Gamma_eta| < 0.75 and |e+_eta| < 0.75 and |e-_eta| < 0.75
341         eventCutArray[3]="8000011"; photonCutArray[3] = "012000093603007200000000"; ElecCutarray[3] = "90475400236102623710"; MesonCutarray[3] = "01034035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.5 and |Gamma_eta| < 0.60 and |e+_eta| < 0.60 and |e-_eta| < 0.60  
342         
343 } else if ( trainConfig == 20 ) {
344         
345         eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011
346         eventCutArray[1]="8000012"; photonCutArray[1] = "032000093603007200000000"; ElecCutarray[1] = "90475400239102623710"; MesonCutarray[1] = "01033035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.6 and |Gamma_eta| < 0.65 and |e+_eta| < 0.65 and |e-_eta| < 0.65 
347         eventCutArray[2]="8000012"; photonCutArray[2] = "042000093603007200000000"; ElecCutarray[2] = "90475400235102623710"; MesonCutarray[2] = "01032035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.7 and |Gamma_eta| < 0.75 and |e+_eta| < 0.75 and |e-_eta| < 0.75
348         eventCutArray[3]="8000012"; photonCutArray[3] = "012000093603007200000000"; ElecCutarray[3] = "90475400236102623710"; MesonCutarray[3] = "01034035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 +  |Y| < 0.5 and |Gamma_eta| < 0.60 and |e+_eta| < 0.60 and |e-_eta| < 0.60  
349         
350 } else if ( trainConfig == 21 ) {
351         
352         eventCutArray[0]="8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400433102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011  +  3Cls ITS
353         eventCutArray[1]="8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400533102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011  +  4Cls ITS
354         eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400633102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011  +  5Cls ITS
355         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400733102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011  +  4Cls ITS no Any
356         
357         
358 } else if ( trainConfig == 22 ) {
359         
360         eventCutArray[0]="8000012"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400433102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011  + 3 ITScls
361         eventCutArray[1]="8000012"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400533102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011  + 4 ITScls
362         eventCutArray[2]="8000012"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400633102623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011  + 5 ITScls
363         eventCutArray[3]="8000012"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400733102623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011  + 4 ITScls no Any
364         
365 } else if ( trainConfig == 23 ) {
366   
367         eventCutArray[0]="8000011"; photonCutArray[0] = "002000493603007200000000"; ElecCutarray[0] = "90475400233102623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + Pt > 0.075
368         eventCutArray[1]="8000011"; photonCutArray[1] = "002000193603007200000000"; ElecCutarray[1] = "90475400233102623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + Pt > 0.100
369         eventCutArray[2]="8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233102633710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + Pt{e} > 0.150
370         eventCutArray[3]="8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233102653710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + Pt{e} > 0.175
371         
372   
373 } else if ( trainConfig == 24 ) {
374   
375         eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400833202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth  + new psiPair cut   0.60, 0.0 0.12
376         eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400133202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kFirst + new psiPair cut + 0.60  0.0 0.12
377         eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut   0.60, 0.0 0.12
378         eventCutArray[3] = "8000011"; photonCutArray[3] = "005000093603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut + gammaR >  10cm  0.60, 0.0 0.12
379   
380 } else if ( trainConfig == 25 ) {
381   
382         eventCutArray[0] = "8000011"; photonCutArray[0] = "008000093603007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut + gammaR >  12.5cm  0.60, 0.0 0.12
383         eventCutArray[1] = "8000011"; photonCutArray[1] = "006000093603007200000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut + gammaR >  20 cm  0.60, 0.0 0.12
384         eventCutArray[2] = "8000011"; photonCutArray[2] = "007000093603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut + gammaR >  35 cm  0.60, 0.0 0.12
385         eventCutArray[3] = "8000011"; photonCutArray[3] = "009000093603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut + gammaR >  7.5 cm  0.60, 0.0 0.12
386         
387 } else if ( trainConfig == 26 ) {
388   
389         eventCutArray[0] = "8000011"; photonCutArray[0] = "008000093603007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut + gammaR >  12.5cm  0.60, 0.0 0.12
390         eventCutArray[1] = "8000011"; photonCutArray[1] = "006000093603007200000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut + gammaR >  20 cm  0.60, 0.0 0.12
391         eventCutArray[2] = "8000011"; photonCutArray[2] = "007000093603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut + gammaR >  35 cm  0.60, 0.0 0.12
392         eventCutArray[3] = "8000011"; photonCutArray[3] = "009000093603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + new psiPair Cut + gammaR >  7.5 cm  0.60, 0.0 0.12
393         
394 } else if ( trainConfig == 27 ) {
395   
396         eventCutArray[0]= "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233202613710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Ptprim > 100 MeV
397         eventCutArray[1]= "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233202633710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Ptprim > 150 MeV
398         eventCutArray[2]= "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233202603710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Ptprim > 75 MeV
399         eventCutArray[3]= "8000011"; photonCutArray[3] = "002000493603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Ptsec  > 75 MeV
400         
401 } else if ( trainConfig == 28 ) {
402   
403         eventCutArray[0] = "8000011"; photonCutArray[0] = "002000193603007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Pt sec  > 100 MeV
404         eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233202653710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Pt prim > 175 GeV
405         eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603001200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Qt > 0.1
406         eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603002200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12 + Qt > 0.07
407         
408 } else if ( trainConfig == 29 ) {
409   
410         eventCutArray[0]= "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90375400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth + dEdx primary electron -5,5                            
411         eventCutArray[1]= "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90575400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth + dEdx primary electron -3,5
412         eventCutArray[2]= "8000011"; photonCutArray[2] = "002000091603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth + dEdx secondary electron -5,5
413         eventCutArray[3]= "8000011"; photonCutArray[3] = "002000092603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth + dEdx secondary electron -3,5
414         
415 } else if ( trainConfig == 30 ) {
416   
417         eventCutArray[0]= "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90435400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth + dEdx pion rejec primary 2.0 sigmas Low and 0 High momentum
418         eventCutArray[1]= "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90425400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth + dEdx pion rejec primary 2.0 sigmas Low and -1 High momentum 
419         eventCutArray[2]= "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90477400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth + dEdx pion rejec primary 0.3 GeV Low and 3.5 High momentum 
420         eventCutArray[3]= "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475200233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth + dEdx pion rejec primary 0.5 GeV Low and 5.0 High momentum*/
421         
422 } else if ( trainConfig == 31 ) {
423   
424         eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093653007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair cut + dEdx pion rejec sec 0.3 GeV Low and 3.5 High momentum
425         eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093601007200000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair cut + dEdx pion rejec sec 0.5 GeV Low and 5.0 High momentum
426         eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093803007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair cut + dEdx pion rejec sec  2.0 sigmas Low and  1 High momentum
427         eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007900000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair cut + Chi2  < 15
428         
429 } else if ( trainConfig == 32 ) {
430   
431         eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007800000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair cut + Chi2  < 20
432         eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007100000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair cut + Chi2  < 50
433         eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031005009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair cut + Alpha < 0.7
434         eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233002623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny   + no psi pair +  weights
435         
436 } else if ( trainConfig == 33 ) {
437   
438         eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233202622710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + DCAxy < 1 cm
439         eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400233202623810"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + 0.015 < InvMass(e+,e-) < 0.050
440         eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400233202623910"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + 0.025 < InvMass(e+,e-) < 0.035
441         eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233202723710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + 100 events background
442         
443 } else if ( trainConfig == 34 ) {
444   
445         eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233201623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + Background method V0 multiplicity
446         eventCutArray[1] = "8000011"; photonCutArray[1] = "032000093603007200000000"; ElecCutarray[1] = "90475400239202623710"; MesonCutarray[1] = "01033035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + |Y| < 0.6 and |Gamma_eta| < 0.65 and |e+_eta| < 0.65 and |e-_eta| < 0.65 
447         eventCutArray[2] = "8000011"; photonCutArray[2] = "042000093603007200000000"; ElecCutarray[2] = "90475400235202623710"; MesonCutarray[2] = "01032035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + |Y| < 0.7 and |Gamma_eta| < 0.75 and |e+_eta| < 0.75 and |e-_eta| < 0.75
448         eventCutArray[3] = "8000011"; photonCutArray[3] = "012000093603007200000000"; ElecCutarray[3] = "90475400236202623710"; MesonCutarray[3] = "01034035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + |Y| < 0.5 and |Gamma_eta| < 0.60 and |e+_eta| < 0.60 and |e-_eta| < 0.60
449         
450 } else if ( trainConfig == 35 ) {
451         
452         eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + new psiPair Cut    0.60, 0.0 0.12
453         eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603003200000000"; ElecCutarray[1] = "90475400233202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  +  Qt < 0.05
454         eventCutArray[2] = "8000011"; photonCutArray[2] = "002000083603007200000000"; ElecCutarray[2] = "90475400233202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + Findable Cls > 0.35 Secondary
455         eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400273202623710"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + Findable Cls > 0.60 primary
456         
457 } else if ( trainConfig == 36 ) {
458   
459         eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400433202623710"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + ITScls >= 3
460         eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400533202623710"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + ITScls >= 4
461         eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400633202623710"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + ITScls >= 5
462         eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233202623710"; MesonCutarray[3] = "01031035000000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny  + No extra smearing
463         
464 } else if ( trainConfig == 37 ) {
465   
466         eventCutArray[0] = "8000011"; photonCutArray[0] = "002000093603007200000000"; ElecCutarray[0] = "90475400233202623700"; MesonCutarray[0] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny    + new psiPair Cut    0.60, 0.0 0.12
467         eventCutArray[1] = "8000011"; photonCutArray[1] = "002000093603007200000000"; ElecCutarray[1] = "90475400833202623700"; MesonCutarray[1] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kBoth   + new psiPair cut    0.60, 0.0 0.12
468         eventCutArray[2] = "8000011"; photonCutArray[2] = "002000093603007200000000"; ElecCutarray[2] = "90475400133202623700"; MesonCutarray[2] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kFirst  + new psiPair Cut    0.60, 0.0 0.12
469         eventCutArray[3] = "8000011"; photonCutArray[3] = "002000093603007200000000"; ElecCutarray[3] = "90475400233002623700"; MesonCutarray[3] = "01031035009000"; //standard cut Pi0 pPb 00-100  //Tracks 2011 + kAny    + no psi pair + no weights
470         
471 }
472         
473
474
475
476
477
478
479    TList *EventCutList = new TList();
480    TList *ConvCutList  = new TList();
481    TList *MesonCutList = new TList();
482    TList *ElecCutList  = new TList();
483
484    TList *HeaderList = new TList();
485    TObjString *Header1 = new TObjString("pi0_1");
486    HeaderList->Add(Header1);
487    TObjString *Header3 = new TObjString("eta_2");
488    HeaderList->Add(Header3);
489    
490    EventCutList->SetOwner(kTRUE);
491    AliConvEventCuts **analysisEventCuts         = new AliConvEventCuts*[numberOfCuts];
492    ConvCutList->SetOwner(kTRUE);
493    AliConversionPhotonCuts **analysisCuts       = new AliConversionPhotonCuts*[numberOfCuts];
494    MesonCutList->SetOwner(kTRUE);
495    AliConversionMesonCuts **analysisMesonCuts   = new AliConversionMesonCuts*[numberOfCuts];
496    ElecCutList->SetOwner(kTRUE);
497    AliDalitzElectronCuts **analysisElecCuts     = new AliDalitzElectronCuts*[numberOfCuts];
498
499
500
501    for(Int_t i = 0; i<numberOfCuts; i++){
502
503       analysisEventCuts[i] = new AliConvEventCuts();
504
505           if (  ( trainConfig >= 1 && trainConfig <= 9 ) || trainConfig == 19  || trainConfig == 21 || trainConfig == 23 || ( trainConfig >= 24 && trainConfig <=36 )  ){
506             
507             if (doWeighting){
508               if (generatorName.CompareTo("DPMJET")==0){
509                analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "Eta_DPMJET_LHC13b2_efix_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A");
510               } else if (generatorName.CompareTo("HIJING")==0){   
511                analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_pPb_5023GeV_MBV0A", "Eta_Hijing_LHC13e7_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A");
512               }
513             }
514           }
515           else if (  ( trainConfig >= 10 && trainConfig <= 18 ) || trainConfig == 20 || trainConfig == 22 ){
516             
517              if (doWeighting){
518                 analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13e7_addSig_pPb_5023GeV_MBV0A", "Eta_Hijing_LHC13e7_addSig_pPb_5023GeV_MBV0A", "","Pi0_Fit_Data_pPb_5023GeV_MBV0A","Eta_Fit_Data_pPb_5023GeV_MBV0A");
519              }
520             
521           }
522    
523    
524        analysisEventCuts[i]->InitializeCutsFromCutString(eventCutArray[i].Data());
525        if (doEtaShiftIndCuts) {
526             analysisEventCuts[i]->DoEtaShift(doEtaShiftIndCuts);
527             analysisEventCuts[i]->SetEtaShift(stringShift);
528        }
529        EventCutList->Add(analysisEventCuts[i]);
530        analysisEventCuts[i]->SetFillCutHistograms("",kFALSE);
531        analysisEventCuts[i]->SetAcceptedHeader(HeaderList);
532        
533        analysisCuts[i] = new AliConversionPhotonCuts();
534        
535        if( ! analysisCuts[i]->InitializeCutsFromCutString(photonCutArray[i].Data()) ) {
536             cout<<"ERROR: analysisCuts [" <<i<<"]"<<endl;
537             return 0;
538        }
539        analysisCuts[i]->SetIsHeavyIon(isHeavyIon);
540        ConvCutList->Add(analysisCuts[i]);
541        analysisCuts[i]->SetFillCutHistograms("",kFALSE);
542                 
543        
544               
545
546       analysisMesonCuts[i] = new AliConversionMesonCuts();
547     
548       if( ! analysisMesonCuts[i]->InitializeCutsFromCutString(MesonCutarray[i].Data()) ) {
549             cout<<"ERROR: analysisMesonCuts [ " <<i<<" ] "<<endl;
550             return 0;
551       }
552       MesonCutList->Add(analysisMesonCuts[i]);
553       analysisMesonCuts[i]->SetFillCutHistograms("");
554       
555
556
557       //TString cutName( Form("%s_%s_%s",photonCutArray[i].Data(),ElecCutarray[i].Data(),MesonCutarray[i].Data() ) );
558
559
560        analysisElecCuts[i] = new AliDalitzElectronCuts();
561        if( !analysisElecCuts[i]->InitializeCutsFromCutString(ElecCutarray[i].Data())) {
562
563             cout<< "ERROR:  analysisElecCuts [ " <<i<<" ] "<<endl;
564             return 0;
565        }
566        ElecCutList->Add(analysisElecCuts[i]);
567        analysisElecCuts[i]->SetFillCutHistograms("",kFALSE,ElecCutarray[i].Data()); 
568        
569
570    }
571
572    task->SetEventCutList(numberOfCuts,EventCutList);
573    task->SetConversionCutList(numberOfCuts,ConvCutList);
574    task->SetMesonCutList(MesonCutList);
575    task->SetElectronCutList(ElecCutList);
576
577    task->SetMoveParticleAccordingToVertex(kTRUE);
578    task->SetProductionVertextoVGamma(kTRUE);
579
580
581    if(enableQAMesonTask) task->SetDoMesonQA(kTRUE);
582    if(enableDoMesonChic) task->SetDoChicAnalysis(kTRUE);
583
584    //connect containers
585    AliAnalysisDataContainer *coutput =
586    mgr->CreateContainer(Form("GammaConvDalitzV1_%i",trainConfig), TList::Class(),
587                            AliAnalysisManager::kOutputContainer,Form("GammaConvV1Dalitz_%i.root",trainConfig));
588
589    mgr->AddTask(task);
590    mgr->ConnectInput(task,0,cinput);
591    mgr->ConnectOutput(task,1,coutput);
592
593    return;
594
595 }