|
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
|
#include "DetectorDescription/Core/interface/DDFilteredView.h"
#include "DetectorDescription/DDCMS/interface/DDFilteredView.h"
#include "Geometry/TrackerNumberingBuilder/plugins/CmsDetConstruction.h"
#include "Geometry/TrackerNumberingBuilder/plugins/ExtractStringFromDDD.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <memory>
template <class FilteredView>
void CmsDetConstruction<FilteredView>::buildSmallDetsforGlued(FilteredView& fv,
GeometricDet* mother,
const std::string& attribute) {
auto det = std::make_unique<GeometricDet>(&fv,
CmsTrackerLevelBuilder<FilteredView>::theCmsTrackerStringToEnum.type(
ExtractStringFromDDD<FilteredView>::getString(attribute, &fv)));
if (det->stereo()) {
uint32_t temp = 1;
det->setGeographicalID(DetId(temp));
} else {
uint32_t temp = 2;
det->setGeographicalID(DetId(temp));
}
mother->addComponent(det.release());
}
template <class FilteredView>
void CmsDetConstruction<FilteredView>::buildSmallDetsforStack(FilteredView& fv,
GeometricDet* mother,
const std::string& attribute) {
auto det = std::make_unique<GeometricDet>(&fv,
CmsTrackerLevelBuilder<FilteredView>::theCmsTrackerStringToEnum.type(
ExtractStringFromDDD<FilteredView>::getString(attribute, &fv)));
if (det->isLowerSensor()) {
uint32_t temp = 1;
det->setGeographicalID(DetId(temp));
} else if (det->isUpperSensor()) {
uint32_t temp = 2;
det->setGeographicalID(DetId(temp));
} else {
edm::LogError("DetConstruction") << " module defined in a Stack but not upper either lower!? ";
}
mother->addComponent(det.release());
}
template <class FilteredView>
void CmsDetConstruction<FilteredView>::buildSmallDetsfor3D(FilteredView& fv,
GeometricDet* mother,
const std::string& attribute) {
auto det = std::make_unique<GeometricDet>(&fv,
CmsTrackerLevelBuilder<FilteredView>::theCmsTrackerStringToEnum.type(
ExtractStringFromDDD<FilteredView>::getString(attribute, &fv)));
if (det->isFirstSensor()) {
uint32_t temp = 1;
det->setGeographicalID(DetId(temp));
} else if (det->isSecondSensor()) {
uint32_t temp = 2;
det->setGeographicalID(DetId(temp));
} else {
edm::LogError("DetConstruction") << " module defined in a 3D module but not first or second sensor!? ";
}
mother->addComponent(det.release());
}
/*
* OLD DD.
* Module with 2 sensors: calculate the sensor local ID, and add the sensor to its mother volume (module).
* Module with 1 sensor: just add the sensor to to its mother volume (ladder).
*/
template <>
void CmsDetConstruction<DDFilteredView>::buildComponent(DDFilteredView& fv,
GeometricDet* mother,
const std::string& attribute) {
// Mother volume
// Module with 2 sensors: the mother volume is the module volume.
// Module with 1 sensor: the mother volume is the ladder volume.
const std::string& myTopologicalNameInXMLs = ExtractStringFromDDD<DDFilteredView>::getString(attribute, &fv);
const GeometricDet::GDEnumType& myTopologicalType =
CmsTrackerLevelBuilder<DDFilteredView>::theCmsTrackerStringToEnum.type(myTopologicalNameInXMLs);
auto det = std::make_unique<GeometricDet>(&fv, myTopologicalType);
const bool isPhase1ModuleWith2Sensors = (myTopologicalType == GeometricDet::mergedDet);
const bool isPhase2ModuleWith2Sensors = (myTopologicalType == GeometricDet::OTPhase2Stack);
const bool isPhase2BarrelModuleWith2Sensors = (myTopologicalType == GeometricDet::ITPhase2Combined);
// CASE A: MODULE HAS 2 SENSORS
if (isPhase1ModuleWith2Sensors || isPhase2ModuleWith2Sensors || isPhase2BarrelModuleWith2Sensors) {
// Go down in hierarchy: from module to sensor
bool dodets = fv.firstChild(); // very important
while (dodets) {
// PHASE 1 (MERGEDDET)
if (isPhase1ModuleWith2Sensors) {
buildSmallDetsforGlued(fv, det.get(), attribute);
}
// PHASE 2 (STACKDET)
else if (isPhase2ModuleWith2Sensors) {
buildSmallDetsforStack(fv, det.get(), attribute);
} else if (isPhase2BarrelModuleWith2Sensors) {
buildSmallDetsfor3D(fv, det.get(), attribute);
}
dodets = fv.nextSibling();
}
fv.parent();
}
// CASE B: MODULE HAS 1 SENSOR: NOTHING SPECIFIC TO DO
// Indeed, we are not going to sort sensors within module, if there is only 1 sensor!
// ALL CASES: add sensor to its mother volume (module or ladder).
mother->addComponent(det.release());
}
/*
* DD4hep.
* Module with 2 sensors: calculate the sensor local ID, and add the sensor to its mother volume (module).
* Module with 1 sensor: just add the sensor to its mother volume (ladder).
*/
template <>
void CmsDetConstruction<cms::DDFilteredView>::buildComponent(cms::DDFilteredView& fv,
GeometricDet* mother,
const std::string& attribute) {
// Mother volume
// Module with 2 sensors: the mother volume is the module volume.
// Module with 1 sensor: the mother volume is the ladder volume.
const std::string& myTopologicalNameInXMLs = ExtractStringFromDDD<cms::DDFilteredView>::getString(attribute, &fv);
const GeometricDet::GDEnumType& myTopologicalType =
CmsTrackerLevelBuilder<cms::DDFilteredView>::theCmsTrackerStringToEnum.type(myTopologicalNameInXMLs);
auto det = std::make_unique<GeometricDet>(&fv, myTopologicalType);
const bool isPhase1ModuleWith2Sensors = (myTopologicalType == GeometricDet::mergedDet);
const bool isPhase2ModuleWith2Sensors = (myTopologicalType == GeometricDet::OTPhase2Stack);
const bool isPhase2BarrelModuleWith2Sensors = (myTopologicalType == GeometricDet::ITPhase2Combined);
// CASE A: MODULE HAS 2 SENSORS
if (isPhase1ModuleWith2Sensors || isPhase2ModuleWith2Sensors || isPhase2BarrelModuleWith2Sensors) {
// Go down in hierarchy: from module to sensor
if (!fv.firstChild()) { // very important
edm::LogError("CmsDetConstruction::buildComponent. Cannot go down to sensor volume.");
return;
}
// This is the sensor hierarchy level
const int sensorHierarchyLevel = fv.level();
// Loop on all siblings (ie, on all sensors)
while (fv.level() == sensorHierarchyLevel) {
// PHASE 1 (MERGEDDET)
if (isPhase1ModuleWith2Sensors) {
buildSmallDetsforGlued(fv, det.get(), attribute);
}
// PHASE 2 (STACKDET)
else if (isPhase2ModuleWith2Sensors) {
buildSmallDetsforStack(fv, det.get(), attribute);
} else if (isPhase2BarrelModuleWith2Sensors) {
buildSmallDetsfor3D(fv, det.get(), attribute);
}
// Go to the next volume in FilteredView.
// NB: If this volume is another sensor of the same module, will stay in the loop.
// Otherwise, it is very important to access the next volume to be treated anyway.
fv.firstChild();
}
}
// CASE B: MODULE HAS 1 SENSOR: NOTHING SPECIFIC TO DO
// Indeed, we are not going to sort sensors within module, if there is only 1 sensor!
else {
// Go to the next volume in FilteredView.
fv.firstChild();
}
// ALL CASES: add sensor to its mother volume (module or ladder).
mother->addComponent(det.release());
}
|