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 // HMPID base class to produce hits
   7 //.
   8 //.
   9 #include <AliHit.h>           //base class
  10 #include <TVector3.h>         //ctor
  11 #include "AliHMPIDDigit.h"    //QdcTot()
  12
  13 #include <TRandom.h>
  14
  15 class AliHMPIDHit : public AliHit //   TObject-AliHit-AliHMPIDHit
  16 {
  17 public:
  18   AliHMPIDHit(                                                                           ):AliHit(     ),fCh(-1),fPid(-1 ),fQ(-1),fLx(0),fLy(0){} //default ctor
  19   AliHMPIDHit(Int_t c,Float_t &e,Int_t pid,Int_t tid,Float_t x,Float_t y,const TVector3 &p):AliHit(0,tid),fCh(c ),fPid(pid),fQ(0 ),fLx(x),fLy(y){e=QdcTot(e);fX=p.X();fY=p.Y();fZ=p.Z();}
  20   AliHMPIDHit(Int_t c,Float_t &e,Int_t pid,Int_t tid,Float_t x,Float_t y                  ):AliHit(     ),fCh(c ),fPid(pid),fQ(0 ),fLx(x),fLy(y){e=QdcTot(e);fTrack=tid;}//manual ctor
  21   AliHMPIDHit(const AliHMPIDHit &h):AliHit(h),fCh(h.fCh),fPid(h.fPid),fQ(h.fQ),fLx(h.fLx),fLy(h.fLy) {}//copy ctor
  22   virtual ~AliHMPIDHit()                                                                                                                                 {}
  23 //framework part
  24          void    Print(Option_t *opt="")const;                                                    //from TObject to print current status
  25          void    Draw (Option_t *opt="");                                                         //from TObject to Draw this hit
  26 //private part
  27          Int_t   Ch     (                               )const{return fCh;                                    }       //Chamber
  28          void    Hit2Sdi(TClonesArray *pSdiLst,Int_t n=1)const;                                                       //add sdigits of this hit to the list
  29          Float_t LorsX  (                               )const{return fLx;                                    }       //hit X position in LORS, [cm]
  30          Float_t LorsY  (                               )const{return fLy;                                    }       //hit Y position in LORS, [cm]
  31          Int_t   Pid    (                               )const{return fPid;                                   }       //PID
  32          Float_t Q      (                               )const{return fQ;                                     }       //total charge, [QDC]
  33   inline Float_t QdcTot (Float_t e                      );                                                            //calculate total charge of the hit
  34          Int_t   Tid    (                               )const{return fTrack;                                 }       //TID
  35          void    SetQ   (Float_t q                      )     {fQ=q;                                          }       //for debugging...
  36 protected:                                                                     //AliHit has fTrack,fX,fY,fZ
  37   Int_t    fCh;                                                                //Chamber
  38   Int_t    fPid;                                                               //PID
  39   Float_t  fQ;                                                                 //total charge [QDC]
  40   Float_t  fLx;                                                                //hit X position in chamber LORS, [cm]
  41   Float_t  fLy;                                                                //hit Y position in chamber LORS, [cm]
  42   ClassDef(AliHMPIDHit,5)                                                      //HMPID hit class
  43 };//class AliHMPIDhit
  44 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  45 Float_t AliHMPIDHit::QdcTot(Float_t e)
  46 {
  47 // Samples total charge of the hit
  48 // Arguments: e- hit energy [GeV] for mip Eloss for photon Etot
  49 //   Returns: total QDC
  50   Int_t pc,px,py;
  51   AliHMPIDParam::Lors2Pad(fLx,fLy,pc,px,py);
  52   if(py<0) fQ=0;
  53  else {
  54   Float_t y=AliHMPIDParam::LorsY(pc,py);
  55   fLy=((y-fLy)>0)?y-0.2:y+0.2;                                                                       //shift to the nearest anod wire
  56
  57   Float_t  x=(fLx > 66.6)? fLx-66.6:fLx;                                                             //sagita is for PC (0-64) and not for chamber
  58   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
  59
  60   Int_t iNele=Int_t(e/26e-9);  if(iNele<1) iNele=1;                                                  //number of electrons created by hit
  61   fQ=0;
  62   for(Int_t i=1;i<=iNele;i++){
  63     Double_t rnd=gRandom->Rndm(); if(rnd==0) rnd=1e-12;                                              //1e-12 is a protection against 0 from rndm
  64     fQ-=qdcEle*TMath::Log(rnd);
  65   }
  66  }
  67   return fQ;
  68 }
  69 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  70
  71 #endif