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

Private Member Functions
void	addCSCGeometry ()

void	addDTGeometry ()

void	addEcalCaloGeometry ()

void	addGEMGeometry ()

void	addHcalCaloGeometryBarrel ()

void	addHcalCaloGeometryEndcap ()

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
edm::ESHandle< CaloGeometry >	m_caloGeom

TGeoMedium *	m_dummyMedium

boost::shared_ptr < FWTGeoRecoGeometry >	m_fwGeometry

edm::ESHandle < GlobalTrackingGeometry >	m_geomRecord

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

std::map< ERecoDet, TGeoMedium * >	m_recoMedium

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.

Member Enumeration Documentation

enum FWTGeoRecoGeometryESProducer::ERecoDet

private

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

Definition at line 30 of file FWTGeoRecoGeometryESProducer.h.

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

Constructor & Destructor Documentation

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

Definition at line 63 of file FWTGeoRecoGeometryESProducer.cc.

References edm::ESProducer::setWhatProduced().

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

FWTGeoRecoGeometryESProducer::~FWTGeoRecoGeometryESProducer ( void )

virtual

Definition at line 69 of file FWTGeoRecoGeometryESProducer.cc.

70 {}

FWTGeoRecoGeometryESProducer::FWTGeoRecoGeometryESProducer ( const FWTGeoRecoGeometryESProducer & )

private

Member Function Documentation

void FWTGeoRecoGeometryESProducer::addCSCGeometry ( void )

private

Definition at line 663 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), CSCDetId, GetDaughter(), GetTopHolder(), kMuonCSC, kMuonRPC, 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* 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 = 0;
 
       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 574 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, DTChamberId, DTLayerId, DTSuperLayerId, end, GetDaughter(), GetTopHolder(), kMuonDT, kMuonRPC, m_geomRecord, alignCSCRings::s, 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 1031 of file FWTGeoRecoGeometryESProducer.cc.

References cond::rpcobgas::detid, DetId::Ecal, EcalBarrel, EcalEndcap, end, CaloCellGeometry::getCorners(), GetDaughter(), getEcalTrans(), GetMedium(), GetTopHolder(), 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 TruncatedPyramid* cell = dynamic_cast<const TruncatedPyramid*> ( m_caloGeom->getGeometry( *it ));
          if (!cell) { printf("ecalBarrel cell not a TruncatedPyramid !!\n"); return; }
 
          TGeoVolume* volume = 0;
          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 TruncatedPyramid* cell = dynamic_cast<const TruncatedPyramid*> (m_caloGeom->getGeometry( *it ));
          if (!cell) { printf("ecalEndcap cell not a TruncatedPyramid !!\n"); continue;}
 
          TGeoVolume* volume = 0;
          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 706 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by produce().

 { 
   
    DetId detId( DetId::Muon, MuonSubdetId::GEM );
    const GEMGeometry* gemGeom = (const GEMGeometry*) m_geomRecord->slaveGeometry( detId );
 
    TGeoVolume* tv =  GetTopHolder("Muon", kMuonRPC);
    TGeoVolume *assembly = GetDaughter(tv, "GEM", 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));
       }
    }
 }

void FWTGeoRecoGeometryESProducer::addHcalCaloGeometryBarrel ( void )

private

Definition at line 797 of file FWTGeoRecoGeometryESProducer.cc.

References trackerHits::c, cond::rpcobgas::detid, end, CaloCellGeometry::etaPos(), CaloCellGeometry::getCorners(), GetDaughter(), GetMedium(), GetTopHolder(), DetId::Hcal, HcalBarrel, HcalDetId::ieta(), kHCal, IdealObliquePrism::localCorners(), m_caloGeom, CaloCellGeometry::param(), CaloCellGeometry::phiPos(), makeMuonMisalignmentScenario::rot, x, detailsBasic3DVector::y, detailsBasic3DVector::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());
 
       const CaloCellGeometry* cellb= m_caloGeom->getGeometry(*it);
 
       const IdealObliquePrism* cell = dynamic_cast<const IdealObliquePrism*> (cellb);
    
       if (!cell) { printf ("HB not olique !!!\n"); continue; }
   
       TGeoVolume* volume = 0;
       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 881 of file FWTGeoRecoGeometryESProducer.cc.

