|
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
|
#include "DataFormats/Math/interface/angle_units.h"
#include "DD4hep/DetFactoryHelper.h"
#include "DetectorDescription/DDCMS/interface/DDPlugins.h"
#include "DetectorDescription/DDCMS/interface/DDutils.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"
//#define EDM_ML_DEBUG
using namespace angle_units::operators;
static long algorithm(dd4hep::Detector& /* description */, cms::DDParsingContext& ctxt, xml_h e) {
cms::DDNamespace ns(ctxt, e, true);
cms::DDAlgoArguments args(ctxt, e);
std::string motherName = args.parentName();
auto const& m_tiltAngle = args.value<double>("tiltAngle"); // Tilt angle
auto const& m_rMin = args.value<double>("rMin"); // Inner radius
auto const& m_rMax = args.value<double>("rMax"); // Outer radius
auto const& m_zoffset = args.value<double>("zoffset"); // Offset in z
auto const& m_xyoffset = args.value<double>("xyoffset"); // Offset in x or y
auto const& m_startCopyNo = args.value<int>("startCopyNo"); // Start copy Number
auto const& m_incrCopyNo = args.value<int>("incrCopyNo"); // Increment copy Number
auto const& m_childName = args.value<std::string>("ChildName"); // Children name
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HGCalGeom") << "Tilt Angle " << m_tiltAngle << " R " << cms::convert2mm(m_rMin) << ":"
<< cms::convert2mm(m_rMax) << " Offset " << cms::convert2mm(m_zoffset) << ":"
<< cms::convert2mm(m_xyoffset) << " Copy " << m_startCopyNo << ":" << m_incrCopyNo
<< " Child " << m_childName;
edm::LogVerbatim("HGCalGeom") << "DDHGCalNoTaperEndcap: NameSpace " << ns.name() << "\tParent " << args.parentName();
#endif
dd4hep::Volume parent = ns.volume(args.parentName());
std::string name = ns.prepend(m_childName);
const int ix[4] = {1, -1, -1, 1};
const int iy[4] = {1, 1, -1, -1};
int copyNo = m_startCopyNo;
for (int i = 0; i < 4; ++i) {
int xQuadrant = ix[i];
int yQuadrant = iy[i];
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HGCalGeom") << "Create quarter " << xQuadrant << ":" << yQuadrant;
#endif
double tiltAngle = m_tiltAngle;
double xphi = xQuadrant * tiltAngle;
double yphi = yQuadrant * tiltAngle;
double theta = 90._deg;
double phiX = 0.0;
double phiY = theta;
double phiZ = 3 * theta;
double offsetZ = m_zoffset;
double offsetXY = m_xyoffset;
double offsetX = xQuadrant * 0.5 * offsetXY;
double offsetY = yQuadrant * 0.5 * offsetXY;
#ifdef EDM_ML_DEBUG
int rowmax(0), column(0);
#endif
while (std::abs(offsetX) < m_rMax) {
#ifdef EDM_ML_DEBUG
column++;
int row(0);
#endif
while (std::abs(offsetY) < m_rMax) {
#ifdef EDM_ML_DEBUG
row++;
#endif
double limit1 = sqrt((offsetX + 0.5 * xQuadrant * offsetXY) * (offsetX + 0.5 * xQuadrant * offsetXY) +
(offsetY + 0.5 * yQuadrant * offsetXY) * (offsetY + 0.5 * yQuadrant * offsetXY));
double limit2 = sqrt((offsetX - 0.5 * xQuadrant * offsetXY) * (offsetX - 0.5 * xQuadrant * offsetXY) +
(offsetY - 0.5 * yQuadrant * offsetXY) * (offsetY - 0.5 * yQuadrant * offsetXY));
// Make sure we do not add supermodules in rMin area
if (limit2 > m_rMin && limit1 < m_rMax) {
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HGCalGeom") << m_childName << " copyNo = " << copyNo << " (" << column << "," << row
<< "): offsetX,Y = " << cms::convert2mm(offsetX) << ","
<< cms::convert2mm(offsetY) << " limit=" << cms::convert2mm(limit1) << ":"
<< cms::convert2mm(limit2) << " rMin, rMax = " << cms::convert2mm(m_rMin) << ","
<< cms::convert2mm(m_rMax);
#endif
dd4hep::Rotation3D rotation = (cms::makeRotation3D(theta, phiX, theta + yphi, phiY, -yphi, phiZ) *
cms::makeRotation3D(theta + xphi, phiX, 90._deg, 90._deg, xphi, 0.0));
dd4hep::Position tran(offsetX, offsetY, offsetZ);
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HGCalGeom") << "Module " << copyNo << ": location = (" << cms::convert2mm(offsetX) << ","
<< cms::convert2mm(offsetY) << "," << cms::convert2mm(offsetZ) << ") Rotation "
<< rotation;
#endif
parent.placeVolume(ns.volume(name), copyNo, dd4hep::Transform3D(rotation, tran));
copyNo += m_incrCopyNo;
} else {
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HGCalGeom") << " (" << column << "," << row << "): offsetX,Y = " << cms::convert2mm(offsetX)
<< "," << cms::convert2mm(offsetY)
<< " is out of limit=" << cms::convert2mm(limit1) << ":"
<< cms::convert2mm(limit2) << " rMin, rMax = " << cms::convert2mm(m_rMin) << ","
<< cms::convert2mm(m_rMax);
#endif
}
yphi += yQuadrant * 2. * tiltAngle;
offsetY += yQuadrant * offsetXY;
}
#ifdef EDM_ML_DEBUG
if (row > rowmax)
rowmax = row;
#endif
xphi += xQuadrant * 2. * tiltAngle;
yphi = yQuadrant * tiltAngle;
offsetY = yQuadrant * 0.5 * offsetXY;
offsetX += xQuadrant * offsetXY;
}
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HGCalGeom") << rowmax << " rows and " << column << " columns in quadrant " << xQuadrant << ":"
<< yQuadrant;
#endif
}
return cms::s_executed;
}
// first argument is the type from the xml file
DECLARE_DDCMS_DETELEMENT(DDCMS_hgcal_DDHGCalNoTaperEndcap, algorithm)
|