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),
41 fTimeResolutionPar0(0),
42 fTimeResolutionPar1(0),
43 fNADCEC(0),//Digitizer
46 fECPrimThreshold(0.) //SDigitizer
50 //Parameters in Digitizer
51 fMeanPhotonElectron = 4400; // electrons per GeV
52 fPinNoise = 0.012; // APD noise in GeV from analysis test beam data
53 fDigitThreshold = 3; // 3 ADC counts not anymore cut in energy: //fPinNoise * 3; // 3 * sigma
54 fTimeNoise = 1.28e-5; // time noise in s
55 fTimeResolutionPar0 = 0.26666; // From F. Blanco: 0.51639^2
56 fTimeResolutionPar1 = 1.45861; // From F. Blanco: 1.20773^2
57 fTimeDelay = 600e-9 ; // 600 ns
59 fNADCEC = (Int_t) TMath::Power(2,16) ; // number of channels in Tower ADC - 65536
63 fB = 1.e+6; // Dynamic range now 2 TeV
64 fECPrimThreshold = 0.05; // GeV // threshold for deposit energy of hit
69 //-----------------------------------------------------------------------------
70 AliEMCALSimParam::AliEMCALSimParam(const AliEMCALSimParam& ):
73 fMeanPhotonElectron(0),
77 fTimeResolutionPar0(0),
78 fTimeResolutionPar1(0),
82 fECPrimThreshold(0.)//SDigitizer
85 AliError("Should not use copy constructor for singleton") ;
91 //-----------------------------------------------------------------------------
92 AliEMCALSimParam * AliEMCALSimParam::GetInstance(){
96 fgSimParam = new AliEMCALSimParam() ;
103 //-----------------------------------------------------------------------------
104 AliEMCALSimParam& AliEMCALSimParam::operator = (const AliEMCALSimParam& simParam)
106 //Assignment operator.
108 if(this != &simParam) {
109 AliError("Should not use operator= for singleton\n") ;
115 //-----------------------------------------------------------------------------
116 void AliEMCALSimParam::Print(Option_t *) const
118 // Print simulation parameters to stdout
120 printf("=== Parameters in Digitizer === \n");
121 printf("\t Electronics noise in EMC (fPinNoise) = %f, (fTimeNoise) = %e\n", fPinNoise, fTimeNoise) ;
122 printf("\t Threshold in EMC (fDigitThreshold) = %d\n", fDigitThreshold) ;
123 printf("\t Time Resolution (fTimeResolutionPar0) = %g\n", fTimeResolutionPar0) ;
124 printf("\t Time Resolution (fTimeResolutionPar1) = %g\n", fTimeResolutionPar1) ;
125 printf("\t Time Delay (fTimeDelay) = %g\n", fTimeDelay) ;
126 printf("\t Mean Photon-Electron (fMeanPhotonElectron) = %d\n", fMeanPhotonElectron) ;
127 printf("\t N channels in EC section ADC (fNADCEC) = %d\n", fNADCEC) ;
131 printf("=== Parameters in SDigitizer === \n");
132 printf("\t sdigitization parameters A = %f\n", fA);
133 printf("\t B = %f\n", fB);
134 printf("\t Threshold for EC Primary assignment = %f\n", fECPrimThreshold);