#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::TrackIsoCalculator	(	reco::TrackCollection const &	trackCollection,
		std::string const &	trackQuality
	)

Definition at line 9 of file TrackIsoCalculator.cc.

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

References duplicaterechits_cfi::trackCollection.

Member Function Documentation

◆ getBkgSubTrackIso()

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.

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

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

Referenced by photonIsolationHIProducer::produce().

◆ getTrackIso()

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.

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