EMCAL/AliEMCALGetter.h

   1 #ifndef ALIEMCALGETTER_H
   2 #define ALIEMCALGETTER_H
   3 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   4  * See cxx source for full Copyright notice                               */
   5
   6 /* $Id$ */
   7
   8 //_________________________________________________________________________
   9 //  A singleton that returns various objects
  10 //  Should be used on the analysis stage to avoid confusing between different
  11 //  branches of reconstruction tree: e.g. reading RecPoints and TS made from
  12 //  another set of RecPoints.
  13 //
  14 //  The objects are retrived from folders.
  15 //*-- Author: Yves Schutz (SUBATECH) & Dmitri Peressounko (RRC KI & SUBATECH)
  16 //
  17
  18
  19 // --- ROOT system ---
  20 #include "TObject.h"
  21 #include "TClonesArray.h"
  22 class TParticle ;
  23 class TTree ;
  24
  25 // --- Standard library ---
  26
  27 // --- AliRoot header files ---
  28 #include "AliConfig.h"
  29
  30 // #include "AliRun.h"
  31 #include "AliEMCALHit.h"
  32 #include "AliEMCALRecParticle.h"
  33 #include "AliEMCALDigitizer.h"
  34 #include "AliEMCALSDigitizer.h"
  35 #include "AliEMCALLoader.h"
  36 class AliEMCAL ;
  37 class AliEMCALGeometry ;
  38 class AliEMCALClusterizer ;
  39 class AliEMCALRecPoint ;
  40 class AliEMCALTowerRecPoint ;
  41 class AliEMCALTrackSegmentMaker ;
  42 class AliEMCALTrackSegment ;
  43 class AliEMCALPID ;
  44 class AliEMCALBeamTestEvent ;
  45
  46
  47 class AliEMCALGetter : public TObject {
  48
  49  public:
  50   AliEMCALGetter(){    // ctor: this is a singleton, the ctor should never be called but cint needs it as public
  51     Fatal("ctor", "AliEMCALGetter is a singleton default ctor not callable") ;
  52   }
  53   AliEMCALGetter(const AliEMCALGetter & obj):TObject(obj) {
  54     // cpy ctor requested by Coding Convention
  55     Fatal("cpy ctor", "not implemented") ;
  56   }
  57
  58   AliEMCALGetter & operator = (const AliEMCALGetter & ) {
  59     // assignement operator requested by coding convention, but not needed
  60     Fatal("operator =", "not implemented") ;
  61     return *this ;
  62   }
  63   virtual ~AliEMCALGetter() ;
  64
  65   //=========== Instantiators ================
  66   static AliEMCALGetter * Instance(const char* headerFile,
  67                                   const char* version = AliConfig::fgkDefaultEventFolderName,
  68                                   Option_t * openingOption = "READ" ) ;
  69   static AliEMCALGetter * Instance() ;
  70
  71   static void Print() ;
  72
  73 //   //=========== General information about run ==============
  74   Bool_t IsLoaded(const TString tree) const { return fLoadingStatus.Contains(tree) ; }
  75   void   SetLoaded(const TString tree) { fLoadingStatus += tree ; }
  76
  77   Int_t  MaxEvent() const ;
  78   Int_t  EventNumber() const ;
  79   Bool_t VersionExists(TString & opt) const ;
  80   UShort_t EventPattern(void) const ;
  81   Float_t  BeamEnergy(void) const ;
  82
  83 //   //========== EMCALGeometry and EMCAL =============
  84   AliEMCAL *         EMCAL() const  ;
  85   AliEMCALGeometry * EMCALGeometry() const ;
  86
  87 //   //========== Methods to read something from file ==========
  88   void   Event(const Int_t event, const char * opt = "HSDRP") ;
  89   void   Track(const Int_t itrack) ;
  90
  91 //   //-----------------now getter's data--------------------------------------
  92 // //  AliEMCALCalibrationDB * CalibrationDB(){return  fcdb; }
  93 // //  void ReadCalibrationDB(const char * name, const char * filename) ;
  94
  95   //=========== Primaries ============
  96 //   TTree *           TreeK(TString filename="") ;
  97   TClonesArray *    Primaries(void)  ;
  98   TParticle * Primary(Int_t index) const ;
  99   Int_t       NPrimaries()const { return fNPrimaries; }
 100   TParticle * Secondary(const TParticle * p, const Int_t index=1) const ;
 101
 102 //   //=========== Hits =================
 103 //   TTree *               TreeH(TString filename="") ;
 104   TClonesArray *  Hits(void) const  ;
 105   AliEMCALHit *    Hit(const Int_t index) const { return dynamic_cast<AliEMCALHit*>(Hits()->At(index) );}
 106   TTree *         TreeH() const ;
 107
 108   //=========== SDigits ==============
 109   TClonesArray *      SDigits() const ;
 110   AliEMCALDigit *      SDigit(const Int_t index) const { return static_cast<AliEMCALDigit *>(SDigits()->At(index)) ;}
 111   TTree *             TreeS() const ;
 112   AliEMCALSDigitizer * SDigitizer() ;
 113
 114   TString             GetSDigitsFileName() const { return EmcalLoader()->GetSDigitsFileName() ; }
 115   Int_t               LoadSDigits(Option_t* opt="") const { return EmcalLoader()->LoadSDigits(opt) ; }
 116   Int_t               LoadSDigitizer(Option_t* opt="") const { return  EmcalLoader()->LoadSDigitizer(opt) ; }
 117   Int_t               WriteSDigits(Option_t* opt="") const { return EmcalLoader()->WriteSDigits(opt) ; }
 118   Int_t               WriteSDigitizer(Option_t* opt="") const {
 119     return  EmcalLoader()->WriteSDigitizer(opt) ; }
 120
 121   //========== Digits ================
 122   TClonesArray * Digits() const ;
 123   AliEMCALDigit * Digit(const Int_t index) const { return static_cast<AliEMCALDigit *>(Digits()->At(index)) ;}
 124   TTree *        TreeD() const ;
 125   AliEMCALDigitizer * Digitizer() ;
 126   TString             GetDigitsFileName() const { return EmcalLoader()->GetDigitsFileName() ; }
 127   Int_t               LoadDigits(Option_t* opt="") const { return EmcalLoader()->LoadDigits(opt) ; }
 128   Int_t               LoadDigitizer(Option_t* opt="") const {
 129     return  EmcalLoader()->LoadDigitizer(opt) ; }
 130   Int_t               WriteDigits(Option_t* opt="") const { return EmcalLoader()->WriteDigits(opt) ; }
 131   Int_t               WriteDigitizer(Option_t* opt="") const {
 132     return  EmcalLoader()->WriteDigitizer(opt) ; }
 133
 134   //========== RecPoints =============
 135   TObjArray *             PRERecPoints() const ;
 136   AliEMCALRecPoint *      PRERecPoint(const Int_t index) const { return static_cast<AliEMCALRecPoint *>(PRERecPoints()->At(index)) ;}
 137   TObjArray *             ECARecPoints() const ;
 138   AliEMCALTowerRecPoint * ECARecPoint(const Int_t index) const { return static_cast<AliEMCALTowerRecPoint *>(ECARecPoints()->At(index)) ;}
 139   TObjArray *             HCARecPoints() const ;
 140   AliEMCALTowerRecPoint * HCARecPoint(const Int_t index) const { return static_cast<AliEMCALTowerRecPoint *>(HCARecPoints()->At(index)) ;}
 141   TTree *                 TreeR() const ;
 142   AliEMCALClusterizer *   Clusterizer()  ;
 143   TString                 GetRecPointsFileName() const { return EmcalLoader()->GetRecPointsFileName() ; }
 144   Int_t                   LoadRecPoints(Option_t* opt="") const { return EmcalLoader()->LoadRecPoints(opt) ; }
 145   Int_t                   LoadClusterizer(Option_t* opt="") const {
 146     return  EmcalLoader()->LoadClusterizer(opt) ; }
 147   Int_t                   WriteRecPoints(Option_t* opt="") const { return EmcalLoader()->WriteRecPoints(opt) ; }
 148   Int_t                   WriteClusterizer(Option_t* opt="") const {
 149     return  EmcalLoader()->WriteClusterizer(opt) ; }
 150
 151   //========== TrackSegments   TClonesArray * TrackSegments(const char * name = 0) {
 152   TClonesArray *           TrackSegments() const ;
 153   AliEMCALTrackSegment *  TrackSegments(const Int_t index) const { return static_cast<AliEMCALTrackSegment *>(TrackSegments()->At(index)) ;}
 154   TTree *               TreeT() const ;
 155   AliEMCALTrackSegmentMaker * TrackSegmentMaker() ;
 156   TString               GetTracksFileName() const { return EmcalLoader()->GetTracksFileName() ; }
 157   Int_t                 LoadTracks(Option_t* opt="") const { return EmcalLoader()->LoadTracks(opt) ; }
 158   Int_t                 LoadTrackSegementMaker(Option_t* opt="") const {
 159     return  EmcalLoader()->LoadTrackSegmentMaker(opt) ; }
 160   Int_t                 WriteTracks(Option_t* opt="") const { return EmcalLoader()->WriteTracks(opt) ; }
 161   Int_t                 WriteTrackSegmentMaker(Option_t* opt="") const {
 162     return  EmcalLoader()->WriteTracker(opt) ; }
 163   //========== RecParticles ===========
 164
 165   TClonesArray *         RecParticles() const ;
 166   AliEMCALRecParticle *   RecPaticles(const Int_t index) const { return static_cast<AliEMCALRecParticle *>(RecParticles()->At(index)) ;}
 167   TTree *               TreeP() const ;
 168   AliEMCALPID * PID() ;
 169   TString               GetRecParticlesFileName() const { return EmcalLoader()->GetRecParticlesFileName() ; }
 170   Int_t                 LoadRecParticles(Option_t* opt="") const { return EmcalLoader()->LoadRecParticles(opt) ; }
 171   Int_t                 LoadPID(Option_t* opt="") const {
 172     return  EmcalLoader()->LoadPID(opt) ; }
 173   Int_t                 WriteRecParticles(Option_t* opt="") const { return EmcalLoader()->WriteRecParticles(opt) ; }
 174   Int_t                 WritePID(Option_t* opt="") const {
 175     return  EmcalLoader()->WritePID(opt) ; }
 176
 177
 178   void SetDebug(Int_t level) {fgDebug = level;} // Set debug level
 179     void PostClusterizer(AliEMCALClusterizer * clu)
 180     const{EmcalLoader()->PostClusterizer(clu) ; }
 181   void PostPID(AliEMCALPID * pid)
 182     const{EmcalLoader()->PostPID(pid) ; }
 183   void PostTrackSegmentMaker(AliEMCALTrackSegmentMaker * tr)
 184     const{EmcalLoader()->PostTrackSegmentMaker(tr) ; }
 185   void PostSDigitizer (AliEMCALSDigitizer * sdigitizer)
 186     const {EmcalLoader()->PostSDigitizer(sdigitizer);}
 187   void PostDigitizer (AliEMCALDigitizer * digitizer)
 188     const {EmcalLoader()->PostDigitizer(dynamic_cast<AliDigitizer *>(digitizer));}
 189
 190   TString Version() const  { return EmcalLoader()->GetTitle() ; }
 191   AliEMCALLoader * EmcalLoader() const { return  fgEmcalLoader ; }
 192
 193 private:
 194
 195   AliEMCALGetter(const char* headerFile,
 196                 const char* version = AliConfig::fgkDefaultEventFolderName,
 197                 Option_t * openingOption = "READ") ;
 198
 199   Int_t ReadTreeD(void) ;
 200   Int_t ReadTreeH(void) ;
 201   Int_t ReadTreeR(void) ;
 202   Int_t ReadTreeT(void) ;
 203   Int_t ReadTreeS(void) ;
 204   Int_t ReadTreeP(void) ;
 205
 206
 207   void ReadPrimaries(void) ;
 208
 209 private:
 210
 211 //   static TFile * fgFile;           //!
 212
 213 //  AliEMCALBeamTestEvent * fBTE ;           //! Header if BeamTest Event
 214
 215   static Int_t          fgDebug ;             //! Debug level
 216
 217   TString           fLoadingStatus ;     //! tells which trees are loaded
 218   Int_t             fNPrimaries ;        //! # of primaries
 219   TClonesArray *    fPrimaries ;         //! list of lists of primaries
 220
 221 //  AliEMCALCalibrationDB * fcdb ;       //!
 222
 223   static AliEMCALLoader * fgEmcalLoader ; // pointer to EMCAL Loader
 224   static AliEMCALGetter * fgObjGetter; // pointer to the unique instance of the singleton
 225
 226   enum EDataTypes{kHits,kSDigits,kDigits,kRecPoints,kTracks,kNDataTypes};
 227
 228
 229   ClassDef(AliEMCALGetter,2)  // Algorithm class that provides methods to retrieve objects from a list knowing the index
 230
 231 };
 232
 233 #endif // AliEMCALGETTER_H