FWRecoGeometryESProducer Class Reference

#include <FWRecoGeometryESProducer.h>

Inheritance diagram for FWRecoGeometryESProducer:

Public Member Functions
	FWRecoGeometryESProducer (const edm::ParameterSet &)
	FWRecoGeometryESProducer (const FWRecoGeometryESProducer &)=delete
const FWRecoGeometryESProducer &	operator= (const FWRecoGeometryESProducer &)=delete
std::unique_ptr< FWRecoGeometry >	produce (const FWRecoGeometryRecord &)
	~FWRecoGeometryESProducer (void) override
Public Member Functions inherited from edm::ESProducer
	ESProducer ()
	ESProducer (const ESProducer &)=delete
ESProxyIndex const *	getTokenIndices (unsigned int iIndex) const
ESRecordIndex const *	getTokenRecordIndices (unsigned int iIndex) const
bool	hasMayConsumes () const noexcept
size_t	numberOfTokenIndices (unsigned int iIndex) const
ESProducer const &	operator= (const ESProducer &)=delete
SerialTaskQueueChain &	queue ()
template<typename Record >
std::optional< std::vector < ESProxyIndex > >	updateFromMayConsumes (unsigned int iIndex, const Record &iRecord) const
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &) final
	~ESProducer () noexcept(false) override
Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()
	ESProxyFactoryProducer (const ESProxyFactoryProducer &)=delete
const ESProxyFactoryProducer &	operator= (const ESProxyFactoryProducer &)=delete
	~ESProxyFactoryProducer () noexcept(false) override
Public Member Functions inherited from edm::eventsetup::DataProxyProvider
void	createKeyedProxies (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)
	DataProxyProvider ()
	DataProxyProvider (const DataProxyProvider &)=delete
const ComponentDescription &	description () const
void	fillRecordsNotAllowingConcurrentIOVs (std::set< EventSetupRecordKey > &recordsNotAllowingConcurrentIOVs) const
virtual void	initConcurrentIOVs (EventSetupRecordKey const &key, unsigned int nConcurrentIOVs)
bool	isUsingRecord (const EventSetupRecordKey &key) const
KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey, unsigned int iovIndex=0)
const DataProxyProvider &	operator= (const DataProxyProvider &)=delete
void	setAppendToDataLabel (const edm::ParameterSet &)
void	setDescription (const ComponentDescription &iDescription)
std::set< EventSetupRecordKey >	usingRecords () const
virtual	~DataProxyProvider () noexcept(false)

Private Member Functions
void	ADD_PIXEL_TOPOLOGY (unsigned int rawid, const GeomDet *detUnit, FWRecoGeometry &)
void	addCaloGeometry (FWRecoGeometry &)
void	addCSCGeometry (FWRecoGeometry &)
void	addDTGeometry (FWRecoGeometry &)
void	addFTLGeometry (FWRecoGeometry &)
void	addGEMGeometry (FWRecoGeometry &)
void	addME0Geometry (FWRecoGeometry &)
void	addPixelBarrelGeometry (FWRecoGeometry &)
void	addPixelForwardGeometry (FWRecoGeometry &)
void	addRPCGeometry (FWRecoGeometry &)
void	addTECGeometry (FWRecoGeometry &)
void	addTIBGeometry (FWRecoGeometry &)
void	addTIDGeometry (FWRecoGeometry &)
void	addTOBGeometry (FWRecoGeometry &)
void	fillPoints (unsigned int id, std::vector< GlobalPoint >::const_iterator begin, std::vector< GlobalPoint >::const_iterator end, FWRecoGeometry &)
void	fillShapeAndPlacement (unsigned int id, const GeomDet *det, FWRecoGeometry &)
unsigned int	insert_id (unsigned int id, FWRecoGeometry &)
void	writeTrackerParametersXML (FWRecoGeometry &)

Private Attributes
bool	m_calo
const CaloGeometry *	m_caloGeom = nullptr
edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	m_caloGeomToken
unsigned int	m_current
const FastTimeGeometry *	m_ftlBarrelGeom = nullptr
edm::ESGetToken < FastTimeGeometry, IdealGeometryRecord >	m_ftlBarrelGeomToken
const FastTimeGeometry *	m_ftlEndcapGeom = nullptr
edm::ESGetToken < FastTimeGeometry, IdealGeometryRecord >	m_ftlEndcapGeomToken
bool	m_gem
bool	m_muon
bool	m_timing
bool	m_tracker
const TrackerGeometry *	m_trackerGeom = nullptr
const GlobalTrackingGeometry *	m_trackingGeom = nullptr
edm::ESGetToken < GlobalTrackingGeometry, GlobalTrackingGeometryRecord >	m_trackingGeomToken

