#include <FWTGeoRecoGeometryESProducer.h>

Inheritance diagram for FWTGeoRecoGeometryESProducer:

Public Member Functions
	FWTGeoRecoGeometryESProducer (const edm::ParameterSet &)

std::unique_ptr < FWTGeoRecoGeometry >	produce (const FWTGeoRecoGeometryRecord &)

	~FWTGeoRecoGeometryESProducer (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 Types
enum	ERecoDet { kDummy, kSiPixel, kSiStrip, kMuonDT, kMuonRPC, kMuonCSC, kMuonGEM, kMuonME0, kECal, kHCal, kCaloTower, kHGCE, kHGCH }

Private Member Functions
void	addCaloTowerGeometry ()

void	addCSCGeometry ()

void	addDTGeometry ()

void	addEcalCaloGeometry ()

void	addGEMGeometry ()

void	addHcalCaloGeometryBarrel ()

void	addHcalCaloGeometryEndcap ()

void	addHcalCaloGeometryForward ()

void	addHcalCaloGeometryOuter ()

void	addME0Geometry ()

void	addPixelBarrelGeometry ()

void	addPixelForwardGeometry ()

void	addRPCGeometry ()

void	addTECGeometry ()

void	addTIBGeometry ()

void	addTIDGeometry ()

void	addTOBGeometry ()

TGeoManager *	createManager (int level)

TGeoShape *	createShape (const GeomDet *det)

TGeoVolume *	createVolume (const std::string &name, const GeomDet *det, ERecoDet=kDummy)

	FWTGeoRecoGeometryESProducer (const FWTGeoRecoGeometryESProducer &)

TGeoVolume *	GetDaughter (TGeoVolume mother, const char prefix, ERecoDet cidx, int id)

TGeoVolume *	GetDaughter (TGeoVolume mother, const char prefix, ERecoDet cidx)

TGeoMedium *	GetMedium (ERecoDet)

TGeoVolume *	GetTopHolder (const char *prefix, ERecoDet cidx)

const FWTGeoRecoGeometryESProducer &	operator= (const FWTGeoRecoGeometryESProducer &)

Private Attributes
bool	m_calo

const CaloGeometry *	m_caloGeom = nullptr

edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	m_caloGeomToken

TGeoMedium *	m_dummyMedium = nullptr

bool	m_gem

bool	m_muon

std::map< std::string, TGeoShape * >	m_nameToShape

std::map< ERecoDet, TGeoMedium * >	m_recoMedium

std::map< TGeoShape , TGeoVolume >	m_shapeToVolume

bool	m_tracker

const TrackerGeometry *	m_trackerGeom = nullptr

const TrackerTopology *	m_trackerTopology = nullptr

edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	m_trackerTopologyToken

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 31 of file FWTGeoRecoGeometryESProducer.h.

Member Enumeration Documentation

enum FWTGeoRecoGeometryESProducer::ERecoDet

private

Enumerator
kDummy
kSiPixel
kSiStrip
kMuonDT
kMuonRPC
kMuonCSC
kMuonGEM
kMuonME0
kECal
kHCal
kCaloTower
kHGCE
kHGCH

Definition at line 32 of file FWTGeoRecoGeometryESProducer.h.

                 {
     kDummy,
     kSiPixel,
     kSiStrip,
     kMuonDT,
     kMuonRPC,
     kMuonCSC,
     kMuonGEM,
     kMuonME0,
     kECal,
     kHCal,
     kCaloTower,
     kHGCE,
     kHGCH
   };

Constructor & Destructor Documentation

FWTGeoRecoGeometryESProducer::FWTGeoRecoGeometryESProducer ( const edm::ParameterSet & pset )

Definition at line 63 of file FWTGeoRecoGeometryESProducer.cc.

References edm::ParameterSet::getUntrackedParameter(), m_calo, m_caloGeomToken, m_gem, m_muon, m_tracker, m_trackerTopologyToken, m_trackingGeomToken, and edm::ESProducer::setWhatProduced().

                                                                                       : m_dummyMedium(nullptr) {
   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);
 
   if (m_muon)
     m_gem = true;
 
   auto cc = setWhatProduced(this);
   if (m_tracker || m_muon || m_gem) {
     m_trackingGeomToken = cc.consumes();
   }
   if (m_tracker) {
     m_trackerTopologyToken = cc.consumes();
   }
   if (m_calo) {
     m_caloGeomToken = cc.consumes();
   }
 }

FWTGeoRecoGeometryESProducer::~FWTGeoRecoGeometryESProducer ( void )

override

Definition at line 84 of file FWTGeoRecoGeometryESProducer.cc.

84 {}

FWTGeoRecoGeometryESProducer::FWTGeoRecoGeometryESProducer ( const FWTGeoRecoGeometryESProducer & )

private

Member Function Documentation

void FWTGeoRecoGeometryESProducer::addCaloTowerGeometry ( )

private

Definition at line 1074 of file FWTGeoRecoGeometryESProducer.cc.

References c, DetId::Calo, PVValHelper::dz, dataset::end, CaloCellGeometry::etaPos(), CaloCellGeometry::getCorners(), GetDaughter(), CaloGeometry::getGeometry(), GetMedium(), GetTopHolder(), CaloGeometry::getValidDetIds(), CaloTowerDetId::ieta(), kCaloTower, IdealObliquePrism::localCorners(), m_caloGeom, CaloCellGeometry::param(), CaloCellGeometry::phiPos(), gpuVertexFinder::printf(), makeMuonMisalignmentScenario::rot, CaloTowerDetId::SubdetId, and CaloTowerDetId::zside().

Referenced by produce().

                                                         {
   CaloVolMap caloShapeMapP;
   CaloVolMap caloShapeMapN;
 
   TGeoVolume* tv = GetTopHolder("CaloTower", kCaloTower);
   TGeoVolume* assembly = GetDaughter(tv, "CaloTower", kCaloTower);
 
   std::vector<DetId> vid = m_caloGeom->getValidDetIds(DetId::Calo, CaloTowerDetId::SubdetId);
   for (std::vector<DetId>::const_iterator it = vid.begin(), end = vid.end(); it != end; ++it) {
     CaloTowerDetId detid = CaloTowerDetId(it->rawId());
     const CaloCellGeometry* cellb = (m_caloGeom->getGeometry(*it)).get();
     const IdealObliquePrism* cell = dynamic_cast<const IdealObliquePrism*>(cellb);
     if (!cell) {
       printf("EC not oblique \n");
       continue;
     }
     TGeoVolume* volume = nullptr;
     CaloVolMap& caloShapeMap = (cell->etaPos() > 0) ? caloShapeMapP : caloShapeMapN;
     CaloVolMap::iterator volIt = caloShapeMap.find(cell->param());
     if (volIt == caloShapeMap.end()) {
       IdealObliquePrism::Pt3DVec lc(8);
       IdealObliquePrism::Pt3D ref;
       IdealObliquePrism::localCorners(lc, cell->param(), ref);
       HepGeom::Vector3D<float> lCenter;
       for (int c = 0; c < 8; ++c)
         lCenter += lc[c];
       lCenter *= 0.125;
 
       static const int arrP[] = {3, 2, 1, 0, 7, 6, 5, 4};
       static const int arrN[] = {7, 6, 5, 4, 3, 2, 1, 0};
       const int* arr = (detid.ieta() > 0) ? &arrP[0] : &arrN[0];
 
       double points[16];
       for (int c = 0; c < 8; ++c) {
         points[c * 2 + 0] = lc[arr[c]].x() - lCenter.x();
         points[c * 2 + 1] = lc[arr[c]].y() - lCenter.y();
       }
 
       float dz = (lc[4].z() - lc[0].z()) * 0.5;
       TGeoShape* solid = new TGeoArb8(dz, &points[0]);
       volume = new TGeoVolume("ecal oblique prism", solid, GetMedium(kCaloTower));
       caloShapeMap[cell->param()] = volume;
     } else {
       volume = volIt->second;
     }
 
     HepGeom::Vector3D<float> gCenter;
     CaloCellGeometry::CornersVec const& gc = cell->getCorners();
     for (int c = 0; c < 8; ++c) {
       gCenter += HepGeom::Vector3D<float>(gc[c].x(), gc[c].y(), gc[c].z());
     }
     gCenter *= 0.125;
 
     TGeoTranslation gtr(gCenter.x(), gCenter.y(), gCenter.z());
     TGeoRotation rot;
     rot.SetAngles(cell->phiPos() * TMath::RadToDeg(), 0, 0);
 
     TGeoVolume* holder = GetDaughter(assembly, "side", kCaloTower, detid.zside());
     holder = GetDaughter(holder, "ieta", kCaloTower, detid.ieta());
     std::stringstream nname;
     nname << detid;
     AddLeafNode(holder, volume, nname.str().c_str(), new TGeoCombiTrans(gtr, rot));
   }
 }

