---

HcalGeometry Class Reference

#include <Geometry/HcalTowerAlgo/interface/HcalGeometry.h>

Inheritance diagram for HcalGeometry:

List of all members.

Public Member Functions
virtual DetId	getClosestCell (const GlobalPoint &r) const
virtual std::vector< DetId > const &	getValidDetIds (DetId::Detector det, int subdet) const
	Get a list of valid detector ids (for the given subdetector).
	HcalGeometry (const HcalTopology *topology)
virtual	~HcalGeometry ()
	The HcalGeometry will delete all its cell geometries at destruction time.
Private Member Functions
int	etaRing (HcalSubdetector bc, double abseta) const
	helper methods for getClosestCell
int	phiBin (double phi, int etaring) const
Private Attributes
DetId::Detector	lastReqDet_
int	lastReqSubdet_
const HcalTopology *	theTopology

---

Detailed Description

Definition at line 8 of file HcalGeometry.h.

---

Constructor & Destructor Documentation

HcalGeometry::HcalGeometry

(

const HcalTopology *

topology

)

[explicit]

Definition at line 5 of file HcalGeometry.cc.

: theTopology(topology),
  lastReqDet_(DetId::Detector(0)), 
  lastReqSubdet_(0) 
{
}

HcalGeometry::~HcalGeometry

(

)

[virtual]

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

Definition at line 13 of file HcalGeometry.cc.

{
}

---

Member Function Documentation

int HcalGeometry::etaRing	(	HcalSubdetector	bc,
		double	abseta
	)			const [private]

helper methods for getClosestCell

Definition at line 109 of file HcalGeometry.cc.

References HcalTopology::firstHFRing(), HcalForward, HcalTopology::lastHERing(), HcalTopology::lastHFRing(), theHBHEEtaBounds, theHFEtaBounds, and theTopology.

Referenced by getClosestCell().

{
  int etaring;
  if( bc == HcalForward ) {
    for(etaring = theTopology->firstHFRing();
        etaring <= theTopology->lastHFRing(); ++etaring)
    {
      if(theHFEtaBounds[etaring-theTopology->firstHFRing()+1] > abseta) break;
    }
  }
  else
  {
    for(etaring = 1;
        etaring <= theTopology->lastHERing(); ++etaring)
    {
      if(theHBHEEtaBounds[etaring] >= abseta) break;
    }
  }

  return etaring;
}

DetId HcalGeometry::getClosestCell

(

const GlobalPoint &

r

)

const [virtual]

Reimplemented from CaloSubdetectorGeometry.

Definition at line 39 of file HcalGeometry.cc.

References PV3DBase< T, PVType, FrameType >::eta(), etaRing(), CaloSubdetectorGeometry::getGeometry(), CaloCellGeometry::getPosition(), HcalBarrel, HcalEmpty, HcalEndcap, HcalForward, HcalTopology::incrementDepth(), HcalTopology::lastHBRing(), HcalTopology::lastHERing(), HcalTopology::lastHFRing(), PV3DBase< T, PVType, FrameType >::mag(), PV3DBase< T, PVType, FrameType >::phi(), phiBin(), funct::pow(), radius(), funct::sqrt(), theHBHEEtaBounds, theTopology, PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), 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());
  
  // figure out subdetector, giving preference to HE in HE/HF overlap region
  HcalSubdetector bc= HcalEmpty;
  if( abseta <= theHBHEEtaBounds[theTopology->lastHBRing()] )
  {
    bc = HcalBarrel;
  }
  else if( abseta <= theHBHEEtaBounds[theTopology->lastHERing()] ) 
  {
    bc = HcalEndcap;
  }
  else
  {
    bc = HcalForward;
  }

  if (bc == HcalForward) {
    static const double z_long=1100.0;
    static const double z_short=1120.0;
    // determine front-face eta
    double radius=sqrt(pow(r.x(),2)+pow(r.y(),2));
    double trueAeta=asinh(z_long/radius);
    // find eta bin
    int etaring = etaRing(bc, trueAeta);
    if (etaring>theTopology->lastHFRing()) etaring=theTopology->lastHFRing(); 
  
    int phibin = phiBin(r.phi(), etaring);

    // add a sign to the etaring
    int etabin = (r.z() > 0) ? etaring : -etaring;
    HcalDetId bestId(bc,etabin,phibin,((fabs(r.z())>=z_short)?(2):(1)));
    return bestId;
  } else {

    // find eta bin
    int etaring = etaRing(bc, abseta);
    
    int phibin = phiBin(r.phi(), etaring);
    
    // add a sign to the etaring
    int etabin = (r.z() > 0) ? etaring : -etaring;
    
    //Now do depth if required
    int dbin = 1;
    double pointradius=r.mag();
    double dradius=99999.;
    HcalDetId currentId(bc, etabin, phibin, dbin);
    HcalDetId bestId;
    for(  ; currentId != HcalDetId(); theTopology->incrementDepth(currentId))
      {    
        const CaloCellGeometry * cell = getGeometry(currentId);
        assert(cell != 0);
        double radius=cell->getPosition().mag();
        if(fabs(pointradius-radius)<dradius) 
          {
            bestId = currentId;
            dradius=fabs(pointradius-radius);
          }
      }
    
    return bestId;
  }
}

std::vector< DetId > const & HcalGeometry::getValidDetIds	(	DetId::Detector	det,
		int	subdet
	)			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 18 of file HcalGeometry.cc.

References begin, CaloSubdetectorGeometry::cellGeometries(), i, id, lastReqDet_, lastReqSubdet_, CaloSubdetectorGeometry::m_validIds, size, and python::multivaluedict::sort().

                                                                                         {
  if (lastReqDet_!=det || lastReqSubdet_!=subdet) {
    lastReqDet_=det;
    lastReqSubdet_=subdet;
    m_validIds.clear();
  }
  if (m_validIds.empty()) {
    m_validIds.reserve(cellGeometries().size());
    CaloSubdetectorGeometry::CellCont::const_iterator i;
    for (i=cellGeometries().begin(); i!=cellGeometries().end(); i++) {
      DetId id(i->first);
      if (id.det()==det && id.subdetId()==subdet) 
        m_validIds.push_back(id);
    }
      std::sort(m_validIds.begin(),m_validIds.end());
  }

  return m_validIds;
}

int HcalGeometry::phiBin	(	double	phi,
		int	etaring
	)			const [private]

Definition at line 132 of file HcalGeometry.cc.

References HcalTopology::firstHFQuadPhiRing(), HcalTopology::nPhiBins(), and theTopology.

Referenced by getClosestCell().

{
  double twopi = M_PI+M_PI;
  //put phi in correct range (0->2pi)
  if(phi<0.0) phi += twopi;
  int nphibins = theTopology->nPhiBins(etaring);
  int phibin= static_cast<int>(phi/twopi*nphibins)+1;
  int iphi;

  // rings 40 and 41 are offset wrt the other phi numbering
  //  1        1         1         2
  //  ------------------------------
  //  72       36        36        1
  if(etaring >= theTopology->firstHFQuadPhiRing())
  {
    phi+=(twopi/36); //shift by half tower.    
    phibin=static_cast<int>(phi/twopi*nphibins);
    if (phibin==0) phibin=18;
    iphi=phibin*4-1; // 71,3,5,
  } else {
    // convert to the convention of numbering 1,3,5, in 36 phi bins
    iphi=(phibin-1)*(72/nphibins) + 1;
  }

  return iphi;
}

---

Member Data Documentation

DetId::Detector HcalGeometry::lastReqDet_ [mutable, private]

Definition at line 24 of file HcalGeometry.h.

Referenced by getValidDetIds().

int HcalGeometry::lastReqSubdet_ [mutable, private]

Definition at line 25 of file HcalGeometry.h.

Referenced by getValidDetIds().

const HcalTopology* HcalGeometry::theTopology [private]

Definition at line 23 of file HcalGeometry.h.

Referenced by etaRing(), getClosestCell(), and phiBin().

---

The documentation for this class was generated from the following files:

• Geometry/HcalTowerAlgo/interface/HcalGeometry.h
• Geometry/HcalTowerAlgo/src/HcalGeometry.cc

---