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 AliESDtrackCuts* CreatedNdPtTrackCuts(Int_t cutMode=1, Bool_t fieldOn = kTRUE, Bool_t hists = kTRUE)
7 AliESDtrackCuts* esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts");
10 esdTrackCuts->DefineHistograms(1);
12 Double_t cov1, cov2, cov3, cov4, cov5;
14 Double_t maxDCAtoVertex, maxDCAtoVertexXY, maxDCAtoVertexZ;
15 Double_t minNClustersTPC;
16 Double_t maxChi2PerClusterTPC;
17 Double_t minPt, maxPt;
19 // default cuts for ITS+TPC
29 maxChi2PerClusterTPC = 3.5;
31 esdTrackCuts->SetMaxCovDiagonalElements(cov1, cov2, cov3, cov4, cov5);
32 esdTrackCuts->SetMinNsigmaToVertex(nSigma);
33 esdTrackCuts->SetRequireSigmaToVertex(kTRUE);
34 esdTrackCuts->SetRequireTPCRefit(kTRUE);
35 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
36 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
37 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
39 TString tag("Global tracking");
42 // TPC-only cuts (vertex n sigma cut)
45 // beta cuts (still under investigation)
55 maxChi2PerClusterTPC = 3.5;
57 esdTrackCuts->SetMaxCovDiagonalElements(cov1, cov2, cov3, cov4, cov5);
58 esdTrackCuts->SetMinNsigmaToVertex(nSigma);
59 esdTrackCuts->SetRequireSigmaToVertex(kTRUE);
60 esdTrackCuts->SetRequireTPCRefit(kFALSE);
61 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
62 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
63 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
65 TString tag = "TPC-only tracking";
68 // TPC-only cuts (vertex maxDCAtoVertex cut)
71 // beta cuts (still under investigation)
72 maxDCAtoVertex = 3.0; // cm
74 maxChi2PerClusterTPC = 3.5;
76 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
77 esdTrackCuts->SetMaxDCAToVertex(maxDCAtoVertex);
78 esdTrackCuts->SetRequireTPCRefit(kFALSE);
79 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
80 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
81 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
83 TString tag = "TPC-only tracking";
86 // TPC-only no vertex cuts
89 // beta cuts (still under investigation)
91 maxChi2PerClusterTPC = 3.5;
93 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
94 esdTrackCuts->SetRequireTPCRefit(kFALSE);
95 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
96 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
97 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
99 TString tag = "TPC-only tracking";
102 // TPC-only no cuts at all
106 // beta cuts (still under investigation)
107 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
108 esdTrackCuts->SetRequireTPCRefit(kFALSE);
109 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
111 TString tag = "TPC-only tracking";
114 // TPC-only no kink removal no chi2
117 // beta cuts (still under investigation)
118 minNClustersTPC = 50;
119 //maxChi2PerClusterTPC = 3.5;
121 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
122 esdTrackCuts->SetRequireTPCRefit(kFALSE);
123 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
124 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
125 //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
127 TString tag = "TPC-only tracking";
130 // TPC-only no kink removal
133 // beta cuts (still under investigation)
134 minNClustersTPC = 50;
135 maxChi2PerClusterTPC = 3.5;
137 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
138 esdTrackCuts->SetRequireTPCRefit(kFALSE);
139 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
140 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
141 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
143 TString tag = "TPC-only tracking";
146 // TPC-only no kink removal no minNClustersTPC
149 // beta cuts (still under investigation)
150 //minNClustersTPC = 50;
151 maxChi2PerClusterTPC = 3.5;
153 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
154 esdTrackCuts->SetRequireTPCRefit(kFALSE);
155 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
156 //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
157 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
159 TString tag = "TPC-only tracking";
161 // TPC-only no kink removal no minNClustersTPC
164 // beta cuts (still under investigation)
165 //minNClustersTPC = 50;
166 maxChi2PerClusterTPC = 3.5;
167 maxDCAtoVertex = 3.0; // cm
169 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
170 esdTrackCuts->SetMaxDCAToVertex(maxDCAtoVertex);
171 esdTrackCuts->SetRequireTPCRefit(kFALSE);
172 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
173 //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
174 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
176 TString tag = "TPC-only tracking";
179 // TPC-only no kink removal no minNClustersTPC no maxChi2PerClusterTPC
182 // beta cuts (still under investigation)
183 //minNClustersTPC = 50;
184 //maxChi2PerClusterTPC = 3.5;
185 maxDCAtoVertex = 3.0; // cm
187 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
188 esdTrackCuts->SetMaxDCAToVertex(maxDCAtoVertex);
189 esdTrackCuts->SetRequireTPCRefit(kFALSE);
190 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
191 //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
192 //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
194 TString tag = "TPC-only tracking";
197 // TPC-only (loose cuts, absolute DCA cut)
200 // beta cuts (still under investigation)
201 minNClustersTPC = 50;
202 maxChi2PerClusterTPC = 4.0;
203 maxDCAtoVertex = 2.8; // cm
207 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
208 esdTrackCuts->SetRequireTPCRefit(kFALSE);
209 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
210 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
211 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
212 esdTrackCuts->SetMaxDCAToVertex(maxDCAtoVertex);
213 esdTrackCuts->SetPtRange(minPt,maxPt);
215 TString tag = "TPC-only tracking";
219 // TPC-only (loose cuts, no DCA cut)
222 // beta cuts (still under investigation)
223 minNClustersTPC = 50;
224 maxChi2PerClusterTPC = 4.0;
225 maxDCAtoVertexXY = 1.e10; // cm
226 maxDCAtoVertexZ = 1.e10; // cm
230 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
231 esdTrackCuts->SetRequireTPCRefit(kFALSE);
232 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
233 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
234 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
235 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
236 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
237 esdTrackCuts->SetDCAToVertex2D(kTRUE);
238 esdTrackCuts->SetPtRange(minPt,maxPt);
240 TString tag = "TPC-only tracking";
243 // TPC-only (standard cuts, no DCA cut)
246 // beta cuts (still under investigation)
247 minNClustersTPC = 96;
248 maxChi2PerClusterTPC = 3.5;
249 maxDCAtoVertexXY = 1.e10; // cm
250 maxDCAtoVertexZ = 1.e10; // cm
254 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
255 esdTrackCuts->SetRequireTPCRefit(kFALSE);
256 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
257 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
258 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
259 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
260 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
261 esdTrackCuts->SetDCAToVertex2D(kTRUE);
262 esdTrackCuts->SetPtRange(minPt,maxPt);
264 TString tag = "TPC-only tracking";
267 // TPC-only (tight cuts, no DCA cut)
270 // beta cuts (still under investigation)
271 minNClustersTPC = 120;
272 maxChi2PerClusterTPC = 3.5;
273 maxDCAtoVertexXY = 1.e10; // cm
274 maxDCAtoVertexZ = 1.e10; // cm
278 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
279 esdTrackCuts->SetRequireTPCRefit(kFALSE);
280 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
281 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
282 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
283 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
284 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
285 esdTrackCuts->SetDCAToVertex2D(kTRUE);
286 esdTrackCuts->SetPtRange(minPt,maxPt);
288 TString tag = "TPC-only tracking";
291 // TPC-only (loose cuts, no pt cut)
294 // beta cuts (still under investigation)
295 minNClustersTPC = 50;
296 maxChi2PerClusterTPC = 4.0;
297 maxDCAtoVertexXY = 1.e10; // cm
298 maxDCAtoVertexZ = 1.e10; // cm
302 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
303 esdTrackCuts->SetRequireTPCRefit(kFALSE);
304 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
305 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
306 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
307 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
308 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
309 esdTrackCuts->SetDCAToVertex2D(kTRUE);
310 esdTrackCuts->SetPtRange(minPt,maxPt);
312 TString tag = "TPC-only tracking";
315 // TPC-only (standard cuts, no pt cut)
318 // beta cuts (still under investigation)
319 minNClustersTPC = 96;
320 maxChi2PerClusterTPC = 3.5;
321 maxDCAtoVertexXY = 1.e10; // cm
322 maxDCAtoVertexZ = 1.e10; // cm
326 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
327 esdTrackCuts->SetRequireTPCRefit(kFALSE);
328 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
329 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
330 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
331 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
332 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
333 esdTrackCuts->SetDCAToVertex2D(kTRUE);
334 esdTrackCuts->SetPtRange(minPt,maxPt);
336 TString tag = "TPC-only tracking";
339 // TPC-only (tight cuts, no pt cuts)
342 // beta cuts (still under investigation)
343 minNClustersTPC = 120;
344 maxChi2PerClusterTPC = 3.5;
345 maxDCAtoVertexXY = 1.e10; // cm
346 maxDCAtoVertexZ = 1.e10; // cm
350 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
351 esdTrackCuts->SetRequireTPCRefit(kFALSE);
352 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
353 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
354 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
355 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
356 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
357 esdTrackCuts->SetDCAToVertex2D(kTRUE);
358 esdTrackCuts->SetPtRange(minPt,maxPt);
360 TString tag = "TPC-only tracking";
362 // TPC-only (loose cuts)
365 // beta cuts (still under investigation)
366 minNClustersTPC = 50;
367 maxChi2PerClusterTPC = 4.0;
368 //maxDCAtoVertexXY = 2.4; // cm
369 //maxDCAtoVertexZ = 3.2; // cm
370 maxDCAtoVertexXY = 1.6; // cm
371 maxDCAtoVertexZ = 2.1; // cm
375 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
376 esdTrackCuts->SetRequireTPCRefit(kFALSE);
377 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
378 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
379 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
380 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
381 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
382 esdTrackCuts->SetDCAToVertex2D(kTRUE);
383 esdTrackCuts->SetPtRange(minPt,maxPt);
385 TString tag = "TPC-only tracking";
388 // TPC-only (standard cuts)
391 // beta cuts (still under investigation)
392 minNClustersTPC = 96;
393 maxChi2PerClusterTPC = 3.5;
394 //maxDCAtoVertexXY = 2.4; // cm
395 //maxDCAtoVertexZ = 3.2; // cm
396 maxDCAtoVertexXY = 1.4; // cm
397 maxDCAtoVertexZ = 1.8; // cm
401 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
402 esdTrackCuts->SetRequireTPCRefit(kFALSE);
403 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
404 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
405 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
406 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
407 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
408 esdTrackCuts->SetDCAToVertex2D(kTRUE);
409 esdTrackCuts->SetPtRange(minPt,maxPt);
411 TString tag = "TPC-only tracking";
414 // TPC-only (tight cuts)
417 // beta cuts (still under investigation)
418 minNClustersTPC = 120;
419 maxChi2PerClusterTPC = 3.0;
420 //maxDCAtoVertexXY = 2.4; // cm
421 //maxDCAtoVertexZ = 3.2; // cm
422 maxDCAtoVertexXY = 1.4; // cm
423 maxDCAtoVertexZ = 1.8; // cm
427 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
428 esdTrackCuts->SetRequireTPCRefit(kFALSE);
429 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
430 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
431 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
432 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
433 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
434 esdTrackCuts->SetDCAToVertex2D(kTRUE);
435 esdTrackCuts->SetPtRange(minPt,maxPt);
437 TString tag = "TPC-only tracking";
440 // TPC-only (arb. cuts, kink cuts included)
443 // beta cuts (still under investigation)
444 minNClustersTPC = 50;
445 maxChi2PerClusterTPC = 1.e10;
446 maxDCAtoVertexXY = 3.0; // cm
447 maxDCAtoVertexZ = 3.0; // cm
451 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
452 esdTrackCuts->SetRequireTPCRefit(kFALSE);
453 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
454 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
455 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
456 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
457 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
458 esdTrackCuts->SetDCAToVertex2D(kTRUE);
459 esdTrackCuts->SetPtRange(minPt,maxPt);
461 TString tag = "TPC-only tracking";
464 // TPC-only (arb. cuts, kink cuts excluded)
467 // beta cuts (still under investigation)
468 minNClustersTPC = 50;
469 maxChi2PerClusterTPC = 1.e10;
470 maxDCAtoVertexXY = 3.0; // cm
471 maxDCAtoVertexZ = 3.0; // cm
475 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
476 esdTrackCuts->SetRequireTPCRefit(kFALSE);
477 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
478 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
479 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
480 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
481 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
482 esdTrackCuts->SetDCAToVertex2D(kTRUE);
483 esdTrackCuts->SetPtRange(minPt,maxPt);
485 TString tag = "TPC-only tracking";
488 // TPC-only (arb. cuts, kink cuts excluded, no chi2, no DCA)
491 // beta cuts (still under investigation)
492 minNClustersTPC = 50;
493 maxChi2PerClusterTPC = 1.e10;
494 maxDCAtoVertexXY = 1.e10; // cm
495 maxDCAtoVertexZ = 1.e10; // cm
499 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
500 esdTrackCuts->SetRequireTPCRefit(kFALSE);
501 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
502 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
503 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
504 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
505 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
506 esdTrackCuts->SetDCAToVertex2D(kTRUE);
507 esdTrackCuts->SetPtRange(minPt,maxPt);
509 TString tag = "TPC-only tracking";
512 // TPC-only + pt cut + eta cut
515 // beta cuts (still under investigation)
516 //minNClustersTPC = 50;
517 minNClustersTPC = 70;
518 maxChi2PerClusterTPC = 4.0;
519 maxDCAtoVertexXY = 2.4; // cm
520 maxDCAtoVertexZ = 3.2; // cm
524 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
525 esdTrackCuts->SetRequireTPCRefit(kFALSE);
526 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
527 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
528 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
529 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
530 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
531 esdTrackCuts->SetDCAToVertex2D(kTRUE);
532 //esdTrackCuts->SetPtRange(minPt,maxPt);
533 //esdTrackCuts->SetEtaRange(minEta,maxEta);
535 TString tag = "TPC-only tracking";
538 // TPC-only (no pt cut, no eta cut)
541 // beta cuts (still under investigation)
542 minNClustersTPC = 50;
543 maxChi2PerClusterTPC = 4.0;
544 maxDCAtoVertexXY = 2.4; // cm
545 maxDCAtoVertexZ = 3.2; // cm
549 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
550 esdTrackCuts->SetRequireTPCRefit(kFALSE);
551 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
552 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
553 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
554 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
555 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
556 esdTrackCuts->SetDCAToVertex2D(kTRUE);
557 esdTrackCuts->SetPtRange(minPt,maxPt);
559 TString tag = "TPC-only tracking";
563 // systematic errors DCA cut studies
568 // beta cuts (still under investigation)
569 minNClustersTPC = 50;
570 maxChi2PerClusterTPC = 4.0;
571 maxDCAtoVertexXY = 1.4; // cm
572 maxDCAtoVertexZ = 2.2; // cm
576 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
577 esdTrackCuts->SetRequireTPCRefit(kFALSE);
578 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
579 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
580 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
581 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
582 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
583 esdTrackCuts->SetDCAToVertex2D(kTRUE);
584 esdTrackCuts->SetPtRange(minPt,maxPt);
586 TString tag = "TPC-only tracking";
591 // beta cuts (still under investigation)
592 minNClustersTPC = 50;
593 maxChi2PerClusterTPC = 4.0;
594 maxDCAtoVertexXY = 1.6; // cm
595 maxDCAtoVertexZ = 2.4; // cm
599 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
600 esdTrackCuts->SetRequireTPCRefit(kFALSE);
601 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
602 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
603 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
604 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
605 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
606 esdTrackCuts->SetDCAToVertex2D(kTRUE);
607 esdTrackCuts->SetPtRange(minPt,maxPt);
609 TString tag = "TPC-only tracking";
613 // systematic errors cut studies
618 // beta cuts (still under investigation)
619 minNClustersTPC = 50;
620 maxChi2PerClusterTPC = 4.0;
621 maxDCAtoVertexXY = 1.8; // cm
622 maxDCAtoVertexZ = 2.6; // cm
626 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
627 esdTrackCuts->SetRequireTPCRefit(kFALSE);
628 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
629 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
630 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
631 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
632 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
633 esdTrackCuts->SetDCAToVertex2D(kTRUE);
634 esdTrackCuts->SetPtRange(minPt,maxPt);
636 TString tag = "TPC-only tracking";
641 // beta cuts (still under investigation)
642 minNClustersTPC = 50;
643 maxChi2PerClusterTPC = 4.0;
644 maxDCAtoVertexXY = 2.0; // cm
645 maxDCAtoVertexZ = 2.8; // cm
649 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
650 esdTrackCuts->SetRequireTPCRefit(kFALSE);
651 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
652 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
653 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
654 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
655 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
656 esdTrackCuts->SetDCAToVertex2D(kTRUE);
657 esdTrackCuts->SetPtRange(minPt,maxPt);
659 TString tag = "TPC-only tracking";
664 // beta cuts (still under investigation)
665 minNClustersTPC = 50;
666 maxChi2PerClusterTPC = 4.0;
667 maxDCAtoVertexXY = 2.2; // cm
668 maxDCAtoVertexZ = 3.0; // cm
672 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
673 esdTrackCuts->SetRequireTPCRefit(kFALSE);
674 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
675 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
676 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
677 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
678 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
679 esdTrackCuts->SetDCAToVertex2D(kTRUE);
680 esdTrackCuts->SetPtRange(minPt,maxPt);
682 TString tag = "TPC-only tracking";
687 // beta cuts (still under investigation)
688 minNClustersTPC = 50;
689 maxChi2PerClusterTPC = 4.0;
690 maxDCAtoVertexXY = 2.4; // cm
691 maxDCAtoVertexZ = 3.2; // cm
695 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
696 esdTrackCuts->SetRequireTPCRefit(kFALSE);
697 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
698 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
699 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
700 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
701 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
702 esdTrackCuts->SetDCAToVertex2D(kTRUE);
703 esdTrackCuts->SetPtRange(minPt,maxPt);
705 TString tag = "TPC-only tracking";
710 // beta cuts (still under investigation)
711 minNClustersTPC = 50;
712 maxChi2PerClusterTPC = 4.0;
713 maxDCAtoVertexXY = 2.6; // cm
714 maxDCAtoVertexZ = 3.4; // cm
718 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
719 esdTrackCuts->SetRequireTPCRefit(kFALSE);
720 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
721 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
722 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
723 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
724 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
725 esdTrackCuts->SetDCAToVertex2D(kTRUE);
726 esdTrackCuts->SetPtRange(minPt,maxPt);
728 TString tag = "TPC-only tracking";
734 // beta cuts (still under investigation)
735 minNClustersTPC = 50;
736 maxChi2PerClusterTPC = 4.0;
737 maxDCAtoVertexXY = 2.8; // cm
738 maxDCAtoVertexZ = 3.6; // cm
742 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
743 esdTrackCuts->SetRequireTPCRefit(kFALSE);
744 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
745 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
746 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
747 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
748 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
749 esdTrackCuts->SetDCAToVertex2D(kTRUE);
750 esdTrackCuts->SetPtRange(minPt,maxPt);
752 TString tag = "TPC-only tracking";
757 // beta cuts (still under investigation)
758 minNClustersTPC = 50;
759 maxChi2PerClusterTPC = 4.0;
760 maxDCAtoVertexXY = 3.0; // cm
761 maxDCAtoVertexZ = 3.8; // cm
765 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
766 esdTrackCuts->SetRequireTPCRefit(kFALSE);
767 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
768 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
769 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
770 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
771 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
772 esdTrackCuts->SetDCAToVertex2D(kTRUE);
773 esdTrackCuts->SetPtRange(minPt,maxPt);
775 TString tag = "TPC-only tracking";
780 // beta cuts (still under investigation)
781 minNClustersTPC = 50;
782 maxChi2PerClusterTPC = 4.0;
783 maxDCAtoVertexXY = 3.2; // cm
784 maxDCAtoVertexZ = 4.0; // cm
788 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
789 esdTrackCuts->SetRequireTPCRefit(kFALSE);
790 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
791 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
792 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
793 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
794 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
795 esdTrackCuts->SetDCAToVertex2D(kTRUE);
796 esdTrackCuts->SetPtRange(minPt,maxPt);
798 TString tag = "TPC-only tracking";
803 // beta cuts (still under investigation)
804 minNClustersTPC = 50;
805 maxChi2PerClusterTPC = 4.0;
806 maxDCAtoVertexXY = 3.4; // cm
807 maxDCAtoVertexZ = 4.2; // cm
811 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
812 esdTrackCuts->SetRequireTPCRefit(kFALSE);
813 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
814 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
815 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
816 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
817 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
818 esdTrackCuts->SetDCAToVertex2D(kTRUE);
819 esdTrackCuts->SetPtRange(minPt,maxPt);
821 TString tag = "TPC-only tracking";
825 // cut stability systematics
830 // beta cuts (still under investigation)
831 minNClustersTPC = 70;
832 maxChi2PerClusterTPC = 4.0;
833 maxDCAtoVertexXY = 2.4; // cm
834 maxDCAtoVertexZ = 3.2; // cm
838 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
839 esdTrackCuts->SetRequireTPCRefit(kFALSE);
840 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
841 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
842 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
843 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
844 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
845 esdTrackCuts->SetDCAToVertex2D(kTRUE);
846 esdTrackCuts->SetPtRange(minPt,maxPt);
848 TString tag = "TPC-only tracking";
853 // beta cuts (still under investigation)
854 minNClustersTPC = 90;
855 maxChi2PerClusterTPC = 4.0;
856 maxDCAtoVertexXY = 2.4; // cm
857 maxDCAtoVertexZ = 3.2; // cm
861 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
862 esdTrackCuts->SetRequireTPCRefit(kFALSE);
863 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
864 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
865 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
866 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
867 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
868 esdTrackCuts->SetDCAToVertex2D(kTRUE);
869 esdTrackCuts->SetPtRange(minPt,maxPt);
871 TString tag = "TPC-only tracking";
876 // beta cuts (still under investigation)
877 minNClustersTPC = 50;
878 maxChi2PerClusterTPC = 3.0;
879 maxDCAtoVertexXY = 2.4; // cm
880 maxDCAtoVertexZ = 3.2; // cm
884 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
885 esdTrackCuts->SetRequireTPCRefit(kFALSE);
886 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
887 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
888 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
889 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
890 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
891 esdTrackCuts->SetDCAToVertex2D(kTRUE);
892 esdTrackCuts->SetPtRange(minPt,maxPt);
894 TString tag = "TPC-only tracking";
899 // beta cuts (still under investigation)
900 minNClustersTPC = 50;
901 maxChi2PerClusterTPC = 5.0;
902 maxDCAtoVertexXY = 2.4; // cm
903 maxDCAtoVertexZ = 3.2; // cm
907 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
908 esdTrackCuts->SetRequireTPCRefit(kFALSE);
909 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
910 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
911 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
912 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
913 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
914 esdTrackCuts->SetDCAToVertex2D(kTRUE);
915 esdTrackCuts->SetPtRange(minPt,maxPt);
917 TString tag = "TPC-only tracking";
922 // beta cuts (still under investigation)
923 minNClustersTPC = 50;
924 maxChi2PerClusterTPC = 4.0;
925 maxDCAtoVertexXY = 1.4; // cm
926 maxDCAtoVertexZ = 2.2; // cm
930 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
931 esdTrackCuts->SetRequireTPCRefit(kFALSE);
932 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
933 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
934 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
935 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
936 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
937 esdTrackCuts->SetDCAToVertex2D(kTRUE);
938 esdTrackCuts->SetPtRange(minPt,maxPt);
940 TString tag = "TPC-only tracking";
945 // beta cuts (still under investigation)
946 minNClustersTPC = 50;
947 maxChi2PerClusterTPC = 4.0;
948 maxDCAtoVertexXY = 3.4; // cm
949 maxDCAtoVertexZ = 4.2; // cm
953 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
954 esdTrackCuts->SetRequireTPCRefit(kFALSE);
955 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
956 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
957 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
958 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
959 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
960 esdTrackCuts->SetDCAToVertex2D(kTRUE);
961 esdTrackCuts->SetPtRange(minPt,maxPt);
963 TString tag = "TPC-only tracking";
968 // beta cuts (still under investigation)
969 minNClustersTPC = 50;
970 maxChi2PerClusterTPC = 4.0;
971 maxDCAtoVertexXY = 2.4; // cm
972 maxDCAtoVertexZ = 3.2; // cm
976 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
977 esdTrackCuts->SetRequireTPCRefit(kFALSE);
978 //esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
979 esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
980 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
981 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
982 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
983 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
984 esdTrackCuts->SetDCAToVertex2D(kTRUE);
985 esdTrackCuts->SetPtRange(minPt,maxPt);
987 TString tag = "TPC-only tracking";
992 // beta cuts (still under investigation)
993 minNClustersTPC = 50;
994 maxChi2PerClusterTPC = 4.0;
995 //maxDCAtoVertexXY = 2.4; // cm
996 //maxDCAtoVertexZ = 3.2; // cm
1000 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1001 esdTrackCuts->SetRequireTPCRefit(kFALSE);
1002 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1003 esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
1004 esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
1005 //esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1006 //esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1007 //esdTrackCuts->SetDCAToVertex2D(kTRUE);
1008 //esdTrackCuts->SetPtRange(minPt,maxPt);
1009 //esdTrackCuts->SetEtaRange(minEta,maxEta);
1011 TString tag = "TPC-only tracking";
1014 // TPC-only + pt cut + eta cut
1017 // beta cuts (still under investigation)
1018 //minNClustersTPC = 50;
1019 //maxChi2PerClusterTPC = 4.0;
1020 maxDCAtoVertexXY = 2.4; // cm
1021 maxDCAtoVertexZ = 3.2; // cm
1025 esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1026 esdTrackCuts->SetRequireTPCRefit(kFALSE);
1027 esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1028 //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
1029 //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
1030 esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1031 esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1032 esdTrackCuts->SetDCAToVertex2D(kTRUE);
1033 //esdTrackCuts->SetPtRange(minPt,maxPt);
1034 //esdTrackCuts->SetEtaRange(minEta,maxEta);
1036 TString tag = "TPC-only tracking";
1044 // cuts for data without field
1048 tag += " without field";
1051 Printf("Created track cuts for: %s", tag.Data());
1053 return esdTrackCuts;