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

virtual const std::vector < DetId > &	getValidDetIds (DetId::Detector det=DetId::Detector(0), int subdet=0) const
	Get a list of valid detector ids (for the given subdetector) More...

	HcalDDDGeometry (const HcalTopology &theTopo)

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

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

virtual	~HcalDDDGeometry ()
	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 ()

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

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

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

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_

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

Detailed Description

Definition at line 15 of file HcalDDDGeometry.h.

Member Typedef Documentation

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

Definition at line 19 of file HcalDDDGeometry.h.

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

Definition at line 20 of file HcalDDDGeometry.h.

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

Definition at line 22 of file HcalDDDGeometry.h.

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

Definition at line 21 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 ( )

virtual

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

protectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 212 of file HcalDDDGeometry.cc.

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

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

void HcalDDDGeometry::fillDetIds ( ) const

private

Definition at line 26 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::getValidDetIds(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, i, m_emptyIds, m_filledDetIds, m_hbIds, m_heIds, m_hfIds, m_hoIds, s_fillLock, and python.multivaluedict::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.push_back( id ) ;
       }
       else
       {
      if( id.subdetId() == HcalEndcap )
      { 
         m_heIds.push_back( id ) ;
      }
      else
      {
         if( id.subdetId() == HcalOuter )
         { 
            m_hoIds.push_back( id ) ;
         }
         else
         {
            if( id.subdetId() == HcalForward )
            { 
           m_hfIds.push_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

virtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 86 of file HcalDDDGeometry.cc.

References constexpr, PV3DBase< T, PVType, FrameType >::eta(), etaMax_, hcalCells_, HcalForward, i, LogDebug, M_PI, PV3DBase< T, PVType, FrameType >::mag(), phi, PV3DBase< T, PVType, FrameType >::phi(), CosmicsPD_Skims::radius, detailsBasic3DVector::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 (unsigned int i=0; i<hcalCells_.size(); i++) {
       if (abseta >=hcalCells_[i].etaMin() && abseta <=hcalCells_[i].etaMax()) {
     HcalSubdetector bc = hcalCells_[i].detType();
     int etaring = hcalCells_[i].etaBin();
     int phibin  = 0;
     if (hcalCells_[i].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/deg)+hcalCells_[i].phiOffset()+
                      0.5*hcalCells_[i].phiBinWidth())/
                     hcalCells_[i].phiBinWidth());
       if (phibin == 0) phibin = hcalCells_[i].nPhiBins();
       phibin = phibin*4 - 1; 
     } else {
       phibin = static_cast<int>(((phi/deg)+hcalCells_[i].phiOffset())/
                     hcalCells_[i].phiBinWidth()) + 1;
       // convert to the convention of numbering 1,3,5, in 36 phi bins
       phibin = (phibin-1)*(hcalCells_[i].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 (hcalCells_[i].depthType()) rz = radius;
       if (rz < hcalCells_[i].depthMax()) {
         dbin   = hcalCells_[i].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

virtual

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

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

Definition at line 150 of file HcalDDDGeometry.cc.

References etaMax_, hcalCells_, i, LogDebug, and pileupDistInMC::num.

Referenced by HcalDDDGeometryLoader::fill().

                                                                 {
 
   hcalCells_.insert(hcalCells_.end(), cells.begin(), cells.end());
   int num = static_cast<int>(hcalCells_.size());
   for (unsigned int i=0; i<cells.size(); i++) {
     if (cells[i].etaMax() > etaMax_ ) etaMax_ = cells[i].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
	)

virtual

Implements CaloSubdetectorGeometry.

Definition at line 165 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::addValidID(), CaloSubdetectorGeometry::cornersMgr(), DetId::det(), HcalTopology::detId2denseId(), cropTnPTrees::din, DetId::Hcal, HcalBarrel, HcalEndcap, HcalOuter, cmsHarvester::index, m_hbCellVec, m_heCellVec, m_hfCellVec, m_hoCellVec, and topo_.

 {
 
   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 ) ;
       }
     }
   }
   addValidID( detId ) ;
 }