#include <FWTGeoRecoGeometryESProducer.h>

Inheritance diagram for FWTGeoRecoGeometryESProducer:

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

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

virtual	~FWTGeoRecoGeometryESProducer (void)

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

virtual	~ESProducer ()

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

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

virtual	~ESProxyFactoryProducer ()

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

Private Member Functions
void	addCSCGeometry ()

void	addDTGeometry ()

void	addEcalCaloGeometry ()

void	addGEMGeometry ()

void	addPixelBarrelGeometry ()

void	addPixelForwardGeometry ()

void	addRPCGeometry ()

void	addTECGeometry ()

void	addTIBGeometry ()

void	addTIDGeometry ()

void	addTOBGeometry ()

TGeoManager *	createManager (int level)

TGeoMaterial *	createMaterial (const std::string &name)

TGeoShape *	createShape (const GeomDet *det)

TGeoVolume *	createVolume (const std::string &name, const GeomDet *det, const std::string &matname="Air")

	FWTGeoRecoGeometryESProducer (const FWTGeoRecoGeometryESProducer &)

TGeoVolume *	getCalloCellVolume (const CaloCellGeometry *cell)

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

TGeoVolume *	getIdealObliquePrismVolume (const CaloCellGeometry *cell)

TGeoVolume *	getIdealZPrismVolume (const CaloCellGeometry *cell)

TGeoVolume *	getTruncatedPyramidVolume (const CaloCellGeometry *cell)

const FWTGeoRecoGeometryESProducer &	operator= (const FWTGeoRecoGeometryESProducer &)

Private Attributes
edm::ESHandle< CaloGeometry >	m_caloGeom

TGeoMedium *	m_dummyMedium

boost::shared_ptr < FWTGeoRecoGeometry >	m_fwGeometry

edm::ESHandle < GlobalTrackingGeometry >	m_geomRecord

std::map< std::string, TGeoMaterial * >	m_nameToMaterial

std::map< std::string, TGeoMedium * >	m_nameToMedium

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

std::map< TGeoShape , TGeoVolume >	m_shapeToVolume

const TrackerGeometry *	m_trackerGeom

Additional Inherited Members
Public Types inherited from edm::eventsetup::DataProxyProvider
typedef std::vector< std::pair < DataKey, boost::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
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::auto_ptr< TFactory > iFactory, const std::string &iLabel=std::string())

virtual void	registerFactoryWithKey (const eventsetup::EventSetupRecordKey &iRecord, std::auto_ptr< eventsetup::ProxyFactoryBase > &iFactory, const std::string &iLabel=std::string())

virtual void	registerProxies (const eventsetup::EventSetupRecordKey &iRecord, KeyedProxies &aProxyList)
	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 28 of file FWTGeoRecoGeometryESProducer.h.

Constructor & Destructor Documentation

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

Definition at line 51 of file FWTGeoRecoGeometryESProducer.cc.

References edm::ESProducer::setWhatProduced().

                                                                                     :
 m_dummyMedium(0)
 {
   setWhatProduced( this );
 }

FWTGeoRecoGeometryESProducer::~FWTGeoRecoGeometryESProducer ( void )

virtual

Definition at line 57 of file FWTGeoRecoGeometryESProducer.cc.

58 {}

FWTGeoRecoGeometryESProducer::FWTGeoRecoGeometryESProducer ( const FWTGeoRecoGeometryESProducer & )

private

Member Function Documentation

void FWTGeoRecoGeometryESProducer::addCSCGeometry ( void )

private