void FWTGeoRecoGeometryESProducer::addCSCGeometry ( )

private

Definition at line 635 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), CSCDetId, TrackingGeometry::dets(), GetDaughter(), GetTopHolder(), kMuonCSC, kMuonRPC, phase1PixelTopology::layer, m_trackingGeom, alignCSCRings::s, GlobalTrackingGeometry::slaveGeometry(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

                                                   {
   if (!m_trackingGeom->slaveGeometry(CSCDetId()))
     throw cms::Exception("FatalError") << "Cannnot find CSCGeometry\n";
 
   TGeoVolume* tv = GetTopHolder("Muon", kMuonRPC);
   TGeoVolume* assembly = GetDaughter(tv, "CSC", kMuonCSC);
 
   auto const& cscGeom = m_trackingGeom->slaveGeometry(CSCDetId())->dets();
   for (auto it = cscGeom.begin(), itEnd = cscGeom.end(); it != itEnd; ++it) {
     unsigned int rawid = (*it)->geographicalId();
     CSCDetId detId(rawid);
     std::stringstream s;
     s << "CSC" << detId;
     std::string name = s.str();
 
     TGeoVolume* child = nullptr;
 
     if (auto chamber = dynamic_cast<const CSCChamber*>(*it))
       child = createVolume(name, chamber, kMuonCSC);
     else if (auto* layer = dynamic_cast<const CSCLayer*>(*it))
       child = createVolume(name, layer, kMuonCSC);
 
     if (child) {
       TGeoVolume* holder = GetDaughter(assembly, "Endcap", kMuonCSC, detId.endcap());
       holder = GetDaughter(holder, "Station", kMuonCSC, detId.station());
       holder = GetDaughter(holder, "Ring", kMuonCSC, detId.ring());
       holder = GetDaughter(holder, "Chamber", kMuonCSC, detId.chamber());
 
       //   holder->AddNode(child, 1,  createPlacement( *it ));
       AddLeafNode(holder, child, name.c_str(), createPlacement(*it));
     }
   }
 }

void FWTGeoRecoGeometryESProducer::addDTGeometry ( )

private

Definition at line 561 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), TrackingGeometry::dets(), DTChamberId, DTLayerId, DTSuperLayerId, dataset::end, GetDaughter(), GetTopHolder(), kMuonDT, kMuonRPC, phase1PixelTopology::layer, m_trackingGeom, alignCSCRings::s, GlobalTrackingGeometry::slaveGeometry(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

                                                  {
   TGeoVolume* tv = GetTopHolder("Muon", kMuonRPC);
   TGeoVolume* assemblyTop = GetDaughter(tv, "DT", kMuonDT);
 
   //
   // DT chambers geometry
   //
   {
     TGeoVolume* assembly = GetDaughter(assemblyTop, "DTChamber", kMuonDT);
     auto const& dtChamberGeom = m_trackingGeom->slaveGeometry(DTChamberId())->dets();
     for (auto it = dtChamberGeom.begin(), end = dtChamberGeom.end(); it != end; ++it) {
       if (auto chamber = dynamic_cast<const DTChamber*>(*it)) {
         DTChamberId detid = chamber->geographicalId();
         std::stringstream s;
         s << detid;
         std::string name = s.str();
 
         TGeoVolume* child = createVolume(name, chamber, kMuonDT);
         TGeoVolume* holder = GetDaughter(assembly, "Wheel", kMuonDT, detid.wheel());
         holder = GetDaughter(holder, "Station", kMuonDT, detid.station());
         holder = GetDaughter(holder, "Sector", kMuonDT, detid.sector());
 
         AddLeafNode(holder, child, name.c_str(), createPlacement(chamber));
       }
     }
   }
 
   // Fill in DT super layer parameters
   {
     TGeoVolume* assembly = GetDaughter(assemblyTop, "DTSuperLayer", kMuonDT);
     auto const& dtSuperLayerGeom = m_trackingGeom->slaveGeometry(DTSuperLayerId())->dets();
     for (auto it = dtSuperLayerGeom.begin(), end = dtSuperLayerGeom.end(); it != end; ++it) {
       if (auto* superlayer = dynamic_cast<const DTSuperLayer*>(*it)) {
         DTSuperLayerId detid(DetId(superlayer->geographicalId()));
         std::stringstream s;
         s << detid;
         std::string name = s.str();
 
         TGeoVolume* child = createVolume(name, superlayer, kMuonDT);
 
         TGeoVolume* holder = GetDaughter(assembly, "Wheel", kMuonDT, detid.wheel());
         holder = GetDaughter(holder, "Station", kMuonDT, detid.station());
         holder = GetDaughter(holder, "Sector", kMuonDT, detid.sector());
         holder = GetDaughter(holder, "SuperLayer", kMuonDT, detid.superlayer());
         AddLeafNode(holder, child, name.c_str(), createPlacement(superlayer));
       }
     }
   }
   // Fill in DT layer parameters
   {
     TGeoVolume* assembly = GetDaughter(assemblyTop, "DTLayer", kMuonDT);
     auto const& dtLayerGeom = m_trackingGeom->slaveGeometry(DTLayerId())->dets();
     for (auto it = dtLayerGeom.begin(), end = dtLayerGeom.end(); it != end; ++it) {
       if (auto layer = dynamic_cast<const DTLayer*>(*it)) {
         DTLayerId detid(DetId(layer->geographicalId()));
 
         std::stringstream s;
         s << detid;
         std::string name = s.str();
 
         TGeoVolume* child = createVolume(name, layer, kMuonDT);
 
         TGeoVolume* holder = GetDaughter(assembly, "Wheel", kMuonDT, detid.wheel());
         holder = GetDaughter(holder, "Station", kMuonDT, detid.station());
         holder = GetDaughter(holder, "Sector", kMuonDT, detid.sector());
         holder = GetDaughter(holder, "SuperLayer", kMuonDT, detid.superlayer());
         holder = GetDaughter(holder, "Layer", kMuonDT, detid.layer());
         AddLeafNode(holder, child, name.c_str(), createPlacement(layer));
       }
     }
   }
 }

void FWTGeoRecoGeometryESProducer::addEcalCaloGeometry ( void )

private

Definition at line 1201 of file FWTGeoRecoGeometryESProducer.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, dataset::end, CaloCellGeometry::getCorners(), GetDaughter(), getEcalTrans(), CaloGeometry::getGeometry(), GetMedium(), GetTopHolder(), CaloGeometry::getValidDetIds(), EBDetId::ieta(), EEDetId::ix(), kECal, m_caloGeom, makeEcalShape(), CaloCellGeometry::param(), gpuVertexFinder::printf(), EBDetId::zside(), and EEDetId::zside().

Referenced by produce().

                                                            {
   TGeoVolume* tv = GetTopHolder("ECal", kECal);
   CaloVolMap caloShapeMap;
 
   {
     TGeoVolume* assembly = GetDaughter(tv, "ECalBarrel", kECal);
 
     std::vector<DetId> vid = m_caloGeom->getValidDetIds(DetId::Ecal, EcalSubdetector::EcalBarrel);
     for (std::vector<DetId>::const_iterator it = vid.begin(), end = vid.end(); it != end; ++it) {
       EBDetId detid(*it);
       const CaloCellGeometry* cellb = (m_caloGeom->getGeometry(*it)).get();
       const TruncatedPyramid* cell = dynamic_cast<const TruncatedPyramid*>(cellb);
       if (!cell) {
         printf("ecalBarrel cell not a TruncatedPyramid !!\n");
         return;
       }
 
       TGeoVolume* volume = nullptr;
       CaloVolMap::iterator volIt = caloShapeMap.find(cell->param());
       if (volIt == caloShapeMap.end()) {
         volume = new TGeoVolume("EE TruncatedPyramid", makeEcalShape(cell), GetMedium(kECal));
         caloShapeMap[cell->param()] = volume;
       } else {
         volume = volIt->second;
       }
       TGeoHMatrix* mtx = getEcalTrans(cell->getCorners());
       TGeoVolume* holder = GetDaughter(assembly, "side", kECal, detid.zside());
       holder = GetDaughter(holder, "ieta", kECal, detid.ieta());
       std::stringstream nname;
       nname << detid;
       AddLeafNode(holder, volume, nname.str().c_str(), mtx);
     }
   }
 
   {
     TGeoVolume* assembly = GetDaughter(tv, "ECalEndcap", kECal);
 
     std::vector<DetId> vid = m_caloGeom->getValidDetIds(DetId::Ecal, EcalSubdetector::EcalEndcap);
     for (std::vector<DetId>::const_iterator it = vid.begin(), end = vid.end(); it != end; ++it) {
       EEDetId detid(*it);
       const CaloCellGeometry* cellb = (m_caloGeom->getGeometry(*it)).get();
       const TruncatedPyramid* cell = dynamic_cast<const TruncatedPyramid*>(cellb);
       if (!cell) {
         printf("ecalEndcap cell not a TruncatedPyramid !!\n");
         continue;
       }
 
       TGeoVolume* volume = nullptr;
       CaloVolMap::iterator volIt = caloShapeMap.find(cell->param());
       if (volIt == caloShapeMap.end()) {
         volume = new TGeoVolume("EE TruncatedPyramid", makeEcalShape(cell), GetMedium(kECal));
         caloShapeMap[cell->param()] = volume;
       } else {
         volume = volIt->second;
       }
       TGeoHMatrix* mtx = getEcalTrans(cell->getCorners());
       TGeoVolume* holder = GetDaughter(assembly, "side", kECal, detid.zside());
       holder = GetDaughter(holder, "ix", kECal, detid.ix());
       std::stringstream nname;
       nname << detid;
       AddLeafNode(holder, volume, nname.str().c_str(), mtx);
     }
   }
 }

void FWTGeoRecoGeometryESProducer::addGEMGeometry ( )

private

Definition at line 671 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), dataset::end, GEMGeometry::etaPartitions(), cppFunctionSkipper::exception, MuonSubdetId::GEM, GeomDet::geographicalId(), GetDaughter(), GetTopHolder(), kMuonGEM, kMuonRPC, m_trackingGeom, DetId::Muon, alignCSCRings::s, GlobalTrackingGeometry::slaveGeometry(), AlCaHLTBitMon_QueryRunRegistry::string, GEMGeometry::superChambers(), and cms::Exception::what().

