#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

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 std::shared_ptr< const CaloCellGeometry >	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	~CaloSubdetectorGeometry ()
	The base class DOES assume that it owns the CaloCellGeometry objects. More...

Protected Member Functions
const CaloCellGeometry *	getGeometryRawPtr (uint32_t index) const override

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

virtual std::shared_ptr< const CaloCellGeometry >	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

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

Definition at line 21 of file HcalDDDGeometry.h.

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

Definition at line 22 of file HcalDDDGeometry.h.

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

Definition at line 24 of file HcalDDDGeometry.h.

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

Definition at line 23 of file HcalDDDGeometry.h.

Constructor & Destructor Documentation

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

override

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

Definition at line 21 of file HcalDDDGeometry.cc.

21 {}

Member Function Documentation

void HcalDDDGeometry::fillDetIds ( ) const

private

Definition at line 23 of file HcalDDDGeometry.cc.

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

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

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

overridevirtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 69 of file HcalDDDGeometry.cc.

References PV3DBase< T, PVType, FrameType >::eta(), etaMax_, hcalCells_, HcalForward, M_PI, PV3DBase< T, PVType, FrameType >::mag(), phi, PV3DBase< T, PVType, FrameType >::phi(), CosmicsPD_Skims::radius, z, and PV3DBase< T, PVType, FrameType >::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;
 }

const CaloCellGeometry * 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 ((nullptr == cell || nullptr == cell->param()) ? nullptr : cell);
 }

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

void HcalDDDGeometry::increaseReserve ( unsigned int extra )

private

Definition at line 206 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::m_validIds.

Referenced by HcalDDDGeometryLoader::fill().

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

CaloSubdetectorGeometry::m_validIds

std::vector< DetId > m_validIds

Definition: CaloSubdetectorGeometry.h:115

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

Definition at line 134 of file HcalDDDGeometry.cc.

References etaMax_, hcalCells_, and pileupDistInMC::num.

Referenced by HcalDDDGeometryLoader::fill().

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

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(), and newCellImpl().

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

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 CaloSubdetectorGeometry::m_validIds, and newCellImpl().

Referenced by HcalDDDGeometryLoader::makeCell().

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

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(), DetId::det(), HcalTopology::detId2denseId(), cropTnPTrees::din, IdealZPrism::EM, 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);
       }
     }
   }
 }