Additional Inherited Members
Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)
Protected Types inherited from edm::ESProxyFactoryProducer
using	EventSetupRecordKey = eventsetup::EventSetupRecordKey
Protected Types inherited from edm::eventsetup::DataProxyProvider
using	KeyedProxiesVector = std::vector< std::pair< DataKey, std::shared_ptr< DataProxy >>>
Protected Member Functions inherited from edm::ESProducer
template<typename T >
auto	setWhatProduced (T *iThis, const es::Label &iLabel={})
template<typename T >
auto	setWhatProduced (T *iThis, const char *iLabel)
template<typename T >
auto	setWhatProduced (T *iThis, const std::string &iLabel)
template<typename T , typename TDecorator >
auto	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel={})
template<typename T , typename TReturn , typename TRecord >
auto	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const es::Label &iLabel={})
template<typename T , typename TReturn , typename TRecord , typename TArg >
ESConsumesCollectorT< TRecord >	setWhatProduced (T *iThis, TReturn(T::*iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel={})
void	usesResources (std::vector< std::string > const &)
Protected Member Functions inherited from edm::ESProxyFactoryProducer
template<class TFactory >
void	registerFactory (std::unique_ptr< TFactory > iFactory, const std::string &iLabel=std::string())
virtual void	registerFactoryWithKey (const EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ProxyFactoryBase > iFactory, const std::string &iLabel=std::string())
KeyedProxiesVector	registerProxies (const EventSetupRecordKey &, unsigned int iovIndex) override
Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
template<class T >
void	usingRecord ()
void	usingRecordWithKey (const EventSetupRecordKey &key)

Detailed Description

Definition at line 25 of file FWRecoGeometryESProducer.h.

Constructor & Destructor Documentation

FWRecoGeometryESProducer::FWRecoGeometryESProducer

(

const edm::ParameterSet &

pset

)

Definition at line 104 of file FWRecoGeometryESProducer.cc.

References edm::ParameterSet::getUntrackedParameter(), m_calo, m_caloGeomToken, m_ftlBarrelGeomToken, m_ftlEndcapGeomToken, m_gem, m_muon, m_timing, m_tracker, m_trackingGeomToken, or, and edm::ESProducer::setWhatProduced().

                                                                              : m_current(-1) {
  m_tracker = pset.getUntrackedParameter<bool>("Tracker", true);
  m_muon = pset.getUntrackedParameter<bool>("Muon", true);
  m_gem = pset.getUntrackedParameter<bool>("GEM", false);
  m_calo = pset.getUntrackedParameter<bool>("Calo", true);
  m_timing = pset.getUntrackedParameter<bool>("Timing", false);
  auto cc = setWhatProduced(this);

  if (m_muon)
    m_gem = true;
  if (m_tracker or m_muon or m_gem) {
    m_trackingGeomToken = cc.consumes();
  }
  if (m_timing) {
    m_ftlBarrelGeomToken = cc.consumes(edm::ESInputTag{"", "FastTimeBarrel"});
    m_ftlEndcapGeomToken = cc.consumes(edm::ESInputTag{"", "SFBX"});
  }
  if (m_calo) {
    m_caloGeomToken = cc.consumes();
  }
}

FWRecoGeometryESProducer::~FWRecoGeometryESProducer ( void  )

override

Definition at line 126 of file FWRecoGeometryESProducer.cc.

{}

FWRecoGeometryESProducer::FWRecoGeometryESProducer ( const FWRecoGeometryESProducer &  )

delete

Member Function Documentation

void FWRecoGeometryESProducer::ADD_PIXEL_TOPOLOGY ( unsigned int  rawid, const GeomDet *  detUnit, FWRecoGeometry &  fwRecoGeometry  )

private

Definition at line 42 of file FWRecoGeometryESProducer.cc.

References validate-o2o-wbm::f, FWRecoGeometry::idToName, PixelTopology::isItBigPixelInX(), PixelTopology::localX(), PixelTopology::localY(), PixelTopology::pitch(), and PixelGeomDetUnit::specificTopology().

Referenced by addPixelBarrelGeometry(), and addPixelForwardGeometry().

                                                                                  {
  const PixelGeomDetUnit* det = dynamic_cast<const PixelGeomDetUnit*>(detUnit);
  if (det) {
    const PixelTopology* topo = &det->specificTopology();

    std::pair<float, float> pitch = topo->pitch();
    fwRecoGeometry.idToName[rawid].topology[0] = pitch.first;
    fwRecoGeometry.idToName[rawid].topology[1] = pitch.second;

    fwRecoGeometry.idToName[rawid].topology[2] = topo->localX(0.f);  // offsetX
    fwRecoGeometry.idToName[rawid].topology[3] = topo->localY(0.f);  // offsetY

    // big pixels layout
    fwRecoGeometry.idToName[rawid].topology[4] = topo->isItBigPixelInX(80) ? 0 : 1;
  }
}

void FWRecoGeometryESProducer::addCaloGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 489 of file FWRecoGeometryESProducer.cc.

References EZArrayFL< T >::begin(), EZArrayFL< T >::end(), dataset::end, fillPoints(), DetId::Forward, ForwardEmpty, relativeConstraints::geom, CaloGeometry::getGeometry(), hgcal::RecHitTools::getLayerWithOffset(), HGCalGeometry::getNewCorners(), HGCalGeometry::getPosition(), hgcal::RecHitTools::getSiThickIndex(), CaloGeometry::getSubdetectorGeometry(), CaloGeometry::getValidDetIds(), DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, mps_fire::i, gpuClustering::id, FWRecoGeometry::idToName, insert_id(), hgcal::RecHitTools::isSilicon(), hgcal::RecHitTools::lastLayerEE(), m_caloGeom, hgcal::RecHitTools::setGeometry(), parallelization::uint(), and hgcal::RecHitTools::zside().

Referenced by produce().

                                                                             {
  std::vector<DetId> vid = m_caloGeom->getValidDetIds();  // Calo
  for (std::vector<DetId>::const_iterator it = vid.begin(), end = vid.end(); it != end; ++it) {
    unsigned int id = insert_id(it->rawId(), fwRecoGeometry);
    if (!((DetId::Forward == it->det()) || (DetId::HGCalEE == it->det()) || (DetId::HGCalHSi == it->det()) ||
          (DetId::HGCalHSc == it->det()))) {
      const CaloCellGeometry::CornersVec& cor = m_caloGeom->getGeometry(*it)->getCorners();
      fillPoints(id, cor.begin(), cor.end(), fwRecoGeometry);
    } else {
      DetId::Detector det = it->det();
      int subdet = (((DetId::HGCalEE == det) || (DetId::HGCalHSi == det) || (DetId::HGCalHSc == det)) ? ForwardEmpty
                                                                                                      : it->subdetId());
      const HGCalGeometry* geom = dynamic_cast<const HGCalGeometry*>(m_caloGeom->getSubdetectorGeometry(det, subdet));
      hgcal::RecHitTools rhtools;
      rhtools.setGeometry(*m_caloGeom);
      const auto cor = geom->getNewCorners(*it);

      // roll = yaw = pitch = 0
      fwRecoGeometry.idToName[id].matrix[0] = 1.0;
      fwRecoGeometry.idToName[id].matrix[4] = 1.0;
      fwRecoGeometry.idToName[id].matrix[8] = 1.0;

      // corners of the front face
      for (uint i = 0; i < (cor.size() - 1); ++i) {
        fwRecoGeometry.idToName[id].points[i * 3 + 0] = cor[i].x();
        fwRecoGeometry.idToName[id].points[i * 3 + 1] = cor[i].y();
        fwRecoGeometry.idToName[id].points[i * 3 + 2] = cor[i].z();
      }

      // center
      auto center = geom->getPosition(*it);
      fwRecoGeometry.idToName[id].points[(cor.size() - 1) * 3 + 0] = center.x();
      fwRecoGeometry.idToName[id].points[(cor.size() - 1) * 3 + 1] = center.y();
      fwRecoGeometry.idToName[id].points[(cor.size() - 1) * 3 + 2] = center.z();

      // thickness
      fwRecoGeometry.idToName[id].shape[3] = cor[cor.size() - 1].z();

      // total points
      fwRecoGeometry.idToName[id].topology[0] = cor.size() - 1;

      // Layer with Offset
      fwRecoGeometry.idToName[id].topology[1] = rhtools.getLayerWithOffset(it->rawId());

      // Zside, +/- 1
      fwRecoGeometry.idToName[id].topology[2] = rhtools.zside(it->rawId());

      // Is Silicon
      fwRecoGeometry.idToName[id].topology[3] = rhtools.isSilicon(it->rawId());

      // Silicon index
      fwRecoGeometry.idToName[id].topology[4] =
          rhtools.isSilicon(it->rawId()) ? rhtools.getSiThickIndex(it->rawId()) : -1.;

      // Last EE layer
      fwRecoGeometry.idToName[id].topology[5] = rhtools.lastLayerEE();
    }
  }
}

void FWRecoGeometryESProducer::addCSCGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 175 of file FWRecoGeometryESProducer.cc.

References CSCRadialStripTopology::angularWidth(), CSCRadialStripTopology::centreToIntersection(), CSCGeometry::chambers(), dataset::end, fillShapeAndPlacement(), GeomDet::geographicalId(), CSCLayer::geometry(), FWRecoGeometry::idToName, insert_id(), phase1PixelTopology::layer, CSCChamber::layers(), m_trackingGeom, DetId::Muon, CSCRadialStripTopology::phiOfOneEdge(), GlobalTrackingGeometry::slaveGeometry(), OffsetRadialStripTopology::stripOffset(), CSCLayerGeometry::topology(), CSCWireTopology::wireAngle(), CSCWireTopology::wireSpacing(), CSCLayerGeometry::wireTopology(), CSCRadialStripTopology::yAxisOrientation(), and CSCRadialStripTopology::yCentreOfStripPlane().

Referenced by produce().

                                                                            {
  DetId detId(DetId::Muon, 2);
  const CSCGeometry* cscGeometry = static_cast<const CSCGeometry*>(m_trackingGeom->slaveGeometry(detId));
  for (auto it = cscGeometry->chambers().begin(), end = cscGeometry->chambers().end(); it != end; ++it) {
    const CSCChamber* chamber = *it;

    if (chamber) {
      unsigned int rawid = chamber->geographicalId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, chamber, fwRecoGeometry);
      //
      // CSC layers geometry
      //
      for (std::vector<const CSCLayer*>::const_iterator lit = chamber->layers().begin(), lend = chamber->layers().end();
           lit != lend;
           ++lit) {
        const CSCLayer* layer = *lit;

        if (layer) {
          unsigned int rawid = layer->geographicalId();
          unsigned int current = insert_id(rawid, fwRecoGeometry);
          fillShapeAndPlacement(current, layer, fwRecoGeometry);

          const CSCStripTopology* stripTopology = layer->geometry()->topology();
          fwRecoGeometry.idToName[current].topology[0] = stripTopology->yAxisOrientation();
          fwRecoGeometry.idToName[current].topology[1] = stripTopology->centreToIntersection();
          fwRecoGeometry.idToName[current].topology[2] = stripTopology->yCentreOfStripPlane();
          fwRecoGeometry.idToName[current].topology[3] = stripTopology->phiOfOneEdge();
          fwRecoGeometry.idToName[current].topology[4] = stripTopology->stripOffset();
          fwRecoGeometry.idToName[current].topology[5] = stripTopology->angularWidth();

          const CSCWireTopology* wireTopology = layer->geometry()->wireTopology();
          fwRecoGeometry.idToName[current].topology[6] = wireTopology->wireSpacing();
          fwRecoGeometry.idToName[current].topology[7] = wireTopology->wireAngle();
        }
      }
    }
  }
}