Definition at line 608 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), CSCDetId, GetDaughter(), m_fwGeometry, m_geomRecord, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    if(! m_geomRecord->slaveGeometry( CSCDetId()))
       throw cms::Exception( "FatalError" ) << "Cannnot find CSCGeometry\n";
 
    
    TGeoVolume *assembly = new TGeoVolumeAssembly("CSC");
 
    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 = 0;
 
       if( auto chamber = dynamic_cast<const CSCChamber*>(*it))
          child = createVolume( name, chamber );
       else if( auto * layer = dynamic_cast<const CSCLayer*>(*it))
          child = createVolume( name, layer );
 
 
 
       if (child) {
          TGeoVolume* holder  = GetDaughter(assembly, "Endcap", detId.endcap());
          holder = GetDaughter(holder, "Station", detId.station());
          holder = GetDaughter(holder, "Ring", detId.ring());
          holder = GetDaughter(holder, "Chamber", detId.chamber());
       
          child->SetLineColor( kBlue );
          //   holder->AddNode(child, 1,  createPlacement( *it ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement(*it));
       }
    }
 
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

void FWTGeoRecoGeometryESProducer::addDTGeometry ( void )

private

Definition at line 514 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, DTChamberId, DTLayerId, DTSuperLayerId, end, GetDaughter(), m_fwGeometry, m_geomRecord, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    //
    // DT chambers geometry
    //
    TGeoVolume *assembly = new TGeoVolumeAssembly( "DT");
 
    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 );
          TGeoVolume* holder  = GetDaughter(assembly, "Wheel", detid.wheel());
          holder = GetDaughter(holder, "Station", detid.station());
          holder = GetDaughter(holder, "Sector", detid.sector());
    
          //   holder->AddNode( child, 1, createPlacement( chamber ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement( chamber));
          child->SetLineColor( kRed );
       }
    }
 
    // Fill in DT super layer parameters
    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))
       {
          //         DetId detidx = superlayer->geographicalId();
          DTSuperLayerId detid( DetId(superlayer->geographicalId()));
          std::stringstream s;
          s << detid;
          std::string name = s.str();
       
          TGeoVolume* child = createVolume( name, superlayer );
 
          TGeoVolume* holder  = GetDaughter(assembly, "Wheel", detid.wheel());
          holder = GetDaughter(holder, "Station", detid.station());
          holder = GetDaughter(holder, "Sector", detid.sector());
          holder = GetDaughter(holder, "SuperLayer", detid.superlayer());
          // holder = GetDaughter(holder, "Layer", detid.layer());
 
          //holder->AddNode( child, 1, createPlacement( superlayer ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement( superlayer));
 
          child->SetLineColor( kBlue );
       }
    }
    
    // Fill in DT layer parameters
    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 );
 
          TGeoVolume* holder  = GetDaughter(assembly, "Wheel", detid.wheel());
          holder = GetDaughter(holder, "Station", detid.station());
          holder = GetDaughter(holder, "Sector", detid.sector());
          holder = GetDaughter(holder, "SuperLayer", detid.superlayer());
          holder = GetDaughter(holder, "Layer", detid.layer());
 
          //         holder->AddNode( child, 1, createPlacement( layer ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement( layer));
          child->SetLineColor( kBlue );
       }
    } 
 
 
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

void FWTGeoRecoGeometryESProducer::addEcalCaloGeometry ( void )

private

Definition at line 901 of file FWTGeoRecoGeometryESProducer.cc.

References cond::rpcobgas::detid, DetId::Ecal, EcalBarrel, EcalEndcap, end, getCalloCellVolume(), GetDaughter(), EBDetId::ism(), m_caloGeom, m_fwGeometry, EEDetId::sc(), and EEDetId::zside().

 {
    if (1)
    {
       TGeoVolume *assembly = new TGeoVolumeAssembly("EcalBarrel");
       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* cell = m_caloGeom->getGeometry( *it );
          TGeoVolume* holder  = GetDaughter(assembly, "ism", detid.ism());
          holder->AddNode( getCalloCellVolume(cell), 1, createCaloPlacement(cell));
 
       }
       m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
    }
 
 
    if (1) {
       TGeoVolume *assembly = new TGeoVolumeAssembly("EcalEndcap");
       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* cell = m_caloGeom->getGeometry( *it );
 
          TGeoVolume* holder  = GetDaughter(assembly, "zside", detid.zside());
          holder  = GetDaughter(holder, "sc", detid.sc());
          holder->AddNode( getCalloCellVolume(cell), 1, createCaloPlacement(cell));
       }
       m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
    }
 }