Referenced by produce().

                                                   {
   try {
     DetId detId(DetId::Muon, MuonSubdetId::GEM);
     const GEMGeometry* gemGeom = static_cast<const GEMGeometry*>(m_trackingGeom->slaveGeometry(detId));
 
     TGeoVolume* tv = GetTopHolder("Muon", kMuonRPC);
     TGeoVolume* assemblyTop = GetDaughter(tv, "GEM", kMuonGEM);
 
     {
       TGeoVolume* assembly = GetDaughter(assemblyTop, "GEMSuperChambers", kMuonGEM);
       for (auto it = gemGeom->superChambers().begin(), end = gemGeom->superChambers().end(); it != end; ++it) {
         const GEMSuperChamber* sc = (*it);
         if (sc) {
           GEMDetId detid = sc->geographicalId();
           std::stringstream s;
           s << detid;
           std::string name = s.str();
 
           TGeoVolume* child = createVolume(name, sc, kMuonGEM);
 
           TGeoVolume* holder = GetDaughter(assembly, "SuperChamber Region", kMuonGEM, detid.region());
           holder = GetDaughter(holder, "Ring", kMuonGEM, detid.ring());
           holder = GetDaughter(holder, "Station", kMuonGEM, detid.station());
           holder = GetDaughter(holder, "Chamber", kMuonGEM, detid.chamber());
 
           AddLeafNode(holder, child, name.c_str(), createPlacement(*it));
         }
       }
     }
 
     {
       TGeoVolume* assembly = GetDaughter(assemblyTop, "GEMetaPartitions", kMuonGEM);
       for (auto it = gemGeom->etaPartitions().begin(), end = gemGeom->etaPartitions().end(); it != end; ++it) {
         const GEMEtaPartition* roll = (*it);
         if (roll) {
           GEMDetId detid = roll->geographicalId();
           std::stringstream s;
           s << detid;
           std::string name = s.str();
 
           TGeoVolume* child = createVolume(name, roll, kMuonGEM);
 
           TGeoVolume* holder = GetDaughter(assembly, "ROLL Region", kMuonGEM, detid.region());
           holder = GetDaughter(holder, "Ring", kMuonGEM, detid.ring());
           holder = GetDaughter(holder, "Station", kMuonGEM, detid.station());
           holder = GetDaughter(holder, "Layer", kMuonGEM, detid.layer());
           holder = GetDaughter(holder, "Chamber", kMuonGEM, detid.chamber());
 
           AddLeafNode(holder, child, name.c_str(), createPlacement(*it));
         }
       }
     }
   } catch (cms::Exception& exception) {
     edm::LogInfo("FWRecoGeometry") << "failed to produce GEM geometry " << exception.what() << std::endl;
   }
 }

void FWTGeoRecoGeometryESProducer::addHcalCaloGeometryBarrel ( void )

private

Definition at line 792 of file FWTGeoRecoGeometryESProducer.cc.

References c, PVValHelper::dz, dataset::end, CaloCellGeometry::etaPos(), CaloCellGeometry::getCorners(), GetDaughter(), CaloGeometry::getGeometry(), GetMedium(), GetTopHolder(), CaloGeometry::getValidDetIds(), DetId::Hcal, HcalBarrel, HcalDetId::ieta(), kHCal, IdealObliquePrism::localCorners(), m_caloGeom, CaloCellGeometry::param(), CaloCellGeometry::phiPos(), gpuVertexFinder::printf(), makeMuonMisalignmentScenario::rot, x, y, z, and HcalDetId::zside().

