#include <HGCalGeometry.h>

Inheritance diagram for HGCalGeometry:

Public Types
enum	{ k_NumberOfParametersPerShape = 12 }

enum	{ k_NumberOfShapes = 50 }

typedef HGCalGeometryRecord	AlignedRecord

typedef CaloCellGeometry::CCGFloat	CCGFloat

typedef std::vector< FlatTrd >	CellVec

typedef std::vector< GlobalPoint >	CornersVec

typedef std::set< DetId >	DetIdSet

typedef PHGCalRcd	PGeometryRecord

typedef CaloCellGeometry::Pt3D	Pt3D

typedef CaloCellGeometry::Pt3DVec	Pt3DVec

Public Types inherited from CaloSubdetectorGeometry
typedef CaloCellGeometry::CCGFloat	CCGFloat

typedef std::vector< std::shared_ptr< const CaloCellGeometry > >	CellSet

typedef std::set< DetId >	DetIdSet

typedef std::vector< CCGFloat >	DimVec

typedef std::vector< unsigned int >	IVec

typedef CaloCellGeometry::ParMgr	ParMgr

typedef CaloCellGeometry::ParVec	ParVec

typedef CaloCellGeometry::ParVecVec	ParVecVec

typedef std::vector< CCGFloat >	TrVec

Public Member Functions
std::string	cellElement () const

virtual void	fillNamedParams (DDFilteredView fv)

DetIdSet	getCells (const GlobalPoint &r, double dR) const override
	Get a list of all cells within a dR of the given cell. More...

DetId	getClosestCell (const GlobalPoint &r) const override

CornersVec	getCorners (const DetId &id) const
	Returns the corner points of this cell's volume. More...

std::shared_ptr< const CaloCellGeometry >	getGeometry (const DetId &id) const override
	Get the cell geometry of a given detector id. Should return false if not found. More...

GlobalPoint	getPosition (const DetId &id) const

void	getSummary (CaloSubdetectorGeometry::TrVec &trVector, CaloSubdetectorGeometry::IVec &iVector, CaloSubdetectorGeometry::DimVec &dimVector, CaloSubdetectorGeometry::IVec &dinsVector) const override

const std::vector< DetId > &	getValidDetIds (DetId::Detector det=DetId::Detector(0), int subdet=0) const override
	Get a list of valid detector ids (for the given subdetector) More...

const std::vector< DetId > &	getValidGeomDetIds (void) const

	HGCalGeometry (const HGCalTopology &topology)

void	initializeParms () override

void	localCorners (Pt3DVec &lc, const CCGFloat *pv, unsigned int i, Pt3D &ref)

void	newCell (const GlobalPoint &f1, const GlobalPoint &f2, const GlobalPoint &f3, const CCGFloat *parm, const DetId &detId) override

bool	present (const DetId &id) const override
	is this detid present in the geometry? More...

void	sortDetIds ()

const HGCalTopology &	topology () const

	~HGCalGeometry () override

Public Member Functions inherited from CaloSubdetectorGeometry
void	allocateCorners (CaloCellGeometry::CornersVec::size_type n)

void	allocatePar (ParVec::size_type n, unsigned int m)

	CaloSubdetectorGeometry ()

	CaloSubdetectorGeometry (const CaloSubdetectorGeometry &)=delete
	avoid copies More...

CaloCellGeometry::CornersMgr *	cornersMgr ()

CCGFloat	deltaEta (const DetId &detId) const

CCGFloat	deltaPhi (const DetId &detId) const

virtual void	fillDefaultNamedParameters () const

virtual CellSet	getCellSet (const GlobalPoint &r, double dR) const

virtual unsigned int	numberOfParametersPerShape () const

virtual unsigned int	numberOfShapes () const

virtual unsigned int	numberOfTransformParms () const

CaloSubdetectorGeometry &	operator= (const CaloSubdetectorGeometry &)=delete

ParMgr *	parMgr ()

const ParMgr *	parMgrConst () const

ParVecVec &	parVecVec ()

const ParVecVec &	parVecVec () const

virtual	~CaloSubdetectorGeometry ()
	The base class DOES assume that it owns the CaloCellGeometry objects. More...

Static Public Member Functions
static std::string	dbString ()

static std::string	producerTag ()

Protected Member Functions
void	addValidID (const DetId &id)

std::shared_ptr< const CaloCellGeometry >	cellGeomPtr (uint32_t index) const override

unsigned int	getClosestCellIndex (const GlobalPoint &r) const

const CaloCellGeometry *	getGeometryRawPtr (uint32_t index) const override

unsigned int	indexFor (const DetId &id) const override

unsigned int	sizeForDenseIndex () const

Protected Member Functions inherited from CaloSubdetectorGeometry
void	addValidID (const DetId &id)

virtual unsigned int	sizeForDenseIndex (const DetId &id) const

