PatRec adapted to new IO
[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
14 #include "AliRICHConst.h"
15 #include "AliRICHChamber.h"
16
17 #include "AliRICHDigit.h"
18 #include "AliRICHSDigit.h"
19 #include "AliRICHRawCluster.h"
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 }
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 }
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 }
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 }
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 }
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 }
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
200 //__________________AliRICH_________________________________________________________________________
201 //__________________________________________________________________________________________________
202 //__________________________________________________________________________________________________
203 class AliRICHParam;
204 class AliRICHSDigit;
205
206 class AliRICH : public AliDetector 
207 {
208 public:
209             AliRICH();                                            
210             AliRICH(const char *name, const char *title);         
211             AliRICH(const AliRICH& RICH):AliDetector(RICH) {;}   
212   virtual  ~AliRICH();                                            
213           
214   AliRICH&  operator=(const AliRICH&)                 {return *this;}
215   virtual Int_t   IsVersion()                                            const =0;            
216           void    Hits2SDigits();                                                                                 //virtual
217           void    SDigits2Digits();                                                                               //virtual
218           void    Digits2Reco();                                                                                  //virtual
219   
220   inline  void    CreateHits();    
221   inline  void    CreateSdigits();  
222   inline  void    CreateDigits();  
223   inline  void    CreateClusters();  
224   inline  void    AddHit(Int_t track, Int_t *vol, Float_t *hits);                                                 //virtual
225   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);
226   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);     
227           void    AddCluster()                                                                                {;}     
228           void    ResetHits()     {AliDetector::ResetHits();fNcerenkovs=0;if(fCerenkovs)fCerenkovs->Clear();fNspecials=0;if(fSpecials)fSpecials->Clear();}  //virtual
229           void    ResetSdigits()  {fNsdigits=0;  if(fSdigits)  fSdigits ->Clear();}                                 
230           void    ResetDigits()   {if(fDigitsNew)for(int i=0;i<kNCH;i++){fDigitsNew->At(i)->Clear();fNdigitsNew[i]=0;}}
231           void    ResetClusters() {if(fClusters) for(int i=0;i<kNCH;i++){fClusters ->At(i)->Clear();fNclusters[i]=0;}}
232                   //Hits provided by AliDetector
233   TClonesArray*   Sdigits()             const{return fSdigits;}
234   TClonesArray*   Digits(Int_t iC)      const{if(fDigitsNew) return (TClonesArray *)fDigitsNew->At(iC-1);else return 0;}
235   TClonesArray*   Clusters(Int_t iC)    const{if(fClusters)  return (TClonesArray *)fClusters->At(iC-1);else return 0;}
236           
237   AliRICHChamber* C(Int_t iC)           const{return (AliRICHChamber*)fChambers->At(iC-1);}
238   AliRICHParam*   Param()               const{return fpParam;}
239   
240   AliRICHhit*     FirstHit(Int_t iTrkN)      {return (AliRICHhit*)AliDetector::FirstHit(iTrkN);}                   //virtual
241   AliRICHhit*     NextHit()                  {return (AliRICHhit*)AliDetector::NextHit();}                         //virtual 
242   
243           void    CreateChambers();         
244           void    CreateMaterials(); //virtual
245   virtual void    BuildGeometry();   //virtual
246   virtual void    CreateGeometry();  //virtual
247           Float_t AbsoCH4(Float_t x);
248           Float_t Fresnel(Float_t ene,Float_t pdoti, Bool_t pola);
249   
250   virtual void    StepManager()=0;
251           void    GenerateFeedbacks(Int_t iChamber,Float_t eloss);
252           void    Print(Option_t *option)const;//virtual
253           void    MakeBranch(Option_t *opt=" ");
254           void    SetTreeAddress();//virtual
255 // OLD staff OLD staff    
256   inline  void    AddCerenkov(Int_t track, Int_t *vol, Float_t *cerenkovs);
257   inline  void    AddSpecialOld(Int_t *);      
258   inline  void    AddDigitOld(Int_t id, Int_t *tracks, Int_t *charges, Int_t *digits);          
259   inline  void    AddClusterOld(Int_t iChamber, const AliRICHRawCluster& cluster);
260           void    AddRecHit1D(Int_t id, Float_t* rechit, Float_t* photons, Int_t* padsx, Int_t* padsy);
261           void    AddRecHit3D(Int_t id, Float_t* rechit, Float_t omega, Float_t theta, Float_t phi);
262           
263   inline  void    CreateCerenkovsOld();  
264   inline  void    CreateSpecialsOld();   
265   inline  void    CreateDigitsOld();    
266   inline  void    CreateRawClustersOld();  
267   inline  void    CreateRecos1Old();  
268   inline  void    CreateRecos3Old();  
269           void    ResetDigitsOld()  {if(fDchambers)  for(int i=0;i<kNCH;i++){fDchambers->At(i)->Clear();fNdch[i]=0;}}   //virtual
270           void    ResetSpecialsOld(){fNspecials=0; if(fSpecials) fSpecials->Clear();}   
271           void    ResetRawClusters(){if(fRawClusters)for(int i=0;i<kNCH;i++){fRawClusters->At(i)->Clear();fNrawch[i]=0;}}
272           void    ResetRecHits1D()  {if(fRecHits1D)  for(int i=0;i<kNCH;i++){fRecHits1D  ->At(i)->Clear();fNrechits1D[i]=0;}}
273           void    ResetRecHits3D()  {if(fRecHits3D)  for(int i=0;i<kNCH;i++){fRecHits3D  ->At(i)->Clear();fNrechits3D[i]=0;}}
274   
275   TClonesArray*   DigitsOld(Int_t iC)   const{if(fDchambers) return (TClonesArray *)fDchambers->At(iC-1);else return 0;}
276   TClonesArray*   ClustersOld(Int_t iC) const{if(fRawClusters)return (TClonesArray *)fRawClusters->At(iC-1);else return 0;}
277   TClonesArray*   Specials()            const{return fSpecials;}
278   TClonesArray*   Cerenkovs()           const{return fCerenkovs;}
279   
280   
281             
282   AliRICHChamber& Chamber(Int_t id)      {return *((AliRICHChamber *) (*fChambers)[id]);}  
283   TObjArray     *Dchambers()                const{return fDchambers;}
284   TObjArray     *RecHits3D()                const{return fRecHits3D;}
285   TObjArray     *RecHits1D()                const{return fRecHits1D;}
286   Int_t         *Ndch()                          {return fNdch;}
287   Int_t         *Nrechits1D()                    {return fNrechits1D;} 
288   Int_t         *Nrechits3D()                    {return fNrechits3D;} 
289   TClonesArray  *DigitsAddress(Int_t id)    const{return ((TClonesArray *) (*fDchambers)[id]);}
290   TClonesArray  *RecHitsAddress1D(Int_t id) const{return ((TClonesArray *) (*fRecHits1D)[id]);}
291   TClonesArray  *RecHitsAddress3D(Int_t id) const{return ((TClonesArray *) (*fRecHits3D)[id]);}
292   TClonesArray  *RawClustAddress(Int_t id)  const{return ((TClonesArray *) (*fRawClusters)[id]);}    
293 //          Int_t DistancetoPrimitive(Int_t /*px*/, Int_t /*py*/)      {return 9999;}
294     
295 protected:  
296   AliRICHParam         *fpParam;             //main RICH parametrization     
297   TObjArray            *fChambers;           //list of RICH chambers
298                 //fHits and fDigits belong to AliDetector
299   TClonesArray         *fSdigits;            //! List of sdigits  
300   Int_t                 fNsdigits;           //Current number of sdigits
301   TObjArray            *fDigitsNew;          //! Each chamber holds it's one lists of digits
302   Int_t                 fNdigitsNew[kNCH];   //Array of current numbers of digits
303   TObjArray            *fClusters;           //! Each chamber holds it's one lists of clusters 
304   Int_t                 fNclusters[kNCH];    //Array of current numbers of raw clusters
305
306   TClonesArray         *fCerenkovs;          //! ??? List of cerenkovs
307   Int_t                 fNcerenkovs;         //??? Current number of cerenkovs
308   TClonesArray         *fSpecials;           //! ??? List of specials
309   Int_t                 fNspecials;          //??? Current number of specials  
310   TObjArray            *fDchambers;          //! Array of lists of digits  
311   Int_t                 fNdch[kNCH];         //Array of current numbers of digits
312   TObjArray            *fRawClusters;        //! Array of lists of raw clusters 
313   Int_t                 fNrawch[kNCH];       //Array of current numbers of raw clusters
314   TObjArray            *fRecHits1D;          //!List of rec. hits
315   Int_t                 fNrechits1D[kNCH];   //Array of current numbers of rec hits 1D
316   TObjArray            *fRecHits3D;          //!List of rec. hits
317   Int_t                 fNrechits3D[kNCH];   //Array of current numbers of rec hits 3D 
318   Int_t fCkovNumber;                         // Number of Cerenkov photons
319   Int_t fFreonProd;                          // Cerenkovs produced in freon
320   Int_t fFeedbacks;                          // Number of feedback photons
321     
322   ClassDef(AliRICH,3)                        //Main RICH class 
323 };//class AliRICH  
324 //__________________________________________________________________________________________________
325 void AliRICH::AddHit(Int_t track, Int_t *vol, Float_t *hits)
326 {//Adds the current hit to the RICH hits list
327   TClonesArray &tmp=*fHits;
328   new(tmp[fNhits++])AliRICHhit(fIshunt,track,vol,hits);
329 }
330 //__________________________________________________________________________________________________
331 void AliRICH::AddSdigit(Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT0,Int_t iT1,Int_t iT2)
332 {//Adds the current Sdigit to the RICH list of Sdigits   
333   TClonesArray &tmp=*fSdigits;
334   new(tmp[fNsdigits++])AliRICHdigit(iC,iX,iY,iAdc,iT0,iT1,iT2);
335
336 //__________________________________________________________________________________________________
337 void AliRICH::AddDigit(Int_t iC,Int_t iX,Int_t iY,Int_t iAdc,Int_t iT0,Int_t iT1,Int_t iT2)
338 {//Adds the current digit to the corresponding RICH list of digits (individual list per chamber)
339   TClonesArray &tmp=*((TClonesArray*)fDigitsNew->At(iC-1));
340   new(tmp[fNdigitsNew[iC-1]++]) AliRICHdigit(iC,iX,iY,iAdc,iT0,iT1,iT2);
341 }
342 //__________________________________________________________________________________________________
343 void AliRICH::CreateHits()
344 {
345   if(fHits) return;
346   if(GetDebug())Info("CreateHits","creating hits container.");
347   fHits=new TClonesArray("AliRICHhit",10000);   fNhits=0;
348 }
349 //__________________________________________________________________________________________________
350 void AliRICH::CreateSdigits()
351 {
352   if(fSdigits) return;
353   if(GetDebug())Info("CreateSdigits","creating sdigits container.");
354   fSdigits=new TClonesArray("AliRICHdigit",10000); fNsdigits=0;
355 }
356 //__________________________________________________________________________________________________
357 void AliRICH::CreateDigits()
358 {
359   if(fDigitsNew) return;
360   if(GetDebug())Info("CreateDigits","creating digits containers.");
361   fDigitsNew = new TObjArray(kNCH);  
362   for(Int_t i=0;i<kNCH;i++) {fDigitsNew->AddAt(new TClonesArray("AliRICHdigit",10000), i); fNdigitsNew[i]=0;}
363 }
364 //__________________________________________________________________________________________________
365 void AliRICH::CreateClusters()
366 {
367   if(fClusters) return;
368   if(GetDebug())Info("CreateClusters","creating clusters containers.");
369   fClusters = new TObjArray(kNCH);  
370   for(Int_t i=0;i<kNCH;i++) {fClusters->AddAt(new TClonesArray("AliRICHcluster",10000), i); fNclusters[i]=0;}
371 }
372
373
374
375 //__________________________________________________________________________________________________
376 void AliRICH::CreateCerenkovsOld()
377 {
378   if(fCerenkovs) return;
379   if(GetDebug())Info("CreateCerenkovs","creating cerenkovs container.");
380   fCerenkovs=new TClonesArray("AliRICHCerenkov",10000);   fNcerenkovs=0;
381 }
382 //__________________________________________________________________________________________________
383 void AliRICH::CreateSpecialsOld()
384 {
385   if(fSpecials) return;
386   if(GetDebug())Info("CreateSpecialsOld","creating SDigits special container.");
387   fSpecials=new TClonesArray("AliRICHSDigit",100000); fNspecials=0;
388 }
389 //__________________________________________________________________________________________________
390 void AliRICH::CreateDigitsOld()
391 {
392   if(fDchambers) return;
393   if(GetDebug())Info("CreateDigitsOld","creating digits containers.");
394   fDchambers = new TObjArray(kNCH);  
395   for(Int_t i=0;i<kNCH;i++) fDchambers->AddAt(new TClonesArray("AliRICHDigit",10000), i); 
396 }
397 //__________________________________________________________________________________________________
398 void AliRICH::CreateRawClustersOld()
399 {
400   if(fRawClusters) return;
401   if(GetDebug())Info("CreateClustersOld","creating clusters containers.");
402   fRawClusters = new TObjArray(kNCH);
403   for(Int_t i=0; i<kNCH ;i++) fRawClusters->AddAt(new TClonesArray("AliRICHRawCluster",10000), i); 
404 }
405 //__________________________________________________________________________________________________
406 void AliRICH::CreateRecos1Old()
407 {
408   if(fRecHits1D) return;
409   if(GetDebug())Info("CreateRecos1DOld","creating recos 1 containers.");
410   fRecHits1D = new TObjArray(kNCH);
411   for(Int_t i=0; i<kNCH ;i++)  fRecHits1D->AddAt(new TClonesArray("AliRICHRecHit1D",1000), i);
412 }
413 //__________________________________________________________________________________________________
414 void AliRICH::CreateRecos3Old()
415 {
416   if(fRecHits3D) return;
417   if(GetDebug())Info("CreateRecos3DOld","creating recos 3 containers.");
418   fRecHits3D = new TObjArray(kNCH);
419   for(Int_t i=0; i<kNCH ;i++)  fRecHits3D->AddAt(new TClonesArray("AliRICHRecHit3D",1000), i);
420 }
421 //__________________________________________________________________________________________________
422 void AliRICH::AddCerenkov(Int_t track, Int_t *vol, Float_t *cerenkovs)
423 {//Adds the current RICH cerenkov hit to the Cerenkovs list   
424   TClonesArray &tmp=*fCerenkovs;
425   new(tmp[fNcerenkovs++]) AliRICHCerenkov(fIshunt,track,vol,cerenkovs);
426 }
427 //__________________________________________________________________________________________________
428 void AliRICH::AddSpecialOld(Int_t *aiSDigit)
429 {// Adds the current Sdigit to the RICH list of Specials
430   TClonesArray &lSDigits = *fSpecials;
431   new(lSDigits[fNspecials++]) AliRICHSDigit(aiSDigit);
432 }
433 //__________________________________________________________________________________________________
434 void AliRICH::AddDigitOld(Int_t iChamber, Int_t *tracks, Int_t *charges, Int_t *digits)
435 {// Adds the current digit to the RICH list of S digits
436    TClonesArray &ldigits = *((TClonesArray*)fDchambers->At(iChamber-1));
437    new(ldigits[fNdch[iChamber-1]++]) AliRICHDigit(tracks,charges,digits);
438 }
439 //__________________________________________________________________________________________________
440 void AliRICH::AddClusterOld(Int_t iChamber, const AliRICHRawCluster& c)
441 {// Add a RICH raw cluster to the list   
442   TClonesArray &tmp= *((TClonesArray*)fRawClusters->At(iChamber-1));
443   new(tmp[fNrawch[iChamber-1]++]) AliRICHRawCluster(c);
444 }
445 //__________________________________________________________________________________________________
446 #endif//#ifndef AliRICH_h