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
  32   void SetPt(const float& lo, const float& hi);
  33   void SetRapidity(const float& lo, const float& hi);
  34   void SetCharge(const int& ch);
  35   void SetPidProbElectron(const float& lo, const float& hi);
  36   void SetPidProbPion(const float& lo, const float& hi);
  37   void SetPidProbKaon(const float& lo, const float& hi);
  38   void SetPidProbProton(const float& lo, const float& hi);
  39   void SetPidProbMuon(const float& lo, const float& hi);
  40   void SetLabel(const bool& flag);
  41   void SetStatus(const long& w);
  42   void SetminTPCclsF(const short& s);
  43   void SetminITScls(const int& s);
  44   void SetRemoveKinks(const bool& flag);
  45   void SetMostProbablePion();
  46   void SetMostProbableKaon();
  47   void SetMostProbableProton();
  48   void SetNoMostProbable();
  49
  50  private:   // here are the quantities I want to cut on...
  51
  52   int               fCharge;             // particle charge
  53   float             fPt[2];              // bounds for transverse momentum
  54   float             fRapidity[2];        // bounds for rapidity
  55   float             fPidProbElectron[2]; // bounds for electron probability
  56   float             fPidProbPion[2];     // bounds for pion probability
  57   float             fPidProbKaon[2];     // bounds for kaon probability
  58   float             fPidProbProton[2];   // bounds for proton probability
  59   float             fPidProbMuon[2];     // bounds for muon probability
  60   bool              fLabel;              // if true label<0 will not pass throught
  61   long              fStatus;             // staus flag
  62   short             fminTPCclsF;         // min number of findable clusters in the TPC
  63   int               fminITScls;          // min number of clusters assigned in the ITS
  64   long              fNTracksPassed;      // passed tracks count
  65   long              fNTracksFailed;      // failed tracks count
  66   bool              fRemoveKinks;        // if true particles with any kink label will not pass
  67   int               fMostProbable;       // this particle type is required to be most probable
  68
  69   float PidFractionElectron(float mom) const;
  70   float PidFractionPion(float mom) const;
  71   float PidFractionKaon(float mom) const;
  72   float PidFractionProton(float mom) const;
  73
  74 #ifdef __ROOT__
  75   ClassDef(AliFemtoESDTrackCut, 1)
  76 #endif
  77     };
  78
  79
  80 inline void AliFemtoESDTrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
  81 inline void AliFemtoESDTrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
  82 inline void AliFemtoESDTrackCut::SetCharge(const int& ch){fCharge = ch;}
  83 inline void AliFemtoESDTrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
  84 inline void AliFemtoESDTrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
  85 inline void AliFemtoESDTrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
  86 inline void AliFemtoESDTrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
  87 inline void AliFemtoESDTrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
  88 inline void AliFemtoESDTrackCut::SetLabel(const bool& flag){fLabel=flag;}
  89 inline void AliFemtoESDTrackCut::SetStatus(const long& status){fStatus=status;}
  90 inline void AliFemtoESDTrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;}
  91 inline void AliFemtoESDTrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;}
  92 inline void AliFemtoESDTrackCut::SetMostProbablePion() { fMostProbable = 2; }
  93 inline void AliFemtoESDTrackCut::SetMostProbableKaon() { fMostProbable = 3; }
  94 inline void AliFemtoESDTrackCut::SetMostProbableProton() { fMostProbable = 4; }
  95 inline void AliFemtoESDTrackCut::SetNoMostProbable() { fMostProbable = 0; }
  96
  97 #endif
  98