CMS 3D CMS Logo
DDCmsMTDConstruction.cc
Go to the documentation of this file.
1 #include "Geometry/MTDNumberingBuilder/plugins/DDCmsMTDConstruction.h"
2 
3 #include <utility>
4 #include <sstream>
5 
6 #include "DetectorDescription/Core/interface/DDFilteredView.h"
7 #include "DetectorDescription/Core/interface/DDCompactView.h"
8 #include "DetectorDescription/DDCMS/interface/DDFilteredView.h"
9 #include "DetectorDescription/DDCMS/interface/DDCompactView.h"
10 #include "DetectorDescription/DDCMS/interface/DDSpecParRegistry.h"
11 #include "Geometry/MTDNumberingBuilder/interface/GeometricTimingDet.h"
12 #include "Geometry/MTDNumberingBuilder/plugins/ExtractStringFromDD.h"
13 #include "Geometry/MTDNumberingBuilder/plugins/CmsMTDConstruction.h"
14 
15 #include "DataFormats/ForwardDetId/interface/BTLDetId.h"
16 #include "DataFormats/ForwardDetId/interface/ETLDetId.h"
17 
18 //#define EDM_ML_DEBUG
19 
20 class DDNameFilter : public DDFilter {
21 public:
22  void add(const std::string& add) { allowed_.emplace_back(add); }
23  void veto(const std::string& veto) { veto_.emplace_back(veto); }
24 
25  bool accept(const DDExpandedView& ev) const final {
26  for (const auto& test : veto_) {
27  if (ev.logicalPart().name().fullname().find(test) != std::string::npos)
28  return false;
29  }
30  for (const auto& test : allowed_) {
31  if (ev.logicalPart().name().fullname().find(test) != std::string::npos)
32  return true;
33  }
34  return false;
35  }
36 
37 private:
38  std::vector<std::string> allowed_;
39  std::vector<std::string> veto_;
40 };
41 
42 std::unique_ptr<GeometricTimingDet> DDCmsMTDConstruction::construct(const DDCompactView& cpv) {
43  std::string attribute{"CMSCutsRegion"};
44  DDNameFilter filter;
45  filter.add("mtd:");
46  filter.add("btl:");
47  filter.add("etl:");
48  filter.veto("service");
49  filter.veto("support");
50  filter.veto("FSide");
51  filter.veto("BSide");
52  filter.veto("LSide");
53  filter.veto("RSide");
54  filter.veto("Between");
55  filter.veto("SupportPlate");
56  filter.veto("Shield");
57  filter.veto("ThermalScreen");
58  filter.veto("Aluminium_Disc");
59  filter.veto("MIC6_Aluminium_Disc");
60  filter.veto("ThermalPad");
61  filter.veto("AlN");
62  filter.veto("LairdFilm");
63  filter.veto("ETROC");
64  filter.veto("SensorModule");
65  filter.veto("DiscSector");
66  filter.veto("LGAD_Substrate");
67 
68  DDFilteredView fv(cpv, filter);
69 
70  CmsMTDStringToEnum theCmsMTDStringToEnum;
71 
72  auto check_root = theCmsMTDStringToEnum.type(ExtractStringFromDD<DDFilteredView>::getString(attribute, &fv));
73  if (check_root != GeometricTimingDet::MTD) {
74  fv.firstChild();
75  auto check_child = theCmsMTDStringToEnum.type(ExtractStringFromDD<DDFilteredView>::getString(attribute, &fv));
76  if (check_child != GeometricTimingDet::MTD) {
77  throw cms::Exception("Configuration") << " The first child of the DDFilteredView is not what is expected \n"
78  << ExtractStringFromDD<DDFilteredView>::getString(attribute, &fv) << "\n";
79  }
80  fv.parent();
81  }
82 
83 #ifdef EDM_ML_DEBUG
84  edm::LogInfo("MTDNumbering") << "Top level node = " << fv.name();
85 #endif
86 
87  auto mtd = std::make_unique<GeometricTimingDet>(&fv, GeometricTimingDet::MTD);
88  size_t limit = 0;
89  CmsMTDConstruction<DDFilteredView> theCmsMTDConstruction;
90 
91  std::vector<GeometricTimingDet*> subdet;
92  std::vector<GeometricTimingDet*> layer;
93 
94  do {
95  GeometricTimingDet::GeometricTimingEnumType fullNode = theCmsMTDStringToEnum.type(fv.name());
96  GeometricTimingDet::GeometricTimingEnumType thisNode =
97  theCmsMTDStringToEnum.type(fv.name().substr(0, CmsMTDStringToEnum::kModStrLen));
98  size_t num = fv.geoHistory().size();
99 
100 #ifdef EDM_ML_DEBUG
101  edm::LogVerbatim("MTDNumbering") << "Module = " << fv.name() << " fullNode = " << fullNode
102  << " thisNode = " << thisNode;
103 #endif
104 
105  if (fullNode == GeometricTimingDet::BTL || fullNode == GeometricTimingDet::ETL) {
106  limit = 0;
107 
108  // define subdetectors as GeometricTimingDet components
109 
110  subdet.emplace_back(theCmsMTDConstruction.buildSubdet(fv));
111  }
112  if (fullNode == GeometricTimingDet::BTLLayer || fullNode == GeometricTimingDet::ETLDisc) {
113  layer.emplace_back(theCmsMTDConstruction.buildLayer(fv));
114 #ifdef EDM_ML_DEBUG
115  edm::LogVerbatim("DD4hep_MTDNumbering") << "Number of layers: " << layer.size();
116 #endif
117  }
118  //
119  // the level chosen corresponds to wafers for D50 and previous scenarios
120  //
121  if ((thisNode == GeometricTimingDet::BTLModule) && limit == 0) {
122  limit = num + 1;
123  }
124  //
125  // workaround to make old and TDR structure to cohexist until needed
126  //
127  else if ((thisNode == GeometricTimingDet::ETLModule) && limit == 0) {
128  if (theCmsMTDConstruction.isETLtdr(fv)) {
129  limit = num;
130  } else {
131  limit = num + 1;
132  }
133  }
134  if (num != limit && limit > 0) {
135  continue;
136  }
137  if (thisNode == GeometricTimingDet::BTLModule) {
138 #ifdef EDM_ML_DEBUG
139  edm::LogVerbatim("MTDNumbering") << "Registered in GeometricTimingDet as type " << thisNode;
140 #endif
141  theCmsMTDConstruction.buildBTLModule(fv, layer.back());
142  limit = num;
143  } else if (thisNode == GeometricTimingDet::ETLModule) {
144 #ifdef EDM_ML_DEBUG
145  edm::LogVerbatim("MTDNumbering") << "Registered in GeometricTimingDet as type " << thisNode;
146 #endif
147  theCmsMTDConstruction.buildETLModule(fv, layer.back());
148  limit = num;
149  }
150  } while (fv.next());
151 
152  // sort GeometricTimingDet
153 
154 #ifdef EDM_ML_DEBUG
155  auto comp = mtd->deepComponents();
156  std::stringstream before(std::stringstream::in | std::stringstream::out);
157  for (const auto& it : comp) {
158  before << "ORDER1 " << it->geographicalId().rawId() << " " << it->type() << " " << it->translation().z() << "\n";
159  }
160  edm::LogVerbatim("MTDNumbering") << "GeometricTimingDet order before sorting \n" << before.str();
161 #endif
162 
163  for (size_t index = 0; index < layer.size(); index++) {
164  GeometricTimingDet::ConstGeometricTimingDetContainer& icomp = layer[index]->components();
165  std::stable_sort(icomp.begin(), icomp.end(), CmsMTDConstruction<DDFilteredView>::mtdOrderZ);
166  std::stable_sort(icomp.begin(), icomp.end(), CmsMTDConstruction<DDFilteredView>::mtdOrderRR);
167  if (index > 0) {
168  std::stable_sort(icomp.begin(), icomp.end(), CmsMTDConstruction<DDFilteredView>::mtdOrderPhi);
169  }
170  }
171 
172  // Add layers to subdetectors:
173  // first BTL (one layer only)
174 
175  subdet[0]->addComponent(layer[0]);
176 
177  // then ETL (number of layers depend on preTDR design or not)
178 
179  if (layer.size() == kNLayerPreTDR) {
180  subdet[1]->addComponent(layer[1]);
181  subdet[2]->addComponent(layer[2]);
182  } else if (layer.size() == kNLayerTDR) {
183  subdet[1]->addComponent(layer[1]);
184  subdet[1]->addComponent(layer[2]);
185  subdet[2]->addComponent(layer[3]);
186  subdet[2]->addComponent(layer[4]);
187  } else {
188  throw cms::Exception("DD4hep_MTDNumbering") << "Wrong number of layers: " << layer.size();
189  }
190 
191  // Add subdetectors to MTD
192 
193  mtd.get()->addComponents(subdet);
194 
195 #ifdef EDM_ML_DEBUG
196  comp.clear();
197  comp = mtd->deepComponents();
198  std::stringstream after(std::stringstream::in | std::stringstream::out);
199  for (const auto& it : comp) {
200  after << "ORDER2 " << it->geographicalId().rawId() << " " << it->type() << " " << it->translation().z() << "\n";
201  }
202  edm::LogVerbatim("MTDNumbering") << "GeometricTimingDet order after sorting \n" << after.str();
203 #endif
204 
205  return mtd;
206 }
207 
208 std::unique_ptr<GeometricTimingDet> DDCmsMTDConstruction::construct(const cms::DDCompactView& cpv) {
209  cms::DDFilteredView fv(cpv.detector(), cpv.detector()->worldVolume());
210 
211  fv.next(0);
212  edm::LogVerbatim("MTDNumbering") << fv.path();
213  auto mtd = std::make_unique<GeometricTimingDet>(&fv, GeometricTimingDet::MTD);
214 
215  cms::DDSpecParRefs ref;
216  const cms::DDSpecParRegistry& mypar = cpv.specpars();
217  std::string attribute("MtdDDStructure");
218  mypar.filter(ref, attribute, "BarrelTimingLayer");
219  mypar.filter(ref, attribute, "EndcapTimingLayer");
220  fv.mergedSpecifics(ref);
221 
222 #ifdef EDM_ML_DEBUG
223  edm::LogVerbatim("MTDNumbering") << "Active filters using " << attribute << ":";
224  fv.printFilter();
225  edm::LogVerbatim("Geometry").log([&ref](auto& log) {
226  log << "Filtered DD SpecPar Registry size: " << ref.size() << "\n";
227  for (const auto& t : ref) {
228  log << "\nRegExps { ";
229  for (const auto& ki : t->paths)
230  log << ki << " ";
231  log << "};\n ";
232  for (const auto& kl : t->spars) {
233  log << kl.first << " = ";
234  for (const auto& kil : kl.second) {
235  log << kil << " ";
236  }
237  log << "\n ";
238  }
239  }
240  });
241 #endif
242 
243  bool doSubdet = fv.firstChild();
244  edm::LogVerbatim("MTDNumbering") << fv.path();
245 
246  CmsMTDStringToEnum theCmsMTDStringToEnum;
247 
248  CmsMTDConstruction<cms::DDFilteredView> theCmsMTDConstruction;
249 
250  std::vector<GeometricTimingDet*> subdet;
251  std::vector<GeometricTimingDet*> layer;
252 
253  while (doSubdet) {
254  std::string nodeName(fv.name());
255  GeometricTimingDet::GeometricTimingEnumType fullNode = theCmsMTDStringToEnum.type(nodeName);
256  GeometricTimingDet::GeometricTimingEnumType thisNode =
257  theCmsMTDStringToEnum.type(nodeName.substr(0, CmsMTDStringToEnum::kModStrLen));
258 
259 #ifdef EDM_ML_DEBUG
260  edm::LogVerbatim("MTDNumbering") << fv.path();
261  edm::LogVerbatim("MTDNumbering") << "Module = " << fv.name() << " fullNode = " << fullNode
262  << " thisNode = " << thisNode;
263 #endif
264 
265  if (fullNode == GeometricTimingDet::BTL || fullNode == GeometricTimingDet::ETL) {
266  // define subdetectors as GeometricTimingDet components
267 
268  subdet.emplace_back(theCmsMTDConstruction.buildSubdet(fv));
269  }
270  if (fullNode == GeometricTimingDet::BTLLayer || fullNode == GeometricTimingDet::ETLDisc) {
271  layer.emplace_back(theCmsMTDConstruction.buildLayer(fv));
272 #ifdef EDM_ML_DEBUG
273  edm::LogVerbatim("DD4hep_MTDNumbering") << "Number of layers: " << layer.size();
274 #endif
275  }
276  if (thisNode == GeometricTimingDet::BTLModule) {
277 #ifdef EDM_ML_DEBUG
278  edm::LogVerbatim("DD4hep_MTDNumbering") << "Registered in GeometricTimingDet as type " << thisNode;
279 #endif
280  theCmsMTDConstruction.buildBTLModule(fv, layer.back());
281  } else if (thisNode == GeometricTimingDet::ETLModule) {
282 #ifdef EDM_ML_DEBUG
283  edm::LogVerbatim("DD4hep_MTDNumbering") << "Registered in GeometricTimingDet as type " << thisNode;
284 #endif
285  theCmsMTDConstruction.buildETLModule(fv, layer.back());
286  }
287 
288  doSubdet = fv.firstChild();
289  }
290 
291  // sort GeometricTimingDet
292 
293 #ifdef EDM_ML_DEBUG
294  auto comp = mtd->deepComponents();
295  std::stringstream before(std::stringstream::in | std::stringstream::out);
296  for (const auto& it : comp) {
297  before << "ORDER1 " << it->geographicalId().rawId() << " " << it->type() << " " << it->translation().z() << "\n";
298  }
299  edm::LogVerbatim("DD4hep_MTDNumbering") << "GeometricTimingDet order before sorting \n" << before.str();
300 #endif
301 
302  for (size_t index = 0; index < layer.size(); index++) {
303  GeometricTimingDet::ConstGeometricTimingDetContainer& icomp = layer[index]->components();
304  std::stable_sort(icomp.begin(), icomp.end(), CmsMTDConstruction<cms::DDFilteredView>::mtdOrderZ);
305  std::stable_sort(icomp.begin(), icomp.end(), CmsMTDConstruction<cms::DDFilteredView>::mtdOrderRR);
306  if (index > 0) {
307  std::stable_sort(icomp.begin(), icomp.end(), CmsMTDConstruction<cms::DDFilteredView>::mtdOrderPhi);
308  }
309  }
310 
311  // Add layers to subdetectors:
312  // first BTL (one layer only)
313 
314  subdet[0]->addComponent(layer[0]);
315 
316  // then ETL (number of layers depend on preTDR design or not)
317 
318  if (layer.size() == kNLayerPreTDR) {
319  subdet[1]->addComponent(layer[1]);
320  subdet[2]->addComponent(layer[2]);
321  } else if (layer.size() == kNLayerTDR) {
322  subdet[1]->addComponent(layer[1]);
323  subdet[1]->addComponent(layer[2]);
324  subdet[2]->addComponent(layer[3]);
325  subdet[2]->addComponent(layer[4]);
326  } else {
327  throw cms::Exception("DD4hep_MTDNumbering") << "Wrong number of layers: " << layer.size();
328  }
329 
330  // Add subdetectors to MTD
331 
332  mtd.get()->addComponents(subdet);
333 
334 #ifdef EDM_ML_DEBUG
335  comp.clear();
336  comp = mtd->deepComponents();
337  std::stringstream after(std::stringstream::in | std::stringstream::out);
338  for (const auto& it : comp) {
339  after << "ORDER2 " << it->geographicalId().rawId() << " " << it->type() << " " << it->translation().z() << "\n";
340  }
341  edm::LogVerbatim("DD4hep_MTDNumbering") << "GeometricTimingDet order after sorting \n" << after.str();
342 #endif
343 
344  return mtd;
345 }
BTLDetId.h
DDNameFilter::veto
void veto(const std::string &veto)
Definition: DDCmsMTDConstruction.cc:23
DDCmsMTDConstruction::construct
static std::unique_ptr< GeometricTimingDet > construct(const DDCompactView &cpv)
Definition: DDCmsMTDConstruction.cc:42
CmsMTDStringToEnum::kModStrLen
static constexpr size_t kModStrLen
Definition: CmsMTDStringToEnum.h:12
CmsMTDConstruction::buildETLModule
void buildETLModule(FilteredView &, GeometricTimingDet *)
edm::LogInfo
Definition: MessageLogger.h:254
cms::DDSpecParRegistry
Definition: DDSpecParRegistry.h:32
CmsMTDConstruction::buildBTLModule
void buildBTLModule(FilteredView &, GeometricTimingDet *)
GeometricTimingDet::ConstGeometricTimingDetContainer
std::vector< GeometricTimingDet const * > ConstGeometricTimingDetContainer
Definition: GeometricTimingDet.h:31
DDNameFilter::allowed_
std::vector< std::string > allowed_
Definition: DDCmsMTDConstruction.cc:38
ExtractStringFromDD
Definition: ExtractStringFromDD.h:12
GeometricTimingDet::BTLModule
Definition: GeometricTimingDet.h:56
cms::DDFilteredView
Definition: DDFilteredView.h:65
DDFilteredView::parent
bool parent()
set the current node to the parent node ...
Definition: DDFilteredView.cc:161
ETLDetId.h
DDCompactView.h
AlCaHLTBitMon_QueryRunRegistry.comp
comp
Definition: AlCaHLTBitMon_QueryRunRegistry.py:249
ExtractStringFromDD.h
test
Definition: SmallWORMDict.h:13
DDNameFilter::add
void add(const std::string &add)
Definition: DDCmsMTDConstruction.cc:22
DDFilteredView::firstChild
bool firstChild()
set the current node to the first child ...
Definition: DDFilteredView.cc:86
GeometricTimingDet::ETLDisc
Definition: GeometricTimingDet.h:60
DDFilteredView.h
DDFilter
A Filter accepts or rejects a DDExpandedNode based on a user-coded decision rule.
Definition: DDFilter.h:15
DDCmsMTDConstruction::kNLayerPreTDR
static constexpr size_t kNLayerPreTDR
Definition: DDCmsMTDConstruction.h:29
DDNameFilter::accept
bool accept(const DDExpandedView &ev) const final
true, if the DDExpandedNode fulfills the filter criteria
Definition: DDCmsMTDConstruction.cc:25
DDFilteredView::next
bool next()
set current node to the next node in the filtered tree
Definition: DDFilteredView.cc:67
DDCompactView
Compact representation of the geometrical detector hierarchy.
Definition: DDCompactView.h:80
DDCmsMTDConstruction.h
DDFilteredView.h
DDCompactView.h
DDExpandedView
Provides an exploded view of the detector (tree-view)
Definition: DDExpandedView.h:41
CmsMTDConstruction::buildLayer
GeometricTimingDet * buildLayer(FilteredView &)
Definition: CmsMTDConstruction.cc:196
OrderedSet.t
t
Definition: OrderedSet.py:90
ALCARECOTkAlBeamHalo_cff.filter
filter
Definition: ALCARECOTkAlBeamHalo_cff.py:27
cms::DDSpecParRefs
std::vector< const DDSpecPar * > DDSpecParRefs
Definition: DDSpecParRegistry.h:30
DDNameFilter::veto_
std::vector< std::string > veto_
Definition: DDCmsMTDConstruction.cc:39
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
DDCmsMTDConstruction::kNLayerTDR
static constexpr size_t kNLayerTDR
Definition: DDCmsMTDConstruction.h:30
DDSpecParRegistry.h
GeometricTimingDet::ETL
Definition: GeometricTimingDet.h:59
GeometricTimingDet::MTD
Definition: GeometricTimingDet.h:52
CmsMTDConstruction::isETLtdr
bool isETLtdr(FilteredView &)
Definition: CmsMTDConstruction.cc:224
recoMuon::in
Definition: RecoMuonEnumerators.h:6
CmsMTDStringToEnum::type
GeometricTimingDet::GeometricTimingEnumType type(std::string const &) const
Definition: CmsMTDStringToEnum.cc:18
edm::LogVerbatim::log
LogVerbatim & log(F &&iF)
Definition: MessageLogger.h:327
cms::DDFilteredView::next
bool next(int)
set current node to the next node in the filtered tree
Definition: DDFilteredView.cc:407
edm::LogVerbatim
Definition: MessageLogger.h:297
DDFilteredView::geoHistory
const DDGeoHistory & geoHistory() const
The list of ancestors up to the root-node of the current node.
Definition: DDFilteredView.cc:30
EgammaValidation_cff.num
num
Definition: EgammaValidation_cff.py:34
cms::DDCompactView::specpars
DDSpecParRegistry const & specpars() const
Definition: DDCompactView.h:33
CmsMTDStringToEnum
Definition: CmsMTDStringToEnum.h:10
ExtractStringFromDD::getString
static std::string getString(const std::string &, FilteredView *)
GeometricTimingDet::BTLLayer
Definition: GeometricTimingDet.h:54
GeometricTimingDet.h
remoteMonitoring_LED_IterMethod_cfg.limit
limit
Definition: remoteMonitoring_LED_IterMethod_cfg.py:427
cms::DDSpecParRegistry::filter
void filter(DDSpecParRefs &, const std::string &, const std::string &="") const
Definition: DDSpecparRegistry.cc:98
ev
bool ev
Definition: Hydjet2Hadronizer.cc:95
Exception
Definition: hltDiff.cc:246
GeometricTimingDet::BTL
Definition: GeometricTimingDet.h:53
DDNameFilter
Definition: DDCmsMTDConstruction.cc:20
cms::DDCompactView
Definition: DDCompactView.h:29
dqm-mbProfile.log
log
Definition: dqm-mbProfile.py:17
AlignmentPI::index
index
Definition: AlignmentPayloadInspectorHelper.h:46
DDFilteredView
Definition: DDFilteredView.h:20
MillePedeFileConverter_cfg.out
out
Definition: MillePedeFileConverter_cfg.py:31
CmsMTDConstruction.h
mtd_digitizer
Definition: MTDDigitizer.h:35
GeometricTimingDet::GeometricTimingEnumType
enum GeometricTimingDet::GTDEnumType GeometricTimingEnumType
CmsMTDConstruction
Definition: CmsMTDConstruction.h:13
CmsMTDConstruction::buildSubdet
GeometricTimingDet * buildSubdet(FilteredView &)
Definition: CmsMTDConstruction.cc:178
GeometricTimingDet::ETLModule
Definition: GeometricTimingDet.h:62
cms::DDDetector::worldVolume
Volume worldVolume() const
Handle to the world volume containing everything.
Definition: DDDetector.cc:41
cms::DDCompactView::detector
const cms::DDDetector * detector() const
Definition: DDCompactView.h:32
DDFilteredView::name
const std::string & name() const
The name of a logical-part of the current node in the filtered-view.
Definition: DDFilteredView.cc:18
Generated for CMSSW Reference Manual by doxygen 1.8.16