|
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
|
#include "Geometry/CaloGeometry/interface/TruncatedPyramid.h"
#include "Geometry/EcalAlgo/interface/EcalBarrelGeometry.h"
#include "Geometry/CaloEventSetup/interface/CaloGeometryLoader.h"
#include "DetectorDescription/DDCMS/interface/DDFilteredView.h"
#include "DetectorDescription/Core/interface/DDFilteredView.h"
#include "DetectorDescription/Core/interface/DDutils.h"
typedef CaloGeometryLoader<EcalBarrelGeometry> EcalBGL;
template <>
void EcalBGL::fillGeom(EcalBarrelGeometry* geom,
const EcalBGL::ParmVec& vv,
const HepGeom::Transform3D& tr,
const DetId& id,
const double& scale);
template <>
void EcalBGL::fillNamedParams(const DDFilteredView& fv, EcalBarrelGeometry* geom);
template <>
void EcalBGL::fillNamedParams(const cms::DDFilteredView& fv, EcalBarrelGeometry* geom);
#include "Geometry/CaloEventSetup/interface/CaloGeometryLoader.icc"
template class CaloGeometryLoader<EcalBarrelGeometry>;
typedef CaloCellGeometry::CCGFloat CCGFloat;
template <>
void EcalBGL::fillGeom(EcalBarrelGeometry* geom,
const EcalBGL::ParmVec& vv,
const HepGeom::Transform3D& tr,
const DetId& id,
const double& scale) {
static constexpr uint32_t maxSize = 11;
std::vector<CCGFloat> pv;
unsigned int size = (vv.size() > maxSize) ? maxSize : vv.size();
unsigned int ioff = (vv.size() > maxSize) ? (vv.size() - maxSize) : 0;
pv.reserve(size);
for (unsigned int i(0); i != size; ++i) {
const CCGFloat factor(1 == i || 2 == i || 6 == i || 10 == i ? 1 : static_cast<CCGFloat>(scale));
pv.emplace_back(factor * vv[i + ioff]);
}
std::vector<GlobalPoint> corners(8);
TruncatedPyramid::createCorners(pv, tr, corners);
const CCGFloat* parmPtr(CaloCellGeometry::getParmPtr(pv, geom->parMgr(), geom->parVecVec()));
const GlobalPoint front(0.25 * (corners[0].x() + corners[1].x() + corners[2].x() + corners[3].x()),
0.25 * (corners[0].y() + corners[1].y() + corners[2].y() + corners[3].y()),
0.25 * (corners[0].z() + corners[1].z() + corners[2].z() + corners[3].z()));
const GlobalPoint back(0.25 * (corners[4].x() + corners[5].x() + corners[6].x() + corners[7].x()),
0.25 * (corners[4].y() + corners[5].y() + corners[6].y() + corners[7].y()),
0.25 * (corners[4].z() + corners[5].z() + corners[6].z() + corners[7].z()));
geom->newCell(front, back, corners[0], parmPtr, id);
}
template <>
void EcalBGL::fillNamedParams(const DDFilteredView& _fv, EcalBarrelGeometry* geom) {
DDFilteredView fv = _fv;
bool doSubDets = fv.firstChild();
while (doSubDets) {
DDsvalues_type sv(fv.mergedSpecifics());
//nxtalPhi
DDValue valnPhi("nxtalPhi");
if (DDfetch(&sv, valnPhi)) {
const std::vector<double>& fvec = valnPhi.doubles();
// this parameter can only appear once
assert(fvec.size() == 1);
geom->setNumXtalsPhiDirection(static_cast<int>(fvec[0]));
} else
continue;
DDValue valnEta("nxtalEta");
if (DDfetch(&sv, valnEta)) {
const std::vector<double>& fmvec = valnEta.doubles();
// there can only be one such value
assert(fmvec.size() == 1);
geom->setNumXtalsEtaDirection(static_cast<int>(fmvec[0]));
} else
// once we find nxtalPhi, the rest must also be defined
assert(1 == 0);
//EtaBaskets
DDValue valEtaB("EtaBaskets");
if (DDfetch(&sv, valEtaB)) {
const std::vector<double>& ebvec = valEtaB.doubles();
assert(!ebvec.empty());
std::vector<int> EtaBaskets;
EtaBaskets.reserve(ebvec.size());
for (const auto& ebv : ebvec)
EtaBaskets.emplace_back(static_cast<int>(ebv));
geom->setEtaBaskets(EtaBaskets);
} else
// once we find nxtalPhi, the rest must also be defined
assert(1 == 0);
//PhiBaskets
DDValue valPhi("PhiBaskets");
if (DDfetch(&sv, valPhi)) {
const std::vector<double>& pvec = valPhi.doubles();
assert(!pvec.empty());
geom->setBasketSizeInPhi(static_cast<int>(pvec[0]));
} else
// once we find nxtalPhi, the rest must also be defined
assert(1 == 0);
break;
doSubDets = fv.nextSibling(); // go to next layer
}
}
template <>
void EcalBGL::fillNamedParams(const cms::DDFilteredView& fv, EcalBarrelGeometry* geom) {
const std::string specName = "ecal_eb";
//nxtalEta
std::vector<double> tempD = fv.get<std::vector<double> >(specName, "nxtalEta");
assert(tempD.size() == 1);
geom->setNumXtalsEtaDirection(static_cast<int>(tempD[0]));
//nxtalPhi
tempD = fv.get<std::vector<double> >(specName, "nxtalPhi");
assert(tempD.size() == 1);
geom->setNumXtalsPhiDirection(static_cast<int>(tempD[0]));
//EtaBaskets
tempD = fv.get<std::vector<double> >(specName, "EtaBaskets");
assert(!tempD.empty());
geom->setEtaBaskets(dbl_to_int(tempD));
//PhiBaskets
tempD = fv.get<std::vector<double> >(specName, "PhiBaskets");
assert(!tempD.empty());
geom->setBasketSizeInPhi(static_cast<int>(tempD[0]));
}
|