void FWRecoGeometryESProducer::addDTGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 215 of file FWRecoGeometryESProducer.cc.

References DTTopology::cellHeight(), DTTopology::cellLenght(), DTTopology::cellWidth(), DTGeometry::chambers(), DTTopology::channels(), dataset::end, fillShapeAndPlacement(), DTTopology::firstChannel(), GeomDet::geographicalId(), DTLayer::id(), FWRecoGeometry::idToName, insert_id(), DTTopology::lastChannel(), phase1PixelTopology::layer, DTGeometry::layers(), m_trackingGeom, DetId::Muon, DetId::rawId(), GlobalTrackingGeometry::slaveGeometry(), DTLayer::specificTopology(), and GeomDet::surface().

Referenced by produce().

                                                                           {
  DetId detId(DetId::Muon, 1);
  const DTGeometry* dtGeometry = static_cast<const DTGeometry*>(m_trackingGeom->slaveGeometry(detId));

  //
  // DT chambers geometry
  //
  for (auto it = dtGeometry->chambers().begin(), end = dtGeometry->chambers().end(); it != end; ++it) {
    const DTChamber* chamber = *it;

    if (chamber) {
      unsigned int rawid = chamber->geographicalId().rawId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, chamber, fwRecoGeometry);
    }
  }

  // Fill in DT layer parameters
  for (auto it = dtGeometry->layers().begin(), end = dtGeometry->layers().end(); it != end; ++it) {
    const DTLayer* layer = *it;

    if (layer) {
      unsigned int rawid = layer->id().rawId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, layer, fwRecoGeometry);

      const DTTopology& topo = layer->specificTopology();
      const BoundPlane& surf = layer->surface();
      // Topology W/H/L:
      fwRecoGeometry.idToName[current].topology[0] = topo.cellWidth();
      fwRecoGeometry.idToName[current].topology[1] = topo.cellHeight();
      fwRecoGeometry.idToName[current].topology[2] = topo.cellLenght();
      fwRecoGeometry.idToName[current].topology[3] = topo.firstChannel();
      fwRecoGeometry.idToName[current].topology[4] = topo.lastChannel();
      fwRecoGeometry.idToName[current].topology[5] = topo.channels();

      // Bounds W/H/L:
      fwRecoGeometry.idToName[current].topology[6] = surf.bounds().width();
      fwRecoGeometry.idToName[current].topology[7] = surf.bounds().thickness();
      fwRecoGeometry.idToName[current].topology[8] = surf.bounds().length();
    }
  }
}

