#include <RecHitTools.h>

Public Member Functions
std::pair< uint32_t, uint32_t >	firstAndLastLayer (DetId::Detector det, int subdet) const

unsigned int	firstLayerBH () const

std::pair< int, int >	getCell (const DetId &) const

int	getCellType (const DetId &id) const

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

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

const CaloGeometry *	getGeometry () const

int	getGeometryType () const

unsigned int	getLayer (DetId::Detector type, bool nose=false) 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

GlobalPoint	getPositionLayer (int layer, bool nose=false) 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::pair< float, float >	getScintDEtaDPhi (const DetId &) const

int	getScintMaxIphi () const

int	getSiThickIndex (const DetId &) const

std::float_t	getSiThickness (const DetId &) const

const CaloSubdetectorGeometry *	getSubdetectorGeometry (const DetId &id) const

std::pair< int, int >	getWafer (const DetId &) const

bool	isHalfCell (const DetId &) const

bool	isOnlySilicon (const unsigned int layer) const

bool	isScintillator (const DetId &) const

bool	isSilicon (const DetId &) const

unsigned int	lastLayer (bool nose=false) const

unsigned int	lastLayerBH () const

unsigned int	lastLayerEE (bool nose=false) const

unsigned int	lastLayerFH () const

bool	maskCell (const DetId &id, int corners=3) const

unsigned int	maxNumberOfWafersPerLayer (bool nose=false) const

	RecHitTools ()

void	setGeometry (CaloGeometry const &)

int	zside (const DetId &id) const

	~RecHitTools ()

Private Attributes
unsigned int	bhFirstLayer_

unsigned int	bhLastLayer_

int	bhMaxIphi_

unsigned int	bhOffset_

unsigned int	eeOffset_

unsigned int	fhLastLayer_

unsigned int	fhOffset_

const CaloGeometry *	geom_

int	geometryType_

unsigned int	maxNumberOfWafersNose_

unsigned int	maxNumberOfWafersPerLayer_

unsigned int	noseLastLayer_

Detailed Description

Definition at line 23 of file RecHitTools.h.

Constructor & Destructor Documentation

◆ RecHitTools()

hgcal::RecHitTools::RecHitTools ( )

inline

Definition at line 25 of file RecHitTools.h.

         : geom_(nullptr),
           eeOffset_(0),
           fhOffset_(0),
           bhFirstLayer_(0),
           bhOffset_(0),
           fhLastLayer_(0),
           noseLastLayer_(0),
           geometryType_(0) {}

◆ ~RecHitTools()

hgcal::RecHitTools::~RecHitTools ( )

inline

Definition at line 34 of file RecHitTools.h.

34 {}

Member Function Documentation

◆ firstAndLastLayer()

std::pair< uint32_t, uint32_t > RecHitTools::firstAndLastLayer	(	DetId::Detector	det,
		int	subdet
	)		const

Definition at line 521 of file RecHitTools.cc.

References bhFirstLayer_, bhLastLayer_, eeOffset_, fhLastLayer_, fhOffset_, DetId::Forward, HFNose, DetId::HGCalEE, DetId::HGCalHSi, HGCEE, HGCHEF, and noseLastLayer_.

                                                                                               {
   if ((det == DetId::HGCalEE) || ((det == DetId::Forward) && (subdet == HGCEE))) {
     return std::make_pair(eeOffset_ + 1, fhOffset_);
   } else if ((det == DetId::HGCalHSi) || ((det == DetId::Forward) && (subdet == HGCHEF))) {
     return std::make_pair(fhOffset_ + 1, fhLastLayer_);
   } else if ((det == DetId::Forward) && (subdet == HFNose)) {
     return std::make_pair(1, noseLastLayer_);
   } else {
     return std::make_pair(bhFirstLayer_, bhLastLayer_);
   }
 }

◆ firstLayerBH()

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

inline

Definition at line 78 of file RecHitTools.h.

References bhFirstLayer_.

Referenced by HGCalClusteringAlgoBase::getEventSetup(), HGCal3DClustering::makeClusters(), and HGCalDepthPreClusterer::makePreClusters().

78 { return bhFirstLayer_; }

hgcal::RecHitTools::bhFirstLayer_

unsigned int bhFirstLayer_

Definition: RecHitTools.h:91

◆ getCell()

std::pair< int, int > RecHitTools::getCell ( const DetId & id ) const

Definition at line 408 of file RecHitTools.cc.

References HGCalDetId::cell(), HFNoseDetId::cellU(), HGCSiliconDetId::cellU(), HFNoseDetId::cellV(), HGCSiliconDetId::cellV(), TauDecayModes::dec, DetId::Forward, HFNose, DetId::HGCalEE, DetId::HGCalHSi, and SiStripPI::max.

                                                             {
   int cellU = std::numeric_limits<int>::max();
   int cellV = 0;
   if ((id.det() == DetId::HGCalEE) || (id.det() == DetId::HGCalHSi)) {
     cellU = HGCSiliconDetId(id).cellU();
     cellV = HGCSiliconDetId(id).cellV();
   } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
     cellU = HFNoseDetId(id).cellU();
     cellV = HFNoseDetId(id).cellV();
   } else if (id.det() == DetId::Forward) {
     cellU = HGCalDetId(id).cell();
   } else {
     edm::LogError("getCell::InvalidSiliconDetid")
         << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal silicon!";
   }
   return std::pair<int, int>(cellU, cellV);
 }

◆ getCellType()

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

Definition at line 439 of file RecHitTools.cc.

References cms::cuda::assert(), DetId::Forward, geom_, getLayerWithOffset(), getSiThickIndex(), CaloGeometry::getSubdetectorGeometry(), HFNose, isScintillator(), lastLayerEE(), Calorimetry_cff::thickness, and funct::true.

Referenced by TICLDumper::analyze().

                                                   {
   auto layer_number = getLayerWithOffset(id);
   auto thickness = getSiThickIndex(id);
   auto geomNose =
       static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HFNose));
   auto isNose = geomNose ? true : false;
   auto isEELayer = (layer_number <= lastLayerEE(isNose));
   auto isScint = isScintillator(id);
   int layerType = -1;
 
   if (isScint) {
     layerType = CE_H_SCINT;
   }
   if (isEELayer) {
     if (thickness == 0) {
       layerType = CE_E_120;
     } else if (thickness == 1) {
       layerType = CE_E_200;
     } else if (thickness == 2) {
       layerType = CE_E_300;
     }
   } else {
     if (thickness == 0) {
       layerType = CE_H_120;
     } else if (thickness == 1) {
       layerType = CE_H_200;
     } else if (thickness == 2) {
       layerType = CE_H_300;
     }
   }
   assert(layerType != -1);
   return layerType;
 }

◆ getEta() [1/2]

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

Definition at line 486 of file RecHitTools.cc.

References PV3DBase< T, PVType, FrameType >::eta(), and position.

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

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

◆ getEta() [2/2]

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

Definition at line 491 of file RecHitTools.cc.

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

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

◆ getGeometry()

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

inline

Definition at line 75 of file RecHitTools.h.

References geom_.

75 { return geom_; };

hgcal::RecHitTools::geom_

const CaloGeometry * geom_

Definition: RecHitTools.h:90

◆ getGeometryType()

int hgcal::RecHitTools::getGeometryType ( ) const

inline

Definition at line 86 of file RecHitTools.h.

References geometryType_.

Referenced by RealisticSimClusterMapper::buildClusters().

86 { return geometryType_; }

hgcal::RecHitTools::geometryType_

int geometryType_

Definition: RecHitTools.h:93

◆ getLayer() [1/3]

unsigned int RecHitTools::getLayer	(	DetId::Detector	type,
		bool	nose = `false`
	)		const

Definition at line 318 of file RecHitTools.cc.

References DetId::Forward, ForwardEmpty, geom_, CaloGeometry::getSubdetectorGeometry(), HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, and hgcalTBTopologyTester_cfi::layers.

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

                                                                             {
   int layer;
   switch (type) {
     case (DetId::HGCalEE): {
       auto geomEE =
           static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(type, ForwardSubdetector::ForwardEmpty));
       layer = (geomEE->topology().dddConstants()).layers(true);
       break;
     }
     case (DetId::HGCalHSi): {
       auto geomFH =
           static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(type, ForwardSubdetector::ForwardEmpty));
       layer = (geomFH->topology().dddConstants()).layers(true);
       break;
     }
     case (DetId::HGCalHSc): {
       auto geomBH =
           static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(type, ForwardSubdetector::ForwardEmpty));
       layer = (geomBH->topology().dddConstants()).layers(true);
       break;
     }
     case (DetId::Forward): {
       if (nose) {
         auto geomHFN = static_cast<const HGCalGeometry*>(
             geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HFNose));
         layer = (geomHFN->topology().dddConstants()).layers(true);
       } else {
         auto geomEE = static_cast<const HGCalGeometry*>(
             geom_->getSubdetectorGeometry(DetId::HGCalEE, ForwardSubdetector::ForwardEmpty));
         layer = (geomEE->topology().dddConstants()).layers(true);
         auto geomFH = static_cast<const HGCalGeometry*>(
             geom_->getSubdetectorGeometry(DetId::HGCalHSi, ForwardSubdetector::ForwardEmpty));
         layer += (geomFH->topology().dddConstants()).layers(true);
       }
       break;
     }
     default:
       layer = 0;
   }
   return (unsigned int)(layer);
 }

