PWG4/macros/CreateTrackCutsPWG4.C

   1 AliESDtrackCuts *CreateTrackCutsPWG4(Int_t cutMode) {
   2
   3   //
   4   // Macro to create track cuts for PWG4 Jet analysis
   5   // User can select a specific set by indicating cutMode
   6   // cutMode has 8 digits: first 4 digits additional cuts, last 4 digits standard cuts
   7   //                       additional cuts are variations of standard cuts (used for hybrid track selection and QA)
   8   // Numbering starts from 1000 for standard and additional cut numbers
   9
  10   AliESDtrackCuts *trackCuts  = new AliESDtrackCuts("AliESDtrackCuts");
  11
  12   TString tag;
  13
  14   Int_t mod = 10000;
  15
  16   Bool_t bStdCutsDefined = kFALSE;
  17
  18
  19   //_____________________________________________________________________
  20   //                     STANDARD CUTS
  21
  22   //Get standard cuts: last 4 digits of cutMode
  23   Int_t stdCutMode = cutMode%mod;
  24
  25   if(stdCutMode == 1000) {
  26
  27     bStdCutsDefined = kTRUE;
  28
  29     trackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kTRUE,1);
  30     trackCuts->SetMinNCrossedRowsTPC(120);
  31     trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.8);
  32     trackCuts->SetMaxChi2PerClusterITS(36);
  33     trackCuts->SetMaxFractionSharedTPCClusters(0.4);
  34
  35     trackCuts->SetEtaRange(-0.9,0.9);
  36     trackCuts->SetPtRange(0.15, 1e10);
  37
  38     tag = "Global track RAA analysis QM2011 + Chi2ITS<36";
  39
  40   }
  41
  42   if(stdCutMode == 1001) {
  43
  44     bStdCutsDefined = kTRUE;
  45
  46     // TPC
  47     trackCuts->SetMinNClustersTPC(90);
  48     trackCuts->SetMaxChi2PerClusterTPC(4);
  49     trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
  50     trackCuts->SetAcceptKinkDaughters(kFALSE);
  51     trackCuts->SetRequireTPCRefit(kTRUE);
  52     trackCuts->SetMaxFractionSharedTPCClusters(0.4);
  53     // ITS
  54     trackCuts->SetRequireITSRefit(kTRUE);
  55     //accept secondaries
  56     trackCuts->SetMaxDCAToVertexXY(2.4);
  57     trackCuts->SetMaxDCAToVertexZ(3.2);
  58     trackCuts->SetDCAToVertex2D(kTRUE);
  59     //reject fakes
  60     trackCuts->SetMaxChi2PerClusterITS(36);
  61
  62     trackCuts->SetRequireSigmaToVertex(kFALSE);
  63
  64     trackCuts->SetEtaRange(-0.9,0.9);
  65     trackCuts->SetPtRange(0.15, 100.);
  66
  67
  68     tag = "Global tracks jet analysis with ITSrefit and NclsIter1=90, noSPD requirement";
  69
  70   }
  71
  72   if(stdCutMode == 1002) {
  73
  74     bStdCutsDefined = kTRUE;
  75
  76     // TPC
  77     trackCuts->SetMinNClustersTPC(80);
  78     trackCuts->SetMaxChi2PerClusterTPC(4);
  79     trackCuts->SetAcceptKinkDaughters(kFALSE);
  80     trackCuts->SetRequireTPCRefit(kTRUE);
  81     trackCuts->SetMaxFractionSharedTPCClusters(0.4);
  82     // ITS
  83     trackCuts->SetRequireITSRefit(kTRUE);
  84     //accept secondaries
  85     trackCuts->SetMaxDCAToVertexXY(2.4);
  86     trackCuts->SetMaxDCAToVertexZ(3.2);
  87     trackCuts->SetDCAToVertex2D(kTRUE);
  88     //reject fakes
  89     trackCuts->SetMaxChi2PerClusterITS(36);
  90
  91     trackCuts->SetRequireSigmaToVertex(kFALSE);
  92
  93     trackCuts->SetEtaRange(-0.9,0.9);
  94     trackCuts->SetPtRange(0.15, 100.);
  95
  96
  97     tag = "Global tracks jet analysis with ITSrefit and Ncls=80, noSPD requirement";
  98
  99   }
 100
 101   if(stdCutMode == 1003) {
 102
 103     bStdCutsDefined = kTRUE;
 104
 105     // tight global tracks
 106     trackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kFALSE,1);
 107     trackCuts->SetMinNClustersTPC(0);
 108     trackCuts->SetMinNCrossedRowsTPC(120);
 109     trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.1);// essentially switches it off
 110     trackCuts->SetMaxDCAToVertexXY(2.4);
 111     trackCuts->SetMaxDCAToVertexZ(3.2);
 112     trackCuts->SetDCAToVertex2D(kTRUE);
 113     trackCuts->SetMaxChi2PerClusterITS(36);
 114     trackCuts->SetMaxFractionSharedTPCClusters(0.4);
 115
 116     tag = "Global tracks ITSTPC2010 + NCrossedRows + loose ITS";
 117
 118   }
 119
 120   if(stdCutMode == 2000) {
 121
 122     bStdCutsDefined = kTRUE;
 123
 124     // TPC
 125     trackCuts->SetMinNClustersTPC(90);
 126     trackCuts->SetMaxChi2PerClusterTPC(4);
 127     trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
 128     trackCuts->SetAcceptKinkDaughters(kFALSE);
 129     trackCuts->SetRequireTPCRefit(kTRUE);
 130     trackCuts->SetMaxFractionSharedTPCClusters(0.4);
 131     //accept secondaries
 132     trackCuts->SetMaxDCAToVertexXY(2.4);
 133     trackCuts->SetMaxDCAToVertexZ(3.2);
 134     trackCuts->SetDCAToVertex2D(kTRUE);
 135
 136     trackCuts->SetRequireSigmaToVertex(kFALSE);
 137
 138     trackCuts->SetEtaRange(-0.9,0.9);
 139     trackCuts->SetPtRange(0.15, 100.);
 140
 141
 142     tag = "Global tracks jet analysis, loose cuts, NClsIter1=90, no ITS requirements";
 143
 144   }
 145
 146   if(stdCutMode == 2001) {
 147
 148     bStdCutsDefined = kTRUE;
 149
 150     trackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
 151     trackCuts->SetMinNClustersTPC(70);
 152
 153     trackCuts->SetEtaRange(-0.9,0.9);
 154     trackCuts->SetPtRange(0.15, 100.);
 155
 156
 157     tag = "TPConly track cuts, loose cuts, NCls=70, no ITS requirements";
 158
 159   }
 160
 161   if(stdCutMode == 2002) {
 162
 163     bStdCutsDefined = kTRUE;
 164
 165     trackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
 166     trackCuts->SetMinNClustersTPC(0);
 167     trackCuts->SetMinNCrossedRowsTPC(120);
 168     trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.1);// essentially switches it off
 169
 170     trackCuts->SetEtaRange(-0.9,0.9);
 171     trackCuts->SetPtRange(0.15, 100.);
 172
 173
 174     tag = "TPConly track cuts, loose cuts, NCls=70, no ITS requirements";
 175
 176   }
 177
 178   if(!bStdCutsDefined) {
 179     printf("last 4 digits do not represent a predefined set of standard cuts. Returning 0\n");
 180     return 0;
 181
 182   }
 183
 184
 185   //_____________________________________________________________________
 186   //                     ADDITIONAL CUTS
 187
 188   //Get additional cut mode: first 4 digits of cutMode
 189   Int_t addCutMode = (int)((float)cutMode/(float)mod);
 190
 191   if(addCutMode == 1000) {
 192
 193     trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
 194
 195     tag += " + additonal: SPD any requirement";
 196
 197   }
 198
 199   if(addCutMode == 1001) {
 200
 201     trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);
 202
 203     tag += " + additional: w/o hits in SPD";
 204
 205   }
 206
 207   if(addCutMode == 1002) {
 208
 209     trackCuts->SetMaxChi2PerClusterITS(1E10);
 210
 211     tag += " + additional: maxITSChi2=1e10";
 212
 213   }
 214
 215   if(addCutMode == 1003) {
 216
 217     trackCuts->SetMinNClustersTPC(0);
 218     trackCuts->SetMinNCrossedRowsTPC(0);
 219     trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.);
 220
 221     tag += " + additional: minClusters=0 minCrossedRows=0 minCrossedRowsOverFindable=0";
 222
 223   }
 224
 225   if(addCutMode == 1004) {
 226
 227     trackCuts->SetRequireITSRefit(kFALSE);
 228
 229     tag += " + additional: ITSrefit=kFALSE";
 230
 231   }
 232
 233   Printf("Created track cuts for: %s", tag.Data());
 234
 235   return trackCuts;
 236
 237 }