References trackerHits::c, cond::rpcobgas::detid, end, GetDaughter(), GetMedium(), GetTopHolder(), DetId::Hcal, HcalEndcap, HcalDetId::ieta(), kHCal, IdealObliquePrism::localCorners(), m_caloGeom, 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 IdealObliquePrism* cell = dynamic_cast<const IdealObliquePrism*> ( m_caloGeom->getGeometry(*it));
    
       if (!cell) { printf ("EC not olique \n"); continue; }
 
       TGeoVolume* volume = 0;
       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::addPixelBarrelGeometry ( void )

private

Definition at line 403 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by produce().

 {
    TGeoVolume* tv =  GetTopHolder("SiPixel", kSiPixel);
    TGeoVolume *assembly = GetDaughter(tv, "PXB", kSiPixel);
 
    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, kSiPixel );
 
        TGeoVolume* holder  = GetDaughter(assembly, "Layer", kSiPixel, xx.layer());
        holder = GetDaughter(holder, "Module", kSiPixel, xx.module());
                                  
        AddLeafNode(holder, child, name.c_str(), createPlacement( *it ));
    }
   
 
 }

void FWTGeoRecoGeometryESProducer::addPixelForwardGeometry ( void )

private

Definition at line 431 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by produce().

 {
    TGeoVolume* tv =  GetTopHolder("SiPixel", kSiPixel);
    TGeoVolume *assembly = GetDaughter(tv, "PXF", kSiPixel);
 
    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, kSiPixel );
 
 
       TGeoVolume* holder  = GetDaughter(assembly, "Side", kSiPixel, detid.side());
       holder = GetDaughter(holder, "Disk", kSiPixel, detid.disk());
       holder = GetDaughter(holder, "Blade", kSiPixel, detid.blade());
       holder = GetDaughter(holder, "Panel", kSiPixel, detid.panel());
    
       // holder->AddNode( child, 1, createPlacement( *it ));
       AddLeafNode(holder, child, name.c_str(),  createPlacement( *it ));
 
    }
   
 }

void FWTGeoRecoGeometryESProducer::addRPCGeometry ( void )

private

Definition at line 744 of file FWTGeoRecoGeometryESProducer.cc.

References createVolume(), cond::rpcobgas::detid, end, GeomDet::geographicalId(), GetDaughter(), GetTopHolder(), kMuonRPC, m_geomRecord, DetId::Muon, RPCGeometry::rolls(), MuonSubdetId::RPC, alignCSCRings::s, 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 542 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by produce().

 {
    TGeoVolume* tv =  GetTopHolder( "SiStrip", kSiStrip);
    TGeoVolume *assembly = GetDaughter(tv, "TEC", kSiStrip);
 
    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, kSiStrip );
 
       TGeoVolume* holder  = GetDaughter(assembly, "Order", kSiStrip, detid.order());
       holder = GetDaughter(holder, "Ring", kSiStrip, detid.ring());
       holder = GetDaughter(holder, "Module", kSiStrip, detid.module());
       AddLeafNode(holder, child, name.c_str(),  createPlacement( *it ));
    }
 }

void FWTGeoRecoGeometryESProducer::addTIBGeometry ( void )

private

Definition at line 464 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by produce().

 {
    TGeoVolume* tv =  GetTopHolder( "SiStrip", kSiStrip);
    TGeoVolume *assembly = GetDaughter(tv,"TIB", kSiStrip);
 
    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, kSiStrip );
 
       TGeoVolume* holder  = GetDaughter(assembly, "Module", kSiStrip, detid.module());
       holder = GetDaughter(holder, "Order", kSiStrip, detid.order());
       holder = GetDaughter(holder, "Side", kSiStrip, detid.side());
       AddLeafNode(holder, child, name.c_str(),  createPlacement( *it ));
    }
 }

void FWTGeoRecoGeometryESProducer::addTIDGeometry ( void )

