TrackIsoCalculator Class Reference

#include <TrackIsoCalculator.h>

Public Member Functions
double	getBkgSubTrackIso (reco::Photon const &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 (reco::Photon const &clus, const double i, const double threshold, const double innerDR=0)
	Return the tracker energy in a cone around the photon. More...
	TrackIsoCalculator (reco::TrackCollection const &trackCollection, std::string const &trackQuality)

Private Attributes
reco::TrackCollection const &	recCollection_
std::string const &	trackQuality_

Detailed Description

Definition at line 8 of file TrackIsoCalculator.h.

Constructor & Destructor Documentation

TrackIsoCalculator::TrackIsoCalculator	(	reco::TrackCollection const &	trackCollection,
		std::string const &	trackQuality
	)

Definition at line 9 of file TrackIsoCalculator.cc.

    : recCollection_{trackCollection}, trackQuality_{trackQuality} {}

Member Function Documentation

double TrackIsoCalculator::getBkgSubTrackIso	(	reco::Photon const &	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 36 of file TrackIsoCalculator.cc.

References funct::abs(), reco::deltaR2(), reco::LeafCandidate::eta(), getTrackIso(), DiDispStaMuonMonitor_cfi::pt, reco::TrackBase::qualityByName(), recCollection_, trackQuality_, and x.

Referenced by photonIsolationHIProducer::produce().

                                                                   {
  double SClusterEta = cluster.eta();
  double totalPt = 0.0;

  for (auto const& recTrack : recCollection_) {
    bool goodtrack = recTrack.quality(reco::TrackBase::qualityByName(trackQuality_));
    if (!goodtrack)
      continue;

    double pt = recTrack.pt();
    if (std::abs(recTrack.eta() - SClusterEta) >= 0.1 * x)
      continue;
    if (reco::deltaR2(cluster, recTrack) < 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	(	reco::Photon const &	clus,
		const double	i,
		const double	threshold,
		const double	innerDR = 0
	)

Return the tracker energy in a cone around the photon.

Definition at line 12 of file TrackIsoCalculator.cc.

References reco::deltaR2(), DiDispStaMuonMonitor_cfi::pt, reco::TrackBase::qualityByName(), recCollection_, and trackQuality_.

Referenced by getBkgSubTrackIso().

                                                             {
  double totalPt = 0;

  for (auto const& recTrack : recCollection_) {
    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

reco::TrackCollection const& TrackIsoCalculator::recCollection_

private

Definition at line 18 of file TrackIsoCalculator.h.

Referenced by getBkgSubTrackIso(), and getTrackIso().

std::string const& TrackIsoCalculator::trackQuality_

private

Definition at line 19 of file TrackIsoCalculator.h.

Referenced by getBkgSubTrackIso(), and getTrackIso().