Private Member Functions
std::shared_ptr< const CaloCellGeometry >	cellGeomPtr (uint32_t index, const GlobalPoint &p) const

Private Attributes
CellVec	m_cellVec

bool	m_halfType

ForwardSubdetector	m_subdet

const HGCalTopology &	m_topology

std::vector< DetId >	m_validGeomIds

Additional Inherited Members
Static Protected Member Functions inherited from CaloSubdetectorGeometry
static CCGFloat	deltaR (const GlobalPoint &p1, const GlobalPoint &p2)

Protected Attributes inherited from CaloSubdetectorGeometry
ParVecVec	m_parVecVec

std::vector< DetId >	m_validIds

Detailed Description

Definition at line 26 of file HGCalGeometry.h.

Member Typedef Documentation

typedef HGCalGeometryRecord HGCalGeometry::AlignedRecord

Definition at line 39 of file HGCalGeometry.h.

typedef CaloCellGeometry::CCGFloat HGCalGeometry::CCGFloat

Definition at line 32 of file HGCalGeometry.h.

typedef std::vector<FlatTrd> HGCalGeometry::CellVec

Definition at line 30 of file HGCalGeometry.h.

typedef std::vector<GlobalPoint> HGCalGeometry::CornersVec

Definition at line 37 of file HGCalGeometry.h.

typedef std::set<DetId> HGCalGeometry::DetIdSet

Definition at line 36 of file HGCalGeometry.h.

typedef PHGCalRcd HGCalGeometry::PGeometryRecord

Definition at line 40 of file HGCalGeometry.h.

typedef CaloCellGeometry::Pt3D HGCalGeometry::Pt3D

Definition at line 33 of file HGCalGeometry.h.

typedef CaloCellGeometry::Pt3DVec HGCalGeometry::Pt3DVec

Definition at line 34 of file HGCalGeometry.h.

Member Enumeration Documentation

anonymous enum

Enumerator
k_NumberOfParametersPerShape

Definition at line 42 of file HGCalGeometry.h.

42 { k_NumberOfParametersPerShape = 12 } ; // FlatTrd

HGCalGeometry::k_NumberOfParametersPerShape

Definition: HGCalGeometry.h:42

anonymous enum

Enumerator
k_NumberOfShapes

Definition at line 43 of file HGCalGeometry.h.

43 { k_NumberOfShapes = 50 } ;

HGCalGeometry::k_NumberOfShapes

Definition: HGCalGeometry.h:43

Constructor & Destructor Documentation

HGCalGeometry::HGCalGeometry ( const HGCalTopology & topology )

Definition at line 23 of file HGCalGeometry.cc.

References gather_cfg::cout, CaloSubdetectorGeometry::m_validIds, topology(), and HGCalTopology::totalGeomModules().

Referenced by dbString().

   : m_topology( topology_ ),
     m_cellVec( topology_.totalGeomModules()),
     m_validGeomIds( topology_.totalGeomModules()),
     m_halfType( topology_.detectorType()),
     m_subdet( topology_.subDetector()) {
   
   m_validIds.reserve( topology().totalModules());
 #ifdef EDM_ML_DEBUG
   std::cout << "Expected total # of Geometry Modules " 
         << topology().totalGeomModules() << std::endl;
 #endif
 }

HGCalGeometry::~HGCalGeometry ( )

override

Definition at line 37 of file HGCalGeometry.cc.

Referenced by dbString().

37 { }

Member Function Documentation

void HGCalGeometry::addValidID ( const DetId & id )

protected

Definition at line 299 of file HGCalGeometry.cc.

Referenced by topology().

                                               {
   edm::LogError("HGCalGeom") << "HGCalGeometry::addValidID is not implemented";
 }

std::string HGCalGeometry::cellElement ( ) const

Definition at line 239 of file HGCalGeometry.cc.

References HGCEE, HGCHEB, HGCHEF, and m_subdet.

Referenced by producerTag().

                                            {
   if      (m_subdet == HGCEE)  return "HGCalEE";
   else if (m_subdet == HGCHEF) return "HGCalHEFront";
   else if (m_subdet == HGCHEB) return "HGCalHEBack";
   else                         return "Unknown";
 }

std::shared_ptr< const CaloCellGeometry > HGCalGeometry::cellGeomPtr ( uint32_t index ) const

overrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 277 of file HGCalGeometry.cc.

References diffTreeTool::index, m_cellVec, and m_validGeomIds.

Referenced by cellGeomPtr(), getGeometry(), getSummary(), and topology().

                                                                                      {
   if ((index >= m_cellVec.size()) || (m_validGeomIds[index].rawId() == 0)) 
     return nullptr;
   static const auto do_not_delete = [](const void*){};
   auto cell = std::shared_ptr<const CaloCellGeometry>(&m_cellVec[index],do_not_delete);
   if (nullptr == cell->param()) return nullptr;
   return cell;
 }

