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 #include "AliDetector.h"
  14 #include "AliRICHConst.h"
  15 #include "AliRICHChamber.h"
  16 static const int kNCH=7;
  17
  18 class TFile;
  19
  20 class AliRICHHit;
  21 class AliRICHSDigit;
  22 class AliRICHRawCluster;
  23 class AliRICHRecHit1D;
  24 class AliRICHRecHit3D;
  25 class AliRICHClusterFinder;
  26 class AliRICHDetect;
  27 class AliRICHChamber;
  28 class AliRICHCerenkov;
  29 class AliSegmentation;
  30 class AliRICHResponse;
  31 class AliRICHEllipse;
  32 class AliRICHGeometry;
  33 class AliRICHMerger;
  34
  35 class AliRICH : public  AliDetector {
  36  public:
  37     AliRICH();
  38     AliRICH(const char *name, const char *title);
  39     AliRICH(const AliRICH& RICH);
  40     virtual       ~AliRICH();
  41     virtual void   AddHit(Int_t track, Int_t *vol, Float_t *hits);
  42     virtual void   AddCerenkov(Int_t track, Int_t *vol, Float_t *cerenkovs);
  43     virtual void   AddSDigit(Int_t *clhits);
  44     virtual void   AddDigits(Int_t id, Int_t *tracks, Int_t *charges, Int_t *digits);
  45     virtual void   AddRawCluster(Int_t id, const AliRICHRawCluster& cluster);
  46     virtual void   AddRecHit1D(Int_t id, Float_t* rechit, Float_t* photons, Int_t* padsx, Int_t* padsy);
  47     virtual void   AddRecHit3D(Int_t id, Float_t* rechit);
  48
  49
  50     virtual void   BuildGeometry();
  51     virtual void   CreateGeometry();
  52     virtual void   CreateMaterials();
  53     virtual Float_t AbsoCH4(Float_t x);
  54     virtual Float_t Fresnel(Float_t ene,Float_t pdoti, Bool_t pola);
  55     virtual void   StepManager();
  56     Int_t          DistancetoPrimitive(Int_t px, Int_t py);
  57     virtual Int_t  IsVersion() const =0;
  58 //
  59     TClonesArray  *SDigits() {return fSDigits;}
  60     TClonesArray  *Cerenkovs() {return fCerenkovs;}
  61     virtual void   MakeBranch(Option_t *opt=" ", char *file=0);
  62     void           SetTreeAddress();
  63     virtual void   ResetHits();
  64     virtual void   ResetDigits();
  65     virtual void   ResetRawClusters();
  66     virtual void   ResetRecHits1D();
  67     virtual void   ResetRecHits3D();
  68     virtual void   FindClusters(Int_t nev,Int_t lastEntry);
  69     virtual void   Digitise(Int_t nev,Int_t flag,Option_t *opt=" ",Text_t *name=" ");
  70     virtual void   SDigits2Digits();
  71 //
  72 // Configuration Methods (per station id)
  73 //
  74 // Set Chamber Segmentation Parameters
  75 // id refers to the station and isec to the cathode plane
  76 // Set Segmentation and Response Model
  77     virtual void   SetGeometryModel(Int_t id, AliRICHGeometry *geometry);
  78     virtual void   SetSegmentationModel(Int_t id, AliSegmentation *segmentation);
  79     virtual void   SetResponseModel(Int_t id, AliRICHResponse *response);
  80 // Set Reconstruction Model
  81     virtual void   SetReconstructionModel(Int_t id, AliRICHClusterFinder *reconstruction);
  82 // Set source debugging level
  83     void SetDebugLevel(Int_t level) {fDebugLevel=level;}
  84 // Set Merger
  85     virtual void   SetMerger(AliRICHMerger* thisMerger) {fMerger=thisMerger;}
  86 // Get source debugging level
  87     Int_t GetDebugLevel() {return fDebugLevel;}
  88 // Response Simulation
  89     virtual Int_t   Hits2SDigits(Float_t xhit,Float_t yhit,Float_t eloss,Int_t id, ResponseType res);
  90 // Return reference to Chamber #id
  91     virtual AliRICHChamber& Chamber(Int_t id) {return *((AliRICHChamber *) (*fChambers)[id]);}
  92 // Retrieve pad hits for a given Hit
  93     virtual AliRICHSDigit* FirstPad(AliRICHHit *hit, TClonesArray *clusters);
  94     virtual AliRICHSDigit* NextPad(TClonesArray *clusters);
  95 // Return pointers to digits
  96     TObjArray            *Dchambers() {return fDchambers;}
  97     Int_t                *Ndch() {return fNdch;}
  98     virtual TClonesArray *DigitsAddress(Int_t id) {return ((TClonesArray *) (*fDchambers)[id]);}
  99 // Return pointers to rec. hits
 100     TObjArray            *RecHits1D() {return fRecHits1D;}
 101     Int_t                *Nrechits1D() {return fNrechits1D;}
 102     virtual TClonesArray *RecHitsAddress1D(Int_t id) {return ((TClonesArray *) (*fRecHits1D)[id]);}
 103     TObjArray            *RecHits3D() {return fRecHits3D;}
 104     Int_t                *Nrechits3D() {return fNrechits3D;}
 105     virtual TClonesArray *RecHitsAddress3D(Int_t id) {return ((TClonesArray *) (*fRecHits3D)[id]);}
 106
 107 // Return pointers to reconstructed clusters
 108     virtual TClonesArray *RawClustAddress(Int_t id) {return ((TClonesArray *) (*fRawClusters)[id]);}
 109 // Assignment operator
 110     AliRICH& operator=(const AliRICH& rhs);
 111
 112
 113  protected:
 114     TObjArray            *fChambers;           // List of Tracking Chambers
 115     Int_t                 fNSDigits;           // Number of clusters
 116     Int_t                 fNcerenkovs;         // Number of cerenkovs
 117     TClonesArray         *fSDigits;            // List of clusters
 118     TObjArray            *fDchambers;          // List of digits
 119     TClonesArray         *fCerenkovs;          // List of cerenkovs
 120     Int_t                 fNdch[kNCH];         // Number of digits
 121     TObjArray            *fRawClusters;        // List of raw clusters
 122     TObjArray            *fRecHits1D;          // List of rec. hits
 123     TObjArray            *fRecHits3D;          // List of rec. hits
 124     Int_t                 fNrawch[kNCH];       // Number of raw clusters
 125     Int_t                 fNrechits1D[kNCH];   // Number of rec hits
 126     Int_t                 fNrechits3D[kNCH];   // Number of rec hits
 127     Int_t                 fDebugLevel;         // Source debugging level
 128
 129     Int_t fCkovNumber;                   // Number of Cerenkov photons
 130     Int_t fCkovQuarz;                    // Cerenkovs crossing quartz
 131     Int_t fCkovGap;                      // Cerenkovs crossing gap
 132     Int_t fCkovCsi;                      // Cerenkovs crossing csi
 133     Int_t fLostRfreo;                    // Cerenkovs reflected in freon
 134     Int_t fLostRquar;                    // Cerenkovs reflected in quartz
 135     Int_t fLostAfreo;                    // Cerenkovs absorbed in freon
 136     Int_t fLostAquarz;                   // Cerenkovs absorbed in quartz
 137     Int_t fLostAmeta;                    // Cerenkovs absorbed in methane
 138     Int_t fLostCsi;                      // Cerenkovs below csi quantum efficiency
 139     Int_t fLostWires;                    // Cerenkovs lost in wires
 140     Int_t fFreonProd;                    // Cerenkovs produced in freon
 141     Float_t fMipx;                       // x coord. of MIP
 142     Float_t fMipy;                       // y coord. of MIP
 143     Int_t fFeedbacks;                    // Number of feedback photons
 144     Int_t fLostFresnel;                  // Cerenkovs lost by Fresnel reflection
 145
 146
 147 // Background eent for event mixing
 148     Text_t *fFileName;           // ! File with background hits
 149     AliRICHMerger *fMerger;   // ! pointer to merger
 150
 151     ClassDef(AliRICH,1)  //Hits manager for set:RICH
 152 };
 153 #endif
 154
 155
 156
 157
 158
 159
 160
 161
 162
 163
 164
 165
 166
 167
 168