JetCorrExtractorT< T > Class Template Reference

#include <JetCorrExtractorT.h>

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

Static Private Member Functions
static double	getCorrection (const T &jet, const reco::JetCorrector *jetCorr)

Detailed Description

template<typename T>
class JetCorrExtractorT< T >

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

NOTE: this "general" implementation is to be used for reco::CaloJets and reco::PFJets, not for pat::Jets

Author

Christian Veelken, LLR

Definition at line 42 of file JetCorrExtractorT.h.

Member Function Documentation

template<typename T>

static double JetCorrExtractorT< T >::getCorrection ( const T &  jet, const reco::JetCorrector *  jetCorr  )

inlinestaticprivate

Definition at line 92 of file JetCorrExtractorT.h.

References reco::JetCorrector::correction().

Referenced by JetReCalibrator.JetReCalibrator::correct().

{ return jetCorr->correction(jet); }

template<typename T>

reco::Candidate::LorentzVector JetCorrExtractorT< T >::operator() ( const T &  rawJet, const reco::JetCorrector *  jetCorr, double  jetCorrEtaMax = 9.9, const reco::Candidate::LorentzVector *const  rawJetP4_specified = 0  ) const

inline

Definition at line 44 of file JetCorrExtractorT.h.

References HLT_2018_cff::jetCorrEtaMax, and jetcorrextractor::sign().

                                                                                  {
    // allow to specify four-vector to be used as "raw" (uncorrected) jet momentum,
    // call 'rawJet.p4()' in case four-vector not specified explicitely
    reco::Candidate::LorentzVector rawJetP4 = (rawJetP4_specified) ? (*rawJetP4_specified) : rawJet.p4();

    double jetCorrFactor = 1.;
    if (fabs(rawJetP4.eta()) < jetCorrEtaMax) {
      jetCorrFactor = getCorrection(rawJet, jetCorr);
    } else {
      reco::Candidate::PolarLorentzVector modJetPolarP4(rawJetP4);
      modJetPolarP4.SetEta(jetcorrextractor::sign(rawJetP4.eta()) * jetCorrEtaMax);

      reco::Candidate::LorentzVector modJetP4(modJetPolarP4);

      T modJet(rawJet);
      modJet.setP4(modJetP4);

      jetCorrFactor = getCorrection(modJet, jetCorr);
      if (jetCorrFactor < 0) {
        edm::LogWarning("JetCorrExtractor") << "Negative jet energy scale correction noticed"
                                            << ".\n";
      }
    }

    reco::Candidate::LorentzVector corrJetP4 = rawJetP4;
    corrJetP4 *= jetCorrFactor;

    return corrJetP4;
  }

template<typename T>

reco::Candidate::LorentzVector JetCorrExtractorT< T >::operator() ( const T &  rawJet, const std::string &  jetCorrLabel, double  jetCorrEtaMax = 9.9, const reco::Candidate::LorentzVector *const  rawJetP4_specified = 0  ) const

inline

Definition at line 78 of file JetCorrExtractorT.h.

                                                                                  {
    edm::LogWarning("JetCorrExtractor")
        << "JetCorrExtractorT<T>::operator(const T&, const std::string&, ...) is deprecated.\n"
        << "Please use JetCorrExtractorT<T>::operator(const T&, const reco::JetCorrector*, ...) instead.\n"
        << "Jet remains uncorrected!";
    reco::Candidate::LorentzVector rawJetP4 = (rawJetP4_specified) ? (*rawJetP4_specified) : rawJet.p4();
    return rawJetP4;
  }