◆ getLayer() [2/3]

unsigned int RecHitTools::getLayer ( ForwardSubdetector type ) const

Definition at line 268 of file RecHitTools.cc.

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

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

◆ getLayer() [3/3]

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

Definition at line 360 of file RecHitTools.cc.

References HcalDetId::depth(), DetId::Forward, DetId::Hcal, HcalEndcap, HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, HGCalDetId::layer(), HFNoseDetId::layer(), HGCScintillatorDetId::layer(), HGCSiliconDetId::layer(), and SiStripPI::max.

                                                         {
   unsigned int layer = std::numeric_limits<unsigned int>::max();
   if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
     layer = HGCSiliconDetId(id).layer();
   } else if (id.det() == DetId::HGCalHSc) {
     layer = HGCScintillatorDetId(id).layer();
   } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
     layer = HFNoseDetId(id).layer();
   } else if (id.det() == DetId::Forward) {
     layer = HGCalDetId(id).layer();
   } else if (id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
     layer = HcalDetId(id).depth();
   }
   return layer;
 }

◆ getLayerWithOffset()

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

Definition at line 376 of file RecHitTools.cc.

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

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), HGCalShowerSeparation::analyze(), HGCalHitCalibration::analyze(), TICLDumper::analyze(), RealisticSimClusterMapper::buildClusters(), TICLTrackstersEdgesValidation::dqmAnalyze(), TrackstersMergeProducerV3::energyRegressionAndID(), TrackstersMergeProducer::energyRegressionAndID(), HGCalHitCalibration::fillWithRecHits(), ticl::ClusterFilterByAlgoAndSizeAndLayerRange::filter(), getCellType(), hgcal::ClusterTools::getLayer(), HGCalRecHitSimpleAlgo::makeRecHit(), HGCal3DClustering::organizeByLayer(), TICLLayerTileProducer::produce(), 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::HGCalHSi || id.det() == DetId::HGCalHSc) {
     // DetId::HGCalHSc hits include the offset w.r.t. EE already
     layer += fhOffset_;
   } else if (id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
     layer += bhOffset_;
   }
   // no need to add offset for HFnose
   return layer;
 }

◆ getPhi() [1/2]

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

Definition at line 497 of file RecHitTools.cc.

References position.

Referenced by HGCalIsoCalculator::setRecHits().

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

◆ getPhi() [2/2]

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

Definition at line 502 of file RecHitTools.cc.

References getPosition(), and position.

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

◆ getPosition()

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

Definition at line 140 of file RecHitTools.cc.

References relativeConstraints::geom, HGCalGeometry::getPosition(), getSubdetectorGeometry(), DetId::Hcal, and position.

Referenced by HGCalShowerSeparation::analyze(), RealisticSimClusterMapper::buildClusters(), HGCalLayerClusterProducer::calculatePosition(), getEta(), getPhi(), getPt(), hgcal::ClusterTools::getWidths(), HGCalIsoCalculator::setRecHits(), and hgcal::EGammaPCAHelper::storeRecHits().

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

◆ getPositionLayer()

GlobalPoint RecHitTools::getPositionLayer	(	int	layer,
		bool	nose = `false`
	)		const

Definition at line 152 of file RecHitTools.cc.

References funct::abs(), fhOffset_, DetId::Forward, geom_, geometryType_, CaloGeometry::getSubdetectorGeometry(), HFNose, and HGCalDDDConstants::waferZ().

Referenced by PositionAtECalEntranceComputer::beginEvent(), ticl::LinkingAlgoByDirectionGeometric::buildLayers(), TICLDumper::buildLayers(), hgcal::EGammaPCAHelper::findZFirstLayer(), ticl::LinkingAlgoByDirectionGeometric::linkTracksters(), TrackstersMergeProducerV3::produce(), SimTrackstersProducer::produce(), and TrackstersMergeProducer::produce().

                                                                     {
   unsigned int lay = std::abs(layer);
   double z(0);
   if (nose) {
     auto geomNose =
         static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HFNose));
     if (geomNose) {
       const HGCalDDDConstants* ddd = &(geomNose->topology().dddConstants());
       if (ddd)
         z = (layer > 0) ? ddd->waferZ(lay, true) : -ddd->waferZ(lay, true);
     }
   } else {
     const HGCalDDDConstants* ddd = get_ddd(geom_, geometryType_, fhOffset_, lay);
     if (geometryType_ == 1) {
       if (lay > fhOffset_)
         lay -= fhOffset_;
     }
     z = (layer > 0) ? ddd->waferZ(lay, true) : -ddd->waferZ(lay, true);
   }
   return GlobalPoint(0, 0, z);
 }

