1 /* $Id: CreateCuts.C,v 1.5 2008/01/11 08:28:52 jgrosseo Exp $ */
3 // this macro creates the track and event cuts used in this analysis
5 // last modified: 2012-02-06
7 // added cut modes 2000-2222 for cut studies
10 AliESDtrackCuts* CreatedNdPtTrackCuts(Int_t cutMode=1, Bool_t fieldOn = kTRUE, Bool_t hists = kTRUE)
12 AliESDtrackCuts* esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts");
15 esdTrackCuts->DefineHistograms(1);
17 Double_t cov1, cov2, cov3, cov4, cov5;
19 Double_t maxDCAtoVertex, maxDCAtoVertexXY, maxDCAtoVertexZ;
20 Double_t minNClustersTPC;
21 Double_t maxChi2PerClusterTPC;
22 Double_t minPt, maxPt;
24 // default cuts for ITS+TPC
34 maxChi2PerClusterTPC = 3.5;
36 esdTrackCuts->SetMaxCovDiagonalElements(cov1, cov2, cov3, cov4, cov5);
37 esdTrackCuts->SetMinNsigmaToVertex(nSigma);
38 esdTrackCuts->SetRequireSigmaToVertex(kTRUE);
39 esdTrackCuts->SetRequireTPCRefit(kTRUE);
40 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
41 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
42 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
44 TString tag("Global tracking");
47 // TPC-only cuts (vertex n sigma cut)
50 // beta cuts (still under investigation)
60 maxChi2PerClusterTPC = 3.5;
62 esdTrackCuts->SetMaxCovDiagonalElements(cov1, cov2, cov3, cov4, cov5);
63 esdTrackCuts->SetMinNsigmaToVertex(nSigma);
64 esdTrackCuts->SetRequireSigmaToVertex(kTRUE);
65 esdTrackCuts->SetRequireTPCRefit(kFALSE);
66 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
67 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
68 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
70 TString tag = "TPC-only tracking";
73 // TPC-only cuts (vertex maxDCAtoVertex cut)
76 // beta cuts (still under investigation)
77 maxDCAtoVertex = 3.0; // cm
79 maxChi2PerClusterTPC = 3.5;
81 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
82 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertex);
83 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertex);
84 esdTrackCuts->SetRequireTPCRefit(kFALSE);
85 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
86 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
87 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
89 TString tag = "TPC-only tracking";
92 // TPC-only no vertex cuts
95 // beta cuts (still under investigation)
97 maxChi2PerClusterTPC = 3.5;
99 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
100 esdTrackCuts->SetRequireTPCRefit(kFALSE);
101 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
102 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
103 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
105 TString tag = "TPC-only tracking";
108 // TPC-only no cuts at all
112 // beta cuts (still under investigation)
113 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
114 esdTrackCuts->SetRequireTPCRefit(kFALSE);
115 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
117 TString tag = "TPC-only tracking";
120 // TPC-only no kink removal no chi2
123 // beta cuts (still under investigation)
124 minNClustersTPC = 50;
125 //maxChi2PerClusterTPC = 3.5;
127 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
128 esdTrackCuts->SetRequireTPCRefit(kFALSE);
129 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
130 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
131 //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
133 TString tag = "TPC-only tracking";
136 // TPC-only no kink removal
139 // beta cuts (still under investigation)
140 minNClustersTPC = 50;
141 maxChi2PerClusterTPC = 3.5;
143 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
144 esdTrackCuts->SetRequireTPCRefit(kFALSE);
145 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
146 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
147 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
149 TString tag = "TPC-only tracking";
152 // TPC-only no kink removal no minNClustersTPC
155 // beta cuts (still under investigation)
156 //minNClustersTPC = 50;
157 maxChi2PerClusterTPC = 3.5;
159 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
160 esdTrackCuts->SetRequireTPCRefit(kFALSE);
161 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
162 //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
163 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
165 TString tag = "TPC-only tracking";
167 // TPC-only no kink removal no minNClustersTPC
170 // beta cuts (still under investigation)
171 //minNClustersTPC = 50;
172 maxChi2PerClusterTPC = 3.5;
173 maxDCAtoVertex = 3.0; // cm
175 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
176 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertex);
177 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertex);
178 esdTrackCuts->SetRequireTPCRefit(kFALSE);
179 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
180 //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
181 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
183 TString tag = "TPC-only tracking";
186 // TPC-only no kink removal no minNClustersTPC no maxChi2PerClusterTPC
189 // beta cuts (still under investigation)
190 //minNClustersTPC = 50;
191 //maxChi2PerClusterTPC = 3.5;
192 maxDCAtoVertex = 3.0; // cm
194 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
195 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertex);
196 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertex);
197 esdTrackCuts->SetRequireTPCRefit(kFALSE);
198 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
199 //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
200 //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
202 TString tag = "TPC-only tracking";
205 // TPC-only (loose cuts, absolute DCA cut)
208 // beta cuts (still under investigation)
209 minNClustersTPC = 50;
210 maxChi2PerClusterTPC = 4.0;
211 maxDCAtoVertex = 2.8; // cm
215 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
216 esdTrackCuts->SetRequireTPCRefit(kFALSE);
217 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
218 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
219 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
220 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertex);
221 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertex);
222 esdTrackCuts->SetPtRange(minPt,maxPt);
224 TString tag = "TPC-only tracking";
228 // TPC-only (loose cuts, no DCA cut)
231 // beta cuts (still under investigation)
232 minNClustersTPC = 50;
233 maxChi2PerClusterTPC = 4.0;
234 maxDCAtoVertexXY = 1.e10; // cm
235 maxDCAtoVertexZ = 1.e10; // cm
239 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
240 esdTrackCuts->SetRequireTPCRefit(kFALSE);
241 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
242 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
243 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
244 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
245 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
246 esdTrackCuts->SetDCAToVertex2D(kTRUE);
247 esdTrackCuts->SetPtRange(minPt,maxPt);
249 TString tag = "TPC-only tracking";
252 // TPC-only (standard cuts, no DCA cut)
255 // beta cuts (still under investigation)
256 minNClustersTPC = 96;
257 maxChi2PerClusterTPC = 3.5;
258 maxDCAtoVertexXY = 1.e10; // cm
259 maxDCAtoVertexZ = 1.e10; // cm
263 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
264 esdTrackCuts->SetRequireTPCRefit(kFALSE);
265 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
266 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
267 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
268 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
269 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
270 esdTrackCuts->SetDCAToVertex2D(kTRUE);
271 esdTrackCuts->SetPtRange(minPt,maxPt);
273 TString tag = "TPC-only tracking";
276 // TPC-only (tight cuts, no DCA cut)
279 // beta cuts (still under investigation)
280 minNClustersTPC = 120;
281 maxChi2PerClusterTPC = 3.5;
282 maxDCAtoVertexXY = 1.e10; // cm
283 maxDCAtoVertexZ = 1.e10; // cm
287 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
288 esdTrackCuts->SetRequireTPCRefit(kFALSE);
289 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
290 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
291 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
292 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
293 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
294 esdTrackCuts->SetDCAToVertex2D(kTRUE);
295 esdTrackCuts->SetPtRange(minPt,maxPt);
297 TString tag = "TPC-only tracking";
300 // TPC-only (loose cuts, no pt cut)
303 // beta cuts (still under investigation)
304 minNClustersTPC = 50;
305 maxChi2PerClusterTPC = 4.0;
306 maxDCAtoVertexXY = 1.e10; // cm
307 maxDCAtoVertexZ = 1.e10; // cm
311 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
312 esdTrackCuts->SetRequireTPCRefit(kFALSE);
313 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
314 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
315 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
316 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
317 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
318 esdTrackCuts->SetDCAToVertex2D(kTRUE);
319 esdTrackCuts->SetPtRange(minPt,maxPt);
321 TString tag = "TPC-only tracking";
324 // TPC-only (standard cuts, no pt cut)
327 // beta cuts (still under investigation)
328 minNClustersTPC = 96;
329 maxChi2PerClusterTPC = 3.5;
330 maxDCAtoVertexXY = 1.e10; // cm
331 maxDCAtoVertexZ = 1.e10; // cm
335 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
336 esdTrackCuts->SetRequireTPCRefit(kFALSE);
337 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
338 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
339 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
340 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
341 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
342 esdTrackCuts->SetDCAToVertex2D(kTRUE);
343 esdTrackCuts->SetPtRange(minPt,maxPt);
345 TString tag = "TPC-only tracking";
348 // TPC-only (tight cuts, no pt cuts)
351 // beta cuts (still under investigation)
352 minNClustersTPC = 120;
353 maxChi2PerClusterTPC = 3.5;
354 maxDCAtoVertexXY = 1.e10; // cm
355 maxDCAtoVertexZ = 1.e10; // cm
359 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
360 esdTrackCuts->SetRequireTPCRefit(kFALSE);
361 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
362 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
363 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
364 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
365 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
366 esdTrackCuts->SetDCAToVertex2D(kTRUE);
367 esdTrackCuts->SetPtRange(minPt,maxPt);
369 TString tag = "TPC-only tracking";
371 // TPC-only (loose cuts)
374 // beta cuts (still under investigation)
375 minNClustersTPC = 50;
376 maxChi2PerClusterTPC = 4.0;
377 //maxDCAtoVertexXY = 2.4; // cm
378 //maxDCAtoVertexZ = 3.2; // cm
379 maxDCAtoVertexXY = 1.6; // cm
380 maxDCAtoVertexZ = 2.1; // cm
384 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
385 esdTrackCuts->SetRequireTPCRefit(kFALSE);
386 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
387 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
388 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
389 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
390 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
391 esdTrackCuts->SetDCAToVertex2D(kTRUE);
392 esdTrackCuts->SetPtRange(minPt,maxPt);
394 TString tag = "TPC-only tracking";
397 // TPC-only (standard cuts)
400 // beta cuts (still under investigation)
401 minNClustersTPC = 96;
402 maxChi2PerClusterTPC = 3.5;
403 //maxDCAtoVertexXY = 2.4; // cm
404 //maxDCAtoVertexZ = 3.2; // cm
405 maxDCAtoVertexXY = 1.4; // cm
406 maxDCAtoVertexZ = 1.8; // cm
410 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
411 esdTrackCuts->SetRequireTPCRefit(kFALSE);
412 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
413 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
414 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
415 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
416 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
417 esdTrackCuts->SetDCAToVertex2D(kTRUE);
418 esdTrackCuts->SetPtRange(minPt,maxPt);
420 TString tag = "TPC-only tracking";
423 // TPC-only (tight cuts)
426 // beta cuts (still under investigation)
427 minNClustersTPC = 120;
428 maxChi2PerClusterTPC = 3.0;
429 //maxDCAtoVertexXY = 2.4; // cm
430 //maxDCAtoVertexZ = 3.2; // cm
431 maxDCAtoVertexXY = 1.4; // cm
432 maxDCAtoVertexZ = 1.8; // cm
436 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
437 esdTrackCuts->SetRequireTPCRefit(kFALSE);
438 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
439 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
440 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
441 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
442 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
443 esdTrackCuts->SetDCAToVertex2D(kTRUE);
444 esdTrackCuts->SetPtRange(minPt,maxPt);
446 TString tag = "TPC-only tracking";
449 // TPC-only (arb. cuts, kink cuts included)
452 // beta cuts (still under investigation)
453 minNClustersTPC = 50;
454 maxChi2PerClusterTPC = 1.e10;
455 maxDCAtoVertexXY = 3.0; // cm
456 maxDCAtoVertexZ = 3.0; // cm
460 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
461 esdTrackCuts->SetRequireTPCRefit(kFALSE);
462 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
463 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
464 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
465 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
466 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
467 esdTrackCuts->SetDCAToVertex2D(kTRUE);
468 esdTrackCuts->SetPtRange(minPt,maxPt);
470 TString tag = "TPC-only tracking";
473 // TPC-only (arb. cuts, kink cuts excluded)
476 // beta cuts (still under investigation)
477 minNClustersTPC = 50;
478 maxChi2PerClusterTPC = 1.e10;
479 maxDCAtoVertexXY = 3.0; // cm
480 maxDCAtoVertexZ = 3.0; // cm
484 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
485 esdTrackCuts->SetRequireTPCRefit(kFALSE);
486 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
487 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
488 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
489 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
490 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
491 esdTrackCuts->SetDCAToVertex2D(kTRUE);
492 esdTrackCuts->SetPtRange(minPt,maxPt);
494 TString tag = "TPC-only tracking";
497 // TPC-only (arb. cuts, kink cuts excluded, no chi2, no DCA)
500 // beta cuts (still under investigation)
501 minNClustersTPC = 50;
502 maxChi2PerClusterTPC = 1.e10;
503 maxDCAtoVertexXY = 1.e10; // cm
504 maxDCAtoVertexZ = 1.e10; // cm
508 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
509 esdTrackCuts->SetRequireTPCRefit(kFALSE);
510 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
511 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
512 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
513 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
514 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
515 esdTrackCuts->SetDCAToVertex2D(kTRUE);
516 esdTrackCuts->SetPtRange(minPt,maxPt);
518 TString tag = "TPC-only tracking";
524 // beta cuts (still under investigation)
525 minNClustersTPC = 70;
526 maxChi2PerClusterTPC = 4.0;
527 maxDCAtoVertexXY = 2.4; // cm
528 maxDCAtoVertexZ = 3.2; // cm
530 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
531 esdTrackCuts->SetRequireTPCRefit(kFALSE);
532 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
533 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
534 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
535 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
536 esdTrackCuts->SetRequireITSRefit(kFALSE);
537 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
538 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
539 esdTrackCuts->SetDCAToVertex2D(kTRUE);
540 //esdTrackCuts->SetPtRange(minPt,maxPt);
541 //esdTrackCuts->SetEtaRange(minEta,maxEta);
543 TString tag = "TPC-only tracking";
546 // TPC-only tight cuts
549 // beta cuts (still under investigation)
550 minNClustersTPC = 70;
551 maxChi2PerClusterTPC = 4.0;
552 maxDCAtoVertexXY = 0.3; // cm
553 maxDCAtoVertexZ = 0.3; // cm
555 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
556 esdTrackCuts->SetRequireTPCRefit(kFALSE);
557 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
558 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
559 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
560 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
561 esdTrackCuts->SetRequireITSRefit(kFALSE);
562 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
563 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
564 esdTrackCuts->SetDCAToVertex2D(kTRUE);
565 //esdTrackCuts->SetPtRange(minPt,maxPt);
566 //esdTrackCuts->SetEtaRange(minEta,maxEta);
568 TString tag = "TPC-only tracking";
572 // TPC (no pt cut, no eta cut)
575 // beta cuts (still under investigation)
576 minNClustersTPC = 70;
577 maxChi2PerClusterTPC = 4.0;
578 maxDCAtoVertexXY = 2.4; // cm
579 maxDCAtoVertexZ = 3.2; // cm
581 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
582 esdTrackCuts->SetRequireTPCRefit(kFALSE);
583 esdTrackCuts->SetRequireTPCStandAlone(kFALSE);
584 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
585 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
586 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
587 esdTrackCuts->SetRequireITSRefit(kFALSE);
588 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
589 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
590 esdTrackCuts->SetDCAToVertex2D(kTRUE);
592 TString tag = "TPC tracking";
595 // TPC-only (no pt cut, no eta cut) updated 2011
598 // beta cuts (still under investigation)
599 //minNClustersTPC = 50;
600 Float_t minNCrossedRowsTPC = 120;
601 Float_t minRatioCrossedRowsOverFindableClustersTPC = 0.8;
602 Float_t maxFractionSharedTPCCluster = 0.4;
603 maxChi2PerClusterTPC = 4.0;
604 maxDCAtoVertexXY = 2.4; // cm
605 maxDCAtoVertexZ = 3.2; // cm
609 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
610 esdTrackCuts->SetRequireTPCRefit(kFALSE);
611 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
612 //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
614 esdTrackCuts->SetMinNCrossedRowsTPC(minNCrossedRowsTPC);
615 esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(minRatioCrossedRowsOverFindableClustersTPC);
616 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
617 esdTrackCuts->SetMaxFractionSharedTPCClusters(maxFractionSharedTPCCluster);
619 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
620 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
621 esdTrackCuts->SetDCAToVertex2D(kTRUE);
622 esdTrackCuts->SetPtRange(minPt,maxPt);
624 TString tag = "TPC-only tracking (2011)";
629 // systematic errors DCA cut studies
634 // beta cuts (still under investigation)
635 minNClustersTPC = 50;
636 maxChi2PerClusterTPC = 4.0;
637 maxDCAtoVertexXY = 1.4; // cm
638 maxDCAtoVertexZ = 2.2; // cm
642 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
643 esdTrackCuts->SetRequireTPCRefit(kFALSE);
644 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
645 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
646 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
647 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
648 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
649 esdTrackCuts->SetDCAToVertex2D(kTRUE);
650 esdTrackCuts->SetPtRange(minPt,maxPt);
652 TString tag = "TPC-only tracking";
657 // beta cuts (still under investigation)
658 minNClustersTPC = 50;
659 maxChi2PerClusterTPC = 4.0;
660 maxDCAtoVertexXY = 1.6; // cm
661 maxDCAtoVertexZ = 2.4; // cm
665 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
666 esdTrackCuts->SetRequireTPCRefit(kFALSE);
667 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
668 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
669 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
670 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
671 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
672 esdTrackCuts->SetDCAToVertex2D(kTRUE);
673 esdTrackCuts->SetPtRange(minPt,maxPt);
675 TString tag = "TPC-only tracking";
679 // systematic errors cut studies
684 // beta cuts (still under investigation)
685 minNClustersTPC = 50;
686 maxChi2PerClusterTPC = 4.0;
687 maxDCAtoVertexXY = 1.8; // cm
688 maxDCAtoVertexZ = 2.6; // cm
692 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
693 esdTrackCuts->SetRequireTPCRefit(kFALSE);
694 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
695 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
696 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
697 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
698 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
699 esdTrackCuts->SetDCAToVertex2D(kTRUE);
700 esdTrackCuts->SetPtRange(minPt,maxPt);
702 TString tag = "TPC-only tracking";
707 // beta cuts (still under investigation)
708 minNClustersTPC = 50;
709 maxChi2PerClusterTPC = 4.0;
710 maxDCAtoVertexXY = 2.0; // cm
711 maxDCAtoVertexZ = 2.8; // cm
715 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
716 esdTrackCuts->SetRequireTPCRefit(kFALSE);
717 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
718 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
719 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
720 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
721 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
722 esdTrackCuts->SetDCAToVertex2D(kTRUE);
723 esdTrackCuts->SetPtRange(minPt,maxPt);
725 TString tag = "TPC-only tracking";
730 // beta cuts (still under investigation)
731 minNClustersTPC = 50;
732 maxChi2PerClusterTPC = 4.0;
733 maxDCAtoVertexXY = 2.2; // cm
734 maxDCAtoVertexZ = 3.0; // cm
738 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
739 esdTrackCuts->SetRequireTPCRefit(kFALSE);
740 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
741 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
742 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
743 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
744 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
745 esdTrackCuts->SetDCAToVertex2D(kTRUE);
746 esdTrackCuts->SetPtRange(minPt,maxPt);
748 TString tag = "TPC-only tracking";
753 // beta cuts (still under investigation)
754 minNClustersTPC = 50;
755 maxChi2PerClusterTPC = 4.0;
756 maxDCAtoVertexXY = 2.4; // cm
757 maxDCAtoVertexZ = 3.2; // cm
761 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
762 esdTrackCuts->SetRequireTPCRefit(kFALSE);
763 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
764 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
765 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
766 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
767 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
768 esdTrackCuts->SetDCAToVertex2D(kTRUE);
769 esdTrackCuts->SetPtRange(minPt,maxPt);
771 TString tag = "TPC-only tracking";
776 // beta cuts (still under investigation)
777 minNClustersTPC = 50;
778 maxChi2PerClusterTPC = 4.0;
779 maxDCAtoVertexXY = 2.6; // cm
780 maxDCAtoVertexZ = 3.4; // cm
784 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
785 esdTrackCuts->SetRequireTPCRefit(kFALSE);
786 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
787 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
788 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
789 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
790 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
791 esdTrackCuts->SetDCAToVertex2D(kTRUE);
792 esdTrackCuts->SetPtRange(minPt,maxPt);
794 TString tag = "TPC-only tracking";
800 // beta cuts (still under investigation)
801 minNClustersTPC = 50;
802 maxChi2PerClusterTPC = 4.0;
803 maxDCAtoVertexXY = 2.8; // cm
804 maxDCAtoVertexZ = 3.6; // cm
808 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
809 esdTrackCuts->SetRequireTPCRefit(kFALSE);
810 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
811 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
812 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
813 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
814 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
815 esdTrackCuts->SetDCAToVertex2D(kTRUE);
816 esdTrackCuts->SetPtRange(minPt,maxPt);
818 TString tag = "TPC-only tracking";
823 // beta cuts (still under investigation)
824 minNClustersTPC = 50;
825 maxChi2PerClusterTPC = 4.0;
826 maxDCAtoVertexXY = 3.0; // cm
827 maxDCAtoVertexZ = 3.8; // cm
831 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
832 esdTrackCuts->SetRequireTPCRefit(kFALSE);
833 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
834 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
835 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
836 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
837 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
838 esdTrackCuts->SetDCAToVertex2D(kTRUE);
839 esdTrackCuts->SetPtRange(minPt,maxPt);
841 TString tag = "TPC-only tracking";
846 // beta cuts (still under investigation)
847 minNClustersTPC = 50;
848 maxChi2PerClusterTPC = 4.0;
849 maxDCAtoVertexXY = 3.2; // cm
850 maxDCAtoVertexZ = 4.0; // cm
854 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
855 esdTrackCuts->SetRequireTPCRefit(kFALSE);
856 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
857 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
858 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
859 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
860 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
861 esdTrackCuts->SetDCAToVertex2D(kTRUE);
862 esdTrackCuts->SetPtRange(minPt,maxPt);
864 TString tag = "TPC-only tracking";
869 // beta cuts (still under investigation)
870 minNClustersTPC = 50;
871 maxChi2PerClusterTPC = 4.0;
872 maxDCAtoVertexXY = 3.4; // cm
873 maxDCAtoVertexZ = 4.2; // cm
877 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
878 esdTrackCuts->SetRequireTPCRefit(kFALSE);
879 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
880 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
881 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
882 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
883 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
884 esdTrackCuts->SetDCAToVertex2D(kTRUE);
885 esdTrackCuts->SetPtRange(minPt,maxPt);
887 TString tag = "TPC-only tracking";
891 // cut stability systematics
896 // beta cuts (still under investigation)
897 minNClustersTPC = 70;
898 maxChi2PerClusterTPC = 4.0;
899 maxDCAtoVertexXY = 2.4; // cm
900 maxDCAtoVertexZ = 3.2; // cm
904 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
905 esdTrackCuts->SetRequireTPCRefit(kFALSE);
906 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
907 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
908 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
909 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
910 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
911 esdTrackCuts->SetDCAToVertex2D(kTRUE);
912 esdTrackCuts->SetPtRange(minPt,maxPt);
914 TString tag = "TPC-only tracking";
919 // beta cuts (still under investigation)
920 minNClustersTPC = 90;
921 maxChi2PerClusterTPC = 4.0;
922 maxDCAtoVertexXY = 2.4; // cm
923 maxDCAtoVertexZ = 3.2; // cm
927 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
928 esdTrackCuts->SetRequireTPCRefit(kFALSE);
929 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
930 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
931 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
932 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
933 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
934 esdTrackCuts->SetDCAToVertex2D(kTRUE);
935 esdTrackCuts->SetPtRange(minPt,maxPt);
937 TString tag = "TPC-only tracking";
942 // beta cuts (still under investigation)
943 minNClustersTPC = 50;
944 maxChi2PerClusterTPC = 3.0;
945 maxDCAtoVertexXY = 2.4; // cm
946 maxDCAtoVertexZ = 3.2; // cm
950 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
951 esdTrackCuts->SetRequireTPCRefit(kFALSE);
952 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
953 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
954 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
955 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
956 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
957 esdTrackCuts->SetDCAToVertex2D(kTRUE);
958 esdTrackCuts->SetPtRange(minPt,maxPt);
960 TString tag = "TPC-only tracking";
965 // beta cuts (still under investigation)
966 minNClustersTPC = 50;
967 maxChi2PerClusterTPC = 5.0;
968 maxDCAtoVertexXY = 2.4; // cm
969 maxDCAtoVertexZ = 3.2; // cm
973 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
974 esdTrackCuts->SetRequireTPCRefit(kFALSE);
975 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
976 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
977 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
978 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
979 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
980 esdTrackCuts->SetDCAToVertex2D(kTRUE);
981 esdTrackCuts->SetPtRange(minPt,maxPt);
983 TString tag = "TPC-only tracking";
988 // beta cuts (still under investigation)
989 minNClustersTPC = 50;
990 maxChi2PerClusterTPC = 4.0;
991 maxDCAtoVertexXY = 1.4; // cm
992 maxDCAtoVertexZ = 2.2; // cm
996 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
997 esdTrackCuts->SetRequireTPCRefit(kFALSE);
998 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
999 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
1000 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
1001 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1002 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1003 esdTrackCuts->SetDCAToVertex2D(kTRUE);
1004 esdTrackCuts->SetPtRange(minPt,maxPt);
1006 TString tag = "TPC-only tracking";
1011 // beta cuts (still under investigation)
1012 minNClustersTPC = 50;
1013 maxChi2PerClusterTPC = 4.0;
1014 maxDCAtoVertexXY = 3.4; // cm
1015 maxDCAtoVertexZ = 4.2; // cm
1019 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1020 esdTrackCuts->SetRequireTPCRefit(kFALSE);
1021 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1022 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
1023 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
1024 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1025 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1026 esdTrackCuts->SetDCAToVertex2D(kTRUE);
1027 esdTrackCuts->SetPtRange(minPt,maxPt);
1029 TString tag = "TPC-only tracking";
1034 // beta cuts (still under investigation)
1035 minNClustersTPC = 50;
1036 maxChi2PerClusterTPC = 4.0;
1037 maxDCAtoVertexXY = 2.4; // cm
1038 maxDCAtoVertexZ = 3.2; // cm
1042 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1043 esdTrackCuts->SetRequireTPCRefit(kFALSE);
1044 //esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1045 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
1046 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
1047 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
1048 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1049 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1050 esdTrackCuts->SetDCAToVertex2D(kTRUE);
1051 esdTrackCuts->SetPtRange(minPt,maxPt);
1053 TString tag = "TPC-only tracking";
1058 // beta cuts (still under investigation)
1059 minNClustersTPC = 50;
1060 maxChi2PerClusterTPC = 4.0;
1061 //maxDCAtoVertexXY = 2.4; // cm
1062 //maxDCAtoVertexZ = 3.2; // cm
1066 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1067 esdTrackCuts->SetRequireTPCRefit(kFALSE);
1068 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1069 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
1070 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
1071 //esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1072 //esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1073 //esdTrackCuts->SetDCAToVertex2D(kTRUE);
1074 //esdTrackCuts->SetPtRange(minPt,maxPt);
1075 //esdTrackCuts->SetEtaRange(minEta,maxEta);
1077 TString tag = "TPC-only tracking";
1080 // TPC-only + pt cut + eta cut
1083 // beta cuts (still under investigation)
1084 //minNClustersTPC = 50;
1085 //maxChi2PerClusterTPC = 4.0;
1086 maxDCAtoVertexXY = 2.4; // cm
1087 maxDCAtoVertexZ = 3.2; // cm
1091 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1092 esdTrackCuts->SetRequireTPCRefit(kFALSE);
1093 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1094 //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
1095 //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
1096 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1097 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1098 esdTrackCuts->SetDCAToVertex2D(kTRUE);
1099 //esdTrackCuts->SetPtRange(minPt,maxPt);
1100 //esdTrackCuts->SetEtaRange(minEta,maxEta);
1102 TString tag = "TPC-only tracking";
1105 // TPC-tracks + SPD point + ITS refit
1109 Double_t maxchi2perTPCcl=4.;
1110 //Double_t maxEtaInAcc=0.8;
1111 Double_t maxdcaxyITSTPC=0.2;
1112 Double_t maxdcazITSTPC=1.e9;
1114 esdTrackCuts->SetMaxDCAToVertexXY(maxdcaxyITSTPC);
1115 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1116 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1117 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1118 esdTrackCuts->SetRequireITSRefit(kTRUE);
1119 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1120 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1121 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1122 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1123 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1124 //esdTrackCuts->SetEtaRange(-maxEtaInAcc,maxEtaInAcc);
1126 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster";
1129 // TPC-tracks + SPD point + ITS refit
1134 Double_t maxchi2perTPCcl=4.;
1135 Double_t maxdcaxyITSTPC=0.2;
1136 Double_t maxdcazITSTPC=1.e9;
1138 esdTrackCuts->SetMaxDCAToVertexXY(maxdcaxyITSTPC);
1139 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1140 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1141 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1142 esdTrackCuts->SetRequireITSRefit(kTRUE);
1143 esdTrackCuts->SetMinNClustersITS(minclsITS);
1144 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1145 //esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1146 esdTrackCuts->SetRequireTPCRefit(kTRUE);
1147 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1148 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1149 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1151 TString tag = "Global tracking: TPC refit + ITS refit + >3 ITS clusters + >=1 SPD cluster";
1155 // TPC-tracks + SPD point + ITS refit + DCAr(pt)
1159 Double_t maxchi2perTPCcl=4.;
1160 Double_t maxdcaxyITSTPC=1.e9;
1161 Double_t maxdcazITSTPC=1.e9;
1163 esdTrackCuts->SetMaxDCAToVertexXY(maxdcaxyITSTPC);
1164 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1165 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1166 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1167 esdTrackCuts->SetRequireITSRefit(kTRUE);
1168 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1169 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1170 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1171 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1172 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1174 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt)";
1178 // TPC-tracks + SPD point + ITS refit + DCAr(pt)
1182 Double_t maxchi2perTPCcl=4.;
1183 Double_t maxdcazITSTPC=1.e9;
1188 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1189 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1190 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1191 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1195 esdTrackCuts->SetRequireITSRefit(kTRUE);
1196 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1198 // primary selection
1200 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1201 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1202 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1204 // 7*(0.0050+0.0060/pt^0.9)
1205 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9");
1207 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt)";
1210 // TPC+ITS combine tracking + DCAr(pt) + DCAz(pt)
1214 Double_t maxchi2perTPCcl=4.;
1215 Double_t maxdcazITSTPC=1.e9;
1220 esdTrackCuts->SetRequireTPCRefit(kTRUE);
1221 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1222 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1223 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1227 esdTrackCuts->SetRequireITSRefit(kTRUE);
1228 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1230 // primary selection
1232 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1233 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1234 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1236 // DCArphi parametrization (LHC10c pass2)
1237 // 7*(0.0026+0.0050/pt^1.01)
1238 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01");
1240 // DCArphi parametrization (LHC10c pass2)
1241 // 7*(0.01+0.011/pt^0.72)
1242 esdTrackCuts->SetMaxDCAToVertexZPtDep("0.07+0.077/pt^0.72");
1244 TString tag = "TPC+ITS combine tracking + DCAr(pt) + DCAz(pt)";
1247 // TPC+ITS combine tracking + DCAr(pt) (2010)
1251 Double_t maxchi2perTPCcl=4.;
1252 Double_t maxdcazITSTPC=2.0;
1257 esdTrackCuts->SetRequireTPCRefit(kTRUE);
1258 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1259 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1260 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1264 esdTrackCuts->SetRequireITSRefit(kTRUE);
1265 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1267 // primary selection
1269 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1270 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1271 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1273 // DCArphi parametrization (LHC10c pass2)
1274 // 7*(0.0026+0.0050/pt^1.01)
1275 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01");
1277 TString tag = "TPC+ITS combine tracking + DCAr(pt) (2010)";
1280 // TPC+ITS combine tracking + DCAr(pt) (2011)
1283 //Int_t minclsTPC=70;
1284 Float_t minNCrossedRowsTPC = 120;
1285 Float_t minRatioCrossedRowsOverFindableClustersTPC = 0.8;
1286 Float_t maxFractionSharedTPCCluster = 0.4;
1287 Double_t maxchi2perTPCcl=4.;
1288 Double_t maxdcazITSTPC=2.0;
1293 esdTrackCuts->SetRequireTPCRefit(kTRUE);
1294 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1296 //esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1297 esdTrackCuts->SetMinNCrossedRowsTPC(minNCrossedRowsTPC);
1298 esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(minRatioCrossedRowsOverFindableClustersTPC);
1299 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1300 esdTrackCuts->SetMaxFractionSharedTPCClusters(maxFractionSharedTPCCluster);
1304 esdTrackCuts->SetRequireITSRefit(kTRUE);
1305 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1307 // primary selection
1309 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1310 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1311 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1313 // DCArphi parametrization (LHC10c pass2)
1314 // 7*(0.0026+0.0050/pt^1.01)
1315 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01");
1317 TString tag = "TPC+ITS combine tracking + DCAr(pt) (2011)";
1320 // TPC-tracks + SPD point + ITS refit + DCAr(pt) 4-sigma
1324 Double_t maxchi2perTPCcl=4.;
1325 Double_t maxdcazITSTPC=1.e9;
1330 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1331 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1332 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1333 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1337 esdTrackCuts->SetRequireITSRefit(kTRUE);
1338 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1340 // primary selection
1342 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1343 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1344 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1346 // 4*(0.0050+0.0060/pt^0.9)
1347 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.02+0.024/pt^0.9");
1349 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) 4-sigma";
1352 // TPC-tracks + SPD point + ITS refit + DCAr(pt) 10-sigma
1356 Double_t maxchi2perTPCcl=4.;
1357 Double_t maxdcazITSTPC=1.e9;
1362 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1363 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1364 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1365 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1369 esdTrackCuts->SetRequireITSRefit(kTRUE);
1370 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1372 // primary selection
1374 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1375 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1376 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1378 // 10*(0.0050+0.0060/pt^0.9)
1379 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.05+0.06/pt^0.9");
1381 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) 10 sigma";
1384 // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 60 TPCclust
1388 Double_t maxchi2perTPCcl=4.;
1389 Double_t maxdcazITSTPC=1.e9;
1394 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1395 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1396 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1397 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1401 esdTrackCuts->SetRequireITSRefit(kTRUE);
1402 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1404 // primary selection
1406 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1407 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1408 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1410 // 7*(0.0050+0.0060/pt^0.9)
1411 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9");
1413 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 60 TPCclust";
1416 // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 80 clusters
1420 Double_t maxchi2perTPCcl=4.;
1421 Double_t maxdcazITSTPC=1.e9;
1426 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1427 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1428 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1429 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1433 esdTrackCuts->SetRequireITSRefit(kTRUE);
1434 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1436 // primary selection
1438 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1439 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1440 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1442 // 7*(0.0050+0.0060/pt^0.9)
1443 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9");
1445 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 80 TPCclust";
1448 // TPC-tracks + SPD point + ITS refit + DCAr(pt) + TPCchi2=3.5
1452 Double_t maxchi2perTPCcl=3.5;
1453 Double_t maxdcazITSTPC=1.e9;
1458 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1459 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1460 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1461 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1465 esdTrackCuts->SetRequireITSRefit(kTRUE);
1466 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1468 // primary selection
1470 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1471 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1472 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1474 // 7*(0.0050+0.0060/pt^0.9)
1475 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9");
1477 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + TPCchi2 3.5";
1480 // TPC-tracks + SPD point + ITS refit + DCAr(pt) + TPCchi2=4.5
1484 Double_t maxchi2perTPCcl=4.5;
1485 Double_t maxdcazITSTPC=1.e9;
1490 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1491 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1492 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1493 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1497 esdTrackCuts->SetRequireITSRefit(kTRUE);
1498 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1500 // primary selection
1502 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1503 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1504 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1506 // 7*(0.0050+0.0060/pt^0.9)
1507 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9");
1509 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + TPCchi2 4.5";
1516 minNClustersTPC = 70;
1517 maxChi2PerClusterTPC = 4.0;
1518 maxDCAtoVertexXY = 1.e9; // cm
1519 maxDCAtoVertexZ = 1.e9; // cm
1521 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1522 esdTrackCuts->SetRequireTPCRefit(kFALSE);
1523 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1524 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1525 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
1526 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
1527 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1528 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1529 esdTrackCuts->SetDCAToVertex2D(kTRUE);
1531 TString tag = "TPC-tracks loose criteria";
1535 // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 50 TPCclust
1539 Double_t maxchi2perTPCcl=4.;
1540 Double_t maxdcazITSTPC=1.e9;
1545 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1546 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1547 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1548 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1552 esdTrackCuts->SetRequireITSRefit(kTRUE);
1553 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1555 // primary selection
1557 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1558 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1559 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1561 // 7*(0.0050+0.0060/pt^0.9)
1562 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9");
1564 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 60 TPCclust";
1567 // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 70 TPCclust + accept kink daughters
1571 Double_t maxchi2perTPCcl=4.;
1572 Double_t maxdcazITSTPC=1.e9;
1577 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1578 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
1579 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1580 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1584 esdTrackCuts->SetRequireITSRefit(kTRUE);
1585 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1587 // primary selection
1589 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1590 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1591 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1593 // 7*(0.0050+0.0060/pt^0.9)
1594 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9");
1596 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 60 TPCclust";
1599 // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 30 TPCclust + accept kink daughters
1603 Double_t maxchi2perTPCcl=4.;
1604 Double_t maxdcazITSTPC=1.e9;
1609 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1610 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
1611 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1612 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1616 esdTrackCuts->SetRequireITSRefit(kTRUE);
1617 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1619 // primary selection
1621 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1622 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1623 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1625 // 7*(0.0050+0.0060/pt^0.9)
1626 esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9");
1628 TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 60 TPCclust";
1631 // Adam Kisiel track selectiion
1635 Double_t maxchi2perTPCcl=4.;
1636 Double_t maxdcazITSTPC=0.25;
1637 Double_t maxdcaxyITSTPC=0.2;
1642 //esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1643 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1644 esdTrackCuts->SetMinNClustersTPC(minclsTPC);
1645 esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl);
1649 esdTrackCuts->SetRequireITSRefit(kTRUE);
1650 //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1652 // primary selection
1654 //esdTrackCuts->SetDCAToVertex2D(kFALSE);
1655 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1656 esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC);
1657 esdTrackCuts->SetMaxDCAToVertexXY(maxdcaxyITSTPC);
1659 // 7*(0.0050+0.0060/pt^0.9)
1660 //esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9");
1662 TString tag = "Adam Kisiel track selection";
1665 // TPC+ITS refit + SPD any
1672 esdTrackCuts->SetRequireTPCRefit(kTRUE);
1673 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1677 esdTrackCuts->SetRequireITSRefit(kTRUE);
1678 esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1681 TString tag = "TPC+ITS refit required - for cut studies";
1691 esdTrackCuts->SetRequireTPCRefit(kTRUE);
1692 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1696 //esdTrackCuts->SetRequireITSRefit(kTRUE);
1697 //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1700 TString tag = "TPC refit required - for cut studies";
1710 esdTrackCuts->SetRequireTPCRefit(kFALSE);
1711 esdTrackCuts->SetRequireITSRefit(kFALSE);
1712 esdTrackCuts->SetRequireTPCStandAlone(kTRUE);
1713 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1717 //esdTrackCuts->SetRequireITSRefit(kTRUE);
1718 //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1721 TString tag = "TPC stand alone - for cut studies";
1731 esdTrackCuts->SetRequireTPCRefit(kTRUE);
1732 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1736 esdTrackCuts->SetRequireITSRefit(kTRUE);
1737 //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1740 TString tag = "TPC+ITS refit and KinkRejection required - for cut studies";
1743 // TPC+ITS refit + TPC DCA rough cuts
1750 maxDCAtoVertexXY = 5.0; // cm
1751 maxDCAtoVertexZ = 5.0; // cm
1753 esdTrackCuts->SetRequireTPCRefit(kTRUE);
1754 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1755 esdTrackCuts->SetDCAToVertex2D(kTRUE);
1756 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1757 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1760 esdTrackCuts->SetRequireITSRefit(kTRUE);
1763 // Only TPC refit and KinksRemoval required
1769 esdTrackCuts->SetRequireTPCRefit(kTRUE);
1770 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1774 esdTrackCuts->SetRequireITSRefit(kFALSE);
1775 //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny);
1778 TString tag = "TPC refit + Kink rejection required - for cut studies";
1781 // TPC+ITS combine tracking + DCAr(pt) (2011)
1782 // used for step-by-step cut studies!
1783 if ((cutMode >= 2000) && (cutMode <= 2222))
1785 //Int_t minclsTPC=70;
1786 Float_t minNCrossedRowsTPC = 120;
1787 Float_t minRatioCrossedRowsOverFindableClustersTPC = 0.8;
1788 Float_t maxFractionSharedTPCCluster = 0.4;
1789 Double_t maxchi2perTPCcl=4.;
1790 Double_t maxdcazITSTPC=2.0;
1791 Double_t maxdaczTPC=3.0;
1792 Double_t maxdcaxyTPC=3.0;
1797 if (cutMode >= 2001) { esdTrackCuts->SetRequireTPCRefit(kTRUE); }
1800 if (cutMode >= 2002) { esdTrackCuts->SetMinNCrossedRowsTPC(minNCrossedRowsTPC); }
1801 if (cutMode >= 2003) { esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(minRatioCrossedRowsOverFindableClustersTPC); }
1802 if (cutMode >= 2004) { esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); }
1803 if (cutMode >= 2005) { esdTrackCuts->SetMaxFractionSharedTPCClusters(maxFractionSharedTPCCluster); }
1804 if (cutMode >= 2006) { esdTrackCuts->SetMaxDCAToVertexZ(maxdaczTPC); }
1805 if (cutMode >= 2007) { esdTrackCuts->SetMaxDCAToVertexXY(maxdcaxyTPC); }
1809 if (cutMode >= 2008) { esdTrackCuts->SetRequireITSRefit(kTRUE); }
1810 if (cutMode >= 2009) { esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); }
1811 if (cutMode >= 2010) { esdTrackCuts->SetMaxChi2PerClusterITS(36.); }
1813 // primary selection
1815 esdTrackCuts->SetDCAToVertex2D(kFALSE);
1816 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1817 if (cutMode >= 2011) { esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); }
1819 // DCArphi parametrization (LHC10c pass2)
1820 // 7*(0.0026+0.0050/pt^1.01)
1821 if (cutMode >= 2012) { esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01"); }
1822 if (cutMode >= 2013) { esdTrackCuts->SetAcceptKinkDaughters(kFALSE); }
1824 if (cutMode >= 2014) { esdTrackCuts->SetMaxChi2TPCConstrainedGlobal(36.); }
1826 TString tag = "for cut/efficiency studies)";
1834 // cuts for data without field
1838 tag += " without field";
1841 Printf("Created track cuts for: %s", tag.Data());
1843 return esdTrackCuts;