#include <RecHitTools.h>

Public Member Functions
unsigned int	getCell (const DetId &) const

float	getEta (const GlobalPoint &position, const float &vertex_z=0.) const

float	getEta (const DetId &id, const float &vertex_z=0.) const

void	getEvent (const edm::Event &)

void	getEventSetup (const edm::EventSetup &)

const CaloGeometry *	getGeometry () const

unsigned int	getLayer (ForwardSubdetector type) const

unsigned int	getLayer (const DetId &) const

unsigned int	getLayerWithOffset (const DetId &) const

float	getPhi (const GlobalPoint &position) const

float	getPhi (const DetId &id) const

GlobalPoint	getPosition (const DetId &id) const

float	getPt (const GlobalPoint &position, const float &hitEnergy, const float &vertex_z=0.) const

float	getPt (const DetId &id, const float &hitEnergy, const float &vertex_z=0.) const

std::float_t	getRadiusToSide (const DetId &) const

std::float_t	getSiThickness (const DetId &) const

unsigned int	getWafer (const DetId &) const

bool	isHalfCell (const DetId &) const

unsigned int	lastLayerEE () const

unsigned int	lastLayerFH () const

unsigned int	maxNumberOfWafersPerLayer () const

	RecHitTools ()

int	zside (const DetId &id) const

	~RecHitTools ()

Private Attributes
unsigned int	bhOffset_

unsigned int	fhOffset_

const CaloGeometry *	geom_

unsigned int	maxNumberOfWafersPerLayer_

Detailed Description

Definition at line 18 of file RecHitTools.h.

Constructor & Destructor Documentation

hgcal::RecHitTools::RecHitTools ( )

inline

Definition at line 20 of file RecHitTools.h.

20 : geom_(nullptr), fhOffset_(0), bhOffset_(0) {}

hgcal::RecHitTools::fhOffset_

unsigned int fhOffset_

Definition: RecHitTools.h:57

hgcal::RecHitTools::bhOffset_

unsigned int bhOffset_

Definition: RecHitTools.h:57

hgcal::RecHitTools::geom_

const CaloGeometry * geom_

Definition: RecHitTools.h:56

hgcal::RecHitTools::~RecHitTools ( )

inline

Definition at line 21 of file RecHitTools.h.

References position, and ecaldqm::zside().

21 {}

Member Function Documentation

unsigned int RecHitTools::getCell ( const DetId & id ) const

Definition at line 201 of file RecHitTools.cc.

References HGCalDetId::cell(), TauDecayModes::dec, DetId::Forward, and SiStripPI::max.

                                                        {
   if( id.det() != DetId::Forward ) {
     edm::LogError("getCell::InvalidSiliconDetid")
       << "det id: " << std::hex << id.rawId() << std::dec << ":" 
       << id.det() << " is not HGCal silicon!";
     return std::numeric_limits<unsigned int>::max();
   }
   const HGCalDetId hid(id);
   unsigned int cell = hid.cell();
   return cell;
 }

float RecHitTools::getEta	(	const GlobalPoint &	position,
		const float &	vertex_z = `0.`
	)		const

Definition at line 225 of file RecHitTools.cc.

References PV3DBase< T, PVType, FrameType >::eta(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by getEta(), getPt(), and HGCalIsoCalculator::setRecHits().

                                                                                   {
   GlobalPoint corrected_position = GlobalPoint(position.x(), position.y(), position.z()-vertex_z);
   return corrected_position.eta();
 }

float RecHitTools::getEta	(	const DetId &	id,
		const float &	vertex_z = `0.`
	)		const

Definition at line 230 of file RecHitTools.cc.

References PVValHelper::eta, getEta(), getPosition(), and position.

                                                                       {
   GlobalPoint position = getPosition(id);
   float eta = getEta(position, vertex_z);
   return eta;
 }

void RecHitTools::getEvent ( const edm::Event & ev )

Definition at line 58 of file RecHitTools.cc.

Referenced by hgcal::ClusterTools::getEvent().

58 {

59 }

void RecHitTools::getEventSetup ( const edm::EventSetup & es )

Definition at line 61 of file RecHitTools.cc.

References bhOffset_, fhOffset_, DetId::Forward, relativeConstraints::geom, geom_, edm::EventSetup::get(), CaloGeometry::getSubdetectorGeometry(), HGCEE, HGCHEF, LayerTriplets::layers(), SiStripPI::max, maxNumberOfWafersPerLayer_, and edm::ESHandle< T >::product().

Referenced by HGCalHitCalibration::analyze(), HGCalShowerSeparation::analyze(), HGCalEgammaIDHelper::eventInit(), hgcal::ClusterTools::getEventSetup(), HGCal3DClustering::getEventSetup(), HGCalImagingAlgo::getEventSetup(), PFHGCalRecHitCreator< DET, Layer, subdet >::importRecHits(), and RealisticSimClusterMapper::update().

                                                        {
 
   edm::ESHandle<CaloGeometry> geom;
   es.get<CaloGeometryRecord>().get(geom);
 
   geom_ = geom.product();
   auto geomEE = dynamic_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward,ForwardSubdetector::HGCEE));
   fhOffset_ = (geomEE->topology().dddConstants()).layers(true);
   unsigned int wmaxEE = 1 + (geomEE->topology().dddConstants()).waferMax();
   auto geomFH = dynamic_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward,ForwardSubdetector::HGCHEF));
   bhOffset_ = fhOffset_ + (geomFH->topology().dddConstants()).layers(true);
   unsigned int wmaxFH = 1 + (geomFH->topology().dddConstants()).waferMax();
   maxNumberOfWafersPerLayer_ = std::max(wmaxEE,wmaxFH);
 
 }