◆ getPt() [1/2]

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

Definition at line 508 of file RecHitTools.cc.

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

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

◆ getPt() [2/2]

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

Definition at line 514 of file RecHitTools.cc.

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

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

◆ getRadiusToSide()

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

Definition at line 246 of file RecHitTools.cc.

References HGCalDDDConstants::cellSizeHex(), TauDecayModes::dec, DetId::Forward, relativeConstraints::geom, getSubdetectorGeometry(), HFNose, DetId::HGCalEE, DetId::HGCalHSi, SiStripPI::max, HFNoseDetId::type(), HGCSiliconDetId::type(), and HGCalDetId::waferType().

                                                              {
   auto geom = getSubdetectorGeometry(id);
   std::float_t size(std::numeric_limits<std::float_t>::max());
   if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
     const HGCSiliconDetId hid(id);
     auto ddd = get_ddd(geom, hid);
     size = ddd->cellSizeHex(hid.type());
   } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
     const HFNoseDetId hid(id);
     auto ddd = get_ddd(geom, hid);
     size = ddd->cellSizeHex(hid.type());
   } else if (id.det() == DetId::Forward) {
     const HGCalDetId hid(id);
     auto ddd = get_ddd(geom, hid);
     size = ddd->cellSizeHex(hid.waferType());
   } else {
     edm::LogError("getRadiusToSide::InvalidSiliconDetid")
         << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal silicon!";
   }
   return size;
 }

◆ getScintDEtaDPhi()

std::pair< float, float > RecHitTools::getScintDEtaDPhi ( const DetId & id ) const

Definition at line 236 of file RecHitTools.cc.

References TauDecayModes::dec, CaloSubdetectorGeometry::getGeometry(), getSubdetectorGeometry(), isScintillator(), and LogDebug.

                                                                          {
   if (!isScintillator(id)) {
     LogDebug("getScintDEtaDPhi::InvalidScintDetid")
         << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal scintillator!";
     return {0.f, 0.f};
   }
   auto cellGeom = getSubdetectorGeometry(id)->getGeometry(id);
   return {cellGeom->etaSpan(), cellGeom->phiSpan()};
 }

◆ getScintMaxIphi()

int hgcal::RecHitTools::getScintMaxIphi ( ) const

inline

Definition at line 85 of file RecHitTools.h.

References bhMaxIphi_.

Referenced by HGCalClusteringAlgoBase::getEventSetup().

85 { return bhMaxIphi_; }

hgcal::RecHitTools::bhMaxIphi_

int bhMaxIphi_

Definition: RecHitTools.h:94

◆ getSiThickIndex()

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

Definition at line 216 of file RecHitTools.cc.

References cms::cuda::assert(), DetId::Forward, getSiThickness(), HFNose, DetId::HGCalEE, DetId::HGCalHSi, Calorimetry_cff::thickness, HFNoseDetId::type(), and HGCSiliconDetId::type().

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), HGCalLayerClusterProducer::calculatePosition(), getCellType(), and hgcal::EGammaPCAHelper::storeRecHits().

                                                       {
   int thickIndex = -1;
   if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
     thickIndex = HGCSiliconDetId(id).type();
   } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
     thickIndex = HFNoseDetId(id).type();
   } else if (id.det() == DetId::Forward) {
     float thickness = getSiThickness(id);
     if (thickness > 99. && thickness < 121.)
       thickIndex = 0;
     else if (thickness > 199. && thickness < 201.)
       thickIndex = 1;
     else if (thickness > 299. && thickness < 301.)
       thickIndex = 2;
     else
       assert(thickIndex > 0 && "ERROR - silicon thickness has a nonsensical value");
   }
   return thickIndex;
 }

◆ getSiThickness()

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

Definition at line 190 of file RecHitTools.cc.

References HGCalDDDConstants::cellThickness(), TauDecayModes::dec, DetId::Forward, relativeConstraints::geom, getSubdetectorGeometry(), HFNose, DetId::HGCalEE, DetId::HGCalHSi, HGCalDetId::layer(), HFNoseDetId::layer(), HGCSiliconDetId::layer(), LogDebug, HGCalDetId::wafer(), HFNoseDetId::waferU(), HGCSiliconDetId::waferU(), HFNoseDetId::waferV(), and HGCSiliconDetId::waferV().

