CMS 3D CMS Logo
ME0GeometryParsFromDD.cc
Go to the documentation of this file.
1 /* Implementation of the ME0GeometryParsFromDD Class
2  * Build the ME0Geometry from the DDD and DD4Hep description
3  *
4  * DD4hep part added to the original old file (DD version) made by M. Maggi (INFN Bari)
5  * Author: Sergio Lo Meo (sergio.lo.meo@cern.ch)
6  * Created: Thu, 25 Feb 2021
7  *
8  */
9 #include "Geometry/GEMGeometryBuilder/src/ME0GeometryParsFromDD.h"
10 #include "DataFormats/MuonDetId/interface/ME0DetId.h"
11 #include "CondFormats/GeometryObjects/interface/RecoIdealGeometry.h"
12 #include "DetectorDescription/Core/interface/DDFilter.h"
13 #include "DetectorDescription/Core/interface/DDFilteredView.h"
14 #include "DetectorDescription/Core/interface/DDSolid.h"
15 #include "Geometry/MuonNumbering/interface/MuonGeometryNumbering.h"
16 #include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
17 #include "Geometry/MuonNumbering/interface/ME0NumberingScheme.h"
18 #include "Geometry/MuonNumbering/interface/MuonGeometryConstants.h"
19 #include "DataFormats/GeometryVector/interface/Basic3DVector.h"
20 #include "DetectorDescription/DDCMS/interface/DDFilteredView.h"
21 #include "DetectorDescription/DDCMS/interface/DDCompactView.h"
22 #include "DetectorDescription/DDCMS/interface/DDSpecParRegistry.h"
23 
24 #include <iostream>
25 #include <algorithm>
26 
27 // DD
28 
29 void ME0GeometryParsFromDD::build(const DDCompactView* cview,
30  const MuonGeometryConstants& muonConstants,
31  RecoIdealGeometry& rgeo) {
32  std::string attribute = "MuStructure";
33  std::string value = "MuonEndCapME0";
34 
35  // Asking only for the MuonME0's
36  DDSpecificsMatchesValueFilter filter{DDValue(attribute, value, 0.0)};
37  DDFilteredView fview(*cview, filter);
38 
39  this->buildGeometry(fview, muonConstants, rgeo);
40 }
41 
42 void ME0GeometryParsFromDD::buildGeometry(DDFilteredView& fv,
43  const MuonGeometryConstants& muonConstants,
44  RecoIdealGeometry& rgeo) {
45  LogDebug("ME0GeometryParsFromDD") << "Building the geometry service";
46  LogDebug("ME0GeometryParsFromDD") << "About to run through the ME0 structure\n"
47  << " First logical part " << fv.logicalPart().name().name();
48 
49  edm::LogVerbatim("ME0GeometryParsFromDD") << "(0) ME0GeometryParsFromDD - DDD ";
50  MuonGeometryNumbering muonDDDNumbering(muonConstants);
51  ME0NumberingScheme me0Numbering(muonConstants);
52 
53  bool doChambers = fv.firstChild();
54  LogDebug("ME0GeometryParsFromDD") << "doChamber = " << doChambers;
55  // loop over superchambers
56  while (doChambers) {
57  // getting chamber id from eta partitions
58  fv.firstChild();
59  fv.firstChild();
60  ME0DetId detIdCh =
61  ME0DetId(me0Numbering.baseNumberToUnitNumber(muonDDDNumbering.geoHistoryToBaseNumber(fv.geoHistory())));
62  // back to chambers
63  fv.parent();
64  fv.parent();
65 
66  buildChamber(fv, detIdCh, rgeo);
67 
68  // loop over chambers
69  // only 1 chamber
70  bool doLayers = fv.firstChild();
71  while (doLayers) {
72  // get layer ID
73  fv.firstChild();
74  ME0DetId detIdLa =
75  ME0DetId(me0Numbering.baseNumberToUnitNumber(muonDDDNumbering.geoHistoryToBaseNumber(fv.geoHistory())));
76  fv.parent();
77  // build layer
78  buildLayer(fv, detIdLa, rgeo);
79 
80  // loop over ME0EtaPartitions
81  bool doEtaPart = fv.firstChild();
82  while (doEtaPart) {
83  ME0DetId detId =
84  ME0DetId(me0Numbering.baseNumberToUnitNumber(muonDDDNumbering.geoHistoryToBaseNumber(fv.geoHistory())));
85  buildEtaPartition(fv, detId, rgeo);
86 
87  doEtaPart = fv.nextSibling();
88  }
89  fv.parent();
90  doLayers = fv.nextSibling();
91  }
92  fv.parent();
93  doChambers = fv.nextSibling();
94  }
95 }
96 
97 void ME0GeometryParsFromDD::buildChamber(DDFilteredView& fv, ME0DetId detId, RecoIdealGeometry& rgeo) {
98  LogDebug("ME0GeometryParsFromDD") << "buildChamber " << fv.logicalPart().name().name() << " " << detId << std::endl;
99 
100  std::vector<double> pars = getDimension(fv);
101  std::vector<double> vtra = getTranslation(fv);
102  std::vector<double> vrot = getRotation(fv);
103  edm::LogVerbatim("ME0GeometryParsFromDD")
104  << "(4) DDD, Chamber DetID " << detId.chamberId().rawId() << " Name " << fv.logicalPart().name().name();
105 
106  rgeo.insert(detId.chamberId().rawId(), vtra, vrot, pars, {fv.logicalPart().name().name()});
107 }
108 
109 void ME0GeometryParsFromDD::buildLayer(DDFilteredView& fv, ME0DetId detId, RecoIdealGeometry& rgeo) {
110  LogDebug("ME0GeometryParsFromDD") << "buildLayer " << fv.logicalPart().name().name() << " " << detId << std::endl;
111 
112  std::vector<double> pars = getDimension(fv);
113  std::vector<double> vtra = getTranslation(fv);
114  std::vector<double> vrot = getRotation(fv);
115 
116  edm::LogVerbatim("ME0GeometryParsFromDD")
117  << "(5) DDD, Layer DetID " << detId.layerId().rawId() << " Name " << fv.logicalPart().name().name();
118  rgeo.insert(detId.layerId().rawId(), vtra, vrot, pars, {fv.logicalPart().name().name()});
119 }
120 
121 void ME0GeometryParsFromDD::buildEtaPartition(DDFilteredView& fv, ME0DetId detId, RecoIdealGeometry& rgeo) {
122  LogDebug("ME0GeometryParsFromDD") << "buildEtaPartition " << fv.logicalPart().name().name() << " " << detId
123  << std::endl;
124 
125  // EtaPartition specific parameter (nstrips and npads)
126  DDValue numbOfStrips("nStrips");
127  DDValue numbOfPads("nPads");
128  const std::vector<const DDsvalues_type*>& specs = fv.specifics();
129  double nStrips = 0., nPads = 0.;
130  for (auto const& is : specs) {
131  if (DDfetch(is, numbOfStrips))
132  nStrips = numbOfStrips.doubles()[0];
133  if (DDfetch(is, numbOfPads))
134  nPads = numbOfPads.doubles()[0];
135  }
136  LogDebug("ME0GeometryParsFromDD") << ((nStrips == 0.) ? ("No nStrips found!!")
137  : ("Number of strips: " + std::to_string(nStrips)));
138  LogDebug("ME0GeometryParsFromDD") << ((nPads == 0.) ? ("No nPads found!!")
139  : ("Number of pads: " + std::to_string(nPads)));
140 
141  std::vector<double> pars = getDimension(fv);
142  pars.emplace_back(nStrips);
143  pars.emplace_back(nPads);
144  std::vector<double> vtra = getTranslation(fv);
145  std::vector<double> vrot = getRotation(fv);
146 
147  edm::LogVerbatim("ME0GeometryParsFromDD")
148  << "(6) DDD, Eta Partion DetID " << detId.rawId() << " Name " << fv.logicalPart().name().name() << " nStrips "
149  << nStrips << " nPads " << nPads;
150 
151  rgeo.insert(detId.rawId(), vtra, vrot, pars, {fv.logicalPart().name().name()});
152 }
153 
154 std::vector<double> ME0GeometryParsFromDD::getDimension(DDFilteredView& fv) {
155  std::vector<double> dpar = fv.logicalPart().solid().parameters();
156  //dpar[4] bottom width is along local X
157  //dpar[8] top width is along local X
158  //dpar[3] thickness is long local Z
159  //dpar[0] length is along local Y
160  LogDebug("ME0GeometryParsFromDD") << "dimension dx1 " << dpar[4] << ", dx2 " << dpar[8] << ", dy " << dpar[0]
161  << ", dz " << dpar[3];
162  edm::LogVerbatim("ME0GeometryParsFromDD")
163  << "(1) DDD, dimension dx1 " << dpar[4] << ", dx2 " << dpar[8] << ", dy " << dpar[0] << ", dz " << dpar[3];
164  return {dpar[4], dpar[8], dpar[0], dpar[3]};
165 }
166 
167 std::vector<double> ME0GeometryParsFromDD::getTranslation(DDFilteredView& fv) {
168  const DDTranslation& tran = fv.translation();
169  edm::LogVerbatim("ME0GeometryParsFromDD")
170  << "(2) DDD, tran vector " << tran.x() << " " << tran.y() << " " << tran.z();
171  return {tran.x(), tran.y(), tran.z()};
172 }
173 
174 std::vector<double> ME0GeometryParsFromDD::getRotation(DDFilteredView& fv) {
175  const DDRotationMatrix& rota = fv.rotation(); //.Inverse();
176  DD3Vector x, y, z;
177  rota.GetComponents(x, y, z);
178  edm::LogVerbatim("ME0GeometryParsFromDD")
179  << "(3) DDD, rot matrix " << x.X() << " " << x.Y() << " " << x.Z() << " " << y.X() << " " << y.Y() << " "
180  << y.Z() << " " << z.X() << " " << z.Y() << " " << z.Z();
181  return {x.X(), x.Y(), x.Z(), y.X(), y.Y(), y.Z(), z.X(), z.Y(), z.Z()};
182 }
183 
184 // DD4HEP
185 
186 void ME0GeometryParsFromDD::build(const cms::DDCompactView* cview,
187  const MuonGeometryConstants& muonConstants,
188  RecoIdealGeometry& rgeo) {
189  std::string attribute = "MuStructure";
190  std::string value = "MuonEndCapME0";
191  const cms::DDFilter filter(attribute, value);
192  cms::DDFilteredView fview(*cview, filter);
193  this->buildGeometry(fview, muonConstants, rgeo);
194 }
195 
196 void ME0GeometryParsFromDD::buildGeometry(cms::DDFilteredView& fv,
197  const MuonGeometryConstants& muonConstants,
198  RecoIdealGeometry& rgeo) {
199  edm::LogVerbatim("ME0GeometryParsFromDD") << "(0) ME0GeometryParsFromDD - DD4HEP ";
200 
201  MuonGeometryNumbering mdddnum(muonConstants);
202  ME0NumberingScheme me0Num(muonConstants);
203 
204  static constexpr uint32_t levelChamber = 7;
205  static constexpr uint32_t levelLayer = 8;
206  uint32_t theLevelPart = muonConstants.getValue("level");
207  uint32_t theSectorLevel = muonConstants.getValue("m0_sector") / theLevelPart;
208 
209  while (fv.firstChild()) {
210  const auto& history = fv.history();
211  MuonBaseNumber num(mdddnum.geoHistoryToBaseNumber(history));
212  ME0DetId detId(me0Num.baseNumberToUnitNumber(num));
213 
214  if (fv.level() == levelChamber) {
215  buildChamber(fv, detId, rgeo);
216  } else if (fv.level() == levelLayer) {
217  buildLayer(fv, detId, rgeo);
218  } else if (history.tags.size() > theSectorLevel) {
219  buildEtaPartition(fv, detId, rgeo);
220  }
221  } // end while
222 } // end buildGeometry
223 
224 void ME0GeometryParsFromDD::buildChamber(cms::DDFilteredView& fv, ME0DetId detId, RecoIdealGeometry& rgeo) {
225  std::string_view name = fv.name();
226  std::vector<double> pars = getDimension(fv);
227  std::vector<double> vtra = getTranslation(fv);
228  std::vector<double> vrot = getRotation(fv);
229 
230  edm::LogVerbatim("ME0GeometryParsFromDD")
231  << "(4) DD4HEP, Chamber DetID " << detId.chamberId().rawId() << " Name " << std::string(name);
232 
233  rgeo.insert(detId.chamberId().rawId(), vtra, vrot, pars, {std::string(name)});
234 }
235 
236 void ME0GeometryParsFromDD::buildLayer(cms::DDFilteredView& fv, ME0DetId detId, RecoIdealGeometry& rgeo) {
237  std::string_view name = fv.name();
238  std::vector<double> pars = getDimension(fv);
239  std::vector<double> vtra = getTranslation(fv);
240  std::vector<double> vrot = getRotation(fv);
241 
242  edm::LogVerbatim("ME0GeometryParsFromDD")
243  << "(5) DD4HEP, Layer DetID " << detId.layerId().rawId() << " Name " << std::string(name);
244  rgeo.insert(detId.layerId().rawId(), vtra, vrot, pars, {std::string(name)});
245 }
246 
247 void ME0GeometryParsFromDD::buildEtaPartition(cms::DDFilteredView& fv, ME0DetId detId, RecoIdealGeometry& rgeo) {
248  auto nStrips = fv.get<double>("nStrips");
249  auto nPads = fv.get<double>("nPads");
250  std::string_view name = fv.name();
251  std::vector<double> pars = getDimension(fv);
252  pars.emplace_back(nStrips);
253  pars.emplace_back(nPads);
254  std::vector<double> vtra = getTranslation(fv);
255  std::vector<double> vrot = getRotation(fv);
256 
257  edm::LogVerbatim("ME0GeometryParsFromDD") << "(6) DD4HEP, Eta Partion DetID " << detId.rawId() << " Name "
258  << std::string(name) << " nStrips " << nStrips << " nPads " << nPads;
259 
260  rgeo.insert(detId.rawId(), vtra, vrot, pars, {std::string(name)});
261 }
262 
263 std::vector<double> ME0GeometryParsFromDD::getDimension(cms::DDFilteredView& fv) {
264  std::vector<double> dpar = fv.parameters();
265 
266  edm::LogVerbatim("ME0GeometryParsFromDD")
267  << "(1) DD4HEP, dimension dx1 " << dpar[0] / dd4hep::mm << ", dx2 " << dpar[1] / dd4hep::mm << ", dy "
268  << dpar[3] / dd4hep::mm << ", dz " << dpar[2] / dd4hep::mm;
269 
270  return {dpar[0] / dd4hep::mm, dpar[1] / dd4hep::mm, dpar[3] / dd4hep::mm, dpar[2] / dd4hep::mm};
271 }
272 
273 std::vector<double> ME0GeometryParsFromDD::getTranslation(cms::DDFilteredView& fv) {
274  std::vector<double> tran(3);
275  tran[0] = static_cast<double>(fv.translation().X()) / dd4hep::mm;
276  tran[1] = static_cast<double>(fv.translation().Y()) / dd4hep::mm;
277  tran[2] = static_cast<double>(fv.translation().Z()) / dd4hep::mm;
278 
279  edm::LogVerbatim("ME0GeometryParsFromDD")
280  << "(2) DD4HEP, tran vector " << tran[0] << " " << tran[1] << " " << tran[2];
281  return {tran[0], tran[1], tran[2]};
282 }
283 
284 std::vector<double> ME0GeometryParsFromDD::getRotation(cms::DDFilteredView& fv) {
285  DDRotationMatrix rota;
286  fv.rot(rota);
287  DD3Vector x, y, z;
288  rota.GetComponents(x, y, z);
289  const std::vector<double> rot = {x.X(), x.Y(), x.Z(), y.X(), y.Y(), y.Z(), z.X(), z.Y(), z.Z()};
290  edm::LogVerbatim("ME0GeometryParsFromDD")
291  << "(3) DD4HEP, rot matrix " << rot[0] << " " << rot[1] << " " << rot[2] << " " << rot[3] << " " << rot[4]
292  << " " << rot[5] << " " << rot[6] << " " << rot[7] << " " << rot[8];
293  return {rot[0], rot[1], rot[2], rot[3], rot[4], rot[5], rot[6], rot[7], rot[8]};
294 }
cms::DDFilteredView::rot
const Double_t * rot() const
The absolute rotation of the current node.
Definition: DDFilteredView.cc:120
ME0GeometryParsFromDD::buildEtaPartition
void buildEtaPartition(DDFilteredView &fv, ME0DetId detId, RecoIdealGeometry &rgeo)
Definition: ME0GeometryParsFromDD.cc:121
MuonGeometryConstants
Definition: MuonGeometryConstants.h:20
DDAxes::y
cms::DDFilteredView::parameters
const std::vector< double > parameters() const
extract shape parameters
Definition: DDFilteredView.cc:536
MuonGeometryNumbering::geoHistoryToBaseNumber
MuonBaseNumber geoHistoryToBaseNumber(const DDGeoHistory &history) const
Definition: MuonGeometryNumbering.cc:38
ME0GeometryParsFromDD::buildLayer
void buildLayer(DDFilteredView &fv, ME0DetId detId, RecoIdealGeometry &rgeo)
Definition: ME0GeometryParsFromDD.cc:109
ME0NumberingScheme::baseNumberToUnitNumber
int baseNumberToUnitNumber(const MuonBaseNumber &) const override
Definition: ME0NumberingScheme.cc:30
DDFilteredView::logicalPart
const DDLogicalPart & logicalPart() const
The logical-part of the current node in the filtered-view.
Definition: DDFilteredView.cc:16
ME0DetId.h
MuonGeometryConstants.h
HistogramManager_cfi.specs
specs
Definition: HistogramManager_cfi.py:83
ME0GeometryParsFromDD::getDimension
std::vector< double > getDimension(DDFilteredView &fv)
Definition: ME0GeometryParsFromDD.cc:154
Basic3DVector.h
ME0GeometryParsFromDD::buildGeometry
void buildGeometry(DDFilteredView &, const MuonGeometryConstants &, RecoIdealGeometry &)
Definition: ME0GeometryParsFromDD.cc:42
MuonGeometryNumbering.h
DDAxes::x
cms::DDFilteredView
Definition: DDFilteredView.h:70
DDFilteredView::parent
bool parent()
set the current node to the parent node ...
Definition: DDFilteredView.cc:161
ME0DetId::layerId
ME0DetId layerId() const
Return the corresponding LayerId (mask eta partition)
Definition: ME0DetId.h:55
cms::DDFilteredView::name
std::string_view name() const
Definition: DDFilteredView.cc:854
cms::DDFilteredView::translation
const Translation translation() const
Definition: DDFilteredView.cc:103
cms::DDFilteredView::get
T get(const std::string &)
extract attribute value
DDFilteredView::firstChild
bool firstChild()
set the current node to the first child ...
Definition: DDFilteredView.cc:86
DDFilteredView.h
cms::DDFilter
Definition: DDFilteredView.h:59
DDTranslation
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DDTranslation
Definition: DDTranslation.h:7
DDCompactView
Compact representation of the geometrical detector hierarchy.
Definition: DDCompactView.h:81
DDBase::name
const N & name() const
Definition: DDBase.h:59
DDAxes::z
cms::DDFilteredView::firstChild
bool firstChild()
set the current node to the first child
Definition: DDFilteredView.cc:268
ME0GeometryParsFromDD::buildChamber
void buildChamber(DDFilteredView &fv, ME0DetId detId, RecoIdealGeometry &rgeo)
Definition: ME0GeometryParsFromDD.cc:97
DDSolid.h
DDFilteredView::nextSibling
bool nextSibling()
set the current node to the next sibling ...
Definition: DDFilteredView.cc:124
ME0NumberingScheme.h
DDfetch
bool DDfetch(const DDsvalues_type *, DDValue &)
helper for retrieving DDValues from DDsvalues_type *.
Definition: DDsvalues.cc:79
DDFilteredView.h
DDCompactView.h
ALCARECOTkAlBeamHalo_cff.filter
filter
Definition: ALCARECOTkAlBeamHalo_cff.py:27
RecoIdealGeometry.h
ME0GeometryParsFromDD::getRotation
std::vector< double > getRotation(DDFilteredView &fv)
Definition: ME0GeometryParsFromDD.cc:174
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
DD3Vector
ROOT::Math::DisplacementVector3D< ROOT::Math::Cartesian3D< double > > DD3Vector
Definition: PGeometricDetBuilder.cc:19
DDSpecParRegistry.h
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:233
cms::DDFilteredView::level
const int level() const
get Iterator level
Definition: DDFilteredView.cc:636
ME0GeometryParsFromDD.h
MuonBaseNumber.h
cms::DDFilteredView::history
const ExpandedNodes & history()
The numbering history of the current node.
Definition: DDFilteredView.cc:683
DDFilter.h
value
Definition: value.py:1
DDName::name
const std::string & name() const
Returns the name.
Definition: DDName.cc:41
MuonGeometryConstants::getValue
int getValue(const std::string &name) const
Definition: MuonGeometryConstants.cc:8
DDFilteredView::geoHistory
const DDGeoHistory & geoHistory() const
The list of ancestors up to the root-node of the current node.
Definition: DDFilteredView.cc:30
ME0DetId
Definition: ME0DetId.h:16
EgammaValidation_cff.num
num
Definition: EgammaValidation_cff.py:34
DDFilteredView::specifics
std::vector< const DDsvalues_type * > specifics() const
Definition: DDFilteredView.cc:32
me0TriggerPseudoDigis_cff.nStrips
nStrips
1.2 is to make the matching window safely the two nearest strips 0.35 is the size of an ME0 chamber i...
Definition: me0TriggerPseudoDigis_cff.py:26
ME0GeometryParsFromDD::build
void build(const DDCompactView *, const MuonGeometryConstants &, RecoIdealGeometry &)
Definition: ME0GeometryParsFromDD.cc:29
ME0NumberingScheme
Definition: ME0NumberingScheme.h:9
ME0DetId::chamberId
ME0DetId chamberId() const
Return the corresponding ChamberId (mask layers)
Definition: ME0DetId.h:53
MuonBaseNumber
Definition: MuonBaseNumber.h:21
DetId::rawId
constexpr uint32_t rawId() const
get the raw id
Definition: DetId.h:57
DDRotationMatrix
ROOT::Math::Rotation3D DDRotationMatrix
A DDRotationMatrix is currently implemented with a ROOT Rotation3D.
Definition: DDRotationMatrix.h:8
DDValue
Definition: DDValue.h:21
edm::LogVerbatim
Log< level::Info, true > LogVerbatim
Definition: MessageLogger.h:128
DDFilteredView::translation
const DDTranslation & translation() const
The absolute translation of the current node.
Definition: DDFilteredView.cc:26
cms::DDCompactView
Definition: DDCompactView.h:31
makeMuonMisalignmentScenario.rot
rot
Definition: makeMuonMisalignmentScenario.py:322
DDSolid::parameters
const std::vector< double > & parameters(void) const
Give the parameters of the solid.
Definition: DDSolid.cc:121
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
DDSpecificsMatchesValueFilter
Definition: DDFilter.h:70
DDFilteredView::rotation
const DDRotationMatrix & rotation() const
The absolute rotation of the current node.
Definition: DDFilteredView.cc:28
MuonGeometryNumbering
Definition: MuonGeometryNumbering.h:24
DDValue::doubles
const std::vector< double > & doubles() const
a reference to the double-valued values stored in the given instance of DDValue
Definition: DDValue.cc:111
DDFilteredView
Definition: DDFilteredView.h:20
RecoIdealGeometry::insert
bool insert(DetId id, const std::vector< double > &trans, const std::vector< double > &rot, const std::vector< double > &pars)
Definition: RecoIdealGeometry.h:33
DDLogicalPart::solid
const DDSolid & solid(void) const
Returns a reference object of the solid being the shape of this LogicalPart.
Definition: DDLogicalPart.cc:120
RecoIdealGeometry
Definition: RecoIdealGeometry.h:28
ME0GeometryParsFromDD::getTranslation
std::vector< double > getTranslation(DDFilteredView &fv)
Definition: ME0GeometryParsFromDD.cc:167
Generated for CMSSW Reference Manual by doxygen 1.8.16