std::shared_ptr< const CaloCellGeometry > HGCalGeometry::cellGeomPtr	(	uint32_t	index,
		const GlobalPoint &	p
	)		const

private

Definition at line 286 of file HGCalGeometry.cc.

References cellGeomPtr(), diffTreeTool::index, m_cellVec, and m_validGeomIds.

                                                                                                              {
   if ((index >= m_cellVec.size()) || (m_validGeomIds[index].rawId() == 0)) 
     return nullptr;
   if (pos == GlobalPoint()) return cellGeomPtr(index);
   auto cell = std::make_shared<FlatTrd>(m_cellVec[index]);
   cell->setPosition(pos);
 #ifdef EDM_ML_DEBUG
 //std::cout << "cellGeomPtr " << newcell << ":" << cell << std::endl;
 #endif
   if (nullptr == cell->param()) return nullptr;
   return cell;
 }

static std::string HGCalGeometry::dbString ( )

inlinestatic

Definition at line 45 of file HGCalGeometry.h.

References connectstrParser::f1, connectstrParser::f2, connectstrParser::f3, getCorners(), getGeometry(), getPosition(), getSummary(), HGCalGeometry(), mps_fire::i, localCorners(), newCell(), present(), MetAnalyzer::pv(), topology(), and ~HGCalGeometry().

Referenced by CaloGeometryDBEP< T, U >::produceAligned().

45 { return "PHGCalRcd" ; }

void HGCalGeometry::fillNamedParams ( DDFilteredView fv )

virtual

Definition at line 39 of file HGCalGeometry.cc.

Referenced by getValidGeomDetIds().

39 {}

HGCalGeometry::DetIdSet HGCalGeometry::getCells	(	const GlobalPoint &	r,
		double	dR
	)		const

overridevirtual

Get a list of all cells within a dR of the given cell.

The default implementation makes a loop over all cell geometries. Cleverer implementations are suggested to use rough conversions between eta/phi and ieta/iphi and test on the boundaries.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 234 of file HGCalGeometry.cc.

Referenced by getValidGeomDetIds().

                                                                                     {
    HGCalGeometry::DetIdSet dss;
    return dss;
 }

DetId HGCalGeometry::getClosestCell ( const GlobalPoint & r ) const

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 198 of file HGCalGeometry.cc.

References HGCalDDDConstants::assignCell(), gather_cfg::cout, HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalTopology::encode(), HGCalDDDConstants::geomMode(), getClosestCellIndex(), HGCalTopology::DecodedDetId::iCell, HGCalTopology::DecodedDetId::iLay, HGCalTopology::DecodedDetId::iSec, HGCalTopology::DecodedDetId::iSubSec, DTRecHitQuality_cfi::local, m_cellVec, m_validGeomIds, HGCalGeometryMode::Square, topology(), HGCalDDDConstants::waferTypeT(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), PV3DBase< T, PVType, FrameType >::z(), and HGCalTopology::DecodedDetId::zside.

Referenced by getValidGeomDetIds().

                                                               {
   unsigned int cellIndex = getClosestCellIndex(r);
   if (cellIndex < m_cellVec.size()) {
     HGCalTopology::DecodedDetId id_ = topology().decode(m_validGeomIds[cellIndex]);
     HepGeom::Point3D<float> local;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       local = m_cellVec[cellIndex].getLocal(r);
     } else {
       if (r.z() > 0) local = HepGeom::Point3D<float>(r.x(),r.y(),0);
       else           local = HepGeom::Point3D<float>(-r.x(),r.y(),0);
     }
     std::pair<int,int> kxy = 
       topology().dddConstants().assignCell(local.x(),local.y(),id_.iLay,
                        id_.iSubSec,true);
     id_.iCell   = kxy.second;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       id_.iSubSec = kxy.first;
     } else {
       id_.iSec    = kxy.first;
       id_.iSubSec = topology().dddConstants().waferTypeT(kxy.first);
       if (id_.iSubSec != 1) id_.iSubSec = -1;
     }
 #ifdef EDM_ML_DEBUG
     std::cout << "getClosestCell: local " << local << " Id " << id_.zside 
           << ":" << id_.iLay << ":" << id_.iSec << ":" << id_.iSubSec
           << ":" << id_.iCell << std::endl;
 #endif
 
     //check if returned cell is valid
     if (id_.iCell>=0) return topology().encode(id_);
   }
 
   //if not valid or out of bounds return a null DetId
   return DetId();
 }

unsigned int HGCalGeometry::getClosestCellIndex ( const GlobalPoint & r ) const

protected

Definition at line 304 of file HGCalGeometry.cc.

