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=" ", const char *file=0);
  62     virtual void   MakeBranchInTreeD(TTree *treeD, const char *file=0);
  63     void           SetTreeAddress();
  64     virtual void   ResetHits();
  65     virtual void   ResetDigits();
  66     virtual void   ResetRawClusters();
  67     virtual void   ResetRecHits1D();
  68     virtual void   ResetRecHits3D();
  69     virtual void   FindClusters(Int_t nev,Int_t lastEntry);
  70     virtual void   Hits2SDigits();
  71     virtual void   SDigits2Digits();
  72     virtual void   SDigits2Digits(Int_t nev, Int_t flag);
  73     virtual void   Digits2Reco();
  74 //
  75 // Configuration Methods (per station id)
  76 //
  77 // Set Chamber Segmentation Parameters
  78 // id refers to the station and isec to the cathode plane
  79 // Set Segmentation and Response Model
  80     virtual void   SetGeometryModel(Int_t id, AliRICHGeometry *geometry);
  81     virtual void   SetSegmentationModel(Int_t id, AliSegmentation *segmentation);
  82     virtual void   SetResponseModel(Int_t id, AliRICHResponse *response);
  83 // Set Reconstruction Model
  84     virtual void   SetReconstructionModel(Int_t id, AliRICHClusterFinder *reconstruction);
  85 // Set source debugging level
  86     void           SetDebugLevel(Int_t level) {fDebugLevel=level;}
  87 // Set Merger
  88     virtual void   SetMerger(AliRICHMerger* thisMerger) {fMerger=thisMerger;}
  89 // Get source debugging level
  90     Int_t          GetDebugLevel() {return fDebugLevel;}
  91 // Response Simulation
  92     virtual Int_t  Hits2SDigits(Float_t xhit,Float_t yhit,Float_t eloss,Int_t id, ResponseType res);
  93 // Return reference to Chamber #id
  94     virtual        AliRICHChamber& Chamber(Int_t id) {return *((AliRICHChamber *) (*fChambers)[id]);}
  95 // Retrieve pad hits for a given Hit
  96     virtual        AliRICHSDigit* FirstPad(AliRICHHit *hit, TClonesArray *clusters);
  97     virtual        AliRICHSDigit* NextPad(TClonesArray *clusters);
  98 // Return pointers to digits
  99     TObjArray     *Dchambers() {return fDchambers;}
 100     Int_t         *Ndch() {return fNdch;}
 101     virtual TClonesArray *DigitsAddress(Int_t id) {return ((TClonesArray *) (*fDchambers)[id]);}
 102 // Return pointers to rec. hits
 103     TObjArray     *RecHits1D() {return fRecHits1D;}
 104     Int_t         *Nrechits1D() {return fNrechits1D;}
 105     virtual TClonesArray *RecHitsAddress1D(Int_t id) {return ((TClonesArray *) (*fRecHits1D)[id]);}
 106     TObjArray     *RecHits3D() {return fRecHits3D;}
 107     Int_t         *Nrechits3D() {return fNrechits3D;}
 108     virtual TClonesArray *RecHitsAddress3D(Int_t id) {return ((TClonesArray *) (*fRecHits3D)[id]);}
 109
 110 // Return pointers to reconstructed clusters
 111     virtual TClonesArray *RawClustAddress(Int_t id) {return ((TClonesArray *) (*fRawClusters)[id]);}
 112 // Assignment operator
 113     AliRICH& operator=(const AliRICH& rhs);
 114
 115 // Analysis routines
 116     // Full events
 117     virtual void DiagnosticsFE(Int_t evNumber1=0,Int_t evNumber2=0);
 118     // Single events
 119     virtual void DiagnosticsSE(Int_t diaglevel,Int_t evNumber1=0,Int_t evNumber2=0);
 120
 121  protected:
 122     TObjArray            *fChambers;           // List of Tracking Chambers
 123     Int_t                 fNSDigits;           // Number of clusters
 124     Int_t                 fNcerenkovs;         // Number of cerenkovs
 125     TClonesArray         *fSDigits;            // List of clusters
 126     TObjArray            *fDchambers;          // List of digits
 127     TClonesArray         *fCerenkovs;          // List of cerenkovs
 128     Int_t                 fNdch[kNCH];         // Number of digits
 129     TObjArray            *fRawClusters;        // List of raw clusters
 130     TObjArray            *fRecHits1D;          // List of rec. hits
 131     TObjArray            *fRecHits3D;          // List of rec. hits
 132     Int_t                 fNrawch[kNCH];       // Number of raw clusters
 133     Int_t                 fNrechits1D[kNCH];   // Number of rec hits
 134     Int_t                 fNrechits3D[kNCH];   // Number of rec hits
 135     Int_t                 fDebugLevel;         // Source debugging level
 136
 137     Int_t fCkovNumber;                         // Number of Cerenkov photons
 138     Int_t fCkovQuarz;                          // Cerenkovs crossing quartz
 139     Int_t fCkovGap;                            // Cerenkovs crossing gap
 140     Int_t fCkovCsi;                            // Cerenkovs crossing csi
 141     Int_t fLostRfreo;                          // Cerenkovs reflected in freon
 142     Int_t fLostRquar;                          // Cerenkovs reflected in quartz
 143     Int_t fLostAfreo;                          // Cerenkovs absorbed in freon
 144     Int_t fLostAquarz;                         // Cerenkovs absorbed in quartz
 145     Int_t fLostAmeta;                          // Cerenkovs absorbed in methane
 146     Int_t fLostCsi;                            // Cerenkovs below csi quantum efficiency
 147     Int_t fLostWires;                          // Cerenkovs lost in wires
 148     Int_t fFreonProd;                          // Cerenkovs produced in freon
 149     Float_t fMipx;                             // x coord. of MIP
 150     Float_t fMipy;                             // y coord. of MIP
 151     Int_t fFeedbacks;                          // Number of feedback photons
 152     Int_t fLostFresnel;                        // Cerenkovs lost by Fresnel reflection
 153
 154
 155 // Background eent for event mixing
 156     Text_t *fFileName;                         // ! File with background hits
 157     AliRICHMerger *fMerger;                    // ! pointer to merger
 158
 159     ClassDef(AliRICH,1)                        //Hits manager for set:RICH
 160 };
 161 #endif
 162