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