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=nullptr) const
reco::Candidate::LorentzVector	operator() (const T &rawJet, const std::string &jetCorrLabel, double jetCorrEtaMax=9.9, const reco::Candidate::LorentzVector *const rawJetP4_specified=nullptr) 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

getCorrection()

template<typename T>

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

inlinestaticprivate

Definition at line 92 of file JetCorrExtractorT.h.

{ return jetCorr->correction(jet); }

References reco::JetCorrector::correction(), and metsig::jet.

Referenced by JetReCalibrator.JetReCalibrator::correct(), and JetCorrExtractorT< T >::operator()().

operator()() [1/2]

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 = nullptr  ) const

inline

Definition at line 44 of file JetCorrExtractorT.h.

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

References JetCorrExtractorT< T >::getCorrection(), HLT_FULL_cff::jetCorrEtaMax, and jetcorrextractor::sign().

operator()() [2/2]

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 = nullptr  ) 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;
  }