References a, b, gather_cfg::cout, PVValHelper::dz, getPosition(), gen::k, m_cellVec, M_PI, PV3DBase< T, PVType, FrameType >::phi(), DetId::rawId(), z, and PV3DBase< T, PVType, FrameType >::z().

Referenced by getClosestCell(), and topology().

                                                                            {
 
   float phip = r.phi();
   float zp   = r.z();
   unsigned int cellIndex =  m_cellVec.size();
   float dzmin(9999), dphimin(9999), dphi10(0.175);
   for (unsigned int k=0; k<m_cellVec.size(); ++k) {
     float dphi = phip-m_cellVec[k].phiPos();
     while (dphi >   M_PI) dphi -= 2*M_PI;
     while (dphi <= -M_PI) dphi += 2*M_PI;
     if (fabs(dphi) < dphi10) {
       float dz = fabs(zp - m_cellVec[k].getPosition().z());
       if (dz < (dzmin+0.001)) {
     dzmin     = dz;
     if (fabs(dphi) < (dphimin+0.01)) {
       cellIndex = k;
       dphimin   = fabs(dphi);
     } else {
       if (cellIndex >= m_cellVec.size()) cellIndex = k;
     }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   std::cout << "getClosestCellIndex::Input " << zp << ":" << phip << " Index "
         << cellIndex;
   if (cellIndex < m_cellVec.size()) 
     std::cout << " Cell z " << m_cellVec[cellIndex].getPosition().z() << ":"
           << dzmin << " phi " << m_cellVec[cellIndex].phiPos() << ":"
           << dphimin;
   std::cout << std::endl;
 
 #endif
   return cellIndex;
 }

HGCalGeometry::CornersVec HGCalGeometry::getCorners ( const DetId & id ) const

Returns the corner points of this cell's volume.

Definition at line 173 of file HGCalGeometry.cc.

References HGCalTopology::dddConstants(), HGCalTopology::decode(), PVValHelper::dx, PVValHelper::dz, mps_fire::i, HGCalTopology::DecodedDetId::iCell, HGCalTopology::DecodedDetId::iLay, indexFor(), HGCalTopology::DecodedDetId::iSec, HGCalTopology::DecodedDetId::iSubSec, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), m_cellVec, HGCalGeometryMode::Square, topology(), and geometryCSVtoXML::xy.

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), and dbString().

                                                                        {
 
   HGCalGeometry::CornersVec co (8, GlobalPoint(0,0,0));
   unsigned int cellIndex =  indexFor(id);
   if (cellIndex <  m_cellVec.size()) {
     HGCalTopology::DecodedDetId id_ = topology().decode(id);
     std::pair<float,float> xy;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       xy = topology().dddConstants().locateCell(id_.iCell,id_.iLay,id_.iSubSec,true);
     } else {
       xy = topology().dddConstants().locateCellHex(id_.iCell,id_.iSec,true);
     }
     float dz = m_cellVec[cellIndex].param()[0];
     float dx = 0.5*m_cellVec[cellIndex].param()[11];
     static const int signx[] = {-1,-1,1,1,-1,-1,1,1};
     static const int signy[] = {-1,1,1,-1,-1,1,1,-1};
     static const int signz[] = {-1,-1,-1,-1,1,1,1,1};
     for (unsigned int i = 0; i != 8; ++i) {
       const HepGeom::Point3D<float> lcoord(xy.first+signx[i]*dx,xy.second+signy[i]*dx,signz[i]*dz);
       co[i] = m_cellVec[cellIndex].getPosition(lcoord);
     }
   }
   return co;
 }

std::shared_ptr< const CaloCellGeometry > HGCalGeometry::getGeometry ( const DetId & id ) const

overridevirtual

Get the cell geometry of a given detector id. Should return false if not found.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 120 of file HGCalGeometry.cc.

References cellGeomPtr(), HGCEEDetId::geometryCell(), HGCHEDetId::geometryCell(), HGCalDetId::geometryCell(), getPosition(), HGCEE, HGCalGeometryMode::Square, and topology().

Referenced by dbString().

                                                                                       {
   if (id == DetId()) return nullptr; // nothing to get
   DetId geoId;
   if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
     geoId = (id.subdetId() == HGCEE ? 
          (DetId)(HGCEEDetId(id).geometryCell()) : 
          (DetId)(HGCHEDetId(id).geometryCell()));
   } else {
     geoId = (DetId)(HGCalDetId(id).geometryCell());
   }
   const uint32_t cellIndex (topology().detId2denseGeomId(geoId));
   const GlobalPoint pos = (id != geoId) ? getPosition(id) : GlobalPoint();
   return cellGeomPtr (cellIndex, pos);
 
 }

const CaloCellGeometry * HGCalGeometry::getGeometryRawPtr ( uint32_t index ) const

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 270 of file HGCalGeometry.cc.

References m_cellVec, and CaloCellGeometry::param().

