RPCGeometryParsFromDD Class Reference

#include <RPCGeometryParsFromDD.h>

Public Member Functions
void	build (const DDCompactView *cview, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rgeo)
void	build (const cms::DDCompactView *cview, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rgeo)
	RPCGeometryParsFromDD ()
	~RPCGeometryParsFromDD ()

Private Member Functions
void	buildGeometry (DDFilteredView &fview, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rgeo)
void	buildGeometry (cms::DDFilteredView &fview, const MuonGeometryConstants &muonConstants, RecoIdealGeometry &rgeo)

Detailed Description

Definition at line 27 of file RPCGeometryParsFromDD.h.

Constructor & Destructor Documentation

RPCGeometryParsFromDD()

RPCGeometryParsFromDD::RPCGeometryParsFromDD

(

)

Definition at line 27 of file RPCGeometryParsFromDD.cc.

{}

~RPCGeometryParsFromDD()

RPCGeometryParsFromDD::~RPCGeometryParsFromDD

(

)

Definition at line 29 of file RPCGeometryParsFromDD.cc.

{}

Member Function Documentation

build() [1/2]

void RPCGeometryParsFromDD::build	(	const DDCompactView *	cview,
		const MuonGeometryConstants &	muonConstants,
		RecoIdealGeometry &	rgeo
	)

Definition at line 32 of file RPCGeometryParsFromDD.cc.

References buildGeometry(), ALCARECOTkAlBeamHalo_cff::filter, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by RPCRecoIdealDBLoader::beginRun().

                                                           {
  const std::string attribute = "ReadOutName";
  const std::string value = "MuonRPCHits";

  // Asking only for the MuonRPC's
  DDSpecificsMatchesValueFilter filter{DDValue(attribute, value, 0.0)};
  DDFilteredView fview(*cview, filter);

  this->buildGeometry(fview, muonConstants, rgeo);
}

build() [2/2]

void RPCGeometryParsFromDD::build	(	const cms::DDCompactView *	cview,
		const MuonGeometryConstants &	muonConstants,
		RecoIdealGeometry &	rgeo
	)

Definition at line 47 of file RPCGeometryParsFromDD.cc.

References buildGeometry(), ALCARECOTkAlBeamHalo_cff::filter, and AlCaHLTBitMon_QueryRunRegistry::string.

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

buildGeometry() [1/2]

void RPCGeometryParsFromDD::buildGeometry ( DDFilteredView &  fview, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rgeo  )

private

Definition at line 59 of file RPCGeometryParsFromDD.cc.

References RPCNumberingScheme::baseNumberToUnitNumber(), gather_cfg::cout, DDfetch(), DDValue::doubles(), DDFilteredView::firstChild(), DDFilteredView::geoHistory(), MuonGeometryNumbering::geoHistoryToBaseNumber(), RecoIdealGeometry::insert(), createfilelist::int, DDFilteredView::logicalPart(), Skims_PA_cff::name, DDName::name(), DDBase< N, C >::name(), DDFilteredView::nextSibling(), me0TriggerPseudoDigis_cff::nStrips, DDSolid::parameters(), DetId::rawId(), DDFilteredView::rotation(), DDLogicalPart::solid(), DDFilteredView::specifics(), HistogramManager_cfi::specs, AlCaHLTBitMon_QueryRunRegistry::string, Calorimetry_cff::thickness, DDFilteredView::translation(), ApeEstimator_cff::width, x, y, and z.

Referenced by build().

                                                                   {
  for (bool doSubDets = fview.firstChild(); doSubDets == true; doSubDets = fview.nextSibling()) {
    // Get the Base Muon Number
    MuonGeometryNumbering mdddnum(muonConstants);
    MuonBaseNumber mbn = mdddnum.geoHistoryToBaseNumber(fview.geoHistory());

    // Get the The Rpc det Id
    RPCNumberingScheme rpcnum(muonConstants);
    const int detid = rpcnum.baseNumberToUnitNumber(mbn);

    RPCDetId rpcid(detid);

    DDValue numbOfStrips("nStrips");

    std::vector<const DDsvalues_type*> specs(fview.specifics());
    int nStrips = 0;
    for (auto& spec : specs) {
      if (DDfetch(spec, numbOfStrips)) {
        nStrips = int(numbOfStrips.doubles()[0]);
      }
    }
    if (nStrips == 0)
      std::cout << "No strip found!!" << std::endl;

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

    const std::string name = fview.logicalPart().name().name();

    edm::LogVerbatim("RPCGeometryParsFromDD")
        << " (1) "
        << "detid: " << detid << " name: " << name << " number of Strips: " << nStrips;

    const std::vector<std::string> strpars = {name};
    const DDTranslation& tran = fview.translation();

    const DDRotationMatrix& rota = fview.rotation();  //.Inverse();
    DD3Vector x, y, z;
    rota.GetComponents(x, y, z);
    std::vector<double> pars;
    if (dpar.size() == 3) {
      const double width = dpar[0];
      const double length = dpar[1];
      const double thickness = dpar[2];
      edm::LogVerbatim("RPCGeometryParsFromDD")
          << " (2) dpar.size() == 3, width: " << width << " length: " << length << " thickness: " << thickness;
      pars = {width, length, thickness, numbOfStrips.doubles()[0]};
    } else {
      const double dpar4 = dpar[4];
      const double dpar8 = dpar[8];
      const double dpar0 = dpar[0];
      pars = {dpar4, dpar8, dpar0, 0.4, numbOfStrips.doubles()[0]};
      edm::LogVerbatim("RPCGeometryParsFromDD")
          << " (3), else, dpar[4]: " << dpar4 << " dpar[8]: " << dpar8 << " dpar[0]: " << dpar0;
    }

    const std::vector<double> vtra = {tran.x(), tran.y(), tran.z()};

    edm::LogVerbatim("RPCGeometryParsFromDD")
        << " (4), tran.x() " << tran.x() << " tran.y(): " << tran.y() << " tran.z(): " << tran.z();

    const std::vector<double> vrot = {x.X(), x.Y(), x.Z(), y.X(), y.Y(), y.Z(), z.X(), z.Y(), z.Z()};
    edm::LogVerbatim("RPCGeometryParsFromDD")
        << " (5), x.X(), x.Y(), x.Z(), y.X(), y.Y(), y.Z(), z.X(), z.Y(), z.Z() " << x.X() << ", " << x.Y() << ", "
        << x.Z() << ", " << y.X() << ", " << y.Y() << ", " << y.Z() << ", " << z.X() << ", " << z.Y() << ", " << z.Z();
    rgeo.insert(rpcid.rawId(), vtra, vrot, pars, strpars);
  }
}