Referenced by produce().

                                                                  {
   TGeoVolume* tv = GetTopHolder("HCal", kHCal);
   TGeoVolume* assembly = GetDaughter(tv, "HCalBarrel", kHCal);
 
   std::vector<DetId> vid = m_caloGeom->getValidDetIds(DetId::Hcal, HcalSubdetector::HcalBarrel);
 
   CaloVolMap caloShapeMapP;
   CaloVolMap caloShapeMapN;
   for (std::vector<DetId>::const_iterator it = vid.begin(), end = vid.end(); it != end; ++it) {
     //HcalDetId detid = HcalDetId(it->rawId());
     HcalDetId detid(*it);
     const CaloCellGeometry* cellb = (m_caloGeom->getGeometry(*it)).get();
     const IdealObliquePrism* cell = dynamic_cast<const IdealObliquePrism*>(cellb);
 
     if (!cell) {
       printf("HB not oblique !!!\n");
       continue;
     }
 
     TGeoVolume* volume = nullptr;
     CaloVolMap& caloShapeMap = (cell->etaPos() > 0) ? caloShapeMapP : caloShapeMapN;
     CaloVolMap::iterator volIt = caloShapeMap.find(cell->param());
     if (volIt == caloShapeMap.end()) {
       // printf("FIREWORKS NEW SHAPE BEGIN eta = %f etaPos = %f, phiPos %f >>>>>> \n", cell->eta(), cell->etaPos(), cell->phiPos());
       IdealObliquePrism::Pt3DVec lc(8);
       IdealObliquePrism::Pt3D ref;
       IdealObliquePrism::localCorners(lc, cell->param(), ref);
       HepGeom::Vector3D<float> lCenter;
       for (int c = 0; c < 8; ++c)
         lCenter += lc[c];
       lCenter *= 0.125;
 
       static const int arr[] = {1, 0, 3, 2, 5, 4, 7, 6};
       double points[16];
       for (int c = 0; c < 8; ++c) {
         if (cell->etaPos() > 0)
           points[c * 2 + 0] = -(lc[arr[c]].z() - lCenter.z());
         else
           points[c * 2 + 0] = (lc[arr[c]].z() - lCenter.z());
         points[c * 2 + 1] = (lc[arr[c]].y() - lCenter.y());
         // printf("AMT xy[%d] <=>[%d] = (%.4f, %.4f) \n", arr[c], c, points[c*2],  points[c*2+1]);
       }
 
       float dz = (lc[4].x() - lc[0].x()) * 0.5;
       TGeoShape* solid = new TGeoArb8(dz, &points[0]);
       volume = new TGeoVolume("hcal oblique prism", solid, GetMedium(kHCal));
       caloShapeMap[cell->param()] = volume;
     } else {
       volume = volIt->second;
     }
 
     HepGeom::Vector3D<float> gCenter;
     CaloCellGeometry::CornersVec const& gc = cell->getCorners();
     for (int c = 0; c < 8; ++c)
       gCenter += HepGeom::Vector3D<float>(gc[c].x(), gc[c].y(), gc[c].z());
     gCenter *= 0.125;
 
     TGeoTranslation gtr(gCenter.x(), gCenter.y(), gCenter.z());
     TGeoRotation rot;
     rot.RotateY(90);
 
     TGeoRotation rotPhi;
     rotPhi.SetAngles(0, -cell->phiPos() * TMath::RadToDeg(), 0);
     rot.MultiplyBy(&rotPhi);
 
     TGeoVolume* holder = GetDaughter(assembly, "side", kHCal, detid.zside());
     holder = GetDaughter(holder, "ieta", kHCal, detid.ieta());
     std::stringstream nname;
     nname << detid;
     AddLeafNode(holder, volume, nname.str().c_str(), new TGeoCombiTrans(gtr, rot));
   }
 
   //  printf("HB map size P = %lu , N = %lu", caloShapeMapP.size(),caloShapeMapN.size() );
 }

void FWTGeoRecoGeometryESProducer::addHcalCaloGeometryEndcap ( void )

private

Definition at line 868 of file FWTGeoRecoGeometryESProducer.cc.

References c, PVValHelper::dz, dataset::end, CaloCellGeometry::etaPos(), CaloCellGeometry::getCorners(), GetDaughter(), CaloGeometry::getGeometry(), GetMedium(), GetTopHolder(), CaloGeometry::getValidDetIds(), DetId::Hcal, HcalEndcap, HcalDetId::ieta(), kHCal, IdealObliquePrism::localCorners(), m_caloGeom, CaloCellGeometry::param(), CaloCellGeometry::phiPos(), gpuVertexFinder::printf(), makeMuonMisalignmentScenario::rot, and HcalDetId::zside().

Referenced by produce().

                                                                  {
   CaloVolMap caloShapeMapP;
   CaloVolMap caloShapeMapN;
 
   TGeoVolume* tv = GetTopHolder("HCal", kHCal);
   TGeoVolume* assembly = GetDaughter(tv, "HCalEndcap", kHCal);
 
   std::vector<DetId> vid = m_caloGeom->getValidDetIds(DetId::Hcal, HcalSubdetector::HcalEndcap);
 
   for (std::vector<DetId>::const_iterator it = vid.begin(), end = vid.end(); it != end; ++it) {
     HcalDetId detid = HcalDetId(it->rawId());
     const CaloCellGeometry* cellb = (m_caloGeom->getGeometry(*it)).get();
     const IdealObliquePrism* cell = dynamic_cast<const IdealObliquePrism*>(cellb);
 
     if (!cell) {
       printf("EC not oblique \n");
       continue;
     }
 
     TGeoVolume* volume = nullptr;
     CaloVolMap& caloShapeMap = (cell->etaPos() > 0) ? caloShapeMapP : caloShapeMapN;
     CaloVolMap::iterator volIt = caloShapeMap.find(cell->param());
     if (volIt == caloShapeMap.end()) {
       IdealObliquePrism::Pt3DVec lc(8);
       IdealObliquePrism::Pt3D ref;
       IdealObliquePrism::localCorners(lc, cell->param(), ref);
       HepGeom::Vector3D<float> lCenter;
       for (int c = 0; c < 8; ++c)
         lCenter += lc[c];
       lCenter *= 0.125;
 
       //for( int c = 0; c < 8; ++c)
       //   printf("lc.push_back(TEveVector(%.4f, %.4f, %.4f));\n", lc[c].x(), lc[c].y(), lc[c].z() );
 
       static const int arrP[] = {3, 2, 1, 0, 7, 6, 5, 4};
       static const int arrN[] = {7, 6, 5, 4, 3, 2, 1, 0};
       const int* arr = (detid.ieta() > 0) ? &arrP[0] : &arrN[0];
 
       double points[16];
       for (int c = 0; c < 8; ++c) {
         points[c * 2 + 0] = lc[arr[c]].x() - lCenter.x();
         points[c * 2 + 1] = lc[arr[c]].y() - lCenter.y();
       }
 
       float dz = (lc[4].z() - lc[0].z()) * 0.5;
       TGeoShape* solid = new TGeoArb8(dz, &points[0]);
       volume = new TGeoVolume("ecal oblique prism", solid, GetMedium(kHCal));
       caloShapeMap[cell->param()] = volume;
     } else {
       volume = volIt->second;
     }
 
     HepGeom::Vector3D<float> gCenter;
     CaloCellGeometry::CornersVec const& gc = cell->getCorners();
     for (int c = 0; c < 8; ++c) {
       gCenter += HepGeom::Vector3D<float>(gc[c].x(), gc[c].y(), gc[c].z());
       //  printf("gc.push_back(TEveVector(%.4f, %.4f, %.4f));\n", gc[c].x(), gc[c].y(),gc[c].z() );
     }
     gCenter *= 0.125;
 
     TGeoTranslation gtr(gCenter.x(), gCenter.y(), gCenter.z());
     TGeoRotation rot;
     rot.SetAngles(cell->phiPos() * TMath::RadToDeg(), 0, 0);
 
     TGeoVolume* holder = GetDaughter(assembly, "side", kHCal, detid.zside());
     holder = GetDaughter(holder, "ieta", kHCal, detid.ieta());
     std::stringstream nname;
     nname << detid;
     AddLeafNode(holder, volume, nname.str().c_str(), new TGeoCombiTrans(gtr, rot));
   }
 
   //   printf("HE map size P = %lu , N = %lu", caloShapeMapP.size(),caloShapeMapN.size() );
 }

void FWTGeoRecoGeometryESProducer::addHcalCaloGeometryForward ( )

private

Definition at line 1008 of file FWTGeoRecoGeometryESProducer.cc.

References c, PVValHelper::dz, dataset::end, CaloCellGeometry::etaPos(), CaloCellGeometry::getCorners(), GetDaughter(), CaloGeometry::getGeometry(), GetMedium(), GetTopHolder(), CaloGeometry::getValidDetIds(), DetId::Hcal, HcalForward, HcalDetId::ieta(), kHCal, IdealZPrism::localCorners(), m_caloGeom, CaloCellGeometry::param(), CaloCellGeometry::phiPos(), gpuVertexFinder::printf(), makeMuonMisalignmentScenario::rot, and HcalDetId::zside().

