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 HLT_2018_cff::trackLabel, and StandaloneTrackMonitor_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 52 of file TrackIsoCalculator.cc.

References reco::deltaR2(), HLT_2018_cff::dEta, reco::LeafCandidate::eta(), HLT_2018_cff::eta2, DiDispStaMuonMonitor_cfi::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 27 of file TrackIsoCalculator.cc.

References reco::deltaR2(), DiDispStaMuonMonitor_cfi::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.