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 */
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.
14 // The objects are retrived from folders.
15 //*-- Author: Yves Schutz (SUBATECH) & Dmitri Peressounko (RRC KI & SUBATECH)
19 // --- ROOT system ---
21 #include "TClonesArray.h"
27 // --- Standard library ---
29 // --- AliRoot header files ---
30 #include "AliConfig.h"
32 // #include "AliRun.h"
33 #include "AliEMCALHit.h"
34 #include "AliEMCALRecParticle.h"
35 #include "AliEMCALDigitizer.h"
36 #include "AliEMCALSDigitizer.h"
37 #include "AliEMCALLoader.h"
39 class AliEMCALGeometry ;
40 class AliEMCALClusterizer ;
41 class AliEMCALRecPoint ;
42 class AliEMCALTrackSegmentMaker ;
43 class AliEMCALTrackSegment ;
45 class AliEMCALBeamTestEvent ;
48 class AliEMCALGetter : public TObject {
52 // ctor: this is a singleton, the ctor should never be called but cint needs it as public
53 Fatal("ctor", "AliEMCALGetter is a singleton default ctor not callable") ;
55 AliEMCALGetter(const AliEMCALGetter & obj):TObject(obj) {
56 // cpy ctor requested by Coding Convention
57 Fatal("cpy ctor", "not implemented") ;
60 AliEMCALGetter & operator = (const AliEMCALGetter & ) {
61 // assignement operator requested by coding convention, but not needed
62 Fatal("operator =", "not implemented") ;
65 virtual ~AliEMCALGetter() ;
67 //=========== Instantiators ================
68 static AliEMCALGetter * Instance(const char* headerFile,
69 const char* version = AliConfig::GetDefaultEventFolderName(),
70 Option_t * openingOption = "READ" ) ;
71 static AliEMCALGetter * Instance() ;
72 void OpenFile(const char* headerFile,
73 const char* version = AliConfig::GetDefaultEventFolderName(),
74 Option_t * openingOption = "READ" );
78 // //=========== General information about run ==============
79 Bool_t IsLoaded(TString tree) const { return fLoadingStatus.Contains(tree) ; }
80 void SetLoaded(TString tree) { fLoadingStatus += tree ; }
82 Int_t MaxEvent() const ;
83 Int_t EventNumber() const ;
84 Bool_t VersionExists(TString & opt) const ;
85 UShort_t EventPattern(void) const ;
86 Float_t BeamEnergy(void) const ;
88 // //========== EMCALGeometry and EMCAL =============
89 AliEMCAL * EMCAL() const ;
90 AliEMCALGeometry * EMCALGeometry() const ;
92 // //========== Methods to read something from file ==========
93 void Event(Int_t event, const char * opt = "HSDRP") ;
94 void Track(Int_t itrack) ;
96 // //-----------------now getter's data--------------------------------------
97 // // AliEMCALCalibrationDB * CalibrationDB(){return fcdb; }
98 // // void ReadCalibrationDB(const char * name, const char * filename) ;
100 //=========== Primaries ============
101 // TTree * TreeK(TString filename="") ;
102 TParticle * Primary(Int_t index) const ;
103 Int_t NPrimaries() const;
104 TParticle * Secondary(const TParticle * p, Int_t index=1) const ;
106 // //=========== Hits =================
107 // TTree * TreeH(TString filename="") ;
108 TClonesArray * Hits(void) const ;
109 AliEMCALHit * Hit(Int_t index) const { return dynamic_cast<AliEMCALHit*>(Hits()->At(index) );}
110 TTree * TreeH() const ;
112 //=========== SDigits ==============
113 TClonesArray * SDigits() const ;
114 AliEMCALDigit * SDigit(Int_t index) const { return static_cast<AliEMCALDigit *>(SDigits()->At(index)) ;}
115 TTree * TreeS() const ;
116 AliEMCALSDigitizer * SDigitizer() ;
118 TString GetSDigitsFileName() const { return EmcalLoader()->GetSDigitsFileName() ; }
119 Int_t LoadSDigits(Option_t* opt="") const { return EmcalLoader()->LoadSDigits(opt) ; }
120 Int_t LoadSDigitizer(Option_t* opt="") const { return EmcalLoader()->LoadSDigitizer(opt) ; }
121 Int_t WriteSDigits(Option_t* opt="") const { return EmcalLoader()->WriteSDigits(opt) ; }
122 Int_t WriteSDigitizer(Option_t* opt="") const {
123 return EmcalLoader()->WriteSDigitizer(opt) ; }
125 //========== Digits ================
126 TClonesArray * Digits() const ;
127 AliEMCALDigit * Digit(Int_t index) const { return static_cast<AliEMCALDigit *>(Digits()->At(index)) ;}
128 TTree * TreeD() const ;
129 AliEMCALDigitizer * Digitizer() ;
130 TString GetDigitsFileName() const { return EmcalLoader()->GetDigitsFileName() ; }
131 Int_t LoadDigits(Option_t* opt="") const { return EmcalLoader()->LoadDigits(opt) ; }
132 Int_t LoadDigitizer(Option_t* opt="") const {
133 return EmcalLoader()->LoadDigitizer(opt) ; }
134 Int_t WriteDigits(Option_t* opt="") const { return EmcalLoader()->WriteDigits(opt) ; }
135 Int_t WriteDigitizer(Option_t* opt="") const {
136 return EmcalLoader()->WriteDigitizer(opt) ; }
138 //========== RecPoints =============
139 TObjArray * ECARecPoints() const;
140 AliEMCALRecPoint * ECARecPoint(Int_t index) const{ return static_cast<AliEMCALRecPoint *>(ECARecPoints()->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) ; }
151 //========== TrackSegments TClonesArray * TrackSegments(const char * name = 0) {
152 TClonesArray * TrackSegments() const ;
153 AliEMCALTrackSegment * TrackSegments(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 ===========
165 TClonesArray * RecParticles() const ;
166 AliEMCALRecParticle * RecParticle(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) ; }
177 //========== Raw ===========
178 Int_t ReadRaw(Int_t event) ;
180 void SetDebug(Int_t level) {fgDebug = level;} // Set debug level
181 void PostClusterizer(AliEMCALClusterizer * clu)
182 const{EmcalLoader()->PostClusterizer(clu) ; }
183 void PostPID(AliEMCALPID * pid)
184 const{EmcalLoader()->PostPID(pid) ; }
185 void PostTrackSegmentMaker(AliEMCALTrackSegmentMaker * tr)
186 const{EmcalLoader()->PostTrackSegmentMaker(tr) ; }
187 void PostSDigitizer (AliEMCALSDigitizer * sdigitizer)
188 const {EmcalLoader()->PostSDigitizer(sdigitizer);}
189 void PostDigitizer (AliEMCALDigitizer * digitizer)
190 const {EmcalLoader()->PostDigitizer(dynamic_cast<AliDigitizer *>(digitizer));}
192 TString Version() const { return EmcalLoader()->GetTitle() ; }
193 AliEMCALLoader * EmcalLoader() const { return fgEmcalLoader; }
198 AliEMCALGetter(const char* headerFile,
199 const char* version = AliConfig::GetDefaultEventFolderName(),
200 Option_t * openingOption = "READ") ;
203 Int_t ReadTreeD(void) ;
204 Int_t ReadTreeH(void) ;
205 Int_t ReadTreeR(void) ;
206 Int_t ReadTreeT(void) ;
207 Int_t ReadTreeS(void) ;
208 Int_t ReadTreeP(void) ;
211 void ReadPrimaries(void) ;
213 void FitRaw(Bool_t lowGainFlag, TGraph * gLowGain, TGraph * gHighGain, TF1* signalF, Int_t & amp, Double_t & time) ;
217 // static TFile * fgFile; //!
219 // AliEMCALBeamTestEvent * fBTE ; //! Header if BeamTest Event
221 static Int_t fgDebug ; //! Debug level
223 TString fLoadingStatus ; //! tells which trees are loaded
224 static TString fVersion; //! stores the current folder name
226 // AliEMCALCalibrationDB * fcdb ; //!
228 static AliEMCALLoader * fgEmcalLoader ; // pointer to EMCAL Loader
229 static AliEMCALGetter * fgObjGetter; // pointer to the unique instance of the singleton
231 enum EDataTypes{kHits,kSDigits,kDigits,kRecPoints,kTracks,kNDataTypes};
234 ClassDef(AliEMCALGetter,5) // Algorithm class that provides methods to retrieve objects from a list knowing the index
238 #endif // AliEMCALGETTER_H