|
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
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
|
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/EcalDetId/interface/EcalSubdetector.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
#include "Geometry/EcalAlgo/interface/EcalBarrelGeometry.h"
#include "Geometry/EcalAlgo/interface/EcalEndcapGeometry.h"
#include "Geometry/EcalAlgo/interface/EcalPreshowerGeometry.h"
#include "Geometry/CaloTopology/interface/CaloTopology.h"
#include "FastSimulation/CalorimeterProperties/interface/Calorimeter.h"
#include "FastSimulation/CalorimeterProperties/interface/PreshowerLayer1Properties.h"
#include "FastSimulation/CalorimeterProperties/interface/PreshowerLayer2Properties.h"
#include "FastSimulation/CalorimeterProperties/interface/ECALBarrelProperties.h"
#include "FastSimulation/CalorimeterProperties/interface/ECALEndcapProperties.h"
#include "FastSimulation/CalorimeterProperties/interface/HCALBarrelProperties.h"
#include "FastSimulation/CalorimeterProperties/interface/HCALEndcapProperties.h"
#include "FastSimulation/CalorimeterProperties/interface/HCALForwardProperties.h"
#include "Geometry/HcalTowerAlgo/interface/HcalHardcodeGeometryLoader.h"
Calorimeter::Calorimeter()
: myPreshowerLayer1Properties_(nullptr),
myPreshowerLayer2Properties_(nullptr),
myECALBarrelProperties_(nullptr),
myECALEndcapProperties_(nullptr),
myHCALBarrelProperties_(nullptr),
myHCALEndcapProperties_(nullptr),
myHCALForwardProperties_(nullptr),
EcalBarrelGeometry_(nullptr),
EcalEndcapGeometry_(nullptr),
HcalGeometry_(nullptr),
PreshowerGeometry_(nullptr) {
;
}
Calorimeter::Calorimeter(const edm::ParameterSet& fastCalo)
: myPreshowerLayer1Properties_(nullptr),
myPreshowerLayer2Properties_(nullptr),
myECALBarrelProperties_(nullptr),
myECALEndcapProperties_(nullptr),
myHCALBarrelProperties_(nullptr),
myHCALEndcapProperties_(nullptr),
myHCALForwardProperties_(nullptr),
EcalBarrelGeometry_(nullptr),
EcalEndcapGeometry_(nullptr),
HcalGeometry_(nullptr),
PreshowerGeometry_(nullptr) {
edm::ParameterSet fastDet = fastCalo.getParameter<edm::ParameterSet>("CalorimeterProperties");
edm::ParameterSet fastDetHF = fastCalo.getParameter<edm::ParameterSet>("ForwardCalorimeterProperties");
myPreshowerLayer1Properties_ = new PreshowerLayer1Properties(fastDet);
myPreshowerLayer2Properties_ = new PreshowerLayer2Properties(fastDet);
myECALBarrelProperties_ = new ECALBarrelProperties(fastDet);
myECALEndcapProperties_ = new ECALEndcapProperties(fastDet);
myHCALBarrelProperties_ = new HCALBarrelProperties(fastDet);
myHCALEndcapProperties_ = new HCALEndcapProperties(fastDet);
myHCALForwardProperties_ = new HCALForwardProperties(fastDetHF);
}
Calorimeter::~Calorimeter() {
if (myPreshowerLayer1Properties_)
delete myPreshowerLayer1Properties_;
if (myPreshowerLayer2Properties_)
delete myPreshowerLayer2Properties_;
if (myECALBarrelProperties_)
delete myECALBarrelProperties_;
if (myECALEndcapProperties_)
delete myECALEndcapProperties_;
if (myHCALBarrelProperties_)
delete myHCALBarrelProperties_;
if (myHCALEndcapProperties_)
delete myHCALEndcapProperties_;
if (myHCALForwardProperties_)
delete myHCALForwardProperties_;
}
const ECALProperties* Calorimeter::ecalProperties(int onEcal) const {
if (onEcal) {
if (onEcal == 1)
return myECALBarrelProperties_;
else
return myECALEndcapProperties_;
} else
return nullptr;
}
const HCALProperties* Calorimeter::hcalProperties(int onHcal) const {
if (onHcal) {
if (onHcal == 1)
return myHCALBarrelProperties_;
else if (onHcal == 2)
return myHCALEndcapProperties_;
else {
return myHCALForwardProperties_;
edm::LogInfo("CalorimeterProperties")
<< " Calorimeter::hcalProperties : set myHCALForwardProperties" << std::endl;
}
} else
return nullptr;
}
const PreshowerLayer1Properties* Calorimeter::layer1Properties(int onLayer1) const {
if (onLayer1)
return myPreshowerLayer1Properties_;
else
return nullptr;
}
const PreshowerLayer2Properties* Calorimeter::layer2Properties(int onLayer2) const {
if (onLayer2)
return myPreshowerLayer2Properties_;
else
return nullptr;
}
void Calorimeter::setupGeometry(const CaloGeometry& pG) {
edm::LogInfo("CalorimeterProperties") << " setupGeometry " << std::endl;
EcalBarrelGeometry_ = dynamic_cast<const EcalBarrelGeometry*>(pG.getSubdetectorGeometry(DetId::Ecal, EcalBarrel));
EcalEndcapGeometry_ = dynamic_cast<const EcalEndcapGeometry*>(pG.getSubdetectorGeometry(DetId::Ecal, EcalEndcap));
HcalGeometry_ = pG.getSubdetectorGeometry(DetId::Hcal, HcalBarrel);
// Takes a lot of time
PreshowerGeometry_ =
dynamic_cast<const EcalPreshowerGeometry*>(pG.getSubdetectorGeometry(DetId::Ecal, EcalPreshower));
}
void Calorimeter::setupTopology(const CaloTopology& theTopology) {
EcalBarrelTopology_ = theTopology.getSubdetectorTopology(DetId::Ecal, EcalBarrel);
EcalEndcapTopology_ = theTopology.getSubdetectorTopology(DetId::Ecal, EcalEndcap);
}
const CaloSubdetectorGeometry* Calorimeter::getEcalGeometry(int subdetn) const {
if (subdetn == 1)
return EcalBarrelGeometry_;
if (subdetn == 2)
return EcalEndcapGeometry_;
if (subdetn == 3)
return PreshowerGeometry_;
edm::LogWarning("Calorimeter") << "Requested an invalid ECAL subdetector geometry: " << subdetn << std::endl;
return nullptr;
}
const CaloSubdetectorTopology* Calorimeter::getEcalTopology(int subdetn) const {
if (subdetn == 1)
return EcalBarrelTopology_;
if (subdetn == 2)
return EcalEndcapTopology_;
edm::LogWarning("Calorimeter") << "Requested an invalid ECAL subdetector topology: " << subdetn << std::endl;
return nullptr;
}
|