RICH/AliRICH.h

   1 #ifndef AliRICH_h
   2 #define AliRICH_h
   3
   4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   5  * See cxx source for full Copyright notice                               */
   6
   7 /* $Id$ */
   8
   9
  10 ////////////////////////////////////////////////
  11 //  Manager and hits classes for set:RICH     //
  12 ////////////////////////////////////////////////
  13
  14 #include <TObjArray.h>
  15 #include <AliDetector.h>
  16 #include <AliRICHConst.h>
  17 #include <AliRICHChamber.h>
  18 static const int kNCH=7;
  19
  20 class TFile;
  21
  22 class AliRICHHit;
  23 class AliRICHSDigit;
  24 class AliRICHRawCluster;
  25 class AliRICHRecHit1D;
  26 class AliRICHRecHit3D;
  27 class AliRICHClusterFinder;
  28 class AliRICHDetect;
  29 class AliRICHChamber;
  30 class AliRICHCerenkov;
  31 class AliSegmentation;
  32 class AliRICHResponse;
  33 class AliRICHGeometry;
  34 class AliRICHMerger;
  35
  36 class AliRICH : public AliDetector
  37 {
  38
  39 enum EDebugBits {kDebugStart=BIT(0),kDebugParam=BIT(1),kDebugHit=BIT(2),kDebugDigit=BIT(3),kDebugReco=BIT(4)}; // Debug flags
  40
  41 public:
  42              AliRICH();                                    // default ctor
  43              AliRICH(const char *name, const char *title); // named ctor
  44              AliRICH(const AliRICH& RICH);                 // copy ctor
  45    virtual  ~AliRICH();                     // dtor
  46 // Pure vituls methods
  47    virtual Int_t  IsVersion() const =0;
  48 // The following staff is defined in AliRICHChamber.cxx
  49    virtual void   AddHit(Int_t track, Int_t *vol, Float_t *hits);
  50    virtual void   AddCerenkov(Int_t track, Int_t *vol, Float_t *cerenkovs);
  51    virtual void   AddSDigit(Int_t *clhits);
  52    virtual void   AddDigits(Int_t id, Int_t *tracks, Int_t *charges, Int_t *digits);
  53    virtual void   AddRawCluster(Int_t id, const AliRICHRawCluster& cluster);
  54    virtual void   AddRecHit1D(Int_t id, Float_t* rechit, Float_t* photons, Int_t* padsx, Int_t* padsy);
  55    virtual void   AddRecHit3D(Int_t id, Float_t* rechit);
  56
  57    virtual void   BuildGeometry();   // TNode ROOT variant for event display
  58    virtual void   CreateGeometry();  // GEANT volumes tree for simulation
  59    virtual void   CreateMaterials(); // GEANT materials definition
  60    virtual Float_t AbsoCH4(Float_t x);
  61    virtual Float_t Fresnel(Float_t ene,Float_t pdoti, Bool_t pola);
  62    virtual void   StepManager();
  63    Int_t          DistancetoPrimitive(Int_t px, Int_t py);
  64    virtual void   MakeBranch(Option_t *opt=" ", const char *file=0);
  65    virtual void   MakeBranchInTreeD(TTree *treeD, const char *file=0);
  66    void           SetTreeAddress();
  67    virtual void   ResetHits();
  68    virtual void   ResetDigits();
  69    virtual void   ResetRawClusters();
  70    virtual void   ResetRecHits1D();
  71    virtual void   ResetRecHits3D();
  72    virtual void   FindClusters(Int_t nev,Int_t lastEntry);
  73 // Converters
  74    virtual void   Hits2SDigits();
  75    virtual Int_t  Hits2SDigits(Float_t xhit,Float_t yhit,Float_t eloss,Int_t id, ResponseType res);
  76    virtual void   SDigits2Digits();
  77    virtual void   SDigits2Digits(Int_t nev, Int_t flag);
  78    virtual void   Digits2Reco();
  79 // Retrieve pad hits for a given Hit
  80     virtual        AliRICHSDigit* FirstPad(AliRICHHit *hit, TClonesArray *clusters);
  81     virtual        AliRICHSDigit* NextPad(TClonesArray *clusters);
  82 // inline methods
  83    TClonesArray  *SDigits()   const{return fSDigits;}
  84    TClonesArray  *Cerenkovs() const{return fCerenkovs;}
  85
  86 // Models for chambers
  87    virtual void     SetGeometryModel(Int_t iChamberN, AliRICHGeometry *pRICHGeo)    {       GetChamber(iChamberN)->SetGeometryModel(pRICHGeo);}
  88    AliRICHGeometry* GetGeometryModel(Int_t iChamberN=0)                        const{return GetChamber(iChamberN)->GetGeometryModel();}
  89
  90    virtual void     SetSegmentationModel(Int_t iChamberN, AliSegmentation *pAliSeg) {       GetChamber(iChamberN)->SetSegmentationModel(pAliSeg);}
  91    AliSegmentation* GetSegmentationModel(Int_t iChamberN=0)                    const{return GetChamber(iChamberN)->GetSegmentationModel();}
  92
  93    virtual void     SetResponseModel(Int_t iChamberN, AliRICHResponse *pRICHRes)    {       GetChamber(iChamberN)->SetResponseModel(pRICHRes);}
  94    AliRICHResponse* GetResponseModel(Int_t iChamberN)                          const{return GetChamber(iChamberN)->GetResponseModel();}
  95
  96    virtual void     SetReconstructionModel(Int_t iChamberN, AliRICHClusterFinder *pRICHReco){GetChamber(iChamberN)->SetReconstructionModel(pRICHReco);}
  97 // Debug staff
  98    void   SetDebugLevel(Int_t level) {fDebugLevel=level;}
  99    Int_t  GetDebugLevel()       const{return fDebugLevel;}
 100
 101    void   SetDebugStart()     {fDebugLevel+=kDebugStart;}        // Controls debug message at the entring point of methods
 102    void ResetDebugStart()     {fDebugLevel-=kDebugStart;}        // Controls debug message at the entring point of methods
 103    Bool_t  IsDebugStart()const{return fDebugLevel&kDebugStart;}  // Controls debug message at the entring point of methods
 104
 105    void   SetDebugParam()     {fDebugLevel+=kDebugParam;}        // Controls debug printout for the parameters
 106    void ResetDebugParam()     {fDebugLevel-=kDebugParam;}        // Controls debug printout for the parameters
 107    Bool_t  IsDebugParam()const{return fDebugLevel&kDebugParam;}  // Controls debug printout for the parameters
 108
 109    void   SetDebugHit()       {fDebugLevel+=kDebugHit;}          // Controls debug printout for hits
 110    void ResetDebugHit()       {fDebugLevel-=kDebugHit;}          // Controls debug printout for hits
 111    Bool_t  IsDebugHit()  const{return fDebugLevel&kDebugHit;}    // Controls debug printout for hits
 112
 113    void   SetDebugDigit()     {fDebugLevel+=kDebugDigit;}        // Controls debug printout for digits
 114    void ResetDebugDigit()     {fDebugLevel-=kDebugDigit;}        // Controls debug printout for digits
 115    Bool_t  IsDebugDigit()const{return fDebugLevel&kDebugDigit;}  // Controls debug printout for digits
 116
 117    void   SetDebugReco()      {fDebugLevel+=kDebugReco;}         // Controls debug printout for reco
 118    void ResetDebugReco()      {fDebugLevel-=kDebugReco;}         // Controls debug printout for reco
 119    Bool_t  IsDebugReco() const{return fDebugLevel&kDebugReco;}   // Controls debug printout for reco
 120
 121
 122     virtual void   SetMerger(AliRICHMerger* thisMerger) {fMerger=thisMerger;}
 123 // Return reference to Chamber #id
 124     virtual        AliRICHChamber& Chamber(Int_t id) {return *((AliRICHChamber *) (*fChambers)[id]);}
 125     AliRICHChamber* GetChamber(Int_t iChamberN)     const{return (AliRICHChamber*) (*fChambers)[iChamberN];}
 126 // Return pointers to digits
 127     TObjArray     *Dchambers() {return fDchambers;}
 128     Int_t         *Ndch() {return fNdch;}
 129     virtual TClonesArray *DigitsAddress(Int_t id) {return ((TClonesArray *) (*fDchambers)[id]);}
 130 // Return pointers to rec. hits
 131     TObjArray     *RecHits1D()                      const{return fRecHits1D;}
 132     Int_t         *Nrechits1D()                          {return fNrechits1D;} // returns array
 133     Int_t         *Nrechits3D()                          {return fNrechits3D;} // returns array
 134    virtual TClonesArray *RecHitsAddress1D(Int_t id) const{return ((TClonesArray *) (*fRecHits1D)[id]);}
 135       TObjArray     *RecHits3D()                    const{return fRecHits3D;}
 136    virtual TClonesArray *RecHitsAddress3D(Int_t id) const{return ((TClonesArray *) (*fRecHits3D)[id]);}
 137    virtual TClonesArray *RawClustAddress(Int_t id)  const{return ((TClonesArray *) (*fRawClusters)[id]);}
 138 // Assignment operator
 139    AliRICH& operator=(const AliRICH& rhs);
 140
 141    virtual void DiagnosticsFE(Int_t evNumber1=0,Int_t evNumber2=0);    // Full events
 142    virtual void DiagnosticsSE(Int_t diaglevel,Int_t evNumber1=0,Int_t evNumber2=0);    // Single events
 143
 144    inline virtual void Print(Option_t *option)const; // Prints debug information
 145
 146 protected:
 147    TObjArray            *fChambers;           // List of RICH chambers aka modules
 148    Int_t                 fNSDigits;           // Number of clusters
 149    Int_t                 fNcerenkovs;         // Number of cerenkovs
 150    TClonesArray         *fSDigits;            // List of clusters
 151    TObjArray            *fDchambers;          // List of digits
 152    TClonesArray         *fCerenkovs;          // List of cerenkovs
 153    Int_t                 fNdch[kNCH];         // Number of digits
 154    TObjArray            *fRawClusters;        // List of raw clusters
 155    TObjArray            *fRecHits1D;          // List of rec. hits
 156    TObjArray            *fRecHits3D;          // List of rec. hits
 157    Int_t                 fNrawch[kNCH];       // Number of raw clusters
 158    Int_t                 fNrechits1D[kNCH];   // Number of rec hits
 159    Int_t                 fNrechits3D[kNCH];   // Number of rec hits
 160    Int_t                 fDebugLevel;         // Source debugging level
 161
 162    Int_t fCkovNumber;                         // Number of Cerenkov photons
 163    Int_t fCkovQuarz;                          // Cerenkovs crossing quartz
 164    Int_t fCkovGap;                            // Cerenkovs crossing gap
 165    Int_t fCkovCsi;                            // Cerenkovs crossing csi
 166    Int_t fLostRfreo;                          // Cerenkovs reflected in freon
 167    Int_t fLostRquar;                          // Cerenkovs reflected in quartz
 168    Int_t fLostAfreo;                          // Cerenkovs absorbed in freon
 169    Int_t fLostAquarz;                         // Cerenkovs absorbed in quartz
 170    Int_t fLostAmeta;                          // Cerenkovs absorbed in methane
 171    Int_t fLostCsi;                            // Cerenkovs below csi quantum efficiency
 172    Int_t fLostWires;                          // Cerenkovs lost in wires
 173    Int_t fFreonProd;                          // Cerenkovs produced in freon
 174    Float_t fMipx;                             // x coord. of MIP
 175    Float_t fMipy;                             // y coord. of MIP
 176    Int_t fFeedbacks;                          // Number of feedback photons
 177    Int_t fLostFresnel;                        // Cerenkovs lost by Fresnel reflection
 178
 179
 180 // Background event for event mixing
 181     Text_t *fFileName;                         //! File with background hits
 182     AliRICHMerger *fMerger;                    //! pointer to merger
 183
 184    ClassDef(AliRICH,1)                        // Main RICH class
 185 };//class AliRICH
 186
 187 inline void AliRICH::Print(Option_t *option)const
 188 {
 189    TObject::Print(option);
 190    if(IsDebugParam()){
 191       GetGeometryModel(0)->Print(option);
 192       GetSegmentationModel(0)->Print(option);
 193       GetResponseModel(0)->Print(option);
 194    }
 195 }//inline void AliRICH::Print(Option_t *option)const
 196
 197 #endif