#include <HcalHitRelabeller.h>

Public Member Functions
double	energyWt (const uint32_t testId) const

	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 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

Definition at line 95 of file HcalHitRelabeller.cc.

References hcalRecHitTable_cff::depth, PVValHelper::eta, HcalDDDRecConstants::getLayer0Wt(), nano_mu_digi_cff::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) && (layer <= 1)) || ((det == 2) && (layer <= 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, hcalRecHitTable_cff::depth, HcalDetId::depth(), hcalRecHitTable_cff::energy, energyWt(), PVValHelper::eta, HcalDetId::ietaAbs(), cuy::ii, HcalDetId::iphi(), nano_mu_digi_cff::layer, neutralDensity_, HcalTestNumbering::packHcalIndex(), phi, DetId::rawId(), relabel(), HcalDetId::subdet(), theRecNumber, HcalTestNumbering::unpackHcalIndex(), z, and HcalDetId::zside().

                                                              {
   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
       uint32_t hid;
       int det, z, depth, eta, phi, layer;
       HcalTestNumbering::unpackHcalIndex(hcalHit.id(), det, z, depth, eta, phi, layer);
       if ((det == 2) && (layer == 2) && (eta == 18))
         depth = 2;
       hid = HcalTestNumbering::packHcalIndex(det, z, depth, eta, phi, layer);
       double wt = (neutralDensity_) ? (energyWt(hid)) : 1.0;
       double energy = wt * (hcalHit.energy());
       hcalHit.setEnergy(energy);
       uint32_t newid = relabel(hid).rawId();
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalSim") << "Hit " << ii << " out of " << hcalHits.size() << " " << std::hex << newid
                                   << std::dec << " E " << energy << " wt " << wt;
 #endif
       hcalHit.setID(newid);
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HcalSim") << "Final Setting::subdet: " << HcalDetId(newid).subdet()
                                   << " z: " << HcalDetId(newid).zside() << " depth: " << HcalDetId(newid).depth()
                                   << " ieta: " << HcalDetId(newid).ietaAbs() << " iphi: " << HcalDetId(newid).iphi()
                                   << " wt " << wt;
       ++ii;
 #endif
     }
   } else {
     edm::LogWarning("HcalSim") << "HcalHitRelabeller: no valid HcalDDDRecConstants";
   }
 }

◆ relabel() [1/2]

DetId HcalHitRelabeller::relabel ( const uint32_t testId ) const

Definition at line 56 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 60 of file HcalHitRelabeller.cc.

References TauDecayModes::dec, hcalRecHitTable_cff::depth, PVValHelper::eta, HcalDDDRecConstants::getHCID(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, nano_mu_digi_cff::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);
 
   sign = (z == 0) ? (-1) : (1);
   HcalDDDRecConstants::HcalID id = theRecNumber->getHCID(det, sign * eta, phi, layer, depth);
   int depth0 = ((det == 1) && (layer == 2) && (eta == 15)) ? depth : id.depth;
 
   if (id.subdet == int(HcalBarrel)) {
     hid = HcalDetId(HcalBarrel, sign * id.eta, id.phi, depth0);
   } else if (id.subdet == int(HcalEndcap)) {
     hid = HcalDetId(HcalEndcap, sign * id.eta, id.phi, depth0);
   } else if (id.subdet == int(HcalOuter)) {
     hid = HcalDetId(HcalOuter, sign * id.eta, id.phi, depth0);
   } else if (id.subdet == int(HcalForward)) {
     hid = HcalDetId(HcalForward, sign * id.eta, id.phi, depth0);
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HcalSim") << "Initial indices:"
                               << "subdet: " << det << " "
                               << "z: " << z << " "
                               << "depth: " << depth << " "
                               << "ieta: " << eta << " "
                               << "iphi: " << phi << " "
                               << "layer: " << layer << " new HcalDetId -> hex.RawID = " << std::hex << hid.rawId()
                               << std::dec << " subdet, z, depth, eta, phi = " << det << " " << z << " " << id.depth
                               << " " << id.eta << " " << id.phi << " ---> " << hid;
 #endif
   return hid;
 }