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