|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
#include "FWCore/ParameterSet/interface/ParameterSet.h"
//This class header
#include "FastSimulation/CalorimeterProperties/interface/ECALBarrelProperties.h"
#include "TMath.h"
#include <iostream>
#include "FWCore/MessageLogger/interface/MessageLogger.h"
ECALBarrelProperties::ECALBarrelProperties(const edm::ParameterSet& fastDet) {
edm::ParameterSet fastDetBarrel = fastDet.getParameter<edm::ParameterSet>("BarrelCalorimeterProperties");
lightColl_ = fastDetBarrel.getParameter<double>("lightColl");
lightCollUnif_ = fastDetBarrel.getParameter<double>("lightCollUnif");
photoStatistics_ = fastDetBarrel.getParameter<double>("photoStatistics");
thickness_ = fastDetBarrel.getParameter<double>("thickness");
interactionLength_ = fastDetBarrel.getParameter<double>("interactionLength");
Aeff_ = fastDetBarrel.getParameter<double>("Aeff");
Zeff_ = fastDetBarrel.getParameter<double>("Zeff");
rho_ = fastDetBarrel.getParameter<double>("rho");
radLenIngcm2_ = fastDetBarrel.getParameter<double>("radLenIngcm2");
radLenIncm_ = fastDetBarrel.getParameter<double>("radLenIncm");
radLenIncm_ = (radLenIncm_ < 0) ? radLenIngcm2_ / rho_ : radLenIncm_;
criticalEnergy_ = fastDetBarrel.getParameter<double>("criticalEnergy");
criticalEnergy_ =
(criticalEnergy_ < 0) ? 2.66E-3 * TMath::Power((radLenIngcm2_ * Zeff_ / Aeff_), 1.1) : criticalEnergy_;
moliereRadius_ = fastDetBarrel.getParameter<double>("moliereRadius");
moliereRadius_ = (moliereRadius_ < 0) ? scaleEnergy_ / criticalEnergy_ * radLenIncm_ : moliereRadius_;
Fs_ = fastDetBarrel.getParameter<double>("Fs");
ehat_ = fastDetBarrel.getParameter<double>("ehat");
resE_ = fastDetBarrel.getParameter<double>("resE");
da_ = fastDetBarrel.getParameter<double>("da");
dp_ = fastDetBarrel.getParameter<double>("dp");
bHom_ = fastDetBarrel.getParameter<bool>("bHom");
bool debug = fastDetBarrel.getParameter<bool>("debug");
if (debug)
edm::LogInfo("ECALProperties") << " ========== Barrel ========= "
<< "\n"
<< " isHom ? " << bHom_ << "\n"
<< " da = " << da_ << " dp = " << dp_ << " lightColl = " << lightColl_ << "\n"
<< " lightCollUnif_ = " << lightCollUnif_ << "\n"
<< " photoStatistics_ = " << photoStatistics_ << " photons/GeV\n"
<< " thickness_ = " << thickness_ << " cm\n"
<< " interactionLength_ = " << interactionLength_ << " cm\n"
<< " Aeff_ = " << Aeff_ << "\n"
<< " Zeff_ = " << Zeff_ << "\n"
<< " rho_ = " << rho_ << " g/cm3\n"
<< " radLenIngcm2_ = " << radLenIngcm2_ << " g*cm2\n"
<< " radLenIncm_ = " << radLenIncm_ << " cm\n"
<< " moliereRadius_ = " << moliereRadius_ << " cm\n"
<< " criticalEnergy_ = " << criticalEnergy_ << " GeV\n"
<< " scaleEnergy_ = " << scaleEnergy_ << " GeV\n"
<< " Fs = " << Fs_ << " ehat = " << ehat_ << " resE = " << resE_ << "\n";
}
|