|
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
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
|
#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;
static long algorithm(Detector& /* description */, cms::DDParsingContext& ctxt, xml_h e) {
cms::DDNamespace ns(ctxt, e, true);
DDAlgoArguments args(ctxt, e);
string mother = args.parentName();
string genMat = args.str("GeneralMaterial"); //General material name
int detectorN = args.integer("DetectorNumber"); //Detector planes
double moduleThick = args.dble("ModuleThick"); //Module thickness
double detTilt = args.dble("DetTilt"); //Tilt of stereo detector
double fullHeight = args.dble("FullHeight"); //Height
double dlTop = args.dble("DlTop"); //Width at top of wafer
double dlBottom = args.dble("DlBottom"); //Width at bottom of wafer
double dlHybrid = args.dble("DlHybrid"); //Width at the hybrid end
bool doComponents = ::toupper(args.str("DoComponents")[0]) != 'N'; //Components to be made
string boxFrameName = args.str("BoxFrameName"); //Top frame name
string boxFrameMat = args.str("BoxFrameMaterial"); // material
double boxFrameHeight = args.dble("BoxFrameHeight"); // height
double boxFrameThick = args.dble("BoxFrameThick"); // thickness
double boxFrameWidth = args.dble("BoxFrameWidth"); // extra width
double bottomFrameHeight = args.dble("BottomFrameHeight"); //Bottom of the frame
double bottomFrameOver = args.dble("BottomFrameOver"); // overlap
double topFrameHeight = args.dble("TopFrameHeight"); //Top of the frame
double topFrameOver = args.dble("TopFrameOver"); // overlap
vector<string> sideFrameName = args.vecStr("SideFrameName"); //Side frame name
string sideFrameMat = args.str("SideFrameMaterial"); // material
double sideFrameWidth = args.dble("SideFrameWidth"); // width
double sideFrameThick = args.dble("SideFrameThick"); // thickness
double sideFrameOver = args.dble("SideFrameOver"); // overlap (wrt wafer)
vector<string> holeFrameName = args.vecStr("HoleFrameName"); //Hole in the frame name
vector<string> holeFrameRot = args.vecStr("HoleFrameRotation"); // Rotation matrix
vector<string> kaptonName = args.vecStr("KaptonName"); //Kapton circuit name
string kaptonMat = args.str("KaptonMaterial"); // material
double kaptonThick = args.dble("KaptonThick"); // thickness
double kaptonOver = args.dble("KaptonOver"); // overlap (wrt Wafer)
vector<string> holeKaptonName = args.vecStr("HoleKaptonName"); //Hole in the kapton circuit name
vector<string> holeKaptonRot = args.vecStr("HoleKaptonRotation"); // Rotation matrix
vector<string> waferName = args.vecStr("WaferName"); //Wafer name
string waferMat = args.str("WaferMaterial"); // material
double sideWidthTop = args.dble("SideWidthTop"); // width on the side Top
double sideWidthBottom = args.dble("SideWidthBottom"); // Bottom
vector<string> activeName = args.vecStr("ActiveName"); //Sensitive name
string activeMat = args.str("ActiveMaterial"); // material
double activeHeight = args.dble("ActiveHeight"); // height
vector<double> waferThick =
args.vecDble("WaferThick"); // wafer thickness (active = wafer - backplane)
string activeRot = args.str("ActiveRotation"); // Rotation matrix
vector<double> backplaneThick = args.vecDble("BackPlaneThick"); // thickness
string hybridName = args.str("HybridName"); //Hybrid name
string hybridMat = args.str("HybridMaterial"); // material
double hybridHeight = args.dble("HybridHeight"); // height
double hybridWidth = args.dble("HybridWidth"); // width
double hybridThick = args.dble("HybridThick"); // thickness
vector<string> pitchName = args.vecStr("PitchName"); //Pitch adapter name
string pitchMat = args.str("PitchMaterial"); // material
double pitchHeight = args.dble("PitchHeight"); // height
double pitchThick = args.dble("PitchThick"); // thickness
double pitchStereoTol = args.dble("PitchStereoTolerance"); // tolerance in dimensions of the stereo
string coolName = args.str("CoolInsertName"); // Cool insert name
string coolMat = args.str("CoolInsertMaterial"); // material
double coolHeight = args.dble("CoolInsertHeight"); // height
double coolThick = args.dble("CoolInsertThick"); // thickness
double coolWidth = args.dble("CoolInsertWidth"); // width
edm::LogVerbatim("TIDGeom") << "Parent " << mother << " General Material " << genMat << " Detector Planes "
<< detectorN;
edm::LogVerbatim("TIDGeom") << "ModuleThick " << moduleThick << " Detector Tilt " << convertRadToDeg(detTilt)
<< " Height " << fullHeight << " dl(Top) " << dlTop << " dl(Bottom) " << dlBottom
<< " dl(Hybrid) " << dlHybrid << " doComponents " << doComponents;
edm::LogVerbatim("TIDGeom") << "" << boxFrameName << " Material " << boxFrameMat << " Thickness " << boxFrameThick
<< " width " << boxFrameWidth << " height " << boxFrameHeight
<< " Extra Height at Bottom " << bottomFrameHeight << " Overlap " << bottomFrameOver;
for (int i = 0; i < detectorN; i++)
edm::LogVerbatim("TIDGeom") << sideFrameName[i] << " Material " << sideFrameMat << " Width " << sideFrameWidth
<< " Thickness " << sideFrameThick << " Overlap " << sideFrameOver << " Hole "
<< holeFrameName[i];
for (int i = 0; i < detectorN; i++)
edm::LogVerbatim("TIDGeom") << kaptonName[i] << " Material " << kaptonMat << " Thickness " << kaptonThick
<< " Overlap " << kaptonOver << " Hole " << holeKaptonName[i];
edm::LogVerbatim("TIDGeom") << "Wafer Material " << waferMat << " Side Width Top " << sideWidthTop
<< " Side Width Bottom " << sideWidthBottom;
for (int i = 0; i < detectorN; i++)
edm::LogVerbatim("TIDGeom") << "\twaferName[" << i << "] = " << waferName[i];
edm::LogVerbatim("TIDGeom") << "Active Material " << activeMat << " Height " << activeHeight << " rotated by "
<< activeRot;
for (int i = 0; i < detectorN; i++)
edm::LogVerbatim("TIDGeom") << " translated by (0," << -0.5 * backplaneThick[i] << ",0)\tactiveName[" << i
<< "] = " << activeName[i] << " of thickness " << waferThick[i] - backplaneThick[i];
edm::LogVerbatim("TIDGeom") << "" << hybridName << " Material " << hybridMat << " Height " << hybridHeight
<< " Width " << hybridWidth << " Thickness " << hybridThick;
edm::LogVerbatim("TIDGeom") << "Pitch Adapter Material " << pitchMat << " Height " << pitchHeight << " Thickness "
<< pitchThick;
for (int i = 0; i < detectorN; i++)
edm::LogVerbatim("TIDGeom") << "\tpitchName[" << i << "] = " << pitchName[i];
edm::LogVerbatim("TIDGeom") << "Cool Element Material " << coolMat << " Height " << coolHeight << " Thickness "
<< coolThick << " Width " << coolWidth;
string name = mother;
double sidfr = sideFrameWidth - sideFrameOver; // width of side frame on the sides of module
double botfr; // width of side frame at the the bottom of the modules
double topfr; // width of side frame at the the top of the modules
double kaptonHeight;
if (dlHybrid > dlTop) {
// ring 1, ring 2
topfr = topFrameHeight - pitchHeight - topFrameOver;
botfr = bottomFrameHeight - bottomFrameOver;
kaptonHeight = fullHeight + botfr;
} else {
// ring 3
topfr = topFrameHeight - topFrameOver;
botfr = bottomFrameHeight - bottomFrameOver - pitchHeight;
kaptonHeight = fullHeight + topfr;
}
double sideFrameHeight = fullHeight + pitchHeight + botfr + topfr;
double kaptonWidth = sidfr + kaptonOver;
double dxbot = 0.5 * dlBottom + sidfr;
double dxtop = 0.5 * dlTop + sidfr;
double dxtopenv, dxbotenv; // top/bot width of the module envelope trap
// Envelope
if (dlHybrid > dlTop) {
// ring 1, ring 2
dxtopenv = dxbot + (dxtop - dxbot) * (fullHeight + pitchHeight + topfr + hybridHeight) / fullHeight;
dxbotenv = dxtop - (dxtop - dxbot) * (fullHeight + botfr) / fullHeight;
} else {
// ring 3
dxtopenv = dxbot + (dxtop - dxbot) * (fullHeight + topfr) / fullHeight;
dxbotenv = dxbot;
}
double bl1 = dxbotenv;
double bl2 = dxtopenv;
double h1 = 0.5 * moduleThick;
double dx, dy;
double dz = 0.5 * (boxFrameHeight + sideFrameHeight);
Solid solid = ns.addSolidNS(ns.prepend(name), Trap(dz, 0, 0, h1, bl1, bl1, 0, h1, bl2, bl2, 0));
/* Volume module = */ ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(genMat)));
edm::LogVerbatim("TIDGeom") << solid.name() << " Trap made of " << genMat << " of dimensions " << dz << ", 0, 0, "
<< h1 << ", " << bl1 << ", " << bl1 << ", 0, " << h1 << ", " << bl2 << ", " << bl2
<< ", 0";
if (doComponents) {
//Box frame
name = boxFrameName;
dx = 0.5 * boxFrameWidth;
dy = 0.5 * boxFrameThick;
dz = 0.5 * boxFrameHeight;
solid = ns.addSolidNS(ns.prepend(name), Box(dx, dy, dz));
edm::LogVerbatim("TIDGeom") << solid.name() << " Box made of " << boxFrameMat << " of dimensions " << dx << ", "
<< dy << ", " << dz;
/* Volume boxFrame = */ ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(boxFrameMat)));
// Hybrid
name = hybridName;
dx = 0.5 * hybridWidth;
dy = 0.5 * hybridThick;
dz = 0.5 * hybridHeight;
solid = ns.addSolidNS(ns.prepend(name), Box(dx, dy, dz));
edm::LogVerbatim("TIDGeom") << solid.name() << " Box made of " << hybridMat << " of dimensions " << dx << ", " << dy
<< ", " << dz;
/* Volume hybrid = */ ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(hybridMat)));
// Cool Insert
name = coolName;
dx = 0.5 * coolWidth;
dy = 0.5 * coolThick;
dz = 0.5 * coolHeight;
solid = ns.addSolidNS(ns.prepend(name), Box(dx, dy, dz));
edm::LogVerbatim("TIDGeom") << solid.name() << " Box made of " << coolMat << " of dimensions " << dx << ", " << dy
<< ", " << dz;
/* Volume cool = */ ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(coolMat)));
// Loop over detectors to be placed
for (int k = 0; k < detectorN; k++) {
double bbl1, bbl2; // perhaps useless (bl1 enough)
// Frame Sides
name = sideFrameName[k];
if (dlHybrid > dlTop) {
// ring 1, ring 2
bbl1 = dxtop - (dxtop - dxbot) * (fullHeight + botfr) / fullHeight;
bbl2 = dxbot + (dxtop - dxbot) * (fullHeight + pitchHeight + topfr) / fullHeight;
} else {
// ring 3
bbl1 = dxtop - (dxtop - dxbot) * (fullHeight + pitchHeight + botfr) / fullHeight;
bbl2 = dxbot + (dxtop - dxbot) * (fullHeight + topfr) / fullHeight;
}
h1 = 0.5 * sideFrameThick;
dz = 0.5 * sideFrameHeight;
solid = ns.addSolidNS(ns.prepend(name), Trap(dz, 0., 0., h1, bbl1, bbl1, 0., h1, bbl2, bbl2, 0.));
edm::LogVerbatim("TIDGeom") << solid.name() << " Trap made of " << sideFrameMat << " of dimensions " << dz
<< ", 0, 0, " << h1 << ", " << bbl1 << ", " << bbl1 << ", 0, " << h1 << ", " << bbl2
<< ", " << bbl2 << ", 0";
Volume sideFrame = ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(sideFrameMat)));
std::string rotstr, rotns;
Rotation3D rot;
// Hole in the frame below the wafer
name = holeFrameName[k];
double xpos, zpos;
dz = fullHeight - bottomFrameOver - topFrameOver;
bbl1 = dxbot - sideFrameWidth + bottomFrameOver * (dxtop - dxbot) / fullHeight;
bbl2 = dxtop - sideFrameWidth - topFrameOver * (dxtop - dxbot) / fullHeight;
if (dlHybrid > dlTop) {
// ring 1, ring 2
zpos = -(topFrameHeight + 0.5 * dz - 0.5 * sideFrameHeight);
} else {
// ring 3
zpos = bottomFrameHeight + 0.5 * dz - 0.5 * sideFrameHeight;
}
dz /= 2.;
solid = ns.addSolidNS(ns.prepend(name), Trap(dz, 0, 0, h1, bbl1, bbl1, 0, h1, bbl2, bbl2, 0));
edm::LogVerbatim("TIDGeom") << solid.name() << " Trap made of " << genMat << " of dimensions " << dz << ", 0, 0, "
<< h1 << ", " << bbl1 << ", " << bbl1 << ", 0, " << h1 << ", " << bbl2 << ", " << bbl2
<< ", 0";
Volume holeFrame = ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(genMat)));
rot = ns.rotation(holeFrameRot[k]);
sideFrame.placeVolume(holeFrame, 1, Transform3D(rot, Position(0e0, 0e0, zpos))); // copyNr=1
edm::LogVerbatim("TIDGeom") << holeFrame.name() << " number 1 positioned in " << sideFrame.name() << " at (0,0,"
<< zpos << ") with no rotation";
// Kapton circuit
double kaptonExtraHeight = 0; // kapton extra height in the stereo
if (dlHybrid > dlTop) {
// ring 1, ring 2
bbl1 = dxtop - (dxtop - dxbot) * (fullHeight + botfr) / fullHeight;
if (k == 1) {
kaptonExtraHeight = dlTop * sin(detTilt) - fullHeight * (1 - cos(detTilt));
kaptonExtraHeight = 0.5 * fabs(kaptonExtraHeight);
bbl2 = dxbot + (dxtop - dxbot) * (fullHeight + kaptonExtraHeight) / fullHeight;
} else {
bbl2 = dxtop;
}
} else {
// ring 3
bbl2 = dxbot + (dxtop - dxbot) * (fullHeight + topfr) / fullHeight;
if (k == 1) {
kaptonExtraHeight = dlBottom * sin(detTilt) - fullHeight * (1 - cos(detTilt));
kaptonExtraHeight = 0.5 * fabs(kaptonExtraHeight);
bbl1 = dxtop - (dxtop - dxbot) * (fullHeight + kaptonExtraHeight) / fullHeight;
} else {
bbl1 = dxbot;
}
}
h1 = 0.5 * kaptonThick;
dz = 0.5 * (kaptonHeight + kaptonExtraHeight);
// For the stereo create the uncut solid, the solid to be removed and then the subtraction solid
if (k == 1) {
// Uncut solid
Solid solidUncut =
ns.addSolidNS(ns.prepend(kaptonName[k] + "Uncut"), Trap(dz, 0., 0., h1, bbl1, bbl1, 0., h1, bbl2, bbl2, 0));
// Piece to be cut
dz = (dlHybrid > dlTop) ? 0.5 * dlTop : 0.5 * dlBottom;
h1 = 0.5 * kaptonThick;
bbl1 = fabs(dz * sin(detTilt));
bbl2 = bbl1 * 0.000001;
double thet = atan((bbl1 - bbl2) / (2 * dz));
Solid solidCut =
ns.addSolidNS(ns.prepend(kaptonName[k] + "Cut"), Trap(dz, thet, 0., h1, bbl1, bbl1, 0., h1, bbl2, bbl2, 0));
// Subtraction Solid
name = kaptonName[k];
rot = ns.rotation("tidmodpar:9PYX");
xpos = -0.5 * fullHeight * sin(detTilt);
zpos = 0.5 * kaptonHeight - bbl2;
solid = ns.addSolidNS(ns.prepend(name),
SubtractionSolid(solidUncut, solidCut, Transform3D(rot, Position(xpos, 0.0, zpos))));
} else {
name = kaptonName[k];
solid = ns.addSolidNS(ns.prepend(name), Trap(dz, 0., 0., h1, bbl1, bbl1, 0., h1, bbl2, bbl2, 0.));
}
Volume kapton = ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(kaptonMat)));
edm::LogVerbatim("TIDGeom") << solid.name() << " SUBTRACTION SOLID Trap made of " << kaptonMat
<< " of dimensions " << dz << ", 0, 0, " << h1 << ", " << bbl1 << ", " << bbl1
<< ", 0, " << h1 << ", " << bbl2 << ", " << bbl2 << ", 0";
// Hole in the kapton below the wafer
name = holeKaptonName[k];
dz = fullHeight - kaptonOver;
xpos = 0;
if (dlHybrid > dlTop) {
// ring 1, ring 2
bbl1 = dxbot - kaptonWidth + kaptonOver * (dxtop - dxbot) / fullHeight;
bbl2 = dxtop - kaptonWidth;
zpos = 0.5 * (kaptonHeight - kaptonExtraHeight - dz);
if (k == 1) {
zpos -= 0.5 * kaptonOver * (1 - cos(detTilt));
xpos = -0.5 * kaptonOver * sin(detTilt);
}
} else {
// ring 3
bbl1 = dxbot - kaptonWidth;
bbl2 = dxtop - kaptonWidth - kaptonOver * (dxtop - dxbot) / fullHeight;
zpos = -0.5 * (kaptonHeight - kaptonExtraHeight - dz);
}
dz /= 2.;
solid = ns.addSolidNS(ns.prepend(name), Trap(dz, 0., 0., h1, bbl1, bbl1, 0., h1, bbl2, bbl2, 0.));
edm::LogVerbatim("TIDGeom") << solid.name() << " Trap made of " << genMat << " of dimensions " << dz << ", 0, 0, "
<< h1 << ", " << bbl1 << ", " << bbl1 << ", 0, " << h1 << ", " << bbl2 << ", " << bbl2
<< ", 0";
Volume holeKapton = ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(genMat)));
rot = ns.rotation(holeKaptonRot[k]);
kapton.placeVolume(holeKapton, 1, Transform3D(rot, Position(xpos, 0.0, zpos)));
edm::LogVerbatim("TIDGeom") << holeKapton.name() << " number 1 positioned in " << kapton.name() << " at (0,0,"
<< zpos << ") with no rotation";
// Wafer
name = waferName[k];
if (k == 0 && dlHybrid < dlTop) {
bl1 = 0.5 * dlTop;
bl2 = 0.5 * dlBottom;
} else {
bl1 = 0.5 * dlBottom;
bl2 = 0.5 * dlTop;
}
h1 = 0.5 * waferThick[k];
dz = 0.5 * fullHeight;
solid = ns.addSolidNS(ns.prepend(name), Trap(dz, 0, 0, h1, bl1, bl1, 0, h1, bl2, bl2, 0));
edm::LogVerbatim("TIDGeom") << solid.name() << " Trap made of " << waferMat << " of dimensions " << dz
<< ", 0, 0, " << h1 << ", " << bl1 << ", " << bl1 << ", 0, " << h1 << ", " << bl2
<< ", " << bl2 << ", 0";
Volume wafer = ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(waferMat)));
// Active
name = activeName[k];
if (k == 0 && dlHybrid < dlTop) {
bl1 -= sideWidthTop;
bl2 -= sideWidthBottom;
} else {
bl1 -= sideWidthBottom;
bl2 -= sideWidthTop;
}
dz = 0.5 * (waferThick[k] - backplaneThick[k]); // inactive backplane
h1 = 0.5 * activeHeight;
solid = ns.addSolidNS(ns.prepend(name), Trap(dz, 0, 0, h1, bl2, bl1, 0, h1, bl2, bl1, 0));
edm::LogVerbatim("TIDGeom") << solid.name() << " Trap made of " << activeMat << " of dimensions " << dz
<< ", 0, 0, " << h1 << ", " << bl2 << ", " << bl1 << ", 0, " << h1 << ", " << bl2
<< ", " << bl1 << ", 0";
Volume active = ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(activeMat)));
rot = ns.rotation(activeRot);
Position tran(0.0, -0.5 * backplaneThick[k], 0.0); // from the definition of the wafer local axes
wafer.placeVolume(active, 1, Transform3D(rot, tran)); // inactive backplane copyNr=1
edm::LogVerbatim("TIDGeom") << "DDTIDModuleAlgo test: " << active.name() << " number 1 positioned in "
<< wafer.name() << " at " << tran << " with " << rot;
//Pitch Adapter
name = pitchName[k];
if (dlHybrid > dlTop) {
dz = 0.5 * dlTop;
} else {
dz = 0.5 * dlBottom;
}
if (k == 0) {
dx = dz;
dy = 0.5 * pitchThick;
dz = 0.5 * pitchHeight;
solid = ns.addSolidNS(ns.prepend(name), Box(dx, dy, dz));
edm::LogVerbatim("TIDGeom") << solid.name() << " Box made of " << pitchMat << " of dimensions"
<< " " << dx << ", " << dy << ", " << dz;
} else {
h1 = 0.5 * pitchThick;
bl1 = 0.5 * pitchHeight + 0.5 * dz * sin(detTilt);
bl2 = 0.5 * pitchHeight - 0.5 * dz * sin(detTilt);
dz -= 0.5 * pitchStereoTol;
bl1 -= pitchStereoTol;
bl2 -= pitchStereoTol;
double thet = atan((bl1 - bl2) / (2. * dz));
solid = ns.addSolidNS(ns.prepend(name), Trap(dz, thet, 0, h1, bl1, bl1, 0, h1, bl2, bl2, 0));
edm::LogVerbatim("TIDGeom") << solid.name() << " Trap made of " << pitchMat << " of "
<< "dimensions " << dz << ", " << convertRadToDeg(thet) << ", 0, " << h1 << ", "
<< bl1 << ", " << bl1 << ", 0, " << h1 << ", " << bl2 << ", " << bl2 << ", 0";
}
/* Volume pa = */ ns.addVolumeNS(Volume(ns.prepend(name), solid, ns.material(pitchMat)));
}
}
edm::LogVerbatim("TIDGeom") << "<<== End of DDTIDModuleAlgo construction ...";
return 1;
}
// first argument is the type from the xml file
DECLARE_DDCMS_DETELEMENT(DDCMS_track_DDTIDModuleAlgo, algorithm)
|