void FWRecoGeometryESProducer::addFTLGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 549 of file FWRecoGeometryESProducer.cc.

References fillPoints(), FastTimeGeometry::getCorners(), FastTimeGeometry::getValidDetIds(), insert_id(), m_ftlBarrelGeom, and m_ftlEndcapGeom.

Referenced by produce().

                                                                            {
  // do the barrel
  for (const auto& detid : m_ftlBarrelGeom->getValidDetIds()) {
    unsigned int id = insert_id(detid.rawId(), fwRecoGeometry);
    const auto& cor = m_ftlBarrelGeom->getCorners(detid);
    fillPoints(id, cor.begin(), cor.end(), fwRecoGeometry);
  }
  // do the endcap
  for (const auto& detid : m_ftlEndcapGeom->getValidDetIds()) {
    unsigned int id = insert_id(detid.rawId(), fwRecoGeometry);
    const auto& cor = m_ftlEndcapGeom->getCorners(detid);
    fillPoints(id, cor.begin(), cor.end(), fwRecoGeometry);
  }
}

void FWRecoGeometryESProducer::addGEMGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 288 of file FWRecoGeometryESProducer.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, GEMGeometry::chambers(), alignCSCRings::e, GEMGeometry::etaPartitions(), cppFunctionSkipper::exception, FWRecoGeometry::extraDet, fillShapeAndPlacement(), gpuClustering::id, FWRecoGeometry::idToName, insert_id(), m_trackingGeom, DetId::Muon, StripTopology::nstrips(), StripTopology::pitch(), GlobalTrackingGeometry::slaveGeometry(), StripTopology::stripLength(), GEMGeometry::superChambers(), and cms::Exception::what().

Referenced by produce().

                                                                            {
  //
  // GEM geometry
  //

  try {
    DetId detId(DetId::Muon, 4);
    const GEMGeometry* gemGeom = static_cast<const GEMGeometry*>(m_trackingGeom->slaveGeometry(detId));

    // add in superChambers - gem Segments are based on superChambers
    for (auto sc : gemGeom->superChambers()) {
      if (sc) {
        unsigned int rawid = sc->geographicalId().rawId();
        unsigned int current = insert_id(rawid, fwRecoGeometry);
        fillShapeAndPlacement(current, sc, fwRecoGeometry);
      }
    }
    // add in chambers
    for (auto ch : gemGeom->chambers()) {
      if (ch) {
        unsigned int rawid = ch->geographicalId().rawId();
        unsigned int current = insert_id(rawid, fwRecoGeometry);
        fillShapeAndPlacement(current, ch, fwRecoGeometry);
      }
    }
    // add in etaPartitions - gem rechits are based on etaPartitions
    for (auto roll : gemGeom->etaPartitions()) {
      if (roll) {
        unsigned int rawid = roll->geographicalId().rawId();
        unsigned int current = insert_id(rawid, fwRecoGeometry);
        fillShapeAndPlacement(current, roll, fwRecoGeometry);

        const StripTopology& topo = roll->specificTopology();
        fwRecoGeometry.idToName[current].topology[0] = topo.nstrips();
        fwRecoGeometry.idToName[current].topology[1] = topo.stripLength();
        fwRecoGeometry.idToName[current].topology[2] = topo.pitch();

        float height = topo.stripLength() / 2;
        LocalPoint lTop(0., height, 0.);
        LocalPoint lBottom(0., -height, 0.);
        fwRecoGeometry.idToName[current].topology[3] = roll->localPitch(lTop);
        fwRecoGeometry.idToName[current].topology[4] = roll->localPitch(lBottom);
        fwRecoGeometry.idToName[current].topology[5] = roll->npads();
      }
    }

    fwRecoGeometry.extraDet.Add(new TNamed("GEM", "GEM muon detector"));
    try {
      GEMDetId id(1, 1, 2, 1, 1, 1);
      m_trackingGeom->slaveGeometry(id);
      fwRecoGeometry.extraDet.Add(new TNamed("GE2", "GEM endcap station 2"));
    } catch (std::runtime_error& e) {
      std::cerr << e.what() << std::endl;
    }

  } catch (cms::Exception& exception) {
    edm::LogError("FWRecoGeometry") << " GEM geometry not found " << exception.what() << std::endl;
  }
}

