#include <FWTGeoRecoGeometryESProducer.h>

Inheritance diagram for FWTGeoRecoGeometryESProducer:

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

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

	~FWTGeoRecoGeometryESProducer (void) override

Public Member Functions inherited from edm::ESProducer
	ESProducer ()

	~ESProducer ()(false) override

Public Member Functions inherited from edm::ESProxyFactoryProducer
	ESProxyFactoryProducer ()

void	newInterval (const eventsetup::EventSetupRecordKey &iRecordType, const ValidityInterval &iInterval) override
	overrides DataProxyProvider method More...

	~ESProxyFactoryProducer () noexcept(false) override

Public Member Functions inherited from edm::eventsetup::DataProxyProvider
	DataProxyProvider ()

const ComponentDescription &	description () const

bool	isUsingRecord (const EventSetupRecordKey &) const

const KeyedProxies &	keyedProxies (const EventSetupRecordKey &iRecordKey) const

void	resetProxies (const EventSetupRecordKey &iRecordType)

void	resetProxiesIfTransient (const EventSetupRecordKey &iRecordType)

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

TGeoMedium *	m_dummyMedium

std::shared_ptr< FWTGeoRecoGeometry >	m_fwGeometry

edm::ESHandle< GlobalTrackingGeometry >	m_geomRecord

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

const TrackerTopology *	m_trackerTopology

Additional Inherited Members
Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair< DataKey, edm::propagate_const< std::shared_ptr< DataProxy > > > >	KeyedProxies

typedef std::vector< EventSetupRecordKey >	Keys

typedef std::map< EventSetupRecordKey, KeyedProxies >	RecordProxies

Static Public Member Functions inherited from edm::eventsetup::DataProxyProvider
static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::ESProducer
	ESProducer (const ESProducer &)=delete

ESProducer const &	operator= (const ESProducer &)=delete

template<typename T >
void	setWhatProduced (T *iThis, const es::Label &iLabel=es::Label())

template<typename T >
void	setWhatProduced (T iThis, const char iLabel)

template<typename T >
void	setWhatProduced (T *iThis, const std::string &iLabel)

template<typename T , typename TDecorator >
void	setWhatProduced (T *iThis, const TDecorator &iDec, const es::Label &iLabel=es::Label())

template<typename T , typename TReturn , typename TRecord >
void	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const es::Label &iLabel=es::Label())

template<typename T , typename TReturn , typename TRecord , typename TArg >
void	setWhatProduced (T iThis, TReturn(T::iMethod)(const TRecord &), const TArg &iDec, const es::Label &iLabel=es::Label())

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 eventsetup::EventSetupRecordKey &iRecord, std::unique_ptr< eventsetup::ProxyFactoryBase > iFactory, const std::string &iLabel=std::string())

void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList) override
	override DataProxyProvider method More...

Protected Member Functions inherited from edm::eventsetup::DataProxyProvider
void	eraseAll (const EventSetupRecordKey &iRecordKey)
	deletes all the Proxies in aStream More...

void	invalidateProxies (const EventSetupRecordKey &iRecordKey)

template<class T >
void	usingRecord ()

void	usingRecordWithKey (const EventSetupRecordKey &)

Detailed Description

Definition at line 29 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 31 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 64 of file FWTGeoRecoGeometryESProducer.cc.

References edm::ParameterSet::getUntrackedParameter(), m_calo, m_muon, m_tracker, edm::ESProducer::setWhatProduced(), and funct::true.

   : m_dummyMedium(nullptr)
 {
   m_tracker = pset.getUntrackedParameter<bool>( "Tracker", true );
   m_muon = pset.getUntrackedParameter<bool>( "Muon", true );
   m_calo = pset.getUntrackedParameter<bool>( "Calo", true );
   
   setWhatProduced( this );
 }

FWTGeoRecoGeometryESProducer::~FWTGeoRecoGeometryESProducer ( void )

override

FWTGeoRecoGeometryESProducer::FWTGeoRecoGeometryESProducer ( const FWTGeoRecoGeometryESProducer & )

private

Member Function Documentation

