#include <HcalHitRelabeller.h>

Public Member Functions
	HcalHitRelabeller (bool nd=false)

void	process (std::vector< PCaloHit > &hcalHits)

DetId	relabel (const uint32_t testId) const

void	setGeometry (const HcalDDDRecConstants *&)

Static Public Member Functions
static DetId	relabel (const uint32_t testId, const HcalDDDRecConstants *theRecNumber)

Private Member Functions
double	energyWt (const uint32_t testId) const

Private Attributes
bool	neutralDensity_

const HcalDDDRecConstants *	theRecNumber

Detailed Description

Definition at line 10 of file HcalHitRelabeller.h.

Constructor & Destructor Documentation

◆ HcalHitRelabeller()

HcalHitRelabeller::HcalHitRelabeller ( bool nd = false )

Definition at line 9 of file HcalHitRelabeller.cc.

References neutralDensity_.

                                             : theRecNumber(nullptr), neutralDensity_(nd) {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "HcalHitRelabeller initialized with"
                               << " neutralDensity " << neutralDensity_;
 #endif
 }

Member Function Documentation

◆ energyWt()

double HcalHitRelabeller::energyWt ( const uint32_t testId ) const

private

Definition at line 88 of file HcalHitRelabeller.cc.

References LEDCalibrationChannels::depth, PVValHelper::eta, HcalDDDRecConstants::getLayer0Wt(), phase1PixelTopology::layer, phi, theRecNumber, HcalTestNumbering::unpackHcalIndex(), z, and ecaldqm::zside().

Referenced by process().

                                                               {
   int det, z, depth, eta, phi, layer;
   HcalTestNumbering::unpackHcalIndex(testId, det, z, depth, eta, phi, layer);
   int zside = (z == 0) ? (-1) : (1);
   double wt = (((det == 1) || (det == 2)) && (depth == 1)) ? theRecNumber->getLayer0Wt(det, phi, zside) : 1.0;
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "EnergyWT::det: " << det << " z: " << z << ":" << zside << " depth: " << depth
                               << " ieta: " << eta << " iphi: " << phi << " layer: " << layer << " wt " << wt;
 #endif
   return wt;
 }

◆ process()

void HcalHitRelabeller::process ( std::vector< PCaloHit > & hcalHits )

Definition at line 16 of file HcalHitRelabeller.cc.

References TauDecayModes::dec, HCALHighEnergyHPDFilter_cfi::energy, energyWt(), cuy::ii, neutralDensity_, DetId::rawId(), relabel(), and theRecNumber.

                                                              {
   if (theRecNumber) {
 #ifdef EDM_ML_DEBUG
     int ii(0);
 #endif
     for (auto& hcalHit : hcalHits) {
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalSim") << "Hit[" << ii << "] " << std::hex << hcalHit.id() << std::dec << " Neutral density "
                                   << neutralDensity_;
 #endif
       double energy = (hcalHit.energy());
       if (neutralDensity_) {
         energy *= (energyWt(hcalHit.id()));
         hcalHit.setEnergy(energy);
       }
       DetId newid = relabel(hcalHit.id());
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalSim") << "Hit " << ii << " out of " << hcalHits.size() << " " << std::hex << newid.rawId()
                                   << std::dec << " E " << energy;
 #endif
       hcalHit.setID(newid.rawId());
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalSim") << "Modified Hit " << HcalDetId(hcalHit.id());
       ++ii;
 #endif
     }
   } else {
     edm::LogWarning("HcalSim") << "HcalHitRelabeller: no valid HcalDDDRecConstants";
   }
 }

◆ relabel() [1/2]

DetId HcalHitRelabeller::relabel ( const uint32_t testId ) const

Definition at line 49 of file HcalHitRelabeller.cc.

References theRecNumber.

Referenced by SimHitsValidationHcal::analyze(), HcalTestSimHitID::analyze(), HcalSimHitAnalysis::analyze(), HcalSimHitsValidation::analyze(), HcalTestNumberingTester::analyze(), HcalGeomCheck::analyze(), HcalRecHitsValidation::analyze(), HcalSimHitStudy::analyzeHits(), HcalSimHitCheck::analyzeHits(), HcalTestNS::compare(), GlobalHitsAnalyzer::fillHCal(), CaloParticleDebugger::fillSimHits(), CaloTruthAccumulator::fillSimHits(), process(), and HcalDigisValidation::reco().

                                                             {
   return HcalHitRelabeller::relabel(testId, theRecNumber);
 }

◆ relabel() [2/2]

DetId HcalHitRelabeller::relabel	(	const uint32_t	testId,
		const HcalDDDRecConstants *	theRecNumber
	)

static

Definition at line 53 of file HcalHitRelabeller.cc.

References TauDecayModes::dec, LEDCalibrationChannels::depth, PVValHelper::eta, HcalDDDRecConstants::getHCID(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, phase1PixelTopology::layer, phi, DetId::rawId(), Validation_hcalonly_cfi::sign, theRecNumber, HcalTestNumbering::unpackHcalIndex(), and z.

                                                                                                {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "Enter HcalHitRelabeller::relabel";
 #endif
   HcalDetId hid;
   int det, z, depth, eta, phi, layer, sign;
   HcalTestNumbering::unpackHcalIndex(testId, det, z, depth, eta, phi, layer);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "det: " << det << " "
                               << "z: " << z << " "
                               << "depth: " << depth << " "
                               << "ieta: " << eta << " "
                               << "iphi: " << phi << " "
                               << "layer: " << layer;
 #endif
   sign = (z == 0) ? (-1) : (1);
   HcalDDDRecConstants::HcalID id = theRecNumber->getHCID(det, sign * eta, phi, layer, depth);
 
   if (id.subdet == int(HcalBarrel)) {
     hid = HcalDetId(HcalBarrel, sign * id.eta, id.phi, id.depth);
   } else if (id.subdet == int(HcalEndcap)) {
     hid = HcalDetId(HcalEndcap, sign * id.eta, id.phi, id.depth);
   } else if (id.subdet == int(HcalOuter)) {
     hid = HcalDetId(HcalOuter, sign * id.eta, id.phi, id.depth);
   } else if (id.subdet == int(HcalForward)) {
     hid = HcalDetId(HcalForward, sign * id.eta, id.phi, id.depth);
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << " new HcalDetId -> hex.RawID = " << std::hex << hid.rawId() << std::dec
                               << " det, z, depth, eta, phi = " << det << " " << z << " " << id.depth << " " << id.eta
                               << " " << id.phi << " ---> " << hid;
 #endif
   return hid;
 }