void FWTGeoRecoGeometryESProducer::addGEMGeometry ( void )

private

Definition at line 652 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, end, GEMGeometry::etaPartitions(), MuonSubdetId::GEM, GeomDet::geographicalId(), GetDaughter(), m_fwGeometry, m_geomRecord, DetId::Muon, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 { 
    TGeoVolume *assembly = new TGeoVolumeAssembly("GEM");
 
    DetId detId( DetId::Muon, MuonSubdetId::GEM );
    const GEMGeometry* gemGeom = (const GEMGeometry*) m_geomRecord->slaveGeometry( detId );
    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 );
 
          TGeoVolume* holder  = GetDaughter(assembly, "ROLL Region", detid.region());
          holder = GetDaughter(holder, "Ring", detid.ring());
          holder = GetDaughter(holder, "Station", detid.station()); 
          holder = GetDaughter(holder, "Layer", detid.layer()); 
          holder = GetDaughter(holder, "Chamber", detid.chamber()); 
 
          AddLeafNode(holder, child, name.c_str(),  createPlacement(*it));
 
          child->SetLineColor( kYellow );     
       }
    }
 
    printf("ADD GEM!!!\n");
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

void FWTGeoRecoGeometryESProducer::addPixelBarrelGeometry ( void )

private

Definition at line 330 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, TrackerGeometry::detsPXB(), end, GetDaughter(), PXBDetId::layer(), m_fwGeometry, m_trackerGeom, PXBDetId::module(), DetId::rawId(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    TGeoVolume *assembly = new TGeoVolumeAssembly("PXB");
 
    for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsPXB().begin(),
            end = m_trackerGeom->detsPXB().end();
         it != end; ++it)
    {
        DetId detid = ( *it )->geographicalId();
        unsigned int rawid = detid.rawId();
 
        PXBDetId xx(rawid);
        std::string name = Form("PXB Ly:%d, Md:%d Ld:%d ", xx.layer(), xx.module(), xx.layer());
        TGeoVolume* child = createVolume( name, *it );
        child->SetLineColor( kGreen );
 
        TGeoVolume* holder  = GetDaughter(assembly, "Layer", xx.layer());
        holder = GetDaughter(holder, "Module", xx.module());
                                  
        AddLeafNode(holder, child, name.c_str(), createPlacement( *it ));
    }
   
 
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

void FWTGeoRecoGeometryESProducer::addPixelForwardGeometry ( void )

private

Definition at line 359 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, TrackerGeometry::detsPXF(), end, GetDaughter(), m_fwGeometry, m_trackerGeom, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    TGeoVolume *assembly = new TGeoVolumeAssembly("PXF");
    for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsPXF().begin(),
            end = m_trackerGeom->detsPXF().end();
         it != end; ++it )
    {
       PXFDetId detid = ( *it )->geographicalId();
 
       std::stringstream s;
       s << detid;
       std::string name = s.str();
 
       TGeoVolume* child = createVolume( name, *it );
 
 
       TGeoVolume* holder  = GetDaughter(assembly, "Side", detid.side());
       holder = GetDaughter(holder, "Disk", detid.disk());
       holder = GetDaughter(holder, "Blade", detid.blade());
       holder = GetDaughter(holder, "Panel", detid.panel());
    
       // holder->AddNode( child, 1, createPlacement( *it ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement( *it ));
       child->SetLineColor( kGreen );
 
    }
   
 
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

void FWTGeoRecoGeometryESProducer::addRPCGeometry ( void )

private

