ME0GeometryBuilder Class Reference

#include <ME0GeometryBuilder.h>

Public Member Functions
ME0Geometry *	build (const DDCompactView *cview, const MuonGeometryConstants &muonConstants)
ME0Geometry *	build (const cms::DDCompactView *cview, const MuonGeometryConstants &muonConstants)
	ME0GeometryBuilder ()
	~ME0GeometryBuilder ()

Private Types
typedef ReferenceCountingPointer< BoundPlane >	ME0BoundPlane

Private Member Functions
ME0BoundPlane	boundPlane (const DDFilteredView &fv, Bounds *bounds, bool isOddChamber) const
ME0BoundPlane	boundPlane (const cms::DDFilteredView &fv, Bounds *bounds, bool isOddChamber) const
ME0Chamber *	buildChamber (DDFilteredView &fv, ME0DetId detId) const
ME0Chamber *	buildChamber (cms::DDFilteredView &fv, ME0DetId detId) const
ME0EtaPartition *	buildEtaPartition (DDFilteredView &fv, ME0DetId detId) const
ME0EtaPartition *	buildEtaPartition (cms::DDFilteredView &fv, ME0DetId detId) const
ME0Geometry *	buildGeometry (DDFilteredView &fview, const MuonGeometryConstants &muonConstants)
ME0Geometry *	buildGeometry (cms::DDFilteredView &fview, const MuonGeometryConstants &muonConstants)
ME0Layer *	buildLayer (DDFilteredView &fv, ME0DetId detId) const
ME0Layer *	buildLayer (cms::DDFilteredView &fv, ME0DetId detId) const

Private Attributes
std::map< ME0DetId, std::vector< ME0DetId > >	chids

Static Private Attributes
static constexpr double	k_ScaleFromDD4hep = (1.0 / dd4hep::cm)

Detailed Description

Definition at line 32 of file ME0GeometryBuilder.h.

Member Typedef Documentation

ME0BoundPlane

typedef ReferenceCountingPointer<BoundPlane> ME0GeometryBuilder::ME0BoundPlane

private

Definition at line 46 of file ME0GeometryBuilder.h.

Constructor & Destructor Documentation

ME0GeometryBuilder()

ME0GeometryBuilder::ME0GeometryBuilder

(

)

Definition at line 37 of file ME0GeometryBuilder.cc.

{}

~ME0GeometryBuilder()

ME0GeometryBuilder::~ME0GeometryBuilder

(

)

Definition at line 39 of file ME0GeometryBuilder.cc.

{}

Member Function Documentation

boundPlane() [1/2]

ME0GeometryBuilder::ME0BoundPlane ME0GeometryBuilder::boundPlane ( const DDFilteredView &  fv, Bounds *  bounds, bool  isOddChamber  ) const

private

Definition at line 335 of file ME0GeometryBuilder.cc.

References angle_units::operators::convertMmToCm(), nano_mu_digi_cff::float, TkRotation< T >::rotateAxes(), idealTransformation::rotation, DDFilteredView::rotation(), DDFilteredView::translation(), and x.

                                                                                          {
  // extract the position
  const DDTranslation& trans(fv.translation());
  const Surface::PositionType posResult(
      float(convertMmToCm(trans.x())), float(convertMmToCm(trans.y())), float(convertMmToCm(trans.z())));

  const DDRotationMatrix& rotation = fv.rotation();
  DD3Vector x, y, z;
  rotation.GetComponents(x, y, z);

  Surface::RotationType rotResult(float(x.X()),
                                  float(x.Y()),
                                  float(x.Z()),
                                  float(y.X()),
                                  float(y.Y()),
                                  float(y.Z()),
                                  float(z.X()),
                                  float(z.Y()),
                                  float(z.Z()));

  //Change of axes for the forward
  Basic3DVector<float> newX(1., 0., 0.);
  Basic3DVector<float> newY(0., 0., 1.);
  Basic3DVector<float> newZ(0., 1., 0.);
  newY *= -1;

  rotResult.rotateAxes(newX, newY, newZ);

  return ME0BoundPlane(new BoundPlane(posResult, rotResult, bounds));
}

