Bring AliFemto up to date with latest code developements
[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemtoUser / AliFemtoESDTrackCut.h
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 ///////////////////////////////////////////////////////////////////////////
10
11 #ifndef ALIFEMTOESDTRACKCUT_H
12 #define ALIFEMTOESDTRACKCUT_H
13
14 //#ifndef StMaker_H
15 //#include "StMaker.h"
16 //#endif
17
18 #include "AliFemtoTrackCut.h"
19
20 class AliFemtoESDTrackCut : public AliFemtoTrackCut 
21 {
22
23  public:
24   AliFemtoESDTrackCut();
25   virtual ~AliFemtoESDTrackCut();
26
27   virtual bool Pass(const AliFemtoTrack* aTrack);
28
29   virtual AliFemtoString Report();
30   virtual TList *ListSettings();
31   virtual AliFemtoParticleType Type(){return hbtTrack;}
32
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 SetStatus(const long& w);
43   void SetminTPCclsF(const short& s);
44   void SetminTPCncls(const short& s);
45   void SetminITScls(const int& s);
46   void SetRemoveKinks(const bool& flag);
47   void SetMaxITSChiNdof(const float& maxchi);
48   void SetMaxTPCChiNdof(const float& maxchi);
49   void SetMaxSigmaToVertex(const float& maxsig);
50   void SetMostProbablePion();
51   void SetMostProbableKaon();
52   void SetMostProbableProton();
53   void SetNoMostProbable(); 
54
55  private:   // here are the quantities I want to cut on...
56
57   int               fCharge;             // particle charge
58   float             fPt[2];              // bounds for transverse momentum
59   float             fRapidity[2];        // bounds for rapidity
60   float             fPidProbElectron[2]; // bounds for electron probability
61   float             fPidProbPion[2];     // bounds for pion probability
62   float             fPidProbKaon[2];     // bounds for kaon probability
63   float             fPidProbProton[2];   // bounds for proton probability
64   float             fPidProbMuon[2];     // bounds for muon probability 
65   bool              fLabel;              // if true label<0 will not pass throught 
66   long              fStatus;             // staus flag
67   short             fminTPCclsF;         // min number of findable clusters in the TPC
68   short             fminTPCncls;         // min number of clusters in the TPC
69   int               fminITScls;          // min number of clusters assigned in the ITS 
70   float             fMaxITSchiNdof;      // maximum allowed chi2/ndof for ITS clusters
71   float             fMaxTPCchiNdof;      // maximum allowed chi2/ndof for TPC clusters
72   float             fMaxSigmaToVertex;   // maximum allowed sigma to primary vertex
73   long              fNTracksPassed;      // passed tracks count
74   long              fNTracksFailed;      // failed tracks count
75   bool              fRemoveKinks;        // if true particles with any kink label will not pass
76   int               fMostProbable;       // this particle type is required to be most probable
77
78   float PidFractionElectron(float mom) const;
79   float PidFractionPion(float mom) const;
80   float PidFractionKaon(float mom) const;
81   float PidFractionProton(float mom) const;
82
83 #ifdef __ROOT__ 
84   ClassDef(AliFemtoESDTrackCut, 1)
85 #endif
86     };
87
88
89 inline void AliFemtoESDTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
90 inline void AliFemtoESDTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
91 inline void AliFemtoESDTrackCut::SetCharge(const int& ch){fCharge = ch;}
92 inline void AliFemtoESDTrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
93 inline void AliFemtoESDTrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
94 inline void AliFemtoESDTrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
95 inline void AliFemtoESDTrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
96 inline void AliFemtoESDTrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
97 inline void AliFemtoESDTrackCut::SetLabel(const bool& flag){fLabel=flag;}
98 inline void AliFemtoESDTrackCut::SetStatus(const long& status){fStatus=status;}
99 inline void AliFemtoESDTrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;}
100 inline void AliFemtoESDTrackCut::SetminTPCncls(const short& s){fminTPCncls=s;}
101 inline void AliFemtoESDTrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;}
102 inline void AliFemtoESDTrackCut::SetMostProbablePion() { fMostProbable = 2; }
103 inline void AliFemtoESDTrackCut::SetMostProbableKaon() { fMostProbable = 3; }
104 inline void AliFemtoESDTrackCut::SetMostProbableProton() { fMostProbable = 4; }
105 inline void AliFemtoESDTrackCut::SetNoMostProbable() { fMostProbable = 0; }
106 inline void AliFemtoESDTrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; }
107 inline void AliFemtoESDTrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; }
108 inline void AliFemtoESDTrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }
109
110 #endif
111