const CaloGeometry* hgcal::RecHitTools::getGeometry ( ) const

inline

Definition at line 51 of file RecHitTools.h.

Referenced by HGCalImagingAlgo::computeThreshold().

51 {return geom_;};

hgcal::RecHitTools::geom_

const CaloGeometry * geom_

Definition: RecHitTools.h:56

unsigned int RecHitTools::getLayer ( ForwardSubdetector type ) const

Definition at line 134 of file RecHitTools.cc.

References DetId::Forward, ForwardEmpty, geom_, CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalEndcap, HGCEE, HGCHEB, HGCHEF, and LayerTriplets::layers().

Referenced by RealisticSimClusterMapper::buildClusters(), and getLayerWithOffset().

                                                                       {
 
   int layer;
   switch (type) {
     case(ForwardSubdetector::HGCEE): {
       auto geomEE = dynamic_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward,ForwardSubdetector::HGCEE));
       layer       = (geomEE->topology().dddConstants()).layers(true);
       break;
     }
     case (ForwardSubdetector::HGCHEF): {
       auto geomFH = dynamic_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward,ForwardSubdetector::HGCHEF));
       layer       = (geomFH->topology().dddConstants()).layers(true);
       break;
     }
     case (ForwardSubdetector::HGCHEB): {
       auto geomBH = dynamic_cast<const HcalGeometry*>(geom_->getSubdetectorGeometry(DetId::Hcal,HcalSubdetector::HcalEndcap));
       layer       = (geomBH->topology().dddConstants())->getMaxDepth(1);
       break;
     }
     case (ForwardSubdetector::ForwardEmpty): {
       auto geomEE = dynamic_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward,ForwardSubdetector::HGCEE));
       layer       = (geomEE->topology().dddConstants()).layers(true);
       auto geomFH = dynamic_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward,ForwardSubdetector::HGCHEF));
       layer      += (geomFH->topology().dddConstants()).layers(true);
       auto geomBH = dynamic_cast<const HcalGeometry*>(geom_->getSubdetectorGeometry(DetId::Hcal,HcalSubdetector::HcalEndcap));
       layer      += (geomBH->topology().dddConstants())->getMaxDepth(1);
       break;
     }
     default: layer = 0;
   }
   return (unsigned int)(layer);
 }

unsigned int RecHitTools::getLayer ( const DetId & id ) const

Definition at line 167 of file RecHitTools.cc.

References HcalDetId::depth(), DetId::Forward, DetId::Hcal, HcalEndcap, HGCalDetId::layer(), and SiStripPI::max.

                                                         {
   unsigned int layer = std::numeric_limits<unsigned int>::max();
   if( id.det() == DetId::Forward) {
     const HGCalDetId hid(id);
     layer = hid.layer();
   } else if( id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
     const HcalDetId hcid(id);
     layer = hcid.depth();
   }
   return layer;
 }

unsigned int RecHitTools::getLayerWithOffset ( const DetId & id ) const

Definition at line 179 of file RecHitTools.cc.

References bhOffset_, fhOffset_, DetId::Forward, getLayer(), DetId::Hcal, HcalEndcap, and HGCHEF.