void FWTGeoRecoGeometryESProducer::addCaloTowerGeometry ( )

private

Definition at line 1173 of file FWTGeoRecoGeometryESProducer.cc.

References EnergyCorrector::c, DetId::Calo, PVValHelper::dz, end, CaloCellGeometry::etaPos(), CaloCellGeometry::getCorners(), GetDaughter(), CaloGeometry::getGeometry(), GetMedium(), GetTopHolder(), CaloGeometry::getValidDetIds(), CaloTowerDetId::ieta(), kCaloTower, IdealObliquePrism::localCorners(), m_caloGeom, CaloCellGeometry::param(), CaloCellGeometry::phiPos(), hiPixelPairStep_cff::points, 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 ( void )

private

Definition at line 686 of file FWTGeoRecoGeometryESProducer.cc.

References relativeConstraints::chamber, createVolume(), TrackingGeometry::dets(), GetDaughter(), GetTopHolder(), kMuonCSC, kMuonRPC, m_geomRecord, dataset::name, alignCSCRings::s, GlobalTrackingGeometry::slaveGeometry(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    if(! m_geomRecord->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_geomRecord->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 ( void )

private

Definition at line 597 of file FWTGeoRecoGeometryESProducer.cc.

References relativeConstraints::chamber, createVolume(), TrackingGeometry::dets(), end, GetDaughter(), GetTopHolder(), kMuonDT, kMuonRPC, m_geomRecord, dataset::name, 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_geomRecord->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_geomRecord->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_geomRecord->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 1310 of file FWTGeoRecoGeometryESProducer.cc.

References DetId::Ecal, EcalBarrel, EcalEndcap, end, CaloCellGeometry::getCorners(), GetDaughter(), getEcalTrans(), CaloGeometry::getGeometry(), GetMedium(), GetTopHolder(), CaloGeometry::getValidDetIds(), EBDetId::ieta(), EEDetId::ix(), kECal, m_caloGeom, makeEcalShape(), CaloCellGeometry::param(), 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 ( void )

private

Definition at line 729 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by produce().

 { 
    try {
       DetId detId( DetId::Muon, MuonSubdetId::GEM );
       const GEMGeometry* gemGeom = (const GEMGeometry*) m_geomRecord->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 879 of file FWTGeoRecoGeometryESProducer.cc.

References EnergyCorrector::c, PVValHelper::dz, 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(), hiPixelPairStep_cff::points, 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 957 of file FWTGeoRecoGeometryESProducer.cc.

References EnergyCorrector::c, PVValHelper::dz, 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(), hiPixelPairStep_cff::points, 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 1105 of file FWTGeoRecoGeometryESProducer.cc.

References EnergyCorrector::c, PVValHelper::dz, 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(), hiPixelPairStep_cff::points, 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 1037 of file FWTGeoRecoGeometryESProducer.cc.

References EnergyCorrector::c, PVValHelper::dz, 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(), hiPixelPairStep_cff::points, 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 ( void )

private

Definition at line 834 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), ME0Geometry::etaPartitions(), cppFunctionSkipper::exception, GetDaughter(), GetTopHolder(), kMuonME0, m_geomRecord, DetId::Muon, dataset::name, 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 = (const ME0Geometry*) m_geomRecord->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 ( void )

private

Definition at line 454 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), TrackerGeometry::detsPXB(), GetDaughter(), GetTopHolder(), kSiPixel, PVValHelper::ladder, m_trackerGeom, m_trackerTopology, dataset::name, 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 ( void )

private

Definition at line 476 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), TrackerGeometry::detsPXF(), GetDaughter(), GetTopHolder(), kSiPixel, m_trackerGeom, m_trackerTopology, dataset::name, 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 ( void )

private

Definition at line 800 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), end, GeomDet::geographicalId(), GetDaughter(), GetTopHolder(), kMuonRPC, m_geomRecord, DetId::Muon, dataset::name, 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 = (const RPCGeometry*) m_geomRecord->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 ( void )

private