boundPlane() [2/2]

ME0GeometryBuilder::ME0BoundPlane ME0GeometryBuilder::boundPlane ( const cms::DDFilteredView &  fv, Bounds *  bounds, bool  isOddChamber  ) const

private

Definition at line 490 of file ME0GeometryBuilder.cc.

References cms::DDFilteredView::rot(), TkRotation< T >::rotateAxes(), idealTransformation::rotation, cms::DDFilteredView::trans(), and x.

                                                                                          {
  // extract the position
  const Double_t* trans = fv.trans();
  Surface::PositionType posResult(
      k_ScaleFromDD4hep * trans[0], k_ScaleFromDD4hep * trans[1], k_ScaleFromDD4hep * trans[2]);

  // now the rotation
  DDRotationMatrix rotation;
  fv.rot(rotation);
  DD3Vector x, y, z;
  rotation.GetComponents(x, y, z);
  Surface::RotationType rotResult(float(x.X()),
                                  float(x.Y()),
                                  float(x.Z()),
                                  float(y.X()),
                                  float(y.Y()),
                                  float(y.Z()),
                                  float(z.X()),
                                  float(z.Y()),
                                  float(z.Z()));

  //Change of axes for the forward
  Basic3DVector<float> newX(1., 0., 0.);
  Basic3DVector<float> newY(0., 0., 1.);
  Basic3DVector<float> newZ(0., 1., 0.);
  newY *= -1;

  rotResult.rotateAxes(newX, newY, newZ);

  return ME0BoundPlane(new BoundPlane(posResult, rotResult, bounds));
}

build() [1/2]

ME0Geometry * ME0GeometryBuilder::build	(	const DDCompactView *	cview,
		const MuonGeometryConstants &	muonConstants
	)

Definition at line 41 of file ME0GeometryBuilder.cc.

References ALCARECOTkAlBeamHalo_cff::filter, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ME0GeometryESModule::produce().

                                                                                                             {
  std::string attribute = "MuStructure";
  std::string value = "MuonEndCapME0";
  DDSpecificsMatchesValueFilter filter{DDValue(attribute, value, 0.0)};
  DDFilteredView fview(*cview, filter);
  return this->buildGeometry(fview, muonConstants);
}

build() [2/2]

ME0Geometry * ME0GeometryBuilder::build	(	const cms::DDCompactView *	cview,
		const MuonGeometryConstants &	muonConstants
	)

Definition at line 50 of file ME0GeometryBuilder.cc.

References ALCARECOTkAlBeamHalo_cff::filter, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                                {
  std::string attribute = "MuStructure";
  std::string value = "MuonEndCapME0";
  const cms::DDFilter filter(attribute, value);
  cms::DDFilteredView fview(*cview, filter);
  return this->buildGeometry(fview, muonConstants);
}

buildChamber() [1/2]

ME0Chamber * ME0GeometryBuilder::buildChamber ( DDFilteredView &  fv, ME0DetId  detId  ) const

private

Definition at line 229 of file ME0GeometryBuilder.cc.

