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