PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoMJTrackCut.h

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