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 41 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 93 of file JetCorrExtractorT.h.

References reco::JetCorrector::correction().

Referenced by JetCorrExtractorT< T >::operator()().

  {
    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 45 of file JetCorrExtractorT.h.

References JetCorrExtractorT< T >::getCorrection(), HLT_25ns15e33_v4_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 79 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;
  }