]>
Commit | Line | Data |
---|---|---|
1ee39b3a | 1 | /************************************************************************** |
2 | * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * | |
3 | * * | |
4 | * Author: The ALICE Off-line Project. * | |
5 | * Contributors are mentioned in the code where appropriate. * | |
6 | * * | |
7 | * Permission to use, copy, modify and distribute this software and its * | |
8 | * documentation strictly for non-commercial purposes is hereby granted * | |
9 | * without fee, provided that the above copyright notice appears in all * | |
10 | * copies and that both the copyright notice and this permission notice * | |
11 | * appear in the supporting documentation. The authors make no claims * | |
12 | * about the suitability of this software for any purpose. It is * | |
13 | * provided "as is" without express or implied warranty. * | |
14 | **************************************************************************/ | |
15 | ||
16 | /* $Id: AliTRDv0Info.cxx 27496 2008-07-22 08:35:45Z cblume $ */ | |
17 | ||
18 | //////////////////////////////////////////////////////////////////////////// | |
19 | // // | |
20 | // Reconstruction QA // | |
21 | // // | |
22 | // Gathers all information necessary for reference data selection about // | |
23 | // the track and (in case) its corresponding V0. // | |
24 | // Carries out the selection of electrons (from gamma conversions), // | |
25 | // pions (from K0s decays) and protons (from Lambda and Anti-Lambda // | |
26 | // decays) by cuts specific for the respective decay and particle // | |
27 | // species. // | |
28 | // (M.Heide, 2009/10/06) // | |
29 | // // | |
30 | // Authors: // | |
31 | // Alex Bercuci <A.Bercuci@gsi.de> // | |
32 | // Alex Wilk <wilka@uni-muenster.de> // | |
33 | // Markus Heide <mheide@uni-muenster.de> // | |
34 | // // | |
35 | //////////////////////////////////////////////////////////////////////////// | |
0ed6f095 | 36 | |
1ee39b3a | 37 | #include "TMath.h" |
0ed6f095 | 38 | #include "TDatabasePDG.h" |
1ee39b3a | 39 | |
40 | #include "AliESDtrack.h" | |
41 | #include "AliESDv0.h" | |
1ee39b3a | 42 | #include "AliLog.h" |
0ed6f095 | 43 | #include "TVector3.h" |
44 | #include "AliKFParticle.h" | |
45 | #include "AliKFVertex.h" | |
1ee39b3a | 46 | |
47 | #include "AliTRDv0Info.h" | |
48 | #include "AliTRDtrackInfo.h" | |
49 | #include "AliTRDtrackInfo.h" | |
50 | ||
51 | ClassImp(AliTRDv0Info) | |
52 | ||
53 | //_________________________________________________ | |
54 | AliTRDv0Info::AliTRDv0Info() | |
55 | : TObject() | |
3d19c1b0 | 56 | ,fQuality(0) |
1ee39b3a | 57 | ,fDCA(10) |
58 | ,fPointingAngle(10) | |
59 | ,fOpenAngle(10) | |
60 | ,fPsiPair(99) | |
61 | ,fMagField(0) | |
62 | ,fRadius(0) | |
1ee39b3a | 63 | ,fV0Momentum(0) |
1ee39b3a | 64 | ,fNindex(0) |
65 | ,fPindex(0) | |
0ed6f095 | 66 | ,fInputEvent(NULL) |
67 | ,fPrimaryVertex(NULL) | |
64d57299 | 68 | ,fTrackP(NULL) |
69 | ,fTrackN(NULL) | |
1ee39b3a | 70 | { |
71 | // | |
72 | // Default constructor | |
73 | // | |
74 | ||
1ee39b3a | 75 | memset(fDetPID, 0, 2*kNDaughters*kNDetectors*AliPID::kSPECIES*sizeof(Float_t)); |
8bc8ea55 | 76 | memset(fComPID, 0, 2*kNDaughters*AliPID::kSPECIES*sizeof(Float_t)); |
61cfa442 | 77 | memset(fInvMass, 0, kNDecays*sizeof(Double_t)); |
0ed6f095 | 78 | memset(fArmenteros, 0, kNDecays*sizeof(Bool_t)); |
61cfa442 | 79 | memset(fTPCdEdx, 0, kNDaughters*sizeof(Float_t)); |
0ed6f095 | 80 | memset(fChi2ndf, 0, kNDecays*sizeof(Double_t)); |
1ee39b3a | 81 | |
82 | ///////////////////////////////////////////////////////////////////////////// | |
83 | //Set Cut values: First specify decay in brackets, then the actual cut value! | |
84 | ///////////////////////////////////////////////////////////////////////////// | |
85 | ||
86 | //Upper limit for distance of closest approach of two daughter tracks : | |
8bc8ea55 | 87 | fUpDCA[kGamma] = 1000.; |
88 | fUpDCA[kK0s] = 0.08; | |
89 | fUpDCA[kLambda] = 0.2; | |
90 | fUpDCA[kAntiLambda] = 0.2; | |
1ee39b3a | 91 | |
92 | //Upper limit for pointing angle (= angle between between vector from primary to secondary vertex and reconstructed momentum of V0 mother particle) : | |
93 | fUpPointingAngle[kGamma] = 0.03; | |
94 | fUpPointingAngle[kK0s] = 0.03; | |
8bc8ea55 | 95 | fUpPointingAngle[kLambda] = 0.04; |
96 | fUpPointingAngle[kAntiLambda] = 0.04; | |
1ee39b3a | 97 | |
98 | //Upper limit for invariant mass of V0 mother : | |
8bc8ea55 | 99 | fUpInvMass[kGamma][0] = 0.05;// second pair of brackets is for momentum bin: 0: below mother momentm of 2.5 GeV |
1ee39b3a | 100 | fUpInvMass[kGamma][1] = 0.07;//1: above 2.5 GeV |
8bc8ea55 | 101 | fUpInvMass[kK0s][0] = fUpInvMass[kK0s][1] = 0.50265; |
102 | fUpInvMass[kLambda][0] = fUpInvMass[kLambda][1] = 1.1207; | |
103 | fUpInvMass[kAntiLambda][0] = fUpInvMass[kAntiLambda][1] = 1.1207; | |
1ee39b3a | 104 | |
105 | //Lower limit for invariant mass of V0 mother : | |
106 | fDownInvMass[kGamma] = -1.; | |
8bc8ea55 | 107 | fDownInvMass[kK0s] = 0.49265; |
108 | fDownInvMass[kLambda] = 1.107; | |
109 | fDownInvMass[kAntiLambda] = 1.107; | |
1ee39b3a | 110 | |
0ed6f095 | 111 | //Upper limit for KF Chi2/NDF value; |
112 | fUpChi2ndf[kGamma] = 10000.;//7.; | |
113 | fUpChi2ndf[kK0s] = 10000.;//5.; | |
114 | fUpChi2ndf[kLambda] = 10000.;//5.; | |
115 | fUpChi2ndf[kAntiLambda] = 10000.;//5.; | |
116 | ||
1ee39b3a | 117 | //Lower limit for distance from secondary vertex to primary vertex in x-y plane : |
8bc8ea55 | 118 | fDownRadius[kGamma] = 6.; |
1ee39b3a | 119 | fDownRadius[kK0s] = 0.; |
8bc8ea55 | 120 | fDownRadius[kLambda] = 0.; |
121 | fDownRadius[kAntiLambda] = 0.; | |
1ee39b3a | 122 | |
123 | //Upper limit for distance from secondary vertex to primary vertex in x-y plane : | |
124 | fUpRadius[kGamma] = 1000.; | |
8bc8ea55 | 125 | fUpRadius[kK0s] = 20.; |
1ee39b3a | 126 | fUpRadius[kLambda] = 1000.; |
127 | fUpRadius[kAntiLambda] = 1000.; | |
128 | ||
129 | //Upper limit for opening angle between two daughter tracks (characteristically near zero for conversions) : | |
130 | fUpOpenAngle[kGamma] = 0.1; | |
131 | fUpOpenAngle[kK0s] = 3.15; | |
132 | fUpOpenAngle[kLambda] = 3.15; | |
133 | fUpOpenAngle[kAntiLambda] = 3.15; | |
134 | ||
135 | //Upper limit for angle between daughter momentum plane and plane perpendicular to magnetic field (characteristically around zero for conversions) : | |
8bc8ea55 | 136 | fUpPsiPair[kGamma] = 0.05; |
1ee39b3a | 137 | fUpPsiPair[kK0s] = 1.6; |
138 | fUpPsiPair[kLambda] = 1.6; | |
139 | fUpPsiPair[kAntiLambda] = 1.6; | |
140 | ||
141 | //Lower limit for likelihood value of TPC PID : | |
0ed6f095 | 142 | fDownTPCPIDneg[AliPID::kElectron] = 0.; |
143 | fDownTPCPIDpos[AliPID::kElectron] = 0.; | |
144 | ||
145 | fDownTPCPIDneg[AliPID::kMuon] = 0.; | |
146 | fDownTPCPIDpos[AliPID::kMuon] = 0.; | |
147 | ||
148 | fDownTPCPIDneg[AliPID::kPion] = 0.; | |
149 | fDownTPCPIDpos[AliPID::kPion] = 0.; | |
150 | ||
151 | fDownTPCPIDneg[AliPID::kKaon] = 0.; | |
152 | fDownTPCPIDpos[AliPID::kKaon] = 0.; | |
153 | ||
154 | fDownTPCPIDneg[AliPID::kProton] = 0.; | |
155 | fDownTPCPIDpos[AliPID::kProton] = 0.; | |
156 | ||
157 | //Lower limit for likelihood value of combined PID : | |
158 | fDownComPIDneg[AliPID::kElectron] = 0.; | |
159 | fDownComPIDpos[AliPID::kElectron] = 0.; | |
160 | ||
161 | fDownComPIDneg[AliPID::kMuon] = 0.; | |
162 | fDownComPIDpos[AliPID::kMuon] = 0.; | |
163 | ||
164 | fDownComPIDneg[AliPID::kPion] = 0.; | |
165 | fDownComPIDpos[AliPID::kPion] = 0.; | |
166 | ||
167 | fDownComPIDneg[AliPID::kKaon] = 0.; | |
168 | fDownComPIDpos[AliPID::kKaon] = 0.; | |
169 | ||
170 | fDownComPIDneg[AliPID::kProton] = 0.; | |
171 | fDownComPIDpos[AliPID::kProton] = 0.; | |
172 | ||
173 | //Lower limit for likelihood value of combined PID for daughter track which doesn't enter reference data (here: pion daughters from Lambda decays: | |
174 | fDownComPIDnegPart[AliPID::kElectron] = 0.; | |
175 | fDownComPIDposPart[AliPID::kElectron] = 0.; | |
176 | ||
177 | fDownComPIDnegPart[AliPID::kMuon] = 0.; | |
178 | fDownComPIDposPart[AliPID::kMuon] = 0.; | |
179 | ||
180 | fDownComPIDnegPart[AliPID::kPion] = 0.; | |
181 | fDownComPIDposPart[AliPID::kPion] = 0.; | |
182 | ||
183 | fDownComPIDnegPart[AliPID::kKaon] = 0.; | |
184 | fDownComPIDposPart[AliPID::kKaon] = 0.; | |
185 | ||
186 | fDownComPIDnegPart[AliPID::kProton] = 0.; | |
187 | fDownComPIDposPart[AliPID::kProton] = 0.; | |
188 | ||
189 | //Parameters for data with well-calibrated PID (after usage of tender): | |
190 | /* //Lower limit for likelihood value of TPC PID : | |
1ee39b3a | 191 | fDownTPCPIDneg[AliPID::kElectron] = 0.21; |
192 | fDownTPCPIDpos[AliPID::kElectron] = 0.21; | |
193 | ||
194 | fDownTPCPIDneg[AliPID::kMuon] = 0.21; | |
195 | fDownTPCPIDpos[AliPID::kMuon] = 0.21; | |
196 | ||
197 | fDownTPCPIDneg[AliPID::kPion] = 0.21; | |
198 | fDownTPCPIDpos[AliPID::kPion] = 0.21; | |
199 | ||
200 | fDownTPCPIDneg[AliPID::kKaon] = 0.21; | |
201 | fDownTPCPIDpos[AliPID::kKaon] = 0.21; | |
202 | ||
203 | fDownTPCPIDneg[AliPID::kProton] = 0.21; | |
204 | fDownTPCPIDpos[AliPID::kProton] = 0.21; | |
8bc8ea55 | 205 | |
0ed6f095 | 206 | //Lower limit for likelihood value of combined PID : |
8bc8ea55 | 207 | fDownComPIDneg[AliPID::kElectron] = 0.21; |
208 | fDownComPIDpos[AliPID::kElectron] = 0.21; | |
209 | ||
210 | fDownComPIDneg[AliPID::kMuon] = 0.21; | |
211 | fDownComPIDpos[AliPID::kMuon] = 0.21; | |
212 | ||
213 | fDownComPIDneg[AliPID::kPion] = 0.9; | |
214 | fDownComPIDpos[AliPID::kPion] = 0.9; | |
215 | ||
216 | fDownComPIDneg[AliPID::kKaon] = 0.21; | |
217 | fDownComPIDpos[AliPID::kKaon] = 0.21; | |
218 | ||
219 | fDownComPIDneg[AliPID::kProton] = 0.9; | |
220 | fDownComPIDpos[AliPID::kProton] = 0.9; | |
221 | ||
222 | //Lower limit for likelihood value of combined PID for daughter track which doesn't enter reference data (here: pion daughters from Lambda decays: | |
223 | fDownComPIDnegPart[AliPID::kElectron] = 0.05; | |
224 | fDownComPIDposPart[AliPID::kElectron] = 0.05; | |
225 | ||
226 | fDownComPIDnegPart[AliPID::kMuon] = 0.05; | |
227 | fDownComPIDposPart[AliPID::kMuon] = 0.05; | |
228 | ||
229 | fDownComPIDnegPart[AliPID::kPion] = 0.05; | |
230 | fDownComPIDposPart[AliPID::kPion] = 0.05; | |
231 | ||
232 | fDownComPIDnegPart[AliPID::kKaon] = 0.05; | |
233 | fDownComPIDposPart[AliPID::kKaon] = 0.05; | |
234 | ||
235 | fDownComPIDnegPart[AliPID::kProton] = 0.05; | |
0ed6f095 | 236 | fDownComPIDposPart[AliPID::kProton] = 0.05;*/ |
3d19c1b0 | 237 | } |
238 | ||
239 | //_________________________________________________ | |
240 | AliTRDv0Info::AliTRDv0Info(const AliTRDv0Info &ref) | |
241 | : TObject() | |
242 | ,fQuality(ref.fQuality) | |
243 | ,fDCA(ref.fDCA) | |
244 | ,fPointingAngle(ref.fPointingAngle) | |
245 | ,fOpenAngle(ref.fOpenAngle) | |
246 | ,fPsiPair(ref.fPsiPair) | |
247 | ,fMagField(ref.fMagField) | |
0ed6f095 | 248 | ,fRadius(ref.fRadius) |
3d19c1b0 | 249 | ,fV0Momentum(ref.fV0Momentum) |
3d19c1b0 | 250 | ,fNindex(ref.fNindex) |
251 | ,fPindex(ref.fPindex) | |
0ed6f095 | 252 | ,fInputEvent(ref.fInputEvent) |
253 | ,fPrimaryVertex(ref.fPrimaryVertex) | |
64d57299 | 254 | ,fTrackP(ref.fTrackP) |
255 | ,fTrackN(ref.fTrackN) | |
3d19c1b0 | 256 | { |
257 | // | |
258 | // Copy constructor | |
259 | // | |
64d57299 | 260 | |
3d19c1b0 | 261 | memcpy(fDetPID, ref.fDetPID, 2*kNDaughters*kNDetectors*AliPID::kSPECIES*sizeof(Float_t)); |
262 | memcpy(fComPID, ref.fComPID, 2*kNDaughters*AliPID::kSPECIES*sizeof(Float_t)); | |
61cfa442 | 263 | memcpy(fInvMass, ref.fInvMass, kNDecays*sizeof(Double_t)); |
0ed6f095 | 264 | memcpy(fArmenteros, ref.fArmenteros, kNDecays*sizeof(Bool_t)); |
265 | memcpy(fChi2ndf, ref.fChi2ndf, kNDecays*sizeof(Double_t)); | |
266 | memcpy(fTPCdEdx, ref.fTPCdEdx, kNDaughters*sizeof(Float_t)); | |
1ee39b3a | 267 | |
3d19c1b0 | 268 | //Upper limit for distance of closest approach of two daughter tracks : |
269 | memcpy(fUpDCA, ref.fUpDCA, kNDecays*sizeof(Float_t)); | |
270 | memcpy(fUpPointingAngle, ref.fUpPointingAngle, kNDecays*sizeof(Float_t)); | |
271 | memcpy(fUpOpenAngle, ref.fUpOpenAngle, kNDecays*sizeof(Float_t)); | |
272 | memcpy(fDownOpenAngle, ref.fDownOpenAngle, kNDecays*sizeof(Float_t)); | |
273 | memcpy(fUpPsiPair, ref.fUpPsiPair, kNDecays*sizeof(Float_t)); | |
274 | memcpy(fDownPsiPair, ref.fDownPsiPair, kNDecays*sizeof(Float_t)); | |
275 | memcpy(fUpInvMass, ref.fUpInvMass, kNDecays*kNMomBins*sizeof(Double_t)); | |
276 | memcpy(fDownInvMass, ref.fDownInvMass, kNDecays*sizeof(Double_t)); | |
0ed6f095 | 277 | memcpy(fUpChi2ndf, ref.fUpChi2ndf, kNDecays*sizeof(Double_t)); |
3d19c1b0 | 278 | memcpy(fUpRadius, ref.fUpRadius, kNDecays*sizeof(Float_t)); |
279 | memcpy(fDownRadius, ref.fDownRadius, kNDecays*sizeof(Float_t)); | |
280 | memcpy(fDownTPCPIDneg, ref.fDownTPCPIDneg, AliPID::kSPECIES*sizeof(Float_t)); | |
281 | memcpy(fDownTPCPIDpos, ref.fDownTPCPIDpos, AliPID::kSPECIES*sizeof(Float_t)); | |
282 | memcpy(fDownComPIDneg, ref.fDownComPIDneg, AliPID::kSPECIES*sizeof(Float_t)); | |
283 | memcpy(fDownComPIDpos, ref.fDownComPIDpos, AliPID::kSPECIES*sizeof(Float_t)); | |
284 | memcpy(fDownComPIDnegPart, ref.fDownComPIDnegPart, AliPID::kSPECIES*sizeof(Float_t)); | |
285 | memcpy(fDownComPIDposPart, ref.fDownComPIDposPart, AliPID::kSPECIES*sizeof(Float_t)); | |
1ee39b3a | 286 | } |
287 | ||
288 | //_________________________________________________ | |
b37d601d | 289 | void AliTRDv0Info::SetV0Info(const AliESDv0 *esdv0) |
64d57299 | 290 | { |
291 | //Gets values of ESDv0 and daughter track properties | |
1ee39b3a | 292 | //See header file for description of variables |
293 | ||
1ee39b3a | 294 | fQuality = Quality(esdv0);//Attributes an Int_t to the V0 due to quality cuts (= 1 if V0 is accepted, other integers depending on cut which excludes the vertex) |
295 | ||
296 | fRadius = Radius(esdv0);//distance from secondary vertex to primary vertex in x-y plane | |
297 | ||
298 | fDCA = esdv0->GetDcaV0Daughters();//distance of closest approach of two daughter tracks | |
299 | ||
300 | fPointingAngle = TMath::ACos(esdv0->GetV0CosineOfPointingAngle());// pointing angle (= angle between between vector from primary to secondary vertex and reconstructed momentum of V0 mother particle) | |
301 | ||
302 | fOpenAngle = OpenAngle(esdv0);//Opening angle between two daughter tracks | |
303 | ||
304 | fPsiPair = PsiPair(esdv0);//Angle between daughter momentum plane and plane perpendicular to magnetic field | |
305 | ||
306 | fV0Momentum = V0Momentum(esdv0);//Reconstructed momentum of the mother particle | |
307 | ||
3d19c1b0 | 308 | //4 decay types : conversions, K0s, Lambda, Anti-Lambda |
1ee39b3a | 309 | //five particle types: electrons, muons, pions, kaons, protons (muons and kaons not involved) |
3d19c1b0 | 310 | for(Int_t idecay(0), part1(-1), part2(-1); idecay < kNDecays; idecay++){ |
0ed6f095 | 311 | |
312 | fArmenteros[idecay]=Armenteros(esdv0, idecay);//Attribute the Armenteros yes/no decision for every decay type | |
3d19c1b0 | 313 | if(idecay == kLambda){ //protons and pions from Lambda |
314 | part1 = AliPID::kProton; | |
315 | part2 = AliPID::kPion; | |
316 | } else if(idecay == kAntiLambda) { //antiprotons and pions from Anti-Lambda | |
317 | part1 = AliPID::kPion; | |
318 | part2 = AliPID::kProton; | |
319 | } else if(idecay == kK0s) {//pions from K0s | |
320 | part1 = part2 = AliPID::kPion; | |
321 | } else if(idecay == kGamma) {//electrons from conversions | |
322 | part1 = part2 = AliPID::kElectron; | |
323 | } | |
324 | fInvMass[idecay] = InvMass(part1, part2, esdv0);//Calculate invariant mass for all of our four supposed decays | |
24b3cfb9 | 325 | |
326 | // Comment out until bug fix is provided | |
327 | // A.Bercuci 14. July 2010 | |
328 | //fChi2ndf[idecay] = KFChi2ndf(part1, part2,idecay); | |
0ed6f095 | 329 | |
1ee39b3a | 330 | } |
3d19c1b0 | 331 | //Gets all likelihood values from TPC, TOF and ITS PID for the fDetPID[kNDaughters][kNDetectors][AliPID::kSPECIES] array |
332 | GetDetectorPID(); | |
333 | //Bayesian combination of likelihoods from TPC and TOF | |
334 | CombinePID(); | |
0ed6f095 | 335 | //TPC dE/dx values for both tracks |
336 | GetTPCdEdx(); | |
337 | ||
1ee39b3a | 338 | } |
339 | //_________________________________________________ | |
b37d601d | 340 | Float_t AliTRDv0Info::V0Momentum(const AliESDv0 *esdv0) const |
1ee39b3a | 341 | { |
342 | // | |
343 | // Reconstructed momentum of V0 mother particle | |
344 | // | |
345 | ||
346 | Double_t mn[3] = {0,0,0}; | |
347 | Double_t mp[3] = {0,0,0}; | |
348 | ||
349 | ||
350 | esdv0->GetNPxPyPz(mn[0],mn[1],mn[2]);//reconstructed cartesian momentum components of negative daughter; | |
351 | esdv0->GetPPxPyPz(mp[0],mp[1],mp[2]);//reconstructed cartesian momentum components of positive daughter; | |
352 | ||
353 | ||
354 | return TMath::Sqrt((mn[0]+mp[0])*(mn[0]+mp[0]) + (mn[1]+mp[1])*(mn[1]+mp[1])+(mn[2]+mp[2])*(mn[2]+mp[2])); | |
355 | } | |
356 | ||
357 | //_________________________________________________ | |
b37d601d | 358 | Double_t AliTRDv0Info::InvMass(Int_t part1, Int_t part2, const AliESDv0 *esdv0) const |
1ee39b3a | 359 | { |
360 | // | |
361 | // Invariant mass of reconstructed V0 mother | |
362 | // | |
363 | ||
364 | const Double_t kpmass[5] = {AliPID::ParticleMass(AliPID::kElectron),AliPID::ParticleMass(AliPID::kMuon),AliPID::ParticleMass(AliPID::kPion),AliPID::ParticleMass(AliPID::kKaon),AliPID::ParticleMass(AliPID::kProton)}; | |
365 | //Masses of electrons, muons, pions, kaons and protons, as implemented in ROOT | |
366 | ||
367 | ||
368 | Double_t mn[3] = {0,0,0}; | |
369 | Double_t mp[3] = {0,0,0}; | |
370 | ||
371 | esdv0->GetNPxPyPz(mn[0],mn[1],mn[2]);//reconstructed cartesian momentum components of negative daughter; | |
372 | esdv0->GetPPxPyPz(mp[0],mp[1],mp[2]);//reconstructed cartesian momentum components of positive daughter; | |
373 | ||
0ed6f095 | 374 | Double_t mass1 = kpmass[part1];//sets supposed rest masses for both daughters: positive |
375 | Double_t mass2 = kpmass[part2];//negative | |
1ee39b3a | 376 | |
377 | //Calculate daughters' energies : | |
378 | Double_t e1 = TMath::Sqrt(mass1*mass1+ | |
379 | mp[0]*mp[0]+ | |
380 | mp[1]*mp[1]+ | |
381 | mp[2]*mp[2]); | |
382 | Double_t e2 = TMath::Sqrt(mass2*mass2+ | |
383 | mn[0]*mn[0]+ | |
384 | mn[1]*mn[1]+ | |
385 | mn[2]*mn[2]); | |
386 | ||
387 | //Sum of daughter momenta : | |
388 | Double_t momsum = | |
389 | (mn[0]+mp[0])*(mn[0]+mp[0])+ | |
390 | (mn[1]+mp[1])*(mn[1]+mp[1])+ | |
391 | (mn[2]+mp[2])*(mn[2]+mp[2]); | |
392 | ||
393 | //invariant mass : | |
394 | Double_t mInv = TMath::Sqrt((e1+e2)*(e1+e2)-momsum); | |
395 | ||
396 | return mInv; | |
397 | ||
398 | } | |
399 | //_________________________________________________ | |
b37d601d | 400 | Float_t AliTRDv0Info::OpenAngle(const AliESDv0 *esdv0) |
64d57299 | 401 | { |
402 | //Opening angle between two daughter tracks | |
1ee39b3a | 403 | Double_t mn[3] = {0,0,0}; |
404 | Double_t mp[3] = {0,0,0}; | |
405 | ||
406 | ||
407 | esdv0->GetNPxPyPz(mn[0],mn[1],mn[2]);//reconstructed cartesian momentum components of negative daughter; | |
408 | esdv0->GetPPxPyPz(mp[0],mp[1],mp[2]);//reconstructed cartesian momentum components of positive daughter; | |
409 | ||
410 | ||
411 | fOpenAngle = TMath::ACos((mp[0]*mn[0] + mp[1]*mn[1] + mp[2]*mn[2])/(TMath::Sqrt(mp[0]*mp[0] + mp[1]*mp[1] + mp[2]*mp[2])*TMath::Sqrt(mn[0]*mn[0] + mn[1]*mn[1] + mn[2]*mn[2]))); | |
412 | ||
413 | return fOpenAngle; | |
414 | } | |
415 | ||
416 | //_________________________________________________ | |
b37d601d | 417 | Float_t AliTRDv0Info::PsiPair(const AliESDv0 *esdv0) |
64d57299 | 418 | { |
419 | //Angle between daughter momentum plane and plane perpendicular to magnetic field | |
1ee39b3a | 420 | Double_t x, y, z; |
421 | esdv0->GetXYZ(x,y,z);//Reconstructed coordinates of V0; to be replaced by Markus Rammler's method in case of conversions! | |
422 | ||
423 | Double_t mn[3] = {0,0,0}; | |
424 | Double_t mp[3] = {0,0,0}; | |
425 | ||
426 | ||
427 | esdv0->GetNPxPyPz(mn[0],mn[1],mn[2]);//reconstructed cartesian momentum components of negative daughter; | |
428 | esdv0->GetPPxPyPz(mp[0],mp[1],mp[2]);//reconstructed cartesian momentum components of positive daughter; | |
429 | ||
430 | ||
431 | Double_t deltat = 1.; | |
432 | deltat = TMath::ATan(mp[2]/(TMath::Sqrt(mp[0]*mp[0] + mp[1]*mp[1])+1.e-13)) - TMath::ATan(mn[2]/(TMath::Sqrt(mn[0]*mn[0] + mn[1]*mn[1])+1.e-13));//difference of angles of the two daughter tracks with z-axis | |
433 | ||
434 | Double_t radiussum = TMath::Sqrt(x*x + y*y) + 50;//radius to which tracks shall be propagated | |
435 | ||
436 | Double_t momPosProp[3]; | |
437 | Double_t momNegProp[3]; | |
438 | ||
439 | AliExternalTrackParam nt(*fTrackN), pt(*fTrackP); | |
440 | ||
441 | fPsiPair = 4.; | |
442 | ||
443 | if(nt.PropagateTo(radiussum,fMagField) == 0)//propagate tracks to the outside | |
444 | fPsiPair = -5.; | |
445 | if(pt.PropagateTo(radiussum,fMagField) == 0) | |
446 | fPsiPair = -5.; | |
447 | pt.GetPxPyPz(momPosProp);//Get momentum vectors of tracks after propagation | |
448 | nt.GetPxPyPz(momNegProp); | |
449 | ||
450 | Double_t pEle = | |
451 | TMath::Sqrt(momNegProp[0]*momNegProp[0]+momNegProp[1]*momNegProp[1]+momNegProp[2]*momNegProp[2]);//absolute momentum value of negative daughter | |
452 | Double_t pPos = | |
453 | TMath::Sqrt(momPosProp[0]*momPosProp[0]+momPosProp[1]*momPosProp[1]+momPosProp[2]*momPosProp[2]);//absolute momentum value of positive daughter | |
454 | ||
455 | Double_t scalarproduct = | |
456 | momPosProp[0]*momNegProp[0]+momPosProp[1]*momNegProp[1]+momPosProp[2]*momNegProp[2];//scalar product of propagated positive and negative daughters' momenta | |
457 | ||
458 | Double_t chipair = TMath::ACos(scalarproduct/(pEle*pPos));//Angle between propagated daughter tracks | |
459 | ||
460 | fPsiPair = TMath::Abs(TMath::ASin(deltat/chipair)); | |
461 | ||
462 | return fPsiPair; | |
463 | ||
464 | } | |
0ed6f095 | 465 | //_________________________________________________ |
466 | Double_t AliTRDv0Info::KFChi2ndf(Int_t part1, Int_t part2,Int_t decay){ | |
467 | //Calculates Kalman filter Chi2/NDF | |
468 | Int_t mothers[4]={22,310,3122,3122}; | |
469 | ||
470 | const Double_t partMass=TDatabasePDG::Instance()->GetParticle(mothers[decay])->Mass(); | |
471 | const Double_t massWidth[4] = {0.001, 0., 0., 0.}; | |
472 | ||
473 | AliKFParticle *kfMother = CreateMotherParticle(fTrackP, fTrackN, part1, part2); | |
474 | ||
475 | // Lambda | |
476 | if(!kfMother) { | |
477 | return kFALSE; | |
478 | } | |
479 | ||
480 | // production vertex is set in the 'CreateMotherParticle' function | |
481 | kfMother->SetMassConstraint(partMass, massWidth[decay]); | |
482 | ||
483 | Double_t chi2ndf = (kfMother->GetChi2()/kfMother->GetNDF()); | |
484 | ||
485 | if(kfMother)delete kfMother; | |
486 | return chi2ndf; | |
487 | } | |
488 | //________________________________________________________________ | |
64d57299 | 489 | AliKFParticle *AliTRDv0Info::CreateMotherParticle(const AliESDtrack *pdaughter, const AliESDtrack *ndaughter, Int_t pspec, Int_t nspec){ |
0ed6f095 | 490 | // |
491 | // Creates a mother particle | |
492 | // | |
493 | AliKFParticle pkfdaughter(*pdaughter, pspec); | |
494 | AliKFParticle nkfdaughter(*ndaughter, nspec); | |
495 | ||
496 | ||
497 | // Create the mother particle | |
498 | AliKFParticle *m = new AliKFParticle(pkfdaughter, nkfdaughter); | |
499 | ||
500 | AliKFVertex improvedVertex = *fPrimaryVertex; | |
501 | improvedVertex += *m; | |
502 | m->SetProductionVertex(improvedVertex); | |
503 | ||
1ee39b3a | 504 | |
0ed6f095 | 505 | return m; |
506 | } | |
1ee39b3a | 507 | //_________________________________________________ |
64d57299 | 508 | Int_t AliTRDv0Info::HasTrack(const AliTRDtrackInfo * const track) const |
3d19c1b0 | 509 | { |
64d57299 | 510 | //Checks if track is a secondary vertex daughter (due to V0 finder) |
1ee39b3a | 511 | |
d80a6a00 | 512 | if(!track) return 0; |
513 | if(!fTrackP->GetID()) return 0; | |
514 | if(!fTrackN->GetID()) return 0; | |
515 | ||
3d19c1b0 | 516 | Int_t trackID(track->GetTrackId());//index of the track |
b9ddd472 | 517 | return HasTrack(trackID); |
518 | } | |
1ee39b3a | 519 | |
b9ddd472 | 520 | //_________________________________________________ |
64d57299 | 521 | Int_t AliTRDv0Info::HasTrack(Int_t trackID) const |
b9ddd472 | 522 | { |
1ee39b3a | 523 | //comparing index of track with indices of pos./neg. V0 daughter : |
b9ddd472 | 524 | if(fNindex==trackID) return -1; |
525 | else if(fPindex==trackID) return 1; | |
526 | else return 0; | |
1ee39b3a | 527 | } |
3d19c1b0 | 528 | |
1ee39b3a | 529 | //_________________________________________________ |
530 | void AliTRDv0Info::GetDetectorPID() | |
64d57299 | 531 | { |
532 | //PID likelihoods from TPC, TOF, and ITS, for all particle species | |
1ee39b3a | 533 | |
534 | fTrackN->GetTPCpid(fDetPID[kNeg][kTPC]); | |
535 | fTrackP->GetTPCpid(fDetPID[kPos][kTPC]); | |
536 | fTrackN->GetTOFpid(fDetPID[kNeg][kTOF]); | |
537 | fTrackP->GetTOFpid(fDetPID[kPos][kTOF]); | |
538 | fTrackN->GetITSpid(fDetPID[kNeg][kITS]); | |
539 | fTrackP->GetITSpid(fDetPID[kPos][kITS]); | |
540 | ||
8bc8ea55 | 541 | Long_t statusN = fTrackN->GetStatus(); |
542 | Long_t statusP = fTrackP->GetStatus(); | |
543 | ||
544 | if(!(statusN & AliESDtrack::kTPCpid)){ | |
545 | for(Int_t iPart = 0; iPart < AliPID::kSPECIES; iPart++){ | |
546 | fDetPID[kNeg][kTPC][iPart] = 0.2; | |
547 | } | |
548 | } | |
549 | if(!(statusN & AliESDtrack::kTOFpid)){ | |
550 | for(Int_t iPart = 0; iPart < AliPID::kSPECIES; iPart++){ | |
551 | fDetPID[kNeg][kTOF][iPart] = 0.2; | |
552 | } | |
553 | ||
554 | } | |
555 | if(!(statusN & AliESDtrack::kITSpid)){ | |
556 | for(Int_t iPart = 0; iPart < AliPID::kSPECIES; iPart++){ | |
557 | fDetPID[kNeg][kITS][iPart] = 0.2; | |
558 | } | |
559 | } | |
560 | if(!(statusP & AliESDtrack::kTPCpid)){ | |
561 | for(Int_t iPart = 0; iPart < AliPID::kSPECIES; iPart++){ | |
562 | fDetPID[kPos][kTPC][iPart] = 0.2; | |
563 | } | |
564 | } | |
565 | if(!(statusP & AliESDtrack::kTOFpid)){ | |
566 | for(Int_t iPart = 0; iPart < AliPID::kSPECIES; iPart++){ | |
567 | fDetPID[kPos][kTOF][iPart] = 0.2; | |
568 | } | |
569 | ||
570 | } | |
571 | if(!(statusP & AliESDtrack::kITSpid)){ | |
572 | for(Int_t iPart = 0; iPart < AliPID::kSPECIES; iPart++){ | |
573 | fDetPID[kPos][kITS][iPart] = 0.2; | |
574 | } | |
575 | } | |
1ee39b3a | 576 | |
8bc8ea55 | 577 | } |
578 | //____________________________________________________________________________________ | |
579 | void AliTRDv0Info::CombinePID() | |
580 | { | |
64d57299 | 581 | //combined bayesian PID from TPC and TOF |
8bc8ea55 | 582 | Double_t partrat[AliPID::kSPECIES] = {0.208, 0.010, 0.662, 0.019, 0.101}; |
583 | ||
584 | for(Int_t iSign = 0; iSign < kNDaughters; iSign++) | |
585 | { | |
586 | for(Int_t iPart = 0; iPart < AliPID::kSPECIES; iPart++) | |
587 | { | |
588 | fComPID[iSign][iPart] = (partrat[iPart]*fDetPID[iSign][kTPC][iPart]*fDetPID[iSign][kTOF][iPart])/((partrat[0]*fDetPID[iSign][kTPC][0]*fDetPID[iSign][kTOF][0])+(partrat[1]*fDetPID[iSign][kTPC][1]*fDetPID[iSign][kTOF][1])+(partrat[2]*fDetPID[iSign][kTPC][2]*fDetPID[iSign][kTOF][2])+(partrat[3]*fDetPID[iSign][kTPC][3]*fDetPID[iSign][kTOF][3])+(partrat[4]*fDetPID[iSign][kTPC][4]*fDetPID[iSign][kTOF][4])); | |
589 | ||
590 | } | |
591 | } | |
592 | } | |
1ee39b3a | 593 | //_________________________________________________ |
d80a6a00 | 594 | Bool_t AliTRDv0Info::GetTPCdEdx() |
0ed6f095 | 595 | { |
64d57299 | 596 | //gets the TPC dE/dx for both daughter tracks |
d80a6a00 | 597 | if(!fTrackP->GetID()) return 0; |
598 | if(!fTrackN->GetID()) return 0; | |
599 | ||
0ed6f095 | 600 | fTPCdEdx[kNeg] = fTrackN->GetTPCsignal(); |
601 | fTPCdEdx[kPos] = fTrackP->GetTPCsignal(); | |
d80a6a00 | 602 | return 1; |
0ed6f095 | 603 | |
604 | } | |
605 | //_________________________________________________ | |
b37d601d | 606 | Bool_t AliTRDv0Info::TPCdEdxCuts(Int_t part, const AliTRDtrackInfo * const track) |
0ed6f095 | 607 | { |
64d57299 | 608 | //applies cuts on TPC dE/dx according to particle species; cutting lines are drawn shifted to the Bethe-Bloch paremeterization |
d80a6a00 | 609 | if(!fTrackP->GetID()) return 0; |
610 | if(!fTrackN->GetID()) return 0; | |
611 | ||
0ed6f095 | 612 | //Bethe-Bloch lines |
613 | Double_t alephParameters[5]; | |
614 | ||
615 | // data | |
616 | alephParameters[0] = 0.0283086; | |
617 | alephParameters[1] = 2.63394e+01; | |
618 | alephParameters[2] = 5.04114e-11; | |
619 | alephParameters[3] = 2.12543e+00; | |
620 | alephParameters[4] = 4.88663e+00; | |
621 | ||
622 | ||
623 | Double_t deposit = 0; | |
624 | Float_t x = 0; | |
625 | if(HasTrack(track) == 1){ | |
626 | x = fTrackP->P(); | |
627 | deposit = fTPCdEdx[kPos]; | |
628 | } | |
629 | else if(HasTrack(track) == -1){ | |
630 | x = fTrackN->P(); | |
631 | deposit = fTPCdEdx[kNeg]; | |
632 | } | |
633 | else{ | |
634 | printf("No track found"); | |
635 | return 0; | |
636 | } | |
637 | if(x < 0.2)return 0; | |
638 | ||
639 | Float_t upLimits[5]={85,1000,50*AliExternalTrackParam::BetheBlochAleph(x/0.13957, alephParameters[0], alephParameters[1], alephParameters[2], alephParameters[3], alephParameters[4])+6,1000,50*AliExternalTrackParam::BetheBlochAleph(x/0.93827, alephParameters[0], alephParameters[1], alephParameters[2], alephParameters[3], alephParameters[4])+10}; | |
640 | Float_t downLimits[5]={62,40,50*AliExternalTrackParam::BetheBlochAleph(x/0.13957, alephParameters[0], alephParameters[1], alephParameters[2], alephParameters[3], alephParameters[4])-6,40,50*AliExternalTrackParam::BetheBlochAleph(x/0.93827, alephParameters[0], alephParameters[1], alephParameters[2], alephParameters[3], alephParameters[4])-11}; | |
641 | ||
642 | ||
643 | if(x < 0.7){ | |
644 | downLimits[4]=90; | |
645 | } | |
646 | if(x < 1.25){ | |
647 | upLimits[0] = 85; | |
648 | } | |
649 | else{ | |
650 | downLimits[0] = 64; | |
651 | } | |
652 | ||
653 | ||
654 | if(deposit < downLimits[part]) | |
655 | return 0; | |
656 | if(deposit > upLimits[part]) | |
657 | return 0; | |
658 | ||
659 | ||
660 | return 1; | |
661 | ||
662 | } | |
663 | //_________________________________________________ | |
b37d601d | 664 | Float_t AliTRDv0Info::Radius(const AliESDv0 *esdv0) |
64d57299 | 665 | { |
666 | //distance from secondary vertex to primary vertex in x-y plane | |
1ee39b3a | 667 | Double_t x, y, z; |
0ed6f095 | 668 | esdv0->GetXYZ(x,y,z); //Reconstructed coordinates of V0 |
1ee39b3a | 669 | fRadius = TMath::Sqrt(x*x + y*y); |
670 | return fRadius; | |
671 | ||
672 | } | |
673 | ||
674 | //_________________________________________________ | |
b37d601d | 675 | Int_t AliTRDv0Info::Quality(const AliESDv0 *const esdv0) |
1ee39b3a | 676 | { |
677 | // | |
678 | // Checking track and V0 quality status in order to exclude vertices based on poor information | |
679 | // | |
680 | ||
681 | Float_t nClsN; | |
682 | nClsN = fTrackN->GetTPCNcls();//number of found clusters in TPC for negative track | |
683 | Float_t nClsFN; | |
684 | nClsFN = fTrackN->GetTPCNclsF();//number of findable clusters in TPC for negative track | |
685 | Float_t nClsP; | |
686 | nClsP = fTrackP->GetTPCNcls();//number of found clusters in TPC for positive track | |
687 | Float_t nClsFP; | |
688 | nClsFP = fTrackP->GetTPCNclsF();//number of findable clusters in TPC for positive track | |
689 | ||
690 | fQuality = 0; | |
691 | ||
692 | ||
0ed6f095 | 693 | if (!(esdv0->GetOnFlyStatus()))//accept only vertices from online V0 finder |
694 | return -1; | |
695 | ||
1ee39b3a | 696 | Float_t clsRatioN; |
697 | Float_t clsRatioP; | |
698 | ||
0ed6f095 | 699 | if((nClsFN < 80) || (nClsFP < 80)) return -2;//reject all V0s where at least one track has less than 80 TPC clusters |
700 | ||
701 | // Chi2 per TPC cluster | |
702 | Int_t nTPCclustersP = fTrackP->GetTPCclusters(0); | |
703 | Int_t nTPCclustersN = fTrackN->GetTPCclusters(0); | |
704 | Float_t chi2perTPCclusterP = fTrackP->GetTPCchi2()/Float_t(nTPCclustersP); | |
705 | Float_t chi2perTPCclusterN = fTrackN->GetTPCchi2()/Float_t(nTPCclustersN); | |
706 | ||
707 | if((chi2perTPCclusterN > 3.5)||(chi2perTPCclusterP > 3.5)) return -3;//reject all V0s where at least one track has a chi2 above 3.5 | |
1ee39b3a | 708 | |
709 | clsRatioN = nClsN/nClsFN; //ratios of found to findable clusters in TPC | |
710 | clsRatioP = nClsP/nClsFP; | |
0ed6f095 | 711 | |
712 | if((clsRatioN < 0.6)||(clsRatioP < 0.6))//exclude tracks with low ratio of found to findable TPC clusters | |
713 | return -4; | |
714 | ||
1ee39b3a | 715 | if (!((fTrackP->GetStatus() & |
716 | AliESDtrack::kTPCrefit)))//accept only vertices in which both tracks have TPC refit | |
0ed6f095 | 717 | return -5; |
1ee39b3a | 718 | if (!((fTrackN->GetStatus() & |
719 | AliESDtrack::kTPCrefit))) | |
0ed6f095 | 720 | return -6; |
1ee39b3a | 721 | if (fTrackP->GetKinkIndex(0)>0 || |
722 | fTrackN->GetKinkIndex(0)>0 )//exclude tracks with kinks | |
e148d6e6 | 723 | return -7; |
0ed6f095 | 724 | |
725 | if(!(V0SignCheck())) | |
726 | return -8; | |
1ee39b3a | 727 | fQuality = 1; |
728 | return fQuality; | |
729 | } | |
0ed6f095 | 730 | //________________________________________________________________ |
731 | Bool_t AliTRDv0Info::V0SignCheck(){ | |
732 | // | |
733 | // Check if v0 daughters really carry opposite charges | |
734 | // | |
735 | ||
736 | Int_t qP = fTrackP->Charge(); | |
737 | Int_t qN = fTrackN->Charge(); | |
738 | ||
739 | if((qP*qN) != -1) return kFALSE; | |
740 | ||
741 | return kTRUE; | |
742 | } | |
743 | //___________________________________________________________________ | |
b37d601d | 744 | Bool_t AliTRDv0Info::Armenteros(const AliESDv0 *esdv0, Int_t decay){ |
0ed6f095 | 745 | // |
746 | // computes the Armenteros variables for given V0 | |
747 | // | |
748 | Double_t mn[3] = {0,0,0}; | |
749 | Double_t mp[3] = {0,0,0}; | |
750 | Double_t mm[3] = {0,0,0}; | |
751 | ||
752 | if(V0SignCheck()){ | |
753 | esdv0->GetNPxPyPz(mn[0],mn[1],mn[2]); //reconstructed cartesian momentum components of negative daughter | |
754 | esdv0->GetPPxPyPz(mp[0],mp[1],mp[2]); //reconstructed cartesian momentum components of positive daughter | |
755 | } | |
756 | else{ | |
757 | esdv0->GetPPxPyPz(mn[0],mn[1],mn[2]); //reconstructed cartesian momentum components of negative daughter | |
758 | esdv0->GetNPxPyPz(mp[0],mp[1],mp[2]); //reconstructed cartesian momentum components of positive daughter | |
759 | } | |
760 | esdv0->GetPxPyPz(mm[0],mm[1],mm[2]); //reconstructed cartesian momentum components of mother | |
761 | ||
762 | TVector3 vecN(mn[0],mn[1],mn[2]); | |
763 | TVector3 vecP(mp[0],mp[1],mp[2]); | |
764 | TVector3 vecM(mm[0],mm[1],mm[2]); | |
765 | ||
766 | Double_t thetaP = acos((vecP * vecM)/(vecP.Mag() * vecM.Mag())); | |
767 | Double_t thetaN = acos((vecN * vecM)/(vecN.Mag() * vecM.Mag())); | |
768 | ||
769 | Double_t alfa = ((vecP.Mag())*cos(thetaP)-(vecN.Mag())*cos(thetaN))/ | |
770 | ((vecP.Mag())*cos(thetaP)+(vecN.Mag())*cos(thetaN)) ; | |
771 | Double_t qt = vecP.Mag()*sin(thetaP); | |
772 | ||
773 | Float_t ap[2]; | |
774 | ap[0] = alfa; | |
775 | ap[1] = qt; | |
776 | ||
64d57299 | 777 | Double_t lCutAP[2];//Lambda/Anti-Lambda cuts |
0ed6f095 | 778 | if(decay == 0){ |
779 | // armenteros cuts | |
780 | const Double_t cutAlpha[2] = {0.35, 0.45}; // [0.35, 0.45] | |
781 | const Double_t cutQT = 0.015; | |
782 | if(TMath::Abs(ap[0]) > cutAlpha[0] && TMath::Abs(ap[0]) < cutAlpha[1]) return kFALSE; | |
783 | ||
784 | if(ap[1] > cutQT) return kFALSE; | |
785 | } | |
786 | ||
787 | else if(decay == 1){ | |
788 | const Double_t cutQT = 0.1075; | |
789 | const Double_t cutAP = 0.22 * TMath::Sqrt( TMath::Abs( (1-ap[0]*ap[0]/(0.92*0.92)) ) ); | |
790 | if(ap[1] < cutQT) return kFALSE; | |
791 | if(ap[1] > cutAP) return kFALSE; | |
792 | } | |
793 | else if(decay == 2){ | |
794 | const Double_t cutQT = 0.03; | |
795 | const Double_t cutAlpha = 0.7; // VERY strong - should supress the overlap with K0 | |
64d57299 | 796 | lCutAP[0] = 1.0 - (ap[0]-0.7 * ap[0]-0.7)*1.1 - 0.87; |
0ed6f095 | 797 | if(TMath::Abs(ap[0]) > cutAlpha) return kFALSE; |
798 | if(ap[1] < cutQT) return kFALSE; | |
64d57299 | 799 | if(ap[1] > lCutAP[0]) return kFALSE; |
0ed6f095 | 800 | |
801 | } | |
802 | else if(decay == 3){ | |
803 | const Double_t cutQT = 0.03; | |
804 | const Double_t cutAlpha = 0.7; // VERY strong - should supress the overlap with K0 | |
64d57299 | 805 | lCutAP[1] = 1.0 - (ap[0]+0.7 * ap[0]+0.7)*1.1 - 0.87; |
0ed6f095 | 806 | if(TMath::Abs(ap[0]) > cutAlpha) return kFALSE; |
807 | if(ap[1] < cutQT) return kFALSE; | |
64d57299 | 808 | if(ap[1] > lCutAP[1]) return kFALSE; |
0ed6f095 | 809 | } |
810 | return kTRUE; | |
811 | } | |
1ee39b3a | 812 | //_________________________________________________ |
3d19c1b0 | 813 | Int_t AliTRDv0Info::GetPID(Int_t ipart, AliTRDtrackInfo *track) |
814 | { | |
0ed6f095 | 815 | // Decides if track is accepted for one of the reference data samples |
816 | Int_t cutCode = -99; | |
3d19c1b0 | 817 | if(!(track)) { |
818 | AliError("No track info"); | |
b9ddd472 | 819 | return -1; |
3d19c1b0 | 820 | } |
821 | if(!HasTrack(track)){ | |
822 | AliDebug(2, "Track not attached to v0."); | |
0ed6f095 | 823 | return -2; |
3d19c1b0 | 824 | } |
3d19c1b0 | 825 | |
826 | //translate ipart to decay (Anti-Lambda will be treated separately) | |
1ee39b3a | 827 | Int_t iDecay = -1; |
3d19c1b0 | 828 | switch(ipart){ |
829 | case AliPID::kElectron: iDecay = kGamma; break; | |
830 | case AliPID::kPion: iDecay = kK0s; break; | |
831 | case AliPID::kProton: iDecay = kLambda; break; | |
832 | default: | |
b2f4ab8d | 833 | AliDebug(1, Form("Hypothesis \"ipart=%d\" not handled", ipart)); |
0ed6f095 | 834 | return -3; |
3d19c1b0 | 835 | } |
1ee39b3a | 836 | |
3d19c1b0 | 837 | //... it fulfills our quality criteria |
0ed6f095 | 838 | if(!(fQuality == 1)) return -4; |
3d19c1b0 | 839 | //... distance of closest approach between daughters is reasonably small |
0ed6f095 | 840 | if((fDCA > fUpDCA[iDecay])) return -5; |
3d19c1b0 | 841 | //... pointing angle between momentum of mother particle and vector from prim. to sec. vertex is small |
0ed6f095 | 842 | if((fPointingAngle > fUpPointingAngle[iDecay])) return -6; |
3d19c1b0 | 843 | //... x-y plane distance of decay point to prim. vertex is bigger than a certain minimum value (for conversions) |
0ed6f095 | 844 | if((fRadius < fDownRadius[iDecay])) return -7; |
3d19c1b0 | 845 | //...or smaller than a maximum value (for K0s) |
0ed6f095 | 846 | if((fRadius > fUpRadius[iDecay])) return -8; |
3d19c1b0 | 847 | //... opening angle is close enough to zero (for conversions) |
0ed6f095 | 848 | if((fOpenAngle > fUpOpenAngle[iDecay])) return -9; |
3d19c1b0 | 849 | //... Psi-pair angle is close enough to zero(for conversions) |
0ed6f095 | 850 | if((TMath::Abs(fPsiPair) > fUpPsiPair[iDecay])) return -10; |
851 | ||
3d19c1b0 | 852 | |
853 | ||
854 | //Mother momentum slots above/below 2.5 GeV | |
855 | Int_t iPSlot(fV0Momentum > 2.5); | |
b9ddd472 | 856 | Int_t trackID(track->GetTrackId()); |
3d19c1b0 | 857 | |
858 | //specific cut criteria : | |
859 | if(ipart == AliPID::kProton) { | |
0ed6f095 | 860 | if((fInvMass[kK0s] < fUpInvMass[kK0s][iPSlot]) && (fInvMass[kK0s] > fDownInvMass[kK0s])) return -11;//explicit exclusion of K0s decays |
861 | ||
862 | if(fOpenAngle < (0.3 - 0.2*fV0Momentum))return -9; | |
863 | ||
864 | ||
865 | ||
3d19c1b0 | 866 | //for proton sample: separate treatment of Lamba and Anti-Lambda decays: |
867 | //for Anti-Lambda: | |
868 | //Combined PID likelihoods high enough for pi+ and anti-proton ; invariant mass calculated postulating these two particle species... | |
0ed6f095 | 869 | //if((fComPID[kNeg][AliPID::kProton] > fDownComPIDneg[AliPID::kProton]) && (fComPID[kPos][AliPID::kPion] > fDownComPIDposPart[AliPID::kPion])) { |
870 | //if((fDetPID[kNeg][kTPC][AliPID::kProton] > fDownTPCPIDneg[AliPID::kProton]) && (fDetPID[kPos][kTPC][AliPID::kPion] > fDownTPCPIDpos[AliPID::kPion])){ | |
871 | if((TPCdEdxCuts(ipart, track))){//momentary solution: direct cut on TPC dE/dx | |
872 | if(fNindex == trackID) {//we're only interested in the anti-proton | |
873 | if(fArmenteros[kAntiLambda]){//Armenteros condition has to be fulfilled | |
874 | if(fChi2ndf[kAntiLambda] < fUpChi2ndf[kAntiLambda]){//Kalman filter Chi2/NDF not allowed to be too large | |
875 | if((fInvMass[kAntiLambda] < fUpInvMass[kAntiLambda][iPSlot]) && (fInvMass[kAntiLambda] > fDownInvMass[kAntiLambda])){ | |
876 | return 1; | |
877 | } else cutCode = -15; | |
878 | } | |
879 | else cutCode =-14; | |
880 | } | |
881 | else cutCode = -13; | |
3d19c1b0 | 882 | } |
883 | } | |
0ed6f095 | 884 | else cutCode = -12; |
3d19c1b0 | 885 | //for Lambda: |
886 | //TPC PID likelihoods high enough for pi- and proton ; invariant mass calculated accordingly | |
0ed6f095 | 887 | //if((fComPID[kNeg][AliPID::kPion] > fDownComPIDnegPart[AliPID::kPion]) && (fComPID[kPos][AliPID::kProton] > fDownComPIDpos[AliPID::kProton])) { |
888 | //if((fDetPID[kNeg][kTPC][AliPID::kPion] > fDownTPCPIDneg[AliPID::kPion]) && (fDetPID[kPos][kTPC][AliPID::kProton] > fDownTPCPIDpos[AliPID::kProton])){ | |
889 | if((TPCdEdxCuts(ipart, track))){//momentary solution: direct TPC dE/dx cuts | |
3d19c1b0 | 890 | if(fPindex == trackID) { |
0ed6f095 | 891 | if(fArmenteros[kLambda]){ |
892 | if(fChi2ndf[kLambda] < fUpChi2ndf[kLambda]){ | |
893 | if((fInvMass[kLambda] < fUpInvMass[kLambda][iPSlot]) && (fInvMass[kLambda] > fDownInvMass[kLambda])){ | |
894 | return 1; | |
895 | } else cutCode = -15; | |
896 | } | |
897 | else cutCode = -14; | |
898 | } | |
899 | else cutCode = -13; | |
3d19c1b0 | 900 | } |
901 | } | |
0ed6f095 | 902 | else cutCode = -12; |
903 | return cutCode; | |
3d19c1b0 | 904 | } |
0ed6f095 | 905 | |
3d19c1b0 | 906 | //for K0s decays: equal TPC PID likelihood criteria for both daughters ; invariant mass calculated postulating two pions |
907 | if(ipart == AliPID::kPion) { | |
0ed6f095 | 908 | |
909 | if(fOpenAngle < (1.0/(fV0Momentum + 0.3) - 0.1)) | |
910 | return -9; | |
911 | ||
3d19c1b0 | 912 | //explicit exclusion of Lambda decays |
0ed6f095 | 913 | if((fInvMass[kLambda] < fUpInvMass[kLambda][iPSlot]) && (fInvMass[kLambda] > fDownInvMass[kLambda])) return -11; |
3d19c1b0 | 914 | //explicit exclusion of Anti-Lambda decays |
0ed6f095 | 915 | if((fInvMass[kAntiLambda] < fUpInvMass[kAntiLambda][iPSlot]) && (fInvMass[kAntiLambda] > fDownInvMass[kAntiLambda])) return -11; |
916 | ||
917 | //if((fDetPID[kNeg][kTPC][ipart] < fDownTPCPIDneg[ipart]) || (fDetPID[kPos][kTPC][ipart] < fDownTPCPIDpos[ipart])) return -12; | |
918 | if(!(TPCdEdxCuts(ipart, track))){//momentary solution: direct TPC dE/dx cuts | |
919 | return -12; | |
920 | } | |
3d19c1b0 | 921 | } |
0ed6f095 | 922 | |
923 | ||
3d19c1b0 | 924 | //for photon conversions: equal combined PID likelihood criteria for both daughters ; invariant mass calculated postulating two electrons |
925 | //No Lambda/K0s exclusion is provided, since these contributions hardly ever interfere with gamma invariant mass! | |
0ed6f095 | 926 | //Float_t momentum(track->GetESDinfo()->GetOuterParam()->P()); |
3d19c1b0 | 927 | if(ipart == AliPID::kElectron) { |
0ed6f095 | 928 | //if(momentum > 1.75) {//since combined PID performs a little worse in simulations than TPC standalone for higher momenta, ONLY TPC PID is used here |
929 | //if((fDetPID[kNeg][kTPC][ipart] < fDownTPCPIDneg[ipart]) || (fDetPID[kPos][kTPC][ipart] < fDownTPCPIDpos[ipart])) return -12; | |
930 | //} else {//for low momenta, combined PID from TOF and TPC is used to get rid of proton contamination | |
931 | //if((fComPID[kNeg][ipart] > fDownComPIDneg[ipart]) && (fComPID[kPos][ipart] > fDownComPIDpos[ipart])) return 1; | |
932 | //} | |
933 | if(!(TPCdEdxCuts(ipart, track))){//momentary solution for direct TPC dE/dx cut | |
934 | return -12; | |
1ee39b3a | 935 | } |
0ed6f095 | 936 | |
937 | } | |
938 | ||
939 | ||
940 | //Armenteros-Polanski cut | |
941 | if(!(fArmenteros[iDecay])) return -13; | |
942 | ||
943 | //Kalman filter Chi2/NDF cut | |
944 | if(fChi2ndf[iDecay] > fUpChi2ndf[iDecay]) return -14; | |
945 | ||
946 | //Invariant mass cut for K0s and photons, assuming two pions/two electrons as daughters: | |
947 | ||
948 | if((fInvMass[iDecay] > fUpInvMass[iDecay][iPSlot]) || (fInvMass[iDecay] < fDownInvMass[iDecay])) { | |
949 | return -15; | |
950 | ||
3d19c1b0 | 951 | } |
0ed6f095 | 952 | |
953 | return 1; | |
1ee39b3a | 954 | } |
3d19c1b0 | 955 | |
956 | ||
1ee39b3a | 957 | //_________________________________________________ |
b9ddd472 | 958 | void AliTRDv0Info::Print(Option_t *opt) const |
1ee39b3a | 959 | { |
64d57299 | 960 | //prints text for debugging etc. |
b9ddd472 | 961 | printf("V0 P[%d] N[%d]\n", fPindex, fNindex); |
962 | printf(" DCA[%5.3f] Radius[%5.3f]\n", fDCA, fRadius); | |
963 | printf(" Angles : Pointing[%5.3f] Open[%5.3f] Psi[%5.3f]\n", fPointingAngle, fOpenAngle, fPsiPair); | |
964 | if(strcmp(opt, "a")!=0) return; | |
965 | printf(" Reconstructed PID\n" | |
966 | " sgn spec ITS TPC TOF COM\n"); | |
967 | for(Int_t idt=0; idt<kNDaughters; idt++){ | |
968 | printf(" %c", idt?'-':'+'); | |
969 | for(Int_t is(0); is<AliPID::kSPECIES; is++){ | |
970 | printf("%s%s%s", is==0?" ":" ", AliPID::ParticleShortName(is), (is==1||is==2)?" ":" "); | |
971 | for(Int_t id(0); id<kNDetectors; id++){ | |
972 | printf("%5.1f ", 1.e2*fDetPID[idt][id][is]); | |
973 | } | |
974 | printf("%5.1f\n", 1.e2*fComPID[idt][is]); | |
975 | } | |
976 | } | |
977 | } | |
1ee39b3a | 978 | |
b9ddd472 | 979 | //_________________________________________________ |
980 | void AliTRDv0Info::SetV0tracks(AliESDtrack *p, AliESDtrack *n) | |
981 | { | |
64d57299 | 982 | //sets the two daughter trex and their indices |
b9ddd472 | 983 | fTrackP = p; fPindex = p->GetID(); |
984 | fTrackN = n; fNindex = n->GetID(); | |
1ee39b3a | 985 | } |
64d57299 | 986 | //_________________________________________________ |
b9ddd472 | 987 | |
64d57299 | 988 | AliESDtrack *AliTRDv0Info::GetV0Daughter(Int_t sign) |
989 | { | |
990 | //Gets positive of negative daughter of decay | |
991 | if(sign>0) | |
992 | return fTrackP; | |
993 | else if(sign < 0) | |
994 | return fTrackN; | |
b9ddd472 | 995 | |
64d57299 | 996 | return 0; |
997 | } |