#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< CaloCellGeometry const * >	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 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 *	cellGeomPtr (uint32_t index) const override

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

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

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

Definition at line 24 of file HcalDDDGeometry.h.

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

Definition at line 26 of file HcalDDDGeometry.h.

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

Definition at line 25 of file HcalDDDGeometry.h.

Constructor & Destructor Documentation

HcalDDDGeometry::HcalDDDGeometry ( const HcalTopology & theTopo )

explicit

Definition at line 10 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.

22 {

23 }

Member Function Documentation

const CaloCellGeometry * HcalDDDGeometry::cellGeomPtr ( uint32_t index ) const

overrideprotectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 224 of file HcalDDDGeometry.cc.

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

 {
   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 index ( din - m_hbCellVec.size() ) ;
       cell = &m_heCellVec[ index ] ;
     }
     else
     {
       if( m_hbCellVec.size() +
       m_heCellVec.size() +
       m_hoCellVec.size() > din )
       {
     const unsigned int index ( din 
                    - m_hbCellVec.size() 
                    - m_heCellVec.size() ) ;
     cell = &m_hoCellVec[ index ] ;
       }
       else
       {
     if( m_hbCellVec.size() +
         m_heCellVec.size() +
         m_hoCellVec.size() +
         m_hfCellVec.size() > din )
     {
       const unsigned int index ( din 
                      - m_hbCellVec.size() 
                      - m_heCellVec.size() 
                      - m_hoCellVec.size() ) ;
       cell = &m_hfCellVec[ index ] ;
     }
       }
     }
   }
   return ( nullptr == cell || nullptr == cell->param() ? nullptr : cell ) ;
 }

void HcalDDDGeometry::fillDetIds ( ) const

private

Definition at line 26 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::getValidDetIds(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, mps_fire::i, triggerObjects_cff::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 86 of file HcalDDDGeometry.cc.

References constexpr, PV3DBase< T, PVType, FrameType >::eta(), etaMax_, hcalCells_, HcalForward, LogDebug, M_PI, PV3DBase< T, PVType, FrameType >::mag(), phi, PV3DBase< T, PVType, FrameType >::phi(), TCMET_cfi::radius, z, and PV3DBase< T, PVType, FrameType >::z().

 {
   constexpr double twopi = M_PI+M_PI;
   constexpr double deg   = M_PI/180.;
 
   // 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());
 
   LogDebug("HCalGeom") << "HcalDDDGeometry::getClosestCell for eta "
                << r.eta() << " phi " << phi/deg << " z " << r.z()
                << " radius " << radius;
   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;
       }
     }
       }
     }
   }
 
   LogDebug("HCalGeom") << "HcalDDDGeometry::getClosestCell " << bestId;
 
   return bestId;
 }

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

References CaloSubdetectorGeometry::m_validIds.

Referenced by HcalDDDGeometryLoader::fill().

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

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

Definition at line 150 of file HcalDDDGeometry.cc.

References etaMax_, hcalCells_, LogDebug, 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();
   }
 
   LogDebug("HCalGeom") << "HcalDDDGeometry::insertCell " << cells.size()
                << " cells inserted == Total " << num
                << " EtaMax = " << etaMax_;
   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 202 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 213 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 165 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::cornersMgr(), DetId::det(), HcalTopology::detId2denseId(), cropTnPTrees::din, IdealZPrism::EM, IdealZPrism::HADR, DetId::Hcal, HcalBarrel, HcalEndcap, HcalOuter, diffTreeTool::index, 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 ) ;
       }
     }
   }
 }

void HcalDDDGeometry::sortValidIds ( )

private

Definition at line 273 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::m_validIds.

Referenced by HcalDDDGeometryLoader::load().

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

Friends And Related Function Documentation

friend class HcalDDDGeometryLoader

friend

Definition at line 21 of file HcalDDDGeometry.h.

Member Data Documentation

double HcalDDDGeometry::etaMax_

private

Definition at line 77 of file HcalDDDGeometry.h.

Referenced by getClosestCell(), and insertCell().

std::vector<HcalCellType> HcalDDDGeometry::hcalCells_

private

Definition at line 69 of file HcalDDDGeometry.h.

Referenced by getClosestCell(), and insertCell().

std::vector<DetId> HcalDDDGeometry::m_emptyIds

mutableprivate

Definition at line 74 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

std::atomic<bool> HcalDDDGeometry::m_filledDetIds

mutableprivate

Definition at line 83 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

HBCellVec HcalDDDGeometry::m_hbCellVec

private

Definition at line 79 of file HcalDDDGeometry.h.

Referenced by cellGeomPtr(), and newCellImpl().

std::vector<DetId> HcalDDDGeometry::m_hbIds

mutableprivate

Definition at line 70 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

HECellVec HcalDDDGeometry::m_heCellVec

private

Definition at line 80 of file HcalDDDGeometry.h.

Referenced by cellGeomPtr(), and newCellImpl().

std::vector<DetId> HcalDDDGeometry::m_heIds

mutableprivate

Definition at line 71 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

HFCellVec HcalDDDGeometry::m_hfCellVec

private

Definition at line 82 of file HcalDDDGeometry.h.

Referenced by cellGeomPtr(), and newCellImpl().

std::vector<DetId> HcalDDDGeometry::m_hfIds

mutableprivate

Definition at line 73 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

HOCellVec HcalDDDGeometry::m_hoCellVec

private

Definition at line 81 of file HcalDDDGeometry.h.

Referenced by cellGeomPtr(), and newCellImpl().

std::vector<DetId> HcalDDDGeometry::m_hoIds

mutableprivate

Definition at line 72 of file HcalDDDGeometry.h.

Referenced by fillDetIds(), and getValidDetIds().

const HcalTopology& HcalDDDGeometry::topo_

private

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

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation