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 CaloGeometry *geometry, 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	(	const CaloGeometry *	geometry,
		const edm::Handle< HBHERecHitCollection >	hbhe,
		const edm::Handle< HFRecHitCollection >	hfLabel,
		const edm::Handle< HORecHitCollection >	hoLabel
	)

Definition at line 14 of file HcalRechitIsoCalculator.cc.

References relativeConstraints::geometry, edm::HandleBase::isValid(), and edm::Handle< T >::product().

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

  geometry_ = geometry;
}

Member Function Documentation

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 68 of file HcalRechitIsoCalculator.cc.

References CaloRecHit::energy(), PVValHelper::eta, PV3DBase< T, PVType, FrameType >::eta(), and HBHERecHit::id().

Referenced by photonIsolationHIProducer::produce().

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

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 36 of file HcalRechitIsoCalculator.cc.

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

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

Member Data Documentation

const HBHERecHitCollection* HcalRechitIsoCalculator::fHBHERecHits_

private

Definition at line 38 of file HcalRechitIsoCalculator.h.

const HFRecHitCollection* HcalRechitIsoCalculator::fHFRecHits_

private

Definition at line 40 of file HcalRechitIsoCalculator.h.

const HORecHitCollection* HcalRechitIsoCalculator::fHORecHits_

private

Definition at line 39 of file HcalRechitIsoCalculator.h.

const CaloGeometry* HcalRechitIsoCalculator::geometry_

private

Definition at line 41 of file HcalRechitIsoCalculator.h.