1 AliESDtrackCuts *CreateTrackCutsPWGJE(Int_t cutMode) {
4 // Macro to create track cuts for PWG 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
10 AliESDtrackCuts *trackCuts = new AliESDtrackCuts("AliESDtrackCuts");
16 Bool_t bStdCutsDefined = kFALSE;
19 //_____________________________________________________________________
22 //Get standard cuts: last 4 digits of cutMode
23 Int_t stdCutMode = cutMode%mod;
25 if(stdCutMode == 1000) {
27 bStdCutsDefined = kTRUE;
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 trackCuts->SetMaxChi2TPCConstrainedGlobal(36);
36 trackCuts->SetEtaRange(-0.9,0.9);
37 trackCuts->SetPtRange(0.15, 1e10);
39 tag = "Global track RAA analysis QM2011 + Chi2ITS<36";
43 if(stdCutMode == 1001) {
45 bStdCutsDefined = kTRUE;
48 trackCuts->SetMinNClustersTPC(90);
49 trackCuts->SetMaxChi2PerClusterTPC(4);
50 trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
51 trackCuts->SetAcceptKinkDaughters(kFALSE);
52 trackCuts->SetRequireTPCRefit(kTRUE);
53 trackCuts->SetMaxFractionSharedTPCClusters(0.4);
55 trackCuts->SetRequireITSRefit(kTRUE);
57 trackCuts->SetMaxDCAToVertexXY(2.4);
58 trackCuts->SetMaxDCAToVertexZ(3.2);
59 trackCuts->SetDCAToVertex2D(kTRUE);
61 trackCuts->SetMaxChi2PerClusterITS(36);
63 trackCuts->SetRequireSigmaToVertex(kFALSE);
65 trackCuts->SetEtaRange(-0.9,0.9);
66 trackCuts->SetPtRange(0.15, 100.);
69 tag = "Global tracks jet analysis with ITSrefit and NclsIter1=90, noSPD requirement";
73 if(stdCutMode == 1002) {
75 bStdCutsDefined = kTRUE;
78 trackCuts->SetMinNClustersTPC(80);
79 trackCuts->SetMaxChi2PerClusterTPC(4);
80 trackCuts->SetAcceptKinkDaughters(kFALSE);
81 trackCuts->SetRequireTPCRefit(kTRUE);
82 trackCuts->SetMaxFractionSharedTPCClusters(0.4);
84 trackCuts->SetRequireITSRefit(kTRUE);
86 trackCuts->SetMaxDCAToVertexXY(2.4);
87 trackCuts->SetMaxDCAToVertexZ(3.2);
88 trackCuts->SetDCAToVertex2D(kTRUE);
90 trackCuts->SetMaxChi2PerClusterITS(36);
92 trackCuts->SetRequireSigmaToVertex(kFALSE);
94 trackCuts->SetEtaRange(-0.9,0.9);
95 trackCuts->SetPtRange(0.15, 100.);
98 tag = "Global tracks jet analysis with ITSrefit and Ncls=80, noSPD requirement";
102 if(stdCutMode == 1003) {
104 bStdCutsDefined = kTRUE;
106 // tight global tracks
107 trackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010(kFALSE,1);
108 trackCuts->SetMinNClustersTPC(0);
109 trackCuts->SetMinNCrossedRowsTPC(120);
110 trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.1);// essentially switches it off
111 trackCuts->SetMaxDCAToVertexXY(2.4);
112 trackCuts->SetMaxDCAToVertexZ(3.2);
113 trackCuts->SetDCAToVertex2D(kTRUE);
114 trackCuts->SetMaxChi2PerClusterITS(36);
115 trackCuts->SetMaxFractionSharedTPCClusters(0.4);
117 tag = "Global tracks ITSTPC2010 + NCrossedRows + loose ITS";
121 if(stdCutMode == 1004) {
123 bStdCutsDefined = kTRUE;
126 trackCuts->SetMinNClustersTPC(70);
127 trackCuts->SetMaxChi2PerClusterTPC(4);
128 trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
129 trackCuts->SetAcceptKinkDaughters(kFALSE);
130 trackCuts->SetRequireTPCRefit(kTRUE);
131 trackCuts->SetMaxFractionSharedTPCClusters(0.4);
133 trackCuts->SetRequireITSRefit(kTRUE);
135 trackCuts->SetMaxDCAToVertexXY(2.4);
136 trackCuts->SetMaxDCAToVertexZ(3.2);
137 trackCuts->SetDCAToVertex2D(kTRUE);
139 trackCuts->SetMaxChi2PerClusterITS(36);
141 trackCuts->SetRequireSigmaToVertex(kFALSE);
143 trackCuts->SetEtaRange(-0.9,0.9);
144 trackCuts->SetPtRange(0.15, 100.);
146 tag = "Global tracks jet analysis with ITSrefit and NclsIter1=70, noSPD requirement";
149 if(stdCutMode == 1005) {
151 bStdCutsDefined = kTRUE;
154 trackCuts->SetMinNClustersTPC(70);
155 trackCuts->SetMaxChi2PerClusterTPC(4);
156 trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
157 trackCuts->SetAcceptKinkDaughters(kFALSE);
158 trackCuts->SetRequireTPCRefit(kTRUE);
159 trackCuts->SetMaxFractionSharedTPCClusters(0.4);
161 trackCuts->SetRequireITSRefit(kTRUE);
163 trackCuts->SetMaxDCAToVertexXY(2.4);
164 trackCuts->SetMaxDCAToVertexZ(3.2);
165 trackCuts->SetDCAToVertex2D(kTRUE);
167 trackCuts->SetMaxChi2PerClusterITS(36);
169 trackCuts->SetRequireSigmaToVertex(kFALSE);
171 trackCuts->SetEtaRange(-0.9,0.9);
172 trackCuts->SetPtRange(0.15, 1E+15.);
174 tag = "Global tracks jet analysis with ITSrefit and NclsIter1=70, noSPD requirement, no upper pt cut";
177 if(stdCutMode == 1006) {
179 bStdCutsDefined = kTRUE;
182 TFormula *f1NClustersTPCLinearPtDep = new TFormula("f1NClustersTPCLinearPtDep","70.+30./20.*x");
183 trackCuts->SetMinNClustersTPCPtDep(f1NClustersTPCLinearPtDep,20.);
184 trackCuts->SetMinNClustersTPC(70);
185 trackCuts->SetMaxChi2PerClusterTPC(4);
186 trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
187 trackCuts->SetAcceptKinkDaughters(kFALSE);
188 trackCuts->SetRequireTPCRefit(kTRUE);
189 trackCuts->SetMaxFractionSharedTPCClusters(0.4);
191 trackCuts->SetRequireITSRefit(kTRUE);
193 trackCuts->SetMaxDCAToVertexXY(2.4);
194 trackCuts->SetMaxDCAToVertexZ(3.2);
195 trackCuts->SetDCAToVertex2D(kTRUE);
197 trackCuts->SetMaxChi2PerClusterITS(36);
198 trackCuts->SetMaxChi2TPCConstrainedGlobal(36);
200 trackCuts->SetRequireSigmaToVertex(kFALSE);
202 trackCuts->SetEtaRange(-0.9,0.9);
203 trackCuts->SetPtRange(0.15, 1E+15.);
205 tag = "Global tracks jet analysis with ITSrefit and NclsIter1=PtDep, noSPD requirement, no upper pt cut, golden chi2";
211 if(stdCutMode == 2000) {
213 bStdCutsDefined = kTRUE;
216 trackCuts->SetMinNClustersTPC(90);
217 trackCuts->SetMaxChi2PerClusterTPC(4);
218 trackCuts->SetRequireTPCStandAlone(kTRUE); //cut on NClustersTPC and chi2TPC Iter1
219 trackCuts->SetAcceptKinkDaughters(kFALSE);
220 trackCuts->SetRequireTPCRefit(kTRUE);
221 trackCuts->SetMaxFractionSharedTPCClusters(0.4);
223 trackCuts->SetMaxDCAToVertexXY(2.4);
224 trackCuts->SetMaxDCAToVertexZ(3.2);
225 trackCuts->SetDCAToVertex2D(kTRUE);
227 trackCuts->SetRequireSigmaToVertex(kFALSE);
229 trackCuts->SetEtaRange(-0.9,0.9);
230 trackCuts->SetPtRange(0.15, 100.);
233 tag = "Global tracks jet analysis, loose cuts, NClsIter1=90, no ITS requirements";
237 if(stdCutMode == 2001) {
239 bStdCutsDefined = kTRUE;
241 trackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
242 trackCuts->SetRequireTPCRefit(kTRUE);
243 trackCuts->SetMinNClustersTPC(70);
245 trackCuts->SetEtaRange(-0.9,0.9);
246 trackCuts->SetPtRange(0.15, 100.);
249 tag = "TPConly track cuts, loose cuts, NCls=70, no ITS requirements";
253 if(stdCutMode == 2002) {
255 bStdCutsDefined = kTRUE;
257 trackCuts = AliESDtrackCuts::GetStandardTPCOnlyTrackCuts();
258 trackCuts->SetRequireTPCRefit(kTRUE);
259 trackCuts->SetMinNClustersTPC(0);
260 trackCuts->SetMinNCrossedRowsTPC(120);
261 trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.1);// essentially switches it off
263 trackCuts->SetEtaRange(-0.9,0.9);
264 trackCuts->SetPtRange(0.15, 100.);
267 tag = "TPConly track cuts, loose cuts, NCls=70, no ITS requirements";
271 if(!bStdCutsDefined) {
272 printf("last 4 digits do not represent a predefined set of standard cuts. Returning 0\n");
278 //_____________________________________________________________________
281 //Get additional cut mode: first 4 digits of cutMode
282 Int_t addCutMode = (int)((float)cutMode/(float)mod);
284 if(addCutMode == 1000) {
286 trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kAny);
288 tag += " + additonal: SPD any requirement";
292 if(addCutMode == 1001) {
294 trackCuts->SetClusterRequirementITS(AliESDtrackCuts::kSPD, AliESDtrackCuts::kNone);
296 tag += " + additional: w/o hits in SPD";
300 if(addCutMode == 1002) {
302 trackCuts->SetMaxChi2PerClusterITS(1E10);
304 tag += " + additional: maxITSChi2=1e10";
308 if(addCutMode == 1003) {
310 trackCuts->SetMinNClustersTPC(0);
311 trackCuts->SetMinNCrossedRowsTPC(0);
312 trackCuts->SetMinRatioCrossedRowsOverFindableClustersTPC(0.);
314 tag += " + additional: minClusters=0 minCrossedRows=0 minCrossedRowsOverFindable=0";
318 if(addCutMode == 1004) {
320 trackCuts->SetRequireITSRefit(kFALSE);
322 tag += " + additional: ITSrefit=kFALSE";
326 Printf("Created track cuts for: %s", tag.Data());