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

References edm::EventSetup::get(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), edm::Handle< T >::product(), and edm::ESHandle< class >::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;

  ESHandle<CaloGeometry> geometryHandle;
  iSetup.get<CaloGeometryRecord>().get(geometryHandle);
  if(geometryHandle.isValid())
    geometry_ = geometryHandle.product();
  else
    geometry_ = nullptr;

}

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

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

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

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

{
   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 35 of file HcalRechitIsoCalculator.h.

const HFRecHitCollection* HcalRechitIsoCalculator::fHFRecHits_

private

Definition at line 37 of file HcalRechitIsoCalculator.h.

const HORecHitCollection* HcalRechitIsoCalculator::fHORecHits_

private

Definition at line 36 of file HcalRechitIsoCalculator.h.

const CaloGeometry* HcalRechitIsoCalculator::geometry_

private

Definition at line 38 of file HcalRechitIsoCalculator.h.