private

Definition at line 491 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by produce().

 {
    TGeoVolume* tv =  GetTopHolder( "SiStrip", kSiStrip);
    TGeoVolume *assembly = GetDaughter( tv, "TID", kSiStrip);
 
    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, kSiStrip );
       TGeoVolume* holder  = GetDaughter(assembly, "Side", kSiStrip, detid.side());
       holder = GetDaughter(holder, "Wheel", kSiStrip, detid.wheel());
       holder = GetDaughter(holder, "Ring", kSiStrip, detid.ring());
       AddLeafNode(holder, child, name.c_str(),  createPlacement( *it ));
    }
 }

void FWTGeoRecoGeometryESProducer::addTOBGeometry ( void )

private

Definition at line 516 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by produce().

 {
    TGeoVolume* tv =  GetTopHolder( "SiStrip", kSiStrip);
    TGeoVolume *assembly = GetDaughter(tv, "TOB",  kSiStrip);
 
    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, kSiStrip );
       TGeoVolume* holder  = GetDaughter(assembly, "Rod", kSiStrip, detid.rodNumber());
       holder = GetDaughter(holder, "Side", kSiStrip, detid.side());
       holder = GetDaughter(holder, "Module", kSiStrip, detid.moduleNumber());
       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 298 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,
		ERecoDet	mid = `kDummy`
	)

private

Create TGeo volume for GeomDet

Definition at line 375 of file FWTGeoRecoGeometryESProducer.cc.

References createShape(), GetMedium(), 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;
    
 
    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 111 of file FWTGeoRecoGeometryESProducer.cc.

References GetMedium().

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

 {
    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(GetMedium(cidx));
       mother->AddNode(res, 1);
    }
 
    return res;
 }

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

private

Definition at line 128 of file FWTGeoRecoGeometryESProducer.cc.

References GetMedium().

 {
    TGeoVolume* res = 0;
    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 155 of file FWTGeoRecoGeometryESProducer.cc.

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

Referenced by addEcalCaloGeometry(), addHcalCaloGeometryBarrel(), addHcalCaloGeometryEndcap(), 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;
    // see TGeoManager::DefaultColors()
    switch (det)
    {
       case kSiPixel:
          name = "SiPixel";
          color = 4;//kGreen-10
          break;
       case kSiStrip:
          name = "SiStrip";
          color = 29; // kOrange+9
          break;
       case kMuonDT:
          name = "MuonDT";
          color = 13; // kBlue -10
          break;
       case kMuonGEM:
          name = "MuonGEM";
          color = 16; //kYellow+1
          break;
       case kMuonRPC:
          name = "MuonRPC";
          color = 10; //kRed-10
          break;
       case kMuonCSC:
          name = "MuonCSC"; // kMagenta
          color = 9;
          break;
       case kECal:
          name = "ECal";
          color = 6;
          break;
       case kHCal:
          name = "HCal";    
          color = 10; // kRed-10
          break;
       case kHGCE:
          name = "HGCEE";
          color = 10;
          break;
       case kHGCH:
          name = "HGCEH";
          color = 6; // kBlue -7
          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 146 of file FWTGeoRecoGeometryESProducer.cc.

References GetDaughter(), and prof2calltree::prefix.

Referenced by addCSCGeometry(), addDTGeometry(), addEcalCaloGeometry(), addGEMGeometry(), addHcalCaloGeometryBarrel(), addHcalCaloGeometryEndcap(), 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

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

Definition at line 226 of file FWTGeoRecoGeometryESProducer.cc.

References addCSCGeometry(), addDTGeometry(), addEcalCaloGeometry(), addGEMGeometry(), addHcalCaloGeometryBarrel(), addHcalCaloGeometryEndcap(), 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_trackerGeom, fff_deleter::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);
 
 
    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();
    
    addHcalCaloGeometryBarrel();
    addHcalCaloGeometryEndcap();
 
    geom->CloseGeometry();
 
    geom->DefaultColors();
    geom->CloseGeometry();
 
    return m_fwGeometry;
 }