1 ///////////////////////////////////////////////////////////////////////////
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 //
9 ///////////////////////////////////////////////////////////////////////////
11 #ifndef ALIFEMTOESDTRACKCUT_H
12 #define ALIFEMTOESDTRACKCUT_H
15 //#include "StMaker.h"
18 #include "AliFemtoTrackCut.h"
20 class AliFemtoESDTrackCut : public AliFemtoTrackCut
24 AliFemtoESDTrackCut();
25 virtual ~AliFemtoESDTrackCut();
27 virtual bool Pass(const AliFemtoTrack* aTrack);
29 virtual AliFemtoString Report();
30 virtual TList *ListSettings();
31 virtual AliFemtoParticleType Type(){return hbtTrack;}
33 void SetPt(const float& lo, const float& hi);
34 void SetRapidity(const float& lo, const float& hi);
35 void SetEta(const float& lo, const float& hi);
36 void SetCharge(const int& ch);
37 void SetPidProbElectron(const float& lo, const float& hi);
38 void SetPidProbPion(const float& lo, const float& hi);
39 void SetPidProbKaon(const float& lo, const float& hi);
40 void SetPidProbProton(const float& lo, const float& hi);
41 void SetPidProbMuon(const float& lo, const float& hi);
42 void SetLabel(const bool& flag);
43 void SetStatus(const long& w);
44 void SetminTPCclsF(const short& s);
45 void SetminTPCncls(const short& s);
46 void SetminITScls(const int& s);
47 void SetRemoveKinks(const bool& flag);
48 void SetMaxITSChiNdof(const float& maxchi);
49 void SetMaxTPCChiNdof(const float& maxchi);
50 void SetMaxSigmaToVertex(const float& maxsig);
51 void SetMaxImpactXY(const float& maximpxy);
52 void SetMaxImpactZ(const float& maximpz);
53 void SetMostProbablePion();
54 void SetMostProbableKaon();
55 void SetMostProbableProton();
56 void SetNoMostProbable();
58 void SetMomRangeTOFpidIs(const float& minp, const float& maxp);
59 void SetMomRangeTPCpidIs(const float& minp, const float& maxp);
60 void SetMomRangeITSpidIs(const float& minp, const float& maxp);
62 private: // here are the quantities I want to cut on...
64 int fCharge; // particle charge
65 float fPt[2]; // bounds for transverse momentum
66 float fRapidity[2]; // bounds for rapidity
67 float fEta[2]; // bounds for pseudorapidity
68 float fPidProbElectron[2]; // bounds for electron probability
69 float fPidProbPion[2]; // bounds for pion probability
70 float fPidProbKaon[2]; // bounds for kaon probability
71 float fPidProbProton[2]; // bounds for proton probability
72 float fPidProbMuon[2]; // bounds for muon probability
73 bool fLabel; // if true label<0 will not pass throught
74 long fStatus; // staus flag
76 short fminTPCclsF; // min number of findable clusters in the TPC
77 short fminTPCncls; // min number of clusters in the TPC
78 int fminITScls; // min number of clusters assigned in the ITS
79 float fMaxITSchiNdof; // maximum allowed chi2/ndof for ITS clusters
80 float fMaxTPCchiNdof; // maximum allowed chi2/ndof for TPC clusters
81 float fMaxSigmaToVertex; // maximum allowed sigma to primary vertex
82 long fNTracksPassed; // passed tracks count
83 long fNTracksFailed; // failed tracks count
84 bool fRemoveKinks; // if true particles with any kink label will not pass
85 int fMostProbable; // this particle type is required to be most probable
87 float fMaxImpactXY; // Max XY impact parameter
88 float fMaxImpactZ; // Max Z impact parameter
90 float fMinPforTOFpid; // momentum from which TOF PID is requested
91 float fMaxPforTOFpid; // momentum till which TOF PID is requested
92 float fMinPforTPCpid; // momentum from which TPC PID is requested
93 float fMaxPforTPCpid; // momentum till which TPC PID is requested
94 float fMinPforITSpid; // momentum from which ITS PID is requested
95 float fMaxPforITSpid; // momentum till which ITS PID is requested
97 float PidFractionElectron(float mom) const;
98 float PidFractionPion(float mom) const;
99 float PidFractionKaon(float mom) const;
100 float PidFractionProton(float mom) const;
102 bool IsPionTPCdEdx(float mom, float dEdx);
103 bool IsKaonTPCdEdx(float mom, float dEdx);
106 ClassDef(AliFemtoESDTrackCut, 1)
111 inline void AliFemtoESDTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
112 inline void AliFemtoESDTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
113 inline void AliFemtoESDTrackCut::SetEta(const float& lo,const float& hi){fEta[0]=lo; fEta[1]=hi;}
114 inline void AliFemtoESDTrackCut::SetCharge(const int& ch){fCharge = ch;}
115 inline void AliFemtoESDTrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
116 inline void AliFemtoESDTrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
117 inline void AliFemtoESDTrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
118 inline void AliFemtoESDTrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
119 inline void AliFemtoESDTrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
120 inline void AliFemtoESDTrackCut::SetLabel(const bool& flag){fLabel=flag;}
121 inline void AliFemtoESDTrackCut::SetStatus(const long& status){fStatus=status;}
122 inline void AliFemtoESDTrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;}
123 inline void AliFemtoESDTrackCut::SetminTPCncls(const short& s){fminTPCncls=s;}
124 inline void AliFemtoESDTrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;}
125 inline void AliFemtoESDTrackCut::SetMostProbablePion() { fMostProbable = 2; }
126 inline void AliFemtoESDTrackCut::SetMostProbableKaon() { fMostProbable = 3; }
127 inline void AliFemtoESDTrackCut::SetMostProbableProton() { fMostProbable = 4; }
128 inline void AliFemtoESDTrackCut::SetNoMostProbable() { fMostProbable = 0; }
129 inline void AliFemtoESDTrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; }
130 inline void AliFemtoESDTrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; }
131 inline void AliFemtoESDTrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }
132 inline void AliFemtoESDTrackCut::SetMaxImpactXY(const float& maximpxy) { fMaxImpactXY = maximpxy; }
133 inline void AliFemtoESDTrackCut::SetMaxImpactZ(const float& maximpz) { fMaxImpactZ = maximpz; }