#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

using	CellMayOwnPtr = CaloCellGeometryMayOwnPtr

using	CellPtr = CaloCellGeometryPtr

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

double	getArea (const DetId &detid) const
	Returns area of a cell. More...

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

DetId	getClosestCellHex (const GlobalPoint &r, bool extend) const

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

CaloCellGeometryMayOwnPtr	getGeometry (const DetId &id) const override
	Get the cell geometry of a given detector id. Should return false if not found. More...

CornersVec	getNewCorners (const DetId &id, bool debug=false) const

GlobalPoint	getPosition (const DetId &id, bool cog, bool debug) const

GlobalPoint	getPosition (const DetId &id, bool debug=false) 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

GlobalPoint	getWaferPosition (const DetId &id) const

	HGCalGeometry (const HGCalTopology &topology)

void	initializeParms () override

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

DetId	neighborZ (const DetId &idin, const GlobalVector &p) const

DetId	neighborZ (const DetId &idin, const MagneticField *bField, int charge, const GlobalVector &momentum) const

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=default

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 bool	valid (const DetId &id) 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 constexpr unsigned int	k_NumberOfParametersPerHex = 3

static constexpr unsigned int	k_NumberOfParametersPerShape = 3

static constexpr unsigned int	k_NumberOfParametersPerTrd = 12

static constexpr unsigned int	k_NumberOfShapes = 100

static constexpr unsigned int	k_NumberOfShapesTrd = 1000

Protected Member Functions
void	addValidID (const DetId &id)

CaloCellGeometryPtr	cellGeomPtr (uint32_t index) const override

unsigned int	getClosestCellIndex (const GlobalPoint &r) const

CaloCellGeometryPtr	getGeometryRawPtr (uint32_t index) const override

unsigned int	indexFor (const DetId &id) const override

virtual unsigned int	sizeForDenseIndex (const DetId &id) const

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
CaloCellGeometryMayOwnPtr	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

const double	twoBysqrt3_

Static Private Attributes
static constexpr double	k_fac1 = 0.5

static constexpr double	k_fac2 = 1.0 / 3.0

static constexpr double	k_half = 0.5

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

Member Typedef Documentation

◆ AlignedRecord

typedef HGCalGeometryRecord HGCalGeometry::AlignedRecord

Definition at line 40 of file HGCalGeometry.h.

◆ CCGFloat

typedef CaloCellGeometry::CCGFloat HGCalGeometry::CCGFloat

Definition at line 33 of file HGCalGeometry.h.

◆ CellVec

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

Definition at line 31 of file HGCalGeometry.h.

◆ CellVec2

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

Definition at line 32 of file HGCalGeometry.h.

◆ CornersVec

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

Definition at line 38 of file HGCalGeometry.h.

◆ DetIdSet

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

Definition at line 37 of file HGCalGeometry.h.

◆ PGeometryRecord

typedef PHGCalRcd HGCalGeometry::PGeometryRecord

Definition at line 41 of file HGCalGeometry.h.

◆ Pt3D

typedef CaloCellGeometry::Pt3D HGCalGeometry::Pt3D

Definition at line 34 of file HGCalGeometry.h.

◆ Pt3DVec

typedef CaloCellGeometry::Pt3DVec HGCalGeometry::Pt3DVec

Definition at line 35 of file HGCalGeometry.h.

Constructor & Destructor Documentation

◆ HGCalGeometry()

HGCalGeometry::HGCalGeometry ( const HGCalTopology & topology )

Definition at line 28 of file HGCalGeometry.cc.

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

     : m_topology(topology_),
       m_validGeomIds(topology_.totalGeomModules()),
       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::~HGCalGeometry ( )

overridedefault

Member Function Documentation

◆ addValidID()

void HGCalGeometry::addValidID ( const DetId & id )

protected

Definition at line 725 of file HGCalGeometry.cc.

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

◆ cellElement()

std::string HGCalGeometry::cellElement ( ) const

Definition at line 640 of file HGCalGeometry.cc.

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

                                            {
   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";
 }

◆ cellGeomPtr() [1/2]

CaloCellGeometryPtr HGCalGeometry::cellGeomPtr ( uint32_t index ) const

overrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 681 of file HGCalGeometry.cc.

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

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

                                                                    {
   if ((index >= m_cellVec.size() && m_det != DetId::HGCalHSc) ||
       (index >= m_cellVec2.size() && m_det == DetId::HGCalHSc) || (m_validGeomIds[index].rawId() == 0))
     return CaloCellGeometryPtr();
   if (m_det == DetId::HGCalHSc) {
     auto cell = &m_cellVec2[index];
     if (nullptr == cell->param())
       return CaloCellGeometryPtr();
     return CaloCellGeometryPtr(cell);
   } else {
     auto cell = &m_cellVec[index];
     if (nullptr == cell->param())
       return CaloCellGeometryPtr(nullptr);
     return CaloCellGeometryPtr(cell);
   }
 }

◆ cellGeomPtr() [2/2]

CaloCellGeometryMayOwnPtr HGCalGeometry::cellGeomPtr	(	uint32_t	index,
		const GlobalPoint &	p
	)		const

private

Definition at line 698 of file HGCalGeometry.cc.

References cellGeomPtr(), DetId::HGCalHSc, m_cellVec, m_cellVec2, m_det, m_validGeomIds, eostools::move(), and nano_mu_digi_cff::rawId.

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

◆ dbString()

static std::string HGCalGeometry::dbString ( )

inlinestatic

Definition at line 49 of file HGCalGeometry.h.

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

49 { return "PHGCalRcd"; }

◆ fillNamedParams()

void HGCalGeometry::fillNamedParams ( DDFilteredView fv )

virtual

Definition at line 45 of file HGCalGeometry.cc.

45 {}

◆ get8Corners()

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

Definition at line 346 of file HGCalGeometry.cc.

References HGCalDDDConstants::cellSizeTrap(), cms::cuda::co, funct::cos(), HGCalTopology::dddConstants(), debugLocate, HGCalTopology::decode(), ALCARECOPPSCalTrackBasedSel_cff::detid, l1ctLayer1_cff::dr, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, getPosition(), DetId::HGCalHSc, mps_fire::i, indexFor(), FlatTrd::k_Cell, FlatTrd::k_dZ, FlatHexagon::k_dZ, k_fac1, k_fac2, k_half, FlatHexagon::k_r, FlatHexagon::k_R, HGCalDDDConstants::localToGlobal8(), HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), m_cellVec, m_cellVec2, m_det, m_topology, SiStripPI::min, FlatTrd::ncorner_, alignCSCRings::r, findQualityFiles::rr, funct::sin(), findQualityFiles::v, HGCalTopology::waferHexagon6(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, and geometryCSVtoXML::zz.

                                                                            {
   unsigned int ncorner = FlatTrd::ncorner_;
   HGCalGeometry::CornersVec co(ncorner, GlobalPoint(0, 0, 0));
   unsigned int cellIndex = indexFor(detid);
   HGCalTopology::DecodedDetId id = m_topology.decode(detid);
   if (cellIndex < m_cellVec2.size() && m_det == DetId::HGCalHSc) {
     GlobalPoint v = getPosition(detid, false);
     int type = std::min(id.iType, 1);
     std::pair<double, double> rr = m_topology.dddConstants().cellSizeTrap(type, id.iSec1);
     float dr = k_half * (rr.second - rr.first);
     float dfi = m_cellVec2[cellIndex].param()[FlatTrd::k_Cell];
     float dz = id.zSide * m_cellVec2[cellIndex].param()[FlatTrd::k_dZ];
     float r = v.perp();
     float fi = v.phi();
     static const int signr[] = {1, 1, -1, -1, 1, 1, -1, -1};
     static const int signf[] = {-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) {
       co[i] = GlobalPoint((r + signr[i] * dr) * cos(fi + signf[i] * dfi),
                           (r + signr[i] * dr) * sin(fi + signf[i] * dfi),
                           (v.z() + signz[i] * dz));
     }
   } else if (cellIndex < m_cellVec.size() && m_det != DetId::HGCalHSc) {
     std::pair<float, float> xy;
     float dx(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};
     if (m_topology.waferHexagon6()) {
       xy = m_topology.dddConstants().locateCellHex(id.iCell1, id.iSec1, true);
       dx = m_cellVec[cellIndex].param()[FlatHexagon::k_r];
       float dz = m_cellVec[cellIndex].param()[FlatHexagon::k_dZ];
       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_cellVec[cellIndex].getPosition(lcoord);
       }
     } else {
       xy = m_topology.dddConstants().locateCell(
           id.zSide, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2, true, false, true, false, debugLocate);
       dx = k_fac2 * m_cellVec[cellIndex].param()[FlatHexagon::k_r];
       float dy = k_fac1 * m_cellVec[cellIndex].param()[FlatHexagon::k_R];
       float dz = -id.zSide * m_cellVec[cellIndex].param()[FlatHexagon::k_dZ];
       float zz = m_topology.dddConstants().waferZ(id.iLay, true);
       for (unsigned int i = 0; i < ncorner; ++i) {
         auto xyglob = m_topology.dddConstants().localToGlobal8(
             id.zSide, id.iLay, id.iSec1, id.iSec2, (xy.first + signx[i] * dx), (xy.second + signy[i] * dy), true, false);
         double xx = id.zSide * xyglob.first;
         co[i] = GlobalPoint(xx, xyglob.second, id.zSide * (zz + signz[i] * dz));
       }
     }
   }
   return co;
 }

◆ getArea()

double HGCalGeometry::getArea ( const DetId & detid ) const

Returns area of a cell.

Definition at line 275 of file HGCalGeometry.cc.

References funct::abs(), custom_jme_cff::area, ALCARECOPPSCalTrackBasedSel_cff::detid, getNewCorners(), mps_fire::i, dqmiolumiharvest::j, create_idmaps::n, and y.

                                                       {
   HGCalGeometry::CornersVec corners = getNewCorners(detid);
   double area(0);
   if (corners.size() > 1) {
     int n = corners.size() - 1;
     int j = n - 1;
     for (int i = 0; i < n; ++i) {
       area += ((corners[j].x() + corners[i].x()) * (corners[i].y() - corners[j].y()));
       j = i;
     }
   }
   return std::abs(0.5 * area);
 }

◆ getCells()

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

                                                                                    {
   HGCalGeometry::DetIdSet dss;
   return dss;
 }

◆ getClosestCell()

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

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 533 of file HGCalGeometry.cc.

References HGCalDDDConstants::assignCell(), HGCalDDDConstants::assignCellHex(), HGCalDDDConstants::assignCellTrap(), HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalTopology::DecodedDetId::det, HGCalTopology::detector(), HGCalTopology::encode(), getClosestCellIndex(), HGCalDDDConstants::getLayer(), DetId::HGCalHSc, DTRecHitClients_cfi::local, m_cellVec, m_cellVec2, m_det, m_topology, m_validGeomIds, alignCSCRings::r, HGCalTopology::tileTrapezoid(), HGCalTopology::waferHexagon6(), HGCalDDDConstants::waferTypeT(), and ecaldqm::zside().

Referenced by neighborZ().

                                                               {
   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]);
     if (id.det == 0)
       id.det = static_cast<int>(m_topology.detector());
     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 (m_topology.waferHexagon6()) {
       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 (m_topology.tileTrapezoid()) {
       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);
       int zside = (r.z() > 0) ? 1 : -1;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCalGeom") << "ZZ " << r.z() << ":" << zside << " Layer " << id.iLay << " Global " << r
                                     << "  Local " << local;
 #endif
       const auto& kxy =
           m_topology.dddConstants().assignCellHex(local.x(), local.y(), zside, id.iLay, true, false, 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.det << ":" << 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();
 }

◆ getClosestCellHex()

DetId HGCalGeometry::getClosestCellHex	(	const GlobalPoint &	r,
		bool	extend
	)		const

Definition at line 591 of file HGCalGeometry.cc.

References HGCalDDDConstants::assignCellHex(), HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalTopology::DecodedDetId::det, HGCalTopology::detector(), HGCalTopology::encode(), getClosestCellIndex(), HGCalDDDConstants::getLayer(), DetId::HGCalHSc, DTRecHitClients_cfi::local, m_cellVec, m_det, m_topology, m_validGeomIds, alignCSCRings::r, HGCalTopology::waferHexagon8(), and ecaldqm::zside().

                                                                               {
   unsigned int cellIndex = getClosestCellIndex(r);
   if (cellIndex < m_cellVec.size() && m_det != DetId::HGCalHSc) {
     HGCalTopology::DecodedDetId id = m_topology.decode(m_validGeomIds[cellIndex]);
     if (id.det == 0)
       id.det = static_cast<int>(m_topology.detector());
     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 (m_topology.waferHexagon8()) {
       id.iLay = m_topology.dddConstants().getLayer(r.z(), true);
       int zside = (r.z() > 0) ? 1 : -1;
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCalGeom") << "ZZ " << r.z() << ":" << zside << " Layer " << id.iLay << " Global " << r
                                     << "  Local " << local;
 #endif
       const auto& kxy =
           m_topology.dddConstants().assignCellHex(local.x(), local.y(), zside, id.iLay, true, extend, 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.det << ":" << 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();
 }

◆ getClosestCellIndex() [1/2]

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

protected

Definition at line 729 of file HGCalGeometry.cc.

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

Referenced by getClosestCell(), and getClosestCellHex().

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

◆ getClosestCellIndex() [2/2]

template<class T >

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

private

Definition at line 734 of file HGCalGeometry.cc.

References funct::abs(), PVValHelper::dz, getPosition(), isotrackApplyRegressor::k, M_PI, alignCSCRings::r, and 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;
 }

◆ getCorners()

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

Returns the corner points of this cell's volume.

Definition at line 289 of file HGCalGeometry.cc.

References HGCalDDDConstants::cellSizeTrap(), cms::cuda::co, funct::cos(), HGCalTopology::dddConstants(), debugLocate, HGCalTopology::decode(), ALCARECOPPSCalTrackBasedSel_cff::detid, l1ctLayer1_cff::dr, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, getPosition(), DetId::HGCalHSc, mps_fire::i, indexFor(), FlatTrd::k_Cell, FlatTrd::k_dZ, FlatHexagon::k_dZ, k_fac1, k_fac2, k_half, FlatHexagon::k_r, FlatHexagon::k_R, HGCalDDDConstants::localToGlobal8(), HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), m_cellVec, m_cellVec2, m_det, m_topology, SiStripPI::min, FlatHexagon::ncorner_, FlatTrd::ncorner_, alignCSCRings::r, findQualityFiles::rr, funct::sin(), findQualityFiles::v, HGCalTopology::waferHexagon6(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, and geometryCSVtoXML::zz.

                                                                           {
   unsigned int ncorner = ((m_det == DetId::HGCalHSc) ? FlatTrd::ncorner_ : FlatHexagon::ncorner_);
   HGCalGeometry::CornersVec co(ncorner, GlobalPoint(0, 0, 0));
   unsigned int cellIndex = indexFor(detid);
   HGCalTopology::DecodedDetId id = m_topology.decode(detid);
   if (cellIndex < m_cellVec2.size() && m_det == DetId::HGCalHSc) {
     GlobalPoint v = getPosition(detid, false);
     int type = std::min(id.iType, 1);
     std::pair<double, double> rr = m_topology.dddConstants().cellSizeTrap(type, id.iSec1);
     float dr = k_half * (rr.second - rr.first);
     float dfi = m_cellVec2[cellIndex].param()[FlatTrd::k_Cell];
     float dz = id.zSide * m_cellVec2[cellIndex].param()[FlatTrd::k_dZ];
     float r = v.perp();
     float fi = v.phi();
     static const int signr[] = {1, 1, -1, -1, 1, 1, -1, -1};
     static const int signf[] = {-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) {
       co[i] = GlobalPoint((r + signr[i] * dr) * cos(fi + signf[i] * dfi),
                           (r + signr[i] * dr) * sin(fi + signf[i] * dfi),
                           (v.z() + signz[i] * dz));
     }
   } else if (cellIndex < m_cellVec.size() && m_det != DetId::HGCalHSc) {
     std::pair<float, float> xy;
     if (m_topology.waferHexagon6()) {
       xy = m_topology.dddConstants().locateCellHex(id.iCell1, id.iSec1, true);
       float dx = m_cellVec[cellIndex].param()[FlatHexagon::k_r];
       float dy = k_half * m_cellVec[cellIndex].param()[FlatHexagon::k_R];
       float dz = m_cellVec[cellIndex].param()[FlatHexagon::k_dZ];
       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);
       }
     } else {
       xy = m_topology.dddConstants().locateCell(
           id.zSide, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2, true, false, true, false, debugLocate);
       float zz = m_topology.dddConstants().waferZ(id.iLay, true);
       float dx = k_fac2 * m_cellVec[cellIndex].param()[FlatHexagon::k_r];
       float dy = k_fac1 * m_cellVec[cellIndex].param()[FlatHexagon::k_R];
       float dz = -id.zSide * m_cellVec[cellIndex].param()[FlatHexagon::k_dZ];
       static const int signx[] = {1, -1, -2, -1, 1, 2, 1, -1, -2, -1, 1, 2};
       static const int signy[] = {1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0};
       static const int signz[] = {-1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1, 1};
       for (unsigned int i = 0; i < ncorner; ++i) {
         auto xyglob = m_topology.dddConstants().localToGlobal8(
             id.zSide, id.iLay, id.iSec1, id.iSec2, (xy.first + signx[i] * dx), (xy.second + signy[i] * dy), true, false);
         double xx = id.zSide * xyglob.first;
         co[i] = GlobalPoint(xx, xyglob.second, id.zSide * (zz + signz[i] * dz));
       }
     }
   }
   return co;
 }

◆ getGeometry()

CaloCellGeometryMayOwnPtr 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 187 of file HGCalGeometry.cc.

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

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

◆ getGeometryDetId()

DetId HGCalGeometry::getGeometryDetId ( DetId detId ) const

private

Definition at line 882 of file HGCalGeometry.cc.

References hcalRecHitTable_cff::detId, HGCalDetId::geometryCell(), HFNoseDetId::geometryCell(), HGCScintillatorDetId::geometryCell(), HGCSiliconDetId::geometryCell(), HGCalTopology::isHFNose(), m_topology, HGCalTopology::tileTrapezoid(), and HGCalTopology::waferHexagon6().

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

                                                        {
   DetId geomId;
   if (m_topology.waferHexagon6()) {
     geomId = static_cast<DetId>(HGCalDetId(detId).geometryCell());
   } else if (m_topology.tileTrapezoid()) {
     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;
 }

◆ getGeometryRawPtr()

CaloCellGeometryPtr HGCalGeometry::getGeometryRawPtr ( uint32_t index ) const

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 666 of file HGCalGeometry.cc.

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

Referenced by present().

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

◆ getNewCorners()

HGCalGeometry::CornersVec HGCalGeometry::getNewCorners	(	const DetId &	id,
		bool	debug = `false`
	)		const

Definition at line 400 of file HGCalGeometry.cc.

References HGCalDDDConstants::cellSizeTrap(), cms::cuda::co, funct::cos(), HGCalTopology::dddConstants(), debug, HGCalTopology::decode(), ALCARECOPPSCalTrackBasedSel_cff::detid, l1ctLayer1_cff::dr, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, getPosition(), DetId::HGCalHSc, mps_fire::i, indexFor(), FlatTrd::k_Cell, FlatTrd::k_dZ, FlatHexagon::k_dZ, k_fac1, k_fac2, k_half, FlatHexagon::k_r, FlatHexagon::k_R, HGCalDDDConstants::localToGlobal8(), HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), m_cellVec, m_cellVec2, m_det, m_topology, SiStripPI::min, alignCSCRings::r, findQualityFiles::rr, funct::sin(), findQualityFiles::v, HGCalTopology::waferHexagon6(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, and geometryCSVtoXML::zz.

Referenced by getArea().

                                                                                          {
   unsigned int ncorner = (m_det == DetId::HGCalHSc) ? 5 : 7;
   HGCalGeometry::CornersVec co(ncorner, GlobalPoint(0, 0, 0));
   unsigned int cellIndex = indexFor(detid);
   HGCalTopology::DecodedDetId id = m_topology.decode(detid);
   if (debug)
     edm::LogVerbatim("HGCalGeom") << "NewCorners for Layer " << id.iLay << " Wafer " << id.iSec1 << ":" << id.iSec2
                                   << " Cell " << id.iCell1 << ":" << id.iCell2;
   if (cellIndex < m_cellVec2.size() && m_det == DetId::HGCalHSc) {
     GlobalPoint v = getPosition(detid, false);
     int type = std::min(id.iType, 1);
     std::pair<double, double> rr = m_topology.dddConstants().cellSizeTrap(type, id.iSec1);
     float dr = k_half * (rr.second - rr.first);
     float dfi = m_cellVec2[cellIndex].param()[FlatTrd::k_Cell];
     float dz = -id.zSide * m_cellVec2[cellIndex].param()[FlatTrd::k_dZ];
     float r = v.perp();
     float fi = v.phi();
     static const int signr[] = {1, 1, -1, -1};
     static const int signf[] = {-1, 1, 1, -1};
     for (unsigned int i = 0; i < ncorner - 1; ++i) {
       co[i] = GlobalPoint(
           (r + signr[i] * dr) * cos(fi + signf[i] * dfi), (r + signr[i] * dr) * sin(fi + signf[i] * dfi), (v.z() + dz));
     }
     // Used to pass downstream the thickness of this cell
     co[ncorner - 1] = GlobalPoint(0, 0, -2 * dz);
   } else if (cellIndex < m_cellVec.size() && m_det != DetId::HGCalHSc) {
     std::pair<float, float> xy;
     float dx = k_fac2 * m_cellVec[cellIndex].param()[FlatHexagon::k_r];
     float dy = k_fac1 * m_cellVec[cellIndex].param()[FlatHexagon::k_R];
     float dz = -id.zSide * m_cellVec[cellIndex].param()[FlatHexagon::k_dZ];
     static const int signx[] = {1, -1, -2, -1, 1, 2};
     static const int signy[] = {1, 1, 0, -1, -1, 0};
 #ifdef EDM_ML_DEBUG
     if (debug)
       edm::LogVerbatim("HGCalGeom") << "kfac " << k_fac1 << ":" << k_fac2 << " dx:dy:dz " << dx << ":" << dy << ":"
                                     << dz;
 #endif
     if (m_topology.waferHexagon6()) {
       xy = m_topology.dddConstants().locateCellHex(id.iCell1, id.iSec1, true);
       for (unsigned int i = 0; i < ncorner - 1; ++i) {
         const HepGeom::Point3D<float> lcoord(xy.first + signx[i] * dx, xy.second + signy[i] * dy, dz);
         co[i] = m_cellVec[cellIndex].getPosition(lcoord);
       }
     } else {
       xy = m_topology.dddConstants().locateCell(
           id.zSide, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2, true, false, true, false, debug);
       float zz = m_topology.dddConstants().waferZ(id.iLay, true);
       for (unsigned int i = 0; i < ncorner; ++i) {
         double xloc = xy.first + signx[i] * dx;
         double yloc = xy.second + signy[i] * dy;
 #ifdef EDM_ML_DEBUG
         if (debug)
           edm::LogVerbatim("HGCalGeom") << "Corner " << i << " x " << xy.first << ":" << xloc << " y " << xy.second
                                         << ":" << yloc << " z " << zz << ":" << id.zSide * (zz + dz);
 #endif
         auto xyglob =
             m_topology.dddConstants().localToGlobal8(id.zSide, id.iLay, id.iSec1, id.iSec2, xloc, yloc, true, debug);
         double xx = id.zSide * xyglob.first;
         co[i] = GlobalPoint(xx, xyglob.second, id.zSide * (zz + dz));
       }
     }
     // Used to pass downstream the thickness of this cell
     co[ncorner - 1] = GlobalPoint(0, 0, -2 * dz);
   }
   return co;
 }

◆ getPosition() [1/2]

GlobalPoint HGCalGeometry::getPosition	(	const DetId &	id,
		bool	cog,
		bool	debug
	)		const

Definition at line 208 of file HGCalGeometry.cc.

References HGCalTopology::dddConstants(), debug, HGCalTopology::decode(), ALCARECOPPSCalTrackBasedSel_cff::detid, DetId::Forward, RooUtil::FileUtil::glob(), indexFor(), HGCalDDDConstants::locateCell(), HGCalDDDConstants::locateCellHex(), m_cellVec, m_cellVec2, m_topology, reco_skim_cfg_mod::maxSize, HGCalTopology::tileTrapezoid(), HGCalTopology::waferHexagon6(), HGCalDDDConstants::waferZ(), geometryCSVtoXML::xx, geometryCSVtoXML::xy, and geometryCSVtoXML::zz.

Referenced by HGCDigitizer::checkPosition(), HGCalRadiationMap::computeRadius(), get8Corners(), getClosestCellIndex(), getCorners(), getGeometry(), HGCalTriggerGeometryV9Imp2::getModulePosition(), HGCalTriggerGeometryV9Imp3::getModulePosition(), getNewCorners(), hgcal::RecHitTools::getPosition(), getPosition(), HGCalTriggerGeometryV9Imp2::getTriggerCellPosition(), HGCalTriggerGeometryV9Imp3::getTriggerCellPosition(), HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::makeRecHit(), and neighborZ().

                                                                                      {
   unsigned int cellIndex = indexFor(detid);
   GlobalPoint glob;
   unsigned int maxSize = (m_topology.tileTrapezoid() ? m_cellVec2.size() : m_cellVec.size());
   if (cellIndex < maxSize) {
     HGCalTopology::DecodedDetId id = m_topology.decode(detid);
     std::pair<float, float> xy;
     if (m_topology.waferHexagon6()) {
       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);
       if (debug)
         edm::LogVerbatim("HGCalGeom") << "getPosition:: index " << cellIndex << " Local " << lcoord.x() << ":"
                                       << lcoord.y() << " ID " << id.iCell1 << ":" << id.iSec1 << " Global " << glob;
     } else if (m_topology.tileTrapezoid()) {
       const HepGeom::Point3D<float> lcoord(0, 0, 0);
       glob = m_cellVec2[cellIndex].getPosition(lcoord);
       if (debug)
         edm::LogVerbatim("HGCalGeom") << "getPositionTrap:: index " << cellIndex << " Local " << lcoord.x() << ":"
                                       << lcoord.y() << " ID " << id.iLay << ":" << id.iSec1 << ":" << id.iCell1
                                       << " Global " << glob;
     } else {
       if (debug) {
         if (detid.det() == DetId::Forward)
           edm::LogVerbatim("HGCalGeom") << "getPosition for " << HFNoseDetId(detid) << " Layer " << id.iLay << " Wafer "
                                         << id.iSec1 << ":" << id.iSec2 << " Cell " << id.iCell1 << ":" << id.iCell2;
         else
           edm::LogVerbatim("HGCalGeom") << "getPosition for " << HGCSiliconDetId(detid) << " Layer " << id.iLay
                                         << " Wafer " << id.iSec1 << ":" << id.iSec2 << " Cell " << id.iCell1 << ":"
                                         << id.iCell2;
       }
       xy = m_topology.dddConstants().locateCell(
           id.zSide, id.iLay, id.iSec1, id.iSec2, id.iCell1, id.iCell2, true, true, false, cog, debug);
       double xx = id.zSide * xy.first;
       double zz = id.zSide * m_topology.dddConstants().waferZ(id.iLay, true);
       glob = GlobalPoint(xx, xy.second, zz);
       if (debug)
         edm::LogVerbatim("HGCalGeom") << "getPositionWafer:: index " << cellIndex << " Local " << xy.first << ":"
                                       << xy.second << " ID " << id.iLay << ":" << id.iSec1 << ":" << id.iSec2 << ":"
                                       << id.iCell1 << ":" << id.iCell2 << " Global " << glob;
     }
   } else {
     edm::LogVerbatim("HGCalGeom") << "Cannot recognize " << std::hex << detid.rawId() << " cellIndex " << cellIndex
                                   << ":" << maxSize << " Type " << m_topology.tileTrapezoid();
   }
   return glob;
 }

◆ getPosition() [2/2]

GlobalPoint HGCalGeometry::getPosition	(	const DetId &	id,
		bool	debug = `false`
	)		const

Definition at line 204 of file HGCalGeometry.cc.

References debug, ALCARECOPPSCalTrackBasedSel_cff::detid, and getPosition().

                                                                            {
   return getPosition(detid, false, debug);
 }

◆ getSummary()

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

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 780 of file HGCalGeometry.cc.

References HGCalParameters::hgtrap::alpha, HGCalParameters::hgtrap::bl, cellGeomPtr(), HGCalParameters::hgtrap::cellSize, HGCalTopology::dddConstants(), hcalRecHitTable_cff::detId, HGCalTopology::detId2denseGeomId(), HGCalParameters::hgtrap::dz, HGCalParameters::firstModule_, HGCalDDDConstants::getModule(), HGCalDDDConstants::getParameter(), HGCalDDDConstants::getTrForm(), HGCalDDDConstants::getTrFormN(), runTauDisplay::gp, HGCalParameters::hgtrap::h, DetId::HGCalHSc, mps_fire::i, HGCalTopology::isHFNose(), FlatTrd::k_Alp1, FlatTrd::k_Alp2, FlatTrd::k_Cell, FlatTrd::k_dX1, FlatTrd::k_dX2, FlatTrd::k_dX3, FlatTrd::k_dX4, FlatTrd::k_dY1, FlatTrd::k_dY2, FlatTrd::k_dZ, FlatHexagon::k_dZ, k_NumberOfParametersPerHex, k_NumberOfParametersPerTrd, k_NumberOfShapes, FlatTrd::k_Phi, FlatHexagon::k_r, FlatHexagon::k_R, FlatTrd::k_Theta, HGCalParameters::lastModule_, HGCalParameters::hgtrform::lay, nano_mu_digi_cff::layer, HGCalDetId::layer(), HFNoseDetId::layer(), HGCScintillatorDetId::layer(), HGCSiliconDetId::layer(), HGCalDDDConstants::layerIndex(), m_det, m_topology, m_validGeomIds, CaloSubdetectorGeometry::numberOfParametersPerShape(), CaloSubdetectorGeometry::numberOfShapes(), CaloSubdetectorGeometry::numberOfTransformParms(), submitPVValidationJobs::params, findQualityFiles::rr, HGCalDDDConstants::sectors(), HGCalTopology::tileTrapezoid(), HGCalParameters::hgtrap::tl, HGCalTopology::totalGeomModules(), groupFilesInBlocks::tt, twoBysqrt3_, HGCalTopology::waferHexagon6(), and HGCalDDDConstants::waferInLayer().

                                                                               {
   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 (m_topology.waferHexagon6()) {
       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[FlatHexagon::k_dZ] = vol.dz;
           params[FlatHexagon::k_r] = vol.cellSize;
           params[FlatHexagon::k_R] = twoBysqrt3_ * params[FlatHexagon::k_r];
           dimVector.insert(dimVector.end(), params.begin(), params.end());
         }
       }
     } else if (m_topology.tileTrapezoid()) {
       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[FlatTrd::k_dZ] = vol.dz;
         params[FlatTrd::k_Theta] = params[FlatTrd::k_Phi] = 0;
         params[FlatTrd::k_dY1] = params[FlatTrd::k_dY2] = vol.h;
         params[FlatTrd::k_dX1] = params[FlatTrd::k_dX3] = vol.bl;
         params[FlatTrd::k_dX2] = params[FlatTrd::k_dX4] = vol.tl;
         params[FlatTrd::k_Alp1] = params[FlatTrd::k_Alp2] = vol.alpha;
         params[FlatTrd::k_Cell] = 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[FlatHexagon::k_dZ] = vol.dz;
           params[FlatHexagon::k_r] = vol.cellSize;
           params[FlatHexagon::k_R] = twoBysqrt3_ * params[FlatHexagon::k_r];
           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 (m_topology.waferHexagon6()) {
       layer = HGCalDetId(detId).layer();
     } else if (m_topology.tileTrapezoid()) {
       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());
       }
     }
   }
 }

◆ getValidDetIds()

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

References CaloSubdetectorGeometry::m_validIds.

                                                                                                               {
     return m_validIds;
   }

◆ getValidGeomDetIds()

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

inline

Definition at line 93 of file HGCalGeometry.h.

References m_validGeomIds.

93 { return m_validGeomIds; }

HGCalGeometry::m_validGeomIds

std::vector< DetId > m_validGeomIds

Definition: HGCalGeometry.h:142

◆ getWaferPosition()

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

Definition at line 256 of file HGCalGeometry.cc.

References TauDecayModes::dec, ALCARECOPPSCalTrackBasedSel_cff::detid, RooUtil::FileUtil::glob(), indexFor(), m_cellVec, m_cellVec2, m_topology, reco_skim_cfg_mod::maxSize, and HGCalTopology::tileTrapezoid().

                                                                     {
   unsigned int cellIndex = indexFor(detid);
   GlobalPoint glob;
   unsigned int maxSize = (m_topology.tileTrapezoid() ? m_cellVec2.size() : m_cellVec.size());
   if (cellIndex < maxSize) {
     const HepGeom::Point3D<float> lcoord(0, 0, 0);
     if (m_topology.tileTrapezoid()) {
       glob = m_cellVec2[cellIndex].getPosition(lcoord);
     } else {
       glob = m_cellVec[cellIndex].getPosition(lcoord);
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "getPositionTrap:: ID " << std::hex << detid.rawId() << std::dec << " index "
                                   << cellIndex << " Global " << glob;
 #endif
   }
   return glob;
 }

◆ indexFor()

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

overrideprotectedvirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 651 of file HGCalGeometry.cc.

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

Referenced by get8Corners(), getCorners(), getNewCorners(), getPosition(), and getWaferPosition().

                                                              {
   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;
 }

◆ initializeParms()

void HGCalGeometry::initializeParms ( )

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 47 of file HGCalGeometry.cc.

47 {}

◆ localCorners()

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

Definition at line 49 of file HGCalGeometry.cc.

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

                                                                                            {
   if (m_det == DetId::HGCalHSc) {
     FlatTrd::localCorners(lc, pv, ref);
   } else {
     FlatHexagon::localCorners(lc, pv, ref);
   }
 }

◆ neighborZ() [1/2]

DetId HGCalGeometry::neighborZ	(	const DetId &	idin,
		const GlobalVector &	p
	)		const

Definition at line 467 of file HGCalGeometry.cc.

References HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalDDDConstants::firstLayer(), getClosestCell(), getPosition(), HGCalDDDConstants::lastLayer(), m_topology, AlCaHLTBitMon_ParallelJobs::p, alignCSCRings::r, HGCalDDDConstants::rangeR(), topology(), findQualityFiles::v, HGCalDDDConstants::waferZ(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), z, and PV3DBase< T, PVType, FrameType >::z().

                                                                                     {
   DetId idnew;
   HGCalTopology::DecodedDetId id = m_topology.decode(idin);
   int lay = ((momentum.z() * id.zSide > 0) ? (id.iLay + 1) : (id.iLay - 1));
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "neighborz1:: ID " << id.iLay << ":" << id.iSec1 << ":" << id.iSec2 << ":"
                                 << id.iCell1 << ":" << id.iCell2 << " New Layer " << lay << " Range "
                                 << m_topology.dddConstants().firstLayer() << ":"
                                 << m_topology.dddConstants().lastLayer(true) << " pz " << momentum.z();
 #endif
   if ((lay >= m_topology.dddConstants().firstLayer()) && (lay <= m_topology.dddConstants().lastLayer(true)) &&
       (momentum.z() != 0.0)) {
     GlobalPoint v = getPosition(idin, false);
     double z = id.zSide * m_topology.dddConstants().waferZ(lay, true);
     double grad = (z - v.z()) / momentum.z();
     GlobalPoint p(v.x() + grad * momentum.x(), v.y() + grad * momentum.y(), z);
     double r = p.perp();
     auto rlimit = topology().dddConstants().rangeR(z, true);
     if (r >= rlimit.first && r <= rlimit.second)
       idnew = getClosestCell(p);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "neighborz1:: Position " << v << " New Z " << z << ":" << grad << " new position "
                                   << p << " r-limit " << rlimit.first << ":" << rlimit.second;
 #endif
   }
   return idnew;
 }

◆ neighborZ() [2/2]

DetId HGCalGeometry::neighborZ	(	const DetId &	idin,
		const MagneticField *	bField,
		int	charge,
		const GlobalVector &	momentum
	)		const

Definition at line 495 of file HGCalGeometry.cc.

References alongMomentum, Calorimetry_cff::bField, Plane::build(), ALCARECOTkAlJpsiMuMu_cff::charge, HGCalTopology::dddConstants(), HGCalTopology::decode(), HGCalDDDConstants::firstLayer(), getClosestCell(), getPosition(), TrajectoryStateOnSurface::globalPosition(), TrajectoryStateOnSurface::isValid(), HGCalDDDConstants::lastLayer(), M_PI, m_topology, AlCaHLTBitMon_ParallelJobs::p, AnalyticalPropagator::propagate(), alignCSCRings::r, HGCalDDDConstants::rangeR(), topology(), findQualityFiles::v, HGCalDDDConstants::waferZ(), z, and PV3DBase< T, PVType, FrameType >::z().

                                                                    {
   DetId idnew;
   HGCalTopology::DecodedDetId id = m_topology.decode(idin);
   int lay = ((momentum.z() * id.zSide > 0) ? (id.iLay + 1) : (id.iLay - 1));
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "neighborz2:: ID " << id.iLay << ":" << id.iSec1 << ":" << id.iSec2 << ":"
                                 << id.iCell1 << ":" << id.iCell2 << " New Layer " << lay << " Range "
                                 << m_topology.dddConstants().firstLayer() << ":"
                                 << m_topology.dddConstants().lastLayer(true) << " pz " << momentum.z();
 #endif
   if ((lay >= m_topology.dddConstants().firstLayer()) && (lay <= m_topology.dddConstants().lastLayer(true)) &&
       (momentum.z() != 0.0)) {
     GlobalPoint v = getPosition(idin, false);
     double z = id.zSide * m_topology.dddConstants().waferZ(lay, true);
     FreeTrajectoryState fts(v, momentum, charge, bField);
     Plane::PlanePointer nPlane = Plane::build(Plane::PositionType(0, 0, z), Plane::RotationType());
     AnalyticalPropagator myAP(bField, alongMomentum, 2 * M_PI);
     TrajectoryStateOnSurface tsos = myAP.propagate(fts, *nPlane);
     GlobalPoint p;
     auto rlimit = topology().dddConstants().rangeR(z, true);
     if (tsos.isValid()) {
       p = tsos.globalPosition();
       double r = p.perp();
       if (r >= rlimit.first && r <= rlimit.second)
         idnew = getClosestCell(p);
     }
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "neighborz2:: Position " << v << " New Z " << z << ":" << charge << ":"
                                   << tsos.isValid() << " new position " << p << " r limits " << rlimit.first << ":"
                                   << rlimit.second;
 #endif
   }
   return idnew;
 }

◆ newCell()

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

overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 57 of file HGCalGeometry.cc.

References HGCalDDDConstants::cellInLayer(), hgcalTBTopologyTester_cfi::cells, CaloSubdetectorGeometry::cornersMgr(), HGCalTopology::dddConstants(), TauDecayModes::dec, HGCalTopology::decode(), DetId::det(), hcalRecHitTable_cff::detId, HGCalTopology::detId2denseGeomId(), HGCalTopology::encode(), validate-o2o-wbm::f1, validate-o2o-wbm::f2, DetId::Forward, getGeometryDetId(), DetId::HGCalHSc, HGCalTopology::isHFNose(), HGCScintillatorDetId::layer(), m_cellVec, m_cellVec2, m_det, m_topology, m_validGeomIds, CaloSubdetectorGeometry::m_validIds, HGCalDDDConstants::numberCellsHexagon(), DetId::rawId(), HGCScintillatorDetId::ring(), HGCScintillatorDetId::setSiPM(), HGCScintillatorDetId::setType(), HGCalTopology::tileTrapezoid(), HGCalDDDConstants::tileType(), findQualityFiles::v, HGCalTopology::valid(), and HGCalTopology::waferHexagon6().

                                                                                                                    {
   DetId geomId = getGeometryDetId(detId);
   int cells(0);
   HGCalTopology::DecodedDetId id = m_topology.decode(detId);
   if (m_topology.waferHexagon6()) {
     cells = m_topology.dddConstants().numberCellsHexagon(id.iSec1);
 #ifdef EDM_ML_DEBUG
     edm::LogVerbatim("HGCalGeom") << "NewCell " << HGCalDetId(detId) << " GEOM " << HGCalDetId(geomId);
 #endif
   } else if (m_topology.tileTrapezoid()) {
     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
     if (detId.det() == DetId::Forward)
       edm::LogVerbatim("HGCalGeom") << "NewCell " << HFNoseDetId(detId) << " GEOM " << HFNoseDetId(geomId);
     else
       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 (m_topology.waferHexagon6()) {
     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 (m_topology.tileTrapezoid()) {
     DetId idc = m_topology.encode(id);
     if (m_topology.valid(idc)) {
       HGCScintillatorDetId hid(idc);
       std::pair<int, int> typm = m_topology.dddConstants().tileType(hid.layer(), hid.ring(), 0);
       if (typm.first >= 0) {
         hid.setType(typm.first);
         hid.setSiPM(typm.second);
         idc = static_cast<DetId>(hid);
       }
       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, id.zSide, true)) {
             m_validIds.emplace_back(idc);
 #ifdef EDM_ML_DEBUG
             ++cellSelect;
             if (idc.det() == DetId::Forward)
               edm::LogVerbatim("HGCalGeom") << "Valid Id [" << u << ", " << v << "] " << HFNoseDetId(idc);
             else
               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 (m_topology.waferHexagon6()) {
     edm::LogVerbatim("HGCalGeom") << "ID: " << HGCalDetId(detId) << " with valid DetId from " << nOld << " to " << nNew;
   } else if (m_topology.tileTrapezoid()) {
     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
 }

◆ present()

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

overridevirtual

is this detid present in the geometry?

Reimplemented from CaloSubdetectorGeometry.

Definition at line 196 of file HGCalGeometry.cc.

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

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

◆ producerTag()

static std::string HGCalGeometry::producerTag ( )

inlinestatic

Definition at line 110 of file HGCalGeometry.h.

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

110 { return "HGCal"; }

◆ sizeForDenseIndex() [1/2]

unsigned int CaloSubdetectorGeometry::sizeForDenseIndex

protected

Definition at line 257 of file CaloSubdetectorGeometry.cc.

                                                                              {
   return CaloGenericDetId(id).sizeForDenseIndexing();
 }

◆ sizeForDenseIndex() [2/2]

unsigned int HGCalGeometry::sizeForDenseIndex ( ) const

protected

Definition at line 664 of file HGCalGeometry.cc.

References m_topology, and HGCalTopology::totalGeomModules().

664 { return m_topology.totalGeomModules(); }

HGCalTopology::totalGeomModules

unsigned int totalGeomModules() const

Definition: HGCalTopology.h:93

HGCalGeometry::m_topology

const HGCalTopology & m_topology

Definition: HGCalGeometry.h:139

◆ sortDetIds()

void HGCalGeometry::sortDetIds ( void )

Definition at line 775 of file HGCalGeometry.cc.

References CaloSubdetectorGeometry::m_validIds, and jetUpdater_cfi::sort.

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

◆ topology()

const HGCalTopology& HGCalGeometry::topology ( ) const

inline

Definition at line 113 of file HGCalGeometry.h.

References m_topology.

Referenced by PFRecHitQTestThresholdInThicknessNormalizedMIPs::beginEvent(), HGCDigitizer::checkPosition(), HGCalTriggerGeometryBase::eeTopology(), HGCalTriggerGeometryBase::fhTopology(), HGCalTriggerGeometryBase::hscTopology(), HGCalTriggerGeometryBase::hsiTopology(), hgcal::RecHitTools::maskCell(), neighborZ(), HGCalTriggerGeometryBase::noseTopology(), HGCalRecHitWorkerSimple::set(), HGCalRadiationMap::setGeometry(), and HGCalUncalibRecHitRecWeightsAlgo< HGCDataFrame >::setGeometry().

113 { return m_topology; }

HGCalGeometry::m_topology

const HGCalTopology & m_topology

Definition: HGCalGeometry.h:139

Member Data Documentation

◆ k_fac1

constexpr double HGCalGeometry::k_fac1 = 0.5

staticprivate

Definition at line 136 of file HGCalGeometry.h.

Referenced by get8Corners(), getCorners(), and getNewCorners().

◆ k_fac2

constexpr double HGCalGeometry::k_fac2 = 1.0 / 3.0

staticprivate

Definition at line 137 of file HGCalGeometry.h.

Referenced by get8Corners(), getCorners(), and getNewCorners().

◆ k_half

constexpr double HGCalGeometry::k_half = 0.5

staticprivate

Definition at line 135 of file HGCalGeometry.h.

Referenced by get8Corners(), getCorners(), and getNewCorners().

◆ k_NumberOfParametersPerHex

constexpr unsigned int HGCalGeometry::k_NumberOfParametersPerHex = 3

static

Definition at line 44 of file HGCalGeometry.h.

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

◆ k_NumberOfParametersPerShape

constexpr unsigned int HGCalGeometry::k_NumberOfParametersPerShape = 3

static

Definition at line 45 of file HGCalGeometry.h.

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

◆ k_NumberOfParametersPerTrd

constexpr unsigned int HGCalGeometry::k_NumberOfParametersPerTrd = 12

static

Definition at line 43 of file HGCalGeometry.h.

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

◆ k_NumberOfShapes

constexpr unsigned int HGCalGeometry::k_NumberOfShapes = 100

static

Definition at line 46 of file HGCalGeometry.h.

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

◆ k_NumberOfShapesTrd

constexpr unsigned int HGCalGeometry::k_NumberOfShapesTrd = 1000

static

Definition at line 47 of file HGCalGeometry.h.

Referenced by HGCalGeometryLoader::build().

◆ m_cellVec

CellVec HGCalGeometry::m_cellVec

private

Definition at line 140 of file HGCalGeometry.h.

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

◆ m_cellVec2

CellVec2 HGCalGeometry::m_cellVec2

private

Definition at line 141 of file HGCalGeometry.h.

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

◆ m_det

DetId::Detector HGCalGeometry::m_det

private

Definition at line 143 of file HGCalGeometry.h.

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

◆ m_subdet

ForwardSubdetector HGCalGeometry::m_subdet

private

Definition at line 144 of file HGCalGeometry.h.

Referenced by cellElement().

◆ m_topology

const HGCalTopology& HGCalGeometry::m_topology

private

Definition at line 139 of file HGCalGeometry.h.

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

◆ m_validGeomIds

std::vector<DetId> HGCalGeometry::m_validGeomIds

private

Definition at line 142 of file HGCalGeometry.h.

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

◆ twoBysqrt3_

const double HGCalGeometry::twoBysqrt3_

private

Definition at line 145 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

Static Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ AlignedRecord

◆ CCGFloat

◆ CellVec

◆ CellVec2

◆ CornersVec

◆ DetIdSet

◆ PGeometryRecord

◆ Pt3D

◆ Pt3DVec

Constructor & Destructor Documentation

◆ HGCalGeometry()

◆ ~HGCalGeometry()

Member Function Documentation

◆ addValidID()

◆ cellElement()

◆ cellGeomPtr() [1/2]

◆ cellGeomPtr() [2/2]

◆ dbString()

◆ fillNamedParams()

◆ get8Corners()

◆ getArea()

◆ getCells()

◆ getClosestCell()

◆ getClosestCellHex()

◆ getClosestCellIndex() [1/2]

◆ getClosestCellIndex() [2/2]

◆ getCorners()

◆ getGeometry()

◆ getGeometryDetId()

◆ getGeometryRawPtr()

◆ getNewCorners()

◆ getPosition() [1/2]

◆ getPosition() [2/2]

◆ getSummary()

◆ getValidDetIds()

◆ getValidGeomDetIds()

◆ getWaferPosition()

◆ indexFor()

◆ initializeParms()

◆ localCorners()

◆ neighborZ() [1/2]

◆ neighborZ() [2/2]

◆ newCell()

◆ present()

◆ producerTag()

◆ sizeForDenseIndex() [1/2]

◆ sizeForDenseIndex() [2/2]

◆ sortDetIds()

◆ topology()

Member Data Documentation

◆ k_fac1

◆ k_fac2

◆ k_half

◆ k_NumberOfParametersPerHex

◆ k_NumberOfParametersPerShape

◆ k_NumberOfParametersPerTrd

◆ k_NumberOfShapes

◆ k_NumberOfShapesTrd

◆ m_cellVec

◆ m_cellVec2

◆ m_det

◆ m_subdet

◆ m_topology

◆ m_validGeomIds

◆ twoBysqrt3_