References B, b, relativeConstraints::chamber, angle_units::operators::convertMmToCm(), hcalRecHitTable_cff::detId, mps_fire::i, dttmaxenums::L, DDFilteredView::logicalPart(), LogTrace, DDName::name(), DDBase< N, C >::name(), DDSolid::parameters(), and DDLogicalPart::solid().

                                                                                     {
  LogTrace("ME0Geometry") << "buildChamber " << fv.logicalPart().name().name() << " " << detId << std::endl;
  DDBooleanSolid solid = (DDBooleanSolid)(fv.logicalPart().solid());

  std::vector<double> dpar = solid.parameters();

  double L = convertMmToCm(dpar[0]);  // length is along local Y
  double T = convertMmToCm(dpar[3]);  // thickness is long local Z
  double b = convertMmToCm(dpar[4]);  // bottom width is along local X
  double B = convertMmToCm(dpar[8]);  // top width is along local X

#ifdef EDM_ML_DEBUG
  LogTrace("ME0Geometry") << " name of logical part = " << fv.logicalPart().name().name() << std::endl;
  LogTrace("ME0Geometry") << " dpar is vector with size = " << dpar.size() << std::endl;
  for (unsigned int i = 0; i < dpar.size(); ++i) {
    LogTrace("ME0Geometry") << " dpar [" << i << "] = " << convertMmToCm(dpar[i]) << " cm " << std::endl;
  }
  LogTrace("ME0Geometry") << "size  b: " << b << "cm, B: " << B << "cm,  L: " << L << "cm, T: " << T << "cm "
                          << std::endl;
#endif

  bool isOdd = false;  // detId.chamber()%2;
  ME0BoundPlane surf(boundPlane(fv, new TrapezoidalPlaneBounds(b, B, L, T), isOdd));
  ME0Chamber* chamber = new ME0Chamber(detId.chamberId(), surf);
  return chamber;
}

buildChamber() [2/2]

ME0Chamber * ME0GeometryBuilder::buildChamber ( cms::DDFilteredView &  fv, ME0DetId  detId  ) const

private

Definition at line 434 of file ME0GeometryBuilder.cc.

References B, b, relativeConstraints::chamber, hcalRecHitTable_cff::detId, dttmaxenums::L, and cms::DDFilteredView::parameters().

                                                                                        {
  std::vector<double> dpar = fv.parameters();

  double L = k_ScaleFromDD4hep * dpar[3];
  double T = k_ScaleFromDD4hep * dpar[2];
  double b = k_ScaleFromDD4hep * dpar[0];
  double B = k_ScaleFromDD4hep * dpar[1];
  bool isOdd = false;
  ME0BoundPlane surf(boundPlane(fv, new TrapezoidalPlaneBounds(b, B, L, T), isOdd));
  ME0Chamber* chamber = new ME0Chamber(detId.chamberId(), surf);

  return chamber;
}

buildEtaPartition() [1/2]

ME0EtaPartition * ME0GeometryBuilder::buildEtaPartition ( DDFilteredView &  fv, ME0DetId  detId  ) const

private

Definition at line 283 of file ME0GeometryBuilder.cc.