buildGeometry() [2/2]

void RPCGeometryParsFromDD::buildGeometry ( cms::DDFilteredView &  fview, const MuonGeometryConstants &  muonConstants, RecoIdealGeometry &  rgeo  )

private

Definition at line 131 of file RPCGeometryParsFromDD.cc.

References RPCNumberingScheme::baseNumberToUnitNumber(), cms::DDFilteredView::firstChild(), MuonGeometryNumbering::geoHistoryToBaseNumber(), cms::DDFilteredView::get(), cms::DDFilteredView::history(), RecoIdealGeometry::insert(), Skims_PA_cff::name, cms::DDFilteredView::name(), me0TriggerPseudoDigis_cff::nStrips, cms::DDFilteredView::parameters(), cms::DDFilteredView::rot(), makeMuonMisalignmentScenario::rot, cms::DDFilteredView::solid(), AlCaHLTBitMon_QueryRunRegistry::string, cms::DDFilteredView::translation(), x, y, and z.

                                                                   {
  while (fview.firstChild()) {
    MuonGeometryNumbering mdddnum(muonConstants);
    RPCNumberingScheme rpcnum(muonConstants);
    int rawidCh = rpcnum.baseNumberToUnitNumber(mdddnum.geoHistoryToBaseNumber(fview.history()));
    RPCDetId rpcid = RPCDetId(rawidCh);

    auto nStrips = fview.get<double>("nStrips");

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

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

    edm::LogVerbatim("RPCGeometryParsFromDD")
        << " (1), detid: " << rawidCh << " name: " << std::string(name) << " number of Strips: " << nStrips;

    const std::vector<std::string> strpars = {std::string(name)};

    std::vector<double> tran(3);
    tran[0] = static_cast<double>(fview.translation().X()) / dd4hep::mm;
    tran[1] = static_cast<double>(fview.translation().Y()) / dd4hep::mm;
    tran[2] = static_cast<double>(fview.translation().Z()) / dd4hep::mm;

    DDRotationMatrix rota;
    fview.rot(rota);
    DD3Vector x, y, z;
    rota.GetComponents(x, y, z);
    const std::vector<double> rot = {x.X(), x.Y(), x.Z(), y.X(), y.Y(), y.Z(), z.X(), z.Y(), z.Z()};

    if (dd4hep::isA<dd4hep::Box>(fview.solid())) {
      const std::vector<double> pars = {
          dpar[0] / dd4hep::mm, dpar[1] / dd4hep::mm, dpar[2] / dd4hep::mm, double(nStrips)};
      edm::LogVerbatim("RPCGeometryParsFromDD")
          << " (2), dd4hep::Box, width: " << dpar[0] / dd4hep::mm << " length: " << dpar[1] / dd4hep::mm
          << " thickness: " << dpar[2] / dd4hep::mm;
      rgeo.insert(rpcid, tran, rot, pars, strpars);
    } else {
      const double ti = 0.4;
      const std::vector<double> pars = {
          dpar[0] / dd4hep::mm, dpar[1] / dd4hep::mm, dpar[3] / dd4hep::mm, ti, double(nStrips)};
      edm::LogVerbatim("RPCGeometryParsFromDD") << " (3), else, dpar[0] (i.e. dpar[4] for DD): " << dpar[0] / dd4hep::mm
                                                << " dpar[1] (i.e. dpar[8] for DD): " << dpar[1] / dd4hep::mm
                                                << " dpar[3] (i.e. dpar[0] for DD): " << dpar[3] / dd4hep::mm;
      rgeo.insert(rpcid, tran, rot, pars, strpars);
    }
    edm::LogVerbatim("RPCGeometryParsFromDD")
        << " (4), tran.x(): " << tran[0] << " tran.y(): " << tran[1] << " tran.z(): " << tran[2];
    edm::LogVerbatim("RPCGeometryParsFromDD")
        << " (5), x.X(), x.Y(), x.Z(), y.X(), y.Y(), y.Z(), z.X(), z.Y(), z.Z(): " << x.X() << ", " << x.Y() << ", "
        << x.Z() << ", " << y.X() << ", " << y.Y() << ", " << y.Z() << ", " << z.X() << ", " << z.Y() << ", " << z.Z();
  }
}