]>
Commit | Line | Data |
---|---|---|
d92ed900 | 1 | /////////////////////////////////////////////////////////////////////////// |
2 | // // | |
3 | // AliFemtoESDTrackCut: A basic track cut that used information from // | |
4 | // ALICE ESD to accept or reject the track. // | |
5 | // Enables the selection on charge, transverse momentum, rapidity, // | |
6 | // pid probabilities, number of ITS and TPC clusters // | |
7 | // Author: Marek Chojnacki (WUT), mchojnacki@knf.pw.edu.pl // | |
8 | // // | |
9 | /////////////////////////////////////////////////////////////////////////// | |
67427ff7 | 10 | |
d92ed900 | 11 | #ifndef ALIFEMTOESDTRACKCUT_H |
12 | #define ALIFEMTOESDTRACKCUT_H | |
67427ff7 | 13 | |
14 | //#ifndef StMaker_H | |
15 | //#include "StMaker.h" | |
16 | //#endif | |
17 | ||
65423af9 | 18 | #include "AliFemtoTrackCut.h" |
67427ff7 | 19 | |
20 | class AliFemtoESDTrackCut : public AliFemtoTrackCut | |
21 | { | |
271f1c86 | 22 | public: |
23 | ||
24 | enum PIDMethodType {knSigma=0, kContour=1}; | |
25 | typedef enum PIDMethodType ReadPIDMethodType; | |
67427ff7 | 26 | |
d92ed900 | 27 | AliFemtoESDTrackCut(); |
cc5faabc | 28 | virtual ~AliFemtoESDTrackCut(); |
67427ff7 | 29 | |
d92ed900 | 30 | virtual bool Pass(const AliFemtoTrack* aTrack); |
67427ff7 | 31 | |
d92ed900 | 32 | virtual AliFemtoString Report(); |
3a74204a | 33 | virtual TList *ListSettings(); |
0b3bd1ac | 34 | virtual AliFemtoParticleType Type(){return hbtTrack;} |
67427ff7 | 35 | |
d92ed900 | 36 | void SetPt(const float& lo, const float& hi); |
37 | void SetRapidity(const float& lo, const float& hi); | |
c5c1b5e8 | 38 | void SetEta(const float& lo, const float& hi); |
cc5faabc | 39 | void SetCharge(const int& ch); |
d92ed900 | 40 | void SetPidProbElectron(const float& lo, const float& hi); |
41 | void SetPidProbPion(const float& lo, const float& hi); | |
42 | void SetPidProbKaon(const float& lo, const float& hi); | |
43 | void SetPidProbProton(const float& lo, const float& hi); | |
44 | void SetPidProbMuon(const float& lo, const float& hi); | |
45 | void SetLabel(const bool& flag); | |
46 | void SetStatus(const long& w); | |
47 | void SetminTPCclsF(const short& s); | |
0b3bd1ac | 48 | void SetminTPCncls(const short& s); |
d92ed900 | 49 | void SetminITScls(const int& s); |
cc5faabc | 50 | void SetRemoveKinks(const bool& flag); |
3ecdef05 | 51 | void SetRemoveITSFake(const bool& flag); |
0b3bd1ac | 52 | void SetMaxITSChiNdof(const float& maxchi); |
53 | void SetMaxTPCChiNdof(const float& maxchi); | |
54 | void SetMaxSigmaToVertex(const float& maxsig); | |
ab74d4d3 | 55 | void SetMaxImpactXY(const float& maximpxy); |
56 | void SetMaxImpactZ(const float& maximpz); | |
b2722e9a | 57 | void SetMaxImpactXYPtDep(const float& maxoff, const float& maxnrm, const float& maxpow); |
cc5faabc | 58 | void SetMostProbablePion(); |
59 | void SetMostProbableKaon(); | |
60 | void SetMostProbableProton(); | |
61 | void SetNoMostProbable(); | |
271f1c86 | 62 | void SetPIDMethod(ReadPIDMethodType newMethod); |
cc5faabc | 63 | |
e3f643f0 | 64 | void SetMomRangeTOFpidIs(const float& minp, const float& maxp); |
65 | void SetMomRangeTPCpidIs(const float& minp, const float& maxp); | |
66 | void SetMomRangeITSpidIs(const float& minp, const float& maxp); | |
67 | ||
d92ed900 | 68 | private: // here are the quantities I want to cut on... |
67427ff7 | 69 | |
d92ed900 | 70 | int fCharge; // particle charge |
71 | float fPt[2]; // bounds for transverse momentum | |
72 | float fRapidity[2]; // bounds for rapidity | |
c5c1b5e8 | 73 | float fEta[2]; // bounds for pseudorapidity |
d92ed900 | 74 | float fPidProbElectron[2]; // bounds for electron probability |
75 | float fPidProbPion[2]; // bounds for pion probability | |
76 | float fPidProbKaon[2]; // bounds for kaon probability | |
77 | float fPidProbProton[2]; // bounds for proton probability | |
78 | float fPidProbMuon[2]; // bounds for muon probability | |
79 | bool fLabel; // if true label<0 will not pass throught | |
80 | long fStatus; // staus flag | |
271f1c86 | 81 | ReadPIDMethodType fPIDMethod; // which PID mehod to use. 0 - nsgima, 1 - contour |
db2d5b4b | 82 | |
d92ed900 | 83 | short fminTPCclsF; // min number of findable clusters in the TPC |
0b3bd1ac | 84 | short fminTPCncls; // min number of clusters in the TPC |
d92ed900 | 85 | int fminITScls; // min number of clusters assigned in the ITS |
0b3bd1ac | 86 | float fMaxITSchiNdof; // maximum allowed chi2/ndof for ITS clusters |
87 | float fMaxTPCchiNdof; // maximum allowed chi2/ndof for TPC clusters | |
88 | float fMaxSigmaToVertex; // maximum allowed sigma to primary vertex | |
d92ed900 | 89 | long fNTracksPassed; // passed tracks count |
90 | long fNTracksFailed; // failed tracks count | |
cc5faabc | 91 | bool fRemoveKinks; // if true particles with any kink label will not pass |
3ecdef05 | 92 | bool fRemoveITSFake; // if true particles with ITS fake flag will not pass |
cc5faabc | 93 | int fMostProbable; // this particle type is required to be most probable |
94 | ||
ab74d4d3 | 95 | float fMaxImpactXY; // Max XY impact parameter |
96 | float fMaxImpactZ; // Max Z impact parameter | |
97 | ||
b2722e9a | 98 | float fMaxImpactXYPtOff; // Max XY DCA Pt dependent offset |
99 | float fMaxImpactXYPtNrm; // Max XY DCA Pt dependent normalization | |
100 | float fMaxImpactXYPtPow; // Max XY DCA Pt dependent power | |
101 | ||
e3f643f0 | 102 | float fMinPforTOFpid; // momentum from which TOF PID is requested |
103 | float fMaxPforTOFpid; // momentum till which TOF PID is requested | |
104 | float fMinPforTPCpid; // momentum from which TPC PID is requested | |
105 | float fMaxPforTPCpid; // momentum till which TPC PID is requested | |
106 | float fMinPforITSpid; // momentum from which ITS PID is requested | |
107 | float fMaxPforITSpid; // momentum till which ITS PID is requested | |
108 | ||
cc5faabc | 109 | float PidFractionElectron(float mom) const; |
110 | float PidFractionPion(float mom) const; | |
111 | float PidFractionKaon(float mom) const; | |
112 | float PidFractionProton(float mom) const; | |
67427ff7 | 113 | |
c5c1b5e8 | 114 | bool IsPionTPCdEdx(float mom, float dEdx); |
115 | bool IsKaonTPCdEdx(float mom, float dEdx); | |
3c483111 | 116 | bool IsProtonTPCdEdx(float mom, float dEdx); |
117 | ||
118 | bool IsPionTOFTime(float mom, float ttof); | |
119 | bool IsKaonTOFTime(float mom, float ttof); | |
120 | bool IsProtonTOFTime(float mom, float ttof); | |
c5c1b5e8 | 121 | |
271f1c86 | 122 | bool IsKaonTPCdEdxNSigma(float mom, float nsigma); |
123 | bool IsKaonTOFNSigma(float mom, float nsigma); | |
124 | bool IsKaonNSigma(float mom, float nsigmaTPC, float nsigmaTOF); | |
125 | bool IsPionNSigma(float mom, float nsigmaTPC, float nsigmaTOF); | |
126 | bool IsProtonNSigma(float mom, float nsigmaTPC, float nsigmaTOF); | |
127 | ||
67427ff7 | 128 | #ifdef __ROOT__ |
129 | ClassDef(AliFemtoESDTrackCut, 1) | |
130 | #endif | |
d92ed900 | 131 | }; |
67427ff7 | 132 | |
133 | ||
134 | inline void AliFemtoESDTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;} | |
135 | inline void AliFemtoESDTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;} | |
c5c1b5e8 | 136 | inline void AliFemtoESDTrackCut::SetEta(const float& lo,const float& hi){fEta[0]=lo; fEta[1]=hi;} |
67427ff7 | 137 | inline void AliFemtoESDTrackCut::SetCharge(const int& ch){fCharge = ch;} |
138 | inline void AliFemtoESDTrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;} | |
139 | inline void AliFemtoESDTrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;} | |
140 | inline void AliFemtoESDTrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;} | |
141 | inline void AliFemtoESDTrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;} | |
142 | inline void AliFemtoESDTrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;} | |
143 | inline void AliFemtoESDTrackCut::SetLabel(const bool& flag){fLabel=flag;} | |
144 | inline void AliFemtoESDTrackCut::SetStatus(const long& status){fStatus=status;} | |
145 | inline void AliFemtoESDTrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;} | |
0b3bd1ac | 146 | inline void AliFemtoESDTrackCut::SetminTPCncls(const short& s){fminTPCncls=s;} |
67427ff7 | 147 | inline void AliFemtoESDTrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;} |
cc5faabc | 148 | inline void AliFemtoESDTrackCut::SetMostProbablePion() { fMostProbable = 2; } |
149 | inline void AliFemtoESDTrackCut::SetMostProbableKaon() { fMostProbable = 3; } | |
150 | inline void AliFemtoESDTrackCut::SetMostProbableProton() { fMostProbable = 4; } | |
151 | inline void AliFemtoESDTrackCut::SetNoMostProbable() { fMostProbable = 0; } | |
0b3bd1ac | 152 | inline void AliFemtoESDTrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; } |
153 | inline void AliFemtoESDTrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; } | |
154 | inline void AliFemtoESDTrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; } | |
ab74d4d3 | 155 | inline void AliFemtoESDTrackCut::SetMaxImpactXY(const float& maximpxy) { fMaxImpactXY = maximpxy; } |
b2722e9a | 156 | inline void AliFemtoESDTrackCut::SetMaxImpactXYPtDep(const float& maxoff, const float& maxnrm, const float& maxpow) { fMaxImpactXYPtOff = maxoff; fMaxImpactXYPtNrm = maxnrm; fMaxImpactXYPtPow = maxpow; } |
ab74d4d3 | 157 | inline void AliFemtoESDTrackCut::SetMaxImpactZ(const float& maximpz) { fMaxImpactZ = maximpz; } |
67427ff7 | 158 | |
159 | #endif | |
160 |