|
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
|
#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"
using namespace angle_units::operators;
//#define EDM_ML_DEBUG
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 of original DDGEMAngular.h
float startAngle = args.value<float>("startAngle");
float stepAngle = args.value<float>("stepAngle");
int invert = args.value<int>("invert");
double rPos = args.value<double>("rPosition");
double zoffset = args.value<double>("zoffset");
int n = args.value<int>("n");
int startCopyNo = args.value<int>("startCopyNo");
int incrCopyNo = args.value<int>("incrCopyNo");
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "DDGEMAngular: Parameters for positioning-- " << n << " copies in steps of "
<< convertRadToDeg(stepAngle) << " from " << convertRadToDeg(startAngle)
<< " (inversion flag " << invert << ") \trPos " << cms::convert2mm(rPos) << " Zoffest "
<< cms::convert2mm(zoffset) << "\tStart and inremental "
<< "copy nos " << startCopyNo << ", " << incrCopyNo;
#endif
std::string childName = args.value<std::string>("ChildName");
childName = ns.prepend(childName);
std::string parentName = args.parentName();
parentName = ns.prepend(parentName);
dd4hep::Volume parent = ns.volume(parentName);
dd4hep::Volume child = ns.volume(childName);
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "DDGEMAngular: Parent " << parentName << "\tChild " << child.name() << "\tNameSpace "
<< ns.name();
#endif
// Now position child in mother *n* times
double phi = startAngle;
int copyNo = startCopyNo;
double thetax = 90.0_deg;
double thetay = invert == 0 ? 0.0 : 180.0_deg;
for (int ii = 0; ii < n; ii++) {
double phiz = phi;
if (phiz >= 2._pi)
phiz -= 2._pi;
double phix = invert == 0 ? (90.0_deg + phiz) : (-90.0_deg + phiz);
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "DDGEMAngular: Creating a rotation " << convertRadToDeg(thetax) << ", "
<< convertRadToDeg(phix) << ", " << convertRadToDeg(thetay) << ", 0, "
<< convertRadToDeg(thetax) << ", " << convertRadToDeg(phiz);
#endif
dd4hep::Rotation3D rotation = cms::makeRotation3D(thetax, phix, thetay, 0., thetax, phiz);
dd4hep::Position tran(rPos * cos(phiz), rPos * sin(phiz), zoffset);
parent.placeVolume(child, copyNo, dd4hep::Transform3D(rotation, tran));
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "DDGEMAngular: " << child.name() << " number " << copyNo << " positioned in "
<< parentName << " at (" << cms::convert2mm(rPos * cos(phiz)) << ","
<< cms::convert2mm(rPos * sin(phiz)) << "," << cms::convert2mm(zoffset) << ") with "
<< rotation;
#endif
phi += stepAngle;
copyNo += incrCopyNo;
}
return cms::s_executed;
}
// first argument is the type from the xml file
DECLARE_DDCMS_DETELEMENT(DDCMS_muon_DDGEMAngular, algorithm)
|