void FWRecoGeometryESProducer::addME0Geometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 348 of file FWRecoGeometryESProducer.cc.

References ME0Geometry::etaPartitions(), cppFunctionSkipper::exception, FWRecoGeometry::extraDet, fillShapeAndPlacement(), FWRecoGeometry::idToName, insert_id(), m_trackingGeom, DetId::Muon, StripTopology::nstrips(), StripTopology::pitch(), GlobalTrackingGeometry::slaveGeometry(), StripTopology::stripLength(), and cms::Exception::what().

Referenced by produce().

                                                                            {
  //
  // ME0 geometry
  //

  DetId detId(DetId::Muon, 5);
  try {
    const ME0Geometry* me0Geom = static_cast<const ME0Geometry*>(m_trackingGeom->slaveGeometry(detId));
    for (auto roll : me0Geom->etaPartitions()) {
      if (roll) {
        unsigned int rawid = roll->geographicalId().rawId();
        unsigned int current = insert_id(rawid, fwRecoGeometry);
        fillShapeAndPlacement(current, roll, fwRecoGeometry);

        const StripTopology& topo = roll->specificTopology();
        fwRecoGeometry.idToName[current].topology[0] = topo.nstrips();
        fwRecoGeometry.idToName[current].topology[1] = topo.stripLength();
        fwRecoGeometry.idToName[current].topology[2] = topo.pitch();

        float height = topo.stripLength() / 2;
        LocalPoint lTop(0., height, 0.);
        LocalPoint lBottom(0., -height, 0.);
        fwRecoGeometry.idToName[current].topology[3] = roll->localPitch(lTop);
        fwRecoGeometry.idToName[current].topology[4] = roll->localPitch(lBottom);
        fwRecoGeometry.idToName[current].topology[5] = roll->npads();
      }
    }
    fwRecoGeometry.extraDet.Add(new TNamed("ME0", "ME0 muon detector"));
  } catch (cms::Exception& exception) {
    edm::LogError("FWRecoGeometry") << " ME0 geometry not found " << exception.what() << std::endl;
  }
}

void FWRecoGeometryESProducer::addPixelBarrelGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 381 of file FWRecoGeometryESProducer.cc.

References ADD_PIXEL_TOPOLOGY(), TrackerGeometry::detsPXB(), dataset::end, fillShapeAndPlacement(), GeomDet::geographicalId(), TrackerGeometry::idToDetUnit(), insert_id(), m_trackerGeom, and DetId::rawId().

Referenced by produce().

                                                                                    {
  for (TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsPXB().begin(),
                                                     end = m_trackerGeom->detsPXB().end();
       it != end;
       ++it) {
    const GeomDet* det = *it;

    if (det) {
      DetId detid = det->geographicalId();
      unsigned int rawid = detid.rawId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, det, fwRecoGeometry);

      ADD_PIXEL_TOPOLOGY(current, m_trackerGeom->idToDetUnit(detid), fwRecoGeometry);
    }
  }
}

void FWRecoGeometryESProducer::addPixelForwardGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 399 of file FWRecoGeometryESProducer.cc.

References ADD_PIXEL_TOPOLOGY(), TrackerGeometry::detsPXF(), dataset::end, fillShapeAndPlacement(), GeomDet::geographicalId(), TrackerGeometry::idToDetUnit(), insert_id(), m_trackerGeom, and DetId::rawId().

Referenced by produce().

                                                                                     {
  for (TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsPXF().begin(),
                                                     end = m_trackerGeom->detsPXF().end();
       it != end;
       ++it) {
    const GeomDet* det = *it;

    if (det) {
      DetId detid = det->geographicalId();
      unsigned int rawid = detid.rawId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, det, fwRecoGeometry);

      ADD_PIXEL_TOPOLOGY(current, m_trackerGeom->idToDetUnit(detid), fwRecoGeometry);
    }
  }
}

void FWRecoGeometryESProducer::addRPCGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 259 of file FWRecoGeometryESProducer.cc.

References EcnaPython_AdcPeg12_S1_10_R170298_1_0_150_Dee0::cerr, alignCSCRings::e, dataset::end, FWRecoGeometry::extraDet, fillShapeAndPlacement(), GeomDet::geographicalId(), gpuClustering::id, FWRecoGeometry::idToName, insert_id(), m_trackingGeom, DetId::Muon, StripTopology::nstrips(), StripTopology::pitch(), DetId::rawId(), RPCGeometry::rolls(), GlobalTrackingGeometry::slaveGeometry(), RPCRoll::specificTopology(), and StripTopology::stripLength().

Referenced by produce().

                                                                            {
  //
  // RPC rolls geometry
  //
  DetId detId(DetId::Muon, 3);
  const RPCGeometry* rpcGeom = static_cast<const RPCGeometry*>(m_trackingGeom->slaveGeometry(detId));
  for (auto it = rpcGeom->rolls().begin(), end = rpcGeom->rolls().end(); it != end; ++it) {
    const RPCRoll* roll = (*it);
    if (roll) {
      unsigned int rawid = roll->geographicalId().rawId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, roll, fwRecoGeometry);

      const StripTopology& topo = roll->specificTopology();
      fwRecoGeometry.idToName[current].topology[0] = topo.nstrips();
      fwRecoGeometry.idToName[current].topology[1] = topo.stripLength();
      fwRecoGeometry.idToName[current].topology[2] = topo.pitch();
    }
  }

  try {
    RPCDetId id(1, 1, 4, 1, 1, 1, 1);
    m_trackingGeom->slaveGeometry(detId);
    fwRecoGeometry.extraDet.Add(new TNamed("RE4", "RPC endcap station 4"));
  } catch (std::runtime_error& e) {
    std::cerr << e.what() << std::endl;
  }
}

void FWRecoGeometryESProducer::addTECGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 471 of file FWRecoGeometryESProducer.cc.