Referenced by produce().

                                                               {
   CaloVolMap caloShapeMapP;
   CaloVolMap caloShapeMapN;
 
   TGeoVolume* tv = GetTopHolder("HCal", kHCal);
   TGeoVolume* assembly = GetDaughter(tv, "HCalForward", kHCal);
 
   std::vector<DetId> vid = m_caloGeom->getValidDetIds(DetId::Hcal, HcalSubdetector::HcalForward);
 
   for (std::vector<DetId>::const_iterator it = vid.begin(), end = vid.end(); it != end; ++it) {
     HcalDetId detid = HcalDetId(it->rawId());
     const CaloCellGeometry* cellb = (m_caloGeom->getGeometry(*it)).get();
     const IdealZPrism* cell = dynamic_cast<const IdealZPrism*>(cellb);
 
     if (!cell) {
       printf("EC not Z prism \n");
       continue;
     }
 
     TGeoVolume* volume = nullptr;
     CaloVolMap& caloShapeMap = (cell->etaPos() > 0) ? caloShapeMapP : caloShapeMapN;
     CaloVolMap::iterator volIt = caloShapeMap.find(cell->param());
     if (volIt == caloShapeMap.end()) {
       IdealZPrism::Pt3DVec lc(8);
       IdealZPrism::Pt3D ref;
       IdealZPrism::localCorners(lc, cell->param(), ref);
       HepGeom::Vector3D<float> lCenter;
       for (int c = 0; c < 8; ++c)
         lCenter += lc[c];
       lCenter *= 0.125;
       static const int arrP[] = {3, 2, 1, 0, 7, 6, 5, 4};
       static const int arrN[] = {7, 6, 5, 4, 3, 2, 1, 0};
       const int* arr = (detid.ieta() > 0) ? &arrP[0] : &arrN[0];
 
       double points[16];
       for (int c = 0; c < 8; ++c) {
         points[c * 2 + 0] = lc[arr[c]].x() - lCenter.x();
         points[c * 2 + 1] = lc[arr[c]].y() - lCenter.y();
       }
 
       float dz = (lc[4].z() - lc[0].z()) * 0.5;
       TGeoShape* solid = new TGeoArb8(dz, &points[0]);
       volume = new TGeoVolume("ecal oblique prism", solid, GetMedium(kHCal));
       caloShapeMap[cell->param()] = volume;
     } else {
       volume = volIt->second;
     }
     HepGeom::Vector3D<float> gCenter;
     CaloCellGeometry::CornersVec const& gc = cell->getCorners();
     for (int c = 0; c < 8; ++c) {
       gCenter += HepGeom::Vector3D<float>(gc[c].x(), gc[c].y(), gc[c].z());
     }
     gCenter *= 0.125;
 
     TGeoTranslation gtr(gCenter.x(), gCenter.y(), gCenter.z());
     TGeoRotation rot;
     rot.SetAngles(cell->phiPos() * TMath::RadToDeg(), 0, 0);
 
     TGeoVolume* holder = GetDaughter(assembly, "side", kHCal, detid.zside());
     holder = GetDaughter(holder, "ieta", kHCal, detid.ieta());
     std::stringstream nname;
     nname << detid;
     AddLeafNode(holder, volume, nname.str().c_str(), new TGeoCombiTrans(gtr, rot));
   }
 }

void FWTGeoRecoGeometryESProducer::addHcalCaloGeometryOuter ( )

private

Definition at line 942 of file FWTGeoRecoGeometryESProducer.cc.

References c, PVValHelper::dz, dataset::end, CaloCellGeometry::etaPos(), CaloCellGeometry::getCorners(), GetDaughter(), CaloGeometry::getGeometry(), GetMedium(), GetTopHolder(), CaloGeometry::getValidDetIds(), DetId::Hcal, HcalOuter, HcalDetId::ieta(), kHCal, IdealObliquePrism::localCorners(), m_caloGeom, CaloCellGeometry::param(), CaloCellGeometry::phiPos(), gpuVertexFinder::printf(), makeMuonMisalignmentScenario::rot, and HcalDetId::zside().

Referenced by produce().

                                                             {
   CaloVolMap caloShapeMapP;
   CaloVolMap caloShapeMapN;
 
   TGeoVolume* tv = GetTopHolder("HCal", kHCal);
   TGeoVolume* assembly = GetDaughter(tv, "HCalOuter", kHCal);
 
   std::vector<DetId> vid = m_caloGeom->getValidDetIds(DetId::Hcal, HcalSubdetector::HcalOuter);
 
   for (std::vector<DetId>::const_iterator it = vid.begin(), end = vid.end(); it != end; ++it) {
     HcalDetId detid = HcalDetId(it->rawId());
     const CaloCellGeometry* cellb = (m_caloGeom->getGeometry(*it)).get();
     const IdealObliquePrism* cell = dynamic_cast<const IdealObliquePrism*>(cellb);
 
     if (!cell) {
       printf("EC not oblique \n");
       continue;
     }
 
     TGeoVolume* volume = nullptr;
     CaloVolMap& caloShapeMap = (cell->etaPos() > 0) ? caloShapeMapP : caloShapeMapN;
     CaloVolMap::iterator volIt = caloShapeMap.find(cell->param());
     if (volIt == caloShapeMap.end()) {
       IdealObliquePrism::Pt3DVec lc(8);
       IdealObliquePrism::Pt3D ref;
       IdealObliquePrism::localCorners(lc, cell->param(), ref);
       HepGeom::Vector3D<float> lCenter;
       for (int c = 0; c < 8; ++c)
         lCenter += lc[c];
       lCenter *= 0.125;
       static const int arrP[] = {3, 2, 1, 0, 7, 6, 5, 4};
       static const int arrN[] = {7, 6, 5, 4, 3, 2, 1, 0};
       const int* arr = (detid.ieta() > 0) ? &arrP[0] : &arrN[0];
 
       double points[16];
       for (int c = 0; c < 8; ++c) {
         points[c * 2 + 0] = lc[arr[c]].x() - lCenter.x();
         points[c * 2 + 1] = lc[arr[c]].y() - lCenter.y();
       }
 
       float dz = (lc[4].z() - lc[0].z()) * 0.5;
       TGeoShape* solid = new TGeoArb8(dz, &points[0]);
       volume = new TGeoVolume("ecal oblique prism", solid, GetMedium(kHCal));
       caloShapeMap[cell->param()] = volume;
     } else {
       volume = volIt->second;
     }
     HepGeom::Vector3D<float> gCenter;
     CaloCellGeometry::CornersVec const& gc = cell->getCorners();
     for (int c = 0; c < 8; ++c) {
       gCenter += HepGeom::Vector3D<float>(gc[c].x(), gc[c].y(), gc[c].z());
     }
     gCenter *= 0.125;
 
     TGeoTranslation gtr(gCenter.x(), gCenter.y(), gCenter.z());
     TGeoRotation rot;
     rot.SetAngles(cell->phiPos() * TMath::RadToDeg(), 0, 0);
 
     TGeoVolume* holder = GetDaughter(assembly, "side", kHCal, detid.zside());
     holder = GetDaughter(holder, "ieta", kHCal, detid.ieta());
     std::stringstream nname;
     nname << detid;
     AddLeafNode(holder, volume, nname.str().c_str(), new TGeoCombiTrans(gtr, rot));
   }
 }

void FWTGeoRecoGeometryESProducer::addME0Geometry ( )

private

Definition at line 758 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), ME0Geometry::etaPartitions(), cppFunctionSkipper::exception, GetDaughter(), GetTopHolder(), kMuonME0, m_trackingGeom, DetId::Muon, alignCSCRings::s, GlobalTrackingGeometry::slaveGeometry(), AlCaHLTBitMon_QueryRunRegistry::string, and cms::Exception::what().

Referenced by produce().

                                                   {
   TGeoVolume* tv = GetTopHolder("Muon", kMuonME0);
   TGeoVolume* assembly = GetDaughter(tv, "ME0", kMuonME0);
 
   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();
         //std::cout << "AMT FWTTTTRecoGeometryES\n" << rawid ;
 
         ME0DetId detid(rawid);
         std::stringstream s;
         s << detid;
         std::string name = s.str();
         TGeoVolume* child = createVolume(name, roll, kMuonME0);
 
         TGeoVolume* holder = GetDaughter(assembly, "Region", kMuonME0, detid.region());
         holder = GetDaughter(holder, "Layer", kMuonME0, detid.layer());
         holder = GetDaughter(holder, "Chamber", kMuonME0, detid.chamber());
         AddLeafNode(holder, child, name.c_str(), createPlacement(roll));
       }
     }
   } catch (cms::Exception& exception) {
     edm::LogInfo("FWRecoGeometry") << "failed to produce ME0 geometry " << exception.what() << std::endl;
   }
 }

void FWTGeoRecoGeometryESProducer::addPixelBarrelGeometry ( )

private

Definition at line 432 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), TrackerGeometry::detsPXB(), GetDaughter(), GetTopHolder(), kSiPixel, PVValHelper::ladder, phase1PixelTopology::layer, m_trackerGeom, m_trackerTopology, callgraph::module, TrackerTopology::pxbLadder(), TrackerTopology::pxbLayer(), TrackerTopology::pxbModule(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

                                                           {
   TGeoVolume* tv = GetTopHolder("SiPixel", kSiPixel);
   TGeoVolume* assembly = GetDaughter(tv, "PXB", kSiPixel);
 
   for (auto it : m_trackerGeom->detsPXB()) {
     DetId detid = it->geographicalId();
     unsigned int layer = m_trackerTopology->pxbLayer(detid);
     unsigned int module = m_trackerTopology->pxbModule(detid);
     unsigned int ladder = m_trackerTopology->pxbLadder(detid);
 
     std::string name = Form("PXB Ly:%d, Md:%d Ld:%d ", layer, module, ladder);
     TGeoVolume* child = createVolume(name, it, kSiPixel);
 
     TGeoVolume* holder = GetDaughter(assembly, "Layer", kSiPixel, layer);
     holder = GetDaughter(holder, "Module", kSiPixel, module);
 
     AddLeafNode(holder, child, name.c_str(), createPlacement(it));
   }
 }

void FWTGeoRecoGeometryESProducer::addPixelForwardGeometry ( )

private

Definition at line 452 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), TrackerGeometry::detsPXF(), GetDaughter(), GetTopHolder(), kSiPixel, m_trackerGeom, m_trackerTopology, TrackerTopology::pxfBlade(), TrackerTopology::pxfDisk(), TrackerTopology::pxfPanel(), TrackerTopology::side(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

                                                            {
   TGeoVolume* tv = GetTopHolder("SiPixel", kSiPixel);
   TGeoVolume* assembly = GetDaughter(tv, "PXF", kSiPixel);
 
   for (auto it : m_trackerGeom->detsPXF()) {
     DetId detid = it->geographicalId();
     unsigned int disk = m_trackerTopology->pxfDisk(detid);
     unsigned int blade = m_trackerTopology->pxfBlade(detid);
     unsigned int panel = m_trackerTopology->pxfPanel(detid);
     unsigned int side = m_trackerTopology->side(detid);
 
     std::string name = Form("PXF D:%d, B:%d, P:%d, S:%d", disk, blade, panel, side);
     TGeoVolume* child = createVolume(name, it, kSiPixel);
 
     TGeoVolume* holder = GetDaughter(assembly, "Side", kSiPixel, side);
     holder = GetDaughter(holder, "Disk", kSiPixel, disk);
     holder = GetDaughter(holder, "Blade", kSiPixel, blade);
     holder = GetDaughter(holder, "Panel", kSiPixel, panel);
 
     AddLeafNode(holder, child, name.c_str(), createPlacement(it));
   }
 }

void FWTGeoRecoGeometryESProducer::addRPCGeometry ( )

private

Definition at line 730 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), dataset::end, GeomDet::geographicalId(), GetDaughter(), GetTopHolder(), kMuonRPC, m_trackingGeom, DetId::Muon, RPCGeometry::rolls(), MuonSubdetId::RPC, alignCSCRings::s, GlobalTrackingGeometry::slaveGeometry(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

                                                   {
   TGeoVolume* tv = GetTopHolder("Muon", kMuonRPC);
   TGeoVolume* assembly = GetDaughter(tv, "RPC", kMuonRPC);
 
   DetId detId(DetId::Muon, MuonSubdetId::RPC);
   const RPCGeometry* rpcGeom = static_cast<const RPCGeometry*>(m_trackingGeom->slaveGeometry(detId));
   for (auto it = rpcGeom->rolls().begin(), end = rpcGeom->rolls().end(); it != end; ++it) {
     RPCRoll const* roll = (*it);
     if (roll) {
       RPCDetId detid = roll->geographicalId();
       std::stringstream s;
       s << detid;
       std::string name = s.str();
 
       TGeoVolume* child = createVolume(name, roll, kMuonRPC);
 
       TGeoVolume* holder = GetDaughter(assembly, "ROLL Region", kMuonRPC, detid.region());
       holder = GetDaughter(holder, "Ring", kMuonRPC, detid.ring());
       holder = GetDaughter(holder, "Station", kMuonRPC, detid.station());
       holder = GetDaughter(holder, "Sector", kMuonRPC, detid.sector());
       holder = GetDaughter(holder, "Layer", kMuonRPC, detid.layer());
       holder = GetDaughter(holder, "Subsector", kMuonRPC, detid.subsector());
 
       AddLeafNode(holder, child, name.c_str(), createPlacement(*it));
     }
   };
 }

void FWTGeoRecoGeometryESProducer::addTECGeometry ( )

private

Definition at line 535 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), TrackerGeometry::detsTEC(), GetDaughter(), GetTopHolder(), kSiStrip, m_trackerGeom, m_trackerTopology, callgraph::module, TrackerTopology::print(), relativeConstraints::ring, AlCaHLTBitMon_QueryRunRegistry::string, TrackerTopology::tecModule(), TrackerTopology::tecOrder(), and TrackerTopology::tecRing().

Referenced by produce().

                                                   {
   TGeoVolume* tv = GetTopHolder("SiStrip", kSiStrip);
   TGeoVolume* assembly = GetDaughter(tv, "TEC", kSiStrip);
 
   for (auto it : m_trackerGeom->detsTEC()) {
     DetId detid = it->geographicalId();
     unsigned int order = m_trackerTopology->tecOrder(detid);
     unsigned int ring = m_trackerTopology->tecRing(detid);
     unsigned int module = m_trackerTopology->tecModule(detid);
 
     std::string name = m_trackerTopology->print(detid);
 
     TGeoVolume* child = createVolume(name, it, kSiStrip);
 
     TGeoVolume* holder = GetDaughter(assembly, "Order", kSiStrip, order);
     holder = GetDaughter(holder, "Ring", kSiStrip, ring);
     holder = GetDaughter(holder, "Module", kSiStrip, module);
     AddLeafNode(holder, child, name.c_str(), createPlacement(it));
   }
 }

void FWTGeoRecoGeometryESProducer::addTIBGeometry ( )

private

Definition at line 475 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), TrackerGeometry::detsTIB(), GetDaughter(), GetTopHolder(), kSiStrip, m_trackerGeom, m_trackerTopology, callgraph::module, TrackerTopology::print(), AlCaHLTBitMon_QueryRunRegistry::string, TrackerTopology::tibModule(), TrackerTopology::tibOrder(), and TrackerTopology::tibSide().

Referenced by produce().

                                                   {
   TGeoVolume* tv = GetTopHolder("SiStrip", kSiStrip);
   TGeoVolume* assembly = GetDaughter(tv, "TIB", kSiStrip);
 
   for (auto it : m_trackerGeom->detsTIB()) {
     DetId detid = it->geographicalId();
     unsigned int module = m_trackerTopology->tibModule(detid);
     unsigned int order = m_trackerTopology->tibOrder(detid);
     unsigned int side = m_trackerTopology->tibSide(detid);
 
     std::string name = m_trackerTopology->print(detid);
 
     TGeoVolume* child = createVolume(name, it, kSiStrip);
     TGeoVolume* holder = GetDaughter(assembly, "Module", kSiStrip, module);
     holder = GetDaughter(holder, "Order", kSiStrip, order);
     holder = GetDaughter(holder, "Side", kSiStrip, side);
     AddLeafNode(holder, child, name.c_str(), createPlacement(it));
   }
 }

