HcalDDDGeometryLoader Class Reference

#include <HcalDDDGeometryLoader.h>

Public Types
typedef CaloSubdetectorGeometry *	ReturnType

Public Member Functions
	HcalDDDGeometryLoader (const HcalDDDRecConstants *hcons)
	HcalDDDGeometryLoader ()=delete
ReturnType	load (const HcalTopology &topo, DetId::Detector, int)
ReturnType	load (const HcalTopology &topo)
	Load all of HCAL. More...
virtual	~HcalDDDGeometryLoader ()

Private Member Functions
void	fill (HcalSubdetector, HcalDDDGeometry *)
void	makeCell (const HcalDetId &, const HcalCellType &, double, double, HcalDDDGeometry *geom) const

Private Attributes
const HcalDDDRecConstants *	hcalConstants_
bool	isBH_

Detailed Description

Note

The Geometry as loaded from DDD

Author

S. Banerjee

Definition at line 21 of file HcalDDDGeometryLoader.h.

Member Typedef Documentation

typedef CaloSubdetectorGeometry* HcalDDDGeometryLoader::ReturnType

Definition at line 26 of file HcalDDDGeometryLoader.h.

Constructor & Destructor Documentation

HcalDDDGeometryLoader::HcalDDDGeometryLoader ( const HcalDDDRecConstants *  hcons )

explicit

Definition at line 15 of file HcalDDDGeometryLoader.cc.

References hcalConstants_, HcalDDDRecConstants::isBH(), and isBH_.

                                                                             : hcalConstants_(hcons) {
  isBH_ = hcalConstants_->isBH();
}

HcalDDDGeometryLoader::~HcalDDDGeometryLoader ( )

virtual

Definition at line 19 of file HcalDDDGeometryLoader.cc.

{}

HcalDDDGeometryLoader::HcalDDDGeometryLoader ( )

delete

Member Function Documentation

void HcalDDDGeometryLoader::fill ( HcalSubdetector  subdet, HcalDDDGeometry *  geom  )

private

helper functions to make all the ids and cells, and put them into the vectors and mpas passed in.

Definition at line 65 of file HcalDDDGeometryLoader.cc.

References HcalDDDRecConstants::HcalCellTypes(), hcalConstants_, mps_fire::i, gpuClustering::id, HcalDDDGeometry::increaseReserve(), HcalDDDGeometry::insertCell(), makeCell(), and phi.

Referenced by load().

                                                                              {
  // start by making the new HcalDetIds
  std::vector<HcalCellType> hcalCells = hcalConstants_->HcalCellTypes(subdet);
  geom->insertCell(hcalCells);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalDDDGeometryLoader::fill gets " << hcalCells.size() << " cells for subdetector "
                               << subdet;
#endif
  // Make the new HcalDetIds and the cells

  std::vector<HcalDetId> hcalIds;
  for (auto& hcalCell : hcalCells) {
    int etaRing = hcalCell.etaBin();
    int iside = hcalCell.zside();
    int depthBin = hcalCell.depthSegment();
    double dphi = hcalCell.phiBinWidth();
    std::vector<std::pair<int, double> > phis = hcalCell.phis();
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HCalGeom") << "HcalDDDGeometryLoader: Subdet " << subdet << " side " << iside << " eta "
                                 << etaRing << " depth " << depthBin << " with " << phis.size() << "modules:";
    size_t i(0);
#endif
    geom->increaseReserve(phis.size());
    for (auto& phi : phis) {
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HCalGeom") << "HcalDDDGeometryLoader::fill Cell " << i << " eta " << iside * etaRing << " phi "
                                   << phi.first << "(" << phi.second / CLHEP::deg << ", " << dphi / CLHEP::deg
                                   << ") depth " << depthBin;
      i++;
#endif
      HcalDetId id(subdet, iside * etaRing, phi.first, depthBin);
      hcalIds.emplace_back(id);
      makeCell(id, hcalCell, phi.second, dphi, geom);
    }
  }

  edm::LogInfo("HCalGeom") << "Number of HCAL DetIds made for " << subdet << " is " << hcalIds.size();
}

HcalDDDGeometryLoader::ReturnType HcalDDDGeometryLoader::load	(	const HcalTopology &	topo,
		DetId::Detector	det,
		int	subdet
	)

Definition at line 21 of file HcalDDDGeometryLoader.cc.

References CaloSubdetectorGeometry::allocateCorners(), CaloSubdetectorGeometry::allocatePar(), CaloSubdetectorGeometry::cornersMgr(), submitPVResolutionJobs::count, fill(), relativeConstraints::geom, HcalBarrel, hcalConstants_, HcalEndcap, HcalForward, HcalOuter, HcalDDDRecConstants::numberOfCells(), CaloSubdetectorGeometry::parMgr(), and HcalDDDGeometry::sortValidIds().

Referenced by HcalDDDGeometryEP::produceAligned().

                                                                          {
  HcalSubdetector hsub = static_cast<HcalSubdetector>(subdet);

  HcalDDDGeometry* geom(new HcalDDDGeometry(topo));

  if (geom->cornersMgr() == nullptr) {
    const unsigned int count(hcalConstants_->numberOfCells(HcalBarrel) + hcalConstants_->numberOfCells(HcalEndcap) +
                             2 * hcalConstants_->numberOfCells(HcalForward) + hcalConstants_->numberOfCells(HcalOuter));
    geom->allocateCorners(count);
  }

  //  if( geom->cornersMgr() == 0 )  geom->allocateCorners( 2592 ) ;

  if (geom->parMgr() == nullptr)
    geom->allocatePar(500, 3);

  fill(hsub, geom);
  //fast insertion of valid ids requires sort at end
  geom->sortValidIds();
  return geom;
}

HcalDDDGeometryLoader::ReturnType HcalDDDGeometryLoader::load

(

const HcalTopology &

topo

)

Load all of HCAL.

Definition at line 45 of file HcalDDDGeometryLoader.cc.

References CaloSubdetectorGeometry::allocateCorners(), CaloSubdetectorGeometry::allocatePar(), CaloSubdetectorGeometry::cornersMgr(), submitPVResolutionJobs::count, fill(), relativeConstraints::geom, HcalBarrel, hcalConstants_, HcalEndcap, HcalForward, HcalOuter, HcalDDDRecConstants::numberOfCells(), CaloSubdetectorGeometry::parMgr(), and HcalDDDGeometry::sortValidIds().

                                                                                    {
  HcalDDDGeometry* geom(new HcalDDDGeometry(topo));

  if (geom->cornersMgr() == nullptr) {
    const unsigned int count(hcalConstants_->numberOfCells(HcalBarrel) + hcalConstants_->numberOfCells(HcalEndcap) +
                             2 * hcalConstants_->numberOfCells(HcalForward) + hcalConstants_->numberOfCells(HcalOuter));
    geom->allocateCorners(count);
  }
  if (geom->parMgr() == nullptr)
    geom->allocatePar(500, 3);

  fill(HcalBarrel, geom);
  fill(HcalEndcap, geom);
  fill(HcalForward, geom);
  fill(HcalOuter, geom);
  //fast insertion of valid ids requires sort at end
  geom->sortValidIds();
  return geom;
}

void HcalDDDGeometryLoader::makeCell ( const HcalDetId &  detId, const HcalCellType &  hcalCell, double  phi, double  dphi, HcalDDDGeometry *  geom  ) const

private

Definition at line 104 of file HcalDDDGeometryLoader.cc.

References funct::abs(), funct::cos(), HcalCellType::depthMax(), HcalCellType::depthMin(), HcalCellType::depthType(), PVValHelper::eta, HLT_FULL_cff::eta1, HLT_FULL_cff::eta2, HcalCellType::etaMax(), HcalCellType::etaMin(), funct::exp(), CaloCellGeometry::getParmPtr(), HcalCellType::halfSize(), HcalBarrel, HcalForward, HcalOuter, PixelPluginsPhase0_cfi::isBarrel, isBH_, HcalDDDGeometry::newCellFast(), CaloSubdetectorGeometry::parMgr(), CaloSubdetectorGeometry::parVecVec(), point, alignCSCRings::r, diffTwoXMLs::r1, diffTwoXMLs::r2, jetcorrextractor::sign(), funct::sin(), HcalDetId::subdet(), funct::tan(), theta(), TrackerMaterial_cfi::thickness, x, y, z, and HcalDetId::zside().

Referenced by fill().

                                                                                                                {
  // the two eta boundaries of the cell
  double eta1 = hcalCell.etaMin();
  double eta2 = hcalCell.etaMax();
  HcalSubdetector subdet = detId.subdet();
  double eta = 0.5 * (eta1 + eta2) * detId.zside();
  double deta = (eta2 - eta1);
  double theta = 2.0 * atan(exp(-eta));

  // barrel vs forward
  bool rzType = hcalCell.depthType();
  bool isBarrel = (subdet == HcalBarrel || subdet == HcalOuter);

  double z, r, thickness;
#ifdef EDM_ML_DEBUG
  double r0, r1, r2;
#endif
  if (rzType) {
    r = hcalCell.depthMin();
    if (isBarrel) {
      z = r * sinh(eta);                                  // sinh(eta) == 1/tan(theta)
      thickness = (hcalCell.depthMax() - r) * cosh(eta);  // cosh(eta) == 1/sin(theta)
#ifdef EDM_ML_DEBUG
      r1 = r;
      r2 = hcalCell.depthMax();
      r0 = 0.5 * (r1 + r2);
#endif
    } else {
      z = r * sinh(eta2);
      thickness = 2. * hcalCell.halfSize();
      r = z / sinh(std::abs(eta));
#ifdef EDM_ML_DEBUG
      r0 = z / sinh(std::abs(eta));
      r1 = z / sinh(std::abs(eta) + 0.5 * deta);
      r2 = z / sinh(std::abs(eta) - 0.5 * deta);
#endif
    }
#ifdef EDM_ML_DEBUG
    edm::LogVerbatim("HCalGeom") << "HcalDDDGeometryLoader::makeCell SubDet " << subdet << " eta = " << eta
                                 << " theta = " << theta << " r = " << r << " thickness = " << thickness << " r0-r2 ("
                                 << r0 << ":" << r1 << ":" << r2 << ")";
#endif
  } else {
    z = hcalCell.depthMin();
    thickness = hcalCell.depthMax() - z;
    if (isBH_)
      z += (0.5 * thickness);
    z *= detId.zside();  // get the sign right.
    r = z * tan(theta);
    thickness /= std::abs(cos(theta));
#ifdef EDM_ML_DEBUG
    r0 = z / sinh(std::abs(eta));
    r1 = z / sinh(std::abs(eta) + 0.5 * deta);
    r2 = z / sinh(std::abs(eta) - 0.5 * deta);
    edm::LogVerbatim("HCalGeom") << "HcalDDDGeometryLoader::makeCell SubDet " << subdet << " eta = " << eta
                                 << " theta = " << theta << " z = " << z << " r = " << r << " thickness = " << thickness
                                 << " r0-r2 (" << r0 << ":" << r1 << ":" << r2 << ")";
#endif
  }

  double x = r * cos(phi);
  double y = r * sin(phi);
  GlobalPoint point(x, y, z);

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalDDDGeometryLoader::makeCell for " << detId << " Point " << point
                               << " deta = " << deta << " dphi = " << dphi << " thickness = " << thickness
                               << " isBarrel = " << isBarrel << " " << rzType;
#endif

  std::vector<CCGFloat> hp;
  hp.reserve(3);

  if (subdet == HcalForward) {
    hp.emplace_back(deta / 2.);
    hp.emplace_back(dphi / 2.);
    hp.emplace_back(thickness / 2.);
  } else {
    const double sign(isBarrel ? 1 : -1);
    hp.emplace_back(deta / 2.);
    hp.emplace_back(dphi / 2.);
    hp.emplace_back(sign * thickness / 2.);
  }
  geom->newCellFast(point, point, point, CaloCellGeometry::getParmPtr(hp, geom->parMgr(), geom->parVecVec()), detId);
}

Member Data Documentation

const HcalDDDRecConstants* HcalDDDGeometryLoader::hcalConstants_

private

Definition at line 40 of file HcalDDDGeometryLoader.h.

Referenced by fill(), HcalDDDGeometryLoader(), and load().

bool HcalDDDGeometryLoader::isBH_

private

Definition at line 42 of file HcalDDDGeometryLoader.h.

Referenced by HcalDDDGeometryLoader(), and makeCell().