References ADD_SISTRIP_TOPOLOGY, TrackerGeometry::detsTEC(), dataset::end, fillShapeAndPlacement(), GeomDet::geographicalId(), TrackerGeometry::idToDet(), insert_id(), m_trackerGeom, and DetId::rawId().

Referenced by produce().

                                                                            {
  for (TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTEC().begin(),
                                                     end = m_trackerGeom->detsTEC().end();
       it != end;
       ++it) {
    const GeomDet* det = *it;

    if (det) {
      DetId detid = det->geographicalId();
      unsigned int rawid = detid.rawId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, det, fwRecoGeometry);

      ADD_SISTRIP_TOPOLOGY(current, m_trackerGeom->idToDet(detid));
    }
  }
}

void FWRecoGeometryESProducer::addTIBGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 417 of file FWRecoGeometryESProducer.cc.

References ADD_SISTRIP_TOPOLOGY, TrackerGeometry::detsTIB(), dataset::end, fillShapeAndPlacement(), GeomDet::geographicalId(), TrackerGeometry::idToDet(), insert_id(), m_trackerGeom, and DetId::rawId().

Referenced by produce().

                                                                            {
  for (TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTIB().begin(),
                                                     end = m_trackerGeom->detsTIB().end();
       it != end;
       ++it) {
    const GeomDet* det = *it;

    if (det) {
      DetId detid = det->geographicalId();
      unsigned int rawid = detid.rawId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, det, fwRecoGeometry);

      ADD_SISTRIP_TOPOLOGY(current, m_trackerGeom->idToDet(detid));
    }
  }
}

void FWRecoGeometryESProducer::addTIDGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 453 of file FWRecoGeometryESProducer.cc.

References ADD_SISTRIP_TOPOLOGY, TrackerGeometry::detsTID(), dataset::end, fillShapeAndPlacement(), GeomDet::geographicalId(), TrackerGeometry::idToDet(), insert_id(), m_trackerGeom, and DetId::rawId().

Referenced by produce().

                                                                            {
  for (TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTID().begin(),
                                                     end = m_trackerGeom->detsTID().end();
       it != end;
       ++it) {
    const GeomDet* det = *it;

    if (det) {
      DetId detid = det->geographicalId();
      unsigned int rawid = detid.rawId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, det, fwRecoGeometry);

      ADD_SISTRIP_TOPOLOGY(current, m_trackerGeom->idToDet(detid));
    }
  }
}

void FWRecoGeometryESProducer::addTOBGeometry ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 435 of file FWRecoGeometryESProducer.cc.

References ADD_SISTRIP_TOPOLOGY, TrackerGeometry::detsTOB(), dataset::end, fillShapeAndPlacement(), GeomDet::geographicalId(), TrackerGeometry::idToDet(), insert_id(), m_trackerGeom, and DetId::rawId().

Referenced by produce().

                                                                            {
  for (TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTOB().begin(),
                                                     end = m_trackerGeom->detsTOB().end();
       it != end;
       ++it) {
    const GeomDet* det = *it;

    if (det) {
      DetId detid = det->geographicalId();
      unsigned int rawid = detid.rawId();
      unsigned int current = insert_id(rawid, fwRecoGeometry);
      fillShapeAndPlacement(current, det, fwRecoGeometry);

      ADD_SISTRIP_TOPOLOGY(current, m_trackerGeom->idToDet(detid));
    }
  }
}

void FWRecoGeometryESProducer::fillPoints ( unsigned int  id, std::vector< GlobalPoint >::const_iterator  begin, std::vector< GlobalPoint >::const_iterator  end, FWRecoGeometry &  fwRecoGeometry  )

private

Definition at line 572 of file FWRecoGeometryESProducer.cc.

References cms::cuda::assert(), dataset::end, mps_fire::i, gpuClustering::id, FWRecoGeometry::idToName, and FWTGeoRecoGeometry::maxPoints_.

Referenced by addCaloGeometry(), and addFTLGeometry().

                                                                          {
  unsigned int index(0);
  for (std::vector<GlobalPoint>::const_iterator i = begin; i != end; ++i) {
    assert(index < FWTGeoRecoGeometry::maxPoints_ - 1);
    fwRecoGeometry.idToName[id].points[index] = i->x();
    fwRecoGeometry.idToName[id].points[++index] = i->y();
    fwRecoGeometry.idToName[id].points[++index] = i->z();
    ++index;
  }
}

void FWRecoGeometryESProducer::fillShapeAndPlacement ( unsigned int  id, const GeomDet *  det, FWRecoGeometry &  fwRecoGeometry  )

private

Shape of GeomDet

Definition at line 587 of file FWRecoGeometryESProducer.cc.

References b, b2, gpuClustering::id, FWRecoGeometry::idToName, GloballyPositioned< T >::position(), GloballyPositioned< T >::rotation(), GeomDet::surface(), PV3DBase< T, PVType, FrameType >::x(), TkRotation< T >::xx(), TkRotation< T >::xy(), TkRotation< T >::xz(), PV3DBase< T, PVType, FrameType >::y(), TkRotation< T >::yx(), TkRotation< T >::yy(), TkRotation< T >::yz(), PV3DBase< T, PVType, FrameType >::z(), TkRotation< T >::zx(), TkRotation< T >::zy(), and TkRotation< T >::zz().

