1 /**************************************************************************
2 * Copyright(c) 2007, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 // Base class for the EMCAL simulation parameters.
23 // --- Root header files ---
25 // --- AliRoot header files ---
26 #include "AliEMCALSimParam.h"
30 ClassImp(AliEMCALSimParam)
32 AliEMCALSimParam * AliEMCALSimParam::fgSimParam = 0 ;
33 //-----------------------------------------------------------------------------
34 AliEMCALSimParam::AliEMCALSimParam() :
37 fMeanPhotonElectron(0),
42 fTimeResolutionPar0(0),
43 fTimeResolutionPar1(0),
44 fNADCEC(0),//Digitizer
47 fECPrimThreshold(0.) //SDigitizer
51 //Parameters in Digitizer
52 fMeanPhotonElectron = 4400; // electrons per GeV
53 fGainFluctuations = 15.; // obtained empiricaly to match beam test, from Paraskevi Ganoti
54 fPinNoise = 0.012; // APD noise in GeV from analysis test beam data
55 fDigitThreshold = 3; // 3 ADC counts not anymore cut in energy: //fPinNoise * 3; // 3 * sigma
56 fTimeNoise = 1.28e-5; // time noise in s
57 fTimeResolutionPar0 = 0.26666; // From F. Blanco: 0.51639^2
58 fTimeResolutionPar1 = 1.45861; // From F. Blanco: 1.20773^2
59 fTimeDelay = 600e-9 ; // 600 ns
61 fNADCEC = (Int_t) TMath::Power(2,16) ; // number of channels in Tower ADC - 65536
65 fB = 1.e+6; // Dynamic range now 2 TeV
66 fECPrimThreshold = 0.05; // GeV // threshold for deposit energy of hit
71 //-----------------------------------------------------------------------------
72 AliEMCALSimParam::AliEMCALSimParam(const AliEMCALSimParam& ):
75 fMeanPhotonElectron(0),
80 fTimeResolutionPar0(0),
81 fTimeResolutionPar1(0),
85 fECPrimThreshold(0.)//SDigitizer
88 AliError("Should not use copy constructor for singleton") ;
94 //-----------------------------------------------------------------------------
95 AliEMCALSimParam * AliEMCALSimParam::GetInstance(){
99 fgSimParam = new AliEMCALSimParam() ;
106 //-----------------------------------------------------------------------------
107 AliEMCALSimParam& AliEMCALSimParam::operator = (const AliEMCALSimParam& simParam)
109 //Assignment operator.
111 if(this != &simParam) {
112 AliError("Should not use operator= for singleton\n") ;
118 //-----------------------------------------------------------------------------
119 void AliEMCALSimParam::Print(Option_t *) const
121 // Print simulation parameters to stdout
123 printf("=== Parameters in Digitizer === \n");
124 printf("\t Electronics noise in EMC (fPinNoise) = %f, (fTimeNoise) = %e\n", fPinNoise, fTimeNoise) ;
125 printf("\t Threshold in EMC (fDigitThreshold) = %d\n", fDigitThreshold) ;
126 printf("\t Time Resolution (fTimeResolutionPar0) = %g\n", fTimeResolutionPar0) ;
127 printf("\t Time Resolution (fTimeResolutionPar1) = %g\n", fTimeResolutionPar1) ;
128 printf("\t Time Delay (fTimeDelay) = %g\n", fTimeDelay) ;
129 printf("\t Mean Photon-Electron (fMeanPhotonElectron) = %d, Gain Fluc. (fGainFluctuations) %2.1f\n", fMeanPhotonElectron,fGainFluctuations) ;
130 printf("\t N channels in EC section ADC (fNADCEC) = %d\n", fNADCEC) ;
134 printf("=== Parameters in SDigitizer === \n");
135 printf("\t sdigitization parameters A = %f\n", fA);
136 printf("\t B = %f\n", fB);
137 printf("\t Threshold for EC Primary assignment = %f\n", fECPrimThreshold);