Referenced by HGCalHitCalibration::analyze(), HGCalHitCalibration::fillWithRecHits(), and getSiThickIndex().

                                                             {
   auto geom = getSubdetectorGeometry(id);
   std::float_t thick(0.37);
   if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
     const HGCSiliconDetId hid(id);
     auto ddd = get_ddd(geom, hid);
     thick = ddd->cellThickness(hid.layer(), hid.waferU(), hid.waferV());
   } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
     const HFNoseDetId hid(id);
     auto ddd = get_ddd(geom, hid);
     thick = ddd->cellThickness(hid.layer(), hid.waferU(), hid.waferV());
   } else if (id.det() == DetId::Forward) {
     const HGCalDetId hid(id);
     auto ddd = get_ddd(geom, hid);
     thick = ddd->cellThickness(hid.layer(), hid.wafer(), 0);
   } else {
     LogDebug("getSiThickness::InvalidSiliconDetid")
         << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal silicon!";
     // It does not make any sense to return 0.37 as thickness for a detector
     // that is not Silicon(does it?). Mark it as 0. to be able to distinguish
     // the case.
     thick = 0.f;
   }
   return thick;
 }

◆ getSubdetectorGeometry()

const CaloSubdetectorGeometry * RecHitTools::getSubdetectorGeometry ( const DetId & id ) const

Definition at line 130 of file RecHitTools.cc.

References ForwardEmpty, relativeConstraints::geom, geom_, CaloGeometry::getSubdetectorGeometry(), DetId::HGCalEE, DetId::HGCalHSc, and DetId::HGCalHSi.

Referenced by getPosition(), getRadiusToSide(), getScintDEtaDPhi(), getSiThickness(), isHalfCell(), and maskCell().

                                                                                         {
   DetId::Detector det = id.det();
   int subdet = (det == DetId::HGCalEE || det == DetId::HGCalHSi || det == DetId::HGCalHSc)
                    ? ForwardSubdetector::ForwardEmpty
                    : id.subdetId();
   auto geom = geom_->getSubdetectorGeometry(det, subdet);
   check_geom(geom);
   return geom;
 }

◆ getWafer()

std::pair< int, int > RecHitTools::getWafer ( const DetId & id ) const

Definition at line 390 of file RecHitTools.cc.

References TauDecayModes::dec, DetId::Forward, HFNose, DetId::HGCalEE, DetId::HGCalHSi, SiStripPI::max, HGCalDetId::wafer(), HGCalWaferIndex::waferU(), HFNoseDetId::waferU(), HGCSiliconDetId::waferU(), HGCalWaferIndex::waferV(), HFNoseDetId::waferV(), and HGCSiliconDetId::waferV().

                                                              {
   int waferU = std::numeric_limits<int>::max();
   int waferV = 0;
   if ((id.det() == DetId::HGCalEE) || (id.det() == DetId::HGCalHSi)) {
     waferU = HGCSiliconDetId(id).waferU();
     waferV = HGCSiliconDetId(id).waferV();
   } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
     waferU = HFNoseDetId(id).waferU();
     waferV = HFNoseDetId(id).waferV();
   } else if (id.det() == DetId::Forward) {
     waferU = HGCalDetId(id).wafer();
   } else {
     edm::LogError("getWafer::InvalidSiliconDetid")
         << "det id: " << std::hex << id.rawId() << std::dec << ":" << id.det() << " is not HGCal silicon!";
   }
   return std::pair<int, int>(waferU, waferV);
 }

◆ isHalfCell()

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

Definition at line 426 of file RecHitTools.cc.

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

                                                   {
   bool ishalf = false;
   if (id.det() == DetId::Forward) {
     HGCalDetId hid(id);
     auto geom = getSubdetectorGeometry(hid);
     auto ddd = get_ddd(geom, hid);
     const int waferType = ddd->waferTypeT(hid.waferType());
     return ddd->isHalfCell(waferType, hid.cell());
   }
   //new geometry is always false
   return ishalf;
 }

◆ isOnlySilicon()

bool RecHitTools::isOnlySilicon ( const unsigned int layer ) const

Definition at line 480 of file RecHitTools.cc.

References bhLastLayer_, and bhOffset_.

                                                               {
   // HFnose TODO
   bool isonlysilicon = (layer % bhLastLayer_) < bhOffset_;
   return isonlysilicon;
 }

◆ isScintillator()

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

Definition at line 478 of file RecHitTools.cc.

References DetId::HGCalHSc.

Referenced by getCellType(), and getScintDEtaDPhi().

478 { return id.det() == DetId::HGCalHSc; }

DetId::HGCalHSc

Definition: DetId.h:34

◆ isSilicon()

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

Definition at line 473 of file RecHitTools.cc.

References DetId::Forward, HFNose, DetId::HGCalEE, and DetId::HGCalHSi.

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), ticl::ClusterFilterByAlgoAndSize::filter(), and ticl::ClusterFilterByAlgoAndSizeAndLayerRange::filter().

                                                  {
   return (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi ||
           (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)));
 }

◆ lastLayer()

unsigned int hgcal::RecHitTools::lastLayer ( bool nose = false ) const

inline

Definition at line 80 of file RecHitTools.h.