Definition at line 692 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, end, GeomDet::geographicalId(), GetDaughter(), m_fwGeometry, m_geomRecord, DetId::Muon, RPCGeometry::rolls(), MuonSubdetId::RPC, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    TGeoVolume *assembly = new TGeoVolumeAssembly("RPC");
 
    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 );
 
          TGeoVolume* holder  = GetDaughter(assembly, "ROLL Region", detid.region());
          holder = GetDaughter(holder, "Ring", detid.ring());
          holder = GetDaughter(holder, "Station", detid.station()); 
          holder = GetDaughter(holder, "Sector", detid.sector()); 
          holder = GetDaughter(holder, "Layer", detid.layer()); 
          holder = GetDaughter(holder, "Subsector", detid.subsector()); 
 
          AddLeafNode(holder, child, name.c_str(),  createPlacement(*it));
 
          child->SetLineColor( kYellow );     
       }
    }
    /*
   //AMT why no chamber ???
   for( auto it = rpcGeom->chambers().begin(),
            end = rpcGeom->chambers().end(); 
         it != end; ++it )
    {
       RPCChamber const* chamber = (*it);
       if( chamber )
       {
          RPCDetId detid = chamber->geographicalId();
          std::stringstream s;
          s << detid;
          std::string name = s.str();
       
          TGeoVolume* child = createVolume( name, chamber );
 
          TGeoVolume* holder  = GetDaughter(assembly, "CHANBER Region", detid.region());
          holder = GetDaughter(holder, "Ring", detid.ring());
          holder = GetDaughter(holder, "Station", detid.station()); 
          holder = GetDaughter(holder, "Sector", detid.sector()); 
          holder = GetDaughter(holder, "Layer", detid.layer()); 
          holder = GetDaughter(holder, "Subsector", detid.subsector()); 
 
          holder->AddNode( child, 1, createPlacement( chamber ));
          child->SetLineColor( kYellow );     
       }
       else printf("NO CHAMBER \n");
    }
    */
 
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

void FWTGeoRecoGeometryESProducer::addTECGeometry ( void )

private

Definition at line 484 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, TrackerGeometry::detsTEC(), end, GetDaughter(), m_fwGeometry, m_trackerGeom, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    TGeoVolume *assembly = new TGeoVolumeAssembly("TEC");
 
    for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTEC().begin(),
            end = m_trackerGeom->detsTEC().end();
         it != end; ++it )
    {
       TECDetId detid = ( *it )->geographicalId();
 
       std::stringstream s;
       s << detid;
       std::string name = s.str();
 
       TGeoVolume* child = createVolume( name, *it );
 
       TGeoVolume* holder  = GetDaughter(assembly, "Order", detid.order());
       holder = GetDaughter(holder, "Ring", detid.ring());
       holder = GetDaughter(holder, "Module", detid.module());
       // holder->AddNode( child, 1, createPlacement( *it ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement( *it ));
       child->SetLineColor( kGreen );
    }
 
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

void FWTGeoRecoGeometryESProducer::addTIBGeometry ( void )

private

Definition at line 394 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, TrackerGeometry::detsTIB(), end, GetDaughter(), m_fwGeometry, m_trackerGeom, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    TGeoVolume *assembly = new TGeoVolumeAssembly("TIB");
    for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTIB().begin(),
            end = m_trackerGeom->detsTIB().end();
         it != end; ++it )
    {
       TIBDetId detid(( *it )->geographicalId());
 
       std::stringstream s;
       s << detid;
       std::string name = s.str();
 
       TGeoVolume* child = createVolume( name, *it );
 
       TGeoVolume* holder  = GetDaughter(assembly, "Module", detid.module());
       holder = GetDaughter(holder, "Order", detid.order());
       holder = GetDaughter(holder, "Side", detid.side());
       // holder->AddNode( child, 1, createPlacement( *it ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement( *it ));
       child->SetLineColor( kGreen );
    }
   
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

void FWTGeoRecoGeometryESProducer::addTIDGeometry ( void )

private