void FWTGeoRecoGeometryESProducer::addTIDGeometry ( )

private

Definition at line 495 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), TrackerGeometry::detsTID(), GetDaughter(), GetTopHolder(), kSiStrip, m_trackerGeom, m_trackerTopology, TrackerTopology::print(), relativeConstraints::ring, AlCaHLTBitMon_QueryRunRegistry::string, TrackerTopology::tidRing(), TrackerTopology::tidSide(), and TrackerTopology::tidWheel().

Referenced by produce().

                                                   {
   TGeoVolume* tv = GetTopHolder("SiStrip", kSiStrip);
   TGeoVolume* assembly = GetDaughter(tv, "TID", kSiStrip);
 
   for (auto it : m_trackerGeom->detsTID()) {
     DetId detid = it->geographicalId();
     unsigned int side = m_trackerTopology->tidSide(detid);
     unsigned int wheel = m_trackerTopology->tidWheel(detid);
     unsigned int ring = m_trackerTopology->tidRing(detid);
 
     std::string name = m_trackerTopology->print(detid);
 
     TGeoVolume* child = createVolume(name, it, kSiStrip);
     TGeoVolume* holder = GetDaughter(assembly, "Side", kSiStrip, side);
     holder = GetDaughter(holder, "Wheel", kSiStrip, wheel);
     holder = GetDaughter(holder, "Ring", kSiStrip, ring);
     AddLeafNode(holder, child, name.c_str(), createPlacement(it));
   }
 }

void FWTGeoRecoGeometryESProducer::addTOBGeometry ( )

private

Definition at line 515 of file FWTGeoRecoGeometryESProducer.cc.

References dd4hep::createPlacement(), createVolume(), TrackerGeometry::detsTOB(), GetDaughter(), GetTopHolder(), kSiStrip, m_trackerGeom, m_trackerTopology, callgraph::module, TrackerTopology::print(), AlCaHLTBitMon_QueryRunRegistry::string, TrackerTopology::tobModule(), TrackerTopology::tobRod(), and TrackerTopology::tobSide().

Referenced by produce().

                                                   {
   TGeoVolume* tv = GetTopHolder("SiStrip", kSiStrip);
   TGeoVolume* assembly = GetDaughter(tv, "TOB", kSiStrip);
 
   for (auto it : m_trackerGeom->detsTOB()) {
     DetId detid = it->geographicalId();
     unsigned int rod = m_trackerTopology->tobRod(detid);
     unsigned int side = m_trackerTopology->tobSide(detid);
     unsigned int module = m_trackerTopology->tobModule(detid);
 
     std::string name = m_trackerTopology->print(detid);
 
     TGeoVolume* child = createVolume(name, it, kSiStrip);
     TGeoVolume* holder = GetDaughter(assembly, "Rod", kSiStrip, rod);
     holder = GetDaughter(holder, "Side", kSiStrip, side);
     holder = GetDaughter(holder, "Module", kSiStrip, module);
     AddLeafNode(holder, child, name.c_str(), createPlacement(it));
   }
 }

TGeoManager* FWTGeoRecoGeometryESProducer::createManager ( int level )

private

TGeoShape * FWTGeoRecoGeometryESProducer::createShape ( const GeomDet * det )

private

Create TGeo shape for GeomDet

Definition at line 346 of file FWTGeoRecoGeometryESProducer.cc.

References b, b2, Surface::bounds(), Bounds::length(), m_nameToShape, mergeVDriftHistosByStation::name, TrapezoidalPlaneBounds::parameters(), alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, GeomDet::surface(), Bounds::thickness(), TrackerMaterial_cfi::thickness, and Bounds::width().

Referenced by createVolume().

                                                                        {
   TGeoShape* shape = nullptr;
 
   // Trapezoidal
   const Bounds* b = &((det->surface()).bounds());
   const TrapezoidalPlaneBounds* b2 = dynamic_cast<const TrapezoidalPlaneBounds*>(b);
   if (b2) {
     std::array<const float, 4> const& par = b2->parameters();
 
     // These parameters are half-lengths, as in CMSIM/GEANT3
     float hBottomEdge = par[0];
     float hTopEdge = par[1];
     float thickness = par[2];
     float apothem = par[3];
 
     std::stringstream s;
     s << "T_" << hBottomEdge << "_" << hTopEdge << "_" << thickness << "_" << apothem;
     std::string name = s.str();
 
     // Do not create identical shape,
     // if one already exists
     shape = m_nameToShape[name];
     if (nullptr == shape) {
       shape = new TGeoTrap(name.c_str(),
                            thickness,    //dz
                            0,            //theta
                            0,            //phi
                            apothem,      //dy1
                            hBottomEdge,  //dx1
                            hTopEdge,     //dx2
                            0,            //alpha1
                            apothem,      //dy2
                            hBottomEdge,  //dx3
                            hTopEdge,     //dx4
                            0);           //alpha2
 
       m_nameToShape[name] = shape;
     }
   }
   if (dynamic_cast<const RectangularPlaneBounds*>(b) != nullptr) {
     // Rectangular
     float length = det->surface().bounds().length();
     float width = det->surface().bounds().width();
     float thickness = det->surface().bounds().thickness();
 
     std::stringstream s;
     s << "R_" << width << "_" << length << "_" << thickness;
     std::string name = s.str();
 
     // Do not create identical shape,
     // if one already exists
     shape = m_nameToShape[name];
     if (nullptr == shape) {
       shape = new TGeoBBox(name.c_str(), width / 2., length / 2., thickness / 2.);  // dx, dy, dz
 
       m_nameToShape[name] = shape;
     }
   }
 
   return shape;
 }

TGeoVolume * FWTGeoRecoGeometryESProducer::createVolume	(	const std::string &	name,
		const GeomDet *	det,
		ERecoDet	mid = `kDummy`
	)

private

Create TGeo volume for GeomDet

Definition at line 409 of file FWTGeoRecoGeometryESProducer.cc.

References createShape(), GetMedium(), and m_shapeToVolume.

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

                                                                                                               {
   TGeoShape* solid = createShape(det);
 
   std::map<TGeoShape*, TGeoVolume*>::iterator vIt = m_shapeToVolume.find(solid);
   if (vIt != m_shapeToVolume.end())
     return vIt->second;
 
   TGeoVolume* volume = new TGeoVolume(name.c_str(), solid, GetMedium(mid));
 
   m_shapeToVolume[solid] = volume;
 
   return volume;
 }

TGeoVolume * FWTGeoRecoGeometryESProducer::GetDaughter	(	TGeoVolume *	mother,
		const char *	prefix,
		ERecoDet	cidx,
		int	id
	)

private

Definition at line 125 of file FWTGeoRecoGeometryESProducer.cc.

References GetMedium().

Referenced by addCaloTowerGeometry(), addCSCGeometry(), addDTGeometry(), addEcalCaloGeometry(), addGEMGeometry(), addHcalCaloGeometryBarrel(), addHcalCaloGeometryEndcap(), addHcalCaloGeometryForward(), addHcalCaloGeometryOuter(), addME0Geometry(), addPixelBarrelGeometry(), addPixelForwardGeometry(), addRPCGeometry(), addTECGeometry(), addTIBGeometry(), addTIDGeometry(), addTOBGeometry(), and GetTopHolder().

                                                                                                                    {
   TGeoVolume* res = nullptr;
   if (mother->GetNdaughters()) {
     TGeoNode* n = mother->FindNode(Form("%s_%d_1", prefix, id));
     if (n)
       res = n->GetVolume();
   }
 
   if (!res) {
     res = new TGeoVolumeAssembly(Form("%s_%d", prefix, id));
     res->SetMedium(GetMedium(cidx));
     mother->AddNode(res, 1);
   }
 
   return res;
 }

