]>
Commit | Line | Data |
---|---|---|
2f614988 | 1 | #ifndef AliRICH_h |
2 | #define AliRICH_h | |
2e5f0f7b | 3 | |
3da30618 | 4 | /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * |
5 | * See cxx source for full Copyright notice */ | |
6 | ||
7 | /* $Id$ */ | |
8 | ||
2e5f0f7b | 9 | |
c11ee6f0 | 10 | #include <TObjArray.h> |
853634d3 | 11 | #include <TClonesArray.h> |
c11ee6f0 | 12 | #include <AliDetector.h> |
853634d3 | 13 | #include <AliHit.h> |
14 | #include <AliDigit.h> | |
dfb4e77d | 15 | #include "AliRICHConst.h" |
16 | #include "AliRICHChamber.h" | |
853634d3 | 17 | |
237c933d | 18 | static const int kNCH=7; |
fe4da5cc | 19 | |
2ab0c725 | 20 | class TFile; |
fe4da5cc | 21 | |
2e5f0f7b | 22 | class AliRICHRawCluster; |
a4622d0f | 23 | class AliRICHRecHit1D; |
24 | class AliRICHRecHit3D; | |
2e5f0f7b | 25 | class AliRICHClusterFinder; |
26 | class AliRICHDetect; | |
27 | class AliRICHChamber; | |
28 | class AliRICHCerenkov; | |
a2f7eaf6 | 29 | class AliSegmentation; |
2e5f0f7b | 30 | class AliRICHResponse; |
237c933d | 31 | class AliRICHGeometry; |
d8a72780 | 32 | class AliRICHMerger; |
ddae0931 | 33 | |
853634d3 | 34 | |
35 | ||
36 | class AliRICHhit : public AliHit | |
37 | { | |
38 | public: | |
39 | inline AliRICHhit(); | |
40 | inline AliRICHhit(Int_t fIshunt, Int_t track, Int_t *vol, Float_t *hits); | |
41 | virtual ~AliRICHhit() {;} | |
42 | ||
43 | Int_t Chamber() {return fChamber;} | |
44 | Float_t Particle() {return fParticle;} | |
45 | Float_t Theta() {return fTheta;} | |
46 | Float_t Phi() {return fPhi;} | |
47 | Float_t Tlength() {return fTlength;} | |
48 | Float_t Eloss() {return fEloss;} | |
49 | Float_t Loss() {return fLoss;} | |
50 | Float_t PHfirst() {return fPHfirst;} | |
51 | Float_t PHlast() {return fPHlast;} | |
52 | Float_t MomX() {return fMomX;} | |
53 | Float_t MomY() {return fMomY;} | |
54 | Float_t MomZ() {return fMomZ;} | |
55 | Float_t CerenkovAngle() {return fMomX;} | |
56 | Float_t MomFreoX() {return fMomX;} | |
57 | Float_t MomFreoY() {return fMomY;} | |
58 | Float_t MomFreoZ() {return fMomZ;} | |
59 | protected: | |
60 | Int_t fChamber; // Chamber number | |
61 | Float_t fParticle; // Geant3 particle type | |
62 | Float_t fTheta ; // Incident theta angle in degrees | |
63 | Float_t fPhi ; // Incident phi angle in degrees | |
64 | Float_t fTlength; // Track length inside the chamber | |
65 | Float_t fEloss; // ionisation energy loss in gas | |
66 | Float_t fPHfirst; // first padhit | |
67 | Float_t fPHlast; // last padhit | |
68 | Float_t fLoss; // did it hit the freon? | |
69 | Float_t fMomX; // x Momentum at photochatode entry point | |
70 | Float_t fMomY; // y Momentum at photochatode entry point | |
71 | Float_t fMomZ; // z Momentum at photochatode entry point | |
72 | Float_t fNPads; // Pads hit | |
73 | Float_t fCerenkovAngle; // Dummy cerenkov angle | |
74 | Float_t fMomFreoX; // x Momentum at freon entry point | |
75 | Float_t fMomFreoY; // y Momentum at freon entry point | |
76 | Float_t fMomFreoZ; // z Momentum at freon entry point | |
77 | ||
78 | ClassDef(AliRICHhit,1) //RICH hit class | |
79 | };//class AliRICHhit | |
80 | //______________________________________________________________________________ | |
81 | AliRICHhit::AliRICHhit() | |
82 | :AliHit() | |
83 | {//default ctor | |
84 | fChamber=-1; | |
85 | fParticle=fTheta=fPhi=fTlength=fEloss=fPHfirst=fPHlast=fLoss=-1; | |
86 | fMomX=fMomY=fMomZ=fNPads=fCerenkovAngle=fMomFreoX=fMomFreoY=fMomFreoZ=-1; | |
87 | }//AliRICHhit::default ctor | |
88 | //______________________________________________________________________________ | |
89 | AliRICHhit::AliRICHhit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits): | |
90 | AliHit(shunt, track) | |
91 | {//ctor | |
92 | fChamber=vol[0]; | |
93 | fParticle=hits[0]; | |
94 | fX=hits[1]; | |
95 | fY=hits[2]; | |
96 | fZ=hits[3]; | |
97 | fTheta=hits[4]; | |
98 | fPhi=hits[5]; | |
99 | fTlength=hits[6]; | |
100 | fEloss=hits[7]; | |
101 | fPHfirst=(Int_t) hits[8]; | |
102 | fPHlast=(Int_t) hits[9]; | |
103 | fLoss=hits[13]; | |
104 | fMomX=hits[14]; | |
105 | fMomY=hits[15]; | |
106 | fMomZ=hits[16]; | |
107 | fNPads=hits[17]; | |
108 | fCerenkovAngle=hits[18]; | |
109 | fMomFreoX=hits[19]; | |
110 | fMomFreoY=hits[20]; | |
111 | fMomFreoZ=hits[21]; | |
112 | }//AliRICHhit::ctor | |
113 | //______________________________________________________________________________ | |
114 | //______________________________________________________________________________ | |
115 | //______________________________________________________________________________ | |
116 | class AliRICHCerenkov: public AliHit | |
117 | { | |
118 | public: | |
119 | inline AliRICHCerenkov(); | |
120 | inline AliRICHCerenkov(Int_t fIshunt, Int_t track, Int_t *vol, Float_t *Cerenkovs); | |
121 | virtual ~AliRICHCerenkov() {;} | |
122 | public: | |
123 | Int_t fChamber; // Chamber number | |
124 | Float_t fTheta ; // Incident theta angle in degrees | |
125 | Float_t fPhi ; // Incident phi angle in degrees | |
126 | Float_t fTlength; // Track length inside the chamber | |
127 | Float_t fEloss; // ionisation energy loss in gas | |
128 | Int_t fPHfirst; // first padhit | |
129 | Int_t fPHlast; // last padhit | |
130 | Int_t fCMother; // index of mother particle | |
131 | Float_t fLoss; // nature of particle loss | |
132 | Float_t fIndex; // Index of photon | |
133 | Float_t fProduction; // Point of production | |
134 | Float_t fMomX; // Local Momentum | |
135 | Float_t fMomY; // Local Momentum | |
136 | Float_t fMomZ; // Local Momentum | |
137 | Float_t fNPads; // Pads hit | |
138 | Float_t fCerenkovAngle; // Cerenkov Angle | |
139 | ||
140 | ClassDef(AliRICHCerenkov,1) //RICH cerenkov class | |
141 | };//class AliRICHCerenkov | |
142 | //______________________________________________________________________________ | |
143 | AliRICHCerenkov::AliRICHCerenkov() | |
144 | {//ctor | |
145 | fChamber=-1; | |
146 | fX=fY=fZ=fTheta=fPhi=fTlength=fEloss=-1; | |
147 | fPHfirst=fPHlast=fCMother=-1; | |
148 | fLoss=fIndex=fProduction=fMomX=fMomY=fMomZ=fNPads=fCerenkovAngle=-1; | |
149 | }//AliRICHCerenkov::ctor | |
150 | //______________________________________________________________________________ | |
151 | AliRICHCerenkov::AliRICHCerenkov(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits) | |
152 | :AliHit(shunt, track) | |
153 | {//ctor | |
154 | fChamber=vol[0]; | |
155 | fX=hits[1]; | |
156 | fY=hits[2]; | |
157 | fZ=hits[3]; | |
158 | fTheta=hits[4]; | |
159 | fPhi=hits[5]; | |
160 | fTlength=hits[6]; | |
161 | fEloss=hits[7]; | |
162 | fPHfirst=(Int_t) hits[8]; | |
163 | fPHlast=(Int_t) hits[9]; | |
164 | fCMother=Int_t(hits[10]); | |
165 | fIndex = hits[11]; | |
166 | fProduction = hits[12]; | |
167 | fLoss=hits[13]; | |
168 | fMomX=hits[14]; | |
169 | fMomY=hits[15]; | |
170 | fMomZ=hits[16]; | |
171 | fNPads=hits[17]; | |
172 | fCerenkovAngle=hits[18]; | |
173 | }//AliRICHCerenkov::ctor | |
174 | //______________________________________________________________________________ | |
175 | //______________________________________________________________________________ | |
176 | //______________________________________________________________________________ | |
177 | class AliRICHdigit :public AliDigit | |
178 | { | |
179 | public: | |
180 | AliRICHdigit() {fPadX=fPadY=fChamber=fAdc=fTracks[0]=fTracks[1]=fTracks[2]=-1;} | |
181 | inline AliRICHdigit(Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT1,Int_t iT2,Int_t iT3); | |
182 | virtual ~AliRICHdigit() {;} | |
183 | Int_t C() const{return fChamber;} | |
184 | Int_t X() const{return fPadX;} | |
185 | Int_t Y() const{return fPadY;} | |
186 | Int_t Adc() const{return fAdc;} | |
187 | protected: | |
188 | Int_t fChamber; //module number | |
189 | Int_t fPadX; //pad number along X | |
190 | Int_t fPadY; //pad number along Y | |
191 | Int_t fAdc; //ADC value | |
192 | ClassDef(AliRICHdigit,1) //RICH digit class | |
193 | };//class AliRICHdigit | |
194 | //______________________________________________________________________________ | |
195 | AliRICHdigit::AliRICHdigit(Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT0,Int_t iT1,Int_t iT2) | |
196 | { | |
197 | fChamber=iC;fPadX=iX;fPadY=iY;fAdc=iAdc; | |
198 | fTracks[0]=iT0;fTracks[1]=iT1;fTracks[2]=iT2; | |
199 | }//AliRICHdigit::ctor | |
200 | //______________________________________________________________________________ | |
201 | //______________________________________________________________________________ | |
202 | //______________________________________________________________________________ | |
3ea9cb08 | 203 | class AliRICHParam; |
204 | class AliRICHSDigit; | |
205 | ||
2f614988 | 206 | class AliRICH : public AliDetector |
207 | { | |
2f614988 | 208 | public: |
853634d3 | 209 | AliRICH(); |
210 | AliRICH(const char *name, const char *title); | |
942194a4 | 211 | AliRICH(const AliRICH& RICH):AliDetector(RICH) {;} |
853634d3 | 212 | virtual ~AliRICH(); |
dfb4e77d | 213 | |
942194a4 | 214 | AliRICH& operator=(const AliRICH&) {return *this;} |
3ea9cb08 | 215 | virtual Int_t IsVersion()const =0; |
216 | AliRICHParam *Param() {return fpParam;} | |
853634d3 | 217 | inline void AddHit(Int_t track, Int_t *vol, Float_t *hits);//virtual |
218 | inline void AddCerenkov(Int_t track, Int_t *vol, Float_t *cerenkovs); | |
219 | inline void AddSDigit(Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT0,Int_t iT1=-1,Int_t iT2=-1); | |
853634d3 | 220 | inline void ResetHits(); //virtual |
221 | inline void ResetSDigits(); //virtual | |
853634d3 | 222 | void Hits2SDigits(); //virtual |
3ea9cb08 | 223 | |
224 | TClonesArray *SDigits() const{return fSDigits;} | |
225 | TClonesArray *Cerenkovs() const{return fCerenkovs;} | |
226 | void CreateChambers(); | |
dfb4e77d | 227 | virtual void CreateMaterials(); //GEANT materials definition |
228 | Float_t AbsoCH4(Float_t x); | |
229 | Float_t Fresnel(Float_t ene,Float_t pdoti, Bool_t pola); | |
230 | virtual void BuildGeometry(); //TNode ROOT variant for event display | |
231 | virtual void CreateGeometry(); //GEANT volumes tree for simulation | |
232 | virtual void StepManager()=0; | |
3ea9cb08 | 233 | |
234 | AliRICHChamber* C(Int_t i) const{return (AliRICHChamber*)fChambers->At(i);}//return pointer to requested chamber | |
235 | AliRICHChamber& Chamber(Int_t id) {return *((AliRICHChamber *) (*fChambers)[id]);} | |
236 | TObjArray* Chambers() const{return fChambers;} | |
237 | ||
238 | void AddDigits(Int_t id, Int_t *tracks, Int_t *charges, Int_t *digits); | |
239 | void AddRawCluster(Int_t id, const AliRICHRawCluster& cluster); | |
240 | void AddRecHit1D(Int_t id, Float_t* rechit, Float_t* photons, Int_t* padsx, Int_t* padsy); | |
241 | void AddRecHit3D(Int_t id, Float_t* rechit, Float_t omega, Float_t theta, Float_t phi); | |
242 | void ResetDigits(); //virtual | |
243 | void ResetRawClusters(); | |
244 | void ResetRecHits1D(); | |
245 | void ResetRecHits3D(); | |
246 | virtual void FindClusters(Int_t nev); | |
247 | Int_t Hits2SDigits(Float_t xhit,Float_t yhit,Float_t eloss,Int_t id, ResponseType res);//kir ????? to be removed | |
248 | virtual void SDigits2Digits(); | |
249 | virtual void Digits2Reco(); | |
942194a4 | 250 | Int_t DistancetoPrimitive(Int_t /*px*/, Int_t /*py*/) {return 9999;} |
2f614988 | 251 | |
dfb4e77d | 252 | virtual void MakeBranch(Option_t *opt=" "); |
253 | virtual void MakeBranchInTreeD(TTree *treeD, const char *file=0); | |
254 | virtual void SetTreeAddress(); | |
2f614988 | 255 | |
853634d3 | 256 | |
2f614988 | 257 | |
853634d3 | 258 | AliRICHSDigit* FirstPad(AliRICHhit *hit, TClonesArray *clusters); |
dfb4e77d | 259 | AliRICHSDigit* NextPad(TClonesArray *clusters); |
2f614988 | 260 | |
c11ee6f0 | 261 | |
3ea9cb08 | 262 | void SetGeometryModel(Int_t iC,AliRICHGeometry *pRICHGeo) {C(iC)->SetGeometryModel(pRICHGeo);} |
263 | void SetSegmentationModel(Int_t iC, AliSegmentation *pAliSeg) {C(iC)->SetSegmentationModel(pAliSeg);} | |
264 | void SetResponseModel(Int_t iC, AliRICHResponse *pRICHRes) {C(iC)->SetResponseModel(pRICHRes);} | |
265 | void SetReconstructionModel(Int_t iC, AliRICHClusterFinder *pRICHReco) {C(iC)->SetReconstructionModel(pRICHReco);} | |
266 | AliRICHGeometry* GetGeometryModel(Int_t iC=0) const{return C(iC)->GetGeometryModel();} | |
267 | AliSegmentation* GetSegmentationModel(Int_t iC=0) const{return C(iC)->GetSegmentationModel();} | |
268 | AliRICHResponse* GetResponseModel(Int_t iC) const{return C(iC)->GetResponseModel();} | |
dfb4e77d | 269 | |
853634d3 | 270 | //kir virtual void SetMerger(AliRICHMerger* thisMerger) {fMerger=thisMerger;} |
dfb4e77d | 271 | |
3ba5db3e | 272 | TObjArray *Dchambers() {return fDchambers;} |
273 | TObjArray *RecHits3D() const{return fRecHits3D;} | |
274 | TObjArray *RecHits1D() const{return fRecHits1D;} | |
275 | Int_t *Ndch() {return fNdch;} | |
276 | Int_t *Nrechits1D() {return fNrechits1D;} | |
277 | Int_t *Nrechits3D() {return fNrechits3D;} | |
3ba5db3e | 278 | TClonesArray *DigitsAddress(Int_t id) {return ((TClonesArray *) (*fDchambers)[id]);} |
279 | TClonesArray *RecHitsAddress1D(Int_t id) const{return ((TClonesArray *) (*fRecHits1D)[id]);} | |
280 | TClonesArray *RecHitsAddress3D(Int_t id) const{return ((TClonesArray *) (*fRecHits3D)[id]);} | |
281 | TClonesArray *RawClustAddress(Int_t id) const{return ((TClonesArray *) (*fRawClusters)[id]);} | |
dfb4e77d | 282 | |
283 | void DiagnosticsFE(Int_t evNumber1=0,Int_t evNumber2=0); // Full events | |
284 | void DiagnosticsSE(Int_t diaglevel,Int_t evNumber1=0,Int_t evNumber2=0); // Single events | |
c11ee6f0 | 285 | |
dfb4e77d | 286 | virtual void Print(Option_t *option)const; // Prints debug information |
2f614988 | 287 | |
288 | protected: | |
3ea9cb08 | 289 | AliRICHParam *fpParam; //main RICH parametrization |
dfb4e77d | 290 | TObjArray *fChambers; //! List of RICH chambers |
853634d3 | 291 | Int_t fNsdigits; //Current number of sdigits |
dfb4e77d | 292 | Int_t fNcerenkovs; //Current number of cerenkovs |
293 | TClonesArray *fSDigits; //! List of sdigits | |
294 | TObjArray *fDchambers; //! Array of lists of digits | |
295 | TClonesArray *fCerenkovs; //! List of cerenkovs | |
296 | Int_t fNdch[kNCH]; //Array of current numbers of digits | |
297 | TObjArray *fRawClusters; // !List of raw clusters | |
298 | TObjArray *fRecHits1D; // !List of rec. hits | |
299 | TObjArray *fRecHits3D; // !List of rec. hits | |
300 | Int_t fNrawch[kNCH]; //Array of current numbers of raw clusters | |
301 | Int_t fNrechits1D[kNCH]; //Array of current numbers of rec hits 1D | |
302 | Int_t fNrechits3D[kNCH]; //Array of current numbers of rec hits 3D | |
303 | ||
304 | Int_t fCkovNumber; // Number of Cerenkov photons | |
dfb4e77d | 305 | Int_t fFreonProd; // Cerenkovs produced in freon |
dfb4e77d | 306 | Int_t fFeedbacks; // Number of feedback photons |
853634d3 | 307 | //kir Int_t fCkovQuarz; // Cerenkovs crossing quartz |
308 | //kir Int_t fCkovGap; // Cerenkovs crossing gap | |
309 | //kir Int_t fCkovCsi; // Cerenkovs crossing csi | |
310 | //kir Int_t fLostRfreo; // Cerenkovs reflected in freon | |
311 | //kir Int_t fLostRquar; // Cerenkovs reflected in quartz | |
312 | //kir Int_t fLostAfreo; // Cerenkovs absorbed in freon | |
313 | //kir Int_t fLostAquarz; // Cerenkovs absorbed in quartz | |
314 | //kir Int_t fLostAmeta; // Cerenkovs absorbed in methane | |
315 | //kir Int_t fLostCsi; // Cerenkovs below csi quantum efficiency | |
316 | //kir Int_t fLostWires; // Cerenkovs lost in wires | |
317 | //kir Float_t fMipx; // x coord. of MIP | |
318 | //kir Float_t fMipy; // y coord. of MIP | |
319 | //kir Int_t fLostFresnel; // Cerenkovs lost by Fresnel reflection | |
d8a72780 | 320 | |
853634d3 | 321 | //kir Text_t *fFileName; //! File with background hits |
322 | //kir AliRICHMerger *fMerger; //! pointer to merger | |
ddae0931 | 323 | |
dfb4e77d | 324 | ClassDef(AliRICH,2) //Main RICH class |
c11ee6f0 | 325 | };//class AliRICH |
853634d3 | 326 | |
327 | //______________________________________________________________________________ | |
328 | void AliRICH::AddHit(Int_t track, Int_t *vol, Float_t *hits) | |
329 | {//Adds the current hit to the RICH hits list | |
330 | TClonesArray &tmp=*fHits; | |
331 | new(tmp[fNhits++])AliRICHhit(fIshunt,track,vol,hits); | |
332 | } | |
333 | //______________________________________________________________________________ | |
334 | void AliRICH::AddCerenkov(Int_t track, Int_t *vol, Float_t *cerenkovs) | |
335 | {//Adds the current RICH cerenkov hit to the Cerenkovs list | |
336 | TClonesArray &tmp=*fCerenkovs; | |
337 | new(tmp[fNcerenkovs++]) AliRICHCerenkov(fIshunt,track,vol,cerenkovs); | |
338 | } | |
339 | //______________________________________________________________________________ | |
340 | void AliRICH::ResetHits() | |
341 | {//Resets hits and cerenkovs | |
342 | AliDetector::ResetHits(); | |
343 | fNcerenkovs = 0; | |
344 | if(fCerenkovs)fCerenkovs->Clear(); | |
345 | } | |
346 | //______________________________________________________________________________ | |
347 | void AliRICH::AddSDigit(Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT0,Int_t iT1,Int_t iT2) | |
348 | {//Adds the current Sdigit to the RICH list of Sdigits | |
349 | TClonesArray &tmp = *fSDigits; | |
350 | new(tmp[fNsdigits++])AliRICHdigit(iC,iX,iY,iAdc,iT0,iT1,iT2); | |
351 | } | |
352 | //______________________________________________________________________________ | |
353 | void AliRICH::ResetSDigits() | |
354 | {//Resets sdigits | |
355 | fNsdigits=0; | |
356 | if(fSDigits)fSDigits->Clear(); | |
357 | } | |
358 | //______________________________________________________________________________ | |
2f614988 | 359 | #endif |