Referenced by present(), and topology().

                                                                              {
   // Modify the RawPtr class
   const CaloCellGeometry* cell(&m_cellVec[index]);
   return (m_cellVec.size() < index ||
       nullptr == cell->param() ? nullptr : cell);
 }

GlobalPoint HGCalGeometry::getPosition ( const DetId & id ) const

Definition at line 150 of file HGCalGeometry.cc.

References gather_cfg::cout, HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalTopology::DecodedDetId::iCell, HGCalTopology::DecodedDetId::iLay, indexFor(), HGCalTopology::DecodedDetId::iSec, HGCalTopology::DecodedDetId::iSubSec, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), m_cellVec, HGCalGeometryMode::Square, topology(), and geometryCSVtoXML::xy.

Referenced by HGCalTBAnalyzer::analyzeDigi(), HGCalTBAnalyzer::analyzeRecHits(), HGCalTriggerGeometryHexImp1::buildTriggerCellsAndModules(), dbString(), HGCalTriggerNtupleHGCDigis::fill(), getClosestCellIndex(), getGeometry(), HGCalTriggerGeometryHexImp2::getModulePosition(), HGCalTriggerGeometryHexLayerBasedImp1::getModulePosition(), hgcal::RecHitTools::getPosition(), HGCalTriggerGeometryHexImp2::getTriggerCellPosition(), HGCalTriggerGeometryHexLayerBasedImp1::getTriggerCellPosition(), and HGCalTriggerGeometryImp1::initialize().

                                                             {
 
   unsigned int cellIndex =  indexFor(id);
   GlobalPoint glob;
   if (cellIndex <  m_cellVec.size()) {
     HGCalTopology::DecodedDetId id_ = topology().decode(id);
     std::pair<float,float> xy;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       xy = topology().dddConstants().locateCell(id_.iCell,id_.iLay,id_.iSubSec,true);
     } else {
       xy = topology().dddConstants().locateCellHex(id_.iCell,id_.iSec,true);
     }
     const HepGeom::Point3D<float> lcoord(xy.first,xy.second,0);
     glob = m_cellVec[cellIndex].getPosition(lcoord);
 #ifdef EDM_ML_DEBUG
     std::cout << "getPosition:: index " << cellIndex << " Local " << lcoord.x()
           << ":" << lcoord.y() << " ID " << id_.iCell << ":" << id_.iLay 
           << " Global " << glob << std::endl;
 #endif
   } 
   return glob;
 }

void HGCalGeometry::getSummary	(	CaloSubdetectorGeometry::TrVec &	trVector,
		CaloSubdetectorGeometry::IVec &	iVector,
		CaloSubdetectorGeometry::DimVec &	dimVector,
		CaloSubdetectorGeometry::IVec &	dinsVector
	)		const

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 354 of file HGCalGeometry.cc.

References HGCalParameters::hgtrap::alpha, HGCalParameters::hgtrap::bl, cellGeomPtr(), HGCalParameters::hgtrap::cellSize, HGCalTopology::dddConstants(), HGCalParameters::hgtrap::dz, HGCalDDDConstants::getModule(), HGCalDDDConstants::getTrForm(), HGCalDDDConstants::getTrFormN(), runTauDisplay::gp, HGCalParameters::hgtrap::h, HGCEE, mps_fire::i, gen::k, k_NumberOfParametersPerShape, k_NumberOfShapes, HGCalParameters::hgtrform::lay, HGCEEDetId::layer(), HGCHEDetId::layer(), m_validGeomIds, CaloSubdetectorGeometry::numberOfParametersPerShape(), CaloSubdetectorGeometry::numberOfShapes(), CaloSubdetectorGeometry::numberOfTransformParms(), findQualityFiles::rr, HGCalDDDConstants::sectors(), HGCalGeometryMode::Square, DetId::subdetId(), HGCalParameters::hgtrap::tl, topology(), HGCalTopology::totalGeomModules(), groupFilesInBlocks::tt, HGCalDDDConstants::volumes(), and HGCalDDDConstants::waferInLayer().

