/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
* See cxx source for full Copyright notice */
-#include <AliDetector.h>
-#include <TClonesArray.h>
+#include <AliDetector.h> //inheritance
+#include <TClonesArray.h>
#include <TObjArray.h>
#include <TVector.h>
#include <TVector3.h>
-#include "AliRICHDigitizer.h"
#include "AliRICHParam.h"
#include "AliRICHCluster.h"
#include "AliRICHHit.h"
//framework part
virtual Int_t IsVersion() const =0; //interface from
virtual void StepManager() =0; //interface from AliMC
- virtual void Hits2SDigits(); //interface from AliSimulation
- virtual void Digits2Raw(); //interface from AliSimulation
- virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* man) const {return new AliRICHDigitizer(man);} //interface from AliSimulation
virtual void SetTreeAddress(); //interface from AliLoader
virtual void MakeBranch(Option_t *opt=" "); //interface from AliLoader
virtual void CreateMaterials(); //interface from AliMC
inline void HitAdd (Int_t c,Int_t tid,TVector3 in,TVector3 out,Double_t e=0); //add new hit
inline void HitsCreate ( ); //create hits container
void HitsPrint (Int_t iEvent=0 )const; //prints hits
+ void HitsQA (Double_t cut=0,Double_t cutele=0,Double_t cutR=999);
TClonesArray* SDigits ( )const{return fSdigits;} //pointer to sdigits list
inline void SDigitAdd (Int_t c,TVector pad,Double_t q,Int_t pid,Int_t tid ) ; //add new sdigit
using AliDetector::Digits;
TClonesArray* Digits (Int_t iC )const{return fDigs ? (TClonesArray *)fDigs->At(iC-1):0;}
inline void DigitAdd (Int_t c,TVector pad,int q,int cfm,int *tid ) ; //add new digit
+ inline void DigitAdd (AliRICHDigit &dif ) ; //add new digit
inline void DigitsCreate ( ) ; //create digits
void DigitsReset ( ) {if(fDigs)for(int i=0;i<kNchambers;i++){fDigs->At(i)->Clear();fNdigs[i]=0;}} //virtual
void DigitsPrint (Int_t iEvent=0 )const; //prints digits
void ClustersReset ( ) {if(fClus)for(int i=0;i<kNchambers;i++){fClus ->At(i)->Clear();fNclus[i]=0;}}
void ClustersPrint (Int_t iEvent=0 )const; //prints a list of clusters for a given event
- void OccupancyPrint(Int_t iEvent=0 )const;
+ void OccupancyPrint(Int_t iEvent=-1 )const;
+ void SummaryOfEvent(Int_t iEvent=0 )const;
AliRICHChamber* C(Int_t iC) const{return fParam->C(iC);} //provides pointer to a given chamber
AliRICHParam* P() const{return fParam;} //provides pointer to a RICH params
void DisplayEvent(Int_t,Int_t)const; //Display event
static Int_t Nparticles(Int_t iPid,Int_t iEventN,AliRunLoader *pRunLoader); //counts total number of particle with iPid
void PrintTracks (Int_t iEvent=0); //prints a list of tracks for a given event
- void CheckPR ()const; //utility-> run staff for stack ??????
void ReadESD(Int_t iEventN, Int_t iChamber)const;
void DrawRing(TVector3 entrance,TVector3 vectorTrack,Double_t thetaCer)const;
for(Int_t i=0;i<kNchambers;i++) {fDigs->AddAt(new TClonesArray("AliRICHDigit",10000), i); fNdigs[i]=0;}
}
//__________________________________________________________________________________________________
+void AliRICH::DigitAdd(AliRICHDigit &dig)
+{
+//special for digit formed from raw
+ TClonesArray &tmp=*((TClonesArray*)fDigs->At(dig.Chamber()-1));
+ new(tmp[fNdigs[dig.Chamber()-1]++])AliRICHDigit(dig);
+}
+//__________________________________________________________________________________________________
void AliRICH::DigitAdd(int c,TVector pad,int q,int cfm,int *tid)
{
TClonesArray &tmp=*((TClonesArray*)fDigs->At(c-1));