removed printout in Terminate
[u/mrichter/AliRoot.git] / PWG0 / dNdPt / macros / CreatedNdPtTrackCuts.C
1 /* $Id: CreateCuts.C,v 1.5 2008/01/11 08:28:52 jgrosseo Exp $ */
2
3 // this macro creates the track and event cuts used in this analysis
4
5 AliESDtrackCuts* CreatedNdPtTrackCuts(Int_t cutMode=1, Bool_t fieldOn = kTRUE, Bool_t hists = kTRUE)
6 {
7   AliESDtrackCuts* esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts");
8
9   if (hists)
10     esdTrackCuts->DefineHistograms(1);
11
12   Double_t cov1, cov2, cov3, cov4, cov5;
13   Double_t nSigma;
14   Double_t maxDCAtoVertex, maxDCAtoVertexXY, maxDCAtoVertexZ;
15   Double_t minNClustersTPC;
16   Double_t maxChi2PerClusterTPC;
17   Double_t minPt, maxPt;
18
19   // default cuts for ITS+TPC
20   if (cutMode == 0) 
21   {
22     cov1 = 2;
23     cov2 = 2;
24     cov3 = 0.5;
25     cov4 = 0.5;
26     cov5 = 2;
27     nSigma = 3;
28     minNClustersTPC = 50;
29     maxChi2PerClusterTPC = 3.5;
30
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);
38
39     TString tag("Global tracking");
40   }
41
42   // TPC-only cuts (vertex n sigma cut)
43   if (cutMode == 1) 
44   {
45     // beta cuts (still under investigation)
46     //cov1 = 4;
47     //cov2 = 4;
48     cov1 = 2;
49     cov2 = 2;
50     cov3 = 0.5;
51     cov4 = 0.5;
52     cov5 = 2;
53     nSigma = 4;
54     minNClustersTPC = 50;
55     maxChi2PerClusterTPC = 3.5;
56
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);
64
65     TString tag = "TPC-only tracking";
66   }
67
68   // TPC-only cuts (vertex maxDCAtoVertex cut)
69   if (cutMode == 2) 
70   {
71     // beta cuts (still under investigation)
72     maxDCAtoVertex = 3.0; // cm
73     minNClustersTPC = 50;
74     maxChi2PerClusterTPC = 3.5;
75
76     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
77     esdTrackCuts->SetMaxDCAToVertex(maxDCAtoVertex);    
78     esdTrackCuts->SetRequireTPCRefit(kFALSE);
79     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
80     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
81     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
82
83     TString tag = "TPC-only tracking";
84   }
85
86   // TPC-only no vertex cuts
87   if (cutMode == 3) 
88   {
89     // beta cuts (still under investigation)
90     minNClustersTPC = 50;
91     maxChi2PerClusterTPC = 3.5;
92
93     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
94     esdTrackCuts->SetRequireTPCRefit(kFALSE);
95     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
96     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
97     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
98
99     TString tag = "TPC-only tracking";
100   }
101
102   // TPC-only no cuts at all 
103   if (cutMode == 4) 
104   {
105
106     // beta cuts (still under investigation)
107     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
108     esdTrackCuts->SetRequireTPCRefit(kFALSE);
109     esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
110
111     TString tag = "TPC-only tracking";
112   }
113
114   // TPC-only no kink removal no chi2 
115   if (cutMode == 5) 
116   {
117     // beta cuts (still under investigation)
118     minNClustersTPC = 50;
119     //maxChi2PerClusterTPC = 3.5;
120
121     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
122     esdTrackCuts->SetRequireTPCRefit(kFALSE);
123     esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
124     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
125     //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
126
127     TString tag = "TPC-only tracking";
128   }
129
130   // TPC-only no kink removal 
131   if (cutMode == 6) 
132   {
133     // beta cuts (still under investigation)
134     minNClustersTPC = 50;
135     maxChi2PerClusterTPC = 3.5;
136
137     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
138     esdTrackCuts->SetRequireTPCRefit(kFALSE);
139     esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
140     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
141     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
142
143     TString tag = "TPC-only tracking";
144   }
145
146   // TPC-only no kink removal no minNClustersTPC 
147   if (cutMode == 7) 
148   {
149     // beta cuts (still under investigation)
150     //minNClustersTPC = 50;
151     maxChi2PerClusterTPC = 3.5;
152
153     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
154     esdTrackCuts->SetRequireTPCRefit(kFALSE);
155     esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
156     //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
157     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
158
159     TString tag = "TPC-only tracking";
160   }
161   // TPC-only no kink removal no minNClustersTPC 
162   if (cutMode == 8) 
163   {
164     // beta cuts (still under investigation)
165     //minNClustersTPC = 50;
166     maxChi2PerClusterTPC = 3.5;
167     maxDCAtoVertex = 3.0; // cm
168
169     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
170     esdTrackCuts->SetMaxDCAToVertex(maxDCAtoVertex);    
171     esdTrackCuts->SetRequireTPCRefit(kFALSE);
172     esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
173     //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
174     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
175
176     TString tag = "TPC-only tracking";
177   }
178
179   // TPC-only no kink removal no minNClustersTPC no maxChi2PerClusterTPC
180   if (cutMode == 9) 
181   {
182     // beta cuts (still under investigation)
183     //minNClustersTPC = 50;
184     //maxChi2PerClusterTPC = 3.5;
185     maxDCAtoVertex = 3.0; // cm
186
187     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
188     esdTrackCuts->SetMaxDCAToVertex(maxDCAtoVertex);    
189     esdTrackCuts->SetRequireTPCRefit(kFALSE);
190     esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
191     //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
192     //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
193
194     TString tag = "TPC-only tracking";
195   }
196
197   // TPC-only (loose cuts, absolute DCA cut) 
198   if (cutMode == 10) 
199   {
200     // beta cuts (still under investigation)
201     minNClustersTPC = 50;
202     maxChi2PerClusterTPC = 4.0;
203     maxDCAtoVertex = 2.8; // cm
204     minPt=0.15;
205     maxPt=1.e10;
206
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);
214
215     TString tag = "TPC-only tracking";
216   }
217
218
219   // TPC-only (loose cuts, no DCA cut) 
220   if (cutMode == 11) 
221   {
222     // beta cuts (still under investigation)
223     minNClustersTPC = 50;
224     maxChi2PerClusterTPC = 4.0;
225     maxDCAtoVertexXY = 1.e10; // cm
226     maxDCAtoVertexZ  = 1.e10; // cm
227     minPt=0.15;
228     maxPt=1.e10;
229
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);
239
240     TString tag = "TPC-only tracking";
241   }
242
243   // TPC-only (standard cuts, no DCA cut) 
244   if (cutMode == 12) 
245   {
246     // beta cuts (still under investigation)
247     minNClustersTPC = 96;
248     maxChi2PerClusterTPC = 3.5;
249     maxDCAtoVertexXY = 1.e10; // cm
250     maxDCAtoVertexZ  = 1.e10; // cm
251     minPt=0.2;
252     maxPt=1.e10;
253
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);
263
264     TString tag = "TPC-only tracking";
265   }
266
267   // TPC-only (tight cuts, no DCA cut) 
268   if (cutMode == 13) 
269   {
270     // beta cuts (still under investigation)
271     minNClustersTPC = 120;
272     maxChi2PerClusterTPC = 3.5;
273     maxDCAtoVertexXY = 1.e10; // cm
274     maxDCAtoVertexZ  = 1.e10; // cm
275     minPt=0.3;
276     maxPt=1.e10;
277
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);
287
288     TString tag = "TPC-only tracking";
289   }
290
291   // TPC-only (loose cuts, no pt cut) 
292   if (cutMode == 14) 
293   {
294     // beta cuts (still under investigation)
295     minNClustersTPC = 50;
296     maxChi2PerClusterTPC = 4.0;
297     maxDCAtoVertexXY = 1.e10; // cm
298     maxDCAtoVertexZ  = 1.e10; // cm
299     minPt=0.0;
300     maxPt=1.e10;
301
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);
311
312     TString tag = "TPC-only tracking";
313   }
314
315   // TPC-only (standard cuts, no pt cut) 
316   if (cutMode == 15) 
317   {
318     // beta cuts (still under investigation)
319     minNClustersTPC = 96;
320     maxChi2PerClusterTPC = 3.5;
321     maxDCAtoVertexXY = 1.e10; // cm
322     maxDCAtoVertexZ  = 1.e10; // cm
323     minPt=0.0;
324     maxPt=1.e10;
325
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);
335
336     TString tag = "TPC-only tracking";
337   }
338
339   // TPC-only (tight cuts, no pt cuts) 
340   if (cutMode == 16) 
341   {
342     // beta cuts (still under investigation)
343     minNClustersTPC = 120;
344     maxChi2PerClusterTPC = 3.5;
345     maxDCAtoVertexXY = 1.e10; // cm
346     maxDCAtoVertexZ  = 1.e10; // cm
347     minPt=0.0;
348     maxPt=1.e10;
349
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);
359
360     TString tag = "TPC-only tracking";
361   }
362   // TPC-only (loose cuts)
363   if (cutMode == 17) 
364   {
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
372     minPt=0.15;
373     maxPt=1.e10;
374
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);
384
385     TString tag = "TPC-only tracking";
386   }
387
388   // TPC-only (standard cuts)
389   if (cutMode == 18) 
390   {
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
398     minPt=0.2;
399     maxPt=1.e10;
400
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);
410
411     TString tag = "TPC-only tracking";
412   }
413
414  // TPC-only (tight cuts)
415   if (cutMode == 19) 
416   {
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
424     minPt=0.3;
425     maxPt=1.e10;
426
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);
436
437     TString tag = "TPC-only tracking";
438   }
439
440   // TPC-only (arb. cuts, kink cuts included)
441   if (cutMode == 20) 
442   {
443     // beta cuts (still under investigation)
444     minNClustersTPC = 50;
445     maxChi2PerClusterTPC = 1.e10;
446     maxDCAtoVertexXY = 3.0; // cm
447     maxDCAtoVertexZ  = 3.0; // cm
448     minPt=0.0;
449     maxPt=1.e10;
450
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);
460
461     TString tag = "TPC-only tracking";
462   }
463
464   // TPC-only (arb. cuts, kink cuts excluded)
465   if (cutMode == 21) 
466   {
467     // beta cuts (still under investigation)
468     minNClustersTPC = 50;
469     maxChi2PerClusterTPC = 1.e10;
470     maxDCAtoVertexXY = 3.0; // cm
471     maxDCAtoVertexZ  = 3.0; // cm
472     minPt=0.0;
473     maxPt=1.e10;
474
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);
484
485     TString tag = "TPC-only tracking";
486   }
487
488   // TPC-only (arb. cuts, kink cuts excluded, no chi2, no DCA)
489   if (cutMode == 22) 
490   {
491     // beta cuts (still under investigation)
492     minNClustersTPC = 50;
493     maxChi2PerClusterTPC = 1.e10;
494     maxDCAtoVertexXY = 1.e10; // cm
495     maxDCAtoVertexZ  = 1.e10; // cm
496     minPt=0.15;
497     maxPt=1.e10;
498
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);
508
509     TString tag = "TPC-only tracking";
510   }
511
512   // TPC-only + pt cut + eta cut 
513   if (cutMode == 23) 
514   {
515     // beta cuts (still under investigation)
516     minNClustersTPC = 50;
517     maxChi2PerClusterTPC = 4.0;
518     maxDCAtoVertexXY = 2.4; // cm
519     maxDCAtoVertexZ  = 3.2; // cm
520     //minPt=0.15;
521     //maxPt=1.e10;
522
523     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
524     esdTrackCuts->SetRequireTPCRefit(kFALSE);
525     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
526     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
527     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
528     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
529     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
530     esdTrackCuts->SetDCAToVertex2D(kTRUE);
531     //esdTrackCuts->SetPtRange(minPt,maxPt);
532     //esdTrackCuts->SetEtaRange(minEta,maxEta);
533
534     TString tag = "TPC-only tracking";
535   }
536
537   // TPC-only (no pt cut, no eta cut)
538   if (cutMode == 24) 
539   {
540     // beta cuts (still under investigation)
541     minNClustersTPC = 50;
542     maxChi2PerClusterTPC = 4.0;
543     maxDCAtoVertexXY = 2.4; // cm
544     maxDCAtoVertexZ  = 3.2; // cm
545     minPt=0.0;
546     maxPt=1.e10;
547
548     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
549     esdTrackCuts->SetRequireTPCRefit(kFALSE);
550     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
551     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
552     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
553     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
554     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
555     esdTrackCuts->SetDCAToVertex2D(kTRUE);
556     esdTrackCuts->SetPtRange(minPt,maxPt);
557
558     TString tag = "TPC-only tracking";
559   }
560
561   //
562   // systematic errors DCA cut studies
563   //
564   // TPC-only
565   if (cutMode == 25) 
566   {
567     // beta cuts (still under investigation)
568     minNClustersTPC = 50;
569     maxChi2PerClusterTPC = 4.0;
570     maxDCAtoVertexXY = 1.4; // cm
571     maxDCAtoVertexZ  = 2.2; // cm
572     minPt=0.0;
573     maxPt=1.e10;
574
575     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
576     esdTrackCuts->SetRequireTPCRefit(kFALSE);
577     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
578     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
579     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
580     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
581     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
582     esdTrackCuts->SetDCAToVertex2D(kTRUE);
583     esdTrackCuts->SetPtRange(minPt,maxPt);
584
585     TString tag = "TPC-only tracking";
586   }
587
588   if (cutMode == 26) 
589   {
590     // beta cuts (still under investigation)
591     minNClustersTPC = 50;
592     maxChi2PerClusterTPC = 4.0;
593     maxDCAtoVertexXY = 1.6; // cm
594     maxDCAtoVertexZ  = 2.4; // cm
595     minPt=0.0;
596     maxPt=1.e10;
597
598     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
599     esdTrackCuts->SetRequireTPCRefit(kFALSE);
600     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
601     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
602     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
603     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
604     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
605     esdTrackCuts->SetDCAToVertex2D(kTRUE);
606     esdTrackCuts->SetPtRange(minPt,maxPt);
607
608     TString tag = "TPC-only tracking";
609   }
610
611   //
612   // systematic errors cut studies
613   //
614   // TPC-only
615   if (cutMode == 27) 
616   {
617     // beta cuts (still under investigation)
618     minNClustersTPC = 50;
619     maxChi2PerClusterTPC = 4.0;
620     maxDCAtoVertexXY = 1.8; // cm
621     maxDCAtoVertexZ  = 2.6; // cm
622     minPt=0.0;
623     maxPt=1.e10;
624
625     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
626     esdTrackCuts->SetRequireTPCRefit(kFALSE);
627     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
628     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
629     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
630     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
631     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
632     esdTrackCuts->SetDCAToVertex2D(kTRUE);
633     esdTrackCuts->SetPtRange(minPt,maxPt);
634
635     TString tag = "TPC-only tracking";
636   }
637
638   if (cutMode == 28) 
639   {
640     // beta cuts (still under investigation)
641     minNClustersTPC = 50;
642     maxChi2PerClusterTPC = 4.0;
643     maxDCAtoVertexXY = 2.0; // cm
644     maxDCAtoVertexZ  = 2.8; // cm
645     minPt=0.0;
646     maxPt=1.e10;
647
648     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
649     esdTrackCuts->SetRequireTPCRefit(kFALSE);
650     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
651     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
652     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
653     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
654     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
655     esdTrackCuts->SetDCAToVertex2D(kTRUE);
656     esdTrackCuts->SetPtRange(minPt,maxPt);
657
658     TString tag = "TPC-only tracking";
659   }
660
661   if (cutMode == 29) 
662   {
663     // beta cuts (still under investigation)
664     minNClustersTPC = 50;
665     maxChi2PerClusterTPC = 4.0;
666     maxDCAtoVertexXY = 2.2; // cm
667     maxDCAtoVertexZ  = 3.0; // cm
668     minPt=0.0;
669     maxPt=1.e10;
670
671     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
672     esdTrackCuts->SetRequireTPCRefit(kFALSE);
673     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
674     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
675     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
676     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
677     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
678     esdTrackCuts->SetDCAToVertex2D(kTRUE);
679     esdTrackCuts->SetPtRange(minPt,maxPt);
680
681     TString tag = "TPC-only tracking";
682   }
683
684   if (cutMode == 30) 
685   {
686     // beta cuts (still under investigation)
687     minNClustersTPC = 50;
688     maxChi2PerClusterTPC = 4.0;
689     maxDCAtoVertexXY = 2.4; // cm
690     maxDCAtoVertexZ  = 3.2; // cm
691     minPt=0.0;
692     maxPt=1.e10;
693
694     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
695     esdTrackCuts->SetRequireTPCRefit(kFALSE);
696     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
697     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
698     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
699     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
700     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
701     esdTrackCuts->SetDCAToVertex2D(kTRUE);
702     esdTrackCuts->SetPtRange(minPt,maxPt);
703
704     TString tag = "TPC-only tracking";
705   }
706
707   if (cutMode == 31) 
708   {
709     // beta cuts (still under investigation)
710     minNClustersTPC = 50;
711     maxChi2PerClusterTPC = 4.0;
712     maxDCAtoVertexXY = 2.6; // cm
713     maxDCAtoVertexZ  = 3.4; // cm
714     minPt=0.0;
715     maxPt=1.e10;
716
717     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
718     esdTrackCuts->SetRequireTPCRefit(kFALSE);
719     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
720     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
721     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
722     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
723     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
724     esdTrackCuts->SetDCAToVertex2D(kTRUE);
725     esdTrackCuts->SetPtRange(minPt,maxPt);
726
727     TString tag = "TPC-only tracking";
728   }
729
730
731   if (cutMode == 32) 
732   {
733     // beta cuts (still under investigation)
734     minNClustersTPC = 50;
735     maxChi2PerClusterTPC = 4.0;
736     maxDCAtoVertexXY = 2.8; // cm
737     maxDCAtoVertexZ  = 3.6; // cm
738     minPt=0.0;
739     maxPt=1.e10;
740
741     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
742     esdTrackCuts->SetRequireTPCRefit(kFALSE);
743     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
744     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
745     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
746     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
747     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
748     esdTrackCuts->SetDCAToVertex2D(kTRUE);
749     esdTrackCuts->SetPtRange(minPt,maxPt);
750
751     TString tag = "TPC-only tracking";
752   }
753
754   if (cutMode == 33) 
755   {
756     // beta cuts (still under investigation)
757     minNClustersTPC = 50;
758     maxChi2PerClusterTPC = 4.0;
759     maxDCAtoVertexXY = 3.0; // cm
760     maxDCAtoVertexZ  = 3.8; // cm
761     minPt=0.0;
762     maxPt=1.e10;
763
764     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
765     esdTrackCuts->SetRequireTPCRefit(kFALSE);
766     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
767     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
768     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
769     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
770     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
771     esdTrackCuts->SetDCAToVertex2D(kTRUE);
772     esdTrackCuts->SetPtRange(minPt,maxPt);
773
774     TString tag = "TPC-only tracking";
775   }
776
777   if (cutMode == 34) 
778   {
779     // beta cuts (still under investigation)
780     minNClustersTPC = 50;
781     maxChi2PerClusterTPC = 4.0;
782     maxDCAtoVertexXY = 3.2; // cm
783     maxDCAtoVertexZ  = 4.0; // cm
784     minPt=0.0;
785     maxPt=1.e10;
786
787     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
788     esdTrackCuts->SetRequireTPCRefit(kFALSE);
789     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
790     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
791     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
792     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
793     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
794     esdTrackCuts->SetDCAToVertex2D(kTRUE);
795     esdTrackCuts->SetPtRange(minPt,maxPt);
796
797     TString tag = "TPC-only tracking";
798   }
799
800   if (cutMode == 35) 
801   {
802     // beta cuts (still under investigation)
803     minNClustersTPC = 50;
804     maxChi2PerClusterTPC = 4.0;
805     maxDCAtoVertexXY = 3.4; // cm
806     maxDCAtoVertexZ  = 4.2; // cm
807     minPt=0.0;
808     maxPt=1.e10;
809
810     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
811     esdTrackCuts->SetRequireTPCRefit(kFALSE);
812     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
813     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
814     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
815     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
816     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
817     esdTrackCuts->SetDCAToVertex2D(kTRUE);
818     esdTrackCuts->SetPtRange(minPt,maxPt);
819
820     TString tag = "TPC-only tracking";
821   }
822
823 //
824 // cut stability systematics
825 //
826
827   if (cutMode == 36) 
828   {
829     // beta cuts (still under investigation)
830     minNClustersTPC = 70;
831     maxChi2PerClusterTPC = 4.0;
832     maxDCAtoVertexXY = 2.4; // cm
833     maxDCAtoVertexZ  = 3.2; // cm
834     minPt=0.0;
835     maxPt=1.e10;
836
837     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
838     esdTrackCuts->SetRequireTPCRefit(kFALSE);
839     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
840     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
841     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
842     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
843     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
844     esdTrackCuts->SetDCAToVertex2D(kTRUE);
845     esdTrackCuts->SetPtRange(minPt,maxPt);
846
847     TString tag = "TPC-only tracking";
848   }
849
850  if (cutMode == 37) 
851   {
852     // beta cuts (still under investigation)
853     minNClustersTPC = 90;
854     maxChi2PerClusterTPC = 4.0;
855     maxDCAtoVertexXY = 2.4; // cm
856     maxDCAtoVertexZ  = 3.2; // cm
857     minPt=0.0;
858     maxPt=1.e10;
859
860     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
861     esdTrackCuts->SetRequireTPCRefit(kFALSE);
862     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
863     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
864     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
865     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
866     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
867     esdTrackCuts->SetDCAToVertex2D(kTRUE);
868     esdTrackCuts->SetPtRange(minPt,maxPt);
869
870     TString tag = "TPC-only tracking";
871   }
872
873   if (cutMode == 38) 
874   {
875     // beta cuts (still under investigation)
876     minNClustersTPC = 50;
877     maxChi2PerClusterTPC = 3.0;
878     maxDCAtoVertexXY = 2.4; // cm
879     maxDCAtoVertexZ  = 3.2; // cm
880     minPt=0.0;
881     maxPt=1.e10;
882
883     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
884     esdTrackCuts->SetRequireTPCRefit(kFALSE);
885     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
886     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
887     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
888     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
889     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
890     esdTrackCuts->SetDCAToVertex2D(kTRUE);
891     esdTrackCuts->SetPtRange(minPt,maxPt);
892
893     TString tag = "TPC-only tracking";
894   }
895
896   if (cutMode == 39) 
897   {
898     // beta cuts (still under investigation)
899     minNClustersTPC = 50;
900     maxChi2PerClusterTPC = 5.0;
901     maxDCAtoVertexXY = 2.4; // cm
902     maxDCAtoVertexZ  = 3.2; // cm
903     minPt=0.0;
904     maxPt=1.e10;
905
906     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
907     esdTrackCuts->SetRequireTPCRefit(kFALSE);
908     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
909     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
910     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
911     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
912     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
913     esdTrackCuts->SetDCAToVertex2D(kTRUE);
914     esdTrackCuts->SetPtRange(minPt,maxPt);
915
916     TString tag = "TPC-only tracking";
917   }
918
919   if (cutMode == 40) 
920   {
921     // beta cuts (still under investigation)
922     minNClustersTPC = 50;
923     maxChi2PerClusterTPC = 4.0;
924     maxDCAtoVertexXY = 1.4; // cm
925     maxDCAtoVertexZ  = 2.2; // cm
926     minPt=0.0;
927     maxPt=1.e10;
928
929     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
930     esdTrackCuts->SetRequireTPCRefit(kFALSE);
931     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
932     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
933     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
934     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
935     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
936     esdTrackCuts->SetDCAToVertex2D(kTRUE);
937     esdTrackCuts->SetPtRange(minPt,maxPt);
938
939     TString tag = "TPC-only tracking";
940   }
941
942   if (cutMode == 41) 
943   {
944     // beta cuts (still under investigation)
945     minNClustersTPC = 50;
946     maxChi2PerClusterTPC = 4.0;
947     maxDCAtoVertexXY = 3.4; // cm
948     maxDCAtoVertexZ  = 4.2; // cm
949     minPt=0.0;
950     maxPt=1.e10;
951
952     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
953     esdTrackCuts->SetRequireTPCRefit(kFALSE);
954     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
955     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
956     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
957     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
958     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
959     esdTrackCuts->SetDCAToVertex2D(kTRUE);
960     esdTrackCuts->SetPtRange(minPt,maxPt);
961
962     TString tag = "TPC-only tracking";
963   }
964
965   if (cutMode == 42) 
966   {
967     // beta cuts (still under investigation)
968     minNClustersTPC = 50;
969     maxChi2PerClusterTPC = 4.0;
970     maxDCAtoVertexXY = 2.4; // cm
971     maxDCAtoVertexZ  = 3.2; // cm
972     minPt=0.0;
973     maxPt=1.e10;
974
975     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
976     esdTrackCuts->SetRequireTPCRefit(kFALSE);
977     //esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
978     esdTrackCuts->SetAcceptKinkDaughters(kTRUE);
979     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
980     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
981     esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
982     esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
983     esdTrackCuts->SetDCAToVertex2D(kTRUE);
984     esdTrackCuts->SetPtRange(minPt,maxPt);
985
986     TString tag = "TPC-only tracking";
987   }
988   // test
989   if (cutMode == 43) 
990   {
991     // beta cuts (still under investigation)
992     minNClustersTPC = 50;
993     maxChi2PerClusterTPC = 4.0;
994     //maxDCAtoVertexXY = 2.4; // cm
995     //maxDCAtoVertexZ  = 3.2; // cm
996     //minPt=0.15;
997     //maxPt=1.e10;
998
999     esdTrackCuts->SetRequireSigmaToVertex(kFALSE);
1000     esdTrackCuts->SetRequireTPCRefit(kFALSE);
1001     esdTrackCuts->SetAcceptKinkDaughters(kFALSE);
1002     esdTrackCuts->SetMinNClustersTPC(minNClustersTPC);
1003     esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC);
1004     //esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY);
1005     //esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ);
1006     //esdTrackCuts->SetDCAToVertex2D(kTRUE);
1007     //esdTrackCuts->SetPtRange(minPt,maxPt);
1008     //esdTrackCuts->SetEtaRange(minEta,maxEta);
1009
1010     TString tag = "TPC-only tracking";
1011   }
1012
1013
1014   // cuts for data without field
1015   if (!fieldOn)
1016   {
1017     cov5 = 1e10;
1018     tag += " without field";
1019   }
1020
1021   Printf("Created track cuts for: %s", tag.Data());
1022
1023   return esdTrackCuts;
1024 }