References B, b, angle_units::operators::convertMmToCm(), DDfetch(), hcalRecHitTable_cff::detId, DDValue::doubles(), mps_fire::i, dttmaxenums::L, DDFilteredView::logicalPart(), LogTrace, GeomDetEnumerators::ME0, Skims_PA_cff::name, DDName::name(), DDBase< N, C >::name(), me0TriggerPseudoDigis_cff::nStrips, DDSolid::parameters(), DDLogicalPart::solid(), DDFilteredView::specifics(), HistogramManager_cfi::specs, AlCaHLTBitMon_QueryRunRegistry::string, submitPVValidationJobs::t, and to_string().

                                                                                               {
  LogTrace("ME0Geometry") << "buildEtaPartition " << fv.logicalPart().name().name() << " " << detId << std::endl;

  // EtaPartition specific parameter (nstrips and npads)
  DDValue numbOfStrips("nStrips");
  DDValue numbOfPads("nPads");
  std::vector<const DDsvalues_type*> specs(fv.specifics());
  double nStrips = 0., nPads = 0.;
  for (const auto& is : specs) {
    if (DDfetch(is, numbOfStrips))
      nStrips = numbOfStrips.doubles()[0];
    if (DDfetch(is, numbOfPads))
      nPads = numbOfPads.doubles()[0];
  }

  LogTrace("ME0Geometry") << ((nStrips == 0.) ? ("No nStrips found!!")
                                              : ("Number of strips: " + std::to_string(nStrips)));
  LogTrace("ME0Geometry") << ((nPads == 0.) ? ("No nPads found!!") : ("Number of pads: " + std::to_string(nPads)));

  // EtaPartition specific parameter (size)
  std::vector<double> dpar = fv.logicalPart().solid().parameters();
  double b = convertMmToCm(dpar[4]);  // half bottom edge
  double B = convertMmToCm(dpar[8]);  // half top edge
  double L = convertMmToCm(dpar[0]);  // half apothem
  double t = convertMmToCm(dpar[3]);  // half thickness

#ifdef EDM_ML_DEBUG
  LogTrace("ME0Geometry") << " name of logical part = " << fv.logicalPart().name().name() << std::endl;
  LogTrace("ME0Geometry") << " dpar is vector with size = " << dpar.size() << std::endl;
  for (unsigned int i = 0; i < dpar.size(); ++i) {
    LogTrace("ME0Geometry") << " dpar [" << i << "] = " << convertMmToCm(dpar[i]) << " cm " << std::endl;
  }
  LogTrace("ME0Geometry") << "size  b: " << b << "cm, B: " << B << "cm,  L: " << L << "cm, t: " << t << "cm "
                          << std::endl;
#endif

  std::vector<float> pars;
  pars.emplace_back(b);
  pars.emplace_back(B);
  pars.emplace_back(L);
  pars.emplace_back(nStrips);
  pars.emplace_back(nPads);

  bool isOdd = false;  // detId.chamber()%2;
  ME0BoundPlane surf(boundPlane(fv, new TrapezoidalPlaneBounds(b, B, L, t), isOdd));
  std::string name = fv.logicalPart().name().name();
  ME0EtaPartitionSpecs* e_p_specs = new ME0EtaPartitionSpecs(GeomDetEnumerators::ME0, name, pars);

  ME0EtaPartition* etaPartition = new ME0EtaPartition(detId, surf, e_p_specs);
  return etaPartition;
}

buildEtaPartition() [2/2]

ME0EtaPartition * ME0GeometryBuilder::buildEtaPartition ( cms::DDFilteredView &  fv, ME0DetId  detId  ) const

private

Definition at line 462 of file ME0GeometryBuilder.cc.

References B, b, hcalRecHitTable_cff::detId, nano_mu_digi_cff::float, dttmaxenums::L, GeomDetEnumerators::ME0, Skims_PA_cff::name, cms::DDFilteredView::name(), me0TriggerPseudoDigis_cff::nStrips, cms::DDFilteredView::parameters(), AlCaHLTBitMon_QueryRunRegistry::string, and submitPVValidationJobs::t.

                                                                                                  {
  //      auto nStrips = fv.get<double>("nStrips"); //it doesn't work
  //    auto nPads = fv.get<double>("nPads"); //it doesn't work

  auto nStrips = 384;  // from GEMSpecs
  auto nPads = 192;    // from GEMSpecs

  std::vector<double> dpar = fv.parameters();

  double b = k_ScaleFromDD4hep * dpar[0];
  double B = k_ScaleFromDD4hep * dpar[1];
  double L = k_ScaleFromDD4hep * dpar[3];
  double t = k_ScaleFromDD4hep * dpar[2];

  const std::vector<float> pars{float(dpar[0]), float(dpar[1]), float(dpar[3]), float(nStrips), float(nPads)};

  bool isOdd = false;
  ME0BoundPlane surf(boundPlane(fv, new TrapezoidalPlaneBounds(b, B, L, t), isOdd));

  std::string_view name = fv.name();

  ME0EtaPartitionSpecs* e_p_specs = new ME0EtaPartitionSpecs(GeomDetEnumerators::ME0, std::string(name), pars);

  ME0EtaPartition* etaPartition = new ME0EtaPartition(detId, surf, e_p_specs);

  return etaPartition;
}

buildGeometry() [1/2]

ME0Geometry * ME0GeometryBuilder::buildGeometry ( DDFilteredView &  fview, const MuonGeometryConstants &  muonConstants  )

private

Definition at line 58 of file ME0GeometryBuilder.cc.