Definition at line 424 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, TrackerGeometry::detsTID(), end, GetDaughter(), m_fwGeometry, m_trackerGeom, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    TGeoVolume *assembly = new TGeoVolumeAssembly("TID");
 
    for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTID().begin(),
            end = m_trackerGeom->detsTID().end();
         it != end; ++it)
    {
       TIDDetId detid = ( *it )->geographicalId();
 
       std::stringstream s;
       s << detid;
       std::string name = s.str();
 
       TGeoVolume* child = createVolume( name, *it );
       TGeoVolume* holder  = GetDaughter(assembly, "Side", detid.side());
       holder = GetDaughter(holder, "Wheel", detid.wheel());
       holder = GetDaughter(holder, "Ring", detid.ring());
       //  holder->AddNode( child, 1, createPlacement( *it ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement( *it ));
    
       child->SetLineColor( kGreen );
    }
 
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

void FWTGeoRecoGeometryESProducer::addTOBGeometry ( void )

private

Definition at line 454 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, TrackerGeometry::detsTOB(), end, GetDaughter(), m_fwGeometry, m_trackerGeom, alignCSCRings::s, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by produce().

 {
    TGeoVolume *assembly = new TGeoVolumeAssembly("TOB");
 
    for( TrackerGeometry::DetContainer::const_iterator it = m_trackerGeom->detsTOB().begin(),
            end = m_trackerGeom->detsTOB().end();
         it != end; ++it )
    {
       TOBDetId detid(( *it )->geographicalId());
 
       std::stringstream s;
       s << detid;
       std::string name = s.str();
 
       TGeoVolume* child = createVolume( name, *it );
       TGeoVolume* holder  = GetDaughter(assembly, "Rod", detid.rodNumber());
       holder = GetDaughter(holder, "Side", detid.side());
       holder = GetDaughter(holder, "Module", detid.moduleNumber());
       //holder->AddNode( child, 1, createPlacement( *it ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement( *it ));
    
       child->SetLineColor( kGreen );
    }
 
    m_fwGeometry->manager()->GetTopVolume()->AddNode( assembly, 1 );
 }

TGeoManager* FWTGeoRecoGeometryESProducer::createManager ( int level )

private

TGeoMaterial * FWTGeoRecoGeometryESProducer::createMaterial ( const std::string & name )

private

Create TGeo material based on its name

Definition at line 303 of file FWTGeoRecoGeometryESProducer.cc.

References m_nameToMaterial, and mergeVDriftHistosByStation::name.

Referenced by createVolume().

 {
   TGeoMaterial *material = m_nameToMaterial[name];
 
   if( material == 0 )
   {
     // FIXME: Do we need to set real parameters of the material?
     material = new TGeoMaterial( name.c_str(),
                                  0, 0, 0 );
     m_nameToMaterial[name] = material;
   }
 
   return material;
 }

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

private

Create TGeo shape for GeomDet

Definition at line 203 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by createVolume().

 {
   TGeoShape* shape = 0;
 
   // 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( 0 == 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) != 0 )
   {
     // 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( 0 == 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,
		const std::string &	material = `"Air"`
	)

private

Create TGeo volume for GeomDet

Definition at line 280 of file FWTGeoRecoGeometryESProducer.cc.

References createMaterial(), createShape(), m_nameToMedium, and m_shapeToVolume.

Referenced by addCSCGeometry(), addDTGeometry(), addGEMGeometry(), 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;
    
 
    TGeoMedium* medium = m_nameToMedium[material];
    if( 0 == medium )
    {
       medium = new TGeoMedium( material.c_str(), 0, createMaterial( material ));
       m_nameToMedium[material] = medium;
    }
    TGeoVolume* volume = new TGeoVolume( name.c_str(),solid, medium);
 
    m_shapeToVolume[solid] = volume;
 
    return volume;
 }

TGeoVolume * FWTGeoRecoGeometryESProducer::getCalloCellVolume ( const CaloCellGeometry * cell )

private

Definition at line 882 of file FWTGeoRecoGeometryESProducer.cc.

References getIdealObliquePrismVolume(), getIdealZPrismVolume(), and getTruncatedPyramidVolume().