Referenced by addCSCGeometry(), addDTGeometry(), addGEMGeometry(), addME0Geometry(), addPixelBarrelGeometry(), addPixelForwardGeometry(), addRPCGeometry(), addTECGeometry(), addTIBGeometry(), addTIDGeometry(), and addTOBGeometry().

                                                                                     {
  // Trapezoidal
  const Bounds* b = &((det->surface()).bounds());
  if (const TrapezoidalPlaneBounds* b2 = dynamic_cast<const TrapezoidalPlaneBounds*>(b)) {
    std::array<const float, 4> const& par = b2->parameters();

    // These parameters are half-lengths, as in CMSIM/GEANT3
    fwRecoGeometry.idToName[id].shape[0] = 1;
    fwRecoGeometry.idToName[id].shape[1] = par[0];  // hBottomEdge
    fwRecoGeometry.idToName[id].shape[2] = par[1];  // hTopEdge
    fwRecoGeometry.idToName[id].shape[3] = par[2];  // thickness
    fwRecoGeometry.idToName[id].shape[4] = par[3];  // apothem
  }
  if (const RectangularPlaneBounds* b2 = dynamic_cast<const RectangularPlaneBounds*>(b)) {
    // Rectangular
    fwRecoGeometry.idToName[id].shape[0] = 2;
    fwRecoGeometry.idToName[id].shape[1] = b2->width() * 0.5;      // half width
    fwRecoGeometry.idToName[id].shape[2] = b2->length() * 0.5;     // half length
    fwRecoGeometry.idToName[id].shape[3] = b2->thickness() * 0.5;  // half thickness
  }

  // Position of the DetUnit's center
  GlobalPoint pos = det->surface().position();
  fwRecoGeometry.idToName[id].translation[0] = pos.x();
  fwRecoGeometry.idToName[id].translation[1] = pos.y();
  fwRecoGeometry.idToName[id].translation[2] = pos.z();

  // Add the coeff of the rotation matrix
  // with a projection on the basis vectors
  TkRotation<float> detRot = det->surface().rotation();
  fwRecoGeometry.idToName[id].matrix[0] = detRot.xx();
  fwRecoGeometry.idToName[id].matrix[1] = detRot.yx();
  fwRecoGeometry.idToName[id].matrix[2] = detRot.zx();
  fwRecoGeometry.idToName[id].matrix[3] = detRot.xy();
  fwRecoGeometry.idToName[id].matrix[4] = detRot.yy();
  fwRecoGeometry.idToName[id].matrix[5] = detRot.zy();
  fwRecoGeometry.idToName[id].matrix[6] = detRot.xz();
  fwRecoGeometry.idToName[id].matrix[7] = detRot.yz();
  fwRecoGeometry.idToName[id].matrix[8] = detRot.zz();
}

unsigned int FWRecoGeometryESProducer::insert_id ( unsigned int  id, FWRecoGeometry &  fwRecoGeometry  )

private

Definition at line 564 of file FWRecoGeometryESProducer.cc.

References FWRecoGeometry::idToName, and m_current.

Referenced by addCaloGeometry(), addCSCGeometry(), addDTGeometry(), addFTLGeometry(), addGEMGeometry(), addME0Geometry(), addPixelBarrelGeometry(), addPixelForwardGeometry(), addRPCGeometry(), addTECGeometry(), addTIBGeometry(), addTIDGeometry(), and addTOBGeometry().

                                                                                                   {
  ++m_current;
  fwRecoGeometry.idToName.push_back(FWRecoGeom::Info());
  fwRecoGeometry.idToName.back().id = rawid;

  return m_current;
}

const FWRecoGeometryESProducer& FWRecoGeometryESProducer::operator= ( const FWRecoGeometryESProducer &  )

delete

std::unique_ptr< FWRecoGeometry > FWRecoGeometryESProducer::produce

(

const FWRecoGeometryRecord &

record

)

Definition at line 128 of file FWRecoGeometryESProducer.cc.

References addCaloGeometry(), addCSCGeometry(), addDTGeometry(), addFTLGeometry(), addGEMGeometry(), addME0Geometry(), addPixelBarrelGeometry(), addPixelForwardGeometry(), addRPCGeometry(), addTECGeometry(), addTIBGeometry(), addTIDGeometry(), addTOBGeometry(), edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::get(), edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), m_calo, m_caloGeom, m_caloGeomToken, m_current, m_ftlBarrelGeom, m_ftlBarrelGeomToken, m_ftlEndcapGeom, m_ftlEndcapGeomToken, m_gem, m_muon, m_timing, m_tracker, m_trackerGeom, m_trackingGeom, m_trackingGeomToken, GlobalTrackingGeometry::slaveGeometry(), edm::swap(), DetId::Tracker, and writeTrackerParametersXML().

                                                                                                  {
  using namespace edm;

  auto fwRecoGeometry = std::make_unique<FWRecoGeometry>();

  if (m_tracker || m_muon || m_gem) {
    m_trackingGeom = &record.get(m_trackingGeomToken);
  }

  if (m_tracker) {
    DetId detId(DetId::Tracker, 0);
    m_trackerGeom = static_cast<const TrackerGeometry*>(m_trackingGeom->slaveGeometry(detId));
    addPixelBarrelGeometry(*fwRecoGeometry);
    addPixelForwardGeometry(*fwRecoGeometry);
    addTIBGeometry(*fwRecoGeometry);
    addTIDGeometry(*fwRecoGeometry);
    addTOBGeometry(*fwRecoGeometry);
    addTECGeometry(*fwRecoGeometry);
    writeTrackerParametersXML(*fwRecoGeometry);
  }
  if (m_muon) {
    addDTGeometry(*fwRecoGeometry);
    addCSCGeometry(*fwRecoGeometry);
    addRPCGeometry(*fwRecoGeometry);
    addME0Geometry(*fwRecoGeometry);
  }
  if (m_gem) {
    addGEMGeometry(*fwRecoGeometry);
  }
  if (m_calo) {
    m_caloGeom = &record.get(m_caloGeomToken);
    addCaloGeometry(*fwRecoGeometry);
  }

  if (m_timing) {
    m_ftlBarrelGeom = &record.getRecord<CaloGeometryRecord>().get(m_ftlBarrelGeomToken);
    m_ftlEndcapGeom = &record.getRecord<CaloGeometryRecord>().get(m_ftlEndcapGeomToken);
    addFTLGeometry(*fwRecoGeometry);
  }

  fwRecoGeometry->idToName.resize(m_current + 1);
  std::vector<FWRecoGeom::Info>(fwRecoGeometry->idToName).swap(fwRecoGeometry->idToName);
  std::sort(fwRecoGeometry->idToName.begin(), fwRecoGeometry->idToName.end());

  return fwRecoGeometry;
}

