]>
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 | ||
c11ee6f0 | 7 | #include <TObjArray.h> |
853634d3 | 8 | #include <TClonesArray.h> |
c021cb15 | 9 | #include <TLorentzVector.h> |
c11ee6f0 | 10 | #include <AliDetector.h> |
853634d3 | 11 | #include <AliHit.h> |
12 | #include <AliDigit.h> | |
543d5224 | 13 | |
dfb4e77d | 14 | #include "AliRICHConst.h" |
15 | #include "AliRICHChamber.h" | |
853634d3 | 16 | |
c60862bf | 17 | #include "AliRICHDigit.h" |
18 | #include "AliRICHSDigit.h" | |
543d5224 | 19 | #include "AliRICHRawCluster.h" |
a4622d0f | 20 | class AliRICHRecHit1D; |
853634d3 | 21 | |
c021cb15 | 22 | //__________________AliRICHhit______________________________________________________________________ |
23 | //__________________________________________________________________________________________________ | |
24 | //__________________________________________________________________________________________________ | |
853634d3 | 25 | class AliRICHhit : public AliHit |
26 | { | |
27 | public: | |
28 | inline AliRICHhit(); | |
29 | inline AliRICHhit(Int_t fIshunt, Int_t track, Int_t *vol, Float_t *hits); | |
c021cb15 | 30 | inline AliRICHhit(Int_t track,Int_t iPID,Int_t iChamber,TLorentzVector x4,Float_t eloss); |
853634d3 | 31 | virtual ~AliRICHhit() {;} |
32 | ||
c60862bf | 33 | Int_t C() {return fChamber;} |
853634d3 | 34 | Int_t Chamber() {return fChamber;} |
c60862bf | 35 | Int_t Pid() {return fPid;} |
36 | Int_t Particle() {return fPid;} | |
853634d3 | 37 | Float_t Theta() {return fTheta;} |
38 | Float_t Phi() {return fPhi;} | |
39 | Float_t Tlength() {return fTlength;} | |
40 | Float_t Eloss() {return fEloss;} | |
41 | Float_t Loss() {return fLoss;} | |
c021cb15 | 42 | Float_t PHfirst() {return fPHfirst;} |
43 | Float_t PHlast() {return fPHlast;} | |
853634d3 | 44 | Float_t MomX() {return fMomX;} |
45 | Float_t MomY() {return fMomY;} | |
46 | Float_t MomZ() {return fMomZ;} | |
47 | Float_t CerenkovAngle() {return fMomX;} | |
48 | Float_t MomFreoX() {return fMomX;} | |
49 | Float_t MomFreoY() {return fMomY;} | |
50 | Float_t MomFreoZ() {return fMomZ;} | |
543d5224 | 51 | void Print(Option_t *option="")const; //virtual |
853634d3 | 52 | protected: |
c021cb15 | 53 | Int_t fChamber; //chamber number |
c60862bf | 54 | Int_t fPid; //particle code |
c021cb15 | 55 | Float_t fTheta,fPhi ; //incident theta phi angles in degrees |
56 | Float_t fTlength; //track length inside the chamber | |
57 | Float_t fEloss; //ionisation energy loss in gas | |
58 | Float_t fPHfirst; //first padhit | |
59 | Float_t fPHlast; //last padhit | |
60 | Float_t fLoss; // did it hit the freon? | |
61 | Float_t fMomX,fMomY,fMomZ; //momentum at photochatode entry point | |
62 | Float_t fNPads; // Pads hit | |
63 | Float_t fCerenkovAngle; // Dummy cerenkov angle | |
64 | Float_t fMomFreoX,fMomFreoY,fMomFreoZ; //momentum at freon entry point | |
65 | ClassDef(AliRICHhit,1) //RICH hit class | |
853634d3 | 66 | };//class AliRICHhit |
c021cb15 | 67 | |
68 | //__________________________________________________________________________________________________ | |
853634d3 | 69 | AliRICHhit::AliRICHhit() |
70 | :AliHit() | |
71 | {//default ctor | |
c60862bf | 72 | fChamber=fPid=kBad; |
c021cb15 | 73 | fTheta=fPhi=fTlength=fEloss=fPHfirst=fPHlast=fLoss=kBad; |
74 | fMomX=fMomY=fMomZ=fNPads=fCerenkovAngle=fMomFreoX=fMomFreoY=fMomFreoZ=kBad; | |
543d5224 | 75 | } |
c021cb15 | 76 | //__________________________________________________________________________________________________ |
77 | AliRICHhit::AliRICHhit(Int_t shunt, Int_t track, Int_t *vol, Float_t *hit): | |
853634d3 | 78 | AliHit(shunt, track) |
79 | {//ctor | |
80 | fChamber=vol[0]; | |
c60862bf | 81 | fPid=(Int_t)hit[0]; |
c021cb15 | 82 | fX=hit[1];fY=hit[2];fZ=hit[3]; |
83 | fTheta=hit[4];fPhi=hit[5]; | |
84 | fTlength=hit[6]; | |
85 | fEloss=hit[7]; | |
86 | fPHfirst=(Int_t)hit[8]; | |
87 | fPHlast=(Int_t)hit[9]; | |
88 | fLoss=hit[13]; | |
89 | fMomX=hit[14];fMomY=hit[15];fMomZ=hit[16]; | |
90 | fNPads=hit[17]; | |
91 | fCerenkovAngle=hit[18]; | |
92 | fMomFreoX=hit[19];fMomFreoY=hit[20];fMomFreoZ=hit[21]; | |
543d5224 | 93 | } |
c021cb15 | 94 | //__________________________________________________________________________________________________ |
95 | AliRICHhit::AliRICHhit(Int_t track,Int_t iPID,Int_t iChamber,TLorentzVector x4,Float_t eloss): | |
96 | AliHit(0, track) | |
97 | {//ctor | |
98 | fChamber=iChamber; | |
c60862bf | 99 | fPid=iPID; |
c021cb15 | 100 | fX=x4.X();fY=x4.Y();fZ=x4.Z(); |
101 | fEloss=eloss; | |
543d5224 | 102 | } |
c021cb15 | 103 | |
104 | //__________________AliRICHCerenkov_________________________________________________________________ | |
105 | //__________________________________________________________________________________________________ | |
106 | //__________________________________________________________________________________________________ | |
853634d3 | 107 | class AliRICHCerenkov: public AliHit |
108 | { | |
109 | public: | |
110 | inline AliRICHCerenkov(); | |
111 | inline AliRICHCerenkov(Int_t fIshunt, Int_t track, Int_t *vol, Float_t *Cerenkovs); | |
112 | virtual ~AliRICHCerenkov() {;} | |
113 | public: | |
c021cb15 | 114 | Int_t fChamber; //chamber number |
115 | Float_t fTheta,fPhi; //incident theta phi angles in degrees | |
116 | Float_t fTlength; //track length inside the chamber | |
117 | Float_t fEloss; //ionisation energy loss in gas | |
118 | Int_t fPHfirst; //first padhit | |
119 | Int_t fPHlast; //last padhit | |
120 | Int_t fCMother; //index of mother particle | |
121 | Float_t fLoss; //nature of particle loss | |
122 | Float_t fIndex; //index of photon | |
123 | Float_t fProduction; //point of production | |
124 | Float_t fMomX,fMomY,fMomZ; //local Momentum | |
853634d3 | 125 | Float_t fNPads; // Pads hit |
126 | Float_t fCerenkovAngle; // Cerenkov Angle | |
127 | ||
128 | ClassDef(AliRICHCerenkov,1) //RICH cerenkov class | |
129 | };//class AliRICHCerenkov | |
c021cb15 | 130 | |
131 | //__________________________________________________________________________________________________ | |
853634d3 | 132 | AliRICHCerenkov::AliRICHCerenkov() |
133 | {//ctor | |
c021cb15 | 134 | fChamber=kBad; |
135 | fX=fY=fZ=fTheta=fPhi=fTlength=fEloss=kBad; | |
136 | fPHfirst=fPHlast=fCMother=kBad; | |
137 | fLoss=fIndex=fProduction=fMomX=fMomY=fMomZ=fNPads=fCerenkovAngle=kBad; | |
543d5224 | 138 | } |
c021cb15 | 139 | //__________________________________________________________________________________________________ |
853634d3 | 140 | AliRICHCerenkov::AliRICHCerenkov(Int_t shunt, Int_t track, Int_t *vol, Float_t *hits) |
141 | :AliHit(shunt, track) | |
142 | {//ctor | |
143 | fChamber=vol[0]; | |
c021cb15 | 144 | fX=hits[1];fY=hits[2];fZ=hits[3]; |
145 | fTheta=hits[4];fPhi=hits[5]; | |
853634d3 | 146 | fTlength=hits[6]; |
147 | fEloss=hits[7]; | |
c021cb15 | 148 | fPHfirst=(Int_t)hits[8];fPHlast=(Int_t)hits[9]; |
853634d3 | 149 | fCMother=Int_t(hits[10]); |
150 | fIndex = hits[11]; | |
151 | fProduction = hits[12]; | |
152 | fLoss=hits[13]; | |
c021cb15 | 153 | fMomX=hits[14];fMomY=hits[15];fMomZ=hits[16]; |
853634d3 | 154 | fNPads=hits[17]; |
155 | fCerenkovAngle=hits[18]; | |
543d5224 | 156 | } |
c021cb15 | 157 | |
158 | //__________________AliRICHdigit____________________________________________________________________ | |
159 | //__________________________________________________________________________________________________ | |
160 | //__________________________________________________________________________________________________ | |
853634d3 | 161 | class AliRICHdigit :public AliDigit |
162 | { | |
163 | public: | |
c60862bf | 164 | AliRICHdigit() {fPadX=fPadY=fChamber=fQdc=fTracks[0]=fTracks[1]=fTracks[2]=kBad;} |
165 | inline AliRICHdigit(Int_t iC,Int_t iX,Int_t iY,Int_t iQdc,Int_t iT1,Int_t iT2,Int_t iT3); | |
853634d3 | 166 | virtual ~AliRICHdigit() {;} |
c60862bf | 167 | inline Int_t Compare(const TObject *pObj) const;//virtual |
168 | Bool_t IsSortable() const{return kTRUE;}//virtual | |
169 | Int_t C() const{return fChamber;} | |
170 | Int_t X() const{return fPadX;} | |
171 | Int_t Y() const{return fPadY;} | |
172 | Int_t Id() const{return fChamber*1000000+fPadX*1000+fPadY;} | |
173 | Int_t Qdc() const{return fQdc;} | |
174 | Int_t T(Int_t i) const{return fTracks[i];} | |
175 | void Print(Option_t *option)const; //virtual | |
853634d3 | 176 | protected: |
177 | Int_t fChamber; //module number | |
c021cb15 | 178 | Int_t fPadX; //pad number along X |
179 | Int_t fPadY; //pad number along Y | |
c60862bf | 180 | Int_t fQdc; //ADC value |
853634d3 | 181 | ClassDef(AliRICHdigit,1) //RICH digit class |
182 | };//class AliRICHdigit | |
c021cb15 | 183 | //__________________________________________________________________________________________________ |
c60862bf | 184 | AliRICHdigit::AliRICHdigit(Int_t iC,Int_t iX,Int_t iY,Int_t iQdc,Int_t iT0,Int_t iT1,Int_t iT2) |
853634d3 | 185 | { |
c60862bf | 186 | fChamber=iC;fPadX=iX;fPadY=iY;fQdc=iQdc; |
853634d3 | 187 | fTracks[0]=iT0;fTracks[1]=iT1;fTracks[2]=iT2; |
543d5224 | 188 | } |
c60862bf | 189 | //__________________________________________________________________________________________________ |
190 | Int_t AliRICHdigit::Compare(const TObject *pObj)const | |
191 | { | |
192 | if(Id()==((AliRICHdigit*)pObj)->Id()) | |
193 | return 0; | |
194 | else if(Id()>((AliRICHdigit*)pObj)->Id()) | |
195 | return 1; | |
196 | else | |
197 | return -1; | |
543d5224 | 198 | } |
c021cb15 | 199 | //__________________AliRICH_________________________________________________________________________ |
200 | //__________________________________________________________________________________________________ | |
201 | //__________________________________________________________________________________________________ | |
3ea9cb08 | 202 | class AliRICHParam; |
203 | class AliRICHSDigit; | |
204 | ||
2f614988 | 205 | class AliRICH : public AliDetector |
206 | { | |
2f614988 | 207 | public: |
853634d3 | 208 | AliRICH(); |
209 | AliRICH(const char *name, const char *title); | |
942194a4 | 210 | AliRICH(const AliRICH& RICH):AliDetector(RICH) {;} |
853634d3 | 211 | virtual ~AliRICH(); |
dfb4e77d | 212 | |
942194a4 | 213 | AliRICH& operator=(const AliRICH&) {return *this;} |
c60862bf | 214 | virtual Int_t IsVersion() const =0; |
215 | void Hits2SDigits(); //virtual | |
216 | void SDigits2Digits(); //virtual | |
217 | void Digits2Reco(); //virtual | |
c021cb15 | 218 | |
c60862bf | 219 | inline void CreateHits(); |
220 | inline void CreateSdigits(); | |
221 | inline void CreateDigits(); | |
222 | inline void CreateClusters(); | |
223 | inline void AddHit(Int_t track, Int_t *vol, Float_t *hits); //virtual | |
c60862bf | 224 | inline void AddSdigit(Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT0,Int_t iT1=kBad,Int_t iT2=kBad); |
543d5224 | 225 | inline void AddDigit (Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT0,Int_t iT1=kBad,Int_t iT2=kBad); |
c60862bf | 226 | void AddCluster() {;} |
543d5224 | 227 | void ResetHits() {AliDetector::ResetHits();fNcerenkovs=0;if(fCerenkovs)fCerenkovs->Clear();fNspecials=0;if(fSpecials)fSpecials->Clear();} //virtual |
228 | void ResetSdigits() {fNsdigits=0; if(fSdigits) fSdigits ->Clear();} | |
229 | void ResetDigits() {if(fDigitsNew)for(int i=0;i<kNCH;i++){fDigitsNew->At(i)->Clear();fNdigitsNew[i]=0;}} | |
230 | void ResetClusters() {if(fClusters) for(int i=0;i<kNCH;i++){fClusters ->At(i)->Clear();fNclusters[i]=0;}} | |
c60862bf | 231 | //Hits provided by AliDetector |
232 | TClonesArray* Sdigits() const{return fSdigits;} | |
543d5224 | 233 | TClonesArray* Digits(Int_t iC) const{if(fDigitsNew) return (TClonesArray *)fDigitsNew->At(iC-1);else return 0;} |
c60862bf | 234 | TClonesArray* Clusters(Int_t iC) const{if(fClusters) return (TClonesArray *)fClusters->At(iC-1);else return 0;} |
235 | ||
236 | AliRICHChamber* C(Int_t iC) const{return (AliRICHChamber*)fChambers->At(iC-1);} | |
237 | AliRICHParam* Param() const{return fpParam;} | |
c021cb15 | 238 | |
c60862bf | 239 | AliRICHhit* FirstHit(Int_t iTrkN) {return (AliRICHhit*)AliDetector::FirstHit(iTrkN);} //virtual |
240 | AliRICHhit* NextHit() {return (AliRICHhit*)AliDetector::NextHit();} //virtual | |
241 | ||
3ea9cb08 | 242 | void CreateChambers(); |
c60862bf | 243 | void CreateMaterials(); //virtual |
244 | virtual void BuildGeometry(); //virtual | |
245 | virtual void CreateGeometry(); //virtual | |
dfb4e77d | 246 | Float_t AbsoCH4(Float_t x); |
247 | Float_t Fresnel(Float_t ene,Float_t pdoti, Bool_t pola); | |
3ea9cb08 | 248 | |
c021cb15 | 249 | virtual void StepManager()=0; |
c60862bf | 250 | void GenerateFeedbacks(Int_t iChamber,Float_t eloss); |
251 | void Print(Option_t *option)const;//virtual | |
252 | void MakeBranch(Option_t *opt=" "); | |
253 | void SetTreeAddress();//virtual | |
543d5224 | 254 | // OLD staff OLD staff |
255 | inline void AddCerenkov(Int_t track, Int_t *vol, Float_t *cerenkovs); | |
256 | inline void AddSpecialOld(Int_t *); | |
257 | inline void AddDigitOld(Int_t id, Int_t *tracks, Int_t *charges, Int_t *digits); | |
258 | inline void AddClusterOld(Int_t iChamber, const AliRICHRawCluster& cluster); | |
259 | void AddRecHit1D(Int_t id, Float_t* rechit, Float_t* photons, Int_t* padsx, Int_t* padsy); | |
543d5224 | 260 | |
261 | inline void CreateCerenkovsOld(); | |
262 | inline void CreateSpecialsOld(); | |
263 | inline void CreateDigitsOld(); | |
264 | inline void CreateRawClustersOld(); | |
265 | inline void CreateRecos1Old(); | |
543d5224 | 266 | void ResetDigitsOld() {if(fDchambers) for(int i=0;i<kNCH;i++){fDchambers->At(i)->Clear();fNdch[i]=0;}} //virtual |
267 | void ResetSpecialsOld(){fNspecials=0; if(fSpecials) fSpecials->Clear();} | |
268 | void ResetRawClusters(){if(fRawClusters)for(int i=0;i<kNCH;i++){fRawClusters->At(i)->Clear();fNrawch[i]=0;}} | |
269 | void ResetRecHits1D() {if(fRecHits1D) for(int i=0;i<kNCH;i++){fRecHits1D ->At(i)->Clear();fNrechits1D[i]=0;}} | |
543d5224 | 270 | |
271 | TClonesArray* DigitsOld(Int_t iC) const{if(fDchambers) return (TClonesArray *)fDchambers->At(iC-1);else return 0;} | |
272 | TClonesArray* ClustersOld(Int_t iC) const{if(fRawClusters)return (TClonesArray *)fRawClusters->At(iC-1);else return 0;} | |
273 | TClonesArray* Specials() const{return fSpecials;} | |
274 | TClonesArray* Cerenkovs() const{return fCerenkovs;} | |
c60862bf | 275 | |
3ea9cb08 | 276 | |
c60862bf | 277 | |
278 | AliRICHChamber& Chamber(Int_t id) {return *((AliRICHChamber *) (*fChambers)[id]);} | |
c60862bf | 279 | TObjArray *Dchambers() const{return fDchambers;} |
3ba5db3e | 280 | TObjArray *RecHits1D() const{return fRecHits1D;} |
281 | Int_t *Ndch() {return fNdch;} | |
282 | Int_t *Nrechits1D() {return fNrechits1D;} | |
c60862bf | 283 | TClonesArray *DigitsAddress(Int_t id) const{return ((TClonesArray *) (*fDchambers)[id]);} |
3ba5db3e | 284 | TClonesArray *RecHitsAddress1D(Int_t id) const{return ((TClonesArray *) (*fRecHits1D)[id]);} |
3ba5db3e | 285 | TClonesArray *RawClustAddress(Int_t id) const{return ((TClonesArray *) (*fRawClusters)[id]);} |
c60862bf | 286 | // Int_t DistancetoPrimitive(Int_t /*px*/, Int_t /*py*/) {return 9999;} |
2f614988 | 287 | |
c60862bf | 288 | protected: |
c021cb15 | 289 | AliRICHParam *fpParam; //main RICH parametrization |
290 | TObjArray *fChambers; //list of RICH chambers | |
c60862bf | 291 | //fHits and fDigits belong to AliDetector |
292 | TClonesArray *fSdigits; //! List of sdigits | |
853634d3 | 293 | Int_t fNsdigits; //Current number of sdigits |
c60862bf | 294 | TObjArray *fDigitsNew; //! Each chamber holds it's one lists of digits |
295 | Int_t fNdigitsNew[kNCH]; //Array of current numbers of digits | |
296 | TObjArray *fClusters; //! Each chamber holds it's one lists of clusters | |
297 | Int_t fNclusters[kNCH]; //Array of current numbers of raw clusters | |
298 | ||
299 | TClonesArray *fCerenkovs; //! ??? List of cerenkovs | |
300 | Int_t fNcerenkovs; //??? Current number of cerenkovs | |
301 | TClonesArray *fSpecials; //! ??? List of specials | |
302 | Int_t fNspecials; //??? Current number of specials | |
303 | TObjArray *fDchambers; //! Array of lists of digits | |
dfb4e77d | 304 | Int_t fNdch[kNCH]; //Array of current numbers of digits |
c60862bf | 305 | TObjArray *fRawClusters; //! Array of lists of raw clusters |
dfb4e77d | 306 | Int_t fNrawch[kNCH]; //Array of current numbers of raw clusters |
c60862bf | 307 | TObjArray *fRecHits1D; //!List of rec. hits |
dfb4e77d | 308 | Int_t fNrechits1D[kNCH]; //Array of current numbers of rec hits 1D |
dfb4e77d | 309 | Int_t fCkovNumber; // Number of Cerenkov photons |
dfb4e77d | 310 | Int_t fFreonProd; // Cerenkovs produced in freon |
dfb4e77d | 311 | Int_t fFeedbacks; // Number of feedback photons |
ddae0931 | 312 | |
c60862bf | 313 | ClassDef(AliRICH,3) //Main RICH class |
314 | };//class AliRICH | |
c021cb15 | 315 | //__________________________________________________________________________________________________ |
853634d3 | 316 | void AliRICH::AddHit(Int_t track, Int_t *vol, Float_t *hits) |
317 | {//Adds the current hit to the RICH hits list | |
318 | TClonesArray &tmp=*fHits; | |
319 | new(tmp[fNhits++])AliRICHhit(fIshunt,track,vol,hits); | |
543d5224 | 320 | } |
c60862bf | 321 | //__________________________________________________________________________________________________ |
322 | void AliRICH::AddSdigit(Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT0,Int_t iT1,Int_t iT2) | |
323 | {//Adds the current Sdigit to the RICH list of Sdigits | |
324 | TClonesArray &tmp=*fSdigits; | |
325 | new(tmp[fNsdigits++])AliRICHdigit(iC,iX,iY,iAdc,iT0,iT1,iT2); | |
326 | } | |
327 | //__________________________________________________________________________________________________ | |
c60862bf | 328 | void AliRICH::AddDigit(Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT0,Int_t iT1,Int_t iT2) |
329 | {//Adds the current digit to the corresponding RICH list of digits (individual list per chamber) | |
330 | TClonesArray &tmp=*((TClonesArray*)fDigitsNew->At(iC-1)); | |
331 | new(tmp[fNdigitsNew[iC-1]++]) AliRICHdigit(iC,iX,iY,iAdc,iT0,iT1,iT2); | |
543d5224 | 332 | } |
c60862bf | 333 | //__________________________________________________________________________________________________ |
334 | void AliRICH::CreateHits() | |
335 | { | |
336 | if(fHits) return; | |
337 | if(GetDebug())Info("CreateHits","creating hits container."); | |
338 | fHits=new TClonesArray("AliRICHhit",10000); fNhits=0; | |
543d5224 | 339 | } |
c60862bf | 340 | //__________________________________________________________________________________________________ |
341 | void AliRICH::CreateSdigits() | |
342 | { | |
343 | if(fSdigits) return; | |
344 | if(GetDebug())Info("CreateSdigits","creating sdigits container."); | |
345 | fSdigits=new TClonesArray("AliRICHdigit",10000); fNsdigits=0; | |
543d5224 | 346 | } |
c60862bf | 347 | //__________________________________________________________________________________________________ |
348 | void AliRICH::CreateDigits() | |
349 | { | |
350 | if(fDigitsNew) return; | |
351 | if(GetDebug())Info("CreateDigits","creating digits containers."); | |
352 | fDigitsNew = new TObjArray(kNCH); | |
353 | for(Int_t i=0;i<kNCH;i++) {fDigitsNew->AddAt(new TClonesArray("AliRICHdigit",10000), i); fNdigitsNew[i]=0;} | |
543d5224 | 354 | } |
c60862bf | 355 | //__________________________________________________________________________________________________ |
356 | void AliRICH::CreateClusters() | |
357 | { | |
358 | if(fClusters) return; | |
359 | if(GetDebug())Info("CreateClusters","creating clusters containers."); | |
360 | fClusters = new TObjArray(kNCH); | |
361 | for(Int_t i=0;i<kNCH;i++) {fClusters->AddAt(new TClonesArray("AliRICHcluster",10000), i); fNclusters[i]=0;} | |
543d5224 | 362 | } |
c60862bf | 363 | |
364 | ||
365 | ||
366 | //__________________________________________________________________________________________________ | |
367 | void AliRICH::CreateCerenkovsOld() | |
368 | { | |
369 | if(fCerenkovs) return; | |
370 | if(GetDebug())Info("CreateCerenkovs","creating cerenkovs container."); | |
371 | fCerenkovs=new TClonesArray("AliRICHCerenkov",10000); fNcerenkovs=0; | |
543d5224 | 372 | } |
c60862bf | 373 | //__________________________________________________________________________________________________ |
374 | void AliRICH::CreateSpecialsOld() | |
375 | { | |
376 | if(fSpecials) return; | |
377 | if(GetDebug())Info("CreateSpecialsOld","creating SDigits special container."); | |
378 | fSpecials=new TClonesArray("AliRICHSDigit",100000); fNspecials=0; | |
543d5224 | 379 | } |
c60862bf | 380 | //__________________________________________________________________________________________________ |
381 | void AliRICH::CreateDigitsOld() | |
382 | { | |
383 | if(fDchambers) return; | |
384 | if(GetDebug())Info("CreateDigitsOld","creating digits containers."); | |
385 | fDchambers = new TObjArray(kNCH); | |
386 | for(Int_t i=0;i<kNCH;i++) fDchambers->AddAt(new TClonesArray("AliRICHDigit",10000), i); | |
543d5224 | 387 | } |
c60862bf | 388 | //__________________________________________________________________________________________________ |
389 | void AliRICH::CreateRawClustersOld() | |
390 | { | |
391 | if(fRawClusters) return; | |
392 | if(GetDebug())Info("CreateClustersOld","creating clusters containers."); | |
393 | fRawClusters = new TObjArray(kNCH); | |
394 | for(Int_t i=0; i<kNCH ;i++) fRawClusters->AddAt(new TClonesArray("AliRICHRawCluster",10000), i); | |
543d5224 | 395 | } |
c60862bf | 396 | //__________________________________________________________________________________________________ |
397 | void AliRICH::CreateRecos1Old() | |
398 | { | |
399 | if(fRecHits1D) return; | |
400 | if(GetDebug())Info("CreateRecos1DOld","creating recos 1 containers."); | |
401 | fRecHits1D = new TObjArray(kNCH); | |
402 | for(Int_t i=0; i<kNCH ;i++) fRecHits1D->AddAt(new TClonesArray("AliRICHRecHit1D",1000), i); | |
403 | } | |
404 | //__________________________________________________________________________________________________ | |
853634d3 | 405 | void AliRICH::AddCerenkov(Int_t track, Int_t *vol, Float_t *cerenkovs) |
406 | {//Adds the current RICH cerenkov hit to the Cerenkovs list | |
407 | TClonesArray &tmp=*fCerenkovs; | |
408 | new(tmp[fNcerenkovs++]) AliRICHCerenkov(fIshunt,track,vol,cerenkovs); | |
409 | } | |
c021cb15 | 410 | //__________________________________________________________________________________________________ |
543d5224 | 411 | void AliRICH::AddSpecialOld(Int_t *aiSDigit) |
c60862bf | 412 | {// Adds the current Sdigit to the RICH list of Specials |
413 | TClonesArray &lSDigits = *fSpecials; | |
543d5224 | 414 | new(lSDigits[fNspecials++]) AliRICHSDigit(aiSDigit); |
415 | } | |
c021cb15 | 416 | //__________________________________________________________________________________________________ |
543d5224 | 417 | void AliRICH::AddDigitOld(Int_t iChamber, Int_t *tracks, Int_t *charges, Int_t *digits) |
c60862bf | 418 | {// Adds the current digit to the RICH list of S digits |
543d5224 | 419 | TClonesArray &ldigits = *((TClonesArray*)fDchambers->At(iChamber-1)); |
420 | new(ldigits[fNdch[iChamber-1]++]) AliRICHDigit(tracks,charges,digits); | |
421 | } | |
422 | //__________________________________________________________________________________________________ | |
423 | void AliRICH::AddClusterOld(Int_t iChamber, const AliRICHRawCluster& c) | |
424 | {// Add a RICH raw cluster to the list | |
425 | TClonesArray &tmp= *((TClonesArray*)fRawClusters->At(iChamber-1)); | |
426 | new(tmp[fNrawch[iChamber-1]++]) AliRICHRawCluster(c); | |
427 | } | |
c021cb15 | 428 | //__________________________________________________________________________________________________ |
c60862bf | 429 | #endif//#ifndef AliRICH_h |