HcalNumberingFromDDD Class Reference

#include <HcalNumberingFromDDD.h>

Classes
struct	HcalID

Public Member Functions
	HcalNumberingFromDDD (const HcalDDDSimConstants *hcons)
HcalID	unitID (int det, const math::XYZVectorD &pos, int depth, int lay=-1) const
HcalID	unitID (double eta, double phi, int depth=1, int lay=-1) const
HcalID	unitID (int det, double etaR, double phi, int depth, int lay=-1) const
HcalID	unitID (int det, int zside, int depth, int etaR, int phi, int lay=-1) const
	~HcalNumberingFromDDD ()

Private Attributes
const HcalDDDSimConstants *	hcalConstants

Detailed Description

Definition at line 16 of file HcalNumberingFromDDD.h.

Constructor & Destructor Documentation

HcalNumberingFromDDD()

HcalNumberingFromDDD::HcalNumberingFromDDD

(

const HcalDDDSimConstants *

hcons

)

Definition at line 18 of file HcalNumberingFromDDD.cc.

                                                                           : hcalConstants(hcons) {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "Creating HcalNumberingFromDDD";
#endif
}

~HcalNumberingFromDDD()

HcalNumberingFromDDD::~HcalNumberingFromDDD

(

)

Definition at line 24 of file HcalNumberingFromDDD.cc.

                                            {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "Deleting HcalNumberingFromDDD";
#endif
}

Member Function Documentation

unitID() [1/4]

HcalNumberingFromDDD::HcalID HcalNumberingFromDDD::unitID	(	int	det,
		const math::XYZVectorD &	pos,
		int	depth,
		int	lay = -1
	)		const

Definition at line 30 of file HcalNumberingFromDDD.cc.

References Matriplex::atan2(), hcalRecHitTable_cff::depth, HcalDDDSimConstants::getEtaHO(), HcalBarrel, hcalConstants, HcalEndcap, HcalForward, CrabHelper::log, WZElectronSkims53X_cff::max, SiStripPI::min, point, funct::sin(), mathSSE::sqrt(), and funct::tan().

Referenced by HcalDumpGeometry::dumpTouch(), and unitID().

                                                                         {
  double hx = point.x();
  double hy = point.y();
  double hz = point.z();
  double hR = sqrt(hx * hx + hy * hy + hz * hz);
  double htheta = (hR == 0. ? 0. : acos(std::max(std::min(hz / hR, 1.0), -1.0)));
  double hsintheta = sin(htheta);
  double hphi = (hR * hsintheta == 0. ? 0. : atan2(hy, hx));
  double heta = (fabs(hsintheta) == 1. ? 0. : -log(fabs(tan(htheta / 2.))));

  int hsubdet = 0;
  double etaR;

  //First eta index
  if (det == 5) {  // Forward HCal
    hsubdet = static_cast<int>(HcalForward);
    hR = sqrt(hx * hx + hy * hy);
    etaR = (heta >= 0. ? hR : -hR);
  } else {  // Barrel or Endcap
    etaR = heta;
    if (det == 3) {
      hsubdet = static_cast<int>(HcalBarrel);
      etaR = hcalConstants->getEtaHO(heta, hx, hy, hz);
    } else {
      hsubdet = static_cast<int>(HcalEndcap);
    }
  }

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalNumberingFromDDD: point = " << point << " det " << det << ":" << hsubdet
                               << " eta/R " << etaR << " phi " << hphi << " depth " << depth << " layer " << lay;
#endif
  return unitID(hsubdet, etaR, hphi, depth, lay);
}

unitID() [2/4]

HcalNumberingFromDDD::HcalID HcalNumberingFromDDD::unitID	(	double	eta,
		double	phi,
		int	depth = 1,
		int	lay = -1
	)		const

Definition at line 68 of file HcalNumberingFromDDD.cc.

References hcalRecHitTable_cff::depth, PVValHelper::eta, HcalDDDSimConstants::getDetEta(), hcalConstants, and unitID().

                                                                                                       {
  std::pair<int, double> detEta = hcalConstants->getDetEta(eta, depth);
  return unitID(detEta.first, detEta.second, fi, depth, lay);
}

unitID() [3/4]

HcalNumberingFromDDD::HcalID HcalNumberingFromDDD::unitID	(	int	det,
		double	etaR,
		double	phi,
		int	depth,
		int	lay = -1
	)		const

Definition at line 73 of file HcalNumberingFromDDD.cc.

References hcalRecHitTable_cff::depth, HcalDDDSimConstants::getEta(), HcalDDDSimConstants::getPhiCons(), hcalConstants, hcalRecHitTable_cff::ieta, createfilelist::int, hcalRecHitTable_cff::iphi, nphi, phi, unitID(), and ecaldqm::zside().

                                                                                                                  {
  double hetaR = fabs(etaR);
  int ieta = hcalConstants->getEta(det, lay, hetaR);
  std::pair<double, double> ficons = hcalConstants->getPhiCons(det, ieta);

  int nphi = int((2._pi + 0.1 * ficons.second) / ficons.second);
  int zside = etaR > 0 ? 1 : 0;
  double hphi = phi + ficons.first;
  if (hphi < 0)
    hphi += (2._pi);
  int iphi = int(hphi / ficons.second) + 1;
  if (iphi > nphi)
    iphi = 1;

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalNumberingFromDDD: etaR = " << etaR << " : " << zside << "/" << ieta << " phi "
                               << hphi << " : " << iphi;
#endif
  return unitID(det, zside, depth, ieta, iphi, lay);
}

unitID() [4/4]

HcalNumberingFromDDD::HcalID HcalNumberingFromDDD::unitID	(	int	det,
		int	zside,
		int	depth,
		int	etaR,
		int	phi,
		int	lay = -1
	)		const

Definition at line 94 of file HcalNumberingFromDDD.cc.

References hcalRecHitTable_cff::depth, HcalDDDSimConstants::getEtaDepth(), HcalBarrel, hcalConstants, HcalOuter, phi, HcalDDDSimConstants::phiNumber(), createJobs::tmp, HcalDDDSimConstants::unitPhi(), units(), and ecaldqm::zside().

                                                                     {
  if (det == static_cast<int>(HcalBarrel) && lay > 17)
    det = static_cast<int>(HcalOuter);

  int units = hcalConstants->unitPhi(det, etaR);
  int iphi_skip = hcalConstants->phiNumber(phi, units);

  if ((lay == 1) && (etaR == 16))
    etaR = 15;

  std::pair<int, int> etaDepth = hcalConstants->getEtaDepth(det, etaR, iphi_skip, zside, depth, lay);

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalNumberingFromDDD: phi units= " << units << "  iphi_skip= " << iphi_skip;
#endif
  HcalNumberingFromDDD::HcalID tmp(det, zside, etaDepth.second, etaDepth.first, phi, iphi_skip, lay);

#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalNumberingFromDDD: det = " << det << " " << tmp.subdet << " zside = " << tmp.zside
                               << " depth = " << tmp.depth << " eta/R = " << tmp.etaR << " phi = " << tmp.phi
                               << " layer = " << tmp.lay;
#endif
  return tmp;
}

Member Data Documentation

hcalConstants

const HcalDDDSimConstants* HcalNumberingFromDDD::hcalConstants

private

Definition at line 33 of file HcalNumberingFromDDD.h.

Referenced by unitID().