TrackIsoCalculator Class Reference

#include <TrackIsoCalculator.h>

Public Member Functions
double	getBkgSubTrackIso (const reco::Photon clus, const double i, const double threshold, const double innerDR=0)
	Return the background-subtracted tracker energy in a cone around the photon. More...
double	getTrackIso (const reco::Photon clus, const double i, const double threshold, const double innerDR=0)
	Return the tracker energy in a cone around the photon. More...
	TrackIsoCalculator (const edm::Event &iEvent, const edm::EventSetup &iSetup, const edm::Handle< reco::TrackCollection > trackLabel, const std::string trackQuality_)

Private Attributes
edm::Handle< reco::TrackCollection >	recCollection
std::string	trackQuality_

Detailed Description

Definition at line 22 of file TrackIsoCalculator.h.

Constructor & Destructor Documentation

TrackIsoCalculator::TrackIsoCalculator	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup,
		const edm::Handle< reco::TrackCollection >	trackLabel,
		const std::string	trackQuality_
	)

Definition at line 19 of file TrackIsoCalculator.cc.

References MuonErrorMatrixAdjuster_cfi::trackLabel, and photonIsolationHIProducer_cfi::trackQuality.

{
  recCollection = trackLabel;
  trackQuality_ = trackQuality;
}

Member Function Documentation

double TrackIsoCalculator::getBkgSubTrackIso	(	const reco::Photon	clus,
		const double	i,
		const double	threshold,
		const double	innerDR = 0
	)

Return the background-subtracted tracker energy in a cone around the photon.

Definition at line 48 of file TrackIsoCalculator.cc.

References reco::deltaR2(), particleFlow_cfi::dEta, reco::LeafCandidate::eta(), EnergyCorrector::pt, and reco::TrackBase::qualityByName().

Referenced by photonIsolationHIProducer::produce().

{
  double SClusterEta = cluster.eta();
  double TotalPt = 0;

  TotalPt = 0;

  for(reco::TrackCollection::const_iterator
    recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
  {
    bool goodtrack = recTrack->quality(reco::TrackBase::qualityByName(trackQuality_));
    if(!goodtrack) continue;

    double pt = recTrack->pt();
    double eta2 = recTrack->eta();
    double dEta = fabs(eta2-SClusterEta);
    double dR2 = reco::deltaR2(cluster, *recTrack);
    if(dEta >= 0.1 * x)
      continue;
    if(dR2 < innerDR*innerDR)
      continue;

    if(pt > threshold)
      TotalPt = TotalPt + pt;
  }

  double Tx = getTrackIso(cluster,x,threshold,innerDR);
  double CTx = (Tx - TotalPt / 40.0 * x)*(1/(1-x/40.));

  return CTx;
}

double TrackIsoCalculator::getTrackIso	(	const reco::Photon	clus,
		const double	i,
		const double	threshold,
		const double	innerDR = 0
	)

Return the tracker energy in a cone around the photon.

Definition at line 25 of file TrackIsoCalculator.cc.

References reco::deltaR2(), EnergyCorrector::pt, and reco::TrackBase::qualityByName().

{
  double TotalPt = 0;

  for(reco::TrackCollection::const_iterator
    recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
  {
    bool goodtrack = recTrack->quality(reco::TrackBase::qualityByName(trackQuality_));
    if(!goodtrack) continue;

    double pt = recTrack->pt();
    double dR2 = reco::deltaR2(cluster, *recTrack);
    if(dR2 >= (0.01 * x*x))
      continue;
    if(dR2 < innerDR*innerDR)
      continue;
    if(pt > threshold)
      TotalPt = TotalPt + pt;
  }

  return TotalPt;
}

Member Data Documentation

edm::Handle<reco::TrackCollection> TrackIsoCalculator::recCollection

private

Definition at line 35 of file TrackIsoCalculator.h.

std::string TrackIsoCalculator::trackQuality_

private

Definition at line 36 of file TrackIsoCalculator.h.