STEER/STEERBase/AliITSPIDResponse.h

   1 #ifndef ALIITSPIDRESPONSE_H
   2 #define ALIITSPIDRESPONSE_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 //-------------------------------------------------------
   9 //                    ITS PID response class
  10 //
  11 //
  12 //-------------------------------------------------------
  13 //#include <Rtypes.h>
  14 #include <TObject.h>
  15 #include "AliPID.h"
  16
  17 class AliITSPIDResponse : public TObject {
  18
  19 public:
  20   AliITSPIDResponse(Bool_t isMC=kFALSE);
  21   //AliITSPIDResponse(Double_t *param);
  22  ~AliITSPIDResponse() {}
  23
  24  void SetBetheBlochParamsITSTPC(Double_t* param){
  25    for(Int_t iPar=0; iPar<5; iPar++) fBBtpcits[iPar]=param[iPar];
  26  }
  27  void SetBetheBlochParamsITSsa(Double_t* param){
  28    for(Int_t iPar=0; iPar<5; iPar++) fBBsa[iPar]=param[iPar];
  29  }
  30  void SetElectronBetheBlochParamsITSsa(Double_t* param){
  31    for(Int_t iPar=0; iPar<5; iPar++) fBBsaElectron[iPar]=param[iPar];
  32  }
  33
  34  Double_t BetheAleph(Double_t p,Double_t mass) const;
  35  Double_t Bethe(Double_t p, Double_t mass, Bool_t iSA=kFALSE) const;
  36  Double_t GetResolution(Double_t bethe, Int_t nPtsForPid=4, Bool_t isSA=kFALSE) const;
  37  void GetITSProbabilities(Float_t mom, Double_t qclu[4], Double_t condprobfun[AliPID::kSPECIES],Bool_t isMC=kFALSE) const;
  38  Float_t GetNumberOfSigmas(Float_t mom, Float_t signal, AliPID::EParticleType type, Int_t nPtsForPid=4, Bool_t isSA=kFALSE) const {
  39    const Double_t chargeFactor = TMath::Power(AliPID::ParticleCharge(type),2.);
  40    Float_t bethe = Bethe(mom,AliPID::ParticleMassZ(type),isSA)*chargeFactor;
  41    return (signal - bethe)/GetResolution(bethe,nPtsForPid,isSA);
  42  }
  43  Int_t GetParticleIdFromdEdxVsP(Float_t mom, Float_t signal, Bool_t isSA=kFALSE) const;
  44
  45 private:
  46
  47
  48   // Data members for truncated mean method
  49   Float_t  fRes;             // relative dEdx resolution
  50   Double_t fKp1;             // ALEPH BB param 1
  51   Double_t fKp2;             // ALEPH BB param 2
  52   Double_t fKp3;             // ALEPH BB param 3
  53   Double_t fKp4;             // ALEPH BB param 4
  54   Double_t fKp5;             // ALEPH BB param
  55   Double_t  fBBsa[5];         // parameters of BB for SA tracks
  56   Double_t  fBBsaElectron[5];         // parameters of BB for SA tracks
  57   Double_t  fBBtpcits[5];     // parameters of BB for TPC+ITS tracks
  58   Float_t  fResolSA[5];      // resolutions vs. n. of SDD/SSD points
  59   Float_t  fResolTPCITS[5];  // resolutions vs. n. of SDD/SSD points
  60
  61   ClassDef(AliITSPIDResponse,3)   // ITS PID class
  62 };
  63
  64 #endif
  65
  66