|
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
|
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "Geometry/EcalTestBeam/plugins/EcalTBHodoscopeGeometryLoaderFromDDD.h"
#include "Geometry/EcalTestBeam/interface/EcalTBHodoscopeGeometry.h"
#include "DataFormats/EcalDetId/interface/EBDetId.h"
#include "DetectorDescription/Core/interface/DDFilteredView.h"
#include "DetectorDescription/Core/interface/DDFilter.h"
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDSolid.h"
#include "DetectorDescription/Core/interface/DDTransform.h"
#include "DetectorDescription/Core/interface/DDExpandedView.h"
#include "DetectorDescription/Core/interface/DDName.h"
typedef CaloCellGeometry::CCGFloat CCGFloat;
#include <iostream>
#include <utility>
#include <vector>
#include <memory>
std::unique_ptr<CaloSubdetectorGeometry> EcalTBHodoscopeGeometryLoaderFromDDD::load(const DDCompactView* cpv) {
edm::LogVerbatim("EcalTBHodoscopeGeometry")
<< "[EcalTBHodoscopeGeometryLoaderFromDDD]:: start the construction of EcalTBHodoscope";
std::unique_ptr<CaloSubdetectorGeometry> ebg(new EcalTBHodoscopeGeometry());
makeGeometry(cpv, ebg.get());
edm::LogVerbatim("EcalTBHodoscopeGeometry")
<< "[EcalTBHodoscopeGeometryLoaderFromDDD]:: Returning EcalTBHodoscopeGeometry";
return ebg;
}
void EcalTBHodoscopeGeometryLoaderFromDDD::makeGeometry(const DDCompactView* cpv, CaloSubdetectorGeometry* ebg) {
if (ebg->cornersMgr() == nullptr)
ebg->allocateCorners(EBDetId::kSizeForDenseIndexing);
if (ebg->parMgr() == nullptr)
ebg->allocatePar(10, 3);
std::unique_ptr<DDFilter> filter{getDDFilter()};
DDFilteredView fv(*cpv, *filter);
bool doSubDets;
for (doSubDets = fv.firstChild(); doSubDets; doSubDets = fv.nextSibling()) {
#if 0
std::string answer = getDDDString("ReadOutName",&fv);
if (answer != "EcalTBH4BeamHits")
continue;
#endif
const DDSolid& solid = fv.logicalPart().solid();
if (solid.shape() != DDSolidShape::ddbox) {
throw cms::Exception("DDException")
<< std::string(__FILE__)
<< "\n CaloGeometryEcalTBHodoscope::upDate(...): currently only box fiber shapes supported ";
edm::LogWarning("EcalTBHodoscopeGeometry") << "Wrong shape for sensitive volume!" << solid;
}
std::vector<double> pv = solid.parameters();
// use preshower strip as box in space representation
// rotate the box and then move it
DD3Vector x, y, z;
fv.rotation().GetComponents(x, y, z);
CLHEP::Hep3Vector hx(x.X(), x.Y(), x.Z());
CLHEP::Hep3Vector hy(y.X(), y.Y(), y.Z());
CLHEP::Hep3Vector hz(z.X(), z.Y(), z.Z());
CLHEP::HepRotation hrot(hx, hy, hz);
CLHEP::Hep3Vector htran(fv.translation().X(), fv.translation().Y(), fv.translation().Z());
const HepGeom::Transform3D ht3d(hrot, // only scale translation
CaloCellGeometry::k_ScaleFromDDDtoGeant * htran);
const HepGeom::Point3D<float> ctr(ht3d * HepGeom::Point3D<float>(0, 0, 0));
const GlobalPoint refPoint(ctr.x(), ctr.y(), ctr.z());
std::vector<CCGFloat> vv;
vv.reserve(pv.size() + 1);
for (unsigned int i(0); i != pv.size(); ++i) {
vv.emplace_back(CaloCellGeometry::k_ScaleFromDDDtoGeant * pv[i]);
}
vv.emplace_back(0.); // tilt=0 here
const CCGFloat* pP(CaloCellGeometry::getParmPtr(vv, ebg->parMgr(), ebg->parVecVec()));
const DetId detId(getDetIdForDDDNode(fv));
//Adding cell to the Geometry
ebg->newCell(refPoint, refPoint, refPoint, pP, detId);
} // loop over all children
}
unsigned int EcalTBHodoscopeGeometryLoaderFromDDD::getDetIdForDDDNode(const DDFilteredView& fv) {
// perform some consistency checks
// get the parents and grandparents of this node
DDGeoHistory parents = fv.geoHistory();
assert(parents.size() >= 3);
EcalBaseNumber baseNumber;
//baseNumber.setSize(parents.size());
for (unsigned int i = 1; i <= parents.size(); i++) {
baseNumber.addLevel(parents[parents.size() - i].logicalPart().name().name(), parents[parents.size() - i].copyno());
}
return _scheme.getUnitID(baseNumber);
}
DDFilter* EcalTBHodoscopeGeometryLoaderFromDDD::getDDFilter() {
return new DDAndFilter<DDSpecificsMatchesValueFilter, DDSpecificsMatchesValueFilter>(
DDSpecificsMatchesValueFilter{DDValue("SensitiveDetector", "EcalTBH4BeamDetector", 0)},
DDSpecificsMatchesValueFilter{DDValue("ReadOutName", "EcalTBH4BeamHits", 0)});
}
|