]> git.uio.no Git - u/mrichter/AliRoot.git/blame - RICH/AliRICH.h
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
More tools for QA
[u/mrichter/AliRoot.git] / RICH / AliRICH.h
CommitLineData
2f614988 1#ifndef AliRICH_h
2#define AliRICH_h
3da30618 3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
5
0fe8fa07 6#include <AliDetector.h>
853634d3 7#include <TClonesArray.h>
1c349ae5 8#include <TObjArray.h>
9#include <TVector.h>
d4c94996 10#include <TVector3.h>
1c349ae5 11
3582c1f9 12#include "AliRICHDigitizer.h"
cb8189e7 13#include "AliRICHParam.h"
0fe8fa07 14#include "AliRICHCluster.h"
15#include "AliRICHHit.h"
d4c94996 16
c021cb15 17//__________________AliRICH_________________________________________________________________________
e42a7b46 18class AliESD;
3ea9cb08 19
2f614988 20class AliRICH : public AliDetector
21{
2f614988 22public:
a25b3368 23//ctor & dtor
853634d3 24 AliRICH();
af3d25a6 25 AliRICH(const char *name, const char *title);
d4c94996 26 AliRICH(const AliRICH& RICH):AliDetector(RICH) {;} //copy ctor
853634d3 27 virtual ~AliRICH();
dfb4e77d 28
942194a4 29 AliRICH& operator=(const AliRICH&) {return *this;}
d4c94996 30//framework part
e42a7b46 31 virtual Int_t IsVersion() const =0; //interface from
32 virtual void StepManager() =0; //interface from AliMC
33 virtual void Hits2SDigits(); //interface from AliSimulation
d3eb6079 34 virtual void Digits2Raw(); //interface from AliSimulation
e42a7b46 35 virtual AliDigitizer* CreateDigitizer(AliRunDigitizer* man) const {return new AliRICHDigitizer(man);} //interface from AliSimulation
e42a7b46 36 virtual void SetTreeAddress(); //interface from AliLoader
37 virtual void MakeBranch(Option_t *opt=" "); //interface from AliLoader
38 virtual void CreateMaterials(); //interface from AliMC
39 virtual void CreateGeometry(); //interface from AliMC
0fe8fa07 40 virtual void BuildGeometry(); //interface
41 virtual void Print(Option_t *option="") const; //prints current RICH status
42//private part
d3eb6079 43 void GeomPadPanelFrame()const; //defines PPF geometry
44 void GeomAmpGap() const; //defines gap geometry + anod wires
45 void GeomRadiators() const; //defines radiators geometry
46 void GeomSandBox() const; //defines sandbox geometry
47 void GeomRadioSrc() const; //defines radio source geometry
48 void GeomAerogel() const; //defines aerogel geometry
a25b3368 49 static Float_t Fresnel(Float_t ene,Float_t pdoti, Bool_t pola); //deals with Fresnel absorption
0fe8fa07 50
51 AliRICHHit* Hit (Int_t tid )const; //first hit of given TID
52 inline void HitAdd (Int_t c,Int_t tid,TVector3 in,TVector3 out,Double_t e=0); //add new hit
53 inline void HitsCreate ( ); //create hits container
54 void HitsPrint (Int_t iEvent=0 )const; //prints hits
e53ffca5 55 void HitsQA (Double_t cut=0,Double_t cutele=0,Double_t cutR=999);
0fe8fa07 56
57 TClonesArray* SDigits ( )const{return fSdigits;} //pointer to sdigits list
58 inline void SDigitAdd (Int_t c,TVector pad,Double_t q,Int_t pid,Int_t tid ) ; //add new sdigit
59 inline void SDigitsCreate ( ) ; //create sdigits container
60 void SDigitsReset ( ) {fNsdigits=0; if(fSdigits) fSdigits ->Clear();} //clean a list of sdigits
61 void SDigitsPrint (Int_t iEvent=0 )const; //prints sdigits
62
63 using AliDetector::Digits;
64 TClonesArray* Digits (Int_t iC )const{return fDigs ? (TClonesArray *)fDigs->At(iC-1):0;}
65 inline void DigitAdd (Int_t c,TVector pad,int q,int cfm,int *tid ) ; //add new digit
66 inline void DigitsCreate ( ) ; //create digits
67 void DigitsReset ( ) {if(fDigs)for(int i=0;i<kNchambers;i++){fDigs->At(i)->Clear();fNdigs[i]=0;}} //virtual
68 void DigitsPrint (Int_t iEvent=0 )const; //prints digits
69
70 TClonesArray* Clusters (Int_t iC )const{if(fClus) return (TClonesArray *)fClus->At(iC-1);else return 0;}
71 inline void ClusterAdd (AliRICHCluster &cl ) ; //add new cluster
72 inline void ClustersCreate( ) ; //create clusters container
73 void ClustersReset ( ) {if(fClus)for(int i=0;i<kNchambers;i++){fClus ->At(i)->Clear();fNclus[i]=0;}}
74 void ClustersPrint (Int_t iEvent=0 )const; //prints a list of clusters for a given event
c23b3bc0 75
e53ffca5 76 void OccupancyPrint(Int_t iEvent=-1 )const;
77 void SummaryOfEvent(Int_t iEvent=0 )const;
0fe8fa07 78
79 AliRICHChamber* C(Int_t iC) const{return fParam->C(iC);} //provides pointer to a given chamber
80 AliRICHParam* P() const{return fParam;} //provides pointer to a RICH params
a25b3368 81 AliRICH* R() {return this;} //provides pointer to RICH main object
82 TVector Counters() const{return fCounters;} //provides a set of counters
83 void ControlPlots(); //creates ~/RCP.root with a set of QA plots
0fe8fa07 84 void Display()const; //Display event
85 void DisplayEvent(Int_t,Int_t)const; //Display event
86 static Int_t Nparticles(Int_t iPid,Int_t iEventN,AliRunLoader *pRunLoader); //counts total number of particle with iPid
a25b3368 87 void PrintTracks (Int_t iEvent=0); //prints a list of tracks for a given event
d3eb6079 88 void CheckPR ()const; //utility-> run staff for stack ??????
56030c8d 89 void ReadESD(Int_t iEventN, Int_t iChamber)const;
90 void DrawRing(TVector3 entrance,TVector3 vectorTrack,Double_t thetaCer)const;
91
c60862bf 92protected:
d3eb6079 93 enum EMedia {kAir=1,kRoha,kSiO2,kC6F14,kCH4,kCsI,kGridCu,kOpSiO2,kGap,kAl,kGlass,kCu,kW,kSteel,kPerpex,kSr90,kMylar,kGel,kReflector};
a25b3368 94 enum ECounters {kStepManager=0,kCerProdTot,kCerProdRad,kCerKillTot,kCerKillRad,kCerKillRef,kEleProdTot};
0fe8fa07 95 AliRICHParam *fParam; //main RICH parametrization
a25b3368 96 //fHits and fDigits belong to AliDetector
97 TClonesArray *fSdigits; //! list of sdigits
98 Int_t fNsdigits; //! current number of sdigits
84365c70 99
0fe8fa07 100 TObjArray *fDigs; //! each chamber holds it's one lists of digits
101 Int_t fNdigs[7]; //! array of current numbers of digits
84365c70 102
0fe8fa07 103 TObjArray *fClus; //! each chamber holds it's one lists of clusters
104 Int_t fNclus[7]; //! array of current numbers of raw clusters
84365c70 105
a25b3368 106 TVector fCounters; //Particle history counters, explanation in StepManager()
107
0fe8fa07 108 ClassDef(AliRICH,9) //Main RICH class
c60862bf 109};//class AliRICH
ed3ceb24 110
c1863e3c 111//__________________________________________________________________________________________________
0fe8fa07 112void AliRICH::HitsCreate()
c60862bf 113{
114 if(fHits) return;
998b831f 115 AliDebug(1,"creating hits container.");
0fe8fa07 116 fHits=new TClonesArray("AliRICHHit",10000); fNhits=0;
543d5224 117}
c60862bf 118//__________________________________________________________________________________________________
0fe8fa07 119void AliRICH::HitAdd(Int_t c,Int_t tid,TVector3 i3,TVector3 o3,Double_t eloss)
a25b3368 120{
121//add new RICH hit to the list of hits
122 TClonesArray &tmp=*fHits;
0fe8fa07 123 new(tmp[fNhits++])AliRICHHit(c,tid,i3,o3,eloss);
a25b3368 124}//AddHit()
125//__________________________________________________________________________________________________
0fe8fa07 126void AliRICH::SDigitsCreate()
127{
128 if(fSdigits) return;
129 AliDebug(1,"creating sdigits container.");
130 fSdigits=new TClonesArray("AliRICHDigit",10000); fNsdigits=0;
131}
132//__________________________________________________________________________________________________
133void AliRICH::SDigitAdd(Int_t c,TVector pad,Double_t q,Int_t pid,Int_t tid)
e42a7b46 134{
135 Int_t cfm;
ed3ceb24 136 switch(pid){
e42a7b46 137 case 50000050: cfm=1000000;break;//cerenkov
138 case 50000051: cfm=1000; break;//feedback
139 default: cfm=1; break;//mip
ed3ceb24 140 }
0fe8fa07 141 TClonesArray &tmp=*fSdigits; new(tmp[fNsdigits++])AliRICHDigit(c,pad,q,cfm,tid,-1,-1);
e42a7b46 142}
d4c94996 143//__________________________________________________________________________________________________
0fe8fa07 144void AliRICH::DigitsCreate()
a25b3368 145{
0fe8fa07 146 if(fDigs) return;
147 AliDebug(1,"creating digits containers.");
148 fDigs = new TObjArray(kNchambers);
149 for(Int_t i=0;i<kNchambers;i++) {fDigs->AddAt(new TClonesArray("AliRICHDigit",10000), i); fNdigs[i]=0;}
150}
151//__________________________________________________________________________________________________
152void AliRICH::DigitAdd(int c,TVector pad,int q,int cfm,int *tid)
153{
154 TClonesArray &tmp=*((TClonesArray*)fDigs->At(c-1));
155 new(tmp[fNdigs[c-1]++])AliRICHDigit(c,pad,q,cfm,tid[0],tid[1],tid[2]);
a25b3368 156}
157//__________________________________________________________________________________________________
0fe8fa07 158void AliRICH::ClustersCreate()
159{
160 if(fClus) return;
161 AliDebug(1,"creating clusters containers.");
162 fClus = new TObjArray(kNchambers);
163 for(Int_t i=0;i<kNchambers;i++) {fClus->AddAt(new TClonesArray("AliRICHCluster",10000), i); fNclus[i]=0;}
164}
165//__________________________________________________________________________________________________
166void AliRICH::ClusterAdd(AliRICHCluster &cl)
a25b3368 167{
0fe8fa07 168 Int_t c=cl.C()-1;TClonesArray &tmp=*((TClonesArray*)fClus->At(c));
169 new(tmp[fNclus[c]++])AliRICHCluster(cl);
a25b3368 170}
0fe8fa07 171//__________________________________________________________________________________________________
c60862bf 172#endif//#ifndef AliRICH_h