4d7454cebb378c0880d4b6b0811dad4ae3b9c3a5
[u/mrichter/AliRoot.git] / PWG2 / FEMTOSCOPY / AliFemto / AliFemtoAODTrackCut.h
1 ///////////////////////////////////////////////////////////////////////////
2 //                                                                       //
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                   //
8 //                                                                       //
9 ///////////////////////////////////////////////////////////////////////////
10
11 #ifndef ALIFEMTOAODTRACKCUT_H
12 #define ALIFEMTOAODTRACKCUT_H
13
14 //#ifndef StMaker_H
15 //#include "StMaker.h"
16 //#endif
17
18 #include "AliFemtoTrackCut.h"
19
20 class AliFemtoAODTrackCut : public AliFemtoTrackCut 
21 {
22
23  public:
24   AliFemtoAODTrackCut();
25   virtual ~AliFemtoAODTrackCut();
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 SetMaxChiNdof(const float& maxchi);
43   void SetMaxSigmaToVertex(const float& maxsig);
44   void SetMostProbablePion();
45   void SetMostProbableKaon();
46   void SetMostProbableProton();
47   void SetNoMostProbable(); 
48
49  private:   // here are the quantities I want to cut on...
50
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
65
66   float PidFractionElectron(float mom) const;
67   float PidFractionPion(float mom) const;
68   float PidFractionKaon(float mom) const;
69   float PidFractionProton(float mom) const;
70
71 #ifdef __ROOT__ 
72   ClassDef(AliFemtoAODTrackCut, 1)
73 #endif
74     };
75
76
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; }
92
93 #endif
94