HcalRechitIsoCalculator Class Reference

#include <HcalRechitIsoCalculator.h>

Public Member Functions
double	getBkgSubHcalRechitIso (const reco::SuperClusterRef clus, const double i, const double threshold, const double innerR=0.0)
	Return the background-subtracted hcal rechit energy in a cone around the SC. More...
double	getHcalRechitIso (const reco::SuperClusterRef clus, const double i, const double threshold, const double innerR=0.0)
	Return the hcal rechit energy in a cone around the SC. More...
	HcalRechitIsoCalculator (const edm::Event &iEvent, const edm::EventSetup &iSetup, const edm::Handle< HBHERecHitCollection > hbhe, const edm::Handle< HFRecHitCollection > hfLabel, const edm::Handle< HORecHitCollection > hoLabel)

Private Attributes
const HBHERecHitCollection *	fHBHERecHits_
const HFRecHitCollection *	fHFRecHits_
const HORecHitCollection *	fHORecHits_
const CaloGeometry *	geometry_

Detailed Description

Definition at line 19 of file HcalRechitIsoCalculator.h.

Constructor & Destructor Documentation

HcalRechitIsoCalculator()

HcalRechitIsoCalculator::HcalRechitIsoCalculator	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const edm::Handle< HBHERecHitCollection >	hbhe,
		const edm::Handle< HFRecHitCollection >	hfLabel,
		const edm::Handle< HORecHitCollection >	hoLabel
	)

Definition at line 16 of file HcalRechitIsoCalculator.cc.

                                                                                         {
  if (hf.isValid())
    fHFRecHits_ = hf.product();
  else
    fHFRecHits_ = nullptr;
 
  if (ho.isValid())
    fHORecHits_ = ho.product();
  else
    fHORecHits_ = nullptr;
 
  if (hbhe.isValid())
    fHBHERecHits_ = hbhe.product();
  else
    fHBHERecHits_ = nullptr;
 
  ESHandle<CaloGeometry> geometryHandle;
  iSetup.get<CaloGeometryRecord>().get(geometryHandle);
  if (geometryHandle.isValid())
    geometry_ = geometryHandle.product();
  else
    geometry_ = nullptr;
}

References edm::EventSetup::get(), edm::get(), photonIsolationHIProducer_cfi::hbhe, es_hardcode_cfi::hf, es_hardcode_cfi::ho, edm::ESHandleBase::isValid(), and edm::ESHandle< T >::product().

Member Function Documentation

getBkgSubHcalRechitIso()

double HcalRechitIsoCalculator::getBkgSubHcalRechitIso	(	const reco::SuperClusterRef	clus,
		const double	i,
		const double	threshold,
		const double	innerR = 0.0
	)

Return the background-subtracted hcal rechit energy in a cone around the SC.

Definition at line 76 of file HcalRechitIsoCalculator.cc.

                                                                            {
  if (!fHBHERecHits_) {
    return -100;
  }
 
  double SClusterEta = cluster->eta();
  double TotalEt = 0;
 
  for (size_t index = 0; index < fHBHERecHits_->size(); index++) {
    const HBHERecHit &rechit = (*fHBHERecHits_)[index];
    const DetId &detid = rechit.id();
    const GlobalPoint &hitpoint = geometry_->getPosition(detid);
    double eta = hitpoint.eta();
    double dEta = fabs(eta - SClusterEta);
 
    if (dEta < x * 0.1) {
      double et = rechit.energy() / cosh(eta);
      if (et < threshold)
        et = 0;
      TotalEt += et;
    }
  }
 
  double Rx = getHcalRechitIso(cluster, x, threshold, innerR);
  double CRx = (Rx - TotalEt * (0.01 * x * x - innerR * innerR) / (2 * 2 * 0.1 * x)) * (1 / (1 - x / 40.));
 
  return CRx;
}

References HLT_2018_cff::dEta, CaloRecHit::energy(), EgHLTOffHistBins_cfi::et, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), HBHERecHit::id(), and remoteMonitoring_LED_IterMethod_cfg::threshold.

Referenced by photonIsolationHIProducer::produce().

getHcalRechitIso()

double HcalRechitIsoCalculator::getHcalRechitIso	(	const reco::SuperClusterRef	clus,
		const double	i,
		const double	threshold,
		const double	innerR = 0.0
	)

Return the hcal rechit energy in a cone around the SC.

Definition at line 44 of file HcalRechitIsoCalculator.cc.

                                                                      {
  if (!fHBHERecHits_) {
    return -100;
  }
 
  double TotalEt = 0;
 
  for (size_t index = 0; index < fHBHERecHits_->size(); index++) {
    const HBHERecHit &rechit = (*fHBHERecHits_)[index];
    const DetId &detid = rechit.id();
    const GlobalPoint &hitpoint = geometry_->getPosition(detid);
    double eta = hitpoint.eta();
 
    double dR2 = reco::deltaR2(*cluster, hitpoint);
    // veto inner cone///////////////
    if (dR2 < innerR * innerR)
      continue;
    if (dR2 < (x * x * 0.01)) {
      double et = rechit.energy() / cosh(eta);
      if (et < threshold)
        et = 0;
      TotalEt += et;
    }
  }
 
  return TotalEt;
}

References reco::deltaR2(), CaloRecHit::energy(), EgHLTOffHistBins_cfi::et, PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), HBHERecHit::id(), and remoteMonitoring_LED_IterMethod_cfg::threshold.

Member Data Documentation

fHBHERecHits_

const HBHERecHitCollection* HcalRechitIsoCalculator::fHBHERecHits_

private

Definition at line 39 of file HcalRechitIsoCalculator.h.

fHFRecHits_

const HFRecHitCollection* HcalRechitIsoCalculator::fHFRecHits_

private

Definition at line 41 of file HcalRechitIsoCalculator.h.

fHORecHits_

const HORecHitCollection* HcalRechitIsoCalculator::fHORecHits_

private

Definition at line 40 of file HcalRechitIsoCalculator.h.

geometry_

const CaloGeometry* HcalRechitIsoCalculator::geometry_

private

Definition at line 42 of file HcalRechitIsoCalculator.h.