References ME0Layer::add(), ME0Chamber::add(), ME0NumberingScheme::baseNumberToUnitNumber(), angle_units::operators::convertMmToCm(), hcalRecHitTable_cff::detId, DDFilteredView::firstChild(), DDFilteredView::geoHistory(), MuonGeometryNumbering::geoHistoryToBaseNumber(), relativeConstraints::geometry, mps_fire::i, DDFilteredView::logicalPart(), LogTrace, DDName::name(), DDBase< N, C >::name(), DDFilteredView::nextSibling(), DDSolid::parameters(), DDFilteredView::parent(), DetId::rawId(), nano_mu_digi_cff::rawId, and DDLogicalPart::solid().

                                                                                                             {
  ME0Geometry* geometry = new ME0Geometry();
  MuonGeometryNumbering mdddnum(muonConstants);
  ME0NumberingScheme me0Num(muonConstants);

  LogTrace("ME0Geometry") << "Building the geometry service";
  LogTrace("ME0Geometry") << "About to run through the ME0 structure\n"
                          << "Top level logical part: " << fv.logicalPart().name().name();

// ==========================================
// ===  Test to understand the structure  ===
// ==========================================
#ifdef EDM_ML_DEBUG
  bool testChambers = fv.firstChild();
  LogTrace("ME0Geometry") << "doChamber = fv.firstChild() = " << testChambers;

  while (testChambers) {
    // to etapartitions
    LogTrace("ME0Geometry") << "to layer " << fv.firstChild();
    LogTrace("ME0Geometry") << "to etapt " << fv.firstChild();
    int rawId = me0Num.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
    ME0DetId detId = ME0DetId(rawId);
    ME0DetId detIdCh = detId.chamberId();
    // back to chambers
    LogTrace("ME0Geometry") << "back to layer " << fv.parent();
    LogTrace("ME0Geometry") << "back to chamb " << fv.parent();

    LogTrace("ME0Geometry") << "In DoChambers Loop :: ME0DetId " << detId << " = " << detId.rawId()
                            << " (which belongs to ME0Chamber " << detIdCh << " = " << detIdCh.rawId() << ")";
    LogTrace("ME0Geometry") << "Second level logical part: " << fv.logicalPart().name().name();
    DDBooleanSolid solid2 = (DDBooleanSolid)(fv.logicalPart().solid());
    std::vector<double> dpar2 = solid2.parameters();
    std::stringstream parameters2;
    for (unsigned int i = 0; i < dpar2.size(); ++i) {
      parameters2 << " dpar[" << i << "]=" << convertMmToCm(dpar2[i]) << "cm ";
    }
    LogTrace("ME0Geometry") << "Second level parameters: vector with size = " << dpar2.size() << " and elements "
                            << parameters2.str();

    bool doLayers = fv.firstChild();

    LogTrace("ME0Geometry") << "doLayer = fv.firstChild() = " << doLayers;
    while (doLayers) {
      // to etapartitions
      LogTrace("ME0Geometry") << "to etapt " << fv.firstChild();
      int rawId = me0Num.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
      ME0DetId detId = ME0DetId(rawId);
      ME0DetId detIdLa = detId.layerId();
      // back to layers
      LogTrace("ME0Geometry") << "back to layer " << fv.parent();
      LogTrace("ME0Geometry") << "In DoLayers Loop :: ME0DetId " << detId << " = " << detId.rawId()
                              << " (which belongs to ME0Layer " << detIdLa << " = " << detIdLa.rawId() << ")";
      LogTrace("ME0Geometry") << "Third level logical part: " << fv.logicalPart().name().name();
      DDBooleanSolid solid3 = (DDBooleanSolid)(fv.logicalPart().solid());
      std::vector<double> dpar3 = solid3.parameters();
      std::stringstream parameters3;
      for (unsigned int i = 0; i < dpar3.size(); ++i) {
        parameters3 << " dpar[" << i << "]=" << convertMmToCm(dpar3[i]) << "cm ";
      }
      LogTrace("ME0Geometry") << "Third level parameters: vector with size = " << dpar3.size() << " and elements "
                              << parameters3.str();
      bool doEtaParts = fv.firstChild();

      LogTrace("ME0Geometry") << "doEtaPart = fv.firstChild() = " << doEtaParts;
      while (doEtaParts) {
        LogTrace("ME0Geometry") << "In DoEtaParts Loop :: ME0DetId " << detId << " = " << detId.rawId();
        LogTrace("ME0Geometry") << "Fourth level logical part: " << fv.logicalPart().name().name();
        DDBooleanSolid solid4 = (DDBooleanSolid)(fv.logicalPart().solid());
        std::vector<double> dpar4 = solid4.parameters();
        std::stringstream parameters4;
        for (unsigned int i = 0; i < dpar4.size(); ++i) {
          parameters4 << " dpar[" << i << "]=" << convertMmToCm(dpar4[i]) << "cm ";
        }
        LogTrace("ME0Geometry") << "Fourth level parameters: vector with size = " << dpar4.size() << " and elements "
                                << parameters4.str();

        doEtaParts = fv.nextSibling();
        LogTrace("ME0Geometry") << "doEtaPart = fv.nextSibling() = " << doEtaParts;
      }
      fv.parent();
      LogTrace("ME0Geometry") << "went back to parent :: name = " << fv.logicalPart().name().name()
                              << " will now ask for nextSibling";
      doLayers = fv.nextSibling();
      LogTrace("ME0Geometry") << "doLayer = fv.nextSibling() = " << doLayers;
    }
    fv.parent();
    LogTrace("ME0Geometry") << "went back to parent :: name = " << fv.logicalPart().name().name()
                            << " will now ask for nextSibling";
    testChambers = fv.nextSibling();
    LogTrace("ME0Geometry") << "doChamber = fv.nextSibling() = " << testChambers;
  }
  fv.parent();
#endif

  // ==========================================
  // === Here the Real ME0 Geometry Builder ===
  // ==========================================
  bool doChambers = fv.firstChild();

  while (doChambers) {
    // to etapartitions and back again to pick up DetId
    fv.firstChild();
    fv.firstChild();

#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ME0Geometry") << "MuonGeometry 1 " << fv.geoHistory() << " Levels "
                                    << mdddnum.geoHistoryToBaseNumber(fv.geoHistory()).getLevels();
#endif
    int rawId = me0Num.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
    ME0DetId detId = ME0DetId(rawId);
    ME0DetId detIdCh = detId.chamberId();

    fv.parent();
    fv.parent();
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ME0Geometry") << "MuonGeometry 2 " << fv.geoHistory() << " Levels "
                                    << mdddnum.geoHistoryToBaseNumber(fv.geoHistory()).getLevels();
#endif
    // build chamber
    ME0Chamber* me0Chamber = buildChamber(fv, detIdCh);
    geometry->add(me0Chamber);

    // loop over layers of the chamber
    bool doLayers = fv.firstChild();

    while (doLayers) {
      // to etapartitions and back again to pick up DetId
      fv.firstChild();
      int rawId = me0Num.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
      ME0DetId detId = ME0DetId(rawId);
      ME0DetId detIdLa = detId.layerId();
      fv.parent();
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("ME0Geometry") << "MuonGeometry 3 " << fv.geoHistory() << " Levels "
                                      << mdddnum.geoHistoryToBaseNumber(fv.geoHistory()).getLevels();
#endif
      // build layer
      ME0Layer* me0Layer = buildLayer(fv, detIdLa);
      me0Chamber->add(me0Layer);
      geometry->add(me0Layer);

      // loop over etapartitions of the layer
      bool doEtaParts = fv.firstChild();

      while (doEtaParts) {
        // pick up DetId
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("ME0Geometry") << "MuonGeometry 4 " << fv.geoHistory() << " Levels "
                                        << mdddnum.geoHistoryToBaseNumber(fv.geoHistory()).getLevels();
#endif
        int rawId = me0Num.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fv.geoHistory()));
        ME0DetId detId = ME0DetId(rawId);

        // build etapartition
        ME0EtaPartition* etaPart = buildEtaPartition(fv, detId);
        me0Layer->add(etaPart);
        geometry->add(etaPart);

        doEtaParts = fv.nextSibling();
      }
      fv.parent();

      doLayers = fv.nextSibling();
    }
    fv.parent();

    doChambers = fv.nextSibling();
  }
  return geometry;
}

