HcalNumberingFromDDD Class Reference

#include <HcalNumberingFromDDD.h>

Classes
struct	HcalID

Public Member Functions
	HcalNumberingFromDDD (HcalDDDSimConstants *hcons)
HcalID	unitID (int det, const CLHEP::Hep3Vector &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
HcalDDDSimConstants *	hcalConstants

Detailed Description

Definition at line 16 of file HcalNumberingFromDDD.h.

Constructor & Destructor Documentation

HcalNumberingFromDDD::HcalNumberingFromDDD

(

HcalDDDSimConstants *

hcons

)

Definition at line 17 of file HcalNumberingFromDDD.cc.

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

HcalNumberingFromDDD::~HcalNumberingFromDDD

(

)

Definition at line 24 of file HcalNumberingFromDDD.cc.

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

Member Function Documentation

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

Definition at line 30 of file HcalNumberingFromDDD.cc.

References gather_cfg::cout, HcalDDDSimConstants::getEtaHO(), HcalBarrel, hcalConstants, HcalEndcap, HcalForward, cmsBatch::log, SiStripPI::max, min(), funct::sin(), mathSSE::sqrt(), and funct::tan().

Referenced by HcalTestAnalysis::fill(), SimG4HcalValidation::fill(), HcalNumberingFromDDD::HcalID::HcalID(), HcalTestAnalysis::qieAnalysis(), unitID(), and SimG4HGCalValidation::update().

                                             {


  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
  std::cout << "HcalNumberingFromDDD: point = " << point
        << " det " << det << ":" << hsubdet << " eta/R " 
        << etaR << " phi " << hphi << " depth " << depth
        << " layer " << lay << std::endl;
#endif
  return unitID(hsubdet,etaR,hphi,depth,lay);
}

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

Definition at line 72 of file HcalNumberingFromDDD.cc.

References HcalDDDSimConstants::getDetEta(), hcalConstants, and unitID().

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

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

Definition at line 80 of file HcalNumberingFromDDD.cc.

References gather_cfg::cout, HcalDDDSimConstants::getEta(), HcalDDDSimConstants::getPhiCons(), hcalConstants, createfilelist::int, 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((CLHEP::twopi+0.1*ficons.second)/ficons.second);
  int    zside = etaR>0 ? 1: 0;
  double hphi  = phi+ficons.first;
  if (hphi < 0)    hphi += CLHEP::twopi;
  int    iphi  = int(hphi/ficons.second) + 1;
  if (iphi > nphi) iphi = 1;

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

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

Definition at line 105 of file HcalNumberingFromDDD.cc.

References gather_cfg::cout, HcalNumberingFromDDD::HcalID::depth, HcalNumberingFromDDD::HcalID::etaR, HcalDDDSimConstants::getEtaDepth(), HcalBarrel, hcalConstants, HcalOuter, HcalNumberingFromDDD::HcalID::lay, HcalNumberingFromDDD::HcalID::phi, HcalDDDSimConstants::phiNumber(), HcalNumberingFromDDD::HcalID::subdet, tmp, HcalDDDSimConstants::unitPhi(), units(), and HcalNumberingFromDDD::HcalID::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);

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

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

#ifdef EDM_ML_DEBUG
  std::cout << "HcalNumberingFromDDD: det = " << det << " " 
        << tmp.subdet << " zside = " << tmp.zside 
        << " depth = " << tmp.depth << " eta/R = " 
        << tmp.etaR << " phi = "   << tmp.phi << " layer = "
        << tmp.lay << std::endl;
#endif
  return tmp;
}

Member Data Documentation

HcalDDDSimConstants* HcalNumberingFromDDD::hcalConstants

private

Definition at line 38 of file HcalNumberingFromDDD.h.

Referenced by unitID().