PWG0/dNdPt/macros/CreateStandardCuts.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* CreateTrackCuts(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 }