PWGCF/FEMTOSCOPY/AliFemtoUser/AliFemtoQATrackCut.h

   1 ///////////////////////////////////////////////////////////////////////////
   2 //                                                                       //
   3 // AliFemtoQATrackCut: 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 AliFemtoQATrackCut_H
  12 #define AliFemtoQATrackCut_H
  13
  14 //#ifndef StMaker_H
  15 //#include "StMaker.h"
  16 //#endif
  17
  18 #include "AliFemtoTrackCut.h"
  19
  20 class AliFemtoQATrackCut : public AliFemtoTrackCut
  21 {
  22
  23  public:
  24   AliFemtoQATrackCut();
  25   virtual ~AliFemtoQATrackCut();
  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 SetStatus(const long& w);
  43   void SetminTPCclsF(const short& s);
  44   void SetminTPCncls(const short& s);
  45   void SetminITScls(const int& s);
  46   void SetminTPCchiNdof(const float& s);
  47   void SetRemoveKinks(const bool& flag);
  48   void SetMaxTPCncls(const short& s);
  49   void SetMaxITSChiNdof(const float& maxchi);
  50   void SetMaxTPCChiNdof(const float& maxchi);
  51   void SetMaxSigmaToVertex(const float& maxsig);
  52   void SetMostProbablePion();
  53   void SetMostProbableKaon();
  54   void SetMostProbableProton();
  55   void SetNoMostProbable();
  56
  57   void SetTPCnclsExclusionSwitch(const bool& excluSwitch);
  58   void SetTPCnclsExclusion(const short& lo, const short& hi);
  59   void SetTPCchiNdofExclusionSwitch(const bool& excluSwitch);
  60   void SetTPCchiNdofExclusion(const float& lo, const float& hi);
  61
  62  private:   // here are the quantities I want to cut on...
  63
  64   int               fCharge;             // particle charge
  65   float             fPt[2];              // bounds for transverse momentum
  66   float             fRapidity[2];        // bounds for rapidity
  67   float             fPidProbElectron[2]; // bounds for electron probability
  68   float             fPidProbPion[2];     // bounds for pion probability
  69   float             fPidProbKaon[2];     // bounds for kaon probability
  70   float             fPidProbProton[2];   // bounds for proton probability
  71   float             fPidProbMuon[2];     // bounds for muon probability
  72   bool              fLabel;              // if true label<0 will not pass throught
  73   long              fStatus;             // staus flag
  74
  75   short             fminTPCclsF;         // min number of findable clusters in the TPC
  76   short             fminTPCncls;         // min number of clusters in the TPC
  77   int               fminITScls;          // min number of clusters assigned in the ITS
  78   float             fminTPCchiNdof;      // min allowed chi2/ndof for TPC clusters
  79   short             fMaxTPCncls;         // maximum allowed clusters in the TPC
  80   float             fMaxITSchiNdof;      // maximum allowed chi2/ndof for ITS clusters
  81   float             fMaxTPCchiNdof;      // maximum allowed chi2/ndof for TPC clusters
  82   float             fMaxSigmaToVertex;   // maximum allowed sigma to primary vertex
  83   long              fNTracksPassed;      // passed tracks count
  84   long              fNTracksFailed;      // failed tracks count
  85   bool              fRemoveKinks;        // if true particles with any kink label will not pass
  86   int               fMostProbable;       // this particle type is required to be most probable
  87
  88   bool    fTPCnclsExclusionSwitch;       // turn on/off TPCncls exclusion zone (true=on)
  89   short   fTPCnclsExclusion[2];          // lower and upper limit of TPCncls QA exclusion zone
  90   bool    fTPCchiNdofExclusionSwitch;        // turn on/off TPCchi exclusion zone (true=on)
  91   float   fTPCchiNdofExclusion[2];           // lower and upper limit of TPCchi QA exclusion zone
  92
  93   float PidFractionElectron(float mom) const;
  94   float PidFractionPion(float mom) const;
  95   float PidFractionKaon(float mom) const;
  96   float PidFractionProton(float mom) const;
  97
  98 #ifdef __ROOT__
  99   ClassDef(AliFemtoQATrackCut, 1)
 100 #endif
 101     };
 102
 103
 104 inline void AliFemtoQATrackCut::SetPt(const float& lo, const float& hi){fPt[0]=lo; fPt[1]=hi;}
 105 inline void AliFemtoQATrackCut::SetRapidity(const float& lo,const float& hi){fRapidity[0]=lo; fRapidity[1]=hi;}
 106 inline void AliFemtoQATrackCut::SetCharge(const int& ch){fCharge = ch;}
 107 inline void AliFemtoQATrackCut::SetPidProbElectron(const float& lo,const float& hi){fPidProbElectron[0]=lo; fPidProbElectron[1]=hi;}
 108 inline void AliFemtoQATrackCut::SetPidProbPion(const float& lo,const float& hi){fPidProbPion[0]=lo; fPidProbPion[1]=hi;}
 109 inline void AliFemtoQATrackCut::SetPidProbKaon(const float& lo,const float& hi){fPidProbKaon[0]=lo; fPidProbKaon[1]=hi;}
 110 inline void AliFemtoQATrackCut::SetPidProbProton(const float& lo,const float& hi){fPidProbProton[0]=lo; fPidProbProton[1]=hi;}
 111 inline void AliFemtoQATrackCut::SetPidProbMuon(const float& lo,const float& hi){fPidProbMuon[0]=lo; fPidProbMuon[1]=hi;}
 112 inline void AliFemtoQATrackCut::SetLabel(const bool& flag){fLabel=flag;}
 113 inline void AliFemtoQATrackCut::SetStatus(const long& status){fStatus=status;}
 114 inline void AliFemtoQATrackCut::SetminTPCclsF(const short& minTPCclsF){fminTPCclsF=minTPCclsF;}
 115 inline void AliFemtoQATrackCut::SetminTPCncls(const short& s){fminTPCncls=s;}
 116 inline void AliFemtoQATrackCut::SetminITScls(const int& minITScls){fminITScls=minITScls;}
 117 inline void AliFemtoQATrackCut::SetminTPCchiNdof(const float& s){fminTPCchiNdof = s;}
 118 inline void AliFemtoQATrackCut::SetMostProbablePion() { fMostProbable = 2; }
 119 inline void AliFemtoQATrackCut::SetMostProbableKaon() { fMostProbable = 3; }
 120 inline void AliFemtoQATrackCut::SetMostProbableProton() { fMostProbable = 4; }
 121 inline void AliFemtoQATrackCut::SetNoMostProbable() { fMostProbable = 0; }
 122 inline void AliFemtoQATrackCut::SetMaxTPCncls(const short& s){fMaxTPCncls=s;}
 123 inline void AliFemtoQATrackCut::SetMaxITSChiNdof(const float& maxchi) { fMaxITSchiNdof = maxchi; }
 124 inline void AliFemtoQATrackCut::SetMaxTPCChiNdof(const float& maxchi) { fMaxTPCchiNdof = maxchi; }
 125 inline void AliFemtoQATrackCut::SetMaxSigmaToVertex(const float& maxsig) { fMaxSigmaToVertex = maxsig; }
 126
 127 inline void AliFemtoQATrackCut::SetTPCnclsExclusionSwitch(const bool& excluSwitch) { fTPCnclsExclusionSwitch = excluSwitch; }
 128 inline void AliFemtoQATrackCut::SetTPCnclsExclusion(const short& lo, const short& hi) {fTPCnclsExclusion[0] = lo; fTPCnclsExclusion[1] = hi;}
 129 inline void AliFemtoQATrackCut::SetTPCchiNdofExclusionSwitch(const bool& excluSwitch) { fTPCchiNdofExclusionSwitch = excluSwitch; }
 130 inline void AliFemtoQATrackCut::SetTPCchiNdofExclusion(const float& lo, const float& hi) {fTPCchiNdofExclusion[0] = lo; fTPCchiNdofExclusion[1] = hi;}
 131
 132 #endif
 133