HMPID/AliHMPIDHit.h

   1 #ifndef AliHMPIDHit_h
   2 #define AliHMPIDHit_h
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 #include <AliHit.h>           //base class
   7 #include <TVector3.h>         //ctor
   8 #include <TRandom.h>          //QdcTot()
   9 #include <TMath.h>
  10
  11 class AliHMPIDHit : public AliHit //   TObject-AliHit-AliHMPIDHit
  12 {
  13 public:
  14   AliHMPIDHit(                                                                             ):AliHit(     ),fCh(-1),fPid(-1  ),fQ(-1),fLorsX(-1),fLorsY(-1) {} //default ctor
  15   AliHMPIDHit(Int_t c,Float_t e,Int_t pid,Int_t tid,Float_t xl,Float_t yl,const TVector3 &p):AliHit(0,tid),fCh(c ),fPid(pid ),fQ(-1),fLorsX(xl),fLorsY(yl) {QdcTot(e);fX=p.X();fY=p.Y();fZ=p.Z();}
  16   AliHMPIDHit(Int_t c,Int_t q  ,                    Float_t xl,Float_t yl                  ):              fCh(c ),fPid(2212),fQ(q ),fLorsX(xl),fLorsY(yl) {fTrack=Int_t(1000*gRandom->Rndm());}//manual ctor
  17   virtual ~AliHMPIDHit()                                                                                                {}
  18 //framework part
  19   void     Print(Option_t *option="")const;                                    //from TObject to print current status
  20 //private part
  21          Int_t   Ch    (         )const{return fCh;                                    }       //Chamber
  22          Float_t LorsX (         )const{return fLorsX;                                 }       //hit X position in LORS, [cm]
  23          Float_t LorsY (         )const{return fLorsY;                                 }       //hit Y position in LORS, [cm]
  24          Int_t   Pid   (         )const{return fPid;                                   }       //PID
  25          Float_t Q     (         )const{return fQ;                                     }       //Eloss for MIP hit or Etot for photon hit, [GeV]
  26   inline Float_t QdcTot(Float_t e);                                                            //calculate total charge of the hit
  27          Int_t   Tid   (         )const{return fTrack;                                 }       //TID
  28
  29 protected:                                                                     //AliHit has fTid,fX,fY,fZ
  30   Int_t    fCh;                                                                //Chamber
  31   Int_t    fPid;                                                               //PID
  32   Float_t  fQ;                                                                 //total charge [QDC]
  33   Float_t  fLorsX;                                                             //hit X position in chamber LORS, [cm]
  34   Float_t  fLorsY;                                                             //hit Y position in chamber LORS, [cm]
  35   ClassDef(AliHMPIDHit,5)                                                      //HMPID hit class
  36 };//class AliHMPIDhit
  37 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  38 Float_t AliHMPIDHit::QdcTot(Float_t e)
  39 {
  40 // Samples total charge of the hit
  41 // Arguments: e- hit energy [GeV]; for mip Eloss for photon Etot
  42 //   Returns: total QDC
  43   Float_t  x=(LorsX() > 66.6)? LorsX()-66.6:LorsX();                                                 //sagita is for PC (0-64) and not for chamber
  44   Float_t  qdcEle=34.06311+0.2337070*x+5.807476e-3*x*x-2.956471e-04*x*x*x+2.310001e-06*x*x*x*x;      //reparametrised from DiMauro
  45
  46   Int_t iNele=Int_t(e/26e-9);  if(iNele<1) iNele=1;                                                  //number of electrons created by hit
  47   for(Int_t i=1;i<=iNele;i++) fQ-=qdcEle*TMath::Log(gRandom->Rndm()+1e-6);                           //1e-6 is a protection against 0 from rndm
  48   return fQ;
  49 }
  50
  51
  52
  53 typedef AliHMPIDHit AliRICHHit; // for backward compatibility
  54
  55 #endif