|
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 "DetectorDescription/DDCMS/interface/DDPlugins.h"
#include "DetectorDescription/DDCMS/interface/DDutils.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DD4hep/DetFactoryHelper.h"
//#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 DDHCalLinearXY.h
int numberX = args.value<int>("NumberX"); //Slices along X
int numberY = args.value<int>("NumberY"); //Slices along Y
double deltaX = args.value<double>("DeltaX"); //Step along X
double deltaY = args.value<double>("DeltaY"); //Step along Y
std::vector<double> centre = args.value<std::vector<double> >("Center");
dd4hep::Volume parent = ns.volume(args.parentName());
std::vector<std::string> childName = args.value<std::vector<std::string> >("Child");
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HCalGeom") << "DDHCalLinearXY: Parent " << parent.name() << " \twith " << childName.size()
<< " children";
unsigned int k(0);
for (auto const& child : childName) {
edm::LogVerbatim("HCalGeom") << "DDHCalLinearXY: Child[" << k << "] = " << child;
++k;
}
edm::LogVerbatim("HCalGeom") << "DDHCalLinearXY: Number along X/Y " << numberX << "/" << numberY
<< "\tDelta along X/Y " << cms::convert2mm(deltaX) << "/" << cms::convert2mm(deltaY)
<< "\tCentre (" << cms::convert2mm(centre[0]) << ", " << cms::convert2mm(centre[1])
<< "," << cms::convert2mm(centre[2]);
#endif
double xoff = centre[0] - (numberX - 1) * 0.5 * deltaX;
double yoff = centre[1] - (numberY - 1) * 0.5 * deltaY;
unsigned int copy = 0;
for (int i = 0; i < numberX; ++i) {
for (int j = 0; j < numberY; ++j) {
unsigned int k = (childName.size() == 1) ? 0 : copy;
++copy;
if (k < childName.size() && (childName[k] != " " && childName[k] != "Null")) {
std::string child = childName[k];
dd4hep::Position tran(xoff + i * deltaX, yoff + j * deltaY, centre[2]);
parent.placeVolume(ns.volume(child), copy, tran);
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HCalGeom") << "DDHCalLinearXY: " << child << " number " << copy << " positioned in "
<< parent.name() << " at (" << cms::convert2mm((xoff + i * deltaX)) << ", "
<< cms::convert2mm((yoff + j * deltaY)) << ", " << cms::convert2mm(centre[2])
<< ") with no rotation";
#endif
} else {
#ifdef EDM_ML_DEBUG
edm::LogVerbatim("HCalGeom") << "DDHCalLinearXY: no child placed for [" << copy << "]";
#endif
}
}
}
return cms::s_executed;
}
// first argument is the type from the xml file
DECLARE_DDCMS_DETELEMENT(DDCMS_hcal_DDHCalLinearXY, algorithm)
|