buildGeometry() [2/2]

ME0Geometry * ME0GeometryBuilder::buildGeometry ( cms::DDFilteredView &  fview, const MuonGeometryConstants &  muonConstants  )

private

Definition at line 370 of file ME0GeometryBuilder.cc.

References ME0NumberingScheme::baseNumberToUnitNumber(), relativeConstraints::chamber, ME0DetId::chamberIdMask_, chambers, ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), TauDecayModes::dec, hcalRecHitTable_cff::detId, cms::DDFilteredView::firstChild(), MuonGeometryNumbering::geoHistoryToBaseNumber(), relativeConstraints::geometry, MuonGeometryConstants::getValue(), cms::DDFilteredView::history(), dqmdumpme::k, ME0DetId::layerIdMask_, hgcalTBTopologyTester_cfi::layers, cms::DDFilteredView::level(), cms::DDFilteredView::name(), and EgammaValidation_cff::num.

                                                                                                                {
  ME0Geometry* geometry = new ME0Geometry();
  MuonGeometryNumbering mdddnum(muonConstants);
  ME0NumberingScheme me0Num(muonConstants);

  static constexpr uint32_t levelChamber = 7;
  static constexpr uint32_t levelLayer = 8;
  uint32_t theLevelPart = muonConstants.getValue("level");
  uint32_t theSectorLevel = muonConstants.getValue("m0_sector") / theLevelPart;
  std::vector<ME0Chamber*> chambers;
  std::vector<ME0Layer*> layers;

  while (fv.firstChild()) {
    const auto& history = fv.history();
    MuonBaseNumber num(mdddnum.geoHistoryToBaseNumber(history));
    ME0DetId detId(me0Num.baseNumberToUnitNumber(num));
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("ME0Geometry") << fv.name() << " with " << history.tags.size() << " Levels and ID " << detId
                                    << " Mask " << std::hex << ME0DetId::chamberIdMask_ << ":" << ME0DetId::layerIdMask_
                                    << std::dec << " and " << ME0DetId(((detId.rawId()) & ME0DetId::chamberIdMask_))
                                    << ":" << ME0DetId(((detId.rawId()) & ME0DetId::layerIdMask_)) << " Sector Level "
                                    << theSectorLevel << ":" << history.tags.size() << ":" << fv.level();
    for (unsigned int k = 0; k < history.tags.size(); ++k)
      edm::LogVerbatim("ME0Geometry") << "[" << k << "] Tag " << history.tags[k] << " Offset " << history.offsets[k]
                                      << " copy " << history.copyNos[k];
#endif

    if (fv.level() == levelChamber) {
      // build chamber
      ME0Chamber* me0Chamber = buildChamber(fv, detId);
      chambers.emplace_back(me0Chamber);
    } else if (fv.level() == levelLayer) {
      // build layer
      ME0Layer* me0Layer = buildLayer(fv, detId);
      layers.emplace_back(me0Layer);
    } else if (history.tags.size() > theSectorLevel) {
      // build first eta partition
      ME0EtaPartition* etaPart = buildEtaPartition(fv, detId);
      geometry->add(etaPart);
    }
  }

  auto const& partitions = geometry->etaPartitions();
  for (auto& layer : layers) {
    uint32_t id0 = ((layer->id().rawId()) & ME0DetId::layerIdMask_);
    for (auto& etaPart : partitions) {
      if (((etaPart->id().rawId()) & ME0DetId::layerIdMask_) == id0) {
        layer->add(etaPart);
      }
    }
    geometry->add(layer);
  }
  for (auto& chamber : chambers) {
    uint32_t id0 = ((chamber->id().rawId()) & ME0DetId::chamberIdMask_);
    for (auto& layer : layers) {
      if (((layer->id().rawId()) & ME0DetId::chamberIdMask_) == id0) {
        chamber->add(layer);
      }
    }
    geometry->add(chamber);
  }
  return geometry;
}

