db11014f6216af290e7d4a32e41438057e65af17
[u/mrichter/AliRoot.git] / 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 "TClonesArray.h"
21 #include "TFolder.h"  
22 #include "TTree.h"
23 #include "TFile.h"
24 class TString ;
25 class TParticle ;
26 class TTask ;
27
28 // --- Standard library ---
29 #include <stdlib.h>
30 #include <iostream.h>
31
32 // --- AliRoot header files ---
33
34 #include "AliRun.h"
35 #include "AliEMCALv1.h" 
36 #include "AliEMCALHit.h" 
37 #include "AliEMCALDigit.h" 
38 #include "AliEMCALDigitizer.h" 
39 #include "AliEMCALSDigitizer.h"
40 #include "AliEMCALTowerRecPoint.h"
41 class AliEMCALGeometry ;
42 class AliEMCALClusterizerv1 ;
43 //class AliEMCALTrackSegment ;
44 //class AliEMCALTrackSegmentMaker ;
45 //class AliEMCALRecParticle ;
46 //class AliEMCALPID ;
47
48 class AliEMCALGetter : public TObject {
49   
50  public:
51   
52   AliEMCALGetter(){ 
53     // ctor: this is a singleton, the ctor should never be called but cint needs it as public
54     cerr << "ERROR: AliPHOGetter is a singleton default ctor not callable" << endl ;
55     abort() ; 
56   } 
57   AliEMCALGetter(const AliEMCALGetter & obj) {
58     // cpy ctor requested by Coding Convention 
59     // but not yet needed
60     abort() ; 
61   } 
62   
63   virtual ~AliEMCALGetter() ; 
64   
65   void CloseFile() ;  
66   const TFolder * Folder(const TString what) const ;
67   const Bool_t HasFailed(void) const {return fFailed ;} 
68   Bool_t PostHits(void ) const ;  
69   Bool_t PostSDigits(      const char * name,  const char * file = 0) const ;  
70   Bool_t PostDigits(       const char * name ) const ;  
71   Bool_t PostRecPoints(    const char * name ) const ;  
72   //Bool_t PostTrackSegments(const char * name) const ;  
73   //Bool_t PostRecParticles( const char * name) const ;  
74
75   Bool_t PostClusterizer( const char * name) const ;  
76   Bool_t PostClusterizer(AliEMCALClusterizerv1 * clu) const ;  
77   Bool_t PostSDigitizer (AliEMCALSDigitizer * sdigitizer) const ;  
78   Bool_t PostSDigitizer ( const char * name, const char * file ) const ;  
79   Bool_t PostDigitizer (AliEMCALDigitizer * digitizer) const ;  
80   Bool_t PostDigitizer  ( const char * name) const ;  
81   //Bool_t PostTrackSegmentMaker(AliEMCALTrackSegmentMaker * tsm) const ;  
82   //Bool_t PostTrackSegmentMaker(const char * name ) const ;  
83   //Bool_t PostPID  (AliEMCALPID * pid) const ;  
84   //Bool_t PostPID  (const char * name ) const ;  
85   //Bool_t PostQA   (void) const ;
86   
87
88   void   Event(const Int_t event, const char * opt = "HSDR") ;    
89   void   Track(Int_t itrack) ;
90
91   //Method to be used when digitizing under AliRunDigitizer, who opens all files etc.
92   void   ReadTreeS(TTree * treeS,Int_t input) ;
93   
94   Int_t  EventNumber()       { return (Int_t) gAlice->GetEvNumber() ; }
95   Int_t  MaxEvent()          { return (Int_t) gAlice->TreeE()->GetEntries() ; }
96   static AliEMCALGetter * GetInstance(const char* headerFile,
97                                      const char* branchTitle = "Default", const Option_t * rw="" ) ; 
98   static AliEMCALGetter *   GetInstance() ; 
99
100   const AliEMCALv1 *         EMCAL()  ;  
101   AliEMCALGeometry * EMCALGeometry() ; 
102    // Alarms
103   //TFolder * Alarms() const { return (TFolder*)(ReturnO("Alarms", 0)) ; }
104   //TObjArray *  Alarms(const char * name ) const { return (TObjArray*)(ReturnO("Alarms", name)) ; }
105
106   // QA Tasks
107   //TTask * QATasks(const char * name = 0) const { return (TTask*)(ReturnT("QATasks", name)) ; }
108
109   // Primaries
110   TClonesArray *  Primaries(void) const { return (TClonesArray*)(ReturnO("Primaries")) ; }
111   
112   
113   // Hits
114   const TClonesArray *  Hits(void) { return static_cast<const TClonesArray*>(ReturnO("Hits")) ; }
115   const AliEMCALHit * Hit(Int_t index)  { return static_cast<const AliEMCALHit*>(Hits()->At(index) );}
116   
117   // SDigits
118   TClonesArray *  SDigits(const char * name = 0, const char * file=0) { 
119     return static_cast<TClonesArray*>(ReturnO("SDigits", name, file)) ; 
120   }
121   const AliEMCALDigit *  SDigit(Int_t index) { return static_cast<const AliEMCALDigit*>(SDigits()->At(index)) ;}
122   
123   AliEMCALSDigitizer *  SDigitizer(const char * name =0) const { 
124     return ((AliEMCALSDigitizer*)(ReturnT("SDigitizer", name))) ; 
125   }
126  
127   // Digits
128    TClonesArray *  Digits(const char * name = 0)const  { 
129     return static_cast<TClonesArray*>(ReturnO("Digits", name)) ; 
130   }
131   const AliEMCALDigit *  Digit(Int_t index) { return static_cast<const AliEMCALDigit *>(Digits()->At(index)) ;}
132   AliEMCALDigitizer *  Digitizer(const char * name =0) const { 
133     return (AliEMCALDigitizer*)(ReturnT("Digitizer", name)) ; 
134   }
135
136   //  RecPoints
137   TObjArray * TowerRecPoints(const char * name = 0) const { 
138     return (TObjArray*)(ReturnO("TowerRecPoints", name)) ; }
139   TObjArray * PreShowerRecPoints(const char * name = 0) const { 
140     return (TObjArray*)(ReturnO("PreShoRecPoints", name)) ; }
141   const AliEMCALTowerRecPoint *  TowerRecPoint(Int_t index) { 
142     return static_cast<const AliEMCALTowerRecPoint *>(TowerRecPoints()->At(index)) ;}
143   const AliEMCALTowerRecPoint *  PreShowerRecPoint(Int_t index) { 
144     return static_cast<const AliEMCALTowerRecPoint *>(PreShowerRecPoints()->At(index)) ;}
145
146     AliEMCALClusterizerv1 * Clusterizer (const char * name =0) const { 
147     return (AliEMCALClusterizerv1*)(ReturnT("Clusterizer", name)) ; }
148
149   // TrackSegments
150   //TClonesArray * TrackSegments(const char * name = 0) const 
151    //                { return (TClonesArray*)(ReturnO("TrackSegments", name)) ; }
152   //AliEMCALTrackSegmentMaker * TrackSegmentMaker (const char * name =0) const 
153     //               { return (AliEMCALTrackSegmentMaker*)(ReturnT("TrackSegmentMaker", name)) ; }
154
155   // RecParticles
156   //TClonesArray * RecParticles(const char * name = 0) const  
157     //               { return (TClonesArray*)(ReturnO("RecParticles", name)) ; }
158     //AliEMCALPID * PID(const char * name =0) const 
159       //             { return (AliEMCALPID*)(ReturnT("PID", name)) ; }
160
161   // Primaries
162   const TParticle *           Primary(Int_t index) const ;
163   const Int_t                 NPrimaries()const { return fNPrimaries; }
164
165   void RemoveTask(TString opt, TString name) const ;
166   void RemoveObjects(TString opt, TString name) const ; 
167   void RemoveSDigits() const ; 
168   void SetDebug(Int_t level) {fDebug = level;} // Set debug level
169
170   AliEMCALGetter & operator = (const AliEMCALGetter & ) {
171     // assignement operator requested by coding convention, but not needed
172     abort() ;
173     return *this ; 
174   }
175   
176   TFolder * SDigitsFolder() { return dynamic_cast<TFolder*>(fSDigitsFolder->FindObject("EMCAL")) ; }
177
178  private:
179
180   AliEMCALGetter(const char* headerFile, const char* branchTitle ="Default", const Option_t * rw ="") ; 
181   void CreateWhiteBoard() const ; 
182   TObject * ReturnO(TString what, TString name=0, TString file=0) const ; 
183   const TTask * ReturnT(TString what,TString name=0) const ; 
184   void DefineBranchTitles(char* branch, char* branchTitle) ;
185   void ReadTreeD() ;
186   void ReadTreeH() ;
187   void ReadTreeR() ;
188   void ReadTreeS(Int_t event) ;
189   //void ReadTreeQA() ;
190   void ReadPrimaries() ;
191
192   TObject ** HitsRef(void) const ;
193   TObject ** SDigitsRef(const char * name, const char * file = 0 ) const;
194   TObject ** DigitsRef (const char * name)   const ;
195   TObject ** TowerRecPointsRef (const char * name) const ;
196   TObject ** PreShoRecPointsRef (const char * name) const ;
197   //TObject ** TrackSegmentsRef(const char * name)   const ;
198   //TObject ** RecParticlesRef (const char * name)   const ;
199   //TObject ** AlarmsRef (void)   const ;
200
201   TObject ** SDigitizerRef (const char * name) const ; 
202   TObject ** DigitizerRef  (const char * name) const ; 
203   TObject ** ClusterizerRef(const char * name) const ; 
204   //TObject ** TSMakerRef    (const char * name) const ; 
205   //TObject ** PIDRef        (const char * name) const ; 
206
207  private:
208
209   static TFile *        fFile ;              //!
210   TString        fHeaderFile ;        //! File in which gAlice lives
211   TString        fBranchTitle ;       //!
212   //TString        fTrackSegmentsTitle ;//! 
213   TString        fRecPointsTitle ;    //!
214   //TString        fRecParticlesTitle ; //!
215   TString        fDigitsTitle ;       //!
216   TString        fSDigitsTitle ;      //!
217
218   Bool_t         fFailed ;            //! true if file is not opened and/or galice not found
219   Int_t          fDebug ;             // Debug level
220
221   Int_t          fNPrimaries ;        //! # of primaries
222   
223   TObjArray *    fPrimaries ;         //! list of lists of primaries-for the case of mixing
224
225   TFolder *      fModuleFolder ;      //!Folder that contains the modules 
226   TFolder *      fHitsFolder ;        //!Folder that contains the Hits 
227   TFolder *      fSDigitsFolder ;     //!Folder that contains the SDigits 
228   TFolder *      fDigitsFolder ;      //!Folder that contains the Digits 
229   TFolder *      fRecoFolder ;        //!Folder that contains the reconstructed objects (RecPoints, TrackSegments, RecParticles) 
230   //TFolder *      fQAFolder ;          //!Folder that contains the QA objects  
231   TFolder *      fTasksFolder ;       //!Folder that contains the Tasks (sdigitizer, digitizer, reconstructioner)
232  
233   static AliEMCALGetter * fgObjGetter; // pointer to the unique instance of the singleton 
234
235   ClassDef(AliEMCALGetter,1)  // Algorithm class that provides methods to retrieve objects from a list knowing the index 
236
237 };
238
239 #endif // AliEMCALGETTER_H