References bhLastLayer_, and noseLastLayer_.

Referenced by HGCalClusteringAlgoBase::getEventSetup(), TICLLayerTileProducer::produce(), and setGeometry().

80 { return (nose ? noseLastLayer_ : bhLastLayer_); }

hgcal::RecHitTools::bhLastLayer_

unsigned int bhLastLayer_

Definition: RecHitTools.h:91

hgcal::RecHitTools::noseLastLayer_

unsigned int noseLastLayer_

Definition: RecHitTools.h:91

◆ lastLayerBH()

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

inline

Definition at line 79 of file RecHitTools.h.

References bhLastLayer_.

Referenced by HGCal3DClustering::getEventSetup(), and hgcal::EGammaPCAHelper::setRecHitTools().

79 { return bhLastLayer_; }

hgcal::RecHitTools::bhLastLayer_

unsigned int bhLastLayer_

Definition: RecHitTools.h:91

◆ lastLayerEE()

unsigned int hgcal::RecHitTools::lastLayerEE ( bool nose = false ) const

inline

Definition at line 76 of file RecHitTools.h.

References fhOffset_, and HFNoseDetId::HFNoseLayerEEmax.

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), RealisticSimClusterMapper::buildClusters(), ticl::LinkingAlgoByDirectionGeometric::buildLayers(), TICLDumper::buildLayers(), hgcal::EGammaPCAHelper::computePCA(), getCellType(), HGCalClusteringAlgoBase::getEventSetup(), hgcal::ClusterTools::getWidths(), ticl::LinkingAlgoByDirectionGeometric::linkTracksters(), HGCal3DClustering::makeClusters(), HGCalDepthPreClusterer::makePreClusters(), TrackstersMergeProducerV3::produce(), SimTrackstersProducer::produce(), and TrackstersMergeProducer::produce().

76 { return (nose ? HFNoseDetId::HFNoseLayerEEmax : fhOffset_); }

hgcal::RecHitTools::fhOffset_

unsigned int fhOffset_

Definition: RecHitTools.h:91

HFNoseDetId::HFNoseLayerEEmax

static const int HFNoseLayerEEmax

Definition: HFNoseDetId.h:29

◆ lastLayerFH()

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

inline

Definition at line 77 of file RecHitTools.h.

References fhLastLayer_.

Referenced by RealisticSimClusterMapper::buildClusters(), and HGCalClusteringAlgoBase::getEventSetup().

77 { return fhLastLayer_; }

hgcal::RecHitTools::fhLastLayer_

unsigned int fhLastLayer_

Definition: RecHitTools.h:91

◆ maskCell()

bool RecHitTools::maskCell	(	const DetId &	id,
		int	corners = `3`
	)		const

Definition at line 533 of file RecHitTools.cc.

References HGCalTopology::dddConstants(), getSubdetectorGeometry(), DetId::Hcal, HGCalDDDConstants::maskCell(), and HGCalGeometry::topology().

                                                              {
   if (id.det() == DetId::Hcal) {
     return false;
   } else {
     auto hg = static_cast<const HGCalGeometry*>(getSubdetectorGeometry(id));
     return hg->topology().dddConstants().maskCell(id, corners);
   }
 }

◆ maxNumberOfWafersPerLayer()

unsigned int hgcal::RecHitTools::maxNumberOfWafersPerLayer ( bool nose = false ) const

inline

Definition at line 82 of file RecHitTools.h.

References maxNumberOfWafersNose_, and maxNumberOfWafersPerLayer_.

                                                                     {
       return (nose ? maxNumberOfWafersNose_ : maxNumberOfWafersPerLayer_);
     }

◆ setGeometry()

void RecHitTools::setGeometry ( CaloGeometry const & geom )

Definition at line 79 of file RecHitTools.cc.

