---

HcalDDDGeometry Class Reference

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

Inheritance diagram for HcalDDDGeometry:

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).
	HcalDDDGeometry ()
int	insertCell (std::vector< HcalCellType::HcalCellType > const &)
virtual	~HcalDDDGeometry ()
	The HcalDDDGeometry will delete all its cell geometries at destruction time.
Private Attributes
double	deg
double	etaMax_
double	firstHFQuadRing_
std::vector < HcalCellType::HcalCellType >	hcalCells_
DetId::Detector	lastReqDet_
int	lastReqSubdet_
double	twopi

---

Detailed Description

Definition at line 10 of file HcalDDDGeometry.h.

---

Constructor & Destructor Documentation

HcalDDDGeometry::HcalDDDGeometry

(

)

[explicit]

Definition at line 6 of file HcalDDDGeometry.cc.

References deg, 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 14 of file HcalDDDGeometry.cc.

{}

---

Member Function Documentation

DetId HcalDDDGeometry::getClosestCell

(

const GlobalPoint &

r

)

const [virtual]

Reimplemented from CaloSubdetectorGeometry.

Definition at line 44 of file HcalDDDGeometry.cc.

References deg, PV3DBase< T, PVType, FrameType >::eta(), etaMax_, gen_hydjet_hard_muon_trigger_cfg::etaMin, firstHFQuadRing_, hcalCells_, HcalForward, i, LogDebug, PV3DBase< T, PVType, FrameType >::mag(), PV3DBase< T, PVType, FrameType >::phi(), phi, radius(), twopi, PV3DBase< T, PVType, FrameType >::z(), and 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());
  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  = static_cast<int>(((phi/deg)+hcalCells_[i].phiOffset())/
                                       hcalCells_[i].phiBinWidth()) + 1;
        // rings 40 and 41 are offset wrt the other phi numbering
        //  1        1         1         2
        //  ------------------------------
        //  72       36        36        1
        if (etaring >= firstHFQuadRing_) {
          ++phibin;
          if (phibin > hcalCells_[i].nPhiBins()) 
            phibin -= hcalCells_[i].nPhiBins();
        }
 
        // convert to the convention of numbering 1,3,5, in 36 phi bins
        // and 1,5,9 in 18 phi bins
        phibin     = (phibin-1)*(hcalCells_[i].nPhiModule()) + 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,
		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 17 of file HcalDDDGeometry.cc.

References begin, CaloSubdetectorGeometry::cellGeometries(), i, id, lastReqDet_, lastReqSubdet_, LogDebug, 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());
  }
  
  LogDebug("HCalGeom") << "HcalDDDGeometry::getValidDetIds: "
                       << m_validIds.size() << " valid IDs found for detector "
                       << det << " Sub-detector " << subdet;
  return m_validIds;
}

int HcalDDDGeometry::insertCell

(

std::vector< HcalCellType::HcalCellType > const &

cells

)

Definition at line 101 of file HcalDDDGeometry.cc.

References gen_hydjet_hard_muon_trigger_cfg::etaMax, etaMax_, hcalCells_, i, LogDebug, and funct::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;
}

---

Member Data Documentation

double HcalDDDGeometry::deg [private]

Definition at line 28 of file HcalDDDGeometry.h.

Referenced by getClosestCell(), and HcalDDDGeometry().

double HcalDDDGeometry::etaMax_ [private]

Definition at line 29 of file HcalDDDGeometry.h.

Referenced by getClosestCell(), and insertCell().

double HcalDDDGeometry::firstHFQuadRing_ [private]

Definition at line 29 of file HcalDDDGeometry.h.

Referenced by getClosestCell().

std::vector<HcalCellType::HcalCellType> HcalDDDGeometry::hcalCells_ [private]

Definition at line 24 of file HcalDDDGeometry.h.

Referenced by getClosestCell(), and insertCell().

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

Definition at line 25 of file HcalDDDGeometry.h.

Referenced by getValidDetIds().

int HcalDDDGeometry::lastReqSubdet_ [mutable, private]

Definition at line 26 of file HcalDDDGeometry.h.

Referenced by getValidDetIds().

double HcalDDDGeometry::twopi [private]

Definition at line 28 of file HcalDDDGeometry.h.

Referenced by getClosestCell(), and HcalDDDGeometry().

---

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

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

---