Referenced by addEcalCaloGeometry().

 {
    if (dynamic_cast<const TruncatedPyramid*> (cell)) 
       return getTruncatedPyramidVolume(cell);
 
    if (dynamic_cast<const IdealZPrism*> (cell))
       return getIdealZPrismVolume(cell);
 
    if (dynamic_cast<const IdealObliquePrism*> (cell)) 
       return getIdealObliquePrismVolume(cell);
 
    return 0;
 }

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

private

Definition at line 118 of file FWTGeoRecoGeometryESProducer.cc.

References m_dummyMedium.

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

 {
    TGeoVolume* res = 0;
    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(m_dummyMedium);
       mother->AddNode(res, 1);
    }
 
    return res;
 }

TGeoVolume * FWTGeoRecoGeometryESProducer::getIdealObliquePrismVolume ( const CaloCellGeometry * cell )

private

Definition at line 819 of file FWTGeoRecoGeometryESProducer.cc.

Referenced by getCalloCellVolume().

 {
    TGeoVolume* volume = 0;
 
    // AMT code unfinished ...
    /*
    CaloVolMap::iterator volIt =  m_caloShapeMap.find(cell->param());
 
    if  (volIt == m_caloShapeMap.end()) {
       printf("FIREWORKS NEW SHAPE BEGIN >>>>>> \n");
      
       double points[24];
       IdealObliquePrism::Pt3DVec co(8);
       IdealObliquePrism::Pt3D ref;
       IdealObliquePrism::localCorners( co, cell->param(), ref);
       
       static const int arr[] = {2, 3, 0, 1, 6, 7, 4,5};
         //static const int arr[] = { 0, 1, 2, 3, 4, 5, 6, 7};
 
       unsigned int idx( 0 ); 
       for( int c = 0; c < 8; ++c, ++idx )
       {
          points[idx*3]     = co[arr[c]].x();
          points[idx*3 + 1] = co[arr[c]].y();
          points[idx*3 + 2] = co[arr[c]].z();
          printf("fw oblique lc [%d] = %.4f %.4f %.4f\n", c, points[idx*3], points[idx*3+1], points[idx*3+2]);
       }  
 
 
        if (cell->param()[2] < 0) {
          printf("negatice height %f\n", cell->param()[2]);
       }
 
       TGeoShape* solid = new TGeoArb8(TMath::Abs(cell->param()[2]), points); 
       volume = new TGeoVolume( "ObliquePrism" ,solid, m_dummyMedium);
       
        
 
       //volume =  m_fwGeometry->manager()->MakeBox( "Z prsim", m_dummyMedium, dz, dz, dz );
 
 
       volume->SetLineColor(kGray);
       m_caloShapeMap[cell->param()]  = volume;
       CaloCellGeometry::CornersVec const & cv = cell->getCorners();
       printf(" cell ------------------- \n");
       for( int c = 0; c < 8; ++c)
          printf("global cell corners [%d] = %.4f %.4f %.4f\n", c, cv[c].x(),cv[c].y(),cv[c].z() );
 
       printf("FIREWORKS NEW SHAPE END >>>>>> \n");
    }
    else {
       volume = volIt->second;
    }
 
    */
    return volume;
 }

TGeoVolume * FWTGeoRecoGeometryESProducer::getIdealZPrismVolume ( const CaloCellGeometry * cell )

private

Definition at line 810 of file FWTGeoRecoGeometryESProducer.cc.

Referenced by getCalloCellVolume().

 {
    TGeoVolume* volume = 0;
 
    return volume;
 }

TGeoVolume * FWTGeoRecoGeometryESProducer::getTruncatedPyramidVolume ( const CaloCellGeometry * cell )

private

Definition at line 770 of file FWTGeoRecoGeometryESProducer.cc.

References trackerHits::c, cmsHarvester::index, m_dummyMedium, and CaloCellGeometry::param().