TGeoVolume * FWTGeoRecoGeometryESProducer::GetDaughter	(	TGeoVolume *	mother,
		const char *	prefix,
		ERecoDet	cidx
	)

private

Definition at line 142 of file FWTGeoRecoGeometryESProducer.cc.

References GetMedium().

                                                                                                            {
   TGeoVolume* res = nullptr;
   if (mother->GetNdaughters()) {
     TGeoNode* n = mother->FindNode(Form("%s_1", prefix));
     if (n)
       res = n->GetVolume();
   }
 
   if (!res) {
     //      printf("GetDau... new holder %s for mother %s \n", mother->GetName(), prefix);
     res = new TGeoVolumeAssembly(prefix);
     res->SetMedium(GetMedium(cidx));
     mother->AddNode(res, 1);
   }
 
   return res;
 }

TGeoMedium * FWTGeoRecoGeometryESProducer::GetMedium ( ERecoDet det )

private

Definition at line 185 of file FWTGeoRecoGeometryESProducer.cc.

References kCaloTower, kECal, kHCal, kHGCE, kHGCH, kMuonCSC, kMuonDT, kMuonGEM, kMuonME0, kMuonRPC, kSiPixel, kSiStrip, m_dummyMedium, m_recoMedium, gpuVertexFinder::printf(), AlCaHLTBitMon_QueryRunRegistry::string, and Phase2TrackerValidateDigi_cff::switch.

Referenced by addCaloTowerGeometry(), addEcalCaloGeometry(), addHcalCaloGeometryBarrel(), addHcalCaloGeometryEndcap(), addHcalCaloGeometryForward(), addHcalCaloGeometryOuter(), createVolume(), and GetDaughter().

                                                                 {
   std::map<ERecoDet, TGeoMedium*>::iterator it = m_recoMedium.find(det);
   if (it != m_recoMedium.end())
     return it->second;
 
   std::string name;
   int color;
 
   switch (det) {
     // TRACKER
     case kSiPixel:
       name = "SiPixel";
       color = GMCol::Green;
       break;
 
     case kSiStrip:
       name = "SiStrip";
       color = GMCol::Gray;
       break;
       // MUON
     case kMuonDT:
       name = "MuonDT";
       color = GMCol::Blue2;
       break;
 
     case kMuonRPC:
       name = "MuonRPC";
       color = GMCol::Red;
       break;
 
     case kMuonGEM:
       name = "MuonGEM";
       color = GMCol::Yellow1;
       break;
 
     case kMuonCSC:
       name = "MuonCSC";
       color = GMCol::Gray;
       break;
 
     case kMuonME0:
       name = "MuonME0";
       color = GMCol::Yellow0;
       break;
 
       // CALO
     case kECal:
       name = "ECal";
       color = GMCol::Blue2;
       break;
     case kHCal:
       name = "HCal";
       color = GMCol::Orange1;
       break;
     case kCaloTower:
       name = "CaloTower";
       color = GMCol::Green;
       break;
     case kHGCE:
       name = "HGCEE";
       color = GMCol::Blue2;
       break;
     case kHGCH:
       name = "HGCEH";
       color = GMCol::Blue1;
       break;
     default:
       printf("invalid medium id \n");
       return m_dummyMedium;
   }
 
   TGeoMaterial* mat = new TGeoMaterial(name.c_str(), 0, 0, 0);
   mat->SetZ(color);
   m_recoMedium[det] = new TGeoMedium(name.c_str(), 0, mat);
   mat->SetFillStyle(3000);  // tansparency 3000-3100
   mat->SetDensity(1);       // disable override of transparency in TGeoManager::DefaultColors()
 
   return m_recoMedium[det];
 }

TGeoVolume * FWTGeoRecoGeometryESProducer::GetTopHolder	(	const char *	prefix,
		ERecoDet	cidx
	)

private

Definition at line 160 of file FWTGeoRecoGeometryESProducer.cc.

References GetDaughter(), and PostProcessorHGCAL_cfi::prefix.

Referenced by addCaloTowerGeometry(), addCSCGeometry(), addDTGeometry(), addEcalCaloGeometry(), addGEMGeometry(), addHcalCaloGeometryBarrel(), addHcalCaloGeometryEndcap(), addHcalCaloGeometryForward(), addHcalCaloGeometryOuter(), addME0Geometry(), addPixelBarrelGeometry(), addPixelForwardGeometry(), addRPCGeometry(), addTECGeometry(), addTIBGeometry(), addTIDGeometry(), and addTOBGeometry().

                                                                                         {
   //   printf("GetTopHolder res = %s \n", prefix);
   TGeoVolume* res = GetDaughter(gGeoManager->GetTopVolume(), prefix, cidx);
   return res;
 }

const FWTGeoRecoGeometryESProducer& FWTGeoRecoGeometryESProducer::operator= ( const FWTGeoRecoGeometryESProducer & )

private

std::unique_ptr< FWTGeoRecoGeometry > FWTGeoRecoGeometryESProducer::produce ( const FWTGeoRecoGeometryRecord & record )

Definition at line 267 of file FWTGeoRecoGeometryESProducer.cc.

References addCaloTowerGeometry(), addCSCGeometry(), addDTGeometry(), addEcalCaloGeometry(), addGEMGeometry(), addHcalCaloGeometryBarrel(), addHcalCaloGeometryEndcap(), addHcalCaloGeometryForward(), addHcalCaloGeometryOuter(), addME0Geometry(), addPixelBarrelGeometry(), addPixelForwardGeometry(), addRPCGeometry(), addTECGeometry(), addTIBGeometry(), addTIDGeometry(), addTOBGeometry(), relativeConstraints::geom, edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::get(), m_calo, m_caloGeom, m_caloGeomToken, m_dummyMedium, m_gem, m_muon, m_tracker, m_trackerGeom, m_trackerTopology, m_trackerTopologyToken, m_trackingGeom, m_trackingGeomToken, GlobalTrackingGeometry::slaveGeometry(), and DetId::Tracker.

                                                                                                               {
   using namespace edm;
 
   auto fwTGeoRecoGeometry = std::make_unique<FWTGeoRecoGeometry>();
 
   if (m_calo) {
     m_caloGeom = &record.get(m_caloGeomToken);
   }
 
   TGeoManager* geom = new TGeoManager("cmsGeo", "CMS Detector");
   if (nullptr == gGeoIdentity) {
     gGeoIdentity = new TGeoIdentity("Identity");
   }
 
   fwTGeoRecoGeometry->manager(geom);
 
   // Default material is Vacuum
   TGeoMaterial* vacuum = new TGeoMaterial("Vacuum", 0, 0, 0);
   m_dummyMedium = new TGeoMedium("reco", 0, vacuum);
 
   TGeoVolume* top = geom->MakeBox("CMS", m_dummyMedium, 270., 270., 120.);
 
   if (nullptr == top) {
     return std::unique_ptr<FWTGeoRecoGeometry>();
   }
   geom->SetTopVolume(top);
   // ROOT chokes unless colors are assigned
   top->SetVisibility(kFALSE);
   top->SetLineColor(kBlue);
 
   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));
 
     m_trackerTopology = &record.get(m_trackerTopologyToken);
 
     addPixelBarrelGeometry();
     addPixelForwardGeometry();
 
     addTIBGeometry();
     addTIDGeometry();
     addTOBGeometry();
     addTECGeometry();
   }
 
   if (m_muon) {
     addDTGeometry();
     addCSCGeometry();
     addRPCGeometry();
     addME0Geometry();
   }
 
   if (m_gem) {
     addGEMGeometry();
   }
 
   if (m_calo) {
     addEcalCaloGeometry();
     addHcalCaloGeometryBarrel();
     addHcalCaloGeometryEndcap();
     addHcalCaloGeometryOuter();
     addHcalCaloGeometryForward();
     addCaloTowerGeometry();
   }
 
   geom->CloseGeometry();
 
   geom->DefaultColors();
   // printf("==== geo manager NNodes = %d \n", geom->GetNNodes());
   geom->CloseGeometry();
 
   return fwTGeoRecoGeometry;
 }