References bhFirstLayer_, bhLastLayer_, bhMaxIphi_, bhOffset_, eeOffset_, fhLastLayer_, fhOffset_, DetId::Forward, ForwardEmpty, relativeConstraints::geom, geom_, geometryType_, CaloGeometry::getSubdetectorGeometry(), DetId::Hcal, HcalEndcap, HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, HGCEE, HGCHEF, lastLayer(), hgcalTBTopologyTester_cfi::layers, SiStripPI::max, maxNumberOfWafersNose_, maxNumberOfWafersPerLayer_, and noseLastLayer_.

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), HeterogeneousHGCalRecHitsValidator::analyze(), HGCalShowerSeparation::analyze(), HGCalHitCalibration::analyze(), PositionAtECalEntranceComputer::beginEvent(), TICLLayerTileProducer::beginRun(), FilteredLayerClustersProducer::beginRun(), TrackstersMergeProducer::beginRun(), TICLDumper::beginRun(), TICLTrackstersEdgesValidation::dqmBeginRun(), HGCalEgammaIDHelper::eventInit(), hgcal::ClusterTools::getEventSetup(), HGCal3DClustering::getEventSetup(), HGCalDepthPreClusterer::getEventSetup(), PFHGCalRecHitCreator< DET, Layer, det, subdet >::importRecHits(), TrackstersMergeProducerV3::produce(), SimTrackstersProducer::produce(), HGCalLayerClusterProducer::produce(), HGCalRecHitAbsAlgo::set(), and RealisticSimClusterMapper::update().

                                                       {
   geom_ = &geom;
   unsigned int wmaxEE(0), wmaxFH(0);
   auto geomEE = static_cast<const HGCalGeometry*>(
       geom_->getSubdetectorGeometry(DetId::HGCalEE, ForwardSubdetector::ForwardEmpty));
   //check if it's the new geometry
   if (geomEE) {
     geometryType_ = 1;
     eeOffset_ = (geomEE->topology().dddConstants()).getLayerOffset();
     wmaxEE = (geomEE->topology().dddConstants()).waferCount(0);
     auto geomFH = static_cast<const HGCalGeometry*>(
         geom_->getSubdetectorGeometry(DetId::HGCalHSi, ForwardSubdetector::ForwardEmpty));
     fhOffset_ = (geomFH->topology().dddConstants()).getLayerOffset();
     wmaxFH = (geomFH->topology().dddConstants()).waferCount(0);
     fhLastLayer_ = fhOffset_ + (geomFH->topology().dddConstants()).lastLayer(true);
     auto geomBH = static_cast<const HGCalGeometry*>(
         geom_->getSubdetectorGeometry(DetId::HGCalHSc, ForwardSubdetector::ForwardEmpty));
     bhOffset_ = (geomBH->topology().dddConstants()).getLayerOffset();
     bhFirstLayer_ = bhOffset_ + (geomBH->topology().dddConstants()).firstLayer();
     bhLastLayer_ = bhOffset_ + (geomBH->topology().dddConstants()).lastLayer(true);
     bhMaxIphi_ = geomBH->topology().dddConstants().maxCells(true);
   } else {
     geometryType_ = 0;
     geomEE =
         static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HGCEE));
     eeOffset_ = (geomEE->topology().dddConstants()).getLayerOffset();
     wmaxEE = 1 + (geomEE->topology().dddConstants()).waferMax();
     auto geomFH =
         static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HGCHEF));
     fhOffset_ = (geomFH->topology().dddConstants()).getLayerOffset();
     fhLastLayer_ = fhOffset_ + (geomFH->topology().dddConstants()).layers(true);
     bhOffset_ = fhLastLayer_;
     bhFirstLayer_ = bhOffset_ + 1;
     wmaxFH = 1 + (geomFH->topology().dddConstants()).waferMax();
     auto geomBH =
         static_cast<const HcalGeometry*>(geom_->getSubdetectorGeometry(DetId::Hcal, HcalSubdetector::HcalEndcap));
     bhLastLayer_ = bhOffset_ + (geomBH->topology().dddConstants())->getMaxDepth(1);
   }
   maxNumberOfWafersPerLayer_ = std::max(wmaxEE, wmaxFH);
   // For nose geometry
   auto geomNose =
       static_cast<const HGCalGeometry*>(geom_->getSubdetectorGeometry(DetId::Forward, ForwardSubdetector::HFNose));
   if (geomNose) {
     maxNumberOfWafersNose_ = (geomNose->topology().dddConstants()).waferCount(0);
     noseLastLayer_ = (geomNose->topology().dddConstants()).layers(true);
   } else {
     maxNumberOfWafersNose_ = 0;
     noseLastLayer_ = 0;
   }
 }

◆ zside()

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

Definition at line 174 of file RecHitTools.cc.

References DetId::Forward, DetId::Hcal, HcalEndcap, HFNose, DetId::HGCalEE, DetId::HGCalHSc, DetId::HGCalHSi, HGCalDetId::zside(), HFNoseDetId::zside(), HGCScintillatorDetId::zside(), HGCSiliconDetId::zside(), and HcalDetId::zside().

Referenced by FWRecoGeometryESProducer::addCaloGeometry(), and TICLLayerTileProducer::produce().

                                             {
   int zside = 0;
   if (id.det() == DetId::HGCalEE || id.det() == DetId::HGCalHSi) {
     zside = HGCSiliconDetId(id).zside();
   } else if (id.det() == DetId::HGCalHSc) {
     zside = HGCScintillatorDetId(id).zside();
   } else if (id.det() == DetId::Forward && id.subdetId() == static_cast<int>(HFNose)) {
     zside = HFNoseDetId(id).zside();
   } else if (id.det() == DetId::Forward) {
     zside = HGCalDetId(id).zside();
   } else if (id.det() == DetId::Hcal && id.subdetId() == HcalEndcap) {
     zside = HcalDetId(id).zside();
   }
   return zside;
 }

