|
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
|
#include "Geometry/CaloGeometry/interface/TruncatedPyramid.h"
#include "Geometry/EcalAlgo/interface/EcalEndcapGeometry.h"
#include "Geometry/CaloEventSetup/interface/CaloGeometryLoader.h"
#include "DetectorDescription/Core/interface/DDFilteredView.h"
typedef CaloGeometryLoader<EcalEndcapGeometry> EcalEGL;
template <>
void EcalEGL::fillGeom(EcalEndcapGeometry* geom,
const EcalEGL::ParmVec& vv,
const HepGeom::Transform3D& tr,
const DetId& id,
const double& scale);
template <>
void EcalEGL::fillNamedParams(const DDFilteredView& fv, EcalEndcapGeometry* geom);
template <>
void EcalEGL::fillNamedParams(const cms::DDFilteredView& fv, EcalEndcapGeometry* geom);
#include "Geometry/CaloEventSetup/interface/CaloGeometryLoader.icc"
template class CaloGeometryLoader<EcalEndcapGeometry>;
typedef CaloCellGeometry::CCGFloat CCGFloat;
template <>
void EcalEGL::fillGeom(EcalEndcapGeometry* geom,
const EcalEGL::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 EcalEGL::fillNamedParams(const DDFilteredView& _fv, EcalEndcapGeometry* geom) {
DDFilteredView fv = _fv;
bool doSubDets = fv.firstChild();
while (doSubDets) {
DDsvalues_type sv(fv.mergedSpecifics());
//ncrys
DDValue valNcrys("ncrys");
if (DDfetch(&sv, valNcrys)) {
const std::vector<double>& fvec = valNcrys.doubles();
// this parameter can only appear once
assert(fvec.size() == 1);
geom->setNumberOfCrystalPerModule(static_cast<int>(fvec[0]));
} else
continue;
//nmods
DDValue valNmods("nmods");
if (DDfetch(&sv, valNmods)) {
const std::vector<double>& fmvec = valNmods.doubles();
// there can only be one such value
assert(fmvec.size() == 1);
geom->setNumberOfModules(static_cast<int>(fmvec[0]));
}
break;
doSubDets = fv.nextSibling(); // go to next layer
}
}
template <>
void EcalEGL::fillNamedParams(const cms::DDFilteredView& fv, EcalEndcapGeometry* geom) {
const std::string specName = "ecal_ee";
//ncrys
std::vector<double> tempD = fv.get<std::vector<double> >(specName, "ncrys");
assert(tempD.size() == 1);
geom->setNumberOfCrystalPerModule(static_cast<int>(tempD[0]));
//nmods
tempD = fv.get<std::vector<double> >(specName, "nmods");
assert(tempD.size() == 1);
geom->setNumberOfModules(static_cast<int>(tempD[0]));
}
|