|
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
|
#include "DD4hep/DetFactoryHelper.h"
#include "DataFormats/Math/interface/angle_units.h"
#include "DetectorDescription/DDCMS/interface/DDPlugins.h"
#include "DetectorDescription/DDCMS/interface/DDutils.h"
#include "FWCore/MessageLogger/interface/MessageLogger.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);
// Header section
std::string idNameSpace = static_cast<std::string>(ns.name()); //Namespace of this and ALL sub-parts
double eta = args.value<double>("Eta"); //Eta at which beam is focussed
double theta = 2.0 * atan(exp(-eta)); //Corresponding theta value
double shiftY = args.value<double>("ShiftX"); //Shift along Y
double shiftX = args.value<double>("ShiftY"); //Shift along X
double zoffset = args.value<double>("Zoffset"); //Offset in z
double dist = args.value<double>("Distance"); //Radial distance
double tilt = args.value<double>("TiltAngle"); //Tilt with respect to y-axis
int copyNumber = args.value<int>("Number"); //Copy Number
std::string childName = args.value<std::string>("ChildName"); //Children name
if (strchr(childName.c_str(), NAMESPACE_SEP) == nullptr)
childName = idNameSpace + childName;
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HCalGeom") << "DDHCalTBZposAlgo: Parameters for positioning-- Eta " << eta << "\tTheta "
<< convertRadToDeg(theta) << "\tShifts " << cms::convert2mm(shiftX) << ", "
<< cms::convert2mm(shiftY) << " along x, y axes; \tZoffest " << cms::convert2mm(zoffset)
<< "\tRadial Distance " << cms::convert2mm(dist) << "\tTilt angle "
<< convertRadToDeg(tilt) << "\tcopyNumber " << copyNumber;
edm::LogVerbatim("HCalGeom") << "DDHCalTBZposAlgo: Parent " << args.parentName() << "\tChild " << childName
<< " NameSpace " << idNameSpace;
#endif
dd4hep::Volume mother = ns.volume(args.parentName());
dd4hep::Volume child = ns.volume(childName);
// Create the rotation matrices
double thetax = 90._deg - theta;
double z = zoffset + dist * tan(thetax);
double x = shiftX - shiftY * sin(tilt);
double y = shiftY * cos(tilt);
dd4hep::Position tran(x, y, z);
dd4hep::Rotation3D rot;
if (tilt != 0) {
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HCalGeom") << "DDHCalZposAlgo: Creating a rotation \t90, " << convertRadToDeg(tilt) << ",90,"
<< (90 + convertRadToDeg(tilt)) << ", 0, 0";
#endif
rot = dd4hep::RotationZ(tilt);
}
mother.placeVolume(child, copyNumber, dd4hep::Transform3D(rot, tran));
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HCalGeom") << "DDHCalTBZposAlgo: " << child.name() << " number " << copyNumber << " positioned in "
<< mother.name() << " at (" << cms::convert2mm(x) << ", " << cms::convert2mm(y) << ", "
<< cms::convert2mm(z) << ") with " << rot;
#endif
return cms::s_executed;
}
// first argument is the type from the xml file
DECLARE_DDCMS_DETELEMENT(DDCMS_hcal_DDHCalTBZposAlgo, algorithm);
|