HcalHitRelabeller Class Reference

#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

(

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

double HcalHitRelabeller::energyWt ( const uint32_t  testId ) const

private

Definition at line 101 of file HcalHitRelabeller.cc.

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

Referenced by process().

                                                              {

  HcalDetId hid;
  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;
}

void HcalHitRelabeller::process

(

std::vector< PCaloHit > &

hcalHits

)

Definition at line 18 of file HcalHitRelabeller.cc.

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

Referenced by HcalDigitizer::accumulateCaloHits(), ConfigBuilder.ConfigBuilder::addExtraStream(), ConfigBuilder.ConfigBuilder::completeInputCommand(), ConfigBuilder.ConfigBuilder::doNotInlineEventContent(), ConfigBuilder.ConfigBuilder.PrintAllModules::leave(), ConfigBuilder.ConfigBuilder::prepare(), ConfigBuilder.ConfigBuilder::prepare_ALCA(), ConfigBuilder.ConfigBuilder::prepare_DQM(), ConfigBuilder.ConfigBuilder::prepare_HLT(), ConfigBuilder.ConfigBuilder::prepare_LHE(), ConfigBuilder.ConfigBuilder::prepare_PATFILTER(), ConfigBuilder.ConfigBuilder::prepare_VALIDATION(), ConfigBuilder.ConfigBuilder::renameHLTprocessInSequence(), ConfigBuilder.ConfigBuilder::renameInputTagsInSequence(), and ConfigBuilder.ConfigBuilder::scheduleSequence().

                                                             {

  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";
  }
  
}

DetId HcalHitRelabeller::relabel

(

const uint32_t

testId

)

const

Definition at line 59 of file HcalHitRelabeller.cc.

References theRecNumber.

Referenced by SimHitsValidationHcal::analyze(), HcalSimHitsValidation::analyze(), HcalHBHEMuonSimAnalyzer::analyze(), HcalRecHitsValidation::analyze(), HGCalSimHitValidation::analyze(), HcalSimHitStudy::analyzeHits(), HcalTestNS::compare(), GlobalHitsAnalyzer::fillHCal(), CaloTruthAccumulator::fillSimHits(), process(), HcalDigisValidation::reco(), HGCalTriggerNtupleHGCTriggerCells::simhits(), and HGCalTriggerNtupleHGCDigis::simhits().

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

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

static

Definition at line 63 of file HcalHitRelabeller.cc.

References TauDecayModes::dec, egammaForCoreTracking_cff::depth, PVValHelper::eta, HcalDDDRecConstants::getHCID(), HcalBarrel, HcalEndcap, HcalForward, HcalOuter, phi, DetId::rawId(), Validation_hcalonly_cfi::sign, 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;
}

void HcalHitRelabeller::setGeometry

(

const HcalDDDRecConstants *&

recNum

)

Definition at line 55 of file HcalHitRelabeller.cc.

References theRecNumber.

Referenced by HcalDigitizer::checkGeometry().

                                                                       {
  theRecNumber = recNum;
}

Member Data Documentation

bool HcalHitRelabeller::neutralDensity_

private

Definition at line 22 of file HcalHitRelabeller.h.

Referenced by HcalHitRelabeller(), and process().

const HcalDDDRecConstants* HcalHitRelabeller::theRecNumber

private

Definition at line 21 of file HcalHitRelabeller.h.

Referenced by energyWt(), process(), relabel(), and setGeometry().