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