Digits::E measured in ADC counts now
[u/mrichter/AliRoot.git] / PHOS / AliPHOSSimParam.cxx
1 /**************************************************************************
2  * Copyright(c) 2007, ALICE Experiment at CERN, All rights reserved.      *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 // Base class for the PHOS simulation parameters.
17 // Do not use in the simulation; use derivative classes instead.
18 // Author: Dmitri Peressounko, RRC KI
19
20 // --- AliRoot header files ---
21 #include "AliPHOSSimParam.h"
22 #include "AliLog.h"
23
24 ClassImp(AliPHOSSimParam)
25
26 AliPHOSSimParam  * AliPHOSSimParam::fgSimParam = 0 ;
27 //-----------------------------------------------------------------------------
28 AliPHOSSimParam::AliPHOSSimParam() :
29   TNamed(),
30   fLightYieldMean(0.),fIntrinsicAPDEfficiency(0.),
31   fLightFactor(0.),fAPDFactor(0.),         
32   fAPDNoise(0.),fEMCDigitThreshold(0.),
33   fEMCADCchannel(0.),fTOFa(0.),fTOFb(0.),
34   fADCpedestalCpv(0.),fADCchanelCpv(0.),
35   fCPVNoise(0.),fCPVDigitThreshold(0.),fNADCcpv(0),
36   fDigitizeE(0)
37 {
38   //Default constructor.
39   for(Int_t i=0; i<10; i++) fDStream[i] = 0 ;
40 }
41
42 //-----------------------------------------------------------------------------
43 AliPHOSSimParam::AliPHOSSimParam(Int_t) :
44   TNamed(),
45   fLightYieldMean(0.),fIntrinsicAPDEfficiency(0.),
46   fLightFactor(0.),fAPDFactor(0.),         
47   fAPDNoise(0.),fEMCDigitThreshold(0.),
48   fEMCADCchannel(0.),fTOFa(0.),fTOFb(0.),
49   fADCpedestalCpv(0.),fADCchanelCpv(0.),
50   fCPVNoise(0.),fCPVDigitThreshold(0.),
51   fNADCcpv(0),
52   fDigitizeE(0)
53 {
54   //Real (private) constructor 
55   //Set default parameters
56
57   //Parameters describing energy deposition and light collection by APD, used in AliPHOSv1
58   //Photoelectron statistics:
59   // The light yield is a poissonian distribution of the number of
60   // photons created in the PbWo4 crystal, calculated using following formula
61   // NumberOfPhotons = EnergyLost * LightYieldMean* APDEfficiency 
62   // LightYieldMean is parameter calculated to be over 47000 photons per GeV
63   // APDEfficiency is 0.02655
64   // k_0 is 0.0045 from Valery Antonenko
65   // The number of electrons created in the APD is
66   // NumberOfElectrons = APDGain * LightYield
67   // The APD Gain is 300
68   fLightYieldMean = 47000;            //Average number of photoelectrons per GeV
69   fIntrinsicAPDEfficiency = 0.02655 ; //APD efficiency including geometric coverage
70 //  fLightYieldAttenuation  = 0.0045 ;  //light attenuation in PWO. Last analysis shows no z-position dependence
71 //                                      //so we removed this dependence from simulations 
72   fLightFactor            = fLightYieldMean * fIntrinsicAPDEfficiency ; //Average number of photons collected by 
73                             //APD per GeV deposited energy
74   fAPDFactor              = (13.418/fLightYieldMean/100.) * 300. ; //factor relating light yield and APD response
75                             //evaluated as (13.418/fLightYieldMean/100) * APDGain ;
76
77
78   //Parameters defining electronic noise calculation and Digits noise thresholds
79   //used in AliPHOSDigitizer
80   fAPDNoise           = 0.004 ;  // [GeV]
81   fEMCDigitThreshold  = 2.5   ;  // [ADC counts]
82   fEMCADCchannel      = 0.005 ;  // [GeV]
83   fTOFa               = 0.5e-9 ; // [sec] constant term
84   fTOFb               = 1.e-9 ;  // [sec/sqrt(GeV)]] stohastic term
85
86   fADCpedestalCpv     = 0.012 ;  // [aux units]
87   fADCchanelCpv       = 0.0012;  // [aux units]    
88   fCPVNoise           = 0.01;    // [aux units]
89   fCPVDigitThreshold  = 0.09 ;   // [aux units]
90   fNADCcpv  =  (Int_t)TMath::Power(2,12) ;
91
92   //Imput streams for merging. If true => this stream contains digits (and thus noise) and not SDigits.
93   for(Int_t i=0; i<10; i++){
94     fDStream[i] = 0 ;
95   }
96   fgSimParam = this ;
97 }
98
99 //-----------------------------------------------------------------------------
100 AliPHOSSimParam::AliPHOSSimParam(const AliPHOSSimParam& ):
101   TNamed(),
102   fLightYieldMean(0.),fIntrinsicAPDEfficiency(0.),
103   fLightFactor(0.),fAPDFactor(0.),         
104   fAPDNoise(0.),fEMCDigitThreshold(0.),
105   fEMCADCchannel(0.),fTOFa(0.),fTOFb(0.),
106   fADCpedestalCpv(0.),fADCchanelCpv(0.),
107   fCPVNoise(0.),fCPVDigitThreshold(0.),fNADCcpv(0),
108   fDigitizeE(0)
109 {
110   //Copy constructor.
111   AliError("Should not use copy constructor for singleton") ;
112   for(Int_t  i=0; i<10; i++){
113     fDStream[i] = 0 ;
114   }
115   fgSimParam = this ;
116 }
117 //-----------------------------------------------------------------------------                                                            
118 AliPHOSSimParam * AliPHOSSimParam::GetInstance(){
119
120   if(!fgSimParam)
121     new AliPHOSSimParam(0) ;
122   return fgSimParam ;
123 }
124 //-----------------------------------------------------------------------------
125 AliPHOSSimParam& AliPHOSSimParam::operator = (const AliPHOSSimParam& simParam)
126 {
127   //Assignment operator.
128
129   if(this != &simParam) {
130     AliError("Should not use operator= for singleton\n") ;
131   }
132
133   return *this;
134 }
135