Initialization of TVector3 data member is corrected
[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   fCellNonLineaityA(0.),fCellNonLineaityB(1.),fCellNonLineaityC(1.),
35   fEMCSubtractPedestals(kFALSE),
36   fGlobalAltroOffset(0),fGlobalAltroThreshold(0),fEMCSampleQualityCut(0),
37   fADCpedestalCpv(0.),fADCchanelCpv(0.),
38   fCPVNoise(0.),fCPVDigitThreshold(0.),fNADCcpv(0),
39   fDigitizeE(0),fCellNonLineaityOn(1)
40 {
41   //Default constructor.
42   for(Int_t i=0; i<10; i++) fDStream[i] = 0 ;
43 }
44
45 //-----------------------------------------------------------------------------
46 AliPHOSSimParam::AliPHOSSimParam(Int_t) :
47   TNamed(),
48   fLightYieldMean(0.),fIntrinsicAPDEfficiency(0.),
49   fLightFactor(0.),fAPDFactor(0.),         
50   fAPDNoise(0.),fEMCDigitThreshold(0.),
51   fEMCADCchannel(0.),fTOFa(0.),fTOFb(0.),
52   fCellNonLineaityA(0.),fCellNonLineaityB(1.),fCellNonLineaityC(1.),
53   fEMCSubtractPedestals(kFALSE),
54   fGlobalAltroOffset(0),fGlobalAltroThreshold(0),fEMCSampleQualityCut(0),
55   fADCpedestalCpv(0.),fADCchanelCpv(0.),
56   fCPVNoise(0.),fCPVDigitThreshold(0.),
57   fNADCcpv(0),
58   fDigitizeE(0),fCellNonLineaityOn(1)
59 {
60   //Real (private) constructor 
61   //Set default parameters
62
63   //Parameters describing energy deposition and light collection by APD, used in AliPHOSv1
64   //Photoelectron statistics:
65   // The light yield is a poissonian distribution of the number of
66   // photons created in the PbWo4 crystal, calculated using following formula
67   // NumberOfPhotons = EnergyLost * LightYieldMean* APDEfficiency 
68   // LightYieldMean is parameter calculated to be over 47000 photons per GeV
69   // APDEfficiency is 0.02655
70   // k_0 is 0.0045 from Valery Antonenko
71   // The number of electrons created in the APD is
72   // NumberOfElectrons = APDGain * LightYield
73   // The APD Gain is 300
74   fLightYieldMean = 47000;            //Average number of photoelectrons per GeV
75   fIntrinsicAPDEfficiency = 0.02655 ; //APD efficiency including geometric coverage
76 //  fLightYieldAttenuation  = 0.0045 ;  //light attenuation in PWO. Last analysis shows no z-position dependence
77 //                                      //so we removed this dependence from simulations 
78   fLightFactor            = fLightYieldMean * fIntrinsicAPDEfficiency ; //Average number of photons collected by 
79                             //APD per GeV deposited energy
80   fAPDFactor              = (13.418/fLightYieldMean/100.) * 300. ; //factor relating light yield and APD response
81                             //evaluated as (13.418/fLightYieldMean/100) * APDGain ;
82
83
84   //Parameters defining electronic noise calculation and Digits noise thresholds
85   //used in AliPHOSDigitizer
86   fAPDNoise           = 0.004 ;  // [GeV]
87   fEMCDigitThreshold  = 2.5   ;  // [ADC counts]
88   fEMCADCchannel      = 0.005 ;  // [GeV]
89   fTOFa               = 0.5e-9 ; // [sec] constant term
90   fTOFb               = 1.e-9 ;  // [sec/sqrt(GeV)]] stohastic term
91   fCellNonLineaityA   = 0.18 ;   //Amp of non-linearity of cell responce
92   fCellNonLineaityB   = 0.109;   //Scale of non-linearity of cell responce
93   fCellNonLineaityC   = 0.976;   //Overall calibration
94
95   fADCpedestalCpv     = 0.012 ;  // [aux units]
96   fADCchanelCpv       = 0.0012;  // [aux units]    
97   fCPVNoise           = 0.01;    // [aux units]
98   fCPVDigitThreshold  = 0.09 ;   // [aux units]
99   fNADCcpv  =  (Int_t)TMath::Power(2,12) ;
100
101   fGlobalAltroOffset = 10;
102   fGlobalAltroThreshold = 5;
103   fEMCSampleQualityCut = 4.;
104
105   //Imput streams for merging. If true => this stream contains digits (and thus noise) and not SDigits.
106   for(Int_t i=0; i<10; i++){
107     fDStream[i] = 0 ;
108   }
109   fgSimParam = this ;
110 }
111
112 //-----------------------------------------------------------------------------
113 AliPHOSSimParam::AliPHOSSimParam(const AliPHOSSimParam& ):
114   TNamed(),
115   fLightYieldMean(0.),fIntrinsicAPDEfficiency(0.),
116   fLightFactor(0.),fAPDFactor(0.),         
117   fAPDNoise(0.),fEMCDigitThreshold(0.),
118   fEMCADCchannel(0.),fTOFa(0.),fTOFb(0.),
119   fCellNonLineaityA(0.),fCellNonLineaityB(1.),fCellNonLineaityC(1.),
120   fEMCSubtractPedestals(kFALSE),
121   fGlobalAltroOffset(0),fGlobalAltroThreshold(0),fEMCSampleQualityCut(1.),
122   fADCpedestalCpv(0.),fADCchanelCpv(0.),
123   fCPVNoise(0.),fCPVDigitThreshold(0.),fNADCcpv(0),
124   fDigitizeE(0),fCellNonLineaityOn(1)
125 {
126   //Copy constructor.
127   AliError("Should not use copy constructor for singleton") ;
128   for(Int_t  i=0; i<10; i++){
129     fDStream[i] = 0 ;
130   }
131   fgSimParam = this ;
132 }
133 //-----------------------------------------------------------------------------                                                            
134 AliPHOSSimParam * AliPHOSSimParam::GetInstance(){
135
136   if(!fgSimParam)
137     new AliPHOSSimParam(0) ;
138   return fgSimParam ;
139 }
140 //-----------------------------------------------------------------------------
141 AliPHOSSimParam& AliPHOSSimParam::operator = (const AliPHOSSimParam& simParam)
142 {
143   //Assignment operator.
144
145   if(this != &simParam) {
146     AliError("Should not use operator= for singleton\n") ;
147   }
148
149   return *this;
150 }
151