File indexing completed on 2024-08-13 05:00:36
0001
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018
0019
0020 #include "CondFormats/GeometryObjects/interface/RecoIdealGeometry.h"
0021 #include "CondFormats/Alignment/interface/DetectorGlobalPosition.h"
0022 #include "CondFormats/Alignment/interface/AlignmentErrorsExtended.h"
0023 #include "CondFormats/AlignmentRecord/interface/GlobalPositionRcd.h"
0024 #include "CondFormats/AlignmentRecord/interface/DTAlignmentRcd.h"
0025 #include "CondFormats/AlignmentRecord/interface/DTAlignmentErrorRcd.h"
0026 #include "CondFormats/AlignmentRecord/interface/DTAlignmentErrorExtendedRcd.h"
0027 #include "Geometry/Records/interface/DTRecoGeometryRcd.h"
0028 #include "DataFormats/GeometrySurface/interface/Plane.h"
0029 #include "DataFormats/GeometrySurface/interface/Bounds.h"
0030 #include "DataFormats/GeometrySurface/interface/RectangularPlaneBounds.h"
0031 #include "Geometry/MuonNumbering/interface/MuonGeometryNumbering.h"
0032 #include "Geometry/Records/interface/MuonNumberingRecord.h"
0033 #include "Geometry/Records/interface/MuonGeometryRecord.h"
0034 #include "Geometry/Records/interface/GeometryFileRcd.h"
0035 #include "Geometry/DTGeometry/interface/DTGeometry.h"
0036 #include "Geometry/MuonNumbering/interface/MuonBaseNumber.h"
0037 #include "Geometry/MuonNumbering/interface/DTNumberingScheme.h"
0038 #include "DTGeometryBuilderFromDD4hep.h"
0039 #include "DD4hep/Detector.h"
0040 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0041
0042 #include <memory>
0043 #include <string>
0044 #include <string_view>
0045
0046 using namespace edm;
0047 using namespace std;
0048
0049 void DTGeometryBuilderFromDD4hep::buildGeometry(cms::DDFilteredView& fview,
0050 DTGeometry& geom,
0051 const MuonGeometryConstants& num) const {
0052 edm::LogVerbatim("DTGeometryBuilder") << "(0) DTGeometryBuilder - DD4hep ";
0053
0054 bool doChamber = fview.firstChild();
0055
0056 while (doChamber) {
0057 DTChamber* chamber = buildChamber(fview, num);
0058
0059 bool doSL = fview.nextSibling();
0060 while (doSL) {
0061 DTSuperLayer* sl = buildSuperLayer(fview, chamber, num);
0062
0063 fview.down();
0064 bool doLayers = fview.sibling();
0065 while (doLayers) {
0066 DTLayer* l = buildLayer(fview, sl, num);
0067 geom.add(l);
0068
0069 doLayers = fview.sibling();
0070 }
0071
0072 geom.add(sl);
0073 doSL = fview.nextSibling();
0074 }
0075 geom.add(chamber);
0076
0077 fview.parent();
0078 doChamber = fview.firstChild();
0079 }
0080 }
0081
0082 DTGeometryBuilderFromDD4hep::RCPPlane DTGeometryBuilderFromDD4hep::plane(const cms::DDFilteredView& fview,
0083 Bounds* bounds) const {
0084 const Double_t* tr = fview.trans();
0085 const Double_t* rot = fview.rot();
0086
0087 return RCPPlane(
0088 new Plane(Surface::PositionType(tr[0] / dd4hep::cm, tr[1] / dd4hep::cm, tr[2] / dd4hep::cm),
0089 Surface::RotationType(rot[0], rot[3], rot[6], rot[1], rot[4], rot[7], rot[2], rot[5], rot[8]),
0090 bounds));
0091 }
0092
0093 DTChamber* DTGeometryBuilderFromDD4hep::buildChamber(cms::DDFilteredView& fview,
0094 const MuonGeometryConstants& muonConstants) const {
0095 MuonGeometryNumbering mdddnum(muonConstants);
0096 DTNumberingScheme dtnum(muonConstants);
0097 int rawid = dtnum.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fview.history()));
0098
0099 DTChamberId detId(rawid);
0100 auto const& par = fview.parameters();
0101
0102 RCPPlane surf(
0103 plane(fview, new RectangularPlaneBounds(par[0] / dd4hep::cm, par[1] / dd4hep::cm, par[2] / dd4hep::cm)));
0104
0105 edm::LogVerbatim("DTGeometryBuilder") << "(1) detId: " << rawid << " par[0]: " << par[0] / dd4hep::cm
0106 << " par[1]: " << par[1] / dd4hep::cm << " par[2]: " << par[2] / dd4hep::cm;
0107
0108 DTChamber* chamber = new DTChamber(detId, surf);
0109
0110 return chamber;
0111 }
0112
0113 DTSuperLayer* DTGeometryBuilderFromDD4hep::buildSuperLayer(cms::DDFilteredView& fview,
0114 DTChamber* chamber,
0115 const MuonGeometryConstants& muonConstants) const {
0116 MuonGeometryNumbering mdddnum(muonConstants);
0117 DTNumberingScheme dtnum(muonConstants);
0118 int rawid = dtnum.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fview.history()));
0119
0120 DTSuperLayerId slId(rawid);
0121
0122 auto const& par = fview.parameters();
0123
0124 RCPPlane surf(
0125 plane(fview, new RectangularPlaneBounds(par[0] / dd4hep::cm, par[1] / dd4hep::cm, par[2] / dd4hep::cm)));
0126
0127 edm::LogVerbatim("DTGeometryBuilder") << "(2) detId: " << rawid << " par[0]: " << par[0] / dd4hep::cm
0128 << " par[1]: " << par[1] / dd4hep::cm << " par[2]: " << par[2] / dd4hep::cm;
0129
0130 DTSuperLayer* slayer = new DTSuperLayer(slId, surf, chamber);
0131
0132 chamber->add(slayer);
0133
0134 return slayer;
0135 }
0136
0137 DTLayer* DTGeometryBuilderFromDD4hep::buildLayer(cms::DDFilteredView& fview,
0138 DTSuperLayer* sl,
0139 const MuonGeometryConstants& muonConstants) const {
0140 MuonGeometryNumbering mdddnum(muonConstants);
0141 DTNumberingScheme dtnum(muonConstants);
0142 int rawid = dtnum.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fview.history()));
0143
0144 DTLayerId layId(rawid);
0145
0146 auto const& par = fview.parameters();
0147
0148 RCPPlane surf(
0149 plane(fview, new RectangularPlaneBounds(par[0] / dd4hep::cm, par[1] / dd4hep::cm, par[2] / dd4hep::cm)));
0150
0151 edm::LogVerbatim("DTGeometryBuilder") << "(3) detId: " << rawid << " par[0]: " << par[0] / dd4hep::cm
0152 << " par[1]: " << par[1] / dd4hep::cm << " par[2]: " << par[2] / dd4hep::cm;
0153
0154 fview.down();
0155 bool doWire = fview.sibling();
0156 int firstWire = fview.volume()->GetNumber();
0157 auto const& wpar = fview.parameters();
0158 float wireLength = wpar[1] / dd4hep::cm;
0159
0160 edm::LogVerbatim("DTGeometryBuilder") << "(4) detId: " << rawid << " wpar[1]: " << wpar[1] / dd4hep::cm
0161 << " firstWire: " << firstWire;
0162
0163 int WCounter = 0;
0164 while (doWire) {
0165 doWire = fview.checkChild();
0166 WCounter++;
0167 }
0168 fview.up();
0169
0170 DTTopology topology(firstWire, WCounter, wireLength);
0171
0172 DTLayerType layerType;
0173
0174 DTLayer* layer = new DTLayer(layId, surf, topology, layerType, sl);
0175
0176 sl->add(layer);
0177 return layer;
0178 }
0179
0180 void DTGeometryBuilderFromDD4hep::build(DTGeometry& geom,
0181 const cms::DDDetector* det,
0182 const MuonGeometryConstants& num,
0183 const dd4hep::SpecParRefs& refs) {
0184 cms::Volume top = det->worldVolume();
0185 cms::DDFilteredView fview(det, top);
0186 fview.mergedSpecifics(refs);
0187 buildGeometry(fview, geom, num);
0188 }