]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGHF/hfe/macros/configs/pPb/ConfigEMCalHFEpA.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Merge branch 'feature-movesplit'
[u/mrichter/AliRoot.git] / PWGHF / hfe / macros / configs / pPb / ConfigEMCalHFEpA.C
CommitLineData
c852fdae 1///*******************************************************
2///Config Description
2f235a8d 3/// 14 July 2014
c852fdae 4///*******************************************************
5
6AliAnalysisTaskEMCalHFEpA* ConfigEMCalHFEpA(
a89620c8 7
8
c852fdae 9
10Bool_t isMC=kFALSE,
11Int_t triggerIndex=0,
12Int_t configIndex=0,
13Int_t centralityIndex=0,
14Bool_t isAOD = kFALSE,
15Bool_t isEMCal = kFALSE,
1dd85087 16Bool_t isTrigger = kFALSE,
c852fdae 17Int_t EMCalThreshould = 0 //0 == EG1, 1 == EG2
18)
19
20{
21///_______________________________________________________________________________________________________________
22///Track selection: Cuts used to ensure a minimum quality level of the tracks selected to perform the analysis
23 AliHFEcuts *hfecuts = new AliHFEcuts("hfeCutsMinBias","HFE Cuts");
24 hfecuts->CreateStandardCuts();
25
26 //TPC Cuts
27 hfecuts->SetTPCmodes(AliHFEextraCuts::kFound, AliHFEextraCuts::kFoundOverFindable);
89483c26 28 if(configIndex==1) hfecuts->SetMinNClustersTPC(90); //Minimum number of clusters on TPC
29 else if(configIndex==2) hfecuts->SetMinNClustersTPC(110);
30 else if(configIndex==3) hfecuts->SetMinNClustersTPC(80);
31 else if(configIndex==4) hfecuts->SetMinNClustersTPC(85);
32 else if(configIndex==5) hfecuts->SetMinNClustersTPC(115);
40e62b79 33 else if(configIndex==6) hfecuts->SetMinNClustersTPC(120);
34 else if(configIndex==80) hfecuts->SetMinNClustersTPC(130);
35 else if(configIndex==81) hfecuts->SetMinNClustersTPC(140);//Minimum number of clusters on TPC
c852fdae 36 else hfecuts->SetMinNClustersTPC(100); //Minimum number of clusters on TPC
37
89483c26 38 if(configIndex==7) hfecuts->SetMinNClustersTPCPID(70);
39 else if (configIndex==8) hfecuts->SetMinNClustersTPCPID(90);
31185d97 40 else if (configIndex==9) hfecuts->SetMinNClustersTPCPID(60);
89483c26 41 else if (configIndex==10) hfecuts->SetMinNClustersTPCPID(65);
42 else if (configIndex==11) hfecuts->SetMinNClustersTPCPID(100);
43 else if (configIndex==12) hfecuts->SetMinNClustersTPCPID(95);
44 else hfecuts->SetMinNClustersTPCPID(80); //Minimum number of clusters for dE/dx
a6f21076 45
c852fdae 46 hfecuts->SetMinRatioTPCclusters(0.6); //Number of clusters (Found/Findable)
47
48 //ITS
31185d97 49 if(configIndex==13) hfecuts->SetCutITSpixel(AliHFEextraCuts::kBoth); //Require at least one cluster on SPD
40e62b79 50 else if(configIndex==82) hfecuts->SetCutITSpixel(AliHFEextraCuts::kFirst);
a6f21076 51 else hfecuts->SetCutITSpixel(AliHFEextraCuts::kAny); //Require at least one cluster on SPD
c852fdae 52 //hfecuts->SetCutITSdrift(AliHFEextraCuts::kAny); //Require at least one cluster on SDD
53 hfecuts->SetCheckITSLayerStatus(kFALSE);
54
89483c26 55 if(configIndex==14) hfecuts->SetMinNClustersITS(2); //Minimum number of clusters on ITS
56 else if(configIndex==15) hfecuts->SetMinNClustersITS(4);
57 else if(configIndex==16) hfecuts->SetMinNClustersITS(1);
58 else if(configIndex==17) hfecuts->SetMinNClustersITS(5);
c852fdae 59 else hfecuts->SetMinNClustersITS(3); //Minimum number of clusters on ITS
60
61 //Additional Cuts
62 hfecuts->SetPtRange(2, 1e6); //Transversal momentum range in GeV/c
628d0da7 63
3d31001a 64 //testing this line for the DCA cut
65 //hfecuts->SetRequireDCAToVertex();
1dd85087 66 //DCA cut included in the analysis 12 March 2014
a4e971f7 67 if(configIndex==18) hfecuts->SetMaxImpactParam(2,3); //changed z to 3
68 else if(configIndex==19) hfecuts->SetMaxImpactParam(0.5,1);//r,z default on AOD 2.4, 3.2
40e62b79 69 else if(configIndex==83) hfecuts->SetMaxImpactParam(0.1,0.2);
628d0da7 70 else hfecuts->SetMaxImpactParam(1,2); //DCA to vertex
c852fdae 71
72 //Event Selection
73 hfecuts->SetVertexRange(10.); //
74 //hfecuts->SetProductionVertex(0,0.3,0,0.3); //
75///_______________________________________________________________________________________________________________
8079a103 76 // new cuts for event selection
77
78 //hfecuts->SetUseCorrelationVertex();
79 //hfecuts->SetSPDVtxResolutionCut();
80 //hfecuts->SetpApileupCut();
c852fdae 81
82///_________________________________________________________________________________________________________________________
83///Task config
84 AliAnalysisTaskEMCalHFEpA *task = new AliAnalysisTaskEMCalHFEpA(Form("HFECuts%d_%d_%d",triggerIndex,configIndex,centralityIndex));
85 printf("task ------------------------ %p\n ", task);
86 task->SetHFECuts(hfecuts);
a89620c8 87 task->SetCorrelationAnalysis(kFALSE);
c852fdae 88 task->SetAODanalysis(isAOD);
89 task->SetEventMixing(kTRUE);
90
31185d97 91 //to separate trigger threshold
c852fdae 92 if(EMCalThreshould==0 && triggerIndex==2) task->SetEMCalTriggerEG1();
93 if(EMCalThreshould==1 && triggerIndex==2) task->SetEMCalTriggerEG2();
94
95 if(isEMCal) task->SetUseEMCal();
1dd85087 96 //Bool_t isTrigger = kFALSE;
97 if(isTrigger) task->SetUseTrigger();
c852fdae 98
3db2410d 99 //task->SetUseTender();
100
101
a4e971f7 102 if(configIndex==104) task->SetdcaCut(2,3);//r,z
103 else if(configIndex==105) task->SetdcaCut(1.5,2.5);//r,z
104 else if(configIndex==106) task->SetdcaCut(0.5,1);//r,z
105 else if(configIndex==107) task->SetdcaCut(0.5,1.5);//r,z
106 else if(configIndex==108) task->SetdcaCut(0.1,0.2);//r,z
107 //else task->SetdcaCut(1,2);//r,z
a89620c8 108
1dd85087 109
89483c26 110 if(configIndex==100){
a89620c8 111 task->SetUseShowerShapeCut(kTRUE);
112 //task->SetM02Cut(0.0,0.3);
113 task->SetM20Cut(0.0,0.3);
114 }
115 task->SetBackground(kTRUE);
116
89483c26 117 //nonHFE cuts
118 if(configIndex==20) task->SetNonHFEmassCut(0.05);
628d0da7 119 else if(configIndex==21) task->SetNonHFEmassCut(0.10);
89483c26 120 else if(configIndex==22) task->SetNonHFEmassCut(0.03);
121 else if(configIndex==23) task->SetNonHFEmassCut(0.18);
122 else if(configIndex==24) task->SetNonHFEmassCut(0.01);
123 else if(configIndex==25) task->SetNonHFEmassCut(0.2);
628d0da7 124 else task->SetNonHFEmassCut(0.15);
c852fdae 125
89483c26 126 if(configIndex==26) task->SetNonHFEangleCut(0.1);
127 if(configIndex==27) task->SetNonHFEangleCut(0.15);
128 if(configIndex==28) task->SetNonHFEangleCut(0.05);
129
2c8e1f6c 130 //partner cuts
131
628d0da7 132 if(configIndex==29) task->SetAdditionalCuts(0.3,80);
2f235a8d 133
2f235a8d 134 else if(configIndex==85) task->SetAdditionalCuts(0.2,80);
135 else if(configIndex==86) task->SetAdditionalCuts(0.4,80);
3d31001a 136 else if(configIndex==87) task->SetAdditionalCuts(0.6,80);
40e62b79 137 else if(configIndex==30) task->SetAdditionalCuts(0.5,80);
138 else if(configIndex==31) task->SetAdditionalCuts(0.7,80);
139 else if(configIndex==32) task->SetAdditionalCuts(0.9,80);
3d31001a 140 else if(configIndex==84) task->SetAdditionalCuts(0.1,80);
40e62b79 141 else if(configIndex==33) task->SetAdditionalCuts(0,60);
142 else if(configIndex==34) task->SetAdditionalCuts(0,70);
143 else if(configIndex==35) task->SetAdditionalCuts(0,90);
144 else if(configIndex==36) task->SetAdditionalCuts(0,100);
2c8e1f6c 145
40e62b79 146 else task->SetAdditionalCuts(0,80);
2c8e1f6c 147
148
89483c26 149 //eta cuts
1dd85087 150 if(configIndex==40) task->SetEtaCut(-0.6,0);
151 else if (configIndex==41) task->SetEtaCut(-0.5,0.1);
152 else if (configIndex==42) task->SetEtaCut(0,0.6);
153 else if (configIndex==43) task->SetEtaCut(-0.1,0.5);
154 else if (configIndex==44) task->SetEtaCut(-0.5,0.5);
155 else if (configIndex==45) task->SetEtaCut(-0.4,0.4);
156 else if (configIndex==46) task->SetEtaCut(-0.3,0.3);
8079a103 157 else task->SetEtaCut(-0.6,0.6);
c852fdae 158
89483c26 159 //track matching cuts
160 if(configIndex==50) task->SetdPhidEtaCut(0.02,0.02);
161 else if (configIndex==51) task->SetdPhidEtaCut(0.03,0.03);
162 else if (configIndex==52) task->SetdPhidEtaCut(0.04,0.04);
a6f21076 163 else task->SetdPhidEtaCut(0.05,0.05);
164
89483c26 165 //E/p Cuts
166 if (configIndex==60) task->SetEoverPCut(0.85,1.2);
a4e971f7 167 else if (configIndex==61) task->SetEoverPCut(0.70,1.2);
168 else if (configIndex==62) task->SetEoverPCut(0.75,1.2);
169 else if (configIndex==63) task->SetEoverPCut(0.76,1.2);
170 else if (configIndex==64) task->SetEoverPCut(0.77,1.2);
171 else if (configIndex==65) task->SetEoverPCut(0.78,1.2);
172 else if (configIndex==66) task->SetEoverPCut(0.79,1.2);
173
174 else if (configIndex==67) task->SetEoverPCut(0.81,1.2);
175 else if (configIndex==68) task->SetEoverPCut(0.82,1.2);
176 else if (configIndex==69) task->SetEoverPCut(0.83,1.2);
177
178 else if (configIndex==88) task->SetEoverPCut(0.84,1.2);
179 else if (configIndex==89) task->SetEoverPCut(0.86,1.2);
180
181 else if (configIndex==90) task->SetEoverPCut(0.85,1.3);
182 else if (configIndex==91) task->SetEoverPCut(0.80,1.3);
183 else if (configIndex==92) task->SetEoverPCut(0.81,1.3);
184 else if (configIndex==93) task->SetEoverPCut(0.78,1.3);
185
186
187 else task->SetEoverPCut(0.80,1.2);
3db2410d 188
189
190 //this line is to set the change on the E/p cut used in the efficiency calculations.
191 task->SetEoverPnsigma(kTRUE);
89483c26 192
c852fdae 193
194 if(centralityIndex==0) task->SetCentrality(0,20);
195 if(centralityIndex==1) task->SetCentrality(20,40);
196 if(centralityIndex==2) task->SetCentrality(40,60);
197 if(centralityIndex==3) task->SetCentrality(60,80);
198 if(centralityIndex==4) task->SetCentrality(80,100);
199 if(centralityIndex==5) task->SetCentrality(0,100);
200///_______________________________________________________________________________________________________________
201
202///_______________________________________________________________________________________________________________
203///Particle identification
204 AliHFEpid *pid = task->GetPID();
205
206//______________________________________
207//In the case of a simulation
208 if(isMC)
209 {
210 pid->SetHasMCData(kTRUE);
211 task->SetMCanalysis();
212 }
213//______________________________________
214
215//______________________________________________________
216//Configure PID
217 //_________________________
218 //TPC PID
219 pid->AddDetector("TPC", 1); //Add TPC PID
220
221 //_________________________
222 //Configure TPC cut
223 //Defaul = -1 to 3 sigmas
224 //Note that it is also possible to define a model instead of a constant
225 //--------->For this change the "cut model"
226
227 Double_t params[4];
228 char *cutmodel;
229 cutmodel = "pol0";
230
89483c26 231 if(configIndex==70) params[0] = -1.5;
232 else if (configIndex==71) params[0] = -0.5;
1dd85087 233 else if (configIndex==72) params[0] = 0;
234 else if (configIndex==73) params[0] = 0.25;
89483c26 235 else if (configIndex==74) params[0] = -1.75;
c852fdae 236 else params[0] = -1;
237
31185d97 238 if(configIndex==75)Double_t max=1.5;
239 else if(configIndex==76)Double_t max=2.0;
240 else if(configIndex==77)Double_t max=2.5;
241 else if(configIndex==78)Double_t max=3.5;
242 else if(configIndex==79)Double_t max=4.0;
243 else Double_t max=3.0;
244
3d31001a 245
246
247 //testing hadron contamination in E/p cut
248 if(configIndex==94){
249 params[0] = 0;
250 task->SetEoverPCut(0.76,1.2);
251 }
252 if(configIndex==95){
253 params[0] = 0;
254 task->SetEoverPCut(0.78,1.2);
255 }
256 if(configIndex==96){
257 params[0] = 0;
258 task->SetEoverPCut(0.80,1.2);
259 }
260 if(configIndex==97){
261 params[0] = 0;
262 task->SetEoverPCut(0.82,1.2);
263 }
264 if(configIndex==98){
265 params[0] = 0;
266 task->SetEoverPCut(0.83,1.2);
267 }
268 if(configIndex==99){
269 params[0] = 0;
270 task->SetEoverPCut(0.84,1.2);
271 }
272 if(configIndex==100){
273 params[0] = 0;
274 task->SetEoverPCut(0.85,1.2);
275 }
276 if(configIndex==101){
277 params[0] = 0;
278 task->SetEoverPCut(0.86,1.2);
279 }
280 if(configIndex==102){
281 params[0] = 0;
6412f1b2 282 task->SetEoverPCut(0.85,1.3);
3d31001a 283 }
284 if(configIndex==103){
285 params[0] = 0;
6412f1b2 286 task->SetEoverPCut(0.86,1.3);
3d31001a 287 }
288
6412f1b2 289 //testing hadron contamination with other TPCsignal : -0.5
290
291 if(configIndex==109){
292 params[0] = -0.5;
293 task->SetEoverPCut(0.76,1.2);
294 }
295 if(configIndex==110){
296 params[0] = -0.5;
297 task->SetEoverPCut(0.78,1.2);
298 }
299 if(configIndex==111){
300 params[0] = -0.5;
301 task->SetEoverPCut(0.80,1.2);
302 }
303 if(configIndex==112){
304 params[0] = -0.5;
305 task->SetEoverPCut(0.82,1.2);
306 }
307 if(configIndex==113){
308 params[0] = -0.5;
309 task->SetEoverPCut(0.83,1.2);
310 }
311 if(configIndex==114){
312 params[0] = -0.5;
313 task->SetEoverPCut(0.84,1.2);
314 }
315 if(configIndex==115){
316 params[0] = -0.5;
317 task->SetEoverPCut(0.85,1.2);
318 }
319 if(configIndex==116){
320 params[0] = -0.5;
321 task->SetEoverPCut(0.86,1.2);
322 }
323 if(configIndex==117){
324 params[0] = -0.5;
325 task->SetEoverPCut(0.85,1.3);
326 }
327 if(configIndex==118){
328 params[0] = -0.5;
329 task->SetEoverPCut(0.86,1.3);
330 }
331
332 //shower shape cut for E/p default
333 if(configIndex==119){
334 task->SetUseShowerShapeCut(kTRUE);
335 //task->SetM02Cut(0.0,0.3);
336 task->SetM20Cut(0.0,0.5);
337 }
338
339
340
341
3d31001a 342
31185d97 343 pid->ConfigureTPCdefaultCut(cutmodel,params,max);
c852fdae 344//_______________________________________________________
31185d97 345
346///_______________________________________________________________________________________________________________
6501a5cf 347/// New configurations for random cuts -- March, 05, 2014 -- Values in the macro "Random_configurations.C"
40e62b79 348 /*
31185d97 349 if (configIndex==80){
350 hfecuts->SetMinNClustersTPC(86);
351 hfecuts->SetMinNClustersTPCPID(76);
6501a5cf 352 hfecuts->SetMinNClustersITS(3);
353 task->SetNonHFEmassCut(0.087);
31185d97 354 task->SetNonHFEangleCut(0.069);
355 task->SetAdditionalCuts(0.152, 91);
356 task->SetdPhidEtaCut(0.019, 0.044);
6501a5cf 357 task->SetEoverPCut(0.798, 1.225);
358 Double_t params[0]=-0.61;
31185d97 359 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
360 }
361 if (configIndex==81){
362 hfecuts->SetMinNClustersTPC(95);
363 hfecuts->SetMinNClustersTPCPID(99);
6501a5cf 364 hfecuts->SetMinNClustersITS(3);
365 task->SetNonHFEmassCut(0.057);
31185d97 366 task->SetNonHFEangleCut(0.135);
367 task->SetAdditionalCuts(0.351, 61);
368 task->SetdPhidEtaCut(0.012, 0.044);
6501a5cf 369 task->SetEoverPCut(0.792, 1.235);
370 Double_t params[0]=-1.17;
31185d97 371 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
372 }
373 if (configIndex==82){
374 hfecuts->SetMinNClustersTPC(117);
375 hfecuts->SetMinNClustersTPCPID(69);
376 hfecuts->SetMinNClustersITS(2);
6501a5cf 377 task->SetNonHFEmassCut(0.054);
31185d97 378 task->SetNonHFEangleCut(0.062);
379 task->SetAdditionalCuts(0.842, 91);
380 task->SetdPhidEtaCut(0.018, 0.033);
6501a5cf 381 task->SetEoverPCut(0.818, 1.212);
382 Double_t params[0]=-1.15;
31185d97 383 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
384 }
385 if (configIndex==83){
386 hfecuts->SetMinNClustersTPC(98);
387 hfecuts->SetMinNClustersTPCPID(93);
6501a5cf 388 hfecuts->SetMinNClustersITS(3);
389 task->SetNonHFEmassCut(0.083);
31185d97 390 task->SetNonHFEangleCut(0.051);
391 task->SetAdditionalCuts(0.415, 83);
392 task->SetdPhidEtaCut(0.047, 0.016);
6501a5cf 393 task->SetEoverPCut(0.826, 1.225);
394 Double_t params[0]=-1.06;
31185d97 395 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
396 }
397 if (configIndex==84){
398 hfecuts->SetMinNClustersTPC(99);
399 hfecuts->SetMinNClustersTPCPID(99);
6501a5cf 400 hfecuts->SetMinNClustersITS(3);
401 task->SetNonHFEmassCut(0.058);
31185d97 402 task->SetNonHFEangleCut(0.145);
403 task->SetAdditionalCuts(0.654, 99);
404 task->SetdPhidEtaCut(0.025, 0.014);
6501a5cf 405 task->SetEoverPCut(0.757, 1.228);
406 Double_t params[0]=-1.29;
31185d97 407 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
408 }
409 if (configIndex==85){
410 hfecuts->SetMinNClustersTPC(85);
411 hfecuts->SetMinNClustersTPCPID(91);
412 hfecuts->SetMinNClustersITS(2);
6501a5cf 413 task->SetNonHFEmassCut(0.167);
31185d97 414 task->SetNonHFEangleCut(0.144);
415 task->SetAdditionalCuts(0.897, 78);
416 task->SetdPhidEtaCut(0.046, 0.043);
6501a5cf 417 task->SetEoverPCut(0.771, 1.238);
418 Double_t params[0]=-1.16;
31185d97 419 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
420 }
421 if (configIndex==86){
422 hfecuts->SetMinNClustersTPC(104);
423 hfecuts->SetMinNClustersTPCPID(75);
424 hfecuts->SetMinNClustersITS(2);
6501a5cf 425 task->SetNonHFEmassCut(0.078);
31185d97 426 task->SetNonHFEangleCut(0.112);
427 task->SetAdditionalCuts(0.036, 93);
428 task->SetdPhidEtaCut(0.019, 0.013);
6501a5cf 429 task->SetEoverPCut(0.824, 1.211);
430 Double_t params[0]=-1.24;
31185d97 431 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
432 }
433 if (configIndex==87){
434 hfecuts->SetMinNClustersTPC(108);
435 hfecuts->SetMinNClustersTPCPID(93);
6501a5cf 436 hfecuts->SetMinNClustersITS(2);
437 task->SetNonHFEmassCut(0.128);
31185d97 438 task->SetNonHFEangleCut(0.140);
439 task->SetAdditionalCuts(0.814, 89);
440 task->SetdPhidEtaCut(0.041, 0.022);
6501a5cf 441 task->SetEoverPCut(0.762, 1.205);
442 Double_t params[0]=-0.93;
31185d97 443 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
444 }
445 if (configIndex==88){
446 hfecuts->SetMinNClustersTPC(80);
447 hfecuts->SetMinNClustersTPCPID(82);
448 hfecuts->SetMinNClustersITS(2);
6501a5cf 449 task->SetNonHFEmassCut(0.064);
31185d97 450 task->SetNonHFEangleCut(0.102);
451 task->SetAdditionalCuts(0.092, 97);
452 task->SetdPhidEtaCut(0.014, 0.031);
6501a5cf 453 task->SetEoverPCut(0.784, 1.216);
454 Double_t params[0]=-1.10;
31185d97 455 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
456 }
457 if (configIndex==89){
458 hfecuts->SetMinNClustersTPC(100);
459 hfecuts->SetMinNClustersTPCPID(66);
460 hfecuts->SetMinNClustersITS(2);
6501a5cf 461 task->SetNonHFEmassCut(0.100);
31185d97 462 task->SetNonHFEangleCut(0.082);
463 task->SetAdditionalCuts(0.339, 76);
464 task->SetdPhidEtaCut(0.040, 0.011);
6501a5cf 465 task->SetEoverPCut(0.811, 1.223);
466 Double_t params[0]=-0.61;
31185d97 467 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
468 }
469 if (configIndex==90){
470 hfecuts->SetMinNClustersTPC(106);
471 hfecuts->SetMinNClustersTPCPID(90);
6501a5cf 472 hfecuts->SetMinNClustersITS(3);
473 task->SetNonHFEmassCut(0.175);
31185d97 474 task->SetNonHFEangleCut(0.098);
475 task->SetAdditionalCuts(0.630, 91);
476 task->SetdPhidEtaCut(0.034, 0.026);
6501a5cf 477 task->SetEoverPCut(0.771, 1.249);
478 Double_t params[0]=-0.64;
31185d97 479 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
480 }
481 if (configIndex==91){
482 hfecuts->SetMinNClustersTPC(84);
483 hfecuts->SetMinNClustersTPCPID(80);
6501a5cf 484 hfecuts->SetMinNClustersITS(2);
485 task->SetNonHFEmassCut(0.083);
31185d97 486 task->SetNonHFEangleCut(0.098);
487 task->SetAdditionalCuts(0.594, 66);
488 task->SetdPhidEtaCut(0.013, 0.043);
6501a5cf 489 task->SetEoverPCut(0.809, 1.246);
490 Double_t params[0]=-0.57;
31185d97 491 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
492 }
493 if (configIndex==92){
494 hfecuts->SetMinNClustersTPC(112);
495 hfecuts->SetMinNClustersTPCPID(80);
6501a5cf 496 hfecuts->SetMinNClustersITS(3);
497 task->SetNonHFEmassCut(0.090);
31185d97 498 task->SetNonHFEangleCut(0.069);
499 task->SetAdditionalCuts(0.803, 91);
500 task->SetdPhidEtaCut(0.043, 0.031);
6501a5cf 501 task->SetEoverPCut(0.805, 1.227);
502 Double_t params[0]=-1.29;
31185d97 503 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
504 }
505 if (configIndex==93){
506 hfecuts->SetMinNClustersTPC(115);
507 hfecuts->SetMinNClustersTPCPID(82);
6501a5cf 508 hfecuts->SetMinNClustersITS(3);
509 task->SetNonHFEmassCut(0.075);
31185d97 510 task->SetNonHFEangleCut(0.113);
511 task->SetAdditionalCuts(0.637, 86);
512 task->SetdPhidEtaCut(0.043, 0.032);
6501a5cf 513 task->SetEoverPCut(0.760, 1.204);
514 Double_t params[0]=-0.84;
31185d97 515 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
516 }
517 if (configIndex==94){
518 hfecuts->SetMinNClustersTPC(81);
519 hfecuts->SetMinNClustersTPCPID(63);
6501a5cf 520 hfecuts->SetMinNClustersITS(3);
521 task->SetNonHFEmassCut(0.105);
31185d97 522 task->SetNonHFEangleCut(0.050);
523 task->SetAdditionalCuts(0.529, 87);
524 task->SetdPhidEtaCut(0.026, 0.047);
6501a5cf 525 task->SetEoverPCut(0.789, 1.200);
526 Double_t params[0]=-1.15;
31185d97 527 pid->ConfigureTPCdefaultCut(cutmodel,params,3.0);
528 }
40e62b79 529 */
31185d97 530
6501a5cf 531
c852fdae 532///_______________________________________________________________________________________________________________
533
534 printf("*************************************\n");
535 printf("Configuring standard Task:\n");
536 pid->PrintStatus();
537 printf("*************************************\n");
538
539 return task;
540}