|
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
|
#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& context, xml_h element) {
cms::DDNamespace ns(context, element, true);
cms::DDAlgoArguments args(context, element);
int n = args.value<int>("n");
int startCopyNo = args.find("startCopyNo") ? args.value<int>("startCopyNo") : 1;
int incrCopyNo = args.find("incrCopyNo") ? args.value<int>("incrCopyNo") : 1;
float startAngle = args.value<float>("startAngle");
float stepAngle = args.value<float>("stepAngle");
float zoffset = args.value<float>("zoffset");
std::string rotns = args.value<std::string>("RotNameSpace");
dd4hep::Volume mother = ns.volume(args.parentName());
std::string childName = args.value<std::string>("ChildName");
childName = ns.prepend(childName);
dd4hep::Volume child = ns.volume(childName);
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "DDMuonAngular: Parameters for positioning:: n " << n << " Start, Step "
<< convertRadToDeg(startAngle) << ", " << convertRadToDeg(stepAngle) << ", zoffset "
<< cms::convert2mm(zoffset) << ", RotNameSpace " << rotns;
edm::LogVerbatim("MuonGeom") << "DDMuonAngular: Parent " << mother.name() << "\tChild " << child.name()
<< " NameSpace " << ns.name();
#endif
double phi = startAngle;
int copyNo = startCopyNo;
for (int i = 0; i < n; ++i) {
double phitmp = phi;
if (phitmp >= 2._pi)
phitmp -= 2._pi;
dd4hep::Rotation3D rotation = cms::makeRotation3D(90._deg, phitmp, 90._deg, 90._deg + phitmp, 0., 0.);
std::string rotstr = ns.nsName(child.name()) + std::to_string(phitmp * 10.);
auto irot = context.rotations.find(ns.prepend(rotstr));
if (irot != context.rotations.end()) {
rotation = ns.rotation(ns.prepend(rotstr));
}
dd4hep::Position tran(0., 0., zoffset);
mother.placeVolume(child, copyNo, dd4hep::Transform3D(rotation, tran));
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("MuonGeom") << "DDMuonAngular:" << child.name() << " number " << copyNo << " positioned in "
<< mother.name() << " at (0,0," << cms::convert2mm(zoffset) << ") with " << rotstr
<< ": " << rotation;
#endif
phi += stepAngle;
copyNo += incrCopyNo;
}
return cms::s_executed;
}
// first argument is the type from the xml file
DECLARE_DDCMS_DETELEMENT(DDCMS_muon_DDMuonAngular, algorithm)
|