Referenced by dbString(), and CaloGeometryDBEP< T, U >::produceAligned().

                                                                    {
 
   unsigned int numberOfCells = topology().totalGeomModules(); // total Geom Modules both sides
   unsigned int numberOfShapes = HGCalGeometry::k_NumberOfShapes;
   unsigned int numberOfParametersPerShape = HGCalGeometry::k_NumberOfParametersPerShape;
 
   trVector.reserve( numberOfCells * numberOfTransformParms());
   iVector.reserve( numberOfCells );
   dimVector.reserve( numberOfShapes * numberOfParametersPerShape );
   dinsVector.reserve( numberOfCells );
   
   if (topology().geomMode() == HGCalGeometryMode::Square) {
     for (unsigned int k=0; k <topology().dddConstants().volumes(); ++k) {
       HGCalParameters::hgtrap vol = topology().dddConstants().getModule(k,false,true);
       ParmVec params( HGCalGeometry::k_NumberOfParametersPerShape, 0 );
       params[0] = vol.dz;
       params[1] = params[2] = 0;
       params[3] = params[7] = vol.h;
       params[4] = params[8] = vol.bl;
       params[5] = params[9] = vol.tl;
       params[6] = params[10]= vol.alpha;
       params[11]= vol.cellSize;
       dimVector.insert( dimVector.end(), params.begin(), params.end());
     }
   } else {
     for (unsigned itr=0; itr<topology().dddConstants().getTrFormN(); ++itr) {
       HGCalParameters::hgtrform mytr = topology().dddConstants().getTrForm(itr);
       int layer  = mytr.lay;
       for (int wafer=0; wafer<topology().dddConstants().sectors(); ++wafer) {
         if (topology().dddConstants().waferInLayer(wafer,layer,true)) {
       HGCalParameters::hgtrap vol = topology().dddConstants().getModule(wafer, true, true);
       ParmVec params( HGCalGeometry::k_NumberOfParametersPerShape, 0 );
       params[0] = vol.dz;
       params[1] = params[2] = 0;
       params[3] = params[7] = vol.h;
       params[4] = params[8] = vol.bl;
       params[5] = params[9] = vol.tl;
       params[6] = params[10]= vol.alpha;
       params[11]= vol.cellSize;
       dimVector.insert( dimVector.end(), params.begin(), params.end());
     }
       }
     }
   }
   
   for (unsigned int i( 0 ); i < numberOfCells; ++i) {
     DetId detId = m_validGeomIds[i];
     int layer = ((detId.subdetId() ==  ForwardSubdetector::HGCEE) ?
          (HGCEEDetId(detId).layer()) :
          (HGCHEDetId(detId).layer()));
     dinsVector.emplace_back( topology().detId2denseGeomId( detId ));
     iVector.emplace_back( layer );
     
     Tr3D tr;
     auto ptr =  cellGeomPtr( i );
     if ( nullptr != ptr ) {
       ptr->getTransform( tr, ( Pt3DVec* ) nullptr );
 
       if( Tr3D() == tr ) { // there is no rotation
     const GlobalPoint& gp( ptr->getPosition()); 
     tr = HepGeom::Translate3D( gp.x(), gp.y(), gp.z());
       }
 
       const CLHEP::Hep3Vector tt( tr.getTranslation());
       trVector.emplace_back( tt.x());
       trVector.emplace_back( tt.y());
       trVector.emplace_back( tt.z());
       if (6 == numberOfTransformParms()) {
     const CLHEP::HepRotation rr( tr.getRotation());
     const ROOT::Math::Transform3D rtr( rr.xx(), rr.xy(), rr.xz(), tt.x(),
                        rr.yx(), rr.yy(), rr.yz(), tt.y(),
                        rr.zx(), rr.zy(), rr.zz(), tt.z());
     ROOT::Math::EulerAngles ea;
     rtr.GetRotation( ea );
     trVector.emplace_back( ea.Phi());
     trVector.emplace_back( ea.Theta());
     trVector.emplace_back( ea.Psi());
       }
     }
   }
 }

const std::vector<DetId>& HGCalGeometry::getValidDetIds	(	DetId::Detector	det = `DetId::Detector(0)`,
		int	subdet = `0`
	)		const

inlineoverridevirtual

Get a list of valid detector ids (for the given subdetector)

Note: The implementation in this class is relevant for SubdetectorGeometries which handle only a single subdetector at a time. It does not look at the det and subdet arguments.

Reimplemented from CaloSubdetectorGeometry.

Definition at line 78 of file HGCalGeometry.h.

References CaloSubdetectorGeometry::m_validIds.

78 { return m_validIds; }

CaloSubdetectorGeometry::m_validIds

std::vector< DetId > m_validIds

Definition: CaloSubdetectorGeometry.h:126

const std::vector<DetId>& HGCalGeometry::getValidGeomDetIds ( void ) const

inline

Definition at line 79 of file HGCalGeometry.h.

References PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, fillNamedParams(), getCells(), getClosestCell(), initializeParms(), m_validGeomIds, and alignCSCRings::r.

79 { return m_validGeomIds; }

HGCalGeometry::m_validGeomIds

std::vector< DetId > m_validGeomIds

Definition: HGCalGeometry.h:122

unsigned int HGCalGeometry::indexFor ( const DetId & id ) const

overrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 246 of file HGCalGeometry.cc.

References gather_cfg::cout, TauDecayModes::dec, HGCalTopology::detId2denseGeomId(), HGCEEDetId::geometryCell(), HGCHEDetId::geometryCell(), HGCalDetId::geometryCell(), HGCEE, m_cellVec, DetId::rawId(), HGCalGeometryMode::Square, and topology().

