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

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 const CaloCellGeometry *	getGeometry (const DetId &id) const
	Get the cell geometry of a given detector id. Should return false if not found. More...

void	getSummary (TrVec &trVector, IVec &iVector, DimVec &dimVector) 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

Private Attributes
double	deg

double	etaMax_

double	firstHFQuadRing_

std::vector< HcalCellType >	hcalCells_

DetId::Detector	lastReqDet_

int	lastReqSubdet_

HBCellVec	m_hbCellVec

HECellVec	m_heCellVec

HFCellVec	m_hfCellVec

HOCellVec	m_hoCellVec

std::vector< DetId >	m_validIds

double	twopi

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

Member Typedef Documentation

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

Definition at line 22 of file HcalDDDGeometry.h.

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

Definition at line 23 of file HcalDDDGeometry.h.

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

Definition at line 25 of file HcalDDDGeometry.h.

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

Definition at line 24 of file HcalDDDGeometry.h.

Constructor & Destructor Documentation

HcalDDDGeometry::HcalDDDGeometry ( )

explicit

Definition at line 11 of file HcalDDDGeometry.cc.

References deg, M_PI, and twopi.

                                  :
    lastReqDet_(DetId::Detector(0)), 
    lastReqSubdet_(0),
    etaMax_(0),
    firstHFQuadRing_(40) ,
    m_hbCellVec ( HcalDetId::kHBSize ) ,
    m_heCellVec ( HcalDetId::kHESize ) ,
    m_hoCellVec ( HcalDetId::kHOSize ) ,
    m_hfCellVec ( HcalDetId::kHFSize ) 
 {
   twopi = M_PI + M_PI;
   deg   = M_PI/180.;
 }

HcalDDDGeometry::~HcalDDDGeometry ( )

virtual

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

Definition at line 26 of file HcalDDDGeometry.cc.

27 {

28 }

Member Function Documentation

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

protectedvirtual

Implements CaloSubdetectorGeometry.

Definition at line 190 of file HcalDDDGeometry.cc.

References cropTnPTrees::din, getHLTprescales::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 ) ;
 }

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

virtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 65 of file HcalDDDGeometry.cc.

References deg, PV3DBase< T, PVType, FrameType >::eta(), etaMax_, jptDQMConfig_cff::etaMin, hcalCells_, HcalForward, i, LogDebug, PV3DBase< T, PVType, FrameType >::mag(), phi, PV3DBase< T, PVType, FrameType >::phi(), CosmicsPD_Skims::radius, twopi, detailsBasic3DVector::z, and PV3DBase< T, PVType, FrameType >::z().

                                                                 {
 
   // 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 DebugLog
   LogDebug("HCalGeom") << "HcalDDDGeometry::getClosestCell for eta "
                << r.eta() << " phi " << phi/deg << " z " << r.z()
                << " radius " << radius;
 #endif
 
   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;
       }
     }
       }
     }
   }
 #ifdef DebugLog
   LogDebug("HCalGeom") << "HcalDDDGeometry::getClosestCell " << bestId;
 #endif
   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 31 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::getValidDetIds(), i, errorMatrix2Lands_multiChannel::id, lastReqDet_, lastReqSubdet_, LogDebug, m_validIds, and python.multivaluedict::sort().

                                                  {
 
   const std::vector<DetId>& baseIds(CaloSubdetectorGeometry::getValidDetIds());
   if (det    == DetId::Detector( 0 ) && subdet == 0) {
     return baseIds ;
   }
    
   if (lastReqDet_  != det || lastReqSubdet_ != subdet ) {
     lastReqDet_     = det    ;
     lastReqSubdet_  = subdet ;
     m_validIds.clear();
     m_validIds.reserve( baseIds.size() ) ;
   }
 
   if (m_validIds.empty() ) {
     for (unsigned int i = 0 ; i != baseIds.size() ; ++i ) {
       const DetId id ( baseIds[i] );
       if (id.det()      == det   && id.subdetId() == subdet ) { 
     m_validIds.push_back( id ) ;
       }
     }
     std::sort(m_validIds.begin(),m_validIds.end());
   }
   
 #ifdef DebugLog
   LogDebug("HCalGeom") << "HcalDDDGeometry::getValidDetIds: "
                << m_validIds.size() << " valid IDs found for detector "
                << det << " Sub-detector " << subdet;
 #endif
   return m_validIds;
 }

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

Definition at line 128 of file HcalDDDGeometry.cc.

References jptDQMConfig_cff::etaMax, etaMax_, hcalCells_, i, and LogDebug.

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();
   }
 #ifdef DebugLog
   LogDebug("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
	)

virtual

Implements CaloSubdetectorGeometry.

Definition at line 144 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::cornersMgr(), CaloGenericDetId::denseIndex(), cropTnPTrees::din, getHLTprescales::index, CaloGenericDetId::isHB(), CaloGenericDetId::isHcal(), CaloGenericDetId::isHE(), CaloGenericDetId::isHO(), m_hbCellVec, m_heCellVec, m_hfCellVec, m_hoCellVec, and m_validIds.

 {
    const CaloGenericDetId cgid ( detId ) ;
 
    const unsigned int din ( cgid.denseIndex() ) ;
 
    assert( cgid.isHcal() ) ;
 
    if( cgid.isHB() )
    {
       m_hbCellVec[ din ] = IdealObliquePrism( f1, cornersMgr(), parm ) ;
    }
    else
    {
       if( cgid.isHE() )
       {
      const unsigned int index ( din - m_hbCellVec.size() ) ;
      m_heCellVec[ index ] = IdealObliquePrism( f1, cornersMgr(), parm ) ;
       }
       else
       {
      if( cgid.isHO() )
      {
         const unsigned int index ( din 
                        - m_hbCellVec.size() 
                        - m_heCellVec.size() ) ;
         m_hoCellVec[ index ] = IdealObliquePrism( f1, cornersMgr(), parm ) ;
      }
      else
      {
         const unsigned int index ( din 
                        - m_hbCellVec.size() 
                        - m_heCellVec.size() 
                        - m_hoCellVec.size() ) ;
         m_hfCellVec[ index ] = IdealZPrism( f1, cornersMgr(), parm ) ;
      }
       }
    }
    m_validIds.push_back( detId ) ;
 }