EcalClusterIsoCalculator Class Reference

#include <EcalClusterIsoCalculator.h>

Public Member Functions
	EcalClusterIsoCalculator (const edm::Event &iEvent, const edm::EventSetup &iSetup, const edm::Handle< reco::BasicClusterCollection > barrel, const edm::Handle< reco::BasicClusterCollection > endcap)
double	getBkgSubEcalClusterIso (const reco::SuperClusterRef clus, const double radius, const double threshold)
	Return the background-subtracted ecal cluster energy in a cone around the SC. More...
double	getEcalClusterIso (const reco::SuperClusterRef clus, const double radius, const double threshold)
	Return the ecal cluster energy in a cone around the SC. More...

Private Attributes
const reco::BasicClusterCollection *	fEBclusters_
const reco::BasicClusterCollection *	fEEclusters_
const CaloGeometry *	geometry_

Detailed Description

Definition at line 17 of file EcalClusterIsoCalculator.h.

Constructor & Destructor Documentation

EcalClusterIsoCalculator()

EcalClusterIsoCalculator::EcalClusterIsoCalculator	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const edm::Handle< reco::BasicClusterCollection >	barrel,
		const edm::Handle< reco::BasicClusterCollection >	endcap
	)

Definition at line 19 of file EcalClusterIsoCalculator.cc.

                                                                                                        {
  if (pEBclusters.isValid())
    fEBclusters_ = pEBclusters.product();
  else
    fEBclusters_ = nullptr;
 
  if (pEEclusters.isValid())
    fEEclusters_ = pEEclusters.product();
  else
    fEEclusters_ = nullptr;
 
  ESHandle<CaloGeometry> geometryHandle;
  iSetup.get<CaloGeometryRecord>().get(geometryHandle);
  if (geometryHandle.isValid())
    geometry_ = geometryHandle.product();
  else
    geometry_ = nullptr;
}

References edm::EventSetup::get(), edm::get(), edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), edm::Handle< T >::product(), and edm::ESHandle< T >::product().

Member Function Documentation

getBkgSubEcalClusterIso()

double EcalClusterIsoCalculator::getBkgSubEcalClusterIso	(	const reco::SuperClusterRef	clus,
		const double	radius,
		const double	threshold
	)

Return the background-subtracted ecal cluster energy in a cone around the SC.

Definition at line 93 of file EcalClusterIsoCalculator.cc.

                                                                                 {
  if (!fEBclusters_) {
    return -100;
  }
 
  if (!fEEclusters_) {
    return -100;
  }
 
  double SClusterEta = cluster->eta();
  double TotalEt = 0;
 
  TotalEt = -cluster->rawEnergy() / cosh(cluster->eta());
 
  for (BasicClusterCollection::const_iterator iclu = fEBclusters_->begin(); iclu != fEBclusters_->end(); ++iclu) {
    const BasicCluster *clu = &(*iclu);
    math::XYZVector ClusPoint(clu->x(), clu->y(), clu->z());
    double eta = ClusPoint.eta();
 
    double dEta = fabs(eta - SClusterEta);
 
    if (dEta < x * 0.1) {
      double et = clu->energy() / cosh(eta);
      if (et < threshold)
        et = 0;
      TotalEt += et;
    }
  }
 
  for (BasicClusterCollection::const_iterator iclu = fEEclusters_->begin(); iclu != fEEclusters_->end(); ++iclu) {
    const BasicCluster *clu = &(*iclu);
    math::XYZVector ClusPoint(clu->x(), clu->y(), clu->z());
    double eta = ClusPoint.eta();
    double dEta = fabs(eta - SClusterEta);
 
    if (dEta < x * 0.1) {
      double et = clu->energy() / cosh(eta);
      if (et < threshold)
        et = 0;
      TotalEt += et;
    }
  }
 
  double Cx = getEcalClusterIso(cluster, x, threshold);
  double CCx = (Cx - TotalEt / 40.0 * x) * (1 / (1 - x / 40.));
 
  return CCx;
}

References HLT_FULL_cff::dEta, reco::CaloCluster::energy(), EgHLTOffHistBins_cfi::et, PVValHelper::eta, remoteMonitoring_LED_IterMethod_cfg::threshold, reco::CaloCluster::x(), reco::CaloCluster::y(), and reco::CaloCluster::z().

Referenced by photonIsolationHIProducer::produce().

getEcalClusterIso()

double EcalClusterIsoCalculator::getEcalClusterIso	(	const reco::SuperClusterRef	clus,
		const double	radius,
		const double	threshold
	)

Return the ecal cluster energy in a cone around the SC.

Definition at line 41 of file EcalClusterIsoCalculator.cc.

                                                                           {
  if (!fEBclusters_) {
    return -100;
  }
 
  if (!fEEclusters_) {
    return -100;
  }
 
  math::XYZVector SClusPoint(cluster->position().x(), cluster->position().y(), cluster->position().z());
 
  double TotalEt = 0;
 
  TotalEt = -cluster->rawEnergy() / cosh(cluster->eta());
 
  // Loop over barrel basic clusters
  for (BasicClusterCollection::const_iterator iclu = fEBclusters_->begin(); iclu != fEBclusters_->end(); ++iclu) {
    const BasicCluster *clu = &(*iclu);
    math::XYZVector ClusPoint(clu->x(), clu->y(), clu->z());
    double eta = ClusPoint.eta();
 
    double dR2 = reco::deltaR2(*clu, *cluster);
 
    if (dR2 < (x * x * 0.01)) {
      double et = clu->energy() / cosh(eta);
      if (et < threshold)
        et = 0;
      TotalEt += et;
    }
  }
 
  for (BasicClusterCollection::const_iterator iclu = fEEclusters_->begin(); iclu != fEEclusters_->end(); ++iclu) {
    const BasicCluster *clu = &(*iclu);
    const GlobalPoint clusPoint(clu->x(), clu->y(), clu->z());
    math::XYZVector ClusPoint(clu->x(), clu->y(), clu->z());
    double eta = ClusPoint.eta();
 
    double dR2 = reco::deltaR2(*clu, *cluster);
 
    if (dR2 < (x * x * 0.01)) {
      double et = clu->energy() / cosh(eta);
      if (et < threshold)
        et = 0;
      TotalEt += et;
    }
  }
 
  return TotalEt;
}

References reco::deltaR2(), reco::CaloCluster::energy(), EgHLTOffHistBins_cfi::et, PVValHelper::eta, remoteMonitoring_LED_IterMethod_cfg::threshold, reco::CaloCluster::x(), reco::CaloCluster::y(), and reco::CaloCluster::z().

Member Data Documentation

fEBclusters_

const reco::BasicClusterCollection* EcalClusterIsoCalculator::fEBclusters_

private

Definition at line 30 of file EcalClusterIsoCalculator.h.

fEEclusters_

const reco::BasicClusterCollection* EcalClusterIsoCalculator::fEEclusters_

private

Definition at line 31 of file EcalClusterIsoCalculator.h.

geometry_

const CaloGeometry* EcalClusterIsoCalculator::geometry_

private

Definition at line 32 of file EcalClusterIsoCalculator.h.