ITS/AliITSPident.h

   1 #ifndef ALIITSPIDENT_H
   2 #define ALIITSPIDENT_H
   3
   4 /////////////////////////////////////////////////////////////////////////
   5 //Class for PID in the ITS                                               //
   6 //The PID is based on the likelihood of all the four ITS' layers,        //
   7 //without using the truncated mean for the dE/dx. The response           //
   8 //functions for each layer are convoluted Landau-Gaussian functions.     //
   9 // Origin: Elena Bruna bruna@to.infn.it, Massimo Masera masera@to.infn.it//
  10 ////////////////////////////////////////////////////////////////////////
  11
  12 #include <TObject.h>
  13 class AliITStrackV2;
  14 class AliITSSteerPid;
  15 class TF1;
  16 class AliITSPidParItem;
  17 class AliITSPident : public TObject{
  18
  19  public:
  20   AliITSPident();
  21   AliITSPident(Double_t mom,Double_t dEdx,AliITSSteerPid *sp,Float_t *Qlay,Float_t *nlay,Float_t priorip=0.066,Float_t priorik=0.103,Float_t prioripi=0.83,Float_t priorie=0.001);
  22
  23   AliITSPident(AliITStrackV2 *trackITS,AliITSSteerPid *sp,Float_t *Qlay,Float_t *nlay,Float_t priorip=0.066,Float_t priorik=0.103,Float_t prioripi=0.83,Float_t priorie=0.001);
  24
  25   virtual ~AliITSPident();
  26   Float_t GetP() const {return fMom;}//local momentum (GeV/c)
  27   Double_t GetdEdx() const {return fdEdx;} //dEdx (ADC units)
  28   Double_t GetCondFunPro(Int_t lay) const {
  29     return fCondFunProLay[lay];
  30   }
  31   Double_t GetProdCondFunPro() const;
  32   Double_t GetCondFunK(Int_t lay) const {
  33     return fCondFunKLay[lay];
  34   }
  35   Double_t GetProdCondFunK() const;
  36   Double_t GetCondFunPi(Int_t lay) const {
  37     return fCondFunPiLay[lay];
  38   }
  39   Double_t GetProdCondFunPi() const;
  40   void PrintParameters() const;
  41   Float_t GetPBayesp()const {return fPBayesp;}
  42   Float_t GetPBayesk()const {return fPBayesk;}
  43   Float_t GetPBayespi()const {return fPBayespi;}
  44   Float_t GetPPriorip() const {return fPPriorip;}
  45   Float_t GetPPriorik() const {return fPPriorik;}
  46   Float_t GetPPrioripi() const {return fPPrioripi;}
  47   Float_t GetPPriorie() const {return fPPriorie;}
  48   void GetNclsPerLayer(Int_t *ncls) const;
  49   static Double_t Langaufun(Double_t *x, Double_t *par);
  50   static Double_t Langaufun2(Double_t *x, Double_t *par);
  51   static Double_t Langaufunnorm(Double_t *x, Double_t *par);
  52   static Double_t Gaus2(Double_t *x, Double_t *par);
  53
  54  private:
  55
  56   void CookFunItsLay(Int_t lay,Int_t opt,Double_t *par,Double_t dedx,Double_t mom,Double_t rangei,Double_t rangef,TString comment);
  57   void CookBayes(Double_t *condfun,Float_t *prior);
  58   Float_t CookCombinedBayes(Double_t condfun[][3],Float_t *prior,Int_t part)const;
  59   Float_t CookProd(Double_t condfun[][3],Int_t part) const;
  60   Float_t CookSum(Double_t condfun[][3],Float_t *prior) const;
  61   AliITSPident(const AliITSPident &ob); // copy constructor
  62   AliITSPident& operator=(const AliITSPident & ob); // ass. op.
  63
  64   Float_t fMom;                   // Particle momentum
  65   Double_t fdEdx;                 // Particle dE/dx
  66   Double_t fCondFunProLay[8];     // one for each silicon layer
  67   Double_t fCondFunKLay[8];       // cond. prob. function kaons per layer
  68   Double_t fCondFunPiLay[8];      // cond. prob. function pions per layer
  69   Float_t fPBayesp;               // Bayes prob.
  70   Float_t fPBayesk;               // Bayes prob. for kaons
  71   Float_t fPBayespi;              // Bayes prob. for pions
  72   Float_t fPPriorip;              // Priori prob.
  73   Float_t fPPriorik;              // Priori prob. for kaons
  74   Float_t fPPrioripi;             // Priori prob. for pions
  75   Float_t fPPriorie;              // Priori prob. for electrons
  76   Int_t fNcls[4];                 // N. of clusters per layer (sdd1,sdd2,ssd1,ssd2)
  77   ClassDef(AliITSPident,2);
  78 };
  79 #endif