Member Data Documentation

◆ bhFirstLayer_

unsigned int hgcal::RecHitTools::bhFirstLayer_

private

Definition at line 91 of file RecHitTools.h.

Referenced by firstAndLastLayer(), firstLayerBH(), and setGeometry().

◆ bhLastLayer_

unsigned int hgcal::RecHitTools::bhLastLayer_

private

Definition at line 91 of file RecHitTools.h.

Referenced by firstAndLastLayer(), isOnlySilicon(), lastLayer(), lastLayerBH(), and setGeometry().

◆ bhMaxIphi_

int hgcal::RecHitTools::bhMaxIphi_

private

Definition at line 94 of file RecHitTools.h.

Referenced by getScintMaxIphi(), and setGeometry().

◆ bhOffset_

unsigned int hgcal::RecHitTools::bhOffset_

private

Definition at line 91 of file RecHitTools.h.

Referenced by getLayerWithOffset(), isOnlySilicon(), and setGeometry().

◆ eeOffset_

unsigned int hgcal::RecHitTools::eeOffset_

private

Definition at line 91 of file RecHitTools.h.

Referenced by firstAndLastLayer(), and setGeometry().

◆ fhLastLayer_

unsigned int hgcal::RecHitTools::fhLastLayer_

private

Definition at line 91 of file RecHitTools.h.

Referenced by firstAndLastLayer(), lastLayerFH(), and setGeometry().

◆ fhOffset_

unsigned int hgcal::RecHitTools::fhOffset_

private

Definition at line 91 of file RecHitTools.h.

Referenced by firstAndLastLayer(), getLayerWithOffset(), getPositionLayer(), lastLayerEE(), and setGeometry().

◆ geom_

const CaloGeometry* hgcal::RecHitTools::geom_

private

Definition at line 90 of file RecHitTools.h.

Referenced by getCellType(), getGeometry(), getLayer(), getPositionLayer(), getSubdetectorGeometry(), and setGeometry().

◆ geometryType_

int hgcal::RecHitTools::geometryType_

private

Definition at line 93 of file RecHitTools.h.

Referenced by getGeometryType(), getPositionLayer(), and setGeometry().

◆ maxNumberOfWafersNose_

unsigned int hgcal::RecHitTools::maxNumberOfWafersNose_

private

Definition at line 92 of file RecHitTools.h.

Referenced by maxNumberOfWafersPerLayer(), and setGeometry().

◆ maxNumberOfWafersPerLayer_

unsigned int hgcal::RecHitTools::maxNumberOfWafersPerLayer_

private

Definition at line 92 of file RecHitTools.h.

Referenced by maxNumberOfWafersPerLayer(), and setGeometry().

◆ noseLastLayer_

unsigned int hgcal::RecHitTools::noseLastLayer_

private

Definition at line 91 of file RecHitTools.h.

Referenced by firstAndLastLayer(), lastLayer(), and setGeometry().

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ RecHitTools()

◆ ~RecHitTools()

Member Function Documentation

◆ firstAndLastLayer()

◆ firstLayerBH()

◆ getCell()

◆ getCellType()

◆ getEta() [1/2]

◆ getEta() [2/2]

◆ getGeometry()

◆ getGeometryType()

◆ getLayer() [1/3]

◆ getLayer() [2/3]

◆ getLayer() [3/3]

◆ getLayerWithOffset()

◆ getPhi() [1/2]

◆ getPhi() [2/2]

◆ getPosition()

◆ getPositionLayer()

◆ getPt() [1/2]

◆ getPt() [2/2]

◆ getRadiusToSide()

◆ getScintDEtaDPhi()

◆ getScintMaxIphi()

◆ getSiThickIndex()

◆ getSiThickness()

◆ getSubdetectorGeometry()

◆ getWafer()

◆ isHalfCell()

◆ isOnlySilicon()

◆ isScintillator()

◆ isSilicon()

◆ lastLayer()

◆ lastLayerBH()

◆ lastLayerEE()

◆ lastLayerFH()

◆ maskCell()

◆ maxNumberOfWafersPerLayer()

◆ setGeometry()

◆ zside()

Member Data Documentation

◆ bhFirstLayer_

◆ bhLastLayer_

◆ bhMaxIphi_

◆ bhOffset_

◆ eeOffset_

◆ fhLastLayer_

◆ fhOffset_

◆ geom_

◆ geometryType_

◆ maxNumberOfWafersNose_

◆ maxNumberOfWafersPerLayer_

◆ noseLastLayer_