void FWRecoGeometryESProducer::writeTrackerParametersXML ( FWRecoGeometry &  fwRecoGeometry )

private

Definition at line 630 of file FWRecoGeometryESProducer.cc.

References edmScanValgrind::buffer, contentValuesFiles::fullPath, edm::FileInPath::fullPath(), TrackerGeometry::isThere(), m_trackerGeom, GeomDetEnumerators::P1PXB, GeomDetEnumerators::P1PXEC, GeomDetEnumerators::P2OTB, GeomDetEnumerators::P2OTEC, GeomDetEnumerators::P2PXB, GeomDetEnumerators::P2PXEC, fed_dqm_sourceclient-live_cfg::path, AlCaHLTBitMon_QueryRunRegistry::string, submitPVValidationJobs::t, and FWRecoGeometry::trackerTopologyXML.

Referenced by produce().

                                                                                       {
  std::string path = "Geometry/TrackerCommonData/data/";
  if (m_trackerGeom->isThere(GeomDetEnumerators::P1PXB) || m_trackerGeom->isThere(GeomDetEnumerators::P1PXEC)) {
    path += "PhaseI/";
  } else if (m_trackerGeom->isThere(GeomDetEnumerators::P2PXB) || m_trackerGeom->isThere(GeomDetEnumerators::P2PXEC) ||
             m_trackerGeom->isThere(GeomDetEnumerators::P2OTB) || m_trackerGeom->isThere(GeomDetEnumerators::P2OTEC)) {
    path += "PhaseII/";
  }
  path += "trackerParameters.xml";
  std::string fullPath = edm::FileInPath(path).fullPath();
  std::ifstream t(fullPath);
  std::stringstream buffer;
  buffer << t.rdbuf();
  fwRecoGeometry.trackerTopologyXML = buffer.str();
}

Member Data Documentation

bool FWRecoGeometryESProducer::m_calo

private

Definition at line 75 of file FWRecoGeometryESProducer.h.

Referenced by FWRecoGeometryESProducer(), and produce().

const CaloGeometry* FWRecoGeometryESProducer::m_caloGeom = nullptr

private

Definition at line 66 of file FWRecoGeometryESProducer.h.

Referenced by addCaloGeometry(), and produce().

edm::ESGetToken<CaloGeometry, CaloGeometryRecord> FWRecoGeometryESProducer::m_caloGeomToken

private

Definition at line 64 of file FWRecoGeometryESProducer.h.

Referenced by FWRecoGeometryESProducer(), and produce().

unsigned int FWRecoGeometryESProducer::m_current

private

Definition at line 71 of file FWRecoGeometryESProducer.h.

Referenced by insert_id(), and produce().

const FastTimeGeometry* FWRecoGeometryESProducer::m_ftlBarrelGeom = nullptr

private

Definition at line 67 of file FWRecoGeometryESProducer.h.

Referenced by addFTLGeometry(), and produce().

edm::ESGetToken<FastTimeGeometry, IdealGeometryRecord> FWRecoGeometryESProducer::m_ftlBarrelGeomToken

private

Definition at line 62 of file FWRecoGeometryESProducer.h.

Referenced by FWRecoGeometryESProducer(), and produce().

const FastTimeGeometry* FWRecoGeometryESProducer::m_ftlEndcapGeom = nullptr

private

Definition at line 68 of file FWRecoGeometryESProducer.h.

Referenced by addFTLGeometry(), and produce().

edm::ESGetToken<FastTimeGeometry, IdealGeometryRecord> FWRecoGeometryESProducer::m_ftlEndcapGeomToken

private

Definition at line 63 of file FWRecoGeometryESProducer.h.

Referenced by FWRecoGeometryESProducer(), and produce().

bool FWRecoGeometryESProducer::m_gem

private

Definition at line 74 of file FWRecoGeometryESProducer.h.

Referenced by FWRecoGeometryESProducer(), and produce().

bool FWRecoGeometryESProducer::m_muon

private

Definition at line 73 of file FWRecoGeometryESProducer.h.

Referenced by FWRecoGeometryESProducer(), and produce().

bool FWRecoGeometryESProducer::m_timing

private

Definition at line 76 of file FWRecoGeometryESProducer.h.

Referenced by FWRecoGeometryESProducer(), and produce().

bool FWRecoGeometryESProducer::m_tracker

private

Definition at line 72 of file FWRecoGeometryESProducer.h.

Referenced by FWRecoGeometryESProducer(), and produce().

const TrackerGeometry* FWRecoGeometryESProducer::m_trackerGeom = nullptr

private

Definition at line 69 of file FWRecoGeometryESProducer.h.

Referenced by addPixelBarrelGeometry(), addPixelForwardGeometry(), addTECGeometry(), addTIBGeometry(), addTIDGeometry(), addTOBGeometry(), produce(), and writeTrackerParametersXML().

const GlobalTrackingGeometry* FWRecoGeometryESProducer::m_trackingGeom = nullptr

private

Definition at line 65 of file FWRecoGeometryESProducer.h.

Referenced by addCSCGeometry(), addDTGeometry(), addGEMGeometry(), addME0Geometry(), addRPCGeometry(), and produce().

edm::ESGetToken<GlobalTrackingGeometry, GlobalTrackingGeometryRecord> FWRecoGeometryESProducer::m_trackingGeomToken

private

Definition at line 61 of file FWRecoGeometryESProducer.h.

Referenced by FWRecoGeometryESProducer(), and produce().