pat::CalculatePtRatioRel Class Reference

#include <CalculatePtRatioRel.h>

Public Member Functions
	CalculatePtRatioRel (float dR2max)
std::array< float, 2 >	computePtRatioRel (const pat::Muon &imuon, const reco::JetTagCollection &bTags, const reco::JetCorrector *correctorL1=nullptr, const reco::JetCorrector *correctorL1L2L3Res=nullptr) const
	~CalculatePtRatioRel ()

Private Attributes
float	dR2max_

Detailed Description

Definition at line 18 of file CalculatePtRatioRel.h.

Constructor & Destructor Documentation

CalculatePtRatioRel()

CalculatePtRatioRel::CalculatePtRatioRel

(

float

dR2max

)

Definition at line 10 of file CalculatePtRatioRel.cc.

: dR2max_(dR2max) {}

~CalculatePtRatioRel()

CalculatePtRatioRel::~CalculatePtRatioRel

(

)

Definition at line 12 of file CalculatePtRatioRel.cc.

{}

Member Function Documentation

computePtRatioRel()

std::array< float, 2 > CalculatePtRatioRel::computePtRatioRel	(	const pat::Muon &	imuon,
		const reco::JetTagCollection &	bTags,
		const reco::JetCorrector *	correctorL1 = nullptr,
		const reco::JetCorrector *	correctorL1L2L3Res = nullptr
	)		const

Definition at line 29 of file CalculatePtRatioRel.cc.

References TtFullHadKinFitProducer_cfi::bTags, reco::JetCorrector::correction(), HLTMuonOfflineAnalyzer_cfi::deltaR2, dR2max_, f, SiStripPI::max, PatBasicFWLiteJetAnalyzer_Selector_cfg::outputs, and electrons_cff::ptRel.

                                                                                                              {
  //Initialise loop variables
  float minDr2 = 9999;
  double jecL1L2L3Res = 1.;
  double jecL1 = 1.;

  // Compute corrected isolation variables
  const float chIso = muon.pfIsolationR04().sumChargedHadronPt;
  const float nIso = muon.pfIsolationR04().sumNeutralHadronEt;
  const float phoIso = muon.pfIsolationR04().sumPhotonEt;
  const float puIso = muon.pfIsolationR04().sumPUPt;
  const float dbCorrectedIsolation = chIso + std::max(nIso + phoIso - .5f * puIso, 0.f);
  const float dbCorrectedRelIso = dbCorrectedIsolation / muon.pt();

  float JetPtRatio = 1. / (1 + dbCorrectedRelIso);
  float JetPtRel = 0.;

  const reco::Candidate::LorentzVector& muP4(muon.p4());
  for (const auto& tagI : bTags) {
    // for each muon with the lepton
    float dr2 = deltaR2(*(tagI.first), muon);
    if (dr2 > minDr2)
      continue;
    minDr2 = dr2;

    reco::Candidate::LorentzVector jetP4(tagI.first->p4());

    if (correctorL1 && correctorL1L2L3Res) {
      jecL1L2L3Res = correctorL1L2L3Res->correction(*(tagI.first));
      jecL1 = correctorL1->correction(*(tagI.first));
    }

    if (minDr2 < dR2max_) {
      if ((jetP4 - muP4).Rho() < 0.0001) {
        JetPtRel = 0.;
        JetPtRatio = 1.;
      } else {
        jetP4 -= muP4 / jecL1;
        jetP4 *= jecL1L2L3Res;
        jetP4 += muP4;

        JetPtRatio = muP4.pt() / jetP4.pt();
        JetPtRel = ptRel(muP4, jetP4);
      }
    }
  }

  if (JetPtRatio > 1.5)
    JetPtRatio = 1.5;

  std::array<float, 2> outputs = {{JetPtRatio, JetPtRel}};
  return outputs;
};

Member Data Documentation

dR2max_

float pat::CalculatePtRatioRel::dR2max_

private

Definition at line 30 of file CalculatePtRatioRel.h.

Referenced by computePtRatioRel().