Referenced by getCorners(), getPosition(), and topology().

                                                           {
   unsigned int cellIndex =  m_cellVec.size();
   if (id != DetId()) {
     DetId geoId;
     if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
       geoId = (id.subdetId() == HGCEE ? 
            (DetId)(HGCEEDetId(id).geometryCell()) : 
            (DetId)(HGCHEDetId(id).geometryCell()));
     } else {
       geoId = (DetId)(HGCalDetId(id).geometryCell());
     }
     cellIndex = topology().detId2denseGeomId(geoId);
 #ifdef EDM_ML_DEBUG
     std::cout << "indexFor " << std::hex << id.rawId() << ":" << geoId.rawId()
           << std::dec << " index " << cellIndex << std::endl;
 #endif
   }
   return cellIndex;
 }

void HGCalGeometry::initializeParms ( )

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 41 of file HGCalGeometry.cc.

Referenced by getValidGeomDetIds().

41 {

42 }

void HGCalGeometry::localCorners	(	Pt3DVec &	lc,
		const CCGFloat *	pv,
		unsigned int	i,
		Pt3D &	ref
	)

Definition at line 44 of file HGCalGeometry.cc.

References FlatTrd::localCorners().

Referenced by dbString().

                                       {
   FlatTrd::localCorners( lc, pv, ref ) ;
 }

void HGCalGeometry::newCell	(	const GlobalPoint &	f1,
		const GlobalPoint &	f2,
		const GlobalPoint &	f3,
		const CCGFloat *	parm,
		const DetId &	detId
	)

overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 51 of file HGCalGeometry.cc.

References CaloSubdetectorGeometry::cornersMgr(), gather_cfg::cout, HGCalTopology::dddConstants(), TauDecayModes::dec, HGCalTopology::decode(), alcaDQMUpload::encode(), HGCalTopology::encode(), HGCEEDetId::geometryCell(), HGCHEDetId::geometryCell(), HGCalDetId::geometryCell(), HGCEE, m_cellVec, m_halfType, m_subdet, m_validGeomIds, CaloSubdetectorGeometry::m_validIds, HGCalDDDConstants::maxCells(), HGCalDDDConstants::numberCellsHexagon(), DetId::rawId(), HGCalGeometryMode::Square, DetId::subdetId(), topology(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by HGCalGeometryLoader::buildGeom(), and dbString().

                                             {
 
   DetId geomId;
   int   cells;
   HGCalTopology::DecodedDetId id = topology().decode(detId);
   if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
     geomId = (detId.subdetId() == HGCEE ? 
           (DetId)(HGCEEDetId(detId).geometryCell()) :
           (DetId)(HGCHEDetId(detId).geometryCell()));
     cells = topology().dddConstants().maxCells(id.iLay,true);
   } else {
     geomId = (DetId)(HGCalDetId(detId).geometryCell());
     cells = topology().dddConstants().numberCellsHexagon(id.iSec);
 #ifdef EDM_ML_DEBUG
     std::cout << "NewCell " << HGCalDetId(detId) << " GEOM " 
           << HGCalDetId(geomId) << std::endl;
 #endif
   }
   const uint32_t cellIndex (topology().detId2denseGeomId(detId));
 
   m_cellVec.at( cellIndex ) = FlatTrd( cornersMgr(), f1, f2, f3, parm ) ;
   m_validGeomIds.at( cellIndex ) = geomId ;
 
 #ifdef EDM_ML_DEBUG
   std::cout << "Store for DetId " << std::hex << detId.rawId() << " GeomId " 
         << geomId.rawId() << std::dec << " Index " << cellIndex
         << " cells " << cells << std::endl;
   unsigned int nOld = m_validIds.size();
 #endif
   for (int cell = 0; cell < cells; ++cell) {
     id.iCell  = cell;
     DetId idc = topology().encode(id);
     if (topology().valid(idc)) {
       m_validIds.emplace_back(idc);
       if ((topology().dddConstants().geomMode() == HGCalGeometryMode::Square) &&
       (!m_halfType)) {
     id.iSubSec = -id.iSubSec;
     m_validIds.emplace_back( topology().encode(id));
     id.iSubSec = -id.iSubSec;
       }
     }
   }
 
 #ifdef EDM_ML_DEBUG
   std::cout << "HGCalGeometry::newCell-> [" << cellIndex << "]"
         << " front:" << f1.x() << '/' << f1.y() << '/' << f1.z() 
             << " back:" <<  f2.x() << '/' << f2.y() << '/' << f2.z()
         << " eta|phi " << m_cellVec[cellIndex].etaPos() << ":"
         << m_cellVec[cellIndex].phiPos() << " id:";
   if (topology().dddConstants().geomMode() != HGCalGeometryMode::Square) {
     std::cout << HGCalDetId(detId);
   } else if (m_subdet == HGCEE) {
     std::cout << HGCEEDetId(detId);
   } else {
     std::cout << HGCHEDetId(detId);
   }
   unsigned int nNew = m_validIds.size();
   std::cout << " with valid DetId from " << nOld << " to " << nNew
         << std::endl; 
   std::cout << "Cell[" << cellIndex << "] " << std::hex << geomId.rawId() 
         << ":"  << m_validGeomIds[cellIndex].rawId() << std::dec
         << std::endl;
 #endif
 }

