|
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
|
#include "DD4hep/DetFactoryHelper.h"
#include "DataFormats/Math/interface/CMSUnits.h"
#include "DetectorDescription/DDCMS/interface/DDPlugins.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
using namespace std;
using namespace dd4hep;
using namespace cms;
using namespace cms_units::operators;
namespace {
long algorithm(Detector& /* description */, cms::DDParsingContext& ctxt, xml_h e) {
DDNamespace ns(ctxt, e, true);
DDAlgoArguments args(ctxt, e);
Volume mother = ns.volume(args.parentName());
string childName = args.value<string>("ChildName");
Volume child = ns.volume(childName);
string parentName = args.parentName();
int n = args.value<int>("N");
int startCopyNo = args.value<int>("StartCopyNo");
int incrCopyNo = args.value<int>("IncrCopyNo");
double rangeAngle = args.value<double>("RangeAngle");
double startAngle = args.value<double>("StartAngle");
double radius = args.value<double>("Radius");
vector<double> center = args.value<vector<double> >("Center");
vector<double> phiAngles = args.value<vector<double> >("phiAngleValues");
vector<double> radiusValues = args.value<vector<double> >("radiusValues");
vector<double> yawAngles = args.value<vector<double> >("yawAngleValues");
bool isZPlus = args.value<int>("IsZPlus") == 1;
double tiltAngle = args.value<double>("TiltAngle");
bool isFlipped = args.value<int>("IsFlipped") == 1;
edm::LogVerbatim("TrackerGeom") << "DDTrackerIrregularRingAlgo debug: Parameters for position"
<< "ing:: n " << n << " Start, Range " << convertRadToDeg(startAngle) << " "
<< convertRadToDeg(rangeAngle) << " Radius " << radius << " Centre " << center[0]
<< ", " << center[1] << ", " << center[2];
edm::LogVerbatim("TrackerGeom") << "DDTrackerIrregularRingAlgo debug: Parent " << parentName << "\tChild "
<< childName << " NameSpace " << ns.name();
Rotation3D flipMatrix, tiltMatrix, phiRotMatrix, globalRotMatrix, phiOwnAxisRotMatrix; // Identity matrix
// flipMatrix calculus
if (isFlipped) {
edm::LogVerbatim("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: "
<< "\t90., 180., "
<< "90., 90., "
<< "180., 0.";
flipMatrix = makeRotation3D(90._deg, 180._deg, 90._deg, 90._deg, 180._deg, 0._deg);
}
// tiltMatrix calculus
if (isZPlus) {
edm::LogVerbatim("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: "
<< "\t90., 90., " << convertRadToDeg(tiltAngle) << ", 180., "
<< convertRadToDeg(90._deg - tiltAngle) << ", 0.";
tiltMatrix = makeRotation3D(90._deg, 90._deg, tiltAngle, 180._deg, 90._deg - tiltAngle, 0._deg);
if (isFlipped) {
tiltMatrix *= flipMatrix;
}
} else {
edm::LogVerbatim("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: "
<< "\t90., 90., " << convertRadToDeg(tiltAngle) << ", 0., "
<< convertRadToDeg(90._deg + tiltAngle) << ", 0.";
tiltMatrix = makeRotation3D(90._deg, 90._deg, tiltAngle, 0._deg, 90._deg + tiltAngle, 0._deg);
if (isFlipped) {
tiltMatrix *= flipMatrix;
}
}
// Loops for all phi values
double theta = 90._deg;
int copy = startCopyNo;
for (int i = 0; i < n; ++i) {
// phiRotMatrix calculus
double phix = convertDegToRad(phiAngles.at(i));
double phix_ownaxis = convertDegToRad(yawAngles.at(i));
radius = radiusValues.at(i);
double phiy = phix + 90._deg;
double phiy_ownaxis = phix_ownaxis + 90._deg;
if (phix_ownaxis != 0.) {
LogDebug("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: "
<< "\t90., " << convertRadToDeg(phix_ownaxis) << ", 90.,"
<< convertRadToDeg(phiy_ownaxis) << ", 0., 0.";
phiOwnAxisRotMatrix = makeRotation3D(theta, phix_ownaxis, theta, phiy_ownaxis, 0., 0.);
}
if (phix != 0.) {
edm::LogVerbatim("TrackerGeom") << "DDTrackerIrregularRingAlgo test: Creating a new rotation: "
<< "\t90., " << convertRadToDeg(phix) << ", 90.," << convertRadToDeg(phiy)
<< ", 0., 0.";
phiRotMatrix = makeRotation3D(theta, phix, theta, phiy, 0., 0.);
}
globalRotMatrix = phiOwnAxisRotMatrix * phiRotMatrix * tiltMatrix;
// translation def
double xpos = radius * cos(phix) + center[0];
double ypos = radius * sin(phix) + center[1];
double zpos = center[2];
Position tran(xpos, ypos, zpos);
// Positions child with respect to parent
mother.placeVolume(child, copy, Transform3D(globalRotMatrix, tran));
edm::LogVerbatim("TrackerGeom") << "DDTrackerIrregularRingAlgo test " << child.data()->GetName() << " number "
<< copy << " positioned in " << mother.data()->GetName() << " at " << tran
<< " with " << globalRotMatrix;
copy += incrCopyNo;
}
return s_executed;
}
} // namespace
// first argument is the type from the xml file
DECLARE_DDCMS_DETELEMENT(DDCMS_track_DDTrackerIrregularRingAlgo, ::algorithm)
|