Change for a realistic time distribution of the raw signal
[u/mrichter/AliRoot.git] / PHOS / AliPHOSGetter.h
1 #ifndef ALIPHOSGETTER_H
2 #define ALIPHOSGETTER_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 class TParticle ;
22 class TTree ; 
23 class TGraph ;
24 class TF1 ;
25
26 // --- Standard library ---
27
28 // --- AliRoot header files ---
29 #include "AliConfig.h" 
30 #include "AliPHOSLoader.h" 
31 #include "AliPHOSHit.h" 
32 #include "AliPHOSDigit.h"
33 #include "AliPHOSEmcRecPoint.h" 
34 #include "AliPHOSCpvRecPoint.h" 
35 #include "AliPHOSTrackSegment.h" 
36 #include "AliPHOSRecParticle.h" 
37 #include "AliPHOSDigitizer.h"
38 #include "AliPHOSSDigitizer.h"
39 class AliPHOS ;  
40 class AliPHOSGeometry ;
41 class AliPHOSClusterizer ;
42 class AliPHOSTrackSegmentMaker ;  
43 class AliPHOSPID ; 
44 class AliPHOSBeamTestEvent ;
45 class AliESD ; 
46
47 class AliPHOSGetter : public TObject {
48   
49 public:  
50   AliPHOSGetter(){    // ctor: this is a singleton, the ctor should never be called but cint needs it as public
51     Fatal("ctor", "AliPHOSGetter is a singleton default ctor not callable") ;
52   } 
53   AliPHOSGetter(const AliPHOSGetter & obj) : TObject(obj) {
54     // cpy ctor requested by Coding Convention 
55     Fatal("cpy ctor", "not implemented") ;
56   } 
57   
58   AliPHOSGetter & operator = (const AliPHOSGetter & ) {
59     // assignement operator requested by coding convention, but not needed
60     Fatal("operator =", "not implemented") ;
61     return *this ; 
62   }
63   virtual ~AliPHOSGetter() ; 
64   
65   //=========== Instantiators ================
66   static AliPHOSGetter * Instance(const char* headerFile,
67                                   const char* version = AliConfig::GetDefaultEventFolderName(),
68                                   Option_t * openingOption = "READ" ) ; 
69   static AliPHOSGetter * Instance() ; 
70   
71   static void Print() ; 
72   
73   //=========== General information about run ==============
74   Bool_t IsLoaded(TString tree) const { return fLoadingStatus.Contains(tree) ; } 
75   void   SetLoaded(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   //========== PHOSGeometry and PHOS ============= 
84   AliPHOS *         PHOS() const  ;  
85   AliPHOSGeometry * PHOSGeometry() const ; 
86   
87   //========== Methods to read something from file ==========
88   void   Event(Int_t event, const char * opt = "HSDRTP") ;    
89   void   Track(Int_t itrack) ;
90   
91   //-----------------now getter's data--------------------------------------
92   //  AliPHOSCalibrationDB * CalibrationDB(){return  fcdb; }
93   //  void ReadCalibrationDB(const char * name, const char * filename) ;
94   
95   //=========== Primaries ============
96   TClonesArray *    Primaries(void) ;
97   TParticle * Primary(Int_t index) const ;
98   Int_t       NPrimaries()const { return fNPrimaries; }
99   TParticle * Secondary(const TParticle * p, Int_t index=1) const ;  
100   
101   //=========== Hits =================
102   TClonesArray *  Hits(void)  ; 
103   AliPHOSHit *    Hit(Int_t index) { return dynamic_cast<AliPHOSHit*>(Hits()->At(index) );}
104   TTree *         TreeH() const ; 
105   
106   //=========== SDigits ==============
107   TClonesArray *      SDigits() ;  
108   AliPHOSDigit *      SDigit(Int_t index) { return static_cast<AliPHOSDigit *>(SDigits()->At(index)) ;} 
109   TTree *             TreeS() const ; 
110   AliPHOSSDigitizer * SDigitizer() ;  
111   
112   TString             GetSDigitsFileName() const { return PhosLoader()->GetSDigitsFileName() ; }  
113   Int_t               LoadSDigits(Option_t* opt="") const { return PhosLoader()->LoadSDigits(opt) ; }
114   Int_t               LoadSDigitizer(Option_t* opt="") const { return  PhosLoader()->LoadSDigitizer(opt) ; }
115   Int_t               WriteSDigits(Option_t* opt="") const  { return PhosLoader()->WriteSDigits(opt) ; }
116   Int_t               WriteSDigitizer(Option_t* opt="") const {
117     return  PhosLoader()->WriteSDigitizer(opt) ; }
118   
119   //========== Digits ================
120   TClonesArray * Digits() ;
121   AliPHOSDigit * Digit(Int_t index) { return static_cast<AliPHOSDigit *>(Digits()->At(index)) ;} 
122   TTree *        TreeD() const ; 
123   AliPHOSDigitizer * Digitizer() ;
124   TString             GetDigitsFileName() const { return PhosLoader()->GetDigitsFileName() ; }  
125   Int_t               LoadDigits(Option_t* opt="") const { return PhosLoader()->LoadDigits(opt) ; }
126   Int_t               LoadDigitizer(Option_t* opt="") const {
127     return  PhosLoader()->LoadDigitizer(opt) ; }
128   Int_t               WriteDigits(Option_t* opt="") const { return PhosLoader()->WriteDigits(opt) ; }
129   Int_t               WriteDigitizer(Option_t* opt="") const {
130     return  PhosLoader()->WriteDigitizer(opt) ; }
131   
132   //========== RecPoints =============
133   TObjArray *           EmcRecPoints() ;
134   AliPHOSEmcRecPoint *  EmcRecPoint(Int_t index) { return static_cast<AliPHOSEmcRecPoint *>(EmcRecPoints()->At(index)) ;} 
135   TObjArray *           CpvRecPoints() ; 
136   AliPHOSCpvRecPoint *  CpvRecPoint(Int_t index) { return static_cast<AliPHOSCpvRecPoint *>(CpvRecPoints()->At(index)) ;} 
137   TTree *               TreeR() const ;
138   AliPHOSClusterizer * Clusterizer() ;
139   TString               GetRecPointsFileName() const { return PhosLoader()->GetRecPointsFileName() ; } 
140   Int_t                 LoadRecPoints(Option_t* opt="") const { return PhosLoader()->LoadRecPoints(opt) ; }
141   Int_t                 LoadClusterizer(Option_t* opt="") const {
142     return  PhosLoader()->LoadClusterizer(opt) ; }
143   Int_t                 WriteRecPoints(Option_t* opt="") const { return PhosLoader()->WriteRecPoints(opt) ; }
144   Int_t                 WriteClusterizer(Option_t* opt="") const {
145     return  PhosLoader()->WriteClusterizer(opt) ; }
146   
147   //========== TrackSegments   TClonesArray * TrackSegments(const char * name = 0) { 
148   TClonesArray *           TrackSegments() ;
149   AliPHOSTrackSegment *  TrackSegment(Int_t index) { return static_cast<AliPHOSTrackSegment *>(TrackSegments()->At(index)) ;} 
150   TTree *               TreeT() const ;
151   AliPHOSTrackSegmentMaker * TrackSegmentMaker() ;
152   TString               GetTracksFileName() const { return PhosLoader()->GetTracksFileName() ; } 
153   Int_t                 LoadTracks(Option_t* opt="") const { return PhosLoader()->LoadTracks(opt) ; }
154   Int_t                 LoadTrackSegementMaker(Option_t* opt="") const {
155     return  PhosLoader()->LoadTrackSegmentMaker(opt) ; }
156   Int_t                 WriteTracks(Option_t* opt="") const { return PhosLoader()->WriteTracks(opt) ; }
157   Int_t                 WriteTrackSegmentMaker(Option_t* opt="") const {
158     return  PhosLoader()->WriteTracker(opt) ; }
159   
160   //========== RecParticles ===========
161   TClonesArray *         RecParticles() ;
162   AliPHOSRecParticle *   RecParticle(Int_t index) { return static_cast<AliPHOSRecParticle *>(RecParticles()->At(index)) ;} 
163   TTree *               TreeP() const ;
164   AliPHOSPID * PID() ;
165   TString               GetRecParticlesFileName() const { return PhosLoader()->GetRecParticlesFileName() ; } 
166   Int_t                 LoadRecParticles(Option_t* opt="") const { return PhosLoader()->LoadRecParticles(opt) ; }
167   Int_t                 LoadPID(Option_t* opt="") const {
168     return  PhosLoader()->LoadPID(opt) ; }
169   Int_t                 WriteRecParticles(Option_t* opt="") const { return PhosLoader()->WriteRecParticles(opt) ; }
170   Int_t                 WritePID(Option_t* opt="") const {
171     return  PhosLoader()->WritePID(opt) ; }
172
173   //========== Raw ===========
174   Int_t ReadRaw(Int_t event) ; 
175
176   void SetDebug(Int_t level) {fgDebug = level;} // Set debug level 
177   void PostClusterizer(AliPHOSClusterizer * clu) 
178     const{PhosLoader()->PostClusterizer(clu) ; }
179   void PostPID(AliPHOSPID * pid) 
180     const{PhosLoader()->PostPID(pid) ; }
181   void PostTrackSegmentMaker(AliPHOSTrackSegmentMaker * tr) 
182     const{PhosLoader()->PostTrackSegmentMaker(tr) ; }
183   void PostSDigitizer (AliPHOSSDigitizer * sdigitizer) 
184     const {PhosLoader()->PostSDigitizer(sdigitizer);}    
185   void PostDigitizer (AliPHOSDigitizer * digitizer)    
186     const {PhosLoader()->PostDigitizer(dynamic_cast<AliDigitizer *>(digitizer));}
187   
188   TString Version() const  { return PhosLoader()->GetTitle() ; } 
189   AliPHOSLoader * PhosLoader() const { return  fgPhosLoader ; }
190   void Reset() ;
191   
192   AliESD * ESD() const { return fESD ; }
193   
194 private:
195   
196   AliPHOSGetter(const char* headerFile,
197                 const char* version = AliConfig::GetDefaultEventFolderName(),
198                 Option_t * openingOption = "READ") ;
199   
200   Int_t ReadTreeD(void) ;
201   Int_t ReadTreeH(void) ;
202   Int_t ReadTreeR(void) ;
203   Int_t ReadTreeT(void) ;
204   Int_t ReadTreeS(void) ;
205   Int_t ReadTreeP(void) ;
206
207   Int_t ReadTreeE(Int_t event) ;    
208   Bool_t OpenESDFile() ;
209   void ReadPrimaries(void) ;
210
211   void FitRaw(Bool_t lowGainFlag, TGraph * gLowGain, TGraph * gHighGain, TF1* signalF, Int_t & amp, Double_t & time) ; 
212
213 private:
214   
215   AliPHOSBeamTestEvent * fBTE ;           //! Header if BeamTest Event
216   static Int_t          fgDebug ;             //! Debug level
217   
218   TString           fLoadingStatus ;     //! tells which trees are loaded
219   Int_t             fNPrimaries ;        //! # of primaries  
220   TClonesArray *    fPrimaries ;         //! list of lists of primaries
221   TFile *           fESDFile ;           //! ESD file
222   TString           fESDFileName ;       //! ESD File Name
223   AliESD *          fESD ;               //! ESD object
224   TTree *           fESDTree ;           //! ESD Tree
225
226   //  AliPHOSCalibrationDB * fcdb ;       //!
227   
228   static AliPHOSLoader * fgPhosLoader ; // the loader for the NewIO
229   static AliPHOSGetter * fgObjGetter; // pointer to the unique instance of the singleton 
230   
231   enum EDataTypes{kHits,kSDigits,kDigits,kRecPoints,kTracks,kNDataTypes};
232   
233   
234   ClassDef(AliPHOSGetter,1)  // Algorithm class that provides methods to retrieve objects from a list knowing the index 
235     
236     };
237
238 #endif // AliPHOSGETTER_H