|
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
63
|
#include "FWCore/ParameterSet/interface/ParameterSet.h"
//This class header
#include "FastSimulation/CalorimeterProperties/interface/ECALEndcapProperties.h"
#include "TMath.h"
#include <iostream>
#include "FWCore/MessageLogger/interface/MessageLogger.h"
ECALEndcapProperties::ECALEndcapProperties(const edm::ParameterSet& fastDet) {
edm::ParameterSet fastDetEndcap = fastDet.getParameter<edm::ParameterSet>("EndcapCalorimeterProperties");
lightColl_ = fastDetEndcap.getParameter<double>("lightColl");
lightCollUnif_ = fastDetEndcap.getParameter<double>("lightCollUnif");
photoStatistics_ = fastDetEndcap.getParameter<double>("photoStatistics");
thickness_ = fastDetEndcap.getParameter<double>("thickness");
interactionLength_ = fastDetEndcap.getParameter<double>("interactionLength");
Aeff_ = fastDetEndcap.getParameter<double>("Aeff");
Zeff_ = fastDetEndcap.getParameter<double>("Zeff");
rho_ = fastDetEndcap.getParameter<double>("rho");
radLenIngcm2_ = fastDetEndcap.getParameter<double>("radLenIngcm2");
// Parameters that might be calculated out of the formulas
radLenIncm_ = fastDetEndcap.getParameter<double>("radLenIncm");
radLenIncm_ = (radLenIncm_ < 0) ? radLenIngcm2_ / rho_ : radLenIncm_;
criticalEnergy_ = fastDetEndcap.getParameter<double>("criticalEnergy");
criticalEnergy_ =
(criticalEnergy_ < 0) ? 2.66E-3 * TMath::Power((radLenIngcm2_ * Zeff_ / Aeff_), 1.1) : criticalEnergy_;
moliereRadius_ = fastDetEndcap.getParameter<double>("moliereRadius");
moliereRadius_ = (moliereRadius_ < 0) ? scaleEnergy_ / criticalEnergy_ * radLenIncm_ : moliereRadius_;
Fs_ = fastDetEndcap.getParameter<double>("Fs");
ehat_ = fastDetEndcap.getParameter<double>("ehat");
resE_ = fastDetEndcap.getParameter<double>("resE");
da_ = fastDetEndcap.getParameter<double>("da");
dp_ = fastDetEndcap.getParameter<double>("dp");
bHom_ = fastDetEndcap.getParameter<bool>("bHom");
bool debug = fastDetEndcap.getParameter<bool>("debug");
if (debug)
edm::LogInfo("ECALProperties") << " ========== Endcap ========= \n"
<< " isHom ? " << bHom_ << "\n"
<< " da = " << da_ << " dp = " << dp_ << " lightColl = " << lightColl_ << "\n"
<< " lightCollUnif_ = " << lightCollUnif_ << "\n"
<< " photoStatistics_ = " << photoStatistics_ << " photons/GeV\n"
<< " thickness_ = " << thickness_ << " in 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";
}
|