flag jets with too high p_T track
[u/mrichter/AliRoot.git] / PWG4 / macros / ConfigTrainGrid.C
CommitLineData
b6b70f4b 1{
2 /*
3
4 // < 1 are private
5 // run >= 100 < 1000 are real data
6 // runs > 1000 are sim data
7 1 = LHC09a1 pt hard 15 - 50 GeV
8 2 = LHC09a2 pt hard 50 -100 GeV
9 3 = LHC09a3 pt hard > 100
10 100 = LHC09d 0.9 GeV real with TPC
11 200 = LHC10b 7 TeV real with TPC
12 202 = LHC10b 7 TeV real with TPC (pass2)
13 300 = LHC10c 900 GeV real with TPC
14 302 = LHC10c 900 GeV pass2 real with TPC
15 400 = LHC10c 7000 GeV real with TPC
16 402 = LHC10c 7000 GeV pass2 real with TPC
17 500 = LHC10d 7000 GeV Real
18 502 = LHC10d 7000 GeV Real pass2
19 801 = LHC10h Pb+Pb real
20 802 = LHC10h Pb+Pb real pass2
21 1012 = LHC10a12 0.9 GeV MC Pythiy D6T
22 1014 = LHC10a14 0.9 GeV MC ojet D6T
23 2001 = LHC10b1 7 TeV Phojet with real OCDB
24 2002 = LHC10b2 7 TeV Pythia Pergiua-0
25 4001 = LHC10d1 pp, Pythia6 Perugia-0, 0.5T, 7000GeV, // misses some runs since i a different direcotory LHC10c9
26 4002 = LHC10d2 pp, Phojet, 0.5T, 7000GeV
27 4004 = LHC10d4 pp, Perugis, 0.5T, 7000GeV
28 5014 = LHC10e14 pp, Jet+Jet different p_T hard bins 0.5T, 7000GeV
29 6006 = LHC10f6 pp, Phojet Nachor Runs LHC10d
30 8012 = 50 80 GeV (LHC10h12)
31 8013 = 80-inf GeV (LHC10h13)
32 8102 = LHC11a2XX (where x = a - j)
33 */
34 Bool_t debugConfig =false;
35 Bool_t ckbConfig = false;
36 Bool_t productionConfig = false; // make this true for productions mode
37 Int_t iVersion = 1; // this is the train version for one datapass
38
39 // bextra == 0 4 plus
40 // bextra == 1 large pass1 split..
41 // bextra == 2 3 plus
b360e4be 42 Int_t bRun = 802; Int_t bExtra = 0; char* cDate = "110818a";
dc6bbd38 43 // Int_t bRun = 8102; Int_t bExtra = 1; char* cDate = "110725a";
b6b70f4b 44 iAODanalysis = 0;
45 // 1 == Read Jets and tracks form the input AOD
46 // needs the jet branchnames set explicitly
47 // 2 == Read tracks from input and jets from output
48
49 // 1 = write the Full AOD for all events
50 // 2 = write the Full AOD for triggered events
51 // 3 = write the deltaAOD for all events
dc6bbd38 52 // iFilterAnalysis = 2;kJetTriggerPtCut = 40;
53 iFilterAnalysis = 3;
54 // iFilterAnalysis = 2;
d3a3f33d 55
b6b70f4b 56 if (kPluginMode.Contains("merge")){
57 // currently merging this one...
dc6bbd38 58 // cDate = "110717a";
65c43de8 59 // bRun = 802; Int_t bExtra = 0; cDate = "110812a";
b6b70f4b 60 }
61 kUseDebug = kFALSE;
62 // this is for testing just one run...
63 // kGridMaxRunsFromList = 1;kUseDebug = kTRUE;
64
65 kUseSysInfo = 100;
66
67 kFillAOD = kTRUE; // fill the aod
68 kGridMergeExclude = ""; kSaveAOD = (1<<0)|(1<<1)|(1<<3); //
69 if(!productionConfig){
70 kUsePAR = kTRUE; //
71 kUseCPAR = kTRUE;
72 kGridFilesPerJob = 20;
73 }
74
75 iPhysicsSelection = 1;
76
77 if(iAODanalysis){
78 // iAODAn
79 if(iAODanalysis == 1){
80 iJETAN = 0;
81 iPWG4Cluster = 0;
82 iJETSUBTRACT = 0;
83 }
84 }
85
86 iJETAN = 3;
87 iDIJETAN = 1; // only run on one JetFinder for the moment
88 iPWG1QASym = 0; // excluded since also on QA train
89 iPWG4TmpSourceSara = 0;
90 iPWG4JetServices = 1;
91 iPWG4Fragmentation = 0;
92 iPWG4JetSpectrum = 7;
93 iPWG4UE = 0; // tmp off awating updates
94 iPWG4LeadingUE = 1;
95 iPWG4CorrectionsUE = 0; // 19.07. OFF awaiting changes by Sara
96 iPWG4PtQAMC = 1;
97 iPWG4PtSpectra = 1;
98 iPWG4PtQATPC = 3;
99 iPWG4PtTrackQA = 1;
100 iPWG4Cosmics = 0;
101 iPWG4JetChem = 1;
102 iPWG4QGSep = 1;
103 iPWG4Minijet = 1;
104 iPWG4ThreeJets = 0; // tmp off mem leak
105 iPWG4KMeans = 1; // Off no feedback
106 iPWG4Cluster = 5;
107 iPWG4Tagged = 1;
108 iPWG4PartCorr = 1;
109 iPWG4CaloQA = 1;
110 iPWG4JetCorr = 0;
111 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
112 iPWG4GammaConv = 1; // TMP OFF for merging
113 iPWG4CaloConv = 0; // 25.08. off: Output size 450MB in memory 03.09 crashes
114 iPWG4omega3pi = 1; //
115 kDeltaAODJetName = ""; // for OTF we do not need to write a delta/saves some extension gymnastics in the tasks
116 kDeltaAODPartCorrName = "";
117 kPluginExecutableCommand = "cat wn.xml; root -b -q "; // dump the file list to stdout for debugging
118 kPluginAliRootVersion = "";
119 kGridRunsPerMaster = 1; // To have the separate configs for eacj run in separate dirs
120
121 if(bRun<100){ // private MC
122 iPWG1QASym = 0;
123 iPWG4TmpSourceSara = 0;
124 iPWG4JetChem = 0;
125 iPWG4UE = 0;
126 iPWG4Cluster = 0;
127 iPWG4PtQAMC = 0;
128 iPWG4PtSpectra = 0;
129 iPWG4PtQATPC = 0;
130 iPWG4Cosmics = 0; // tmp on
131 iPWG4ThreeJets = 0;
132 iPWG4KMeans = 0;
133 iPWG4PartCorr = 0;
134 iPWG4CaloQA = 0;
135 iPWG4CaloConv = 0;
136 iPWG4Tagged = 0;
137 iPWG4JetCorr = 0;
138 iPWG4GammaConv = 0;
139 iPWG4omega3pi = 0;
140 kUseKinefilter = kTRUE;
141 kIsMC = kTRUE;
142 kUseTR = kTRUE;
143 kUseMC = kTRUE;
144 kUseAODMC = kTRUE;
145 // kErrorIgnoreLevel = 3001;
146 // kPluginFastReadOption = kTRUE;
147 }
148
149 if(bRun>=1000){
150 kUseKinefilter = kTRUE;
151 kIsMC = kTRUE;
152 kUseMC = kTRUE;
153 kUseAODMC = kTRUE;
154 kUseTR = kTRUE;
155 iPhysicsSelection = 1;
156 }
157
158
159 if (bRun == 1){
160 kGridRunRange[0] = 0;
161 kGridRunRange[1] = -1;
162 kGridDatadir = "/alice/sim/PDC_08b/LHC09a1/";
163 kGridDataSet = "LHC09a1";
164 kGridRunPattern = "%05d";
165 kGridLocalRunList = "fp_lhc09a1.txt";
166 kTrainName = Form("pwg4train_LHC09a1_%s",cDate);
167 iPhysicsSelection = 0;
168 kHighPtFilterMask = 32;
169 kGridFilesPerJob = 200;
170 }
171 else if (bRun == 2){
172 kGridRunRange[0] = 0;
173 kGridRunRange[1] = -1;
174 kGridDatadir = "/alice/sim/PDC_08b/LHC09a2/";
175 kGridDataSet = "LHC09a2";
176 kGridRunPattern = "%05d";
177 kGridLocalRunList = "fp_lhc09a2.txt";
178 kTrainName = Form("pwg4train_LHC09a2_%s",cDate);
179 iPhysicsSelection = 0;
180 kHighPtFilterMask = 32;
181 kGridFilesPerJob = 80;
182 }
183 else if (bRun == 3){
184 kGridRunRange[0] = 0;
185 kGridRunRange[1] = -1; //
186 kGridDatadir = "/alice/sim/PDC_08b/LHC09a3/";
187 kGridDataSet = "LHC09a3";
188 kGridRunPattern = "%05d";
189 kGridLocalRunList = "fp_lhc09a3.txt";
190 kTrainName = Form("pwg4train_LHC09a3_%s",cDate);
191 iPhysicsSelection = 0;
192 kHighPtFilterMask = 32;
193 kGridFilesPerJob = 80;
194 }
195 else if (bRun == 100){
196 kGridRunRange[0] = 0;
197 kGridRunRange[1] = -1; //
198 kGridDatadir = "/alice/data/2009/LHC09d/";
199 kGridDataSet = "LHC09d";
200 kGridRunPattern = "%09d";
201 kUseKinefilter = kFALSE;
202 kIsMC = kFALSE; // is real!
203 kUseMC = kFALSE;
204 kUseAODMC = kFALSE;
205 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
206 // kGridPassPattern = "pass1"; kGridLocalRunList = "fp_runlist_pass1.txt"; kTrainName = "pwg4train_LHC09d_pass1_100104";
207 kGridPassPattern = "pass6"; kGridLocalRunList = "fp_runlist_pass6.txt"; kTrainName = Form("pwg4train_LHC09d_pass6_%s",cDate);
208 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
209 }
210 else if (bRun == 200){
211 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
212 kGridRunRange[1] = -1; //
213 kGridDatadir = "/alice/data/2010/LHC10b/";
214 kGridDataSet = "LHC10b";
215 kGridRunPattern = "%09d";
216 kUseKinefilter = kFALSE;
217 kIsMC = kFALSE; // is real!
218 kUseMC = kFALSE;
219 kUseAODMC = kFALSE;
220 kGridPassPattern = "pass1"; kGridLocalRunList = "fp_lhc10b_runlist_pass1.txt"; kTrainName = Form("pwg4train_LHC10b_pass1_%s",cDate);
221 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
222 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
223 // stich of tasks not from PWG4JetTasks
224 gROOT->LoadMacro("cleanXML.C");
225 cleanXML();
226 }
227 else if (bRun == 202){
228 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
229 kGridRunRange[1] = -1; //
230 kGridDatadir = "/alice/data/2010/LHC10b/";
231 kGridDataSet = "LHC10b";
232 kGridRunPattern = "%09d";
233 kUseKinefilter = kFALSE;
234 kIsMC = kFALSE; // is real!
235 kUseMC = kFALSE;
236 kUseAODMC = kFALSE;
237 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
238 kGridPassPattern = "pass2"; kGridLocalRunList = "fp_lhc10b_runlist_pass2.txt"; kTrainName = Form("pwg4train_LHC10b_pass2_%s",cDate);
239
240 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
241 // stich of tasks not from PWG4JetTasks
242
243 // switch off tasks with no feedback...
244 iPWG4ThreeJets = 0; // tmp off mem leak
245 iPWG4KMeans = 0; // OFF no FEEDBACK
246 iPWG4Tagged = 0; // OFF crashes on MC
247 iPWG4CaloQA = 0; // OFF not needed on MC
248 iPWG4JetCorr = 0;
249 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
250 iPWG4omega3pi = 0; // OFF no feedback
251
252 gROOT->LoadMacro("cleanXML.C");
253 cleanXML();
254 }
255 else if (bRun == 300){
256 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
257 kGridRunRange[1] = -1; //
258 kGridDatadir = "/alice/data/2010/LHC10c/";
259 kGridDataSet = "LHC10c";
260 kGridRunPattern = "%09d";
261 kUseKinefilter = kFALSE;
262 kIsMC = kFALSE; // is real!
263 kUseMC = kFALSE;
264 kUseAODMC = kFALSE;
265 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
266 kGridPassPattern = "pass1"; kGridLocalRunList = "fp_lhc10c_runlist_900_pass1.txt"; kTrainName = Form("pwg4train_LHC10c_900_pass1_%s",cDate);
267
268 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
269 // stich of tasks not from PWG4JetTasks
270 gROOT->LoadMacro("cleanXML.C");
271 cleanXML();
272 }
273 else if (bRun == 302){
274 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
275 kGridRunRange[1] = -1; //
276 kGridDatadir = "/alice/data/2010/LHC10c/";
277 kGridDataSet = "LHC10c";
278 kGridRunPattern = "%09d";
279 kUseKinefilter = kFALSE;
280 kIsMC = kFALSE; // is real!
281 kUseMC = kFALSE;
282 kUseAODMC = kFALSE;
283 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
284 kGridPassPattern = "pass2"; kGridLocalRunList = "fp_lhc10c_runlist_900_pass2.txt"; kTrainName = Form("pwg4train_LHC10c_900_pass2_%s",cDate);
285
286 // switch off tasks with no feedback...
287 iPWG4ThreeJets = 0; // tmp off mem leak
288 iPWG4KMeans = 0; // OFF no FEEDBACK
289 iPWG4Tagged = 0; // OFF crashes on MC
290 iPWG4CaloQA = 0; // OFF not needed on MC
291 iPWG4JetCorr = 0;
292 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
293 iPWG4omega3pi = 0; // OFF no feedback
294
295
296 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
297 // stich of tasks not from PWG4JetTasks
298 gROOT->LoadMacro("cleanXML.C");
299 cleanXML();
300 }
301 else if (bRun == 400){
302 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
303 kGridRunRange[1] = -1; //
304 kGridDatadir = "/alice/data/2010/LHC10c/";
305 kGridDataSet = "LHC10c";
306 kGridRunPattern = "%09d";
307 kUseKinefilter = kFALSE;
308 kIsMC = kFALSE; // is real!
309 kUseMC = kFALSE;
310 kUseAODMC = kFALSE;
311 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
312 kGridPassPattern = "pass1"; kGridLocalRunList = "fp_lhc10c_runlist_7000_pass1.txt"; kTrainName = Form("pwg4train_LHC10c_7000_pass1_%s",cDate);
313
314 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
315 // stich of tasks not from PWG4JetTasks
316 gROOT->LoadMacro("cleanXML.C");
317 cleanXML();
318 }
319 else if (bRun == 402){
320 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
321 kGridRunRange[1] = -1; //
322 kGridFilesPerJob = 10;
323 kGridDatadir = "/alice/data/2010/LHC10c/";
324 kGridDataSet = "LHC10c";
325 kGridRunPattern = "%09d";
326 kUseKinefilter = kFALSE;
327 kIsMC = kFALSE; // is real!
328 kUseMC = kFALSE;
329 kUseAODMC = kFALSE;
330 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
331 kGridPassPattern = "pass2"; kGridLocalRunList = "fp_lhc10c_runlist_7000_pass2.txt"; kTrainName = Form("pwg4train_LHC10c_7000_pass2_%s",cDate);
332
333 // switch off tasks with no feedback...
334 iPWG4ThreeJets = 0; // tmp off mem leak
335 iPWG4KMeans = 0; // OFF no FEEDBACK
336 iPWG4Tagged = 0; // OFF crashes on MC
337 iPWG4CaloQA = 0; // OFF not needed on MC
338 iPWG4JetCorr = 0;
339 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
340 iPWG4omega3pi = 0; // OFF no feedback
341
342 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
343 // stich of tasks not from PWG4JetTasks
344 gROOT->LoadMacro("cleanXML.C");
345 cleanXML();
346 }
347 else if (bRun == 500){
348 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
349 kGridRunRange[1] = -1; //
350 kGridDatadir = "/alice/data/2010/LHC10d/";
351 kGridDataSet = "LHC10d";
352 kGridRunPattern = "%09d";
353 kUseKinefilter = kFALSE;
354 kIsMC = kFALSE; // is real!
355 kUseMC = kFALSE;
356 kUseAODMC = kFALSE;
357 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
358 kGridPassPattern = "pass1"; kGridLocalRunList = "fp_lhc10d_runlist_7000_pass1.txt"; kTrainName = Form("pwg4train_LHC10d_7000_pass1_%s",cDate);
359
360 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
361 // stich of tasks not from PWG4JetTasks
362 gROOT->LoadMacro("cleanXML.C");
363 cleanXML();
364 }
365 else if (bRun == 502){
366 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
367 kGridRunRange[1] = -1; //
368 kGridDatadir = "/alice/data/2010/LHC10d/";
369 kGridDataSet = "LHC10d";
370 kGridRunPattern = "%09d";
371 kUseKinefilter = kFALSE;
372 kIsMC = kFALSE; // is real!
373 kUseMC = kFALSE;
374 kUseAODMC = kFALSE;
375
376 // switch off tasks with no feedback...
377 iPWG4ThreeJets = 0; // tmp off mem leak
378 iPWG4KMeans = 0; // OFF no FEEDBACK
379 iPWG4Tagged = 0; // OFF crashes on MC
380 iPWG4CaloQA = 0; // OFF not needed on MC
381 iPWG4JetCorr = 0;
382 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
383 iPWG4omega3pi = 0; // OFF no feedback
384
385
386 // swtich off most tasks for Pb+Pb
387 iDIJETAN = 0; // only run on one JetFinder for the moment
388 iPWG4Fragmentation = 0;
389 iPWG4LeadingUE = 0;
390 iPWG4JetChem = 0;
391 iPWG4QGSep = 0;
392 iPWG4Minijet = 0;
393 iPWG4PartCorr = 0;
394 iPWG4GammaConv = 0;
395
396 // running as light a possible
397
398 iPWG4PtQAMC = 1;
399 iPWG4PtSpectra = 1;
400 iPWG4PtQATPC = 1;
401 iPWG4JetSpectrum = 1;
402 iPWG4JetServices = 1; // !!!!!!!!!!!
403 iPWG4Cluster = 1;// not 5....
404 kHighPtFilterMask = 1<<8; // 256 TPC related to SPD
405
406 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
407 kGridPassPattern = "pass2"; kGridLocalRunList = "fp_lhc10d_runlist_7000_pass2.txt"; kTrainName = Form("p4t_10d_7000_p2_%s",cDate);
408 // kDeltaAODJetName = "AliAOD.Jets.root";
409 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
410 // stich of tasks not from PWG4JetTasks
411 gROOT->LoadMacro("cleanXML.C");
412 cleanXML();
413 }
414 else if (bRun == 801){
415 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
416 kGridRunRange[1] = -1; //
417 kGridDatadir = "/alice/data/2010/LHC10h/";
418 kGridDataSet = "LHC10h";
419 kGridRunPattern = "%09d";
420 kUseKinefilter = kFALSE;
421 kIsMC = kFALSE; // is real!
422 kUseMC = kFALSE;
423 kUseAODMC = kFALSE;
424
425 // switch off tasks with no feedback...
426 iPWG4ThreeJets = 0; // tmp off mem leak
427 iPWG4KMeans = 0; // OFF no FEEDBACK
428 iPWG4Tagged = 0; // OFF crashes on MC
429 iPWG4CaloQA = 0; // OFF not needed on MC
430 iPWG4JetCorr = 0;
431 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
432 iPWG4omega3pi = 0; // OFF no feedback
433
434 // swtich off most tasks for Pb+Pb
435 iDIJETAN = 0; // only run on one JetFinder for the moment
436 iPWG4Fragmentation = 0;
437 iPWG4LeadingUE = 0;
438 iPWG4JetChem = 0;
439 iPWG4QGSep = 0;
440 iPWG4Minijet = 0;
441 iPWG4PartCorr = 0;
442 iPWG4GammaConv = 0;
443
444 // running as light a possible
445
446 iPWG4PtQAMC = 1;
447 iPWG4PtSpectra = 1;
448 iPWG4PtQATPC = 1;
449 iPWG4JetSpectrum = 1;
450 iPWG4JetServices = 1; // !!!!!!!!!!!
451 iPWG4Cluster = 1;// not 5....
452 kHighPtFilterMask = 1<<8; // 256 TPC related to SPD
453
454
455 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
456 kGridPassPattern = "pass1_4plus"; kGridLocalRunList = "fp_lhc10h_pass1_4plus.txt";
457 kTrainName = Form("t_10h_1_4p_%s",cDate);
458 if(bExtra==1){
459 kGridPassPattern = "pass1"; kGridLocalRunList = "fp_lhc10h_pass1.txt";
460 // kTrainName = Form("train_pass1_%s",cDate);
461 kTrainName = Form("t_10h_p1_%s",cDate);
462 }
463 else if(bExtra==2){
464 kGridPassPattern = "pass1_plusplusplus"; kGridLocalRunList = "fp_lhc10h_pass1_3plus.txt";
465 // kTrainName = Form("train_LHC10h_pass1_3p_%s",cDate);
466 kTrainName = Form("t_10h_p1_3p_%s",cDate);
467 }
468 kSaveAOD = 1;
469 kDeltaAODJetName = "AliAOD.Jets.root";kSaveAOD = 2;
470
471 if (kPluginMode.Contains("merge")){
472 kSaveAOD = 0; //
473 }
474 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
475 // stich of tasks not from PWG4JetTasks
476 // gROOT->LoadMacro("cleanXML.C");
477 // cleanXML();
478 }
479 else if (bRun == 802){
480 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
481 kGridRunRange[1] = -1; //
482 kGridDatadir = "/alice/data/2010/LHC10h/";
483 kGridDataSet = "LHC10h";
484 kGridExtraAliendirLevel = Form("/%s",cDate);kGridOutdir = "output";
485 kGridRunPattern = "%09d";
486 kUseKinefilter = kFALSE;
487 kIsMC = kFALSE; // is real!
488 kUseMC = kFALSE;
489 kUseAODMC = kFALSE;
490
491 if(iVersion==1){
492 // switch off tasks with no feedback...
493 iPWG4ThreeJets = 0; // tmp off mem leak
494 iPWG4KMeans = 0; // OFF no FEEDBACK
495 iPWG4Tagged = 0; // OFF crashes on MC
496 iPWG4CaloQA = 0; // OFF not needed on MC
497 iPWG4JetCorr = 0;
498 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
499 iPWG4omega3pi = 0; // OFF no feedback
500
501 // swtich off most tasks for Pb+Pb
502 iDIJETAN = 0; // only run on one JetFinder for the moment
503 iPWG4LeadingUE = 0;
504 iPWG4JetChem = 0;
505 iPWG4QGSep = 0;
506 iPWG4Minijet = 0;
507 iPWG4PartCorr = 0;
508 iPWG4GammaConv = 0;
509
510 // running as light a possible
511
512 iPWG4PtQAMC = 0; //
513 iPWG4PtQATPC = 0; //
eb1739b7 514 iPWG4PtSpectra = 0; //
b6b70f4b 515 iPWG4PtTrackQA = 1;
516 iPWG4JetSpectrum = 1;
517 iPWG4JetServices = 1; // !!!!!!!!!!!
518 iPWG4Cluster = 1;// not 5....
d3a3f33d 519 kHighPtFilterMask = 1<<4|1<<8; // Global tracks with SPD requirment global constraitn for the rest
b6b70f4b 520 iPWG4Fragmentation = 1;
521 //
522 }// version1
523
524
525 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
526 kGridPassPattern = "pass2"; // reset for AODs below
527
528 kGridLocalRunList = "fp_lhc10h_pass2c.txt";
529 if(bExtra==1)kGridLocalRunList = "fp_lhc10h_pass2.txt";
530 kTrainName = Form("p4t_10h_pass2_%s",cDate);
531 if (kPluginMode.Contains("merge")){
532 // kGridLocalRunList = "fp_lhc10h_pass2b.txt";
533 kGridLocalRunList = "out_merge_100_1.txt";
534 }
535 if(iAODanalysis)kTrainName = Form("tAOD_10h_pass2_%s",cDate);
536
537 if(iFilterAnalysis==1){
538 kSaveAOD = 1;
539 kGridFilesPerJob = 5;
540 kDeltaAODJetName = "";
541 kFilterAOD = false;
542 }
543 else if(iFilterAnalysis==2){
544 kSaveAOD = 1;
545 kGridFilesPerJob = 20;
546 kDeltaAODJetName = "";
547 kFilterAOD = true;
548 }
549 else if(iFilterAnalysis == 3){
550 kSaveAOD = 2;
551 kGridFilesPerJob = 20;
552 kDeltaAODJetName = "AliAOD.Jets.root";
553 kFilterAOD = true;
554 }
555
556 if(iAODanalysis==0){
557 // select what is running in the spectrum task, list of jet names is filles automatically
558 kGridFilesPerJob = 30;
559 /*
560 // Thnsparses...
21df3cc7 561
b6b70f4b 562 ############# Possible jet branches ###################
21df3cc7 563 1: jetsAOD_UA104_B0_Filter00016_Cut01000
564 2: jetsAOD_UA104_B0_Filter00016_Cut02000
565 3: jetsAOD_UA104_B2_Filter00016_Cut01000
566 4: jetsAOD_UA104_B2_Filter00016_Cut02000
567 5: clustersAOD_KT04_B0_Filter00016_Cut00150_Skip00
568 6: clustersAOD_KT04_B0_Filter00016_Cut00150_Skip00RandomConeSkip00
569 7: clustersAOD_KT04_B0_Filter00016_Cut00150_Skip00RandomCone_random
570 8: clustersAOD_KT04_B0_Filter00016_Cut02000_Skip00
571 9: clustersAOD_KT04_B0_Filter00016_Cut02000_Skip00RandomConeSkip00
572 10: clustersAOD_KT04_B0_Filter00016_Cut02000_Skip00RandomCone_random
573 11: clustersAOD_KT02_B0_Filter00016_Cut00150_Skip00
574 12: clustersAOD_ANTIKT04_B0_Filter00016_Cut00150_Skip02
575 13: clustersAOD_ANTIKT04_B0_Filter00016_Cut00150_Skip02RandomConeSkip02
576 14: clustersAOD_ANTIKT04_B0_Filter00016_Cut00150_Skip02RandomCone_random
577 15: clustersAOD_ANTIKT04_B0_Filter00016_Cut02000_Skip02
578 16: clustersAOD_ANTIKT02_B0_Filter00016_Cut00150_Skip00
579 17: clustersAOD_KT04_B1_Filter00016_Cut00150_Skip00RandomConeSkip00
580 18: clustersAOD_KT04_B1_Filter00016_Cut00150_Skip00RandomCone_random
581 19: clustersAOD_ANTIKT04_B1_Filter00016_Cut00150_Skip02
582 20: clustersAOD_ANTIKT04_B1_Filter00016_Cut00150_Skip02RandomConeSkip02
583 21: clustersAOD_ANTIKT04_B1_Filter00016_Cut00150_Skip02RandomCone_random
584 22: clustersAOD_ANTIKT02_B1_Filter00016_Cut00150_Skip00
585 23: clustersAOD_KT04_B2_Filter00016_Cut00150_Skip00RandomConeSkip00
586 24: clustersAOD_KT04_B2_Filter00016_Cut00150_Skip00RandomCone_random
587 25: clustersAOD_ANTIKT04_B2_Filter00016_Cut00150_Skip02
588 26: clustersAOD_ANTIKT04_B2_Filter00016_Cut00150_Skip02RandomConeSkip02
589 27: clustersAOD_ANTIKT04_B2_Filter00016_Cut00150_Skip02RandomCone_random
590 28: clustersAOD_ANTIKT02_B2_Filter00016_Cut00150_Skip00
591 29: clustersAOD_KT04_B2_Filter00016_Cut02000_Skip00RandomConeSkip00
592 30: clustersAOD_KT04_B2_Filter00016_Cut02000_Skip00RandomCone_random
593 31: clustersAOD_ANTIKT04_B2_Filter00016_Cut02000_Skip02
594 */
2eded560 595 /*
596 ############# Possible jet branches ###################
597 1: jetsAOD_UA104_B0_Filter00272_Cut01000
598 2: jetsAOD_UA104_B0_Filter00272_Cut02000
599 3: jetsAOD_UA104_B2_Filter00272_Cut01000
600 4: jetsAOD_UA104_B2_Filter00272_Cut02000
601 5: clustersAOD_KT04_B0_Filter00272_Cut00150_Skip00
602 6: clustersAOD_KT04_B0_Filter00272_Cut00150_Skip00RandomConeSkip00
603 7: clustersAOD_KT04_B0_Filter00272_Cut00150_Skip00RandomCone_random
604 8: clustersAOD_KT04_B0_Filter00272_Cut01000_Skip00
605 9: clustersAOD_KT04_B0_Filter00272_Cut01000_Skip00RandomConeSkip00
606 10: clustersAOD_KT04_B0_Filter00272_Cut01000_Skip00RandomCone_random
607 11: clustersAOD_KT04_B0_Filter00272_Cut02000_Skip00
608 12: clustersAOD_KT04_B0_Filter00272_Cut02000_Skip00RandomConeSkip00
609 13: clustersAOD_KT04_B0_Filter00272_Cut02000_Skip00RandomCone_random
610 14: clustersAOD_KT02_B0_Filter00272_Cut00150_Skip00
611 15: clustersAOD_ANTIKT04_B0_Filter00272_Cut00150_Skip02
612 16: clustersAOD_ANTIKT04_B0_Filter00272_Cut00150_Skip02RandomConeSkip02
613 17: clustersAOD_ANTIKT04_B0_Filter00272_Cut00150_Skip02RandomCone_random
614 18: clustersAOD_ANTIKT04_B0_Filter00272_Cut01000_Skip02
615 19: clustersAOD_ANTIKT04_B0_Filter00272_Cut02000_Skip02
616 20: clustersAOD_ANTIKT02_B0_Filter00272_Cut00150_Skip00
617 21: clustersAOD_KT04_B1_Filter00272_Cut00150_Skip00RandomConeSkip00
618 22: clustersAOD_KT04_B1_Filter00272_Cut00150_Skip00RandomCone_random
619 23: clustersAOD_ANTIKT04_B1_Filter00272_Cut00150_Skip02
620 24: clustersAOD_ANTIKT04_B1_Filter00272_Cut00150_Skip02RandomConeSkip02
621 25: clustersAOD_ANTIKT04_B1_Filter00272_Cut00150_Skip02RandomCone_random
622 26: clustersAOD_ANTIKT02_B1_Filter00272_Cut00150_Skip00
623 27: clustersAOD_KT04_B2_Filter00272_Cut00150_Skip00RandomConeSkip00
624 28: clustersAOD_KT04_B2_Filter00272_Cut00150_Skip00RandomCone_random
625 29: clustersAOD_ANTIKT04_B2_Filter00272_Cut00150_Skip02
626 30: clustersAOD_ANTIKT04_B2_Filter00272_Cut00150_Skip02RandomConeSkip02
627 31: clustersAOD_ANTIKT04_B2_Filter00272_Cut00150_Skip02RandomCone_random
628 32: clustersAOD_ANTIKT02_B2_Filter00272_Cut00150_Skip00
629 33: clustersAOD_KT04_B2_Filter00272_Cut01000_Skip00RandomConeSkip00
630 34: clustersAOD_KT04_B2_Filter00272_Cut01000_Skip00RandomCone_random
631 35: clustersAOD_ANTIKT04_B2_Filter00272_Cut01000_Skip02
632 36: clustersAOD_KT04_B2_Filter00272_Cut02000_Skip00RandomConeSkip00
633 37: clustersAOD_KT04_B2_Filter00272_Cut02000_Skip00RandomCone_random
634 38: clustersAOD_ANTIKT04_B2_Filter00272_Cut02000_Skip02
635 */
636
637 // in the first map we fill the correlations we want to plot
638 // in the jet back map we associated the branche used for background calculation
639 // to fetch the multiplicity
640
641 // UA1
b6b70f4b 642 kJetMapSpectrum.Add(4,2);
643 kJetBackMapSpectrum.Add(4,8);
644
645 // anti kT 150 MeV
2eded560 646 kJetMapSpectrum.Add(29,15);
647 kJetBackMapSpectrum.Add(29,5);
648 kJetBackMapSpectrum.Add(15,5);
649
650 // anti kT 1000 MeV
651 kJetMapSpectrum.Add(35,18);
652 kJetBackMapSpectrum.Add(35,8);
653 kJetBackMapSpectrum.Add(35,8);
b6b70f4b 654
655 // anti kT 2000 MeV
2eded560 656 kJetMapSpectrum.Add(38,19);
657 kJetBackMapSpectrum.Add(35,11);
658 kJetBackMapSpectrum.Add(19,11);
659
b6b70f4b 660
661 // anti kT 0.2
2eded560 662 kJetMapSpectrum.Add(32,14);
663 kJetBackMapSpectrum.Add(32,5);
664 kJetBackMapSpectrum.Add(14,5);
b6b70f4b 665
666 // random cones
667 kJetMapSpectrum.Add(6,7);
668 kJetBackMapSpectrum.Add(6,5);
669 kJetBackMapSpectrum.Add(7,5);
670
671 kJetMapSpectrum.Add(9,10);
672 kJetBackMapSpectrum.Add(9,8);
673 kJetBackMapSpectrum.Add(10,8);
674
675 }
676 else if (iAODanalysis){
677 kGridFilesPerJob = 20;
678 if(iAODanalysis==1){
679 kGridPassPattern = "pass2/AOD049";
680 iJETAN = 0;
681 iPWG4Cluster = 0;
682 iJETSUBTRACT = 0;
683
684
685 /*
686 reading from AOD043 input
687 1 jetsAOD_UA104_B0_Filter00128_Cut01000 1 0x102c068d0
688 2 jetsAOD_UA104_B2_Filter00128_Cut01000 1 0x102c068d0
689 3 jetsAOD_SISCONE04_B0_Filter00128_Cut00150 1 0x102c068d0
690 4 clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00 1 0x102c068d0
691 5 clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00_random 1 0x102c068d0
692 6 clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00RandomConeSkip00 1 0x102c068d0
693 7 clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00RandomCone_random 1 0x102c068d0
694 8 clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02 1 0x102c068d0
695 9 clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02_random 1 0x102c068d0
696 10 clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02RandomConeSkip02 1 0x102c068d0
697 11 clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02RandomCone_random 1 0x102c068d0
698 12 clustersAOD_ANTIKT02_B0_Filter00128_Cut00150_Skip00 1 0x102c068d0
699 13 clustersAOD_ANTIKT02_B0_Filter00128_Cut00150_Skip00_random 1 0x102c068d0
700 14 jetsAOD_SISCONE04_B1_Filter00128_Cut00150 1 0x102c068d0
701 15 clustersAOD_ANTIKT04_B1_Filter00128_Cut00150_Skip02 1 0x102c068d0
702 16 clustersAOD_ANTIKT02_B1_Filter00128_Cut00150_Skip00 1 0x102c068d0
703 */
704
705 kJetListSpectrum.Add(new TObjString("jetsAOD_UA104_B0_Filter00128_Cut01000"));
706 kJetListSpectrum.Add(new TObjString("jetsAOD_UA104_B2_Filter00128_Cut01000"));
707 kJetListSpectrum.Add(new TObjString("jetsAOD_SISCONE04_B0_Filter00128_Cut00150"));
708 kJetListSpectrum.Add(new TObjString("clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00"));
709 kJetListSpectrum.Add(new TObjString("clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00_random"));
710 kJetListSpectrum.Add(new TObjString("clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00RandomConeSkip00"));
711 kJetListSpectrum.Add(new TObjString("clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00RandomCone_random"));
712 kJetListSpectrum.Add(new TObjString("clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02"));
713 kJetListSpectrum.Add(new TObjString("clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02_random"));
714 kJetListSpectrum.Add(new TObjString("clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02RandomConeSkip02"));
715 kJetListSpectrum.Add(new TObjString("clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02RandomCone_random"));
716 kJetListSpectrum.Add(new TObjString("clustersAOD_ANTIKT02_B0_Filter00128_Cut00150_Skip00"));
717 kJetListSpectrum.Add(new TObjString("clustersAOD_ANTIKT02_B0_Filter00128_Cut00150_Skip00_random"));
718 kJetListSpectrum.Add(new TObjString("jetsAOD_SISCONE04_B1_Filter00128_Cut00150"));
719 kJetListSpectrum.Add(new TObjString("clustersAOD_ANTIKT04_B1_Filter00128_Cut00150_Skip02"));
720 kJetListSpectrum.Add(new TObjString("clustersAOD_ANTIKT02_B1_Filter00128_Cut00150_Skip00"));
721
722 // anti kT 150 MeV
723 kJetMapSpectrum.Add(15,8);
724 kJetBackMapSpectrum.Add(8,4);
725 kJetBackMapSpectrum.Add(15,4);
726
727 // anti kT R = 0.2
728 kJetMapSpectrum.Add(17,12);
729 kJetBackMapSpectrum.Add(17,4);
730 kJetBackMapSpectrum.Add(12,4);
731
732 // random cones
733 kJetMapSpectrum.Add(6,7);
734 kJetBackMapSpectrum.Add(6,4);
735 kJetBackMapSpectrum.Add(7,4);
736
737 kJetMapSpectrum.Add(10,11);
738 kJetBackMapSpectrum.Add(10,4);
739 kJetBackMapSpectrum.Add(11,4);
740 kDeltaAODJetNameInput = "AliAOD.Jets.root";
741 kDeltaAODJetName = "";
742 kHighPtFilterMask = 128; // centrally produced AOD
743 }
744 else if(iAODanalysis==2){
745 /*
746 ############# Possible jet branches ###################
747 1: jetsAOD_UA104_B0_Filter00128_Cut01000
748 2: jetsAOD_UA104_B0_Filter00128_Cut02000
749 3: jetsAOD_UA104_B2_Filter00128_Cut01000
750 4: jetsAOD_UA104_B2_Filter00128_Cut02000
751 5: clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00
752 6: clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00RandomConeSkip00
753 7: clustersAOD_KT04_B0_Filter00128_Cut00150_Skip00RandomCone_random
754 8: clustersAOD_KT04_B0_Filter00128_Cut02000_Skip00
755 9: clustersAOD_KT04_B0_Filter00128_Cut02000_Skip00RandomConeSkip00
756 10: clustersAOD_KT04_B0_Filter00128_Cut02000_Skip00RandomCone_random
757 11: clustersAOD_KT02_B0_Filter00128_Cut00150_Skip00
758 12: clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02
759 13: clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02RandomConeSkip02
760 14: clustersAOD_ANTIKT04_B0_Filter00128_Cut00150_Skip02RandomCone_random
761 15: clustersAOD_ANTIKT04_B0_Filter00128_Cut02000_Skip02
762 16: clustersAOD_ANTIKT02_B0_Filter00128_Cut00150_Skip00
763 17: clustersAOD_KT04_B1_Filter00128_Cut00150_Skip00Skip00RandomConeSkip00
764 18: clustersAOD_KT04_B1_Filter00128_Cut00150_Skip00Skip00RandomCone_random
765 19: clustersAOD_ANTIKT04_B1_Filter00128_Cut00150_Skip02
766 20: clustersAOD_ANTIKT04_B1_Filter00128_Cut00150_Skip02Skip02RandomConeSkip02
767 21: clustersAOD_ANTIKT04_B1_Filter00128_Cut00150_Skip02Skip02RandomCone_random
768 22: clustersAOD_ANTIKT02_B1_Filter00128_Cut00150_Skip00
769 23: clustersAOD_KT04_B2_Filter00128_Cut00150_Skip00Skip00RandomConeSkip00
770 24: clustersAOD_KT04_B2_Filter00128_Cut00150_Skip00Skip00RandomCone_random
771 25: clustersAOD_ANTIKT04_B2_Filter00128_Cut00150_Skip02
772 26: clustersAOD_ANTIKT04_B2_Filter00128_Cut00150_Skip02Skip02RandomConeSkip02
773 27: clustersAOD_ANTIKT04_B2_Filter00128_Cut00150_Skip02Skip02RandomCone_random
774 28: clustersAOD_ANTIKT02_B2_Filter00128_Cut00150_Skip00
775 29: clustersAOD_KT04_B2_Filter00128_Cut02000_Skip00Skip00RandomConeSkip00
776 30: clustersAOD_KT04_B2_Filter00128_Cut02000_Skip00Skip00RandomCone_random
777 31: clustersAOD_ANTIKT04_B2_Filter00128_Cut02000_Skip02
778 */
779
780 // anti kT 150 MeV
781 kJetMapSpectrum.Add(25,12);
782 kJetBackMapSpectrum.Add(25,5);
783 kJetBackMapSpectrum.Add(12,5);
784
785 // anti kT R = 0.2
786 kJetMapSpectrum.Add(28,16);
787 kJetBackMapSpectrum.Add(28,5);
788 kJetBackMapSpectrum.Add(16,5);
789
790 // anti kT 2 GeV
791 kJetMapSpectrum.Add(31,15);
792 kJetBackMapSpectrum.Add(31,8);
793 kJetBackMapSpectrum.Add(15,8);
794
795
796 // random cones
797 kJetMapSpectrum.Add(9,10);
798 kJetBackMapSpectrum.Add(10,5);
799 kJetBackMapSpectrum.Add(9,5);
800
801 kJetMapSpectrum.Add(13,14);
802 kJetBackMapSpectrum.Add(14,5);
803 kJetBackMapSpectrum.Add(13,5);
804 kDeltaAODJetNameInput = "";
805
806
807 kDeltaAODJetName = "AliAOD.Jets.root";kSaveAOD = 2;
808 kHighPtFilterMask = 128; // centrally produced AOD
809 iPWG4Fragmentation = 0; // off for a pass
810
811
812 }
813 }
814 if (kPluginMode.Contains("merge")){
815 kSaveAOD = 0; //
816 }
817 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
818 // stich of tasks not from PWG4JetTasks
819 // gROOT->LoadMacro("cleanXML.C");
820 // cleanXML();
821 }
822 else if (bRun == 1012){
823 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
824 kGridRunRange[1] = -1; //
825 kGridDatadir = "/alice/sim/LHC10a12/";
826 kGridDataSet = "LHC10a12";
827 kGridRunPattern = "%06d";
828 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
829 kGridPassPattern = ""; kGridLocalRunList = "fp_lhc10a12.txt"; kTrainName = Form("pwg4train_LHC10a12_%s",cDate);
830 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
831 kGridFilesPerJob = 250; // only few events in a sim file
832 // stich of tasks not from PWG4JetTasks
833 gROOT->LoadMacro("cleanXML.C");
834 cleanXML();
835 }
836 else if (bRun == 1014){
837 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
838 kGridRunRange[1] = -1; //
839 kGridDatadir = "/alice/sim/LHC10a14/";
840 kGridDataSet = "LHC10a14";
841 kGridRunPattern = "%06d";
842 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
843 kGridPassPattern = ""; kGridLocalRunList = "fp_lhc10a14.txt"; kTrainName = Form("pwg4train_LHC10a14_%s",cDate);
844 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
845 kGridFilesPerJob = 250; // only few events in a sim file
846 // stich of tasks not from PWG4JetTasks
847 gROOT->LoadMacro("cleanXML.C");
848 cleanXML();
849 }
850 else if (bRun == 2001){
851 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
852 kGridRunRange[1] = -1; //
853 kGridDatadir = "/alice/sim/LHC10b1/";
854 kGridDataSet = "LHC10b1";
855 kGridRunPattern = "%06d";
856 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
857 kGridPassPattern = ""; kGridLocalRunList = "fp_lhc10b1.txt"; kTrainName = Form("pwg4train_LHC10b1_%s",cDate);
858 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
859 kGridFilesPerJob = 250; // only few events in a sim file
860 // stich of tasks not from PWG4JetTasks
861 gROOT->LoadMacro("cleanXML.C");
862 cleanXML();
863 }
864 else if (bRun == 2002){
865 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
866 kGridRunRange[1] = -1; //
867 kGridDatadir = "/alice/sim/LHC10b2/";
868 kGridDataSet = "LHC10b2";
869 kGridRunPattern = "%06d";
870 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
871 kGridPassPattern = ""; kGridLocalRunList = "fp_lhc10b2.txt"; kTrainName = Form("pwg4train_LHC10b2_%s",cDate);
872 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
873 kGridFilesPerJob = 250; // only few events in a sim file
874 // stich of tasks not from PWG4JetTasks
875 gROOT->LoadMacro("cleanXML.C");
876 cleanXML();
877 }
878 else if (bRun == 4001){
879 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
880 kGridRunRange[1] = -1; //
881 kGridDatadir = "/alice/sim/LHC10d1/";
882 kGridDataSet = "LHC10d1";
883 kGridRunPattern = "%06d";
884 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
885 kGridPassPattern = ""; kGridLocalRunList = "fp_lhc10d1.txt"; kTrainName = Form("pwg4train_LHC10d1_%s",cDate);
886 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
887 kGridFilesPerJob = 50; // only few events in a sim file
888 // stich of tasks not from PWG4JetTasks
889 gROOT->LoadMacro("cleanXML.C");
890 cleanXML();
891 }
892 else if (bRun == 4002){
893 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
894 kGridRunRange[1] = -1; //
895 kGridDatadir = "/alice/sim/LHC10d2/";
896 kGridDataSet = "LHC10d2";
897 kGridRunPattern = "%06d";
898 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
899 kGridPassPattern = ""; kGridLocalRunList = "fp_lhc10d2.txt"; kTrainName = Form("pwg4train_LHC10d2_%s",cDate);
900 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
901 kGridFilesPerJob = 100; // only few events in a sim file
902 // stich of tasks not from PWG4JetTasks
903 gROOT->LoadMacro("cleanXML.C");
904 cleanXML();
905 }
906 else if (bRun == 4004){
907 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
908 kGridRunRange[1] = -1; //
909 kGridDatadir = "/alice/sim/LHC10d4/";
910 kGridDataSet = "LHC10d4";
911 kGridRunPattern = "%06d";
912 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
913 kGridPassPattern = ""; kGridLocalRunList = "fp_lhc10d4.txt"; kTrainName = Form("pwg4train_LHC10d4_%s",cDate);
914 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
915 kGridFilesPerJob = 100; // only few events in a sim file
916 // stich of tasks not from PWG4JetTasks
917 iPWG4TmpSourceSara = 0;
918 iPWG4UE = 0; //
919 iPWG4ThreeJets = 0; // tmp off mem leak
920 iPWG4KMeans = 0; // OFF no FEEDBACK
921 iPWG4Tagged = 0; // OFF crashes on MC
922 iPWG4CaloQA = 0; // OFF not needed on MC
923 iPWG4JetCorr = 0;
924 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
925 iPWG4CaloConv = 0; // 25.08. off: Output size 03.09 crashes
926 iPWG4omega3pi = 0; // OFF no feedback
927 gROOT->LoadMacro("cleanXML.C");
928 cleanXML();
929 }
930 else if(bRun==5014){
931 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
932 kGridRunRange[1] = -1; //
933 kGridDatadir = "/alice/sim/LHC10e14/120823/";
934 kGridRunPattern = "%d";
935 kGridDataSet = "LHC10e14";
936 kGridLocalRunList = "fp_lhc10e14.txt"; kTrainName = Form("pwg4train_LHC10e14_%s",cDate);
937 kGridPassPattern = "";
938
939 iPWG4TmpSourceSara = 0;
940 iPWG4UE = 0; // off not needed on Jet+Jet
941 iPWG4LeadingUE = 0; // off not needed on Jet +Jet
942 iPWG4CorrectionsUE = 0; // off not needed on Jet +Jet
943 iPWG4Cosmics = 0; // off not needed on Jet +Jet MC
944 iPWG4JetChem = 0; // OFF no FEEDBACK
945 iPWG4ThreeJets = 0; // tmp off mem leak
946 iPWG4KMeans = 0; // OFF no FEEDBACK
947 iPWG4Tagged = 0; // OFF crashes on MC
948 iPWG4CaloQA = 0; // OFF not needed on MC
949 iPWG4JetCorr = 0;
950 iPWG4GammaConv = 0; // TMP OFF cuts not updated not so important for jet+jet
951 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
952 iPWG4CaloConv = 0; // 25.08. off: Output size 03.09 crashes
953 iPWG4omega3pi = 0; // OFF no feedback
954 iPWG1QASym = 0; // excluded since also on QA train
955
956 kDeltaAODJetName = ""; // for OTF we do not need to write a delta/saves some extension gymnastics in the tasks
957
958 kUseKinefilter = kTRUE;
959 kIsMC = kTRUE;
960 kUseTR = kTRUE;
961 kUseMC = kTRUE;
962 kUseAODMC = kTRUE;
963 // gSystem->Exec("cp wn_lhc10b.xml wn.xml");
964 kGridFilesPerJob = 50; // only few events in a sim file
965 gROOT->LoadMacro("cleanXML.C");
966 cleanXML();
967 }
968 else if (bRun == 6006){
969 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
970 kGridRunRange[1] = -1; //
971 kGridDatadir = "/alice/sim/LHC10f6/";
972 kGridDataSet = "LHC10f6";
973 kGridRunPattern = "%06d";
974 kUseKinefilter = kTRUE;
975 kIsMC = kTRUE;
976 kUseMC = kTRUE;
977 kUseAODMC = kTRUE;
978 kGridFilesPerJob = 100;
979 // switch off tasks with no feedback...
980 iPWG4ThreeJets = 0; // tmp off mem leak
981 iPWG4KMeans = 0; // OFF no FEEDBACK
982 iPWG4Tagged = 0; // OFF crashes on MC
983 iPWG4CaloQA = 0; // OFF not needed on MC
984 iPWG4JetCorr = 0;
985 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
986 iPWG4omega3pi = 0; // OFF no feedback
987 iPWG4TmpSourceSara = 0;
988 iPWG4UE = 0; //
989 iPWG4ThreeJets = 0; // tmp off mem leak
990 iPWG4KMeans = 0; // OFF no FEEDBACK
991 iPWG4Tagged = 0; // OFF crashes on MC
992 iPWG4CaloQA = 0; // OFF not needed on MC
993 iPWG4JetCorr = 0;
994 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
995 iPWG4CaloConv = 0; // 25.08. off: Output size 03.09 crashes
996 iPWG4PartCorr = 0; // OFF cjecked back with Gustavo
997 iPWG4omega3pi = 0; // OFF no feedback
998 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
999 kGridPassPattern = ""; kGridLocalRunList = "fp_lhc10f6.txt"; kTrainName = Form("pwg4train_LHC10f6_%s",cDate);
1000 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
1001 // stich of tasks not from PWG4JetTasks
1002 // gROOT->LoadMacro("cleanXML.C");
1003 // cleanXML();
1004 }
1005 else if (bRun == 8012){
1006 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
1007 kGridRunRange[1] = -1; //
1008 kGridDatadir = "/alice/sim/LHC10h12/";
1009 kGridDataSet = "LHC10h12";
1010 kGridRunPattern = "%06d";
1011 kUseKinefilter = kTRUE;
1012 kIsMC = kTRUE;
1013 kUseMC = kTRUE;
1014 kUseAODMC = kTRUE;
1015 kGridFilesPerJob = 50;
1016 // switch off tasks with no feedback...
1017
1018
1019 // switch off tasks with no feedback...
1020 iPWG4ThreeJets = 0; // tmp off mem leak
1021 iPWG4KMeans = 0; // OFF no FEEDBACK
1022 iPWG4Tagged = 0; // OFF crashes on MC
1023 iPWG4CaloQA = 0; // OFF not needed on MC
1024 iPWG4JetCorr = 0;
1025 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
1026 iPWG4omega3pi = 0; // OFF no feedback
1027
1028
1029 // swtich off most tasks for Pb+Pb
1030 iDIJETAN = 0; // only run on one JetFinder for the moment
1031 iPWG4Fragmentation = 0;
1032 iPWG4LeadingUE = 0;
1033 iPWG4JetChem = 0;
1034 iPWG4QGSep = 0;
1035 iPWG4Minijet = 0;
1036 iPWG4PartCorr = 0;
1037 iPWG4GammaConv = 0;
1038
1039 // running as light a possible
1040
1041 iPWG4PtQAMC = 0;
1042 iPWG4PtSpectra = 0;
1043 iPWG4PtQATPC = 0;
1044 iPWG4JetSpectrum = 0;
1045 iPWG4JetServices = 0; // !!!!!!!!!!!
1046 iPWG4Cluster = 1;// not 5....
1047 kHighPtFilterMask = 256;
1048
1049 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
1050 kGridPassPattern = ""; kGridLocalRunList = "fp_lhc10h12.txt"; kTrainName = Form("pwg4train_LHC10h12_%s",cDate);
1051 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
1052 // stich of tasks not from PWG4JetTasks
1053 gROOT->LoadMacro("cleanXML.C");
1054 cleanXML();
1055 }
1056 else if (bRun == 8102){
1057 kGridRunRange[0] = 0; // 0 is a bad number :(, needs a manual commit in alien...
1058 kGridRunRange[1] = -1; //
1059
1060
1061 char a;
1062 if(bExtra<10)a = bExtra+97;
1063
1064 kGridDatadir = Form("/alice/sim/LHC11a2%c/",a);
1065 kGridDataSet = Form("LHC11a2%c/",a);
1066 kGridRunPattern = "%06d";
1067 kUseKinefilter = kTRUE;
1068 kIsMC = kTRUE;
1069 kUseMC = kTRUE;
1070 kUseAODMC = kTRUE;
1071 kGridFilesPerJob = 50;
1072 // switch off tasks with no feedback...
1073
1074
1075 // switch off tasks with no feedback...
1076 iPWG4ThreeJets = 0; // tmp off mem leak
1077 iPWG4KMeans = 0; // OFF no FEEDBACK
1078 iPWG4Tagged = 0; // OFF crashes on MC
1079 iPWG4CaloQA = 0; // OFF not needed on MC
1080 iPWG4JetCorr = 0;
1081 iPWG4JCORRAN = 0; // TMP OFF: Writes a Tree
1082 iPWG4omega3pi = 0; // OFF no feedback
1083
1084
1085 // swtich off most tasks for Pb+Pb
1086 iDIJETAN = 0; // only run on one JetFinder for the moment
1087 iPWG4Fragmentation = 0;
1088 iPWG4LeadingUE = 0;
1089 iPWG4JetChem = 0;
1090 iPWG4QGSep = 0;
1091 iPWG4Minijet = 0;
1092 iPWG4PartCorr = 0;
1093 iPWG4GammaConv = 0;
1094
1095 // running as light a possible
dc6bbd38 1096 iJETSUBTRACT = 0; // no subtraction
1097
b6b70f4b 1098
1099 iPWG4PtQAMC = 1;
1100 iPWG4PtSpectra = 1;
1101 iPWG4PtQATPC = 1;
1102 iPWG4JetSpectrum = 1;
1103 iPWG4JetServices = 1; // !!!!!!!!!!!
1104 iPWG4Cluster = 1;// not 5....
dc6bbd38 1105 kHighPtFilterMask = 1<<4|1<<8;
1106
b6b70f4b 1107
dc6bbd38 1108
1109 if(iFilterAnalysis==1){
1110 kSaveAOD = 1;
1111 kDeltaAODJetName = "";
1112 kFilterAOD = false;
1113 }
1114 else if(iFilterAnalysis==2){
1115 kJetTriggerPtCut = 20; //pt
1116 kSaveAOD = 1;
1117 kDeltaAODJetName = "";
1118 kFilterAOD = true;
1119 }
1120 else if(iFilterAnalysis == 3){
1121 kSaveAOD = 2;
1122 kDeltaAODJetName = "AliAOD.Jets.root";
1123 kFilterAOD = true;
1124 }
1125 kGridFilesPerJob = 100;
1126 /*
1127############# Possible jet branches ###################
1128 1: jetsAOD_UA104_B0_Filter00272_Cut01000
1129 2: jetsAOD_UA104_B0_Filter00272_Cut02000
1130 3: jetsAODMC_UA104_B0_Filter00272_Cut01000
1131 4: jetsAODMC2_UA104_B0_Filter00272_Cut01000
1132 5: clustersAOD_KT06_B0_Filter00272_Cut00150_Skip00
1133 6: clustersAOD_KT04_B0_Filter00272_Cut00150_Skip00
1134 7: clustersAOD_ANTIKT04_B0_Filter00272_Cut00150_Skip02
1135 8: clustersAOD_ANTIKT04_B0_Filter00272_Cut00150_Skip02RandomConeSkip02
1136 9: clustersAOD_ANTIKT04_B0_Filter00272_Cut00150_Skip02RandomCone_random
1137 10: clustersAOD_ANTIKT04_B0_Filter00272_Cut02000_Skip02
1138 11: clustersAOD_ANTIKT02_B0_Filter00272_Cut00150_Skip00
1139 12: clustersAODMC_KT06_B0_Filter00272_Cut00150_Skip00
1140 13: clustersAODMC2_KT06_B0_Filter00272_Cut00150_Skip00
1141 14: clustersAODMC_ANTIKT04_B0_Filter00272_Cut00150_Skip00
1142 15: clustersAODMC2_ANTIKT04_B0_Filter00272_Cut00150_Skip00
1143*/
b6b70f4b 1144 // CLEAN XML FILES LOCALLY AND ON ALIEN WHEN STARTING A NEW PASS!
1145 kGridPassPattern = ""; kGridLocalRunList = "fp_runlist_lhc11a2.txt"; kTrainName = Form("pwg4train_LHC11a2%c_%s",a,cDate);
1146 kGridRunsPerMaster = 1; // Physcicsselection does not support more than on run per job
1147 kGridFilesPerJob = 100;
1148 // stich of tasks not from PWG4JetTasks
1149 gROOT->LoadMacro("cleanXML.C");
1150 cleanXML();
1151 }
1152
1153
1154
1155 if(ckbConfig){
1156 // switch off all but my analyses
1157 iPWG4KMeans = 0;
1158 iPWG1QASym = 0;
1159 iPWG4TmpSourceSara = 0;
1160 iPWG4UE = 0;
1161 iPWG4PtQAMC = 0;
1162 iPWG4PtSpectra = 0;
1163 iPWG4PtQATPC = 0;
1164 iPWG4ThreeJets = 0;
1165 iPWG4KMeans = 0;
1166 iPWG4Tagged = 0;
1167 iPWG4PartCorr = 0;
1168 iPWG4CaloQA = 0;
1169 iPWG4CaloConv = 0;
1170 iPWG4JetCorr = 0;
1171 iPWG4GammaConv = 0;
1172 iPWG4JetChem = 0; // tmp on
1173 iPWG4omega3pi = 0;
1174 kDeltaAODJetName = ""; // for OTF we do not need to write a delta/saves some extension gymnastics in the tasks
1175 kDeltaAODPartCorrName = "";
1176 kUseDebug = kFALSE;
1177 kPluginAliRootVersion = "";
1178 kGridFilesPerJob = 60;
1179 kTrainName = Form("%s_ckb",kTrainName.Data());
1180 }
1181
1182 if(debugConfig){
1183 // debug mode
1184 // kUsePAR = kFALSE; // cannot patch the macro for local test, need to laod FASTjet libs before loading FASTJETA.so
1185 // kUseCPAR = kFALSE;
1186 kUseSysInfo = 100;
1187 kUseDebug = kTRUE;
1188 kGridLocalRunList = "fp_runlist_pass4_debug.txt";
1189 kTrainName = Form("pwg4train_LHC09d_debug",cDate);
1190
1191 // kPluginExecutableCommand = "root -b -q ";
1192 kPluginExecutableCommand = "cat wn.xml; echo \"Root.Stacktrace: yes\" > .rootrc; root -b -q ";
1193 // kPluginExecutableCommand = "echo \"run -b -q pwg4train_LHC09d_debug.C\" > gdb.cmd; echo where >> gdb.cmd; echo quit >> gdb.cmd; echo y >> gdb.cmd; gdb -x gdb.cmd -batch root.exe";
1194 }
1195
1196 if (kPluginMode.Contains("test")){
d3a3f33d 1197 kJetTriggerPtCut = 0.01;
b6b70f4b 1198 if(kAnalysisMode.Contains("grid")){
1199 // kPluginExecutableCommand = "source ~/setup_root.txt; alienroot -b -q";
1200 kPluginExecutableCommand = "root -b -q";
1201 // kPluginExecutableCommand = "valgrind --tool=memcheck --error-limit=no --max-stackframe=3060888 --suppressions=$ROOTSYS/etc/valgrind-root.supp --leak-check=full --num-callers=15 --log-file=valgrind_memcheck.log root.exe -b -q";
1202 kPluginExecutableCommand = "export ALICE_ROOT=./ROOTFILES/;" + kPluginExecutableCommand;
1203 kUseSysInfo = 1;
1204 kUseDebug = kTRUE;
1205 if(bRun==802){
1206 kGridLocalRunList = "fp_lhc10h_anchor.txt";
1207 }
1208 kTrainName = Form("pwg4train_test");
1209 if(iAODanalysis)kNumberOfEvents = 500;
1210 }
1211 else{
1212 // local
1213 if(iAODanalysis)kNumberOfEvents = 2000;
1214 kUseSysInfo = 1;
1215 kUseDebug = kTRUE;
1216 kTrainName = Form("pwg4train_test_local");
1217 kLocalDataList = "local_esd_lhc10d_pass2.txt";
1218 kUsePAR = kFALSE; // cannot patch the macro for local test, need to laod FASTjet libs before loading FASTJETA.so
1219 kUseCPAR = kFALSE;
1220 if(bRun==802){
1221 kLocalDataList = "local_esd_lhc10h.txt";
1222 if(iAODanalysis) kLocalDataList = "local_aod_lhc10h.txt";
1223 }
1224 // iPWG4PtTrackQA = 0;
1225 // iPWG4PtQAMC = 0;
1226 // iPWG4PtSpectra = 0;
1227 // iPWG4PtQATPC = 0;
1228 // iPWG4PtTrackQA = 0;
1229 // iPWG4Cluster = 0;
1230 kUseCPAR = kFALSE;
1231 kUsePAR = kFALSE;
1232 // kNumberOfEvents = 200;
1233 // all OFF
1234 }
1235 }
1236 if(kPluginAliRootVersion.Length()==0){
1237 kPluginExecutableCommand = "export ALICE_ROOT=./ROOTFILES/;" + kPluginExecutableCommand;
1238 }
1239
1240 if (kPluginMode.Contains("merge")){
1241 // swtich of task where macros changed in the meantime
1242 }
1243
1244
1245}