Definition at line 567 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), TrackerGeometry::detsTEC(), GetDaughter(), GetTopHolder(), kSiStrip, m_trackerGeom, m_trackerTopology, dataset::name, 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 ( void )

private

Definition at line 501 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), TrackerGeometry::detsTIB(), GetDaughter(), GetTopHolder(), kSiStrip, m_trackerGeom, m_trackerTopology, dataset::name, 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 ( void )

private

Definition at line 523 of file FWTGeoRecoGeometryESProducer.cc.

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

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 ( void )

private

Definition at line 545 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), TrackerGeometry::detsTOB(), GetDaughter(), GetTopHolder(), kSiStrip, m_trackerGeom, m_trackerTopology, dataset::name, 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 350 of file FWTGeoRecoGeometryESProducer.cc.

References b, Surface::bounds(), Bounds::length(), m_nameToShape, dataset::name, nullptr, TrapezoidalPlaneBounds::parameters(), alignCSCRings::s, AlCaHLTBitMon_QueryRunRegistry::string, GeomDet::surface(), Bounds::thickness(), ApeEstimator_cff::width, 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 427 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 116 of file FWTGeoRecoGeometryESProducer.cc.

References GetMedium(), and gen::n.

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 133 of file FWTGeoRecoGeometryESProducer.cc.

References GetMedium(), and gen::n.

 {
    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 173 of file FWTGeoRecoGeometryESProducer.cc.

References create_public_lumi_plots::color, kCaloTower, kECal, kHCal, kHGCE, kHGCH, kMuonCSC, kMuonDT, kMuonGEM, kMuonME0, kMuonRPC, kSiPixel, kSiStrip, m_dummyMedium, m_recoMedium, dataset::name, and AlCaHLTBitMon_QueryRunRegistry::string.

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 151 of file FWTGeoRecoGeometryESProducer.cc.

References GetDaughter(), and tablePrinter::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::shared_ptr< FWTGeoRecoGeometry > FWTGeoRecoGeometryESProducer::produce ( const FWTGeoRecoGeometryRecord & record )

Definition at line 262 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 >::getRecord(), m_calo, m_caloGeom, m_dummyMedium, m_fwGeometry, m_geomRecord, m_muon, m_tracker, m_trackerGeom, m_trackerTopology, edm::ESHandle< T >::product(), GlobalTrackingGeometry::slaveGeometry(), and DetId::Tracker.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

 {
    using namespace edm;
 
    m_fwGeometry = std::make_shared<FWTGeoRecoGeometry>();
   
    if( m_calo ) {
      edm::ESHandle<CaloGeometry> caloH;
      record.getRecord<CaloGeometryRecord>().get(caloH);
      m_caloGeom = caloH.product();
    }
 
    TGeoManager* geom = new TGeoManager( "cmsGeo", "CMS Detector" );
    if( nullptr == gGeoIdentity )
    {
       gGeoIdentity = new TGeoIdentity( "Identity" );
    }
 
    m_fwGeometry->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::shared_ptr<FWTGeoRecoGeometry>();
    }
    geom->SetTopVolume( top );
    // ROOT chokes unless colors are assigned
    top->SetVisibility( kFALSE );
    top->SetLineColor( kBlue );
 
    if( m_tracker || m_muon )
    {
      record.getRecord<GlobalTrackingGeometryRecord>().get( m_geomRecord );
    }
    
    if( m_tracker )
    {
      DetId detId( DetId::Tracker, 0 );
      m_trackerGeom = (const TrackerGeometry*) m_geomRecord->slaveGeometry( detId );
      
      edm::ESHandle<TrackerTopology> trackerTopologyHandle;
      record.getRecord<TrackerTopologyRcd>().get( trackerTopologyHandle );
      m_trackerTopology = trackerTopologyHandle.product();
      
      addPixelBarrelGeometry();
      addPixelForwardGeometry();
 
      addTIBGeometry();
      addTIDGeometry();
      addTOBGeometry();
      addTECGeometry();
    }
    
    if( m_muon )
    {
      addDTGeometry();
      addCSCGeometry();
      addRPCGeometry();
      addME0Geometry();
      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 m_fwGeometry;
 }