1 ///////////////////////////////////////////////////////////////////////////
3 // AliFemtoAODTrackCut: A basic track cut that used information from //
4 // ALICE AOD 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: Adam Kisiel (WUT, OSU), Adam.Kisiel@cern.ch //
9 ///////////////////////////////////////////////////////////////////////////
11 #ifndef ALIFEMTOAODTRACKCUT_H
12 #define ALIFEMTOAODTRACKCUT_H
15 //#include "StMaker.h"
18 #include "AliFemtoTrackCut.h"
20 class AliFemtoAODTrackCut : public AliFemtoTrackCut
24 AliFemtoAODTrackCut();
25 virtual ~AliFemtoAODTrackCut();
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 SetCharge(const int& ch);
36 void SetPidProbElectron(const float& lo, const float& hi);
37 void SetPidProbPion(const float& lo, const float& hi);
38 void SetPidProbKaon(const float& lo, const float& hi);
39 void SetPidProbProton(const float& lo, const float& hi);
40 void SetPidProbMuon(const float& lo, const float& hi);
41 void SetLabel(const bool& flag);
42 void SetMaxChiNdof(const float& maxchi);
43 void SetMaxSigmaToVertex(const float& maxsig);
44 void SetMostProbablePion();
45 void SetMostProbableKaon();
46 void SetMostProbableProton();
47 void SetNoMostProbable();
49 private: // here are the quantities I want to cut on...
51 int fCharge; // particle charge
52 float fPt[2]; // bounds for transverse momentum
53 float fRapidity[2]; // bounds for rapidity
54 float fPidProbElectron[2]; // bounds for electron probability
55 float fPidProbPion[2]; // bounds for pion probability
56 float fPidProbKaon[2]; // bounds for kaon probability
57 float fPidProbProton[2]; // bounds for proton probability
58 float fPidProbMuon[2]; // bounds for muon probability
59 bool fLabel; // if true label<0 will not pass throught
60 float fMaxchiNdof; // maximum allowed chi2/ndof for TPC clusters
61 float fMaxSigmaToVertex; // maximum allowed sigma to primary vertex
62 long fNTracksPassed; // passed tracks count
63 long fNTracksFailed; // failed tracks count
64 int fMostProbable; // this particle type is required to be most probable
66 float PidFractionElectron(float mom) const;
67 float PidFractionPion(float mom) const;
68 float PidFractionKaon(float mom) const;
69 float PidFractionProton(float mom) const;
72 ClassDef(AliFemtoAODTrackCut, 1)
77 inline void AliFemtoAODTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
78 inline void AliFemtoAODTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
79 inline void AliFemtoAODTrackCut::SetCharge(const int& ch){fCharge = ch;}
80 inline void AliFemtoAODTrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
81 inline void AliFemtoAODTrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
82 inline void AliFemtoAODTrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
83 inline void AliFemtoAODTrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
84 inline void AliFemtoAODTrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
85 inline void AliFemtoAODTrackCut::SetLabel(const bool& flag){fLabel=flag;}
86 inline void AliFemtoAODTrackCut::SetMostProbablePion() { fMostProbable = 2; }
87 inline void AliFemtoAODTrackCut::SetMostProbableKaon() { fMostProbable = 3; }
88 inline void AliFemtoAODTrackCut::SetMostProbableProton() { fMostProbable = 4; }
89 inline void AliFemtoAODTrackCut::SetNoMostProbable() { fMostProbable = 0; }
90 inline void AliFemtoAODTrackCut::SetMaxChiNdof(const float& maxchi) { fMaxchiNdof = maxchi; }
91 inline void AliFemtoAODTrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }