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) {
  edm::LogInfo("HCalGeom") << "Creating HcalNumberingFromDDD";
}

HcalNumberingFromDDD::~HcalNumberingFromDDD

(

)

Definition at line 22 of file HcalNumberingFromDDD.cc.

                                            {
  edm::LogInfo("HCalGeom") << "Deleting HcalNumberingFromDDD";
}

Member Function Documentation

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

Definition at line 26 of file HcalNumberingFromDDD.cc.

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

Referenced by HcalTestAnalysis::fill(), SimG4HcalValidation::fill(), HCalSD::getHitFibreBundle(), HCalSD::getHitPMT(), HCalSD::layerWeight(), HcalTestAnalysis::qieAnalysis(), HCalSD::setDetUnitId(), 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 DebugLog
  LogDebug("HCalGeom") << "HcalNumberingFromDDD: point = " << point << " det "
               << det << ":" << hsubdet << " eta/R " << etaR << " phi "
               << hphi;
#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 67 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 75 of file HcalNumberingFromDDD.cc.

References HcalDDDSimConstants::getEta(), HcalDDDSimConstants::getPhiCons(), hcalConstants, LogDebug, 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 DebugLog
  LogDebug("HCalGeom") << "HcalNumberingFromDDD: etaR = " << etaR << " : "
                       << zside << "/" << ieta << " phi " << hphi << " : "
                       << iphi;
#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 100 of file HcalNumberingFromDDD.cc.

References HcalNumberingFromDDD::HcalID::depth, HcalNumberingFromDDD::HcalID::etaR, HcalDDDSimConstants::getEtaDepth(), HcalBarrel, hcalConstants, HcalOuter, HcalNumberingFromDDD::HcalID::lay, LogDebug, HcalNumberingFromDDD::HcalID::phi, HcalDDDSimConstants::phiNumber(), HcalNumberingFromDDD::HcalID::subdet, tmp, HcalDDDSimConstants::unitPhi(), units(), and HcalNumberingFromDDD::HcalID::zside.

                                             {


  std::pair<int,int> etaDepth = hcalConstants->getEtaDepth(det, etaR, phi, depth, lay);
  if (det == static_cast<int>(HcalBarrel) && etaDepth.second == 4) {
    det = static_cast<int>(HcalOuter);
  }

  int units     = hcalConstants->unitPhi(det, etaDepth.first);
  int iphi_skip = hcalConstants->phiNumber(phi, units);

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

#ifdef DebugLog
  LogDebug("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

HcalDDDSimConstants* HcalNumberingFromDDD::hcalConstants

private

Definition at line 38 of file HcalNumberingFromDDD.h.

Referenced by unitID().