added libraries to AddTask in order to include flow package
[u/mrichter/AliRoot.git] / PWGGA / GammaConv / macros / AddTask_GammaConvV1_PbPb.C
1 void AddTask_GammaConvV1_PbPb(  Int_t trainConfig = 1,  //change different set of cuts
2                               Bool_t isMC   = kFALSE, //run MC 
3                               Int_t enableQAMesonTask = 0, //enable QA in AliAnalysisTaskGammaConvV1
4                               Int_t enableQAPhotonTask = 0, // enable additional QA task
5                               TString fileNameInputForWeighting = "MCSpectraInput.root", // path to file for weigting input
6                               Int_t headerSelectionInt = 0,  // 1 pi0 header, 2 eta header, 3 both (only for "named" boxes)
7                               TString cutnumberAODBranch = "1000000060084000001500000",
8                               TString periodName = "LHC13d2",  //name of the period for added signals and weighting
9                               Bool_t doWeighting = kFALSE  //enable Weighting
10                            ) {
11
12         // ================= Load Librariers =================================
13         gSystem->Load("libCore.so");  
14         gSystem->Load("libTree.so");
15         gSystem->Load("libGeom.so");
16         gSystem->Load("libVMC.so");
17         gSystem->Load("libPhysics.so");
18         gSystem->Load("libMinuit");
19         gSystem->Load("libSTEERBase");
20         gSystem->Load("libESD");
21         gSystem->Load("libAOD");
22         gSystem->Load("libANALYSIS");
23         gSystem->Load("libANALYSISalice");  
24         gSystem->Load("libCDB.so");
25         gSystem->Load("libSTEER.so");
26         gSystem->Load("libSTEERBase.so");
27         gSystem->Load("libTENDER.so");
28         gSystem->Load("libTENDERSupplies.so");
29         gSystem->Load("libPWGflowBase.so");
30         gSystem->Load("libPWGflowTasks.so");
31         gSystem->Load("libPWGGAGammaConv.so");
32
33         Int_t isHeavyIon = 1;
34
35         // ================== GetAnalysisManager ===============================
36         AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
37         if (!mgr) {
38                 Error(Form("AddTask_GammaConvV1_%i",trainConfig), "No analysis manager found.");
39                 return ;
40         }
41
42         // ================== GetInputEventHandler =============================
43         AliVEventHandler *inputHandler=mgr->GetInputEventHandler();
44
45         //========= Add PID Reponse to ANALYSIS manager ====
46         if(!(AliPIDResponse*)mgr->GetTask("PIDResponseTask")){
47                 gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPIDResponse.C");
48                 AddTaskPIDResponse(isMC);
49         }
50
51         //=========  Set Cutnumber for V0Reader ================================
52         TString cutnumberPhoton = "000084001001500000000";
53         TString cutnumberEvent = "1000000";
54         AliAnalysisDataContainer *cinput = mgr->GetCommonInputContainer();
55         //========= Add V0 Reader to  ANALYSIS manager if not yet existent =====
56         if( !(AliV0ReaderV1*)mgr->GetTask("V0ReaderV1") ){
57                 AliV0ReaderV1 *fV0ReaderV1 = new AliV0ReaderV1("V0ReaderV1");
58                 
59                 fV0ReaderV1->SetUseOwnXYZCalculation(kTRUE);
60                 fV0ReaderV1->SetCreateAODs(kFALSE);// AOD Output
61                 fV0ReaderV1->SetUseAODConversionPhoton(kTRUE);
62
63                 if (!mgr) {
64                         Error("AddTask_V0ReaderV1", "No analysis manager found.");
65                         return;
66                 }
67
68                 AliConvEventCuts *fEventCuts=NULL;
69                 if(cutnumberEvent!=""){
70                         fEventCuts= new AliConvEventCuts(cutnumberEvent.Data(),cutnumberEvent.Data());
71                         fEventCuts->SetPreSelectionCutFlag(kTRUE);
72                         if(fEventCuts->InitializeCutsFromCutString(cutnumberEvent.Data())){
73                                 fV0ReaderV1->SetEventCuts(fEventCuts);
74                                 fEventCuts->SetFillCutHistograms("",kTRUE);
75                         }
76                 }
77
78                 
79                 // Set AnalysisCut Number
80                 AliConversionPhotonCuts *fCuts=NULL;
81                 if(cutnumberPhoton!=""){
82                         fCuts= new AliConversionPhotonCuts(cutnumberPhoton.Data(),cutnumberPhoton.Data());
83                         fCuts->SetPreSelectionCutFlag(kTRUE);
84                         fCuts->SetIsHeavyIon(isHeavyIon);
85                         if(fCuts->InitializeCutsFromCutString(cutnumberPhoton.Data())){
86                                 fV0ReaderV1->SetConversionCuts(fCuts);
87                                 fCuts->SetFillCutHistograms("",kTRUE);
88                         }
89                 }
90                 
91                 if(inputHandler->IsA()==AliAODInputHandler::Class()){
92                 // AOD mode
93                         fV0ReaderV1->SetDeltaAODBranchName(Form("GammaConv_%s_gamma",cutnumberAODBranch.Data()));
94                 }
95                 fV0ReaderV1->Init();
96
97                 AliLog::SetGlobalLogLevel(AliLog::kInfo);
98
99                 //connect input V0Reader
100                 mgr->AddTask(fV0ReaderV1);
101                 mgr->ConnectInput(fV0ReaderV1,0,cinput);
102
103         }
104
105         //================================================
106         //========= Add task to the ANALYSIS manager =====
107         //================================================
108         AliAnalysisTaskGammaConvV1 *task=NULL;
109         task= new AliAnalysisTaskGammaConvV1(Form("GammaConvV1_%i",trainConfig));
110         task->SetIsHeavyIon(isHeavyIon);
111         task->SetIsMC(isMC);
112         // Cut Numbers to use in Analysis
113         Int_t numberOfCuts = 5;
114
115         TString *eventCutArray = new TString[numberOfCuts];
116         TString *photonCutArray = new TString[numberOfCuts];
117         TString *mesonCutArray = new TString[numberOfCuts];
118
119         if (trainConfig == 1){ // Standard cuts
120                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
121                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "042092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
122                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "042092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
123                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "042092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
124                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "042092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
125         } else if (trainConfig == 2) { // Standard cuts
126                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
127                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "042092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
128                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "042092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
129                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "042092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
130                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "042092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
131         } else if (trainConfig == 3) { // Standard cuts only added signals
132                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
133                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "042092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
134                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "042092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
135                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "042092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
136                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "042092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
137         } else if (trainConfig == 4) { // Standard cuts only added signals
138                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "042092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
139                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "042092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
140                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "042092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
141                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "042092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 20-40% 
142                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "042092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
143         } else if (trainConfig == 5){ // R-minCut 7.5 cm
144                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "049092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
145                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "049092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
146                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "049092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
147                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "049092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
148                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "049092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
149         } else if (trainConfig == 6) { // R-minCut 7.5 cm
150                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "049092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
151                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "049092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
152                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "049092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
153                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "049092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
154                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "049092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
155         } else if (trainConfig == 7) {// R-minCut 7.5 cm
156                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "049092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
157                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "049092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
158                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "049092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
159                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "049092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
160                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "049092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
161         } else if (trainConfig == 8) { // R-minCut 7.5 cm
162                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "040092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
163                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "049092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
164                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "049092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
165                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "049092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 20-40% 
166                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "049092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%   
167         } else if (trainConfig == 9){ // R-minCut 12.5 cm
168                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "048092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
169                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "048092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
170                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "048092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
171                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "048092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
172                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "048092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
173         } else if (trainConfig == 10) { // R-minCut 12.5 cm
174                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "048092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
175                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "048092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
176                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "048092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
177                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "048092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 40-80%
178                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "048092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%
179         } else if (trainConfig == 11) {// R-minCut 12.5 cm
180                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "048092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 0-5%
181                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "048092970023220000000"; mesonCutArray[ 1] = "01522045009000"; // 5-10%
182                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "048092970023220000000"; mesonCutArray[ 2] = "01522045009000"; // 0-10%
183                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "048092970023220000000"; mesonCutArray[ 3] = "01522045009000"; // 10-20%
184                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "048092970023220000000"; mesonCutArray[ 4] = "01522045009000"; // 0-20%
185         } else if (trainConfig == 12) { // R-minCut 12.5 cm
186                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "040092970023220000000"; mesonCutArray[ 0] = "01522045009000"; // 20-40%
187                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "048092970023220000000"; mesonCutArray[ 1] = "01522065009000"; // 40-60%
188                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "048092970023220000000"; mesonCutArray[ 2] = "01522065009000"; // 60-80%
189                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "048092970023220000000"; mesonCutArray[ 3] = "01522065009000"; // 20-40% 
190                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "048092970023220000000"; mesonCutArray[ 4] = "01522065009000"; // 40-90%      
191         } else  if (trainConfig == 13){ // LHC10h standard, eta 0.65, y = 0.6 
192                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
193                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
194                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
195                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "032092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
196                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "032092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
197         } else if (trainConfig == 14) {  // LHC10h standard, eta 0.65, y = 0.6 
198                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
199                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "032092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
200                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "032092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
201                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "032092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
202                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "032092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
203         } else if (trainConfig == 15) { // LHC10h standard, eta 0.65, y = 0.6  added signals
204                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
205                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
206                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
207                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "032092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
208                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "032092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
209         } else if (trainConfig == 16) { // LHC10h standard, eta 0.65, y = 0.6  added signals
210                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "032092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
211                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "032092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
212                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "032092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
213                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "032092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
214                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "032092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
215         } else  if (trainConfig == 17){ // LHC10h standard, eta 0.65, y = 0.6, photon quality 1 
216                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
217                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
218                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
219                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "032092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
220                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "032092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
221         } else if (trainConfig == 18) {  // LHC10h standard, eta 0.65, y = 0.6, photon quality 1 
222                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
223                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "032092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
224                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "032092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
225                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "032092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
226                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "032092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
227         } else if (trainConfig == 19) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 1  added signal
228                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
229                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
230                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
231                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "032092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
232                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "032092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
233         } else if (trainConfig == 20) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 1 added signal
234                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "032092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
235                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "032092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
236                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "032092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
237                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "032092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
238                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "032092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
239         } else  if (trainConfig == 21){ // LHC10h standard, eta 0.65, y = 0.6, photon quality 2 
240                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
241                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970023220030000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
242                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970023220030000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
243                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "032092970023220030000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
244                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "032092970023220030000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
245         } else if (trainConfig == 22) {  // LHC10h standard, eta 0.65, y = 0.6, photon quality 2 
246                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
247                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "032092970023220030000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
248                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "032092970023220030000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
249                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "032092970023220030000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
250                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "032092970023220030000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
251         } else if (trainConfig == 23) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 2  added signal
252                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
253                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970023220030000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
254                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970023220030000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
255                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "032092970023220030000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
256                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "032092970023220030000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
257         } else if (trainConfig == 24) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 2 added signal
258                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "032092970023220030000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
259                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "032092970023220030000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
260                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "032092970023220030000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
261                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "032092970023220030000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
262                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "032092970023220030000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
263         } else  if (trainConfig == 25){ // LHC10h standard, eta 0.65, y = 0.6, photon quality 3 
264                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
265                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
266                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
267                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "032092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
268                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "032092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
269         } else if (trainConfig == 26) {  // LHC10h standard, eta 0.65, y = 0.6, photon quality 3 
270                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
271                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "032092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
272                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "032092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
273                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "032092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
274                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "032092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
275         } else if (trainConfig == 27) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 3  added signal
276                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
277                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
278                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
279                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "032092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
280                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "032092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
281         } else if (trainConfig == 28) { // LHC10h standard, eta 0.65, y = 0.6, photon quality 3 added signal
282                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "032092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
283                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "032092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
284                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "032092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
285                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "032092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
286                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "032092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
287         } else  if (trainConfig == 29){ // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm 
288                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
289                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "037092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
290                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "037092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
291                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "037092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
292                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "037092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
293         } else if (trainConfig == 30) {  // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm 
294                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
295                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "037092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
296                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "037092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
297                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "037092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
298                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "037092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
299         } else if (trainConfig == 31) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm  added signals
300                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
301                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "037092970023220000000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
302                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "037092970023220000000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
303                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "037092970023220000000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
304                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "037092970023220000000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
305         } else if (trainConfig == 32) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm  added signals
306                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "037092970023220000000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
307                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "037092970023220000000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
308                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "037092970023220000000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
309                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "037092970023220000000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
310                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "037092970023220000000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
311         } else  if (trainConfig == 33){ // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1 
312                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
313                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "037092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
314                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "037092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
315                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "037092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
316                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "037092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
317         } else if (trainConfig == 34) {  // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1  
318                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
319                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "037092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
320                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "037092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
321                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "037092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
322                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "037092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
323         } else if (trainConfig == 35) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1  added signals
324                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
325                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "037092970023220020000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
326                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "037092970023220020000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
327                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "037092970023220020000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
328                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "037092970023220020000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
329         } else if (trainConfig == 36) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 1  added signals
330                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "037092970023220020000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
331                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "037092970023220020000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
332                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "037092970023220020000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
333                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "037092970023220020000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
334                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "037092970023220020000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
335         } else  if (trainConfig == 37){ // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3 
336                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
337                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "037092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
338                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "037092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
339                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "037092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
340                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "037092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
341         } else if (trainConfig == 38) {  // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3  
342                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
343                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "037092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
344                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "037092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
345                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "037092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 40-80%
346                 eventCutArray[ 4] = "5490001"; photonCutArray[ 4] = "037092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
347         } else if (trainConfig == 39) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3  added signals
348                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 0-5%
349                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "037092970023220040000"; mesonCutArray[ 1] = "01523045009000"; // 5-10%
350                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "037092970023220040000"; mesonCutArray[ 2] = "01523045009000"; // 0-10%
351                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "037092970023220040000"; mesonCutArray[ 3] = "01523045009000"; // 10-20%
352                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "037092970023220040000"; mesonCutArray[ 4] = "01523045009000"; // 0-20%
353         } else if (trainConfig == 40) { // LHC10h standard, eta 0.65, y = 0.6, min R = 35 cm, photon quality 3  added signals
354                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "037092970023220040000"; mesonCutArray[ 0] = "01523045009000"; // 20-40%
355                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "037092970023220040000"; mesonCutArray[ 1] = "01523065009000"; // 40-60%
356                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "037092970023220040000"; mesonCutArray[ 2] = "01523065009000"; // 60-80%
357                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "037092970023220040000"; mesonCutArray[ 3] = "01523065009000"; // 20-40% 
358                 eventCutArray[ 4] = "5490002"; photonCutArray[ 4] = "037092970023220040000"; mesonCutArray[ 4] = "01523065009000"; // 40-90%
359         } else if (trainConfig == 41){ // Standard cuts, eta 0.9, only to be run on data :kSemiCentral
360                 eventCutArray[ 0] = "6016001"; photonCutArray[ 0] = "002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; // 0-5%
361                 eventCutArray[ 1] = "6126001"; photonCutArray[ 1] = "002092970023220000000"; mesonCutArray[ 1] = "01525045009000"; // 5-10%
362                 eventCutArray[ 2] = "5016001"; photonCutArray[ 2] = "002092970023220000000"; mesonCutArray[ 2] = "01525045009000"; // 0-10%
363                 eventCutArray[ 3] = "6126001"; photonCutArray[ 3] = "002092970023220000000"; mesonCutArray[ 3] = "01525045009000"; // 10-20%
364                 eventCutArray[ 4] = "5026001"; photonCutArray[ 4] = "002092970023220000000"; mesonCutArray[ 4] = "01525045009000"; // 0-20%
365         } else if (trainConfig == 42) { // Standard cuts, eta 0.9, only to be run on data :kSemiCentral
366                 eventCutArray[ 0] = "5236001"; photonCutArray[ 0] = "002092970023220000000"; mesonCutArray[ 0] = "01525065009000"; 
367                 eventCutArray[ 1] = "5346001"; photonCutArray[ 1] = "002092970023220000000"; mesonCutArray[ 1] = "01525065009000"; 
368                 eventCutArray[ 2] = "5456001"; photonCutArray[ 2] = "002092970023220000000"; mesonCutArray[ 2] = "01525065009000"; 
369                 eventCutArray[ 3] = "5566001"; photonCutArray[ 3] = "002092970023220000000"; mesonCutArray[ 3] = "01525065009000"; 
370                 eventCutArray[ 4] = "5676001"; photonCutArray[ 4] = "002092970023220000000"; mesonCutArray[ 4] = "01525065009000"; 
371         } else if (trainConfig == 43){ // Standard cuts, eta 0.9, only to be run on data
372                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970023220000000"; mesonCutArray[ 0] = "01525045009000"; 
373                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970023220000000"; mesonCutArray[ 1] = "01525065009000"; 
374                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970023220000000"; mesonCutArray[ 2] = "01525065009000"; 
375                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970023220000000"; mesonCutArray[ 3] = "01525065009000"; 
376                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970023220000000"; mesonCutArray[ 4] = "01525065009000"; 
377         } else if ( trainConfig == 44){ // qt elipse cut 0.05
378                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028220000000"; mesonCutArray[ 0] = "01525065009000"; // 0-5%
379                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028220000000"; mesonCutArray[ 1] = "01525065009000"; // 5-10%
380                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028220000000"; mesonCutArray[ 2] = "01525065009000"; // 0-10%
381                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028220000000"; mesonCutArray[ 3] = "01525065009000"; // 10-20%
382                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970028220000000"; mesonCutArray[ 4] = "01525065009000"; // 0-20%
383         } else if ( trainConfig == 45) { // qt elipse cut 0.05
384                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970028220000000"; mesonCutArray[ 0] = "01525065009000"; // 20-40%
385                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970028220000000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
386                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970028220000000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
387                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028220000000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
388                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970028220000000"; mesonCutArray[ 4] = "01525065009000"; // 30-50%  
389         } else if ( trainConfig == 46){ // qt elipse cut 0.05
390                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970028220000000"; mesonCutArray[ 0] = "01525065009000"; // 20-30% 
391                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970028220000000"; mesonCutArray[ 1] = "01525065009000"; // 30-40% 
392                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970028220000000"; mesonCutArray[ 2] = "01525065009000"; // 40-50% 
393                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970028220000000"; mesonCutArray[ 3] = "01525065009000"; // 50-60% 
394                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028220000000"; mesonCutArray[ 4] = "01525065009000"; // 60-70% 
395         } else if ( trainConfig == 47){ // cos(theta_point) cut 0.85
396                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970023220400000"; mesonCutArray[ 0] = "01525065009000"; // 0-5%
397                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970023220400000"; mesonCutArray[ 1] = "01525065009000"; // 5-10%
398                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970023220400000"; mesonCutArray[ 2] = "01525065009000"; // 0-10%
399                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970023220400000"; mesonCutArray[ 3] = "01525065009000"; // 10-20%
400                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970023220400000"; mesonCutArray[ 4] = "01525065009000"; // 0-20%
401         } else if ( trainConfig == 48) { // cos(theta_point) cut 0.85
402                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970023220400000"; mesonCutArray[ 0] = "01525065009000"; // 20-40%
403                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970023220400000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
404                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970023220400000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
405                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970023220400000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
406                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970023220400000"; mesonCutArray[ 4] = "01525065009000"; // 30-50%  
407         } else if ( trainConfig == 49){ // cos(theta_point) cut 0.85
408                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970023220400000"; mesonCutArray[ 0] = "01525065009000"; // 20-30% 
409                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970023220400000"; mesonCutArray[ 1] = "01525065009000"; // 30-40% 
410                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970023220400000"; mesonCutArray[ 2] = "01525065009000"; // 40-50% 
411                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970023220400000"; mesonCutArray[ 3] = "01525065009000"; // 50-60% 
412                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970023220400000"; mesonCutArray[ 4] = "01525065009000"; // 60-70% 
413         } else if ( trainConfig == 50){ // psi pair 2D 0.05
414                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970023260000000"; mesonCutArray[ 0] = "01525065009000"; // 0-5%
415                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970023260000000"; mesonCutArray[ 1] = "01525065009000"; // 5-10%
416                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970023260000000"; mesonCutArray[ 2] = "01525065009000"; // 0-10%
417                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970023260000000"; mesonCutArray[ 3] = "01525065009000"; // 10-20%
418                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970023260000000"; mesonCutArray[ 4] = "01525065009000"; // 0-20%
419         } else if ( trainConfig == 51) { // psi pair 2D 0.05
420                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970023260000000"; mesonCutArray[ 0] = "01525065009000"; // 20-40%
421                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970023260000000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
422                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970023260000000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
423                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970023260000000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
424                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970023260000000"; mesonCutArray[ 4] = "01525065009000"; // 30-50%  
425         } else if ( trainConfig == 52){ // psi pair 2D 0.05
426                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970023260000000"; mesonCutArray[ 0] = "01525065009000"; // 20-30% 
427                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970023260000000"; mesonCutArray[ 1] = "01525065009000"; // 30-40% 
428                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970023260000000"; mesonCutArray[ 2] = "01525065009000"; // 40-50% 
429                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970023260000000"; mesonCutArray[ 3] = "01525065009000"; // 50-60% 
430                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970023260000000"; mesonCutArray[ 4] = "01525065009000"; // 60-70%       
431         } else if ( trainConfig == 53){ // psi pair 2D 0.1
432                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970023250000000"; mesonCutArray[ 0] = "01525065009000"; // 0-5%
433                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970023250000000"; mesonCutArray[ 1] = "01525065009000"; // 5-10%
434                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970023250000000"; mesonCutArray[ 2] = "01525065009000"; // 0-10%
435                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970023250000000"; mesonCutArray[ 3] = "01525065009000"; // 10-20%
436                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970023250000000"; mesonCutArray[ 4] = "01525065009000"; // 0-20%
437         } else if ( trainConfig == 54) { // psi pair 2D 0.1
438                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970023250000000"; mesonCutArray[ 0] = "01525065009000"; // 20-40%
439                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970023250000000"; mesonCutArray[ 1] = "01525065009000"; // 40-60%
440                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970023250000000"; mesonCutArray[ 2] = "01525065009000"; // 60-80%
441                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970023250000000"; mesonCutArray[ 3] = "01525065009000"; // 40-80%
442                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970023250000000"; mesonCutArray[ 4] = "01525065009000"; // 30-50%  
443         } else if ( trainConfig == 55){ // psi pair 2D 0.1
444                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970023250000000"; mesonCutArray[ 0] = "01525065009000"; // 20-30% 
445                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970023250000000"; mesonCutArray[ 1] = "01525065009000"; // 30-40% 
446                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970023250000000"; mesonCutArray[ 2] = "01525065009000"; // 40-50% 
447                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970023250000000"; mesonCutArray[ 3] = "01525065009000"; // 50-60% 
448                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970023250000000"; mesonCutArray[ 4] = "01525065009000"; // 60-70%          
449         } else if ( trainConfig == 56){ // cleaner cuts
450                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
451                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
452                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
453                 eventCutArray[ 3] = "5120001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
454                 eventCutArray[ 4] = "5020001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
455         } else if ( trainConfig == 57){ // cleaner cuts added signal
456                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
457                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
458                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
459                 eventCutArray[ 3] = "5120002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
460                 eventCutArray[ 4] = "5020002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%   
461         } else if ( trainConfig == 58) { // cleaner cuts
462                 eventCutArray[ 0] = "5240001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 20-40%
463                 eventCutArray[ 1] = "5460001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 40-60%
464                 eventCutArray[ 2] = "5680001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 60-80%
465                 eventCutArray[ 3] = "5480001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 40-80%
466                 eventCutArray[ 4] = "5350001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 30-50%  
467         } else if ( trainConfig == 59) { // cleaner cuts added signal
468                 eventCutArray[ 0] = "5240002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 20-40%
469                 eventCutArray[ 1] = "5460002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 40-60%
470                 eventCutArray[ 2] = "5680002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 60-80%
471                 eventCutArray[ 3] = "5480002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 40-80%
472                 eventCutArray[ 4] = "5350002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 30-50%            
473         } else if ( trainConfig == 60){ // cleaner cuts
474                 eventCutArray[ 0] = "5230001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 20-30% 
475                 eventCutArray[ 1] = "5340001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 30-40% 
476                 eventCutArray[ 2] = "5450001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 40-50% 
477                 eventCutArray[ 3] = "5560001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 50-60% 
478                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 60-70%                
479         } else if ( trainConfig == 61){ // cleaner cuts added signal
480                 eventCutArray[ 0] = "5230002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 20-30% 
481                 eventCutArray[ 1] = "5340002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 30-40% 
482                 eventCutArray[ 2] = "5450002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 40-50% 
483                 eventCutArray[ 3] = "5560002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 50-60% 
484                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 60-70%                   
485         } else if ( trainConfig == 62){ // cleaner cuts
486                 eventCutArray[ 0] = "6230001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
487                 eventCutArray[ 1] = "6340001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
488                 eventCutArray[ 2] = "6450001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
489                 eventCutArray[ 3] = "6560001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
490                 eventCutArray[ 4] = "6670001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
491         } else if ( trainConfig == 63){ // cleaner cuts added signal
492                 eventCutArray[ 0] = "6230002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
493                 eventCutArray[ 1] = "6340002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
494                 eventCutArray[ 2] = "6450002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
495                 eventCutArray[ 3] = "6560002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
496                 eventCutArray[ 4] = "6670002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
497         } else if ( trainConfig == 64){ // cleaner cuts
498                 eventCutArray[ 0] = "6780001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
499                 eventCutArray[ 1] = "6890001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
500                 eventCutArray[ 2] = "5670001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
501                 eventCutArray[ 3] = "5780001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
502                 eventCutArray[ 4] = "5890001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
503         } else if ( trainConfig == 65){ // cleaner cuts added signal
504                 eventCutArray[ 0] = "6780002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
505                 eventCutArray[ 1] = "6890002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
506                 eventCutArray[ 2] = "5670002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
507                 eventCutArray[ 3] = "5780002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
508                 eventCutArray[ 4] = "5890002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
509         } else if ( trainConfig == 66){ // cleaner cuts
510                 eventCutArray[ 0] = "7010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
511                 eventCutArray[ 1] = "7120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
512                 eventCutArray[ 2] = "7230001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
513                 eventCutArray[ 3] = "7340001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
514                 eventCutArray[ 4] = "7450001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
515         } else if ( trainConfig == 67){ // cleaner cuts added signal
516                 eventCutArray[ 0] = "7010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
517                 eventCutArray[ 1] = "7120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
518                 eventCutArray[ 2] = "7230002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
519                 eventCutArray[ 3] = "7340002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
520                 eventCutArray[ 4] = "7450002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
521         } else if ( trainConfig == 68){ // cleaner cuts
522                 eventCutArray[ 0] = "7560001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
523                 eventCutArray[ 1] = "7670001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
524                 eventCutArray[ 2] = "7780001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
525                 eventCutArray[ 3] = "7890001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 10-20%
526                 eventCutArray[ 4] = "7090001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 0-20%      
527         } else if ( trainConfig == 69){ // cleaner cuts added signal
528                 eventCutArray[ 0] = "7560002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
529                 eventCutArray[ 1] = "7670002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
530                 eventCutArray[ 2] = "7780002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
531                 eventCutArray[ 3] = "7890002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 10-20%
532                 eventCutArray[ 4] = "7090002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 0-20%                      
533         } else if ( trainConfig == 70){ // variation eta  0.65
534                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "032092970028250400000"; mesonCutArray[ 0]= "01523065000000"; // 0-5%
535                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "032092970028250400000"; mesonCutArray[ 1]= "01523065000000"; // 5-10%
536                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "032092970028250400000"; mesonCutArray[ 2]= "01523065000000"; // 0-10%
537                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "032092970028250400000"; mesonCutArray[ 3]= "01523065000000"; // 20-40%
538                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "032092970028250400000"; mesonCutArray[ 4]= "01523065000000"; // 20-50% 
539         } else if ( trainConfig == 71){ // variation eta  0.65 added signal
540                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "032092970028250400000"; mesonCutArray[ 0]= "01523065000000"; // 0-5%
541                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "032092970028250400000"; mesonCutArray[ 1]= "01523065000000"; // 5-10%
542                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "032092970028250400000"; mesonCutArray[ 2]= "01523065000000"; // 0-10%
543                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "032092970028250400000"; mesonCutArray[ 3]= "01523065000000"; // 20-40%
544                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "032092970028250400000"; mesonCutArray[ 4]= "01523065000000"; // 20-50%             
545         } else if ( trainConfig == 72){ // variation eta  0.75
546                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "042092970028250400000"; mesonCutArray[ 0]= "01522065009000"; // 0-5%
547                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "042092970028250400000"; mesonCutArray[ 1]= "01522065009000"; // 5-10%
548                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "042092970028250400000"; mesonCutArray[ 2]= "01522065009000"; // 0-10%
549                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "042092970028250400000"; mesonCutArray[ 3]= "01522065009000"; // 20-40%
550                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "042092970028250400000"; mesonCutArray[ 4]= "01522065009000"; // 20-50% 
551         } else if ( trainConfig == 73){ // variation eta  0.75 added signal
552                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "042092970028250400000"; mesonCutArray[ 0]= "01522065009000"; // 0-5%
553                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "042092970028250400000"; mesonCutArray[ 1]= "01522065009000"; // 5-10%
554                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "042092970028250400000"; mesonCutArray[ 2]= "01522065009000"; // 0-10%
555                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "042092970028250400000"; mesonCutArray[ 3]= "01522065009000"; // 20-40%
556                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "042092970028250400000"; mesonCutArray[ 4]= "01522065009000"; // 20-50% 
557         } else if ( trainConfig == 74){ // single pt 0.075
558                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002492970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
559                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002492970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
560                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002492970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
561                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002492970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
562                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002492970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
563         } else if ( trainConfig == 75){ // single pt 0.075 added signal
564                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002492970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
565                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002492970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
566                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002492970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
567                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002492970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
568                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002492970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
569         } else if ( trainConfig == 76){ // single pt 0.1
570                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002192970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
571                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002192970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
572                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002192970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
573                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002192970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
574                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002192970028250400000"; mesonCutArray[ 4]= "01525065000000"; //20-50%
575         } else if ( trainConfig == 77){ // single pt 0.1 added signal
576                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002192970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
577                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002192970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
578                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002192970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
579                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002192970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
580                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002192970028250400000"; mesonCutArray[ 4]= "01525065000000"; //20-50%
581         } else if ( trainConfig == 78){ // variation TPC cls 0.7
582                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002062970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
583                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002062970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
584                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002062970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
585                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002062970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
586                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002062970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
587         } else if ( trainConfig == 79){ // variation TPC cls 0.7 added signal
588                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002062970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
589                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002062970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
590                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002062970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
591                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002062970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
592                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002062970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
593         } else if ( trainConfig == 80){ // variation TPC cls 0.35
594                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002082970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
595                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002082970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
596                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002082970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
597                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002082970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
598                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002082970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
599         } else if ( trainConfig == 81){ // variation TPC cls 0.35 added signal
600                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002082970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
601                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002082970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
602                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002082970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
603                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002082970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
604                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002082970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
605         } else if ( trainConfig == 82){ // variation edEdx  -4,5
606                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002093970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
607                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002093970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
608                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002093970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
609                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002093970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
610                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002093970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
611         } else if ( trainConfig == 83){ // variation edEdx  -4,5  added signal
612                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002093970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
613                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002093970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
614                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002093970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
615                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002093970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
616                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002093970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
617         } else if ( trainConfig == 84){ // variation edEdx  -2.5,4
618                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002096970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
619                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002096970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
620                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002096970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
621                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002096970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
622                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002096970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
623         } else if ( trainConfig == 85){ // variation edEdx  -2.5,4  added signal
624                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002096970028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
625                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002096970028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
626                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002096970028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
627                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002096970028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
628                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002096970028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
629         } else if ( trainConfig == 86){ //variation pion p dEdx 0.3-5.
630                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092951028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
631                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092951028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
632                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092951028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
633                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092951028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
634                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092951028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
635         } else if ( trainConfig == 87){ //variation pion p dEdx 0.3-5.  added signal
636                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092951028250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
637                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092951028250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
638                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092951028250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
639                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092951028250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
640                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092951028250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
641         } else if ( trainConfig == 88){ // TOF el. PID -3,5
642                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970038250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
643                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970038250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
644                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970038250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
645                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970038250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
646                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970038250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
647         } else if ( trainConfig == 89){ // TOF el. PID -3,5  added signal
648                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970038250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
649                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970038250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
650                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970038250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
651                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970038250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
652                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970038250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%             
653         } else if ( trainConfig == 90){ // TOF el. PID -2,3
654                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970048250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
655                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970048250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
656                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970048250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
657                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970048250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
658                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970048250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
659         } else if ( trainConfig == 91){ // TOF el. PID -2,3  added signal
660                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970048250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
661                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970048250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
662                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970048250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
663                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970048250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
664                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970048250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
665         } else if ( trainConfig == 92){ // qt 0.03
666                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970029250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
667                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970029250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
668                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970029250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
669                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970029250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
670                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970029250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
671         } else if ( trainConfig == 93){ // qt 0.03  added signal
672                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970029250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
673                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970029250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
674                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970029250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
675                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970029250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
676                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970029250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
677         } else if ( trainConfig == 94){ // qt 0.07 no2D
678                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970022250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
679                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970022250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
680                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970022250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
681                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970022250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
682                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970022250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
683         } else if ( trainConfig         == 95){ // qt 0.07 no2D  added signal
684                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970022250400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
685                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970022250400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
686                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970022250400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
687                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970022250400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
688                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970022250400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
689         } else if ( trainConfig == 96){ // chi2  50.
690                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028150400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
691                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028150400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
692                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028150400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
693                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028150400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
694                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028150400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
695         } else if ( trainConfig == 97){ // chi2  50.  added signal
696                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028150400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
697                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028150400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
698                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028150400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
699                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028150400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
700                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028150400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%             
701         } else if ( trainConfig == 98){ // chi2  20.
702                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028850400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
703                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028850400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
704                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028850400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
705                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028850400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
706                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028850400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
707         } else if ( trainConfig == 99){ // chi2  20.  added signal
708                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028850400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
709                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028850400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
710                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028850400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
711                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028850400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
712                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028850400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
713         } else if ( trainConfig == 100){ // psi pair 0.05
714                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028260400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
715                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028260400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
716                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028260400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
717                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028260400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
718                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028260400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
719         } else if ( trainConfig == 101){ // psi pair 0.05  added signal
720                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028260400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
721                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028260400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
722                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028260400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
723                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028260400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
724                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028260400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
725         } else if ( trainConfig == 102){ // cosPA -1
726                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250000000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
727                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250000000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
728                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250000000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
729                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250000000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
730                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250000000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
731         } else if ( trainConfig == 103){ // cosPA -1  added signal
732                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250000000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
733                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250000000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
734                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250000000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
735                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250000000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
736                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250000000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
737         } else if ( trainConfig == 104){ // variation alpha 0.75
738                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525055000000"; // 0-5%
739                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525055000000"; // 5-10%
740                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525055000000"; // 0-10%
741                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525055000000"; // 20-40%
742                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525055000000"; // 20-50%
743         } else if ( trainConfig == 105){ // variation alpha 0.75  added signal
744                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525055000000"; // 0-5%
745                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525055000000"; // 5-10%
746                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525055000000"; // 0-10%
747                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525055000000"; // 20-40%
748                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525055000000"; // 20-50%     
749         } else if ( trainConfig == 106){ // variation alpha 0.85
750                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525075000000"; // 0-5%
751                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525075000000"; // 5-10%
752                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525075000000"; // 0-10%
753                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525075000000"; // 20-40%
754                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525075000000"; // 20-50%
755         } else if ( trainConfig == 107){ // variation alpha 0.85  added signal
756                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0]= "01525075000000"; // 0-5%
757                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1]= "01525075000000"; // 5-10%
758                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2]= "01525075000000"; // 0-10%
759                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3]= "01525075000000"; // 20-40%
760                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4]= "01525075000000"; // 20-50%
761         } else if ( trainConfig == 108){ // psi pair 0.2
762                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028280400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
763                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028280400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
764                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028280400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
765                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028280400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
766                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028280400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
767         } else if ( trainConfig == 109){ // psi pair 0.2  added signal
768                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028280400000"; mesonCutArray[ 0]= "01525065000000"; // 0-5%
769                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028280400000"; mesonCutArray[ 1]= "01525065000000"; // 5-10%
770                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028280400000"; mesonCutArray[ 2]= "01525065000000"; // 0-10%
771                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028280400000"; mesonCutArray[ 3]= "01525065000000"; // 20-40%
772                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028280400000"; mesonCutArray[ 4]= "01525065000000"; // 20-50%
773         } else if ( trainConfig == 110){ // cleaner cuts
774                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
775                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
776                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
777                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
778                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
779         } else if ( trainConfig == 111){ // cleaner cuts added signal
780                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
781                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
782                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
783                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
784                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%             
785         } else if ( trainConfig == 112){ // cleaner cuts photon Quality 1
786                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
787                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
788                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
789                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
790                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
791         } else if ( trainConfig == 113){ // cleaner cuts added signal photon Quality 1
792                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
793                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
794                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
795                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
796                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%             
797         } else if ( trainConfig == 114){ // cleaner cuts photon Quality 2
798                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250430000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
799                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250430000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
800                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250430000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
801                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250430000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
802                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250430000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
803         } else if ( trainConfig == 115){ // cleaner cuts added signal photon Quality 2
804                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250430000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
805                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250430000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
806                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250430000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
807                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250430000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
808                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250430000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%             
809         } else if ( trainConfig == 116){ // cleaner cuts photon Quality 3
810                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
811                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
812                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
813                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
814                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
815         } else if ( trainConfig == 117){ // cleaner cuts added signal photon Quality 3
816                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
817                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
818                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
819                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
820                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%             
821         } else if ( trainConfig == 118){ // cleaner cuts, min R = 35 cm
822                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "007092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
823                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "007092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
824                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "007092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
825                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "007092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
826                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "007092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
827         } else if ( trainConfig == 119){ // cleaner cuts, min R = 35 cm added signal
828                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "007092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
829                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "007092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
830                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "007092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
831                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "007092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
832                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "007092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%             
833         } else if ( trainConfig == 120){ // cleaner cuts, photon Quality 1, min R = 35 cm
834                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "007092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
835                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "007092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
836                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "007092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
837                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "007092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
838                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "007092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
839         } else if ( trainConfig == 121){ // cleaner cuts added signal, photon Quality 1, min R = 35 cm
840                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "007092970028250420000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
841                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "007092970028250420000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
842                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "007092970028250420000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
843                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "007092970028250420000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
844                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "007092970028250420000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%             
845         } else if ( trainConfig == 122){ // cleaner cuts, photon Quality 3, min R = 35 cm
846                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "007092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
847                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "007092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
848                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "007092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
849                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "007092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
850                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "007092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
851         } else if ( trainConfig == 123){ // cleaner cuts added signal, photon Quality 3, min R = 35 cm
852                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "007092970028250440000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
853                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "007092970028250440000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
854                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "007092970028250440000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
855                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "007092970028250440000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
856                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "007092970028250440000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%             
857         } else if ( trainConfig == 124){ // cleaner cuts, specific centrality trigger selection 
858                 eventCutArray[ 0] = "6013301"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // kCentral 0-5%
859                 eventCutArray[ 1] = "6123301"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // kCentral 5-10%
860                 eventCutArray[ 2] = "5013301"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // kCentral 0-10%
861                 eventCutArray[ 3] = "5013601"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // kSemiCentral 0-10%
862                 eventCutArray[ 4] = "5123601"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // kSemiCentral 10-20%
863         } else if ( trainConfig == 125){ // cleaner cuts, specific centrality trigger selection 
864                 eventCutArray[ 0] = "5233601"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // kSemiCentral 20-30%
865                 eventCutArray[ 1] = "5343601"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // kSemiCentral 30-40%
866                 eventCutArray[ 2] = "5453601"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // kSemiCentral 40-50%
867                 eventCutArray[ 3] = "5243601"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // kSemiCentral 20-40%
868                 eventCutArray[ 4] = "5253601"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // kSemiCentral 20-50%
869         } else if ( trainConfig == 126){ // cleaner cuts, pion line at 2.5 sigma
870                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092370028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
871                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092370028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
872                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092370028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
873                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092370028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
874                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092370028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
875         } else if ( trainConfig == 127){ // cleaner cuts, pion line at 2.5 sigma added signals
876                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092370028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
877                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092370028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
878                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092370028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
879                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092370028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
880                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092370028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
881         } else if ( trainConfig == 128){ // cleaner cuts, pion line at 2.0 sigma
882                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
883                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
884                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
885                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
886                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092570028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
887         } else if ( trainConfig == 129){ // cleaner cuts, pion line at 2.0 sigma added signals
888                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092570028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
889                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092570028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
890                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092570028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
891                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092570028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
892                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092570028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
893         } else if ( trainConfig == 130){ // cleaner cuts, pion line at 2.0, high pt at 1.0
894                 eventCutArray[ 0] = "6010001"; photonCutArray[ 0] = "002092874028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
895                 eventCutArray[ 1] = "6120001"; photonCutArray[ 1] = "002092874028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
896                 eventCutArray[ 2] = "5010001"; photonCutArray[ 2] = "002092874028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
897                 eventCutArray[ 3] = "5240001"; photonCutArray[ 3] = "002092874028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
898                 eventCutArray[ 4] = "5250001"; photonCutArray[ 4] = "002092874028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
899         } else if ( trainConfig == 131){ // cleaner cuts, pion line at 2.0, high pt at 1.0, added signals
900                 eventCutArray[ 0] = "6010002"; photonCutArray[ 0] = "002092874028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
901                 eventCutArray[ 1] = "6120002"; photonCutArray[ 1] = "002092874028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
902                 eventCutArray[ 2] = "5010002"; photonCutArray[ 2] = "002092874028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
903                 eventCutArray[ 3] = "5240002"; photonCutArray[ 3] = "002092874028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
904                 eventCutArray[ 4] = "5250002"; photonCutArray[ 4] = "002092874028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%
905         } else if ( trainConfig == 132){ // cleaner cuts, finer centrality slices
906                 eventCutArray[ 0] = "5120001"; photonCutArray[ 0] = "002092970028250400000"; mesonCutArray[ 0] = "01525065000000"; // 0-5%
907                 eventCutArray[ 1] = "5230001"; photonCutArray[ 1] = "002092970028250400000"; mesonCutArray[ 1] = "01525065000000"; // 5-10%
908                 eventCutArray[ 2] = "5340001"; photonCutArray[ 2] = "002092970028250400000"; mesonCutArray[ 2] = "01525065000000"; // 0-10%
909                 eventCutArray[ 3] = "5450001"; photonCutArray[ 3] = "002092970028250400000"; mesonCutArray[ 3] = "01525065000000"; // 10-20%
910                 eventCutArray[ 4] = "5560001"; photonCutArray[ 4] = "002092970028250400000"; mesonCutArray[ 4] = "01525065000000"; // 0-20%     
911         } else {
912                 Error(Form("GammaConvV1_%i",trainConfig), "wrong trainConfig variable no cuts have been specified for the configuration");
913                 return;
914         }
915
916         TList *EventCutList = new TList();
917         TList *ConvCutList = new TList();
918         TList *MesonCutList = new TList();
919
920         TList *HeaderList = new TList();
921         if (periodName.CompareTo("LHC13d2")==0){
922                 TObjString *Header1 = new TObjString("pi0_1");
923                 HeaderList->Add(Header1);
924         //    TObjString *Header3 = new TObjString("eta_2");
925         //    HeaderList->Add(Header3);
926
927         } else if (periodName.CompareTo("LHC12a17x_fix")==0){
928                 TObjString *Header1 = new TObjString("PARAM");
929                 HeaderList->Add(Header1);
930         } else if (periodName.CompareTo("LHC14a1a")==0){
931                 if (headerSelectionInt == 1){ 
932                         TObjString *Header1 = new TObjString("pi0_1");
933                         HeaderList->Add(Header1);
934                 } else if (headerSelectionInt == 2){
935                         TObjString *Header1 = new TObjString("eta_2");
936                         HeaderList->Add(Header1);
937                 } else {
938                         TObjString *Header1 = new TObjString("pi0_1");
939                         HeaderList->Add(Header1);
940                         TObjString *Header2 = new TObjString("eta_2");
941                         HeaderList->Add(Header2);
942                 }  
943         } else if (periodName.CompareTo("LHC14a1b")==0 || periodName.CompareTo("LHC14a1c")==0){
944                 TObjString *Header1 = new TObjString("BOX");
945                 HeaderList->Add(Header1);
946         }       
947
948         EventCutList->SetOwner(kTRUE);
949         AliConvEventCuts **analysisEventCuts = new AliConvEventCuts*[numberOfCuts];
950         ConvCutList->SetOwner(kTRUE);
951         AliConversionPhotonCuts **analysisCuts = new AliConversionPhotonCuts*[numberOfCuts];
952         MesonCutList->SetOwner(kTRUE);
953         AliConversionMesonCuts **analysisMesonCuts = new AliConversionMesonCuts*[numberOfCuts];
954
955         for(Int_t i = 0; i<numberOfCuts; i++){
956                 analysisEventCuts[i] = new AliConvEventCuts();
957                 if (trainConfig == 1 ||trainConfig == 5 || trainConfig == 9 || trainConfig == 13 || trainConfig == 17 || trainConfig == 21 || trainConfig == 25 || trainConfig == 29 || trainConfig == 33 || trainConfig == 37){ // || trainConfig == 41 
958                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0005TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M");
959                         if (i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0510TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M");
960                         if (i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0010TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M");
961                         if (i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE,fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_1020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M");
962                         if (i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE,fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_0020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M");
963                 } else if (trainConfig == 2 ||trainConfig == 6 || trainConfig == 10 || trainConfig == 14 || trainConfig == 18 || trainConfig == 22 || trainConfig == 26 || trainConfig == 30  || trainConfig == 34 || trainConfig == 38){ // || trainConfig == 42
964                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_2040TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M");
965                         if (i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_4060TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M");
966                         if (i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_6080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M");
967                         if (i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_PbPb_2760GeV_4080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4080V0M");
968                 } else if (trainConfig == 3 ||trainConfig == 7 || trainConfig == 11 || trainConfig == 15 || trainConfig == 19 || trainConfig == 23 || trainConfig == 27 || trainConfig == 31 || trainConfig == 35 || trainConfig == 39 ){ //|| trainConfig == 43 
969                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0005TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M");
970                         if (i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0510TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M");
971                         if (i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0010TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M");
972                         if (i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE,fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_1020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M");
973                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_0020TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M");
974                 } else if (trainConfig == 4 ||trainConfig == 8 || trainConfig == 12 || trainConfig == 16 || trainConfig == 20 || trainConfig == 24 || trainConfig == 28 || trainConfig == 32 || trainConfig == 36 || trainConfig == 40){ // || trainConfig == 44 
975                         if (i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_2040TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M");
976                         if (i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_4060TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M");
977                         if (i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_6080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M");
978                         if (i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kFALSE, kFALSE, fileNameInputForWeighting, "Pi0_Hijing_LHC13d2_addSig_PbPb_2760GeV_4080TPC", "", "","Pi0_Fit_Data_PbPb_2760GeV_4080V0M");
979                 }
980
981                 if (trainConfig == 56 ){
982                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
983                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0005TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0005TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M","Eta_Fit_Data_PbPb_2760GeV_0005V0M");
984                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0510TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0510TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M","Eta_Fit_Data_PbPb_2760GeV_0510V0M");
985                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0010TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0010TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M","Eta_Fit_Data_PbPb_2760GeV_0010V0M");
986                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_1020TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_1020TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M","Eta_Fit_Data_PbPb_2760GeV_1020V0M");
987                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0020TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0020TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M","Eta_Fit_Data_PbPb_2760GeV_0020V0M");
988                         }       
989                 }         
990                 if (trainConfig == 57 ){
991                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
992                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0005TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0005TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M","Eta_Fit_Data_PbPb_2760GeV_0005V0M");
993                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0510TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0510TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M","Eta_Fit_Data_PbPb_2760GeV_0510V0M");
994                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0010TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0010TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M","Eta_Fit_Data_PbPb_2760GeV_0010V0M");
995                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_1020TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_1020TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M","Eta_Fit_Data_PbPb_2760GeV_1020V0M");
996                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0020TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0020TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0020V0M","Eta_Fit_Data_PbPb_2760GeV_0020V0M");
997                         }       
998                 }         
999                 if (trainConfig == 58 ){
1000                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
1001                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_2040TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_2040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M","Eta_Fit_Data_PbPb_2760GeV_2040V0M");
1002                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_4060TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_4060TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M","Eta_Fit_Data_PbPb_2760GeV_4060V0M");
1003                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_6080TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_6080TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M","Eta_Fit_Data_PbPb_2760GeV_6080V0M");
1004                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_4080TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_4080TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M","Eta_Fit_Data_PbPb_2760GeV_1020V0M");
1005                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_3050TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_3050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_3050V0M","Eta_Fit_Data_PbPb_2760GeV_3050V0M");
1006                         }       
1007                 }         
1008                 if (trainConfig == 59 ){
1009                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
1010                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2040TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M","Eta_Fit_Data_PbPb_2760GeV_2040V0M");
1011                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_4060TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_4060TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_4060V0M","Eta_Fit_Data_PbPb_2760GeV_4060V0M");
1012                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_6080TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_6080TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_6080V0M","Eta_Fit_Data_PbPb_2760GeV_6080V0M");
1013                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_4080TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_4080TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_1020V0M","Eta_Fit_Data_PbPb_2760GeV_1020V0M");
1014                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_3050TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_3050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_3050V0M","Eta_Fit_Data_PbPb_2760GeV_3050V0M");
1015                         }       
1016                 }         
1017         
1018                 if (trainConfig == 60 ){
1019                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
1020                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_2030TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_2030TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2030V0M","Eta_Fit_Data_PbPb_2760GeV_2030V0M");
1021                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_3040TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_3040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_3040V0M","Eta_Fit_Data_PbPb_2760GeV_3040V0M");
1022                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_4050TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_4050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_4050V0M","Eta_Fit_Data_PbPb_2760GeV_4050V0M");
1023                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_5060TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_5060TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_5060V0M","Eta_Fit_Data_PbPb_2760GeV_5060V0M");
1024                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_2050TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_2050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2050V0M","Eta_Fit_Data_PbPb_2760GeV_2050V0M");
1025                         }       
1026                 }         
1027                 if (trainConfig == 61 ){
1028                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
1029                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2030TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2030TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2030V0M","Eta_Fit_Data_PbPb_2760GeV_2030V0M");
1030                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_3040TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_3040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_3040V0M","Eta_Fit_Data_PbPb_2760GeV_3040V0M");
1031                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_4050TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_4050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_4050V0M","Eta_Fit_Data_PbPb_2760GeV_4050V0M");
1032                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_5060TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_5060TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_5060V0M","Eta_Fit_Data_PbPb_2760GeV_5060V0M");
1033                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2050TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2050V0M","Eta_Fit_Data_PbPb_2760GeV_2050V0M");
1034                         }       
1035                 }         
1036                 
1037                 if ( trainConfig == 70 || trainConfig == 72  || trainConfig == 74  || trainConfig == 76  || trainConfig == 78  || trainConfig == 80  || trainConfig == 82  || trainConfig == 84 || trainConfig == 86  || trainConfig == 88  || trainConfig == 90 || trainConfig == 92 || trainConfig == 94 || trainConfig == 96  || trainConfig == 98  || trainConfig == 100 || trainConfig == 102  || trainConfig == 104 || trainConfig == 106 || trainConfig == 108 || trainConfig == 110 || trainConfig == 112 || trainConfig == 114 || trainConfig == 116 || trainConfig == 118 || trainConfig == 120 || trainConfig == 122){
1038                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
1039                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0005TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0005TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M","Eta_Fit_Data_PbPb_2760GeV_0005V0M");
1040                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0510TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0510TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M","Eta_Fit_Data_PbPb_2760GeV_0510V0M");
1041                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_0010TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_0010TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M","Eta_Fit_Data_PbPb_2760GeV_0010V0M");
1042                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_2040TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_2040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M","Eta_Fit_Data_PbPb_2760GeV_2040V0M");
1043                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_PbPb_2760GeV_2050TPC",periodName.Data()), Form("Eta_Hijing_%s_PbPb_2760GeV_2050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2050V0M","Eta_Fit_Data_PbPb_2760GeV_2050V0M");
1044                         }       
1045                 } 
1046                 if ( trainConfig == 71 || trainConfig == 73  || trainConfig == 75  || trainConfig == 77  || trainConfig == 79  || trainConfig == 81  || trainConfig == 83  || trainConfig == 85 || trainConfig == 87  || trainConfig == 89  || trainConfig == 91 || trainConfig == 93 || trainConfig == 95 || trainConfig == 97  || trainConfig == 99  || trainConfig == 101 || trainConfig == 103  || trainConfig == 105 || trainConfig == 107 || trainConfig == 109 || trainConfig == 111 || trainConfig == 113 || trainConfig == 115 || trainConfig == 117 || trainConfig == 119 || trainConfig == 121 || trainConfig == 123){
1047                         if (periodName.CompareTo("LHC14a1a") ==0 || periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
1048                                 if ( i == 0 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0005TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0005TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0005V0M","Eta_Fit_Data_PbPb_2760GeV_0005V0M");
1049                                 if ( i == 1 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0510TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0510TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0510V0M","Eta_Fit_Data_PbPb_2760GeV_0510V0M");
1050                                 if ( i == 2 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_0010TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_0010TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_0010V0M","Eta_Fit_Data_PbPb_2760GeV_0010V0M");
1051                                 if ( i == 3 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2040TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2040TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2040V0M","Eta_Fit_Data_PbPb_2760GeV_2040V0M");
1052                                 if ( i == 4 && doWeighting)  analysisEventCuts[i]->SetUseReweightingWithHistogramFromFile(kTRUE, kTRUE, kFALSE,fileNameInputForWeighting, Form("Pi0_Hijing_%s_addSig_PbPb_2760GeV_2050TPC",periodName.Data()), Form("Eta_Hijing_%s_addSig_PbPb_2760GeV_2050TPC",periodName.Data()), "","Pi0_Fit_Data_PbPb_2760GeV_2050V0M","Eta_Fit_Data_PbPb_2760GeV_2050V0M");
1053                         }       
1054                 } 
1055                 
1056                 analysisEventCuts[i]->InitializeCutsFromCutString(eventCutArray[i].Data());
1057                 if (periodName.CompareTo("LHC14a1b") ==0 || periodName.CompareTo("LHC14a1c") ==0 ){
1058                         if (headerSelectionInt == 1) analysisEventCuts[i]->SetAddedSignalPDGCode(111);
1059                         if (headerSelectionInt == 2) analysisEventCuts[i]->SetAddedSignalPDGCode(221);
1060                 }
1061                 EventCutList->Add(analysisEventCuts[i]);                
1062                 if (trainConfig == 37 || trainConfig == 38){
1063                         analysisEventCuts[i]->SelectSpecialTrigger(AliVEvent::kMB, "AliVEvent::kMB" );
1064                 }
1065                 if (trainConfig == 39 || trainConfig == 40){   
1066                         analysisEventCuts[i]->SelectSpecialTrigger(AliVEvent::kCentral,"AliVEvent::kCentral" );
1067                 }   
1068                 if (trainConfig == 41 || trainConfig == 42){   
1069                         analysisEventCuts[i]->SelectSpecialTrigger(AliVEvent::kSemiCentral,"AliVEvent::kSemiCentral" );
1070                 }   
1071                 analysisEventCuts[i]->SetFillCutHistograms("",kFALSE);
1072                 
1073                 analysisCuts[i] = new AliConversionPhotonCuts();
1074                 analysisCuts[i]->InitializeCutsFromCutString(photonCutArray[i].Data());
1075                 ConvCutList->Add(analysisCuts[i]);
1076                 analysisCuts[i]->SetFillCutHistograms("",kFALSE);
1077
1078                 analysisMesonCuts[i] = new AliConversionMesonCuts();
1079                 analysisMesonCuts[i]->InitializeCutsFromCutString(mesonCutArray[i].Data());
1080                 MesonCutList->Add(analysisMesonCuts[i]);
1081                 analysisMesonCuts[i]->SetFillCutHistograms("");
1082                 
1083                 analysisEventCuts[i]->SetAcceptedHeader(HeaderList);
1084         }
1085         
1086         task->SetEventCutList(numberOfCuts,EventCutList);
1087         task->SetConversionCutList(numberOfCuts,ConvCutList);
1088         task->SetMesonCutList(numberOfCuts,MesonCutList);
1089         task->SetMoveParticleAccordingToVertex(kTRUE);
1090         task->SetDoMesonAnalysis(kTRUE);
1091         task->SetDoMesonQA(enableQAMesonTask); //Attention new switch for Pi0 QA
1092         task->SetDoPhotonQA(enableQAPhotonTask);  //Attention new switch small for Photon QA
1093
1094         //connect containers
1095         AliAnalysisDataContainer *coutput =
1096                 mgr->CreateContainer(Form("GammaConvV1_%i",trainConfig), TList::Class(),
1097                                                         AliAnalysisManager::kOutputContainer,Form("GammaConvV1_%i.root",trainConfig));
1098
1099         mgr->AddTask(task);
1100         mgr->ConnectInput(task,0,cinput);
1101         mgr->ConnectOutput(task,1,coutput);
1102
1103         return;
1104
1105 }