1 #ifndef ALIPHOSSIMPARAM_H
2 #define ALIPHOSSIMPARAM_H
3 /* Copyright(c) 2007, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
6 /* $Id: AliPHOSSimParam.h 23530 2008-01-25 06:46:13Z prsnko $ */
8 // Base class for the PHOS simulation parameters.
9 // Do not use in the reconstruction; use derivative classes instead.
14 class AliPHOSSimParam : public TNamed {
19 AliPHOSSimParam(const AliPHOSSimParam& recoParam);
20 AliPHOSSimParam& operator = (const AliPHOSSimParam& recoParam);
21 virtual ~AliPHOSSimParam() {}
23 static AliPHOSSimParam * GetInstance() ;
25 //Parameters used in conversion of deposited energy to APD response [see AliPHOSv1 for details]
26 void SetMeanLightYield(Float_t ly=47000.){fLightYieldMean=ly; //recalculate now dependencies
27 fLightFactor = fLightYieldMean * fIntrinsicAPDEfficiency ;
28 fAPDFactor = (13.418/fLightYieldMean/100.) * 300. ; }
29 void SetAPDEfficiency(Float_t eff=0.02655){fIntrinsicAPDEfficiency=eff;
30 fLightFactor = fLightYieldMean * fIntrinsicAPDEfficiency ;}
31 Float_t GetLightFactor(void) const { return fLightFactor ;}
32 Float_t GetAPDFactor(void) const { return fAPDFactor ;}
34 //Parameters use in EMC noise calculation [see AliPHOSDigitizer for details]
35 Float_t GetAPDNoise() const { return fAPDNoise; } //RMS of APD noise
36 void SetAPDNoise(Float_t noise=0.012){fAPDNoise = noise; }
38 //Parameters to apply non-lineary on cell level
39 Bool_t IsCellNonlinearityOn() const {return fCellNonLineaityOn;}
40 void SetCellNonLinearity(Bool_t on=kTRUE){fCellNonLineaityOn=on;} //default: on=kFALSE
41 Double_t GetCellNonLineairyA(void) const {return fCellNonLineaityA; }
42 Double_t GetCellNonLineairyB(void) const {return fCellNonLineaityB; }
43 void SetCellNonLineairyA(Double_t a=0.30) {fCellNonLineaityA = a; }
44 void SetCellNonLineairyB(Double_t b=0.109){fCellNonLineaityB = b; }
48 Float_t GetEmcDigitsThreshold() const { return fEMCDigitThreshold ; } //Minimal energy to keep digit
49 void SetEMCDigitsThreshold(Float_t thresh=0.01){fEMCDigitThreshold=thresh;}
51 //Parameters for energy digitization [see AliPHOSDigitizer for details]
52 void SetADCchannelW(Float_t width=0.005){fEMCADCchannel=width ;} //EMC channel width
53 Float_t GetADCchannelW(void) const {return fEMCADCchannel; }
55 void SetEDigitizationOn(Bool_t on=kTRUE){fDigitizeE=on ;} //Use digitization in simulation or left it
56 Bool_t IsEDigitizationOn(void)const {return fDigitizeE ;} //for Digits2Raw->Digits procedure
58 //Parameters for EMC TOF resolution [see AliPHOSDigitizer::TimeResolution()]
59 Float_t GetTOFa()const{return fTOFa ;} //constant term
60 Float_t GetTOFb()const{return fTOFb ;} //stohastic term
61 void SetTOFparameters(Float_t a=0.5e-9, Float_t b=1.5e-9){fTOFa=a; fTOFb=b; }
63 //Parameters for CPV noise and digitization [see AliPHOSDigitizer for details]
64 Float_t GetCPVNoise() const {return fCPVNoise ;} //RMS of CPV noise in
65 void SetCPVNoise(Float_t noise=0.01){ fCPVNoise = noise ;} //CPV popugais
67 Float_t GetCpvDigitsThreshold() const {return fCPVDigitThreshold ;} //Minimal energy to keep digit in
68 void SetCpvDigitsThreshold(Float_t thresh=0.09){fCPVDigitThreshold = thresh;} //CPV popugais
70 Float_t GetADCpedestalCpv() const {return fADCpedestalCpv ;} //CPV pedestal value
71 void SetADCpedestalCpv(Float_t ped=0.012){ fADCpedestalCpv=ped ;} //in CPV popugais
73 Float_t GetADCchanelCpv() const {return fADCchanelCpv;} //Price of one ADC channel
74 void SetADCchanelCpv(Float_t w=0.0012) {fADCchanelCpv=w; } //for CPV
76 Int_t GetNADCcpv() const {return fNADCcpv ;} //Max number of channels
77 void SettNADCcpv(Int_t n=12) { fNADCcpv=(Int_t)TMath::Power(2,n) ; } //in CPV ADC
79 //Mark streams for mixing as streams contaning Digits (true) or SDigits (false)
80 //Streamt numbering the same as in StreamManager
81 void SetStreamDigits(Int_t i){if(i<10)fDStream[i]=kTRUE;}
82 Bool_t IsStreamDigits(Int_t i){return fDStream[i]; }
84 //Parameters for RAW embedding
85 void SetEMCSubtractPedestals(Bool_t subtract) { fEMCSubtractPedestals = subtract;}
86 Bool_t EMCSubtractPedestals() const { return fEMCSubtractPedestals; }
88 void SetGlobalAltroOffset(Int_t offset) { fGlobalAltroOffset = offset ; }
89 Int_t GetGlobalAltroOffset() const { return fGlobalAltroOffset ; }
91 void SetGlobalAltroThreshold(Int_t ZSth) { fGlobalAltroThreshold = ZSth; }
92 Int_t GetGlobalAltroThreshold() const { return fGlobalAltroThreshold;}
94 void SetSampleQualityCut(Float_t qcut) { fEMCSampleQualityCut=qcut; }
95 Float_t GetEMCSampleQualityCut() const { return fEMCSampleQualityCut; }
99 AliPHOSSimParam(Int_t i); //True constructor which should be called by GetInstance()
103 //Parameters used in conversion of deposited energy to APD response (AliPHOSv1)
104 Float_t fLightYieldMean ; //Average number of photoelectrons per GeV
105 Float_t fIntrinsicAPDEfficiency; //APD efficiency including geometric coverage
106 Float_t fLightFactor ; //Average number of photons collected by APD per GeV deposited energy
107 Float_t fAPDFactor ; //factor relating light yield and APD response
109 //Parameters used in electronic noise calculation and thresholds for EMC (AliPHOSDigitizer)
110 Float_t fAPDNoise; //RMS of APD noise
111 Float_t fEMCDigitThreshold ; //minimal energy to keep digit
112 Float_t fEMCADCchannel ; //width of ADC channel in GeV
113 Float_t fTOFa ; //constant term of TOF resolution
114 Float_t fTOFb ; //stohastic term of TOF resolution
115 Float_t fCellNonLineaityA ; //Amp of cel non-linearity
116 Float_t fCellNonLineaityB ; //Energy scale of cel non-linearity
118 //Parameters used for RAW embedding
119 Bool_t fEMCSubtractPedestals; // true if pedestal should be subtracted (in non-ZS)
120 Int_t fGlobalAltroOffset ; // Offset used in ALTRO chips in SZ runs
121 Int_t fGlobalAltroThreshold ; // Threshold used in ALTRO chips in SZ runs
122 Float_t fEMCSampleQualityCut; // Cut on pulse shape fit quality
125 Float_t fADCpedestalCpv ; //Pedestal value
126 Float_t fADCchanelCpv ; //ADC channel width
127 Float_t fCPVNoise ; //RMS of CPV noise
128 Float_t fCPVDigitThreshold ; //Minimal energy to keep digit
129 Int_t fNADCcpv ; //Max number of channels in CPV ADC
131 Bool_t fDStream[10] ; //Mark mixing stream contains digits or SDigits
132 Bool_t fDigitizeE ; //Use energy digitization in simulation or left to Digits2Raw()
133 Bool_t fCellNonLineaityOn ; //Model scintillator non-linearity in AliPHOSDigitizer
135 static AliPHOSSimParam * fgSimParam ; // pointer to the unique instance of the class
137 ClassDef(AliPHOSSimParam,2)