#include <HcalDDDGeometry.h>

Inheritance diagram for HcalDDDGeometry:

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	~HcalDDDGeometry ()
	The HcalDDDGeometry will delete all its cell geometries at destruction time. More...

Public Member Functions inherited from CaloSubdetectorGeometry
void	addCell (const DetId &id, CaloCellGeometry *ccg)
	Add a cell to the geometry. More...

void	allocateCorners (CaloCellGeometry::CornersVec::size_type n)

void	allocatePar (ParVec::size_type n, unsigned int m)

	CaloSubdetectorGeometry ()

const CellCont &	cellGeometries () const
	the cells More...

CaloCellGeometry::CornersMgr *	cornersMgr ()

double	deltaEta (const DetId &detId) const

double	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...

Private Attributes
double	deg

double	etaMax_

double	firstHFQuadRing_

std::vector< HcalCellType >	hcalCells_

DetId::Detector	lastReqDet_

int	lastReqSubdet_

std::vector< DetId >	m_validIds

double	twopi

Additional Inherited Members
Public Types inherited from CaloSubdetectorGeometry
typedef std::vector< const CaloCellGeometry * >	CellCont

typedef std::set< DetId >	DetIdSet

typedef std::vector< double >	DimVec

typedef std::vector< unsigned int >	IVec

typedef CaloCellGeometry::ParMgr	ParMgr

typedef CaloCellGeometry::ParVec	ParVec

typedef CaloCellGeometry::ParVecVec	ParVecVec

typedef std::vector< double >	TrVec

Static Protected Member Functions inherited from CaloSubdetectorGeometry
static double	deltaR (const GlobalPoint &p1, const GlobalPoint &p2)

Protected Attributes inherited from CaloSubdetectorGeometry
ParVecVec	m_parVecVec

Detailed Description

Definition at line 10 of file HcalDDDGeometry.h.

Constructor & Destructor Documentation

HcalDDDGeometry::HcalDDDGeometry ( )

explicit

Definition at line 10 of file HcalDDDGeometry.cc.

References deg, M_PI, and twopi.

                                  : lastReqDet_(DetId::Detector(0)), 
                      lastReqSubdet_(0), etaMax_(0),
                      firstHFQuadRing_(40) {
   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 18 of file HcalDDDGeometry.cc.

18 {}

Member Function Documentation

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

virtual

Reimplemented from CaloSubdetectorGeometry.

Definition at line 55 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 21 of file HcalDDDGeometry.cc.

References CaloSubdetectorGeometry::getValidDetIds(), i, 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 118 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;
 }