Referenced by HGCalHitCalibration::analyze(), HGCalShowerSeparation::analyze(), RealisticSimClusterMapper::buildClusters(), HGCalImagingAlgo::computeThreshold(), HGCalHitCalibration::fillWithRecHits(), hgcal::ClusterTools::getLayer(), HGCal3DClustering::organizeByLayer(), HGCalImagingAlgo::populate(), HGCalIsoCalculator::produceHGCalIso(), and hgcal::EGammaPCAHelper::storeRecHits().

                                                                   {
   unsigned int layer = getLayer(id);
   if( id.det() == DetId::Forward && id.subdetId() == HGCHEF ) {
     layer += fhOffset_;
   } else if( id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
     layer += bhOffset_;
   }
   return layer;
 }

float RecHitTools::getPhi ( const GlobalPoint & position ) const

Definition at line 236 of file RecHitTools.cc.

References PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by HGCalIsoCalculator::setRecHits().

                                                            {
   float phi = atan2(position.y(),position.x());
   return phi;
 }

float RecHitTools::getPhi ( const DetId & id ) const

Definition at line 241 of file RecHitTools.cc.

References getPosition(), position, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                {
   GlobalPoint position = getPosition(id);
   float phi = atan2(position.y(),position.x());
   return phi;
 }

GlobalPoint RecHitTools::getPosition ( const DetId & id ) const

Definition at line 77 of file RecHitTools.cc.

References relativeConstraints::geom, geom_, CaloSubdetectorGeometry::getGeometry(), HGCalGeometry::getPosition(), CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, and position.

Referenced by HGCalShowerSeparation::analyze(), RealisticSimClusterMapper::buildClusters(), hgcal::EGammaPCAHelper::findZFirstLayer(), getEta(), getPhi(), getPt(), hgcal::ClusterTools::getWidths(), HGCalImagingAlgo::Hexel::Hexel(), HGCalImagingAlgo::populate(), HGCalIsoCalculator::setRecHits(), and hgcal::EGammaPCAHelper::storeRecHits().

                                                           {
   auto geom = geom_->getSubdetectorGeometry(id);
   check_geom(geom);
   GlobalPoint position;
   if( id.det() == DetId::Hcal ) {
     position = geom->getGeometry(id)->getPosition();
   } else {
     auto hg = dynamic_cast<const HGCalGeometry*>(geom);
     position = hg->getPosition(id);
   }
   return position;
 }

float RecHitTools::getPt	(	const GlobalPoint &	position,
		const float &	hitEnergy,
		const float &	vertex_z = `0.`
	)		const

Definition at line 247 of file RecHitTools.cc.

References PVValHelper::eta, getEta(), and EnergyCorrector::pt.

                                                                                                          {
   float eta = getEta(position, vertex_z);
   float pt = hitEnergy / cosh(eta);
   return pt;
 }

float RecHitTools::getPt	(	const DetId &	id,
		const float &	hitEnergy,
		const float &	vertex_z = `0.`
	)		const

Definition at line 253 of file RecHitTools.cc.

References PVValHelper::eta, getEta(), getPosition(), position, and EnergyCorrector::pt.

                                                                                              {
   GlobalPoint position = getPosition(id);
   float eta = getEta(position, vertex_z);
   float pt = hitEnergy / cosh(eta);
   return pt;
 }

std::float_t RecHitTools::getRadiusToSide ( const DetId & id ) const

Definition at line 119 of file RecHitTools.cc.

References HGCalDDDConstants::cellSizeHex(), TauDecayModes::dec, DetId::Forward, geom_, CaloGeometry::getSubdetectorGeometry(), SiStripPI::max, findQualityFiles::size, and HGCalDetId::waferType().

                                                              {
   auto geom = geom_->getSubdetectorGeometry(id);
   check_geom(geom);
   if( id.det() != DetId::Forward ) {
     edm::LogError("getRadiusToSide::InvalidSiliconDetid")
       << "det id: " << std::hex << id.rawId() << std::dec << ":" 
       << id.det() << " is not HGCal silicon!";
     return std::numeric_limits<std::float_t>::max();
   }
   const HGCalDetId hid(id);
   auto ddd = get_ddd(geom,hid);
   std::float_t size = ddd->cellSizeHex(hid.waferType());
   return size;
 }

std::float_t RecHitTools::getSiThickness ( const DetId & id ) const

Definition at line 103 of file RecHitTools.cc.

