#include <HGCalGeometry.h>

Inheritance diagram for HGCalGeometry:

Public Types
typedef HGCalGeometryRecord	AlignedRecord

typedef CaloCellGeometry::CCGFloat	CCGFloat

typedef std::vector< FlatHexagon >	CellVec

typedef std::vector< FlatTrd >	CellVec2

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)

CornersVec	get8Corners (const DetId &id) const

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

Static Public Attributes
static unsigned int	k_NumberOfParametersPerHex = 3

static unsigned int	k_NumberOfParametersPerShape = 3

static unsigned int	k_NumberOfParametersPerTrd = 12

static unsigned int	k_NumberOfShapes = 100

static unsigned int	k_NumberOfShapesTrd = 1000

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

template<class T >
unsigned int	getClosestCellIndex (const GlobalPoint &r, const std::vector< T > &vec) const

DetId	getGeometryDetId (DetId detId) const

Private Attributes
CellVec	m_cellVec

CellVec2	m_cellVec2

DetId::Detector	m_det

ForwardSubdetector	m_subdet

const HGCalTopology &	m_topology

std::vector< DetId >	m_validGeomIds

HGCalGeometryMode::GeometryMode	mode_

const double	twoBysqrt3_

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 28 of file HGCalGeometry.h.

Member Typedef Documentation

typedef HGCalGeometryRecord HGCalGeometry::AlignedRecord

Definition at line 41 of file HGCalGeometry.h.

typedef CaloCellGeometry::CCGFloat HGCalGeometry::CCGFloat

Definition at line 34 of file HGCalGeometry.h.

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

Definition at line 32 of file HGCalGeometry.h.

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

Definition at line 33 of file HGCalGeometry.h.

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

Definition at line 39 of file HGCalGeometry.h.

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

Definition at line 38 of file HGCalGeometry.h.

typedef PHGCalRcd HGCalGeometry::PGeometryRecord

Definition at line 42 of file HGCalGeometry.h.

typedef CaloCellGeometry::Pt3D HGCalGeometry::Pt3D

Definition at line 35 of file HGCalGeometry.h.

typedef CaloCellGeometry::Pt3DVec HGCalGeometry::Pt3DVec

Definition at line 36 of file HGCalGeometry.h.

Constructor & Destructor Documentation

HGCalGeometry::HGCalGeometry ( const HGCalTopology & topology )

Definition at line 23 of file HGCalGeometry.cc.

References DetId::HGCalHSc, m_cellVec, m_cellVec2, m_det, m_topology, CaloSubdetectorGeometry::m_validIds, HGCalTopology::totalGeomModules(), and HGCalTopology::totalModules().

Referenced by dbString().

   : m_topology( topology_ ),
     m_validGeomIds( topology_.totalGeomModules()),
     mode_( topology_.geomMode()),
     m_det( topology_.detector()),
     m_subdet( topology_.subDetector()),
     twoBysqrt3_(2.0/std::sqrt(3.0)) {
   
   if (m_det == DetId::HGCalHSc) {
     m_cellVec2 = CellVec2(topology_.totalGeomModules());
   } else {
     m_cellVec = CellVec(topology_.totalGeomModules());
   }
   m_validIds.reserve( m_topology.totalModules());
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "Expected total # of Geometry Modules " 
                 << m_topology.totalGeomModules();
 #endif
 }

HGCalGeometry::~HGCalGeometry ( )

override

Definition at line 43 of file HGCalGeometry.cc.

Referenced by dbString().

43 { }

Member Function Documentation

void HGCalGeometry::addValidID ( const DetId & id )

protected

Definition at line 493 of file HGCalGeometry.cc.

Referenced by topology().

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

std::string HGCalGeometry::cellElement ( ) const

Definition at line 414 of file HGCalGeometry.cc.

References DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, HGCEE, HGCHEB, HGCHEF, m_det, and m_subdet.

Referenced by producerTag().

                                            {
   if      (m_subdet == HGCEE  || m_det == DetId::HGCalEE)  return "HGCalEE";
   else if (m_subdet == HGCHEF || m_det == DetId::HGCalHSi) return "HGCalHEFront";
   else if (m_subdet == HGCHEB || m_det == DetId::HGCalHSc) return "HGCalHEBack";
   else                                                     return "Unknown";
 }

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

overrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 453 of file HGCalGeometry.cc.

References DetId::HGCalHSc, m_cellVec, m_cellVec2, m_det, and m_validGeomIds.

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

                                                                                      {
   if ((index >= m_cellVec.size()  && m_det != DetId::HGCalHSc) ||
       (index >= m_cellVec2.size() && m_det == DetId::HGCalHSc) ||
       (m_validGeomIds[index].rawId() == 0)) return nullptr;
   static const auto do_not_delete = [](const void*){};
   if (m_det == DetId::HGCalHSc) {
     auto cell = std::shared_ptr<const CaloCellGeometry>(&m_cellVec2[index],do_not_delete);
     if (nullptr == cell->param()) return nullptr;
     return cell;
   } else {
     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 469 of file HGCalGeometry.cc.

References cellGeomPtr(), DetId::HGCalHSc, m_cellVec, m_cellVec2, m_det, and m_validGeomIds.

                                                                                                              {
   if ((index >= m_cellVec.size()  && m_det != DetId::HGCalHSc) ||
       (index >= m_cellVec2.size() && m_det == DetId::HGCalHSc) ||
       (m_validGeomIds[index].rawId() == 0)) return nullptr;
   if (pos == GlobalPoint()) return cellGeomPtr(index);
   if (m_det == DetId::HGCalHSc) {
     auto cell = std::make_shared<FlatTrd>(m_cellVec2[index]);
     cell->setPosition(pos);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "cellGeomPtr " << index << ":" << cell;
 #endif
     if (nullptr == cell->param()) return nullptr;
     return cell;
   } else {
     auto cell = std::make_shared<FlatHexagon>(m_cellVec[index]);
     cell->setPosition(pos);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "cellGeomPtr " << index << ":" << cell;
 #endif
     if (nullptr == cell->param()) return nullptr;
     return cell;
   }
 }

static std::string HGCalGeometry::dbString ( )

inlinestatic

Definition at line 50 of file HGCalGeometry.h.

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

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

50 { return "PHGCalRcd" ; }

void HGCalGeometry::fillNamedParams ( DDFilteredView fv )

virtual

Definition at line 45 of file HGCalGeometry.cc.

Referenced by getValidGeomDetIds().

45 {}

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

Definition at line 318 of file HGCalGeometry.cc.

References HGCalTopology::dddConstants(), HGCalTopology::decode(), PVValHelper::dx, PVValHelper::dz, getPosition(), HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, DetId::HGCalHSc, mps_fire::i, HGCalTopology::DecodedDetId::iCell1, HGCalTopology::DecodedDetId::iCell2, HGCalTopology::DecodedDetId::iLay, indexFor(), HGCalTopology::DecodedDetId::iSec1, HGCalTopology::DecodedDetId::iSec2, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), HGCalDDDConstants::locateCellTrap(), m_cellVec, m_cellVec2, m_det, m_topology, mode_, FlatTrd::ncorner_, HGCalGeometryMode::Trapezoid, and geometryCSVtoXML::xy.

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

                                                                         {
   
   unsigned int ncorner = FlatTrd::ncorner_;
   HGCalGeometry::CornersVec co (ncorner, GlobalPoint(0,0,0));
   unsigned int cellIndex =  indexFor(id);
   if ((cellIndex <  m_cellVec.size() && m_det != DetId::HGCalHSc) ||
       (cellIndex <  m_cellVec2.size() && m_det == DetId::HGCalHSc)) {
     HGCalTopology::DecodedDetId id_ = m_topology.decode(id);
     std::pair<float,float> xy;
     if ((mode_ == HGCalGeometryMode::Hexagon) ||
     (mode_ == HGCalGeometryMode::HexagonFull)) {
       xy = m_topology.dddConstants().locateCellHex(id_.iCell1,id_.iSec1,true);
     } else if (mode_ == HGCalGeometryMode::Trapezoid) {
       xy = m_topology.dddConstants().locateCellTrap(id_.iLay,id_.iSec1,
                             id_.iCell1,true);
     } else {
       xy = m_topology.dddConstants().locateCell(id_.iLay,id_.iSec1,id_.iSec2,
                         id_.iCell1,id_.iCell2,true,false);
     }
     float dx = ((m_det == DetId::HGCalHSc) ? 0.5*m_cellVec2[cellIndex].param()[11] :
         m_cellVec[cellIndex].param()[1]);
     float dz = ((m_det == DetId::HGCalHSc) ? m_cellVec2[cellIndex].param()[0] :
         m_cellVec[cellIndex].param()[0]);
     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 < ncorner; ++i) {
       const HepGeom::Point3D<float> lcoord(xy.first+signx[i]*dx,
                        xy.second+signy[i]*dx,signz[i]*dz);
       co[i] = ((m_det == DetId::HGCalHSc) ? 
            (m_cellVec2[cellIndex].getPosition(lcoord)) :
            (m_cellVec[cellIndex].getPosition(lcoord)));
     }
   }
   return co;
 }

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 409 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 355 of file HGCalGeometry.cc.

References HGCalDDDConstants::assignCell(), HGCalDDDConstants::assignCellHex(), HGCalDDDConstants::assignCellTrap(), HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalTopology::encode(), getClosestCellIndex(), HGCalDDDConstants::getLayer(), HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, DetId::HGCalHSc, HGCalTopology::DecodedDetId::iCell1, HGCalTopology::DecodedDetId::iCell2, HGCalTopology::DecodedDetId::iLay, HGCalTopology::DecodedDetId::iSec1, HGCalTopology::DecodedDetId::iSec2, HGCalTopology::DecodedDetId::iType, DTRecHitClients_cfi::local, m_cellVec, m_cellVec2, m_det, m_topology, m_validGeomIds, mode_, HGCalGeometryMode::Trapezoid, 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() && m_det != DetId::HGCalHSc) ||
       (cellIndex <  m_cellVec2.size() && m_det == DetId::HGCalHSc)) {
     HGCalTopology::DecodedDetId id_ = m_topology.decode(m_validGeomIds[cellIndex]);
     HepGeom::Point3D<float> local;
     if (r.z() > 0) {
       local    = HepGeom::Point3D<float>(r.x(),r.y(),0);
       id_.zSide = 1;
     } else {
       local = HepGeom::Point3D<float>(-r.x(),r.y(),0);
       id_.zSide =-1;
     }
     if ((mode_ == HGCalGeometryMode::Hexagon) ||
     (mode_ == HGCalGeometryMode::HexagonFull)) {
       const auto & kxy = 
     m_topology.dddConstants().assignCell(local.x(),local.y(),id_.iLay,
                          id_.iType,true);
       id_.iCell1  = kxy.second;
       id_.iSec1   = kxy.first;
       id_.iType   = m_topology.dddConstants().waferTypeT(kxy.first);
       if (id_.iType != 1) id_.iType = -1;
     } else if (mode_ == HGCalGeometryMode::Trapezoid) {
       id_.iLay = m_topology.dddConstants().getLayer(r.z(),true);
       const auto & kxy = 
     m_topology.dddConstants().assignCellTrap(r.x(),r.y(),r.z(),
                          id_.iLay,true);
       id_.iSec1  = kxy[0];
       id_.iCell1 = kxy[1];
       id_.iType  = kxy[2];
     } else {
       id_.iLay = m_topology.dddConstants().getLayer(r.z(),true);
       const auto & kxy = 
     m_topology.dddConstants().assignCellHex(local.x(),local.y(),id_.iLay,
                         true);
       id_.iSec1  = kxy[0]; id_.iSec2  = kxy[1]; id_.iType = kxy[2];
       id_.iCell1 = kxy[3]; id_.iCell2 = kxy[4];
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "getClosestCell: local " << local 
                   << " Id " << id_.zSide << ":" << id_.iLay 
                   << ":" << id_.iSec1 << ":" << id_.iSec2
                   << ":" << id_.iType << ":" << id_.iCell1
                   << ":" << id_.iCell2;
 #endif
 
     //check if returned cell is valid
     if (id_.iCell1>=0) return m_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 497 of file HGCalGeometry.cc.

References DetId::HGCalHSc, m_cellVec, m_cellVec2, and m_det.

Referenced by getClosestCell(), and topology().

                                                                            {
 
   return ((m_det == DetId::HGCalHSc) ? getClosestCellIndex(r,m_cellVec2) : getClosestCellIndex(r,m_cellVec));
 }

template<class T >

unsigned int HGCalGeometry::getClosestCellIndex	(	const GlobalPoint &	r,
		const std::vector< T > &	vec
	)		const

private

Definition at line 503 of file HGCalGeometry.cc.

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

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

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

Returns the corner points of this cell's volume.

Definition at line 272 of file HGCalGeometry.cc.

References HGCalTopology::dddConstants(), HGCalTopology::decode(), PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, DetId::HGCalHSc, mps_fire::i, HGCalTopology::DecodedDetId::iCell1, HGCalTopology::DecodedDetId::iCell2, HGCalTopology::DecodedDetId::iLay, indexFor(), HGCalTopology::DecodedDetId::iSec1, HGCalTopology::DecodedDetId::iSec2, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), HGCalDDDConstants::locateCellTrap(), m_cellVec, m_cellVec2, m_det, m_topology, mode_, FlatTrd::ncorner_, FlatHexagon::ncorner_, HGCalGeometryMode::Trapezoid, and geometryCSVtoXML::xy.

Referenced by dbString().

                                                                        {
 
   unsigned int ncorner = ((m_det == DetId::HGCalHSc) ? FlatTrd::ncorner_ :
               FlatHexagon::ncorner_);
   HGCalGeometry::CornersVec co (ncorner, GlobalPoint(0,0,0));
   unsigned int cellIndex =  indexFor(id);
   if ((cellIndex <  m_cellVec.size() && m_det != DetId::HGCalHSc) ||
       (cellIndex <  m_cellVec2.size() && m_det == DetId::HGCalHSc)) {
     HGCalTopology::DecodedDetId id_ = m_topology.decode(id);
     std::pair<float,float> xy;
     if ((mode_ == HGCalGeometryMode::Hexagon) ||
     (mode_ == HGCalGeometryMode::HexagonFull)) {
       xy = m_topology.dddConstants().locateCellHex(id_.iCell1,id_.iSec1,true);
     } else if (mode_ == HGCalGeometryMode::Trapezoid) {
       xy = m_topology.dddConstants().locateCellTrap(id_.iLay,id_.iSec1,
                             id_.iCell1,true);
     } else {
       xy = m_topology.dddConstants().locateCell(id_.iLay,id_.iSec1,id_.iSec2,
                         id_.iCell1,id_.iCell2,true,false);
     }
     if (m_det == DetId::HGCalHSc) {
       float dx = 0.5*m_cellVec2[cellIndex].param()[11];
       float dz = m_cellVec2[cellIndex].param()[0];
       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 < ncorner; ++i) {
     const HepGeom::Point3D<float> lcoord(xy.first+signx[i]*dx,xy.second+signy[i]*dx,signz[i]*dz);
     co[i] = m_cellVec2[cellIndex].getPosition(lcoord);
       }
     } else {
       float dx = m_cellVec[cellIndex].param()[1];
       float dy = m_cellVec[cellIndex].param()[2];
       float dz = m_cellVec[cellIndex].param()[0];
       static const int signx[] = {0,-1,-1,0,1,1,0,-1,-1,0,1,1};
       static const int signy[] = {-2,-1,1,2,1,-1,-2,-1,1,2,1,-1};
       static const int signz[] = {-1,-1,-1,-1,-1,-1,1,1,1,1,1,1};
       for (unsigned int i = 0; i < ncorner; ++i) {
     const HepGeom::Point3D<float> lcoord(xy.first+signx[i]*dx,xy.second+signy[i]*dy,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 208 of file HGCalGeometry.cc.

References cellGeomPtr(), HGCalTopology::detId2denseGeomId(), getGeometryDetId(), getPosition(), and m_topology.

Referenced by dbString().

                                                                                          {
   if (detId == DetId()) return nullptr; // nothing to get
   DetId geomId = getGeometryDetId(detId);
   const uint32_t cellIndex (m_topology.detId2denseGeomId(geomId));
   const GlobalPoint pos = (detId != geomId) ? getPosition(detId) : GlobalPoint();
   return cellGeomPtr (cellIndex, pos);
 
 }

DetId HGCalGeometry::getGeometryDetId ( DetId detId ) const

private

Definition at line 660 of file HGCalGeometry.cc.

References HGCScintillatorDetId::geometryCell(), HGCalDetId::geometryCell(), HGCSiliconDetId::geometryCell(), HFNoseDetId::geometryCell(), HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, HGCalTopology::isHFNose(), m_topology, mode_, and HGCalGeometryMode::Trapezoid.

Referenced by getGeometry(), indexFor(), newCell(), present(), and topology().

                                                        {
   DetId geomId;
   if ((mode_ == HGCalGeometryMode::Hexagon) || 
       (mode_ == HGCalGeometryMode::HexagonFull)) {
     geomId = static_cast<DetId>(HGCalDetId(detId).geometryCell());
   } else if (mode_ == HGCalGeometryMode::Trapezoid) {
     geomId = static_cast<DetId>(HGCScintillatorDetId(detId).geometryCell());
   } else if (m_topology.isHFNose()) {
     geomId = static_cast<DetId>(HFNoseDetId(detId).geometryCell());
   } else {
     geomId = static_cast<DetId>(HGCSiliconDetId(detId).geometryCell());
   }
   return geomId;
 }

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

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 440 of file HGCalGeometry.cc.

References DetId::HGCalHSc, m_cellVec, m_cellVec2, m_det, and CaloCellGeometry::param().

Referenced by present(), and topology().

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

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

Definition at line 224 of file HGCalGeometry.cc.

References HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, HGCalTopology::DecodedDetId::iCell1, HGCalTopology::DecodedDetId::iCell2, HGCalTopology::DecodedDetId::iLay, indexFor(), HGCalTopology::DecodedDetId::iSec1, HGCalTopology::DecodedDetId::iSec2, HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), m_cellVec, m_cellVec2, m_topology, copyPickMerge_cfg::maxSize, mode_, HGCalGeometryMode::Trapezoid, and geometryCSVtoXML::xy.

                                                             {
 
   unsigned int cellIndex =  indexFor(id);
   GlobalPoint glob;
   unsigned int maxSize = ((mode_ == HGCalGeometryMode::Trapezoid) ? 
               m_cellVec2.size() : m_cellVec.size());
   if (cellIndex <  maxSize) {
     HGCalTopology::DecodedDetId id_ = m_topology.decode(id);
     std::pair<float,float> xy;
     if ((mode_ == HGCalGeometryMode::Hexagon) ||
     (mode_ == HGCalGeometryMode::HexagonFull)) {
       xy = m_topology.dddConstants().locateCellHex(id_.iCell1,id_.iSec1,true);
       const HepGeom::Point3D<float> lcoord(xy.first,xy.second,0);
       glob = m_cellVec[cellIndex].getPosition(lcoord);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCalGeom") << "getPosition:: index " << cellIndex 
                     << " Local " << lcoord.x() << ":" 
                     << lcoord.y() << " ID " << id_.iCell1 
                     << ":" << id_.iSec1 << " Global " << glob;
 #endif
     } else if (mode_ == HGCalGeometryMode::Trapezoid) {
       const HepGeom::Point3D<float> lcoord(0,0,0);
       glob = m_cellVec2[cellIndex].getPosition(lcoord);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCalGeom") << "getPositionTrap:: index " << cellIndex 
                     << " Local " << lcoord.x() << ":" 
                     << lcoord.y() << " ID " << id_.iLay << ":"
                     << id_.iSec1 << ":" << id_.iCell1 
                     << " Global " << glob;
 #endif
     } else {
       xy = m_topology.dddConstants().locateCell(id_.iLay,id_.iSec1,id_.iSec2,
                         id_.iCell1,id_.iCell2,true,false);
       const HepGeom::Point3D<float> lcoord(xy.first,xy.second,0);
       glob = m_cellVec[cellIndex].getPosition(lcoord);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCalGeom") << "getPositionWafer:: index " << cellIndex
                     << " Local " << lcoord.x() << ":" 
                     << lcoord.y() << " ID " << id_.iLay << ":"
                     << id_.iSec1 << ":" << id_.iSec2 << ":"
                     << id_.iCell1 << ":" << id_.iCell2 
                     << " Global " << glob;
 #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 554 of file HGCalGeometry.cc.

References HGCalParameters::hgtrap::alpha, HGCalParameters::hgtrap::bl, cellGeomPtr(), HGCalParameters::hgtrap::cellSize, HGCalTopology::dddConstants(), HGCalTopology::detId2denseGeomId(), HGCalParameters::hgtrap::dz, HGCalParameters::firstModule_, HGCalDDDConstants::getModule(), HGCalDDDConstants::getParameter(), HGCalDDDConstants::getTrForm(), HGCalDDDConstants::getTrFormN(), runTauDisplay::gp, HGCalParameters::hgtrap::h, HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, DetId::HGCalHSc, mps_fire::i, HGCalTopology::isHFNose(), k_NumberOfParametersPerHex, k_NumberOfParametersPerTrd, k_NumberOfShapes, HGCalParameters::lastModule_, HGCalParameters::hgtrform::lay, HGCScintillatorDetId::layer(), HGCalDetId::layer(), HGCSiliconDetId::layer(), HFNoseDetId::layer(), HGCalDDDConstants::layerIndex(), m_det, m_topology, m_validGeomIds, mode_, CaloSubdetectorGeometry::numberOfParametersPerShape(), CaloSubdetectorGeometry::numberOfShapes(), CaloSubdetectorGeometry::numberOfTransformParms(), findQualityFiles::rr, HGCalDDDConstants::sectors(), HGCalParameters::hgtrap::tl, HGCalTopology::totalGeomModules(), HGCalGeometryMode::Trapezoid, groupFilesInBlocks::tt, twoBysqrt3_, and HGCalDDDConstants::waferInLayer().

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

                                                                    {
 
   unsigned int numberOfCells = m_topology.totalGeomModules(); // total Geom Modules both sides
   unsigned int numberOfShapes = k_NumberOfShapes;
   unsigned int numberOfParametersPerShape = 
     ((m_det == DetId::HGCalHSc) ? (unsigned int)(k_NumberOfParametersPerTrd) :
      (unsigned int)(k_NumberOfParametersPerHex));
 
   trVector.reserve( numberOfCells * numberOfTransformParms());
   iVector.reserve( numberOfCells );
   dimVector.reserve( numberOfShapes * numberOfParametersPerShape );
   dinsVector.reserve( numberOfCells );
   
   for (unsigned itr=0; itr<m_topology.dddConstants().getTrFormN(); ++itr) {
     HGCalParameters::hgtrform mytr = m_topology.dddConstants().getTrForm(itr);
     int layer  = mytr.lay;
 
     if ((mode_ == HGCalGeometryMode::Hexagon) ||
     (mode_ == HGCalGeometryMode::HexagonFull)) {
       for (int wafer=0; wafer<m_topology.dddConstants().sectors(); ++wafer) {
         if (m_topology.dddConstants().waferInLayer(wafer,layer,true)) {
       HGCalParameters::hgtrap vol = m_topology.dddConstants().getModule(wafer, true, true);
       ParmVec params( numberOfParametersPerShape, 0 );
       params[0] = vol.dz;
       params[1] = vol.cellSize;
       params[2] = twoBysqrt3_*params[1];
       dimVector.insert( dimVector.end(), params.begin(), params.end());
     }
       }
     } else if (mode_ == HGCalGeometryMode::Trapezoid) {
       int indx = m_topology.dddConstants().layerIndex(layer,true);
       for (int md=m_topology.dddConstants().getParameter()->firstModule_[indx];
        md<=m_topology.dddConstants().getParameter()->lastModule_[indx];
        ++md) {
     HGCalParameters::hgtrap vol = m_topology.dddConstants().getModule(md, true, true);
     ParmVec params( numberOfParametersPerShape, 0 );
     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 (int wafer=0; wafer<m_topology.dddConstants().sectors(); ++wafer) {
         if (m_topology.dddConstants().waferInLayer(wafer,layer,true)) {
       HGCalParameters::hgtrap vol = m_topology.dddConstants().getModule(wafer, true, true);
       ParmVec params( numberOfParametersPerShape, 0 );
       params[0] = vol.dz;
       params[1] = vol.cellSize;
       params[2] = twoBysqrt3_*params[1];
       dimVector.insert( dimVector.end(), params.begin(), params.end());
     }
       }
     }
   }
   
   for (unsigned int i( 0 ); i < numberOfCells; ++i) {
     DetId detId = m_validGeomIds[i];
     int layer(0);
     if ((mode_ == HGCalGeometryMode::Hexagon) ||
     (mode_ == HGCalGeometryMode::HexagonFull)) {
       layer     = HGCalDetId(detId).layer();
     } else if (mode_ == HGCalGeometryMode::Trapezoid) {
       layer     = HGCScintillatorDetId(detId).layer();
     } else if (m_topology.isHFNose()) {
       layer     = HFNoseDetId(detId).layer();
     } else {
       layer     = HGCSiliconDetId(detId).layer();
     }
     dinsVector.emplace_back(m_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 84 of file HGCalGeometry.h.

References CaloSubdetectorGeometry::m_validIds.

84 { 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 85 of file HGCalGeometry.h.

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

85 { return m_validGeomIds; }

HGCalGeometry::m_validGeomIds

std::vector< DetId > m_validGeomIds

Definition: HGCalGeometry.h:131

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

overrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 421 of file HGCalGeometry.cc.

References TauDecayModes::dec, HGCalTopology::detId2denseGeomId(), getGeometryDetId(), DetId::HGCalHSc, m_cellVec, m_cellVec2, m_det, m_topology, and DetId::rawId().

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

                                                              {
   unsigned int cellIndex = ((m_det == DetId::HGCalHSc) ? m_cellVec2.size() :
                 m_cellVec.size());
   if (detId != DetId()) {
     DetId geomId = getGeometryDetId(detId);
     cellIndex = m_topology.detId2denseGeomId(geomId);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "indexFor " << std::hex << detId.rawId() 
                   << ":" << geomId.rawId() << std::dec 
                   << " index " << cellIndex;
 #endif
   }
   return cellIndex;
 }

void HGCalGeometry::initializeParms ( )

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 47 of file HGCalGeometry.cc.

Referenced by getValidGeomDetIds().

47 {

48 }

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

Definition at line 50 of file HGCalGeometry.cc.

References DetId::HGCalHSc, FlatTrd::localCorners(), FlatHexagon::localCorners(), and m_det.

Referenced by dbString().

                                       {
   if (m_det == DetId::HGCalHSc) {
     FlatTrd::localCorners( lc, pv, ref ) ;
   } else {
     FlatHexagon::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 61 of file HGCalGeometry.cc.

References HGCalDDDConstants::cellInLayer(), CaloSubdetectorGeometry::cornersMgr(), HGCalTopology::dddConstants(), TauDecayModes::dec, HGCalTopology::decode(), HGCalTopology::detId2denseGeomId(), HGCalTopology::encode(), getGeometryDetId(), HGCalGeometryMode::Hexagon, HGCalGeometryMode::HexagonFull, DetId::HGCalHSc, HGCalTopology::isHFNose(), m_cellVec, m_cellVec2, m_det, m_topology, m_validGeomIds, CaloSubdetectorGeometry::m_validIds, mode_, HGCalDDDConstants::numberCellsHexagon(), DetId::rawId(), HGCalGeometryMode::Trapezoid, findQualityFiles::v, HGCalTopology::valid(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

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

                                             {
 
   DetId geomId = getGeometryDetId(detId);
   int   cells (0);
   HGCalTopology::DecodedDetId id = m_topology.decode(detId);
   if ((mode_ == HGCalGeometryMode::Hexagon) || 
       (mode_ == HGCalGeometryMode::HexagonFull)) {
     cells = m_topology.dddConstants().numberCellsHexagon(id.iSec1);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "NewCell " << HGCalDetId(detId) 
                   << " GEOM " << HGCalDetId(geomId);
 #endif
   } else if (mode_ == HGCalGeometryMode::Trapezoid) {
     cells  = 1;
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "NewCell " << HGCScintillatorDetId(detId) 
                   << " GEOM " << HGCScintillatorDetId(geomId);
 #endif
   } else {
     cells  = m_topology.dddConstants().numberCellsHexagon(id.iLay,id.iSec1,
                               id.iSec2,false);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "NewCell " << HGCSiliconDetId(detId) 
                   << " GEOM " << HGCSiliconDetId(geomId);
 #endif
   }
   const uint32_t cellIndex (m_topology.detId2denseGeomId(geomId));
 
   if (m_det == DetId::HGCalHSc) {
     m_cellVec2.at( cellIndex ) = FlatTrd( cornersMgr(), f1, f2, f3, parm ) ;
   } else {
     m_cellVec.at( cellIndex ) = FlatHexagon( cornersMgr(), f1, f2, f3, parm ) ;
   }
   m_validGeomIds.at( cellIndex ) = geomId ;
   
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "Store for DetId " << std::hex 
                 << detId.rawId() << " GeomId " 
                 << geomId.rawId() << std::dec << " Index " 
                 << cellIndex << " cells " << cells;
   unsigned int nOld = m_validIds.size();
 #endif
   if ((mode_ == HGCalGeometryMode::Hexagon) || 
       (mode_ == HGCalGeometryMode::HexagonFull)) {
     for (int cell = 0; cell < cells; ++cell) {
       id.iCell1 = cell;
       DetId idc = m_topology.encode(id);
       if (m_topology.valid(idc)) {
     m_validIds.emplace_back(idc);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "Valid Id [" << cell << "] "
                       << HGCalDetId(idc);
 #endif
       }
     }
   } else if (mode_ == HGCalGeometryMode::Trapezoid) {
     DetId idc = m_topology.encode(id);
     if (m_topology.valid(idc)) {
       m_validIds.emplace_back(idc);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCalGeom") << "Valid Id [0] " 
                     << HGCScintillatorDetId(idc);
 #endif
     } else {
       edm::LogWarning("HGCalGeom") << "Check " << HGCScintillatorDetId(idc) 
                    << " from " << HGCScintillatorDetId(detId)
                    << " ERROR ???";
     }
   } else {
 #ifdef EDM_ML_DEBUG
     unsigned int cellAll(0), cellSelect(0);
 #endif
     for (int u=0; u<2*cells; ++u) {
       for (int v=0; v<2*cells; ++v) {
     if (((v-u) < cells) && (u-v) <= cells) {
       id.iCell1 = u; id.iCell2 = v;
       DetId idc = m_topology.encode(id);
 #ifdef EDM_ML_DEBUG
       ++cellAll;
 #endif
       if (m_topology.dddConstants().cellInLayer(id.iSec1,id.iSec2,u,v,
                             id.iLay,true)) {
         m_validIds.emplace_back(idc);
 #ifdef EDM_ML_DEBUG
         ++cellSelect;
         edm::LogVerbatim("HGCalGeom") << "Valid Id [" << u << ", " << v
                       << "] " << HGCSiliconDetId(idc);
 #endif
       }
     }
       }
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "HGCalGeometry keeps " << cellSelect
                   << " out of " << cellAll << " for wafer "
                   << id.iSec1 << ":" << id.iSec2 << " in "
                   << " layer " << id.iLay;
 #endif
   }
 #ifdef EDM_ML_DEBUG
   if (m_det == DetId::HGCalHSc) {
     edm::LogVerbatim("HGCalGeom") << "HGCalGeometry::newCell-> [" << cellIndex 
                   << "]" << " front:" << f1.x() << '/' 
                   << f1.y() << '/' << f1.z() << " back:" 
                   << f2.x() << '/' << f2.y() << '/' << f2.z() 
                   << " eta|phi " 
                   << m_cellVec2[cellIndex].etaPos() << ":"
                   << m_cellVec2[cellIndex].phiPos();
   } else {
     edm::LogVerbatim("HGCalGeom") << "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();
   }
   unsigned int nNew = m_validIds.size();
   if ((mode_ == HGCalGeometryMode::Hexagon) || 
       (mode_ == HGCalGeometryMode::HexagonFull)) {
     edm::LogVerbatim("HGCalGeom") << "ID: " << HGCalDetId(detId) 
                   << " with valid DetId from " << nOld 
                   << " to " << nNew;
   } else if (mode_ == HGCalGeometryMode::Trapezoid) {
     edm::LogVerbatim("HGCalGeom") << "ID: " << HGCScintillatorDetId(detId) 
                   << " with valid DetId from " << nOld 
                   << " to " << nNew;
   } else if (m_topology.isHFNose()) {
     edm::LogVerbatim("HGCalGeom") << "ID: " << HFNoseDetId(detId) 
                   << " with valid DetId from " << nOld 
                   << " to " << nNew;
   } else {
     edm::LogVerbatim("HGCalGeom") << "ID: " << HGCSiliconDetId(detId) 
                   << " with valid DetId from " << nOld 
                   << " to " << nNew;
   }
   edm::LogVerbatim("HGCalGeom") << "Cell[" << cellIndex << "] " << std::hex 
                 << geomId.rawId() << ":"  
                 << m_validGeomIds[cellIndex].rawId() 
                 << std::dec;
 #endif
 }

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

overridevirtual

is this detid present in the geometry?

Reimplemented from CaloSubdetectorGeometry.

Definition at line 217 of file HGCalGeometry.cc.

References HGCalTopology::detId2denseGeomId(), getGeometryDetId(), getGeometryRawPtr(), and m_topology.

Referenced by dbString().

                                                     {
   if (detId == DetId()) return false;
   DetId geomId = getGeometryDetId(detId);
   const uint32_t index (m_topology.detId2denseGeomId(geomId));
   return (nullptr != getGeometryRawPtr(index)) ;
 }

static std::string HGCalGeometry::producerTag ( )

inlinestatic

Definition at line 101 of file HGCalGeometry.h.

References cellElement(), and AlCaHLTBitMon_QueryRunRegistry::string.

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

101 { return "HGCal" ; }

unsigned int HGCalGeometry::sizeForDenseIndex ( ) const

protected

Definition at line 436 of file HGCalGeometry.cc.

References m_topology, and HGCalTopology::totalGeomModules().

Referenced by topology().

                                                     {
   return m_topology.totalGeomModules();
 }

void HGCalGeometry::sortDetIds ( void )

Definition at line 549 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 104 of file HGCalGeometry.h.

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

Referenced by hgc_digi_utils::addCellMetadata(), HGCalRecHitValidation::analyze(), HGCalDigiValidation::analyze(), PFRecHitQTestThresholdInThicknessNormalizedMIPs::beginEvent(), CaloTruthAccumulator::beginLuminosityBlock(), HGCDigitizer::checkPosition(), dbString(), HGCalTriggerGeometryBase::eeTopology(), HGCalTriggerGeometryBase::fhTopology(), CaloParticleDebugger::fillSimHits(), HGCalTriggerGeometryBase::hscTopology(), HGCalTriggerGeometryBase::hsiTopology(), CaloGeometryBuilder::produceAligned(), HGCalRecHitWorkerSimple::set(), and HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::setGeometry().

104 {return m_topology;}

HGCalGeometry::m_topology

const HGCalTopology & m_topology

Definition: HGCalGeometry.h:128

Member Data Documentation

unsigned int HGCalGeometry::k_NumberOfParametersPerHex = 3

static

Definition at line 45 of file HGCalGeometry.h.

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

unsigned int HGCalGeometry::k_NumberOfParametersPerShape = 3

static

Definition at line 46 of file HGCalGeometry.h.

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

unsigned int HGCalGeometry::k_NumberOfParametersPerTrd = 12

static

Definition at line 44 of file HGCalGeometry.h.

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

unsigned int HGCalGeometry::k_NumberOfShapes = 100

static

Definition at line 47 of file HGCalGeometry.h.

Referenced by HGCalGeometryLoader::build(), getSummary(), and CaloGeometryDBEP< T, U >::produceAligned().

unsigned int HGCalGeometry::k_NumberOfShapesTrd = 1000

static

Definition at line 48 of file HGCalGeometry.h.

Referenced by HGCalGeometryLoader::build().

CellVec HGCalGeometry::m_cellVec

private

Definition at line 129 of file HGCalGeometry.h.

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

CellVec2 HGCalGeometry::m_cellVec2

private

Definition at line 130 of file HGCalGeometry.h.

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

DetId::Detector HGCalGeometry::m_det

private

Definition at line 133 of file HGCalGeometry.h.

Referenced by cellElement(), cellGeomPtr(), get8Corners(), getClosestCell(), getClosestCellIndex(), getCorners(), getGeometryRawPtr(), getSummary(), HGCalGeometry(), indexFor(), localCorners(), and newCell().

ForwardSubdetector HGCalGeometry::m_subdet

private

Definition at line 134 of file HGCalGeometry.h.

Referenced by cellElement().

const HGCalTopology& HGCalGeometry::m_topology

private

Definition at line 128 of file HGCalGeometry.h.

Referenced by get8Corners(), getClosestCell(), getCorners(), getGeometry(), getGeometryDetId(), getPosition(), getSummary(), HGCalGeometry(), indexFor(), newCell(), present(), sizeForDenseIndex(), and topology().

std::vector<DetId> HGCalGeometry::m_validGeomIds

private

Definition at line 131 of file HGCalGeometry.h.

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

HGCalGeometryMode::GeometryMode HGCalGeometry::mode_

private

Definition at line 132 of file HGCalGeometry.h.

Referenced by get8Corners(), getClosestCell(), getCorners(), getGeometryDetId(), getPosition(), getSummary(), and newCell().

const double HGCalGeometry::twoBysqrt3_

private

Definition at line 135 of file HGCalGeometry.h.

Referenced by getSummary().

Public Types

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation