]>
Commit | Line | Data |
---|---|---|
4fa29112 | 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 | ||
68f10917 | 5 | // last modified: 2011-03-28 |
6 | // m.l.knichel@gsi.de | |
7 | // added cut modes 200,201: replacing TPCNcluster cut | |
8 | ||
9 | ||
4fa29112 | 10 | AliESDtrackCuts* CreatedNdPtTrackCuts(Int_t cutMode=1, Bool_t fieldOn = kTRUE, Bool_t hists = kTRUE) |
11 | { | |
12 | AliESDtrackCuts* esdTrackCuts = new AliESDtrackCuts("AliESDtrackCuts"); | |
13 | ||
14 | if (hists) | |
15 | esdTrackCuts->DefineHistograms(1); | |
16 | ||
17 | Double_t cov1, cov2, cov3, cov4, cov5; | |
18 | Double_t nSigma; | |
19 | Double_t maxDCAtoVertex, maxDCAtoVertexXY, maxDCAtoVertexZ; | |
20 | Double_t minNClustersTPC; | |
21 | Double_t maxChi2PerClusterTPC; | |
22 | Double_t minPt, maxPt; | |
23 | ||
24 | // default cuts for ITS+TPC | |
25 | if (cutMode == 0) | |
26 | { | |
27 | cov1 = 2; | |
28 | cov2 = 2; | |
29 | cov3 = 0.5; | |
30 | cov4 = 0.5; | |
31 | cov5 = 2; | |
32 | nSigma = 3; | |
33 | minNClustersTPC = 50; | |
34 | maxChi2PerClusterTPC = 3.5; | |
35 | ||
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); | |
43 | ||
44 | TString tag("Global tracking"); | |
45 | } | |
46 | ||
47 | // TPC-only cuts (vertex n sigma cut) | |
48 | if (cutMode == 1) | |
49 | { | |
50 | // beta cuts (still under investigation) | |
51 | //cov1 = 4; | |
52 | //cov2 = 4; | |
53 | cov1 = 2; | |
54 | cov2 = 2; | |
55 | cov3 = 0.5; | |
56 | cov4 = 0.5; | |
57 | cov5 = 2; | |
58 | nSigma = 4; | |
59 | minNClustersTPC = 50; | |
60 | maxChi2PerClusterTPC = 3.5; | |
61 | ||
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); | |
69 | ||
70 | TString tag = "TPC-only tracking"; | |
71 | } | |
72 | ||
73 | // TPC-only cuts (vertex maxDCAtoVertex cut) | |
74 | if (cutMode == 2) | |
75 | { | |
76 | // beta cuts (still under investigation) | |
77 | maxDCAtoVertex = 3.0; // cm | |
78 | minNClustersTPC = 50; | |
79 | maxChi2PerClusterTPC = 3.5; | |
80 | ||
81 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
7f7dd416 | 82 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertex); |
83 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertex); | |
4fa29112 | 84 | esdTrackCuts->SetRequireTPCRefit(kFALSE); |
85 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
86 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
87 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
88 | ||
89 | TString tag = "TPC-only tracking"; | |
90 | } | |
91 | ||
92 | // TPC-only no vertex cuts | |
93 | if (cutMode == 3) | |
94 | { | |
95 | // beta cuts (still under investigation) | |
96 | minNClustersTPC = 50; | |
97 | maxChi2PerClusterTPC = 3.5; | |
98 | ||
99 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
100 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
101 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
102 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
103 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
104 | ||
105 | TString tag = "TPC-only tracking"; | |
106 | } | |
107 | ||
108 | // TPC-only no cuts at all | |
109 | if (cutMode == 4) | |
110 | { | |
111 | ||
112 | // beta cuts (still under investigation) | |
113 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
114 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
115 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
116 | ||
117 | TString tag = "TPC-only tracking"; | |
118 | } | |
119 | ||
120 | // TPC-only no kink removal no chi2 | |
121 | if (cutMode == 5) | |
122 | { | |
123 | // beta cuts (still under investigation) | |
124 | minNClustersTPC = 50; | |
125 | //maxChi2PerClusterTPC = 3.5; | |
126 | ||
127 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
128 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
129 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
130 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
131 | //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
132 | ||
133 | TString tag = "TPC-only tracking"; | |
134 | } | |
135 | ||
136 | // TPC-only no kink removal | |
137 | if (cutMode == 6) | |
138 | { | |
139 | // beta cuts (still under investigation) | |
140 | minNClustersTPC = 50; | |
141 | maxChi2PerClusterTPC = 3.5; | |
142 | ||
143 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
144 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
145 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
146 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
147 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
148 | ||
149 | TString tag = "TPC-only tracking"; | |
150 | } | |
151 | ||
152 | // TPC-only no kink removal no minNClustersTPC | |
153 | if (cutMode == 7) | |
154 | { | |
155 | // beta cuts (still under investigation) | |
156 | //minNClustersTPC = 50; | |
157 | maxChi2PerClusterTPC = 3.5; | |
158 | ||
159 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
160 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
161 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
162 | //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
163 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
164 | ||
165 | TString tag = "TPC-only tracking"; | |
166 | } | |
167 | // TPC-only no kink removal no minNClustersTPC | |
168 | if (cutMode == 8) | |
169 | { | |
170 | // beta cuts (still under investigation) | |
171 | //minNClustersTPC = 50; | |
172 | maxChi2PerClusterTPC = 3.5; | |
173 | maxDCAtoVertex = 3.0; // cm | |
174 | ||
175 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
7f7dd416 | 176 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertex); |
177 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertex); | |
4fa29112 | 178 | esdTrackCuts->SetRequireTPCRefit(kFALSE); |
179 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
180 | //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
181 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
182 | ||
183 | TString tag = "TPC-only tracking"; | |
184 | } | |
185 | ||
186 | // TPC-only no kink removal no minNClustersTPC no maxChi2PerClusterTPC | |
187 | if (cutMode == 9) | |
188 | { | |
189 | // beta cuts (still under investigation) | |
190 | //minNClustersTPC = 50; | |
191 | //maxChi2PerClusterTPC = 3.5; | |
192 | maxDCAtoVertex = 3.0; // cm | |
193 | ||
194 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
7f7dd416 | 195 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertex); |
196 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertex); | |
4fa29112 | 197 | esdTrackCuts->SetRequireTPCRefit(kFALSE); |
198 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
199 | //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
200 | //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
201 | ||
202 | TString tag = "TPC-only tracking"; | |
203 | } | |
204 | ||
205 | // TPC-only (loose cuts, absolute DCA cut) | |
206 | if (cutMode == 10) | |
207 | { | |
208 | // beta cuts (still under investigation) | |
209 | minNClustersTPC = 50; | |
210 | maxChi2PerClusterTPC = 4.0; | |
211 | maxDCAtoVertex = 2.8; // cm | |
212 | minPt=0.15; | |
213 | maxPt=1.e10; | |
214 | ||
215 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
216 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
217 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
218 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
219 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
7f7dd416 | 220 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertex); |
221 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertex); | |
4fa29112 | 222 | esdTrackCuts->SetPtRange(minPt,maxPt); |
223 | ||
224 | TString tag = "TPC-only tracking"; | |
225 | } | |
226 | ||
227 | ||
228 | // TPC-only (loose cuts, no DCA cut) | |
229 | if (cutMode == 11) | |
230 | { | |
231 | // beta cuts (still under investigation) | |
232 | minNClustersTPC = 50; | |
233 | maxChi2PerClusterTPC = 4.0; | |
234 | maxDCAtoVertexXY = 1.e10; // cm | |
235 | maxDCAtoVertexZ = 1.e10; // cm | |
236 | minPt=0.15; | |
237 | maxPt=1.e10; | |
238 | ||
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); | |
248 | ||
249 | TString tag = "TPC-only tracking"; | |
250 | } | |
251 | ||
252 | // TPC-only (standard cuts, no DCA cut) | |
253 | if (cutMode == 12) | |
254 | { | |
255 | // beta cuts (still under investigation) | |
256 | minNClustersTPC = 96; | |
257 | maxChi2PerClusterTPC = 3.5; | |
258 | maxDCAtoVertexXY = 1.e10; // cm | |
259 | maxDCAtoVertexZ = 1.e10; // cm | |
260 | minPt=0.2; | |
261 | maxPt=1.e10; | |
262 | ||
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); | |
272 | ||
273 | TString tag = "TPC-only tracking"; | |
274 | } | |
275 | ||
276 | // TPC-only (tight cuts, no DCA cut) | |
277 | if (cutMode == 13) | |
278 | { | |
279 | // beta cuts (still under investigation) | |
280 | minNClustersTPC = 120; | |
281 | maxChi2PerClusterTPC = 3.5; | |
282 | maxDCAtoVertexXY = 1.e10; // cm | |
283 | maxDCAtoVertexZ = 1.e10; // cm | |
284 | minPt=0.3; | |
285 | maxPt=1.e10; | |
286 | ||
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); | |
296 | ||
297 | TString tag = "TPC-only tracking"; | |
298 | } | |
299 | ||
300 | // TPC-only (loose cuts, no pt cut) | |
301 | if (cutMode == 14) | |
302 | { | |
303 | // beta cuts (still under investigation) | |
304 | minNClustersTPC = 50; | |
305 | maxChi2PerClusterTPC = 4.0; | |
306 | maxDCAtoVertexXY = 1.e10; // cm | |
307 | maxDCAtoVertexZ = 1.e10; // cm | |
308 | minPt=0.0; | |
309 | maxPt=1.e10; | |
310 | ||
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); | |
320 | ||
321 | TString tag = "TPC-only tracking"; | |
322 | } | |
323 | ||
324 | // TPC-only (standard cuts, no pt cut) | |
325 | if (cutMode == 15) | |
326 | { | |
327 | // beta cuts (still under investigation) | |
328 | minNClustersTPC = 96; | |
329 | maxChi2PerClusterTPC = 3.5; | |
330 | maxDCAtoVertexXY = 1.e10; // cm | |
331 | maxDCAtoVertexZ = 1.e10; // cm | |
332 | minPt=0.0; | |
333 | maxPt=1.e10; | |
334 | ||
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); | |
344 | ||
345 | TString tag = "TPC-only tracking"; | |
346 | } | |
347 | ||
348 | // TPC-only (tight cuts, no pt cuts) | |
349 | if (cutMode == 16) | |
350 | { | |
351 | // beta cuts (still under investigation) | |
352 | minNClustersTPC = 120; | |
353 | maxChi2PerClusterTPC = 3.5; | |
354 | maxDCAtoVertexXY = 1.e10; // cm | |
355 | maxDCAtoVertexZ = 1.e10; // cm | |
356 | minPt=0.0; | |
357 | maxPt=1.e10; | |
358 | ||
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); | |
368 | ||
369 | TString tag = "TPC-only tracking"; | |
370 | } | |
371 | // TPC-only (loose cuts) | |
372 | if (cutMode == 17) | |
373 | { | |
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 | |
381 | minPt=0.15; | |
382 | maxPt=1.e10; | |
383 | ||
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); | |
393 | ||
394 | TString tag = "TPC-only tracking"; | |
395 | } | |
396 | ||
397 | // TPC-only (standard cuts) | |
398 | if (cutMode == 18) | |
399 | { | |
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 | |
407 | minPt=0.2; | |
408 | maxPt=1.e10; | |
409 | ||
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); | |
419 | ||
420 | TString tag = "TPC-only tracking"; | |
421 | } | |
422 | ||
423 | // TPC-only (tight cuts) | |
424 | if (cutMode == 19) | |
425 | { | |
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 | |
433 | minPt=0.3; | |
434 | maxPt=1.e10; | |
435 | ||
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); | |
445 | ||
446 | TString tag = "TPC-only tracking"; | |
447 | } | |
448 | ||
449 | // TPC-only (arb. cuts, kink cuts included) | |
450 | if (cutMode == 20) | |
451 | { | |
452 | // beta cuts (still under investigation) | |
453 | minNClustersTPC = 50; | |
454 | maxChi2PerClusterTPC = 1.e10; | |
455 | maxDCAtoVertexXY = 3.0; // cm | |
456 | maxDCAtoVertexZ = 3.0; // cm | |
457 | minPt=0.0; | |
458 | maxPt=1.e10; | |
459 | ||
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); | |
469 | ||
470 | TString tag = "TPC-only tracking"; | |
471 | } | |
472 | ||
473 | // TPC-only (arb. cuts, kink cuts excluded) | |
474 | if (cutMode == 21) | |
475 | { | |
476 | // beta cuts (still under investigation) | |
477 | minNClustersTPC = 50; | |
478 | maxChi2PerClusterTPC = 1.e10; | |
479 | maxDCAtoVertexXY = 3.0; // cm | |
480 | maxDCAtoVertexZ = 3.0; // cm | |
481 | minPt=0.0; | |
482 | maxPt=1.e10; | |
483 | ||
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); | |
493 | ||
494 | TString tag = "TPC-only tracking"; | |
495 | } | |
496 | ||
497 | // TPC-only (arb. cuts, kink cuts excluded, no chi2, no DCA) | |
498 | if (cutMode == 22) | |
499 | { | |
500 | // beta cuts (still under investigation) | |
501 | minNClustersTPC = 50; | |
502 | maxChi2PerClusterTPC = 1.e10; | |
503 | maxDCAtoVertexXY = 1.e10; // cm | |
504 | maxDCAtoVertexZ = 1.e10; // cm | |
505 | minPt=0.15; | |
506 | maxPt=1.e10; | |
507 | ||
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); | |
517 | ||
518 | TString tag = "TPC-only tracking"; | |
519 | } | |
520 | ||
7f7dd416 | 521 | // TPC-only |
4fa29112 | 522 | if (cutMode == 23) |
523 | { | |
524 | // beta cuts (still under investigation) | |
a8ac5525 | 525 | minNClustersTPC = 70; |
4fa29112 | 526 | maxChi2PerClusterTPC = 4.0; |
527 | maxDCAtoVertexXY = 2.4; // cm | |
528 | maxDCAtoVertexZ = 3.2; // cm | |
4fa29112 | 529 | |
530 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
531 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
7f7dd416 | 532 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); |
4fa29112 | 533 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); |
534 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
535 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
0c6e7743 | 536 | esdTrackCuts->SetRequireITSRefit(kFALSE); |
4fa29112 | 537 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); |
538 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
539 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
540 | //esdTrackCuts->SetPtRange(minPt,maxPt); | |
541 | //esdTrackCuts->SetEtaRange(minEta,maxEta); | |
542 | ||
543 | TString tag = "TPC-only tracking"; | |
544 | } | |
68f10917 | 545 | |
546 | ||
4fa29112 | 547 | |
548 | // TPC-only (no pt cut, no eta cut) | |
549 | if (cutMode == 24) | |
550 | { | |
551 | // beta cuts (still under investigation) | |
552 | minNClustersTPC = 50; | |
553 | maxChi2PerClusterTPC = 4.0; | |
554 | maxDCAtoVertexXY = 2.4; // cm | |
555 | maxDCAtoVertexZ = 3.2; // cm | |
556 | minPt=0.0; | |
557 | maxPt=1.e10; | |
558 | ||
559 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
560 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
561 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
562 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
563 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
564 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); | |
565 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
566 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
567 | esdTrackCuts->SetPtRange(minPt,maxPt); | |
568 | ||
569 | TString tag = "TPC-only tracking"; | |
570 | } | |
68f10917 | 571 | |
572 | // TPC-only (no pt cut, no eta cut) updated 2011 | |
573 | if (cutMode == 201) | |
574 | { | |
575 | // beta cuts (still under investigation) | |
576 | //minNClustersTPC = 50; | |
577 | Float_t minNCrossedRowsTPC = 120; | |
578 | Float_t minRatioCrossedRowsOverFindableClustersTPC = 0.8; | |
579 | Float_t maxFractionSharedTPCCluster = 0.4; | |
580 | maxChi2PerClusterTPC = 4.0; | |
581 | maxDCAtoVertexXY = 2.4; // cm | |
582 | maxDCAtoVertexZ = 3.2; // cm | |
583 | minPt=0.0; | |
584 | maxPt=1.e10; | |
585 | ||
586 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
587 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
588 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
589 | //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
590 | ||
591 | esdTrackCuts->SetMinNCrossedRowsTPC(minNCrossedRowsTPC); | |
592 | esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(minRatioCrossedRowsOverFindableClustersTPC); | |
593 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
594 | esdTrackCuts->SetMaxFractionSharedTPCClusters(maxFractionSharedTPCCluster); | |
595 | ||
596 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); | |
597 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
598 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
599 | esdTrackCuts->SetPtRange(minPt,maxPt); | |
600 | ||
601 | TString tag = "TPC-only tracking (2011)"; | |
602 | } | |
603 | ||
4fa29112 | 604 | |
605 | // | |
606 | // systematic errors DCA cut studies | |
607 | // | |
608 | // TPC-only | |
609 | if (cutMode == 25) | |
610 | { | |
611 | // beta cuts (still under investigation) | |
612 | minNClustersTPC = 50; | |
613 | maxChi2PerClusterTPC = 4.0; | |
614 | maxDCAtoVertexXY = 1.4; // cm | |
615 | maxDCAtoVertexZ = 2.2; // cm | |
616 | minPt=0.0; | |
617 | maxPt=1.e10; | |
618 | ||
619 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
620 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
621 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
622 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
623 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
624 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); | |
625 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
626 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
627 | esdTrackCuts->SetPtRange(minPt,maxPt); | |
628 | ||
629 | TString tag = "TPC-only tracking"; | |
630 | } | |
631 | ||
632 | if (cutMode == 26) | |
633 | { | |
634 | // beta cuts (still under investigation) | |
635 | minNClustersTPC = 50; | |
636 | maxChi2PerClusterTPC = 4.0; | |
637 | maxDCAtoVertexXY = 1.6; // cm | |
638 | maxDCAtoVertexZ = 2.4; // cm | |
639 | minPt=0.0; | |
640 | maxPt=1.e10; | |
641 | ||
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); | |
651 | ||
652 | TString tag = "TPC-only tracking"; | |
653 | } | |
654 | ||
655 | // | |
656 | // systematic errors cut studies | |
657 | // | |
658 | // TPC-only | |
659 | if (cutMode == 27) | |
660 | { | |
661 | // beta cuts (still under investigation) | |
662 | minNClustersTPC = 50; | |
663 | maxChi2PerClusterTPC = 4.0; | |
664 | maxDCAtoVertexXY = 1.8; // cm | |
665 | maxDCAtoVertexZ = 2.6; // cm | |
666 | minPt=0.0; | |
667 | maxPt=1.e10; | |
668 | ||
669 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
670 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
671 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
672 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
673 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
674 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); | |
675 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
676 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
677 | esdTrackCuts->SetPtRange(minPt,maxPt); | |
678 | ||
679 | TString tag = "TPC-only tracking"; | |
680 | } | |
681 | ||
682 | if (cutMode == 28) | |
683 | { | |
684 | // beta cuts (still under investigation) | |
685 | minNClustersTPC = 50; | |
686 | maxChi2PerClusterTPC = 4.0; | |
687 | maxDCAtoVertexXY = 2.0; // cm | |
688 | maxDCAtoVertexZ = 2.8; // cm | |
689 | minPt=0.0; | |
690 | maxPt=1.e10; | |
691 | ||
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); | |
701 | ||
702 | TString tag = "TPC-only tracking"; | |
703 | } | |
704 | ||
705 | if (cutMode == 29) | |
706 | { | |
707 | // beta cuts (still under investigation) | |
708 | minNClustersTPC = 50; | |
709 | maxChi2PerClusterTPC = 4.0; | |
710 | maxDCAtoVertexXY = 2.2; // cm | |
711 | maxDCAtoVertexZ = 3.0; // cm | |
712 | minPt=0.0; | |
713 | maxPt=1.e10; | |
714 | ||
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); | |
724 | ||
725 | TString tag = "TPC-only tracking"; | |
726 | } | |
727 | ||
728 | if (cutMode == 30) | |
729 | { | |
730 | // beta cuts (still under investigation) | |
731 | minNClustersTPC = 50; | |
732 | maxChi2PerClusterTPC = 4.0; | |
733 | maxDCAtoVertexXY = 2.4; // cm | |
734 | maxDCAtoVertexZ = 3.2; // cm | |
735 | minPt=0.0; | |
736 | maxPt=1.e10; | |
737 | ||
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); | |
747 | ||
748 | TString tag = "TPC-only tracking"; | |
749 | } | |
750 | ||
751 | if (cutMode == 31) | |
752 | { | |
753 | // beta cuts (still under investigation) | |
754 | minNClustersTPC = 50; | |
755 | maxChi2PerClusterTPC = 4.0; | |
756 | maxDCAtoVertexXY = 2.6; // cm | |
757 | maxDCAtoVertexZ = 3.4; // cm | |
758 | minPt=0.0; | |
759 | maxPt=1.e10; | |
760 | ||
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); | |
770 | ||
771 | TString tag = "TPC-only tracking"; | |
772 | } | |
773 | ||
774 | ||
775 | if (cutMode == 32) | |
776 | { | |
777 | // beta cuts (still under investigation) | |
778 | minNClustersTPC = 50; | |
779 | maxChi2PerClusterTPC = 4.0; | |
780 | maxDCAtoVertexXY = 2.8; // cm | |
781 | maxDCAtoVertexZ = 3.6; // cm | |
782 | minPt=0.0; | |
783 | maxPt=1.e10; | |
784 | ||
785 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
786 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
787 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
788 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
789 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
790 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); | |
791 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
792 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
793 | esdTrackCuts->SetPtRange(minPt,maxPt); | |
794 | ||
795 | TString tag = "TPC-only tracking"; | |
796 | } | |
797 | ||
798 | if (cutMode == 33) | |
799 | { | |
800 | // beta cuts (still under investigation) | |
801 | minNClustersTPC = 50; | |
802 | maxChi2PerClusterTPC = 4.0; | |
803 | maxDCAtoVertexXY = 3.0; // cm | |
804 | maxDCAtoVertexZ = 3.8; // cm | |
805 | minPt=0.0; | |
806 | maxPt=1.e10; | |
807 | ||
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); | |
817 | ||
818 | TString tag = "TPC-only tracking"; | |
819 | } | |
820 | ||
821 | if (cutMode == 34) | |
822 | { | |
823 | // beta cuts (still under investigation) | |
824 | minNClustersTPC = 50; | |
825 | maxChi2PerClusterTPC = 4.0; | |
826 | maxDCAtoVertexXY = 3.2; // cm | |
827 | maxDCAtoVertexZ = 4.0; // cm | |
828 | minPt=0.0; | |
829 | maxPt=1.e10; | |
830 | ||
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); | |
840 | ||
841 | TString tag = "TPC-only tracking"; | |
842 | } | |
843 | ||
844 | if (cutMode == 35) | |
845 | { | |
846 | // beta cuts (still under investigation) | |
847 | minNClustersTPC = 50; | |
848 | maxChi2PerClusterTPC = 4.0; | |
849 | maxDCAtoVertexXY = 3.4; // cm | |
850 | maxDCAtoVertexZ = 4.2; // cm | |
851 | minPt=0.0; | |
852 | maxPt=1.e10; | |
853 | ||
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); | |
863 | ||
864 | TString tag = "TPC-only tracking"; | |
865 | } | |
866 | ||
867 | // | |
868 | // cut stability systematics | |
869 | // | |
870 | ||
871 | if (cutMode == 36) | |
872 | { | |
873 | // beta cuts (still under investigation) | |
874 | minNClustersTPC = 70; | |
875 | maxChi2PerClusterTPC = 4.0; | |
876 | maxDCAtoVertexXY = 2.4; // cm | |
877 | maxDCAtoVertexZ = 3.2; // cm | |
878 | minPt=0.0; | |
879 | maxPt=1.e10; | |
880 | ||
881 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
882 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
883 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
884 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
885 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
886 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); | |
887 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
888 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
889 | esdTrackCuts->SetPtRange(minPt,maxPt); | |
890 | ||
891 | TString tag = "TPC-only tracking"; | |
892 | } | |
893 | ||
894 | if (cutMode == 37) | |
895 | { | |
896 | // beta cuts (still under investigation) | |
897 | minNClustersTPC = 90; | |
898 | maxChi2PerClusterTPC = 4.0; | |
899 | maxDCAtoVertexXY = 2.4; // cm | |
900 | maxDCAtoVertexZ = 3.2; // cm | |
901 | minPt=0.0; | |
902 | maxPt=1.e10; | |
903 | ||
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); | |
913 | ||
914 | TString tag = "TPC-only tracking"; | |
915 | } | |
916 | ||
917 | if (cutMode == 38) | |
918 | { | |
919 | // beta cuts (still under investigation) | |
920 | minNClustersTPC = 50; | |
921 | maxChi2PerClusterTPC = 3.0; | |
922 | maxDCAtoVertexXY = 2.4; // cm | |
923 | maxDCAtoVertexZ = 3.2; // cm | |
924 | minPt=0.0; | |
925 | maxPt=1.e10; | |
926 | ||
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); | |
936 | ||
937 | TString tag = "TPC-only tracking"; | |
938 | } | |
939 | ||
940 | if (cutMode == 39) | |
941 | { | |
942 | // beta cuts (still under investigation) | |
943 | minNClustersTPC = 50; | |
944 | maxChi2PerClusterTPC = 5.0; | |
945 | maxDCAtoVertexXY = 2.4; // cm | |
946 | maxDCAtoVertexZ = 3.2; // cm | |
947 | minPt=0.0; | |
948 | maxPt=1.e10; | |
949 | ||
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); | |
959 | ||
960 | TString tag = "TPC-only tracking"; | |
961 | } | |
962 | ||
963 | if (cutMode == 40) | |
964 | { | |
965 | // beta cuts (still under investigation) | |
966 | minNClustersTPC = 50; | |
967 | maxChi2PerClusterTPC = 4.0; | |
968 | maxDCAtoVertexXY = 1.4; // cm | |
969 | maxDCAtoVertexZ = 2.2; // cm | |
970 | minPt=0.0; | |
971 | maxPt=1.e10; | |
972 | ||
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); | |
982 | ||
983 | TString tag = "TPC-only tracking"; | |
984 | } | |
985 | ||
986 | if (cutMode == 41) | |
987 | { | |
988 | // beta cuts (still under investigation) | |
989 | minNClustersTPC = 50; | |
990 | maxChi2PerClusterTPC = 4.0; | |
991 | maxDCAtoVertexXY = 3.4; // cm | |
992 | maxDCAtoVertexZ = 4.2; // cm | |
993 | minPt=0.0; | |
994 | maxPt=1.e10; | |
995 | ||
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); | |
1005 | ||
1006 | TString tag = "TPC-only tracking"; | |
1007 | } | |
1008 | ||
1009 | if (cutMode == 42) | |
1010 | { | |
1011 | // beta cuts (still under investigation) | |
1012 | minNClustersTPC = 50; | |
1013 | maxChi2PerClusterTPC = 4.0; | |
1014 | maxDCAtoVertexXY = 2.4; // cm | |
1015 | maxDCAtoVertexZ = 3.2; // cm | |
1016 | minPt=0.0; | |
1017 | maxPt=1.e10; | |
1018 | ||
1019 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1020 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
1021 | //esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1022 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
1023 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
1024 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
1025 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); | |
1026 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
1027 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
1028 | esdTrackCuts->SetPtRange(minPt,maxPt); | |
1029 | ||
1030 | TString tag = "TPC-only tracking"; | |
1031 | } | |
1032 | // test | |
1033 | if (cutMode == 43) | |
1034 | { | |
1035 | // beta cuts (still under investigation) | |
1036 | minNClustersTPC = 50; | |
1037 | maxChi2PerClusterTPC = 4.0; | |
1038 | //maxDCAtoVertexXY = 2.4; // cm | |
1039 | //maxDCAtoVertexZ = 3.2; // cm | |
1040 | //minPt=0.15; | |
1041 | //maxPt=1.e10; | |
1042 | ||
1043 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1044 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
1045 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
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); | |
1052 | //esdTrackCuts->SetEtaRange(minEta,maxEta); | |
1053 | ||
1054 | TString tag = "TPC-only tracking"; | |
1055 | } | |
1056 | ||
d269b0e6 | 1057 | // TPC-only + pt cut + eta cut |
1058 | if (cutMode == 45) | |
1059 | { | |
1060 | // beta cuts (still under investigation) | |
1061 | //minNClustersTPC = 50; | |
1062 | //maxChi2PerClusterTPC = 4.0; | |
1063 | maxDCAtoVertexXY = 2.4; // cm | |
1064 | maxDCAtoVertexZ = 3.2; // cm | |
1065 | //minPt=0.15; | |
1066 | //maxPt=1.e10; | |
1067 | ||
1068 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1069 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
1070 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1071 | //esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
1072 | //esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
1073 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); | |
1074 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
1075 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
1076 | //esdTrackCuts->SetPtRange(minPt,maxPt); | |
1077 | //esdTrackCuts->SetEtaRange(minEta,maxEta); | |
1078 | ||
1079 | TString tag = "TPC-only tracking"; | |
1080 | } | |
1081 | ||
7f7dd416 | 1082 | // TPC-tracks + SPD point + ITS refit |
1083 | if (cutMode == 50) | |
1084 | { | |
1085 | Int_t minclsTPC=70; | |
1086 | Double_t maxchi2perTPCcl=4.; | |
1087 | //Double_t maxEtaInAcc=0.8; | |
1088 | Double_t maxdcaxyITSTPC=0.2; | |
1089 | Double_t maxdcazITSTPC=1.e9; | |
1090 | ||
1091 | esdTrackCuts->SetMaxDCAToVertexXY(maxdcaxyITSTPC); | |
1092 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1093 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1094 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1095 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1096 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1097 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1098 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1099 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1100 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1101 | //esdTrackCuts->SetEtaRange(-maxEtaInAcc,maxEtaInAcc); | |
d269b0e6 | 1102 | |
7f7dd416 | 1103 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster"; |
1104 | } | |
d269b0e6 | 1105 | |
7f7dd416 | 1106 | // TPC-tracks + SPD point + ITS refit |
1107 | if (cutMode == 60) | |
1108 | { | |
1109 | Int_t minclsITS=4; | |
1110 | Int_t minclsTPC=70; | |
1111 | Double_t maxchi2perTPCcl=4.; | |
7f7dd416 | 1112 | Double_t maxdcaxyITSTPC=0.2; |
1113 | Double_t maxdcazITSTPC=1.e9; | |
1114 | ||
1115 | esdTrackCuts->SetMaxDCAToVertexXY(maxdcaxyITSTPC); | |
1116 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1117 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1118 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1119 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1120 | esdTrackCuts->SetMinNClustersITS(minclsITS); | |
1121 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1122 | //esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1123 | esdTrackCuts->SetRequireTPCRefit(kTRUE); | |
1124 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1125 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1126 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
d269b0e6 | 1127 | |
7f7dd416 | 1128 | TString tag = "Global tracking: TPC refit + ITS refit + >3 ITS clusters + >=1 SPD cluster"; |
1129 | } | |
4fa29112 | 1130 | |
695facdf | 1131 | /* |
791aaf54 | 1132 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) |
1133 | if (cutMode == 70) | |
1134 | { | |
1135 | Int_t minclsTPC=70; | |
1136 | Double_t maxchi2perTPCcl=4.; | |
791aaf54 | 1137 | Double_t maxdcaxyITSTPC=1.e9; |
1138 | Double_t maxdcazITSTPC=1.e9; | |
1139 | ||
1140 | esdTrackCuts->SetMaxDCAToVertexXY(maxdcaxyITSTPC); | |
1141 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1142 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1143 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1144 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1145 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1146 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1147 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1148 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1149 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
791aaf54 | 1150 | |
1151 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt)"; | |
1152 | } | |
695facdf | 1153 | */ |
1154 | ||
1155 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) | |
1156 | if (cutMode == 70) | |
1157 | { | |
1158 | Int_t minclsTPC=70; | |
1159 | Double_t maxchi2perTPCcl=4.; | |
1160 | Double_t maxdcazITSTPC=1.e9; | |
1161 | ||
1162 | // | |
1163 | // TPC | |
1164 | // | |
1165 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1166 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1167 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1168 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1169 | // | |
1170 | // ITS | |
1171 | // | |
1172 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1173 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1174 | // | |
1175 | // primary selection | |
1176 | // | |
1177 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1178 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1179 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1180 | ||
1181 | // 7*(0.0050+0.0060/pt^0.9) | |
1182 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); | |
1183 | ||
1184 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt)"; | |
1185 | } | |
1186 | ||
ecab7b3d | 1187 | // TPC+ITS combine tracking + DCAr(pt) + DCAz(pt) |
1188 | if (cutMode == 71) | |
1189 | { | |
1190 | Int_t minclsTPC=70; | |
1191 | Double_t maxchi2perTPCcl=4.; | |
68f10917 | 1192 | Double_t maxdcazITSTPC=1.e9; |
ecab7b3d | 1193 | |
1194 | // | |
1195 | // TPC | |
1196 | // | |
1197 | esdTrackCuts->SetRequireTPCRefit(kTRUE); | |
1198 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1199 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1200 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1201 | // | |
1202 | // ITS | |
1203 | // | |
1204 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1205 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1206 | // | |
1207 | // primary selection | |
1208 | // | |
1209 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1210 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
68f10917 | 1211 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); |
ecab7b3d | 1212 | |
1213 | // DCArphi parametrization (LHC10c pass2) | |
1214 | // 7*(0.0026+0.0050/pt^1.01) | |
1215 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01"); | |
1216 | ||
1217 | // DCArphi parametrization (LHC10c pass2) | |
1218 | // 7*(0.01+0.011/pt^0.72) | |
1219 | esdTrackCuts->SetMaxDCAToVertexZPtDep("0.07+0.077/pt^0.72"); | |
1220 | ||
1221 | TString tag = "TPC+ITS combine tracking + DCAr(pt) + DCAz(pt)"; | |
1222 | } | |
1223 | ||
82292e7f | 1224 | // TPC+ITS combine tracking + DCAr(pt) (2010) |
1225 | if (cutMode == 72) | |
1226 | { | |
1227 | Int_t minclsTPC=70; | |
1228 | Double_t maxchi2perTPCcl=4.; | |
1229 | Double_t maxdcazITSTPC=2.0; | |
1230 | ||
1231 | // | |
1232 | // TPC | |
1233 | // | |
1234 | esdTrackCuts->SetRequireTPCRefit(kTRUE); | |
1235 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1236 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1237 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1238 | // | |
1239 | // ITS | |
1240 | // | |
1241 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1242 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1243 | // | |
1244 | // primary selection | |
1245 | // | |
1246 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1247 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1248 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1249 | ||
1250 | // DCArphi parametrization (LHC10c pass2) | |
1251 | // 7*(0.0026+0.0050/pt^1.01) | |
1252 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01"); | |
1253 | ||
1254 | TString tag = "TPC+ITS combine tracking + DCAr(pt) (2010)"; | |
1255 | } | |
1256 | ||
68f10917 | 1257 | // TPC+ITS combine tracking + DCAr(pt) (2011) |
1258 | if (cutMode == 200) | |
1259 | { | |
1260 | //Int_t minclsTPC=70; | |
1261 | Float_t minNCrossedRowsTPC = 120; | |
1262 | Float_t minRatioCrossedRowsOverFindableClustersTPC = 0.8; | |
1263 | Float_t maxFractionSharedTPCCluster = 0.4; | |
1264 | Double_t maxchi2perTPCcl=4.; | |
1265 | Double_t maxdcazITSTPC=2.0; | |
1266 | ||
1267 | // | |
1268 | // TPC | |
1269 | // | |
1270 | esdTrackCuts->SetRequireTPCRefit(kTRUE); | |
1271 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1272 | ||
1273 | //esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1274 | esdTrackCuts->SetMinNCrossedRowsTPC(minNCrossedRowsTPC); | |
1275 | esdTrackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(minRatioCrossedRowsOverFindableClustersTPC); | |
1276 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1277 | esdTrackCuts->SetMaxFractionSharedTPCClusters(maxFractionSharedTPCCluster); | |
1278 | // | |
1279 | // ITS | |
1280 | // | |
1281 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1282 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1283 | // | |
1284 | // primary selection | |
1285 | // | |
1286 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1287 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1288 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1289 | ||
1290 | // DCArphi parametrization (LHC10c pass2) | |
1291 | // 7*(0.0026+0.0050/pt^1.01) | |
1292 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0182+0.0350/pt^1.01"); | |
1293 | ||
1294 | TString tag = "TPC+ITS combine tracking + DCAr(pt) (2011)"; | |
1295 | } | |
1296 | ||
695facdf | 1297 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) 4-sigma |
1298 | if (cutMode == 75) | |
1299 | { | |
1300 | Int_t minclsTPC=70; | |
1301 | Double_t maxchi2perTPCcl=4.; | |
1302 | Double_t maxdcazITSTPC=1.e9; | |
1303 | ||
1304 | // | |
1305 | // TPC | |
1306 | // | |
1307 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1308 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1309 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1310 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1311 | // | |
1312 | // ITS | |
1313 | // | |
1314 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1315 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1316 | // | |
1317 | // primary selection | |
1318 | // | |
1319 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1320 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1321 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1322 | ||
1323 | // 4*(0.0050+0.0060/pt^0.9) | |
1324 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.02+0.024/pt^0.9"); | |
1325 | ||
1326 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) 4-sigma"; | |
1327 | } | |
1328 | ||
1329 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) 10-sigma | |
1330 | if (cutMode == 80) | |
1331 | { | |
1332 | Int_t minclsTPC=70; | |
1333 | Double_t maxchi2perTPCcl=4.; | |
1334 | Double_t maxdcazITSTPC=1.e9; | |
1335 | ||
1336 | // | |
1337 | // TPC | |
1338 | // | |
1339 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1340 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1341 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1342 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1343 | // | |
1344 | // ITS | |
1345 | // | |
1346 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1347 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1348 | // | |
1349 | // primary selection | |
1350 | // | |
1351 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1352 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1353 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1354 | ||
1355 | // 10*(0.0050+0.0060/pt^0.9) | |
1356 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.05+0.06/pt^0.9"); | |
1357 | ||
1358 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) 10 sigma"; | |
1359 | } | |
1360 | ||
1361 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 60 TPCclust | |
1362 | if (cutMode == 85) | |
1363 | { | |
1364 | Int_t minclsTPC=60; | |
1365 | Double_t maxchi2perTPCcl=4.; | |
1366 | Double_t maxdcazITSTPC=1.e9; | |
1367 | ||
1368 | // | |
1369 | // TPC | |
1370 | // | |
1371 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1372 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1373 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1374 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1375 | // | |
1376 | // ITS | |
1377 | // | |
1378 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1379 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1380 | // | |
1381 | // primary selection | |
1382 | // | |
1383 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1384 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1385 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1386 | ||
1387 | // 7*(0.0050+0.0060/pt^0.9) | |
1388 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); | |
1389 | ||
1390 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 60 TPCclust"; | |
1391 | } | |
1392 | ||
1393 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 80 clusters | |
1394 | if (cutMode == 90) | |
1395 | { | |
1396 | Int_t minclsTPC=80; | |
1397 | Double_t maxchi2perTPCcl=4.; | |
1398 | Double_t maxdcazITSTPC=1.e9; | |
1399 | ||
1400 | // | |
1401 | // TPC | |
1402 | // | |
1403 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1404 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1405 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1406 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1407 | // | |
1408 | // ITS | |
1409 | // | |
1410 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1411 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1412 | // | |
1413 | // primary selection | |
1414 | // | |
1415 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1416 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1417 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1418 | ||
1419 | // 7*(0.0050+0.0060/pt^0.9) | |
1420 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); | |
1421 | ||
1422 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 80 TPCclust"; | |
1423 | } | |
1424 | ||
1425 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) + TPCchi2=3.5 | |
1426 | if (cutMode == 95) | |
1427 | { | |
1428 | Int_t minclsTPC=80; | |
1429 | Double_t maxchi2perTPCcl=3.5; | |
1430 | Double_t maxdcazITSTPC=1.e9; | |
1431 | ||
1432 | // | |
1433 | // TPC | |
1434 | // | |
1435 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1436 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1437 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1438 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1439 | // | |
1440 | // ITS | |
1441 | // | |
1442 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1443 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1444 | // | |
1445 | // primary selection | |
1446 | // | |
1447 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1448 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1449 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1450 | ||
1451 | // 7*(0.0050+0.0060/pt^0.9) | |
1452 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); | |
1453 | ||
1454 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + TPCchi2 3.5"; | |
1455 | } | |
1456 | ||
1457 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) + TPCchi2=4.5 | |
1458 | if (cutMode == 100) | |
1459 | { | |
1460 | Int_t minclsTPC=80; | |
1461 | Double_t maxchi2perTPCcl=4.5; | |
1462 | Double_t maxdcazITSTPC=1.e9; | |
1463 | ||
1464 | // | |
1465 | // TPC | |
1466 | // | |
1467 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1468 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1469 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1470 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1471 | // | |
1472 | // ITS | |
1473 | // | |
1474 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1475 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1476 | // | |
1477 | // primary selection | |
1478 | // | |
1479 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1480 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1481 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1482 | ||
1483 | // 7*(0.0050+0.0060/pt^0.9) | |
1484 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); | |
1485 | ||
1486 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + TPCchi2 4.5"; | |
1487 | } | |
1488 | ||
1489 | // TPC-tracks | |
1490 | if (cutMode == 110) | |
1491 | { | |
1492 | ||
1493 | minNClustersTPC = 70; | |
1494 | maxChi2PerClusterTPC = 4.0; | |
1495 | maxDCAtoVertexXY = 1.e9; // cm | |
1496 | maxDCAtoVertexZ = 1.e9; // cm | |
1497 | ||
1498 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1499 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
1500 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1501 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1502 | esdTrackCuts->SetMinNClustersTPC(minNClustersTPC); | |
1503 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxChi2PerClusterTPC); | |
1504 | esdTrackCuts->SetMaxDCAToVertexXY(maxDCAtoVertexXY); | |
1505 | esdTrackCuts->SetMaxDCAToVertexZ(maxDCAtoVertexZ); | |
1506 | esdTrackCuts->SetDCAToVertex2D(kTRUE); | |
1507 | ||
1508 | TString tag = "TPC-tracks loose criteria"; | |
1509 | } | |
791aaf54 | 1510 | |
7a761ceb | 1511 | |
1512 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 50 TPCclust | |
1513 | if (cutMode == 120) | |
1514 | { | |
1515 | Int_t minclsTPC=50; | |
1516 | Double_t maxchi2perTPCcl=4.; | |
1517 | Double_t maxdcazITSTPC=1.e9; | |
1518 | ||
1519 | // | |
1520 | // TPC | |
1521 | // | |
1522 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1523 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1524 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1525 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1526 | // | |
1527 | // ITS | |
1528 | // | |
1529 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1530 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1531 | // | |
1532 | // primary selection | |
1533 | // | |
1534 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1535 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1536 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1537 | ||
1538 | // 7*(0.0050+0.0060/pt^0.9) | |
1539 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); | |
1540 | ||
1541 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 60 TPCclust"; | |
1542 | } | |
1543 | ||
1544 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 70 TPCclust + accept kink daughters | |
1545 | if (cutMode == 130) | |
1546 | { | |
1547 | Int_t minclsTPC=70; | |
1548 | Double_t maxchi2perTPCcl=4.; | |
1549 | Double_t maxdcazITSTPC=1.e9; | |
1550 | ||
1551 | // | |
1552 | // TPC | |
1553 | // | |
1554 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1555 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
1556 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1557 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1558 | // | |
1559 | // ITS | |
1560 | // | |
1561 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1562 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1563 | // | |
1564 | // primary selection | |
1565 | // | |
1566 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1567 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1568 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1569 | ||
1570 | // 7*(0.0050+0.0060/pt^0.9) | |
1571 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); | |
1572 | ||
1573 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 60 TPCclust"; | |
1574 | } | |
1575 | ||
1576 | // TPC-tracks + SPD point + ITS refit + DCAr(pt) + 30 TPCclust + accept kink daughters | |
1577 | if (cutMode == 140) | |
1578 | { | |
1579 | Int_t minclsTPC=30; | |
1580 | Double_t maxchi2perTPCcl=4.; | |
1581 | Double_t maxdcazITSTPC=1.e9; | |
1582 | ||
1583 | // | |
1584 | // TPC | |
1585 | // | |
1586 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1587 | esdTrackCuts->SetAcceptKinkDaughters(kTRUE); | |
1588 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1589 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1590 | // | |
1591 | // ITS | |
1592 | // | |
1593 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1594 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1595 | // | |
1596 | // primary selection | |
1597 | // | |
1598 | esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1599 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1600 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1601 | ||
1602 | // 7*(0.0050+0.0060/pt^0.9) | |
1603 | esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); | |
1604 | ||
1605 | TString tag = "TPC-tracks + ITS refit + >1 SPD cluster + DCAr(Pt) + 60 TPCclust"; | |
1606 | } | |
1607 | ||
1608 | // Adam Kisiel track selectiion | |
1609 | if (cutMode == 150) | |
1610 | { | |
1611 | Int_t minclsTPC=70; | |
1612 | Double_t maxchi2perTPCcl=4.; | |
1613 | Double_t maxdcazITSTPC=0.25; | |
1614 | Double_t maxdcaxyITSTPC=0.2; | |
1615 | ||
1616 | // | |
1617 | // TPC | |
1618 | // | |
1619 | //esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1620 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1621 | esdTrackCuts->SetMinNClustersTPC(minclsTPC); | |
1622 | esdTrackCuts->SetMaxChi2PerClusterTPC(maxchi2perTPCcl); | |
1623 | // | |
1624 | // ITS | |
1625 | // | |
1626 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1627 | //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1628 | // | |
1629 | // primary selection | |
1630 | // | |
1631 | //esdTrackCuts->SetDCAToVertex2D(kFALSE); | |
1632 | esdTrackCuts->SetRequireSigmaToVertex(kFALSE); | |
1633 | esdTrackCuts->SetMaxDCAToVertexZ(maxdcazITSTPC); | |
1634 | esdTrackCuts->SetMaxDCAToVertexXY(maxdcaxyITSTPC); | |
1635 | ||
1636 | // 7*(0.0050+0.0060/pt^0.9) | |
1637 | //esdTrackCuts->SetMaxDCAToVertexXYPtDep("0.0350+0.0420/pt^0.9"); | |
1638 | ||
1639 | TString tag = "Adam Kisiel track selection"; | |
1640 | } | |
1641 | ||
82292e7f | 1642 | // TPC+ITS refit |
1643 | // for cut studies | |
1644 | if (cutMode == 151) | |
1645 | { | |
1646 | // | |
1647 | // TPC | |
1648 | // | |
1649 | esdTrackCuts->SetRequireTPCRefit(kTRUE); | |
1650 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1651 | // | |
1652 | // ITS | |
1653 | // | |
1654 | esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1655 | esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1656 | // | |
1657 | ||
1658 | TString tag = "TPC+ITS refit required - for cut studies"; | |
1659 | } | |
1660 | ||
1661 | // TPC+ITS | |
1662 | // for cut studies | |
1663 | if (cutMode == 152) | |
1664 | { | |
1665 | // | |
1666 | // TPC | |
1667 | // | |
1668 | esdTrackCuts->SetRequireTPCRefit(kTRUE); | |
1669 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1670 | // | |
1671 | // ITS | |
1672 | // | |
1673 | //esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1674 | //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1675 | // | |
1676 | ||
1677 | TString tag = "TPC refit required - for cut studies"; | |
1678 | } | |
1679 | ||
d1c6e451 | 1680 | // TPC |
1681 | // for cut studies | |
1682 | if (cutMode == 153) | |
1683 | { | |
1684 | // | |
1685 | // TPC | |
1686 | // | |
1687 | esdTrackCuts->SetRequireTPCRefit(kFALSE); | |
1688 | esdTrackCuts->SetRequireITSRefit(kFALSE); | |
1689 | esdTrackCuts->SetRequireTPCStandAlone(kTRUE); | |
1690 | esdTrackCuts->SetAcceptKinkDaughters(kFALSE); | |
1691 | // | |
1692 | // ITS | |
1693 | // | |
1694 | //esdTrackCuts->SetRequireITSRefit(kTRUE); | |
1695 | //esdTrackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD,AliESDtrackCuts::kAny); | |
1696 | // | |
1697 | ||
1698 | TString tag = "TPC stand alone - for cut studies"; | |
1699 | } | |
1700 | ||
1701 | ||
82292e7f | 1702 | |
7a761ceb | 1703 | |
1704 | ||
1705 | ||
4fa29112 | 1706 | // cuts for data without field |
1707 | if (!fieldOn) | |
1708 | { | |
1709 | cov5 = 1e10; | |
1710 | tag += " without field"; | |
1711 | } | |
1712 | ||
1713 | Printf("Created track cuts for: %s", tag.Data()); | |
1714 | ||
1715 | return esdTrackCuts; | |
1716 | } |