References TauDecayModes::dec, DetId::Forward, geom_, CaloGeometry::getSubdetectorGeometry(), LogDebug, HGCalDetId::wafer(), and HGCalDDDConstants::waferTypeL().

Referenced by HGCalHitCalibration::analyze(), HGCalImagingAlgo::computeThreshold(), HGCalHitCalibration::fillWithRecHits(), HGCalImagingAlgo::populate(), and hgcal::EGammaPCAHelper::storeRecHits().

                                                             {
   auto geom = geom_->getSubdetectorGeometry(id);
   check_geom(geom);
   if( id.det() != DetId::Forward ) {
     LogDebug("getSiThickness::InvalidSiliconDetid")
       << "det id: " << std::hex << id.rawId() << std::dec << ":" 
       << id.det() << " is not HGCal silicon!";
     return 0.37;
   }
   const HGCalDetId hid(id);
   auto ddd = get_ddd(geom,hid);
   unsigned int wafer = hid.wafer();
   int tidx = ddd->waferTypeL(wafer);
   return idx_to_thickness*tidx;
 }

unsigned int RecHitTools::getWafer ( const DetId & id ) const

Definition at line 189 of file RecHitTools.cc.

References TauDecayModes::dec, DetId::Forward, SiStripPI::max, and HGCalDetId::wafer().

Referenced by HGCalImagingAlgo::computeThreshold(), and HGCalImagingAlgo::populate().

                                                         {
   if( id.det() != DetId::Forward ) {
     edm::LogError("getWafer::InvalidSiliconDetid")
       << "det id: " << std::hex << id.rawId() << std::dec << ":" 
       << id.det() << " is not HGCal silicon!";
     return std::numeric_limits<unsigned int>::max();
   }
   const HGCalDetId hid(id);
   unsigned int wafer = hid.wafer();
   return wafer;
 }

bool RecHitTools::isHalfCell ( const DetId & id ) const

Definition at line 213 of file RecHitTools.cc.

References HGCalDetId::cell(), DetId::Forward, geom_, CaloGeometry::getSubdetectorGeometry(), HGCalDDDConstants::isHalfCell(), HGCalDetId::waferType(), and HGCalDDDConstants::waferTypeT().

Referenced by HGCalImagingAlgo::populate().

                                                   {
   if( id.det() != DetId::Forward ) {
     return false;
   }
   auto geom = geom_->getSubdetectorGeometry(id);
   check_geom(geom);
   const HGCalDetId hid(id);
   auto ddd = get_ddd(geom,hid);
   const int waferType = ddd->waferTypeT(hid.waferType());
   return ddd->isHalfCell(waferType,hid.cell());
 }

unsigned int hgcal::RecHitTools::lastLayerEE ( ) const

inline

Definition at line 52 of file RecHitTools.h.

Referenced by RealisticSimClusterMapper::buildClusters(), hgcal::EGammaPCAHelper::computePCA(), and hgcal::EGammaPCAHelper::findZFirstLayer().

52 {return fhOffset_;}

hgcal::RecHitTools::fhOffset_

unsigned int fhOffset_

Definition: RecHitTools.h:57

unsigned int hgcal::RecHitTools::lastLayerFH ( ) const

inline

Definition at line 53 of file RecHitTools.h.

Referenced by RealisticSimClusterMapper::buildClusters(), and hgcal::EGammaPCAHelper::findZFirstLayer().

53 {return bhOffset_;}

hgcal::RecHitTools::bhOffset_

unsigned int bhOffset_

Definition: RecHitTools.h:57

unsigned int hgcal::RecHitTools::maxNumberOfWafersPerLayer ( ) const

inline

Definition at line 54 of file RecHitTools.h.

54 {return maxNumberOfWafersPerLayer_;}

hgcal::RecHitTools::maxNumberOfWafersPerLayer_

unsigned int maxNumberOfWafersPerLayer_

Definition: RecHitTools.h:57

int RecHitTools::zside ( const DetId & id ) const

Definition at line 91 of file RecHitTools.cc.

References DetId::Forward, DetId::Hcal, HcalEndcap, HGCalDetId::zside(), and HcalDetId::zside().

Referenced by HGCalImagingAlgo::computeThreshold(), and HGCalImagingAlgo::populate().

                                             {
   int zside = 0;
   if( id.det() == DetId::Forward) {
     const HGCalDetId hid(id);
     zside = hid.zside();
   } else if( id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
     const HcalDetId hcid(id);
     zside = hcid.zside();
   }
   return zside;
 }