buildLayer() [1/2]

ME0Layer * ME0GeometryBuilder::buildLayer ( DDFilteredView &  fv, ME0DetId  detId  ) const

private

Definition at line 256 of file ME0GeometryBuilder.cc.

References B, b, angle_units::operators::convertMmToCm(), hcalRecHitTable_cff::detId, mps_fire::i, dttmaxenums::L, DDFilteredView::logicalPart(), LogTrace, DDName::name(), DDBase< N, C >::name(), DDSolid::parameters(), DDLogicalPart::solid(), and submitPVValidationJobs::t.

                                                                                 {
  LogTrace("ME0Geometry") << "buildLayer " << fv.logicalPart().name().name() << " " << detId << std::endl;

  DDBooleanSolid solid = (DDBooleanSolid)(fv.logicalPart().solid());

  std::vector<double> dpar = solid.parameters();
  double L = convertMmToCm(dpar[0]);  // length is along local Y
  double t = convertMmToCm(dpar[3]);  // thickness is long local Z
  double b = convertMmToCm(dpar[4]);  // bottom width is along local X
  double B = convertMmToCm(dpar[8]);  // top width is along local X

#ifdef EDM_ML_DEBUG
  LogTrace("ME0Geometry") << " name of logical part = " << fv.logicalPart().name().name() << std::endl;
  LogTrace("ME0Geometry") << " dpar is vector with size = " << dpar.size() << std::endl;
  for (unsigned int i = 0; i < dpar.size(); ++i) {
    LogTrace("ME0Geometry") << " dpar [" << i << "] = " << convertMmToCm(dpar[i]) << " cm " << std::endl;
  }
  LogTrace("ME0Geometry") << "size  b: " << b << "cm, B: " << B << "cm,  L: " << L << "cm, t: " << t << "cm "
                          << std::endl;
#endif

  bool isOdd = false;  // detId.chamber()%2;
  ME0BoundPlane surf(boundPlane(fv, new TrapezoidalPlaneBounds(b, B, L, t), isOdd));
  ME0Layer* layer = new ME0Layer(detId.layerId(), surf);
  return layer;
}