Referenced by getCalloCellVolume().

 {
    TGeoVolume* volume = 0;
 
    CaloVolMap::iterator volIt =  m_caloShapeMap.find(cell->param());
    if  (volIt == m_caloShapeMap.end()) {
       double points[24];
       const HepGeom::Transform3D idtr;
       std::vector<float> vpar;
       for (int ip =0; ip < 11; ++ip)
          vpar.push_back(cell->param()[ip]);
 
       std::vector<GlobalPoint> co(8);
       TruncatedPyramid::createCorners(vpar, idtr, co);
 
       unsigned int index( 0 ); 
       static const int arr[] = {0, 3, 2, 1, 4, 7, 6, 5};
       for( int c = 0; c < 8; ++c )
       {
          points[index] = co[arr[c]].x();
          points[++index] = co[arr[c]].y();
          points[++index] = co[arr[c]].z();
       }  
 
       TGeoShape* solid = new TGeoArb8(cell->param()[0], points); 
 
 
       volume = new TGeoVolume( "TruncatedPyramid" ,solid, m_dummyMedium);
       volume->SetLineColor(kGray);
       m_caloShapeMap[cell->param()]  = volume;
    }
    else {
       volume = volIt->second;
    }
 
    return volume;
 }

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

private

boost::shared_ptr< FWTGeoRecoGeometry > FWTGeoRecoGeometryESProducer::produce ( const FWTGeoRecoGeometryRecord & record )

Definition at line 136 of file FWTGeoRecoGeometryESProducer.cc.

References addCSCGeometry(), addDTGeometry(), addGEMGeometry(), addPixelBarrelGeometry(), addPixelForwardGeometry(), addRPCGeometry(), addTECGeometry(), addTIBGeometry(), addTIDGeometry(), addTOBGeometry(), cppFunctionSkipper::exception, relativeConstraints::geom, edm::eventsetup::DependentRecordImplementation< RecordT, ListT >::getRecord(), m_caloGeom, m_dummyMedium, m_fwGeometry, m_geomRecord, m_nameToMaterial, m_nameToMedium, m_nameToShape, m_shapeToVolume, m_trackerGeom, fff_deletion::top, DetId::Tracker, and cms::Exception::what().

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

 {
    using namespace edm;
 
   m_fwGeometry = boost::shared_ptr<FWTGeoRecoGeometry>( new FWTGeoRecoGeometry );
 
    record.getRecord<GlobalTrackingGeometryRecord>().get( m_geomRecord );
   
    DetId detId( DetId::Tracker, 0 );
    m_trackerGeom = (const TrackerGeometry*) m_geomRecord->slaveGeometry( detId );
   
    record.getRecord<CaloGeometryRecord>().get( m_caloGeom );
 
    TGeoManager* geom = new TGeoManager( "cmsGeo", "CMS Detector" );
    if( 0 == 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);
    // so is default medium
    TGeoVolume *top = geom->MakeBox( "CMS", m_dummyMedium, 270., 270., 120. );
   
    if( 0 == top )
    {
      return boost::shared_ptr<FWTGeoRecoGeometry>();
    }
    geom->SetTopVolume( top );
    // ROOT chokes unless colors are assigned
    top->SetVisibility( kFALSE );
    top->SetLineColor( kBlue );
    
    addPixelBarrelGeometry();
    addPixelForwardGeometry();
    addTIBGeometry();
    addTIDGeometry();
    addTOBGeometry();
    addTECGeometry();
    addDTGeometry();
    addCSCGeometry();
    addRPCGeometry();
    try {
       addGEMGeometry();
    }
    catch ( cms::Exception& exception ) {
    edm::LogWarning("FWRecoGeometryProducerException")
      << "addGEMGeometry() Exception caught while building GEM geometry: " << exception.what()
      << std::endl; 
    }
 
    //  addEcalCaloGeometry();
    geom->CloseGeometry();
 
    m_nameToShape.clear();
    m_shapeToVolume.clear();
    m_nameToMaterial.clear();
    m_nameToMedium.clear();
 
   return m_fwGeometry;
 }