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