buildLayer() [2/2]

ME0Layer * ME0GeometryBuilder::buildLayer ( cms::DDFilteredView &  fv, ME0DetId  detId  ) const

private

Definition at line 448 of file ME0GeometryBuilder.cc.

References B, b, hcalRecHitTable_cff::detId, dttmaxenums::L, cms::DDFilteredView::parameters(), and submitPVValidationJobs::t.

                                                                                    {
  std::vector<double> dpar = fv.parameters();

  double L = k_ScaleFromDD4hep * dpar[3];
  double t = k_ScaleFromDD4hep * dpar[2];
  double b = k_ScaleFromDD4hep * dpar[0];
  double B = k_ScaleFromDD4hep * dpar[1];
  bool isOdd = false;
  ME0BoundPlane surf(boundPlane(fv, new TrapezoidalPlaneBounds(b, B, L, t), isOdd));
  ME0Layer* layer = new ME0Layer(detId.layerId(), surf);

  return layer;
}

Member Data Documentation

chids

std::map<ME0DetId, std::vector<ME0DetId> > ME0GeometryBuilder::chids

private

Definition at line 44 of file ME0GeometryBuilder.h.

k_ScaleFromDD4hep

constexpr double ME0GeometryBuilder::k_ScaleFromDD4hep = (1.0 / dd4hep::cm)

staticprivate

Definition at line 66 of file ME0GeometryBuilder.h.