bool HGCalGeometry::present ( const DetId & id ) const

overridevirtual

is this detid present in the geometry?

Reimplemented from CaloSubdetectorGeometry.

Definition at line 136 of file HGCalGeometry.cc.

References HGCEEDetId::geometryCell(), HGCHEDetId::geometryCell(), HGCalDetId::geometryCell(), getGeometryRawPtr(), HGCEE, diffTreeTool::index, HGCalGeometryMode::Square, and topology().

Referenced by dbString().

                                                  {
   if (id == DetId()) return false;
   DetId geoId;
   if (topology().dddConstants().geomMode() == HGCalGeometryMode::Square) {
     geoId = (id.subdetId() == HGCEE ? 
          (DetId)(HGCEEDetId(id).geometryCell()) : 
          (DetId)(HGCHEDetId(id).geometryCell()));
   } else {
     geoId = (DetId)(HGCalDetId(id).geometryCell());
   }
   const uint32_t index (topology().detId2denseGeomId(geoId));
   return (nullptr != getGeometryRawPtr(index)) ;
 }

static std::string HGCalGeometry::producerTag ( )

inlinestatic

Definition at line 95 of file HGCalGeometry.h.

References cellElement(), and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by PCaloGeometryBuilder::beginRun(), and CaloGeometryBuilder::produceAligned().

95 { return "HGCal" ; }

unsigned int HGCalGeometry::sizeForDenseIndex ( ) const

protected

Definition at line 266 of file HGCalGeometry.cc.

References topology(), and HGCalTopology::totalGeomModules().

Referenced by topology().

                                                     {
   return topology().totalGeomModules();
 }

void HGCalGeometry::sortDetIds ( void )

Definition at line 349 of file HGCalGeometry.cc.

References CaloSubdetectorGeometry::m_validIds.

Referenced by HGCalGeometryLoader::build(), and topology().

                                      {
   m_validIds.shrink_to_fit();
   std::sort( m_validIds.begin(), m_validIds.end(), rawIdSort());
 }

const HGCalTopology& HGCalGeometry::topology ( ) const

inline

Definition at line 98 of file HGCalGeometry.h.

References addValidID(), cellGeomPtr(), getClosestCellIndex(), getGeometryRawPtr(), diffTreeTool::index, indexFor(), m_topology, AlCaHLTBitMon_ParallelJobs::p, sizeForDenseIndex(), CaloSubdetectorGeometry::sizeForDenseIndex(), and sortDetIds().

Referenced by PFRecHitQTestThresholdInThicknessNormalizedMIPs::beginEvent(), CaloTruthAccumulator::beginLuminosityBlock(), dbString(), HGCalTriggerGeometryBase::eeTopology(), HGCalTriggerGeometryBase::fhTopology(), getClosestCell(), getCorners(), getGeometry(), getPosition(), getSummary(), HGCalGeometry(), indexFor(), newCell(), present(), CaloGeometryBuilder::produceAligned(), HGCalRecHitWorkerSimple::set(), HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::setGeometry(), and sizeForDenseIndex().

98 {return m_topology;}

HGCalGeometry::m_topology

const HGCalTopology & m_topology

Definition: HGCalGeometry.h:119

Member Data Documentation

CellVec HGCalGeometry::m_cellVec

private

Definition at line 121 of file HGCalGeometry.h.

Referenced by cellGeomPtr(), getClosestCell(), getClosestCellIndex(), getCorners(), getGeometryRawPtr(), getPosition(), indexFor(), and newCell().

bool HGCalGeometry::m_halfType

private

Definition at line 123 of file HGCalGeometry.h.

Referenced by newCell().

ForwardSubdetector HGCalGeometry::m_subdet

private

Definition at line 124 of file HGCalGeometry.h.

Referenced by cellElement(), and newCell().

const HGCalTopology& HGCalGeometry::m_topology

private

Definition at line 119 of file HGCalGeometry.h.

Referenced by topology().

std::vector<DetId> HGCalGeometry::m_validGeomIds

private

Definition at line 122 of file HGCalGeometry.h.

Referenced by cellGeomPtr(), getClosestCell(), getSummary(), getValidGeomDetIds(), and newCell().

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation