#include <HcalDDDGeometry.h>

Inheritance diagram for HcalDDDGeometry:

Public Types
typedef std::vector< IdealObliquePrism >	HBCellVec

typedef std::vector< IdealObliquePrism >	HECellVec

typedef std::vector< IdealZPrism >	HFCellVec

typedef std::vector< IdealObliquePrism >	HOCellVec

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
DetId	getClosestCell (const GlobalPoint &r) 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...

	HcalDDDGeometry (const HcalTopology &theTopo)

int	insertCell (std::vector< HcalCellType > const &)

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

	~HcalDDDGeometry () override
	The HcalDDDGeometry will delete all its cell geometries at destruction time. More...

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 DetIdSet	getCells (const GlobalPoint &r, double dR) const
	Get a list of all cells within a dR of the given cell. More...

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

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

virtual void	getSummary (TrVec &trVector, IVec &iVector, DimVec &dimVector, IVec &dinsVector) const

virtual void	initializeParms ()

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	present (const DetId &id) const
	is this detid present in the geometry? More...

virtual bool	valid (const DetId &id) const

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

Protected Member Functions
CaloCellGeometryPtr	getGeometryRawPtr (uint32_t index) const override

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

virtual CellPtr	cellGeomPtr (uint32_t index) const

virtual unsigned int	indexFor (const DetId &id) const

virtual unsigned int	sizeForDenseIndex (const DetId &id) const

Private Member Functions
void	fillDetIds () const

void	increaseReserve (unsigned int extra)

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

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

void	sortValidIds ()

Private Attributes
double	etaMax_

std::vector< HcalCellType >	hcalCells_

std::vector< DetId >	m_emptyIds

std::atomic< bool >	m_filledDetIds

HBCellVec	m_hbCellVec

std::vector< DetId >	m_hbIds

HECellVec	m_heCellVec

std::vector< DetId >	m_heIds

HFCellVec	m_hfCellVec

std::vector< DetId >	m_hfIds

HOCellVec	m_hoCellVec

std::vector< DetId >	m_hoIds

const HcalTopology &	topo_

Friends
class	HcalDDDGeometryLoader

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 17 of file HcalDDDGeometry.h.

Member Typedef Documentation

◆ HBCellVec

typedef std::vector<IdealObliquePrism> HcalDDDGeometry::HBCellVec

Definition at line 21 of file HcalDDDGeometry.h.

◆ HECellVec

typedef std::vector<IdealObliquePrism> HcalDDDGeometry::HECellVec

Definition at line 22 of file HcalDDDGeometry.h.

◆ HFCellVec

typedef std::vector<IdealZPrism> HcalDDDGeometry::HFCellVec

Definition at line 24 of file HcalDDDGeometry.h.

◆ HOCellVec

typedef std::vector<IdealObliquePrism> HcalDDDGeometry::HOCellVec

Definition at line 23 of file HcalDDDGeometry.h.

Constructor & Destructor Documentation

◆ HcalDDDGeometry()

HcalDDDGeometry::HcalDDDGeometry ( const HcalTopology & theTopo )

explicit

Definition at line 12 of file HcalDDDGeometry.cc.

     : topo_(topo),
       etaMax_(0),
       m_hbCellVec(topo.getHBSize()),
       m_heCellVec(topo.getHESize()),
       m_hoCellVec(topo.getHOSize()),
       m_hfCellVec(topo.getHFSize()),
       m_filledDetIds(false) {}

◆ ~HcalDDDGeometry()

HcalDDDGeometry::~HcalDDDGeometry ( )

override

The HcalDDDGeometry will delete all its cell geometries at destruction time.

Definition at line 21 of file HcalDDDGeometry.cc.

21 {}

Member Function Documentation

◆ fillDetIds()

void HcalDDDGeometry::fillDetIds ( ) const

private

Definition at line 23 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::getValidDetIds(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, EcalPhiSymFlatTableProducers_cfi::id, m_emptyIds, m_filledDetIds, m_hbIds, m_heIds, m_hfIds, m_hoIds, s_fillLock, and jetUpdater_cfi::sort.

Referenced by getValidDetIds().

                                        {
   std::lock_guard<std::mutex> guard(s_fillLock);
   if (m_filledDetIds) {
     //another thread already did the work
     return;
   }
   const std::vector<DetId>& baseIds(CaloSubdetectorGeometry::getValidDetIds());
   for (unsigned int i(0); i != baseIds.size(); ++i) {
     const DetId id(baseIds[i]);
     if (id.subdetId() == HcalBarrel) {
       m_hbIds.emplace_back(id);
     } else {
       if (id.subdetId() == HcalEndcap) {
         m_heIds.emplace_back(id);
       } else {
         if (id.subdetId() == HcalOuter) {
           m_hoIds.emplace_back(id);
         } else {
           if (id.subdetId() == HcalForward) {
             m_hfIds.emplace_back(id);
           }
         }
       }
     }
   }
   std::sort(m_hbIds.begin(), m_hbIds.end());
   std::sort(m_heIds.begin(), m_heIds.end());
   std::sort(m_hoIds.begin(), m_hoIds.end());
   std::sort(m_hfIds.begin(), m_hfIds.end());
 
   m_emptyIds.resize(0);
   m_filledDetIds = true;
 }

◆ getClosestCell()

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

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 69 of file HcalDDDGeometry.cc.

References ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), etaMax_, hcalCells_, HcalForward, M_PI, phi, alignCSCRings::r, CosmicsPD_Skims::radius, and z.

                                                                 {
   constexpr double twopi = M_PI + M_PI;
 
   // Now find the closest eta_bin, eta value of a bin i is average
   // of eta[i] and eta[i-1]
   double abseta = fabs(r.eta());
   double phi = r.phi();
   if (phi < 0)
     phi += twopi;
   double radius = r.mag();
   double z = fabs(r.z());
 
 #ifdef EDM_ML_DEBUG
   constexpr double deg = M_PI / 180.;
   edm::LogVerbatim("HCalGeom") << "HcalDDDGeometry::getClosestCell for eta " << r.eta() << " phi " << phi / deg << " z "
                                << r.z() << " radius " << radius;
 #endif
   HcalDetId bestId;
   if (abseta <= etaMax_) {
     for (const auto& hcalCell : hcalCells_) {
       if (abseta >= hcalCell.etaMin() && abseta <= hcalCell.etaMax()) {
         HcalSubdetector bc = hcalCell.detType();
         int etaring = hcalCell.etaBin();
         int phibin = 0;
         if (hcalCell.unitPhi() == 4) {
           // rings 40 and 41 are offset wrt the other phi numbering
           //  1        1         1         2
           //  ------------------------------
           //  72       36        36        1
           phibin =
               static_cast<int>((phi + hcalCell.phiOffset() + 0.5 * hcalCell.phiBinWidth()) / hcalCell.phiBinWidth());
           if (phibin == 0)
             phibin = hcalCell.nPhiBins();
           phibin = phibin * 4 - 1;
         } else {
           phibin = static_cast<int>((phi + hcalCell.phiOffset()) / hcalCell.phiBinWidth()) + 1;
           // convert to the convention of numbering 1,3,5, in 36 phi bins
           phibin = (phibin - 1) * (hcalCell.unitPhi()) + 1;
         }
 
         int dbin = 1;
         int etabin = (r.z() > 0) ? etaring : -etaring;
         if (bc == HcalForward) {
           bestId = HcalDetId(bc, etabin, phibin, dbin);
           break;
         } else {
           double rz = z;
           if (hcalCell.depthType())
             rz = radius;
           if (rz < hcalCell.depthMax()) {
             dbin = hcalCell.depthSegment();
             bestId = HcalDetId(bc, etabin, phibin, dbin);
             break;
           }
         }
       }
     }
   }
 
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDGeometry::getClosestCell " << bestId;
 #endif
   return bestId;
 }

◆ getGeometryRawPtr()

CaloCellGeometryPtr HcalDDDGeometry::getGeometryRawPtr ( uint32_t index ) const

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 187 of file HcalDDDGeometry.cc.

References cropTnPTrees::din, m_hbCellVec, m_heCellVec, m_hfCellVec, m_hoCellVec, and CaloCellGeometry::param().

                                                                          {
   // Modify the RawPtr class
   const CaloCellGeometry* cell(nullptr);
   if (m_hbCellVec.size() > din) {
     cell = (&m_hbCellVec[din]);
   } else if (m_hbCellVec.size() + m_heCellVec.size() > din) {
     const unsigned int ind(din - m_hbCellVec.size());
     cell = (&m_heCellVec[ind]);
   } else if (m_hbCellVec.size() + m_heCellVec.size() + m_hoCellVec.size() > din) {
     const unsigned int ind(din - m_hbCellVec.size() - m_heCellVec.size());
     cell = (&m_hoCellVec[ind]);
   } else if (m_hbCellVec.size() + m_heCellVec.size() + m_hoCellVec.size() + m_hfCellVec.size() > din) {
     const unsigned int ind(din - m_hbCellVec.size() - m_heCellVec.size() - m_hoCellVec.size());
     cell = (&m_hfCellVec[ind]);
   }
 
   return CaloCellGeometryPtr((nullptr == cell || nullptr == cell->param()) ? nullptr : cell);
 }

◆ getValidDetIds()

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

overridevirtual

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 57 of file HcalDDDGeometry.cc.

References fillDetIds(), CaloSubdetectorGeometry::getValidDetIds(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, m_emptyIds, m_filledDetIds, m_hbIds, m_heIds, m_hfIds, and m_hoIds.

                                                                                            {
   if (0 != subdet && not m_filledDetIds)
     fillDetIds();
   return (0 == subdet
               ? CaloSubdetectorGeometry::getValidDetIds()
               : (HcalBarrel == subdet
                      ? m_hbIds
                      : (HcalEndcap == subdet
                             ? m_heIds
                             : (HcalOuter == subdet ? m_hoIds : (HcalForward == subdet ? m_hfIds : m_emptyIds)))));
 }

◆ increaseReserve()

void HcalDDDGeometry::increaseReserve ( unsigned int extra )

private

Definition at line 206 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::m_validIds.

206 { m_validIds.reserve(m_validIds.size() + extra); }

CaloSubdetectorGeometry::m_validIds

std::vector< DetId > m_validIds

Definition: CaloSubdetectorGeometry.h:125

◆ insertCell()

int HcalDDDGeometry::insertCell ( std::vector< HcalCellType > const & cells )

Definition at line 134 of file HcalDDDGeometry.cc.

References hgcalTBTopologyTester_cfi::cells, etaMax_, hcalCells_, and EgammaValidation_cff::num.

                                                                     {
   hcalCells_.insert(hcalCells_.end(), cells.begin(), cells.end());
   int num = static_cast<int>(hcalCells_.size());
   for (const auto& cell : cells) {
     if (cell.etaMax() > etaMax_)
       etaMax_ = cell.etaMax();
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HCalGeom") << "HcalDDDGeometry::insertCell " << cells.size() << " cells inserted == Total " << num
                                << " EtaMax = " << etaMax_;
 #endif
   return num;
 }

◆ newCell()

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

overridevirtual

Implements CaloSubdetectorGeometry.

Definition at line 175 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::addValidID(), hcalRecHitTable_cff::detId, validate-o2o-wbm::f1, validate-o2o-wbm::f2, and newCellImpl().

                                                                                                                    {
   newCellImpl(f1, f2, f3, parm, detId);
   addValidID(detId);
 }

◆ newCellFast()

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

private

Definition at line 181 of file HcalDDDGeometry.cc.

References hcalRecHitTable_cff::detId, validate-o2o-wbm::f1, validate-o2o-wbm::f2, CaloSubdetectorGeometry::m_validIds, and newCellImpl().

                                                                                                                    {
   newCellImpl(f1, f2, f3, parm, detId);
   m_validIds.emplace_back(detId);
 }

◆ newCellImpl()

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

private

Definition at line 148 of file HcalDDDGeometry.cc.

References cms::cuda::assert(), CaloSubdetectorGeometry::cornersMgr(), hcalRecHitTable_cff::detId, HcalTopology::detId2denseId(), cropTnPTrees::din, IdealZPrism::EM, validate-o2o-wbm::f1, IdealZPrism::HADR, DetId::Hcal, HcalBarrel, HcalEndcap, HcalOuter, m_hbCellVec, m_heCellVec, m_hfCellVec, m_hoCellVec, and topo_.

Referenced by newCell(), and newCellFast().

                                                                                                                    {
   assert(detId.det() == DetId::Hcal);
 
   const unsigned int din(topo_.detId2denseId(detId));
 
   HcalDetId hId(detId);
 
   if (hId.subdet() == HcalBarrel) {
     m_hbCellVec[din] = IdealObliquePrism(f1, cornersMgr(), parm);
   } else {
     if (hId.subdet() == HcalEndcap) {
       const unsigned int index(din - m_hbCellVec.size());
       m_heCellVec[index] = IdealObliquePrism(f1, cornersMgr(), parm);
     } else {
       if (hId.subdet() == HcalOuter) {
         const unsigned int index(din - m_hbCellVec.size() - m_heCellVec.size());
         m_hoCellVec[index] = IdealObliquePrism(f1, cornersMgr(), parm);
       } else {  // assuming HcalForward here!
         const unsigned int index(din - m_hbCellVec.size() - m_heCellVec.size() - m_hoCellVec.size());
         m_hfCellVec[index] =
             IdealZPrism(f1, cornersMgr(), parm, hId.depth() == 1 ? IdealZPrism::EM : IdealZPrism::HADR);
       }
     }
   }
 }

◆ sortValidIds()

void HcalDDDGeometry::sortValidIds ( )

private

Definition at line 208 of file HcalDDDGeometry.cc.

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

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

jetUpdater_cfi.sort

sort

Definition: jetUpdater_cfi.py:30

CaloSubdetectorGeometry::m_validIds

std::vector< DetId > m_validIds

Definition: CaloSubdetectorGeometry.h:125

Friends And Related Function Documentation

◆ HcalDDDGeometryLoader

friend class HcalDDDGeometryLoader

friend

Definition at line 19 of file HcalDDDGeometry.h.

Member Data Documentation

◆ etaMax_

double HcalDDDGeometry::etaMax_

private

Definition at line 67 of file HcalDDDGeometry.h.

Referenced by getClosestCell(), and insertCell().

◆ hcalCells_

std::vector<HcalCellType> HcalDDDGeometry::hcalCells_

private

Definition at line 59 of file HcalDDDGeometry.h.

Referenced by getClosestCell(), and insertCell().

◆ m_emptyIds

std::vector<DetId> HcalDDDGeometry::m_emptyIds

mutableprivate

Definition at line 64 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

◆ m_filledDetIds

std::atomic<bool> HcalDDDGeometry::m_filledDetIds

mutableprivate

Definition at line 73 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

◆ m_hbCellVec

HBCellVec HcalDDDGeometry::m_hbCellVec

private

Definition at line 69 of file HcalDDDGeometry.h.

Referenced by getGeometryRawPtr(), and newCellImpl().

◆ m_hbIds

std::vector<DetId> HcalDDDGeometry::m_hbIds

mutableprivate

Definition at line 60 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

◆ m_heCellVec

HECellVec HcalDDDGeometry::m_heCellVec

private

Definition at line 70 of file HcalDDDGeometry.h.

Referenced by getGeometryRawPtr(), and newCellImpl().

◆ m_heIds

std::vector<DetId> HcalDDDGeometry::m_heIds

mutableprivate

Definition at line 61 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

◆ m_hfCellVec

HFCellVec HcalDDDGeometry::m_hfCellVec

private

Definition at line 72 of file HcalDDDGeometry.h.

Referenced by getGeometryRawPtr(), and newCellImpl().

◆ m_hfIds

std::vector<DetId> HcalDDDGeometry::m_hfIds

mutableprivate

Definition at line 63 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

◆ m_hoCellVec

HOCellVec HcalDDDGeometry::m_hoCellVec

private

Definition at line 71 of file HcalDDDGeometry.h.

Referenced by getGeometryRawPtr(), and newCellImpl().

◆ m_hoIds

std::vector<DetId> HcalDDDGeometry::m_hoIds

mutableprivate

Definition at line 62 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

◆ topo_

const HcalTopology& HcalDDDGeometry::topo_

private

Definition at line 66 of file HcalDDDGeometry.h.

Referenced by newCellImpl().

Public Types

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Friends

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ HBCellVec

◆ HECellVec

◆ HFCellVec

◆ HOCellVec

Constructor & Destructor Documentation

◆ HcalDDDGeometry()

◆ ~HcalDDDGeometry()

Member Function Documentation

◆ fillDetIds()

◆ getClosestCell()

◆ getGeometryRawPtr()

◆ getValidDetIds()

◆ increaseReserve()

◆ insertCell()

◆ newCell()

◆ newCellFast()

◆ newCellImpl()

◆ sortValidIds()

Friends And Related Function Documentation

◆ HcalDDDGeometryLoader

Member Data Documentation

◆ etaMax_

◆ hcalCells_

◆ m_emptyIds

◆ m_filledDetIds

◆ m_hbCellVec

◆ m_hbIds

◆ m_heCellVec

◆ m_heIds

◆ m_hfCellVec

◆ m_hfIds

◆ m_hoCellVec

◆ m_hoIds

◆ topo_