PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoKKTrackCut.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 ALIFEMTOKKTRACKCUT_H
  13 #define ALIFEMTOKKTRACKCUT_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 AliFemtoKKTrackCut : public AliFemtoTrackCut
  24 {
  25   public:
  26
  27   enum PIDMethodType {knSigma=0, kContour=1};
  28   typedef enum PIDMethodType ReadPIDMethodType;
  29
  30   AliFemtoKKTrackCut();
  31   virtual ~AliFemtoKKTrackCut();
  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   //ml
  73  void SetNsigmaTPCle250(Double_t);
  74  void SetNsigmaTPC250_400(Double_t);
  75  void SetNsigmaTPC400_450(Double_t);
  76  void SetNsigmaTPC450_500(Double_t);
  77  void SetNsigmaTPCge500(Double_t);
  78
  79  void SetNsigmaTOF500_800(Double_t);
  80  void SetNsigmaTOF800_1000(Double_t);
  81  void SetNsigmaTOFge1000(Double_t);
  82
  83
  84
  85  private:   // here are the quantities I want to cut on...
  86
  87   int               fCharge;             // particle charge
  88   float             fPt[2];              // bounds for transverse momentum
  89   float             fRapidity[2];        // bounds for rapidity
  90   float             fEta[2];             // bounds for pseudorapidity
  91   float             fPidProbElectron[2]; // bounds for electron probability
  92   float             fPidProbPion[2];     // bounds for pion probability
  93   float             fPidProbKaon[2];     // bounds for kaon probability
  94   float             fPidProbProton[2];   // bounds for proton probability
  95   float             fPidProbMuon[2];     // bounds for muon probability
  96
  97   AliESDtrackCuts::ITSClusterRequirement fCutClusterRequirementITS[3];  // detailed ITS cluster requirements for (SPD, SDD, SSD) - from AliESDtrackcuts!
  98   bool              fLabel;              // if true label<0 will not pass throught
  99   long              fStatus;             // staus flag
 100   ReadPIDMethodType fPIDMethod;          // which PID mehod to use. 0 - nsgima, 1 - contour
 101
 102   //ml
 103  Double_t fNsigmaTPCle250;
 104  Double_t fNsigmaTPC250_400;
 105  Double_t fNsigmaTPC400_450;
 106  Double_t fNsigmaTPC450_500;
 107  Double_t fNsigmaTPCge500;
 108
 109  Double_t fNsigmaTOF500_800;
 110  Double_t fNsigmaTOF800_1000;
 111  Double_t fNsigmaTOFge1000;
 112
 113
 114
 115   short             fminTPCclsF;         // min number of findable clusters in the TPC
 116   short             fminTPCncls;         // min number of clusters in the TPC
 117   int               fminITScls;          // min number of clusters assigned in the ITS
 118   float             fMaxITSchiNdof;      // maximum allowed chi2/ndof for ITS clusters
 119   float             fMaxTPCchiNdof;      // maximum allowed chi2/ndof for TPC clusters
 120   float             fMaxSigmaToVertex;   // maximum allowed sigma to primary vertex
 121   long              fNTracksPassed;      // passed tracks count
 122   long              fNTracksFailed;      // failed tracks count
 123   bool              fRemoveKinks;        // if true particles with any kink label will not pass
 124   bool              fRemoveITSFake;      // if true particles with ITS fake flag will not pass
 125   int               fMostProbable;       // this particle type is required to be most probable
 126
 127   float             fMaxImpactXY;        // Max XY impact parameter
 128   float             fMaxImpactZ;         // Max Z impact parameter
 129
 130   float             fMaxImpactXYPtOff;   // Max XY DCA Pt dependent offset
 131   float             fMaxImpactXYPtNrm;   // Max XY DCA Pt dependent normalization
 132   float             fMaxImpactXYPtPow;   // Max XY DCA Pt dependent power
 133
 134   float             fMinPforTOFpid;  // momentum from which TOF PID is requested
 135   float             fMaxPforTOFpid;  // momentum till which TOF PID is requested
 136   float             fMinPforTPCpid;  // momentum from which TPC PID is requested
 137   float             fMaxPforTPCpid;  // momentum till which TPC PID is requested
 138   float             fMinPforITSpid;  // momentum from which ITS PID is requested
 139   float             fMaxPforITSpid;  // momentum till which ITS PID is requested
 140
 141   float PidFractionElectron(float mom) const;
 142   float PidFractionPion(float mom) const;
 143   float PidFractionKaon(float mom) const;
 144   float PidFractionProton(float mom) const;
 145
 146   bool IsPionTPCdEdx(float mom, float dEdx);
 147   bool IsKaonTPCdEdx(float mom, float dEdx);
 148   bool IsProtonTPCdEdx(float mom, float dEdx);
 149
 150   bool IsPionTOFTime(float mom, float ttof);
 151   bool IsKaonTOFTime(float mom, float ttof);
 152   bool IsProtonTOFTime(float mom, float ttof);
 153
 154   bool IsKaonTPCdEdxNSigma(float mom, float nsigma);
 155   bool IsKaonTOFNSigma(float mom, float nsigma);
 156   bool IsKaonNSigma(float mom, float nsigmaTPC, float nsigmaTOF);
 157   bool IsPionNSigma(float mom, float nsigmaTPC, float nsigmaTOF);
 158   bool IsProtonNSigma(float mom, float nsigmaTPC, float nsigmaTOF);
 159
 160   Bool_t CheckITSClusterRequirement(AliESDtrackCuts::ITSClusterRequirement req, Bool_t clusterL1, Bool_t clusterL2); //the same as in AliESDtrackCuts
 161
 162
 163 #ifdef __ROOT__
 164   ClassDef(AliFemtoKKTrackCut, 1)
 165 #endif
 166     };
 167
 168
 169 inline void AliFemtoKKTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
 170 inline void AliFemtoKKTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
 171 inline void AliFemtoKKTrackCut::SetEta(const float& lo,const float& hi){fEta[0]=lo; fEta[1]=hi;}
 172 inline void AliFemtoKKTrackCut::SetCharge(const int& ch){fCharge = ch;}
 173 inline void AliFemtoKKTrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
 174 inline void AliFemtoKKTrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
 175 inline void AliFemtoKKTrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
 176 inline void AliFemtoKKTrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
 177 inline void AliFemtoKKTrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
 178 inline void AliFemtoKKTrackCut::SetLabel(const bool& flag){fLabel=flag;}
 179 inline void AliFemtoKKTrackCut::SetStatus(const long& status){fStatus=status;}
 180 inline void AliFemtoKKTrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;}
 181 inline void AliFemtoKKTrackCut::SetminTPCncls(const short& s){fminTPCncls=s;}
 182 inline void AliFemtoKKTrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;}
 183 inline void AliFemtoKKTrackCut::SetMostProbablePion() { fMostProbable = 2; }
 184 inline void AliFemtoKKTrackCut::SetMostProbableKaon() { fMostProbable = 3; }
 185 inline void AliFemtoKKTrackCut::SetMostProbableProton() { fMostProbable = 4; }
 186 inline void AliFemtoKKTrackCut::SetNoMostProbable() { fMostProbable = 0; }
 187 inline void AliFemtoKKTrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; }
 188 inline void AliFemtoKKTrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; }
 189 inline void AliFemtoKKTrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }
 190 inline void AliFemtoKKTrackCut::SetMaxImpactXY(const float& maximpxy) { fMaxImpactXY = maximpxy; }
 191 inline void AliFemtoKKTrackCut::SetMaxImpactXYPtDep(const float& maxoff, const float& maxnrm, const float& maxpow) { fMaxImpactXYPtOff = maxoff; fMaxImpactXYPtNrm = maxnrm; fMaxImpactXYPtPow = maxpow; }
 192 inline void AliFemtoKKTrackCut::SetMaxImpactZ(const float& maximpz) { fMaxImpactZ = maximpz; }
 193
 194
 195
 196
 197 #endif
 198