PATJetCorrExtractor Class Reference

#include <PATJetCorrExtractor.h>

Public Member Functions
reco::Candidate::LorentzVector	operator() (const pat::Jet rawJet, const reco::JetCorrector *jetCorr, double jetCorrEtaMax=9.9, const reco::Candidate::LorentzVector *const rawJetP4_specified=0) const
reco::Candidate::LorentzVector	operator() (const pat::Jet &jet, const std::string &jetCorrLabel, double jetCorrEtaMax=9.9, const reco::Candidate::LorentzVector *const rawJetP4_specified=0) const

Detailed Description

Retrieve jet energy correction factor for pat::Jets (of either PF-type or Calo-type)

NOTE: this specialization of the "generic" template defined in JetMETCorrections/Type1MET/interface/JetCorrExtractorT.h is to be used for pat::Jets only

Author

Christian Veelken, LLR

Definition at line 51 of file PATJetCorrExtractor.h.

Member Function Documentation

reco::Candidate::LorentzVector PATJetCorrExtractor::operator() ( const pat::Jet  rawJet, const reco::JetCorrector *  jetCorr, double  jetCorrEtaMax = 9.9, const reco::Candidate::LorentzVector *const  rawJetP4_specified = 0  ) const

inline

Definition at line 53 of file PATJetCorrExtractor.h.

References HLT_2018_cff::jetCorrEtaMax.

                                                                                  {
    JetCorrExtractorT<pat::Jet> jetCorrExtractor;
    return jetCorrExtractor(rawJet, jetCorr, jetCorrEtaMax, rawJetP4_specified);
  }

reco::Candidate::LorentzVector PATJetCorrExtractor::operator() ( const pat::Jet &  jet, const std::string &  jetCorrLabel, double  jetCorrEtaMax = 9.9, const reco::Candidate::LorentzVector *const  rawJetP4_specified = 0  ) const

inline

Definition at line 62 of file PATJetCorrExtractor.h.

References pat::Jet::availableJECLevels(), pat::Jet::correctedP4(), Exception, and format_vstring().

                                                                                  {
    reco::Candidate::LorentzVector corrJetP4;

    try {
      corrJetP4 = jet.correctedP4(jetCorrLabel);
      if (rawJetP4_specified != nullptr) {
        //MM: compensate for potential removal of constituents (as muons)
        //similar effect in JetMETCorrection/Type1MET/interface/JetCorrExtractor.h
        reco::Candidate::LorentzVector rawJetP4 = jet.correctedP4("Uncorrected");
        double corrFactor = corrJetP4.pt() / rawJetP4.pt();
        corrJetP4 = (*rawJetP4_specified);
        corrJetP4 *= corrFactor;
        if (corrFactor < 0) {
          edm::LogWarning("PATJetCorrExtractor") << "Negative jet energy scale correction noticed"
                                                 << ".\n";
        }
      }
    } catch (cms::Exception const&) {
      throw cms::Exception("InvalidRequest")
          << "The JEC level " << jetCorrLabel << " does not exist !!\n"
          << "Available levels = " << format_vstring(jet.availableJECLevels()) << ".\n";
    }

    return corrJetP4;
  }