#include <PFJetMETcorrInputProducerT.h>

Inheritance diagram for PFJetMETcorrInputProducerT< T, Textractor >:

Classes
struct	type2BinningEntryType

Public Member Functions
	PFJetMETcorrInputProducerT (const edm::ParameterSet &cfg)

	~PFJetMETcorrInputProducerT () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::Event &evt, const edm::EventSetup &es) override

Private Attributes
double	jetCorrEtaMax_

Textractor	jetCorrExtractor_

edm::InputTag	jetCorrLabel_

edm::InputTag	jetCorrLabelRes_

edm::EDGetTokenT< reco::JetCorrector >	jetCorrResToken_

edm::EDGetTokenT< reco::JetCorrector >	jetCorrToken_

std::string	moduleLabel_

edm::InputTag	offsetCorrLabel_

edm::EDGetTokenT< reco::JetCorrector >	offsetCorrToken_

bool	skipEM_

double	skipEMfractionThreshold_

bool	skipMuons_

StringCutObjectSelector< reco::Candidate > *	skipMuonSelection_

edm::EDGetTokenT< std::vector< T > >	token_

double	type1JetPtThreshold_

std::vector< type2BinningEntryType * >	type2Binning_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

template<typename T, typename Textractor>
class PFJetMETcorrInputProducerT< T, Textractor >

Produce Type 1 + 2 MET corrections corresponding to differences between raw PFJets and PFJets with jet energy corrections (JECs) applied

NOTE: class is templated to that it works with reco::PFJets as well as with pat::Jets of PF-type as input

Authors: Michael Schmitt, Richard Cavanaugh, The University of Florida Florent Lacroix, University of Illinois at Chicago Christian Veelken, LLR

Definition at line 67 of file PFJetMETcorrInputProducerT.h.

Constructor & Destructor Documentation

template<typename T , typename Textractor >

PFJetMETcorrInputProducerT< T, Textractor >::PFJetMETcorrInputProducerT ( const edm::ParameterSet & cfg )

inlineexplicit

Definition at line 69 of file PFJetMETcorrInputProducerT.h.

References edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

       : moduleLabel_(cfg.getParameter<std::string>("@module_label")),
         offsetCorrLabel_(""),
         skipMuonSelection_(nullptr) {
     token_ = consumes<std::vector<T> >(cfg.getParameter<edm::InputTag>("src"));
 
     if (cfg.exists("offsetCorrLabel")) {
       offsetCorrLabel_ = cfg.getParameter<edm::InputTag>("offsetCorrLabel");
       offsetCorrToken_ = consumes<reco::JetCorrector>(offsetCorrLabel_);
     }
     jetCorrLabel_ = cfg.getParameter<edm::InputTag>("jetCorrLabel");        //for MC
     jetCorrLabelRes_ = cfg.getParameter<edm::InputTag>("jetCorrLabelRes");  //for data
     jetCorrToken_ = mayConsume<reco::JetCorrector>(jetCorrLabel_);
     jetCorrResToken_ = mayConsume<reco::JetCorrector>(jetCorrLabelRes_);
 
     jetCorrEtaMax_ = (cfg.exists("jetCorrEtaMax")) ? cfg.getParameter<double>("jetCorrEtaMax") : 9.9;
 
     type1JetPtThreshold_ = cfg.getParameter<double>("type1JetPtThreshold");
 
     skipEM_ = cfg.getParameter<bool>("skipEM");
     if (skipEM_) {
       skipEMfractionThreshold_ = cfg.getParameter<double>("skipEMfractionThreshold");
     }
 
     skipMuons_ = cfg.getParameter<bool>("skipMuons");
     if (skipMuons_) {
       std::string skipMuonSelection_string = cfg.getParameter<std::string>("skipMuonSelection");
       skipMuonSelection_ = new StringCutObjectSelector<reco::Candidate>(skipMuonSelection_string, true);
     }
 
     if (cfg.exists("type2Binning")) {
       typedef std::vector<edm::ParameterSet> vParameterSet;
       vParameterSet cfgType2Binning = cfg.getParameter<vParameterSet>("type2Binning");
       for (vParameterSet::const_iterator cfgType2BinningEntry = cfgType2Binning.begin();
            cfgType2BinningEntry != cfgType2Binning.end();
            ++cfgType2BinningEntry) {
         type2Binning_.push_back(new type2BinningEntryType(*cfgType2BinningEntry));
       }
     } else {
       type2Binning_.push_back(new type2BinningEntryType());
     }
 
     produces<CorrMETData>("type1");
     for (typename std::vector<type2BinningEntryType*>::const_iterator type2BinningEntry = type2Binning_.begin();
          type2BinningEntry != type2Binning_.end();
          ++type2BinningEntry) {
       produces<CorrMETData>((*type2BinningEntry)->getInstanceLabel_full("type2"));
       produces<CorrMETData>((*type2BinningEntry)->getInstanceLabel_full("offset"));
     }
   }

template<typename T , typename Textractor >

PFJetMETcorrInputProducerT< T, Textractor >::~PFJetMETcorrInputProducerT ( )

inlineoverride

Definition at line 119 of file PFJetMETcorrInputProducerT.h.

                                          {
     delete skipMuonSelection_;
 
     for (typename std::vector<type2BinningEntryType*>::const_iterator it = type2Binning_.begin();
          it != type2Binning_.end();
          ++it) {
       delete (*it);
     }
   }

Member Function Documentation

template<typename T , typename Textractor >

static void PFJetMETcorrInputProducerT< T, Textractor >::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

inlinestatic

Definition at line 129 of file PFJetMETcorrInputProducerT.h.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), defaultModuleLabel(), HLT_2018_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                            {
     edm::ParameterSetDescription desc;
     desc.add<edm::InputTag>("src", edm::InputTag("ak4PFJetsCHS"));
     desc.add<edm::InputTag>("offsetCorrLabel", edm::InputTag("ak4PFCHSL1FastjetCorrector"));
     desc.add<edm::InputTag>("jetCorrLabel", edm::InputTag("ak4PFCHSL1FastL2L3Corrector"));
     desc.add<edm::InputTag>("jetCorrLabelRes", edm::InputTag("ak4PFCHSL1FastL2L3ResidualCorrector"));
     desc.add<double>("jetCorrEtaMax", 9.9);
     desc.add<double>("type1JetPtThreshold", 15);
     desc.add<bool>("skipEM", true);
     desc.add<double>("skipEMfractionThreshold", 0.90);
     desc.add<bool>("skipMuons", true);
     desc.add<std::string>("skipMuonSelection", "isGlobalMuon | isStandAloneMuon");
     descriptions.add(defaultModuleLabel<PFJetMETcorrInputProducerT<T, Textractor> >(), desc);
   }

template<typename T , typename Textractor >

void PFJetMETcorrInputProducerT< T, Textractor >::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

inlineoverrideprivate

Definition at line 145 of file PFJetMETcorrInputProducerT.h.

References HLT_2018_cff::cands, edm::Ref< C, T, F >::get(), edm::Event::getByToken(), edm::Ref< C, T, F >::isNonnull(), edm::EventBase::isRealData(), metsig::jet, singleTopDQM_cfi::jets, eostools::move(), amptDefaultParameters_cff::mu, reco::PFCandidate::muonRef(), pfDeepBoostedJetPreprocessParams_cfi::pfcand, edm::Handle< T >::product(), and edm::Event::put().

                                                                 {
     std::unique_ptr<CorrMETData> type1Correction(new CorrMETData());
     for (typename std::vector<type2BinningEntryType*>::iterator type2BinningEntry = type2Binning_.begin();
          type2BinningEntry != type2Binning_.end();
          ++type2BinningEntry) {
       (*type2BinningEntry)->binUnclEnergySum_ = CorrMETData();
       (*type2BinningEntry)->binOffsetEnergySum_ = CorrMETData();
     }
 
     edm::Handle<reco::JetCorrector> jetCorr;
     //automatic switch for residual corrections
     if (evt.isRealData()) {
       jetCorrLabel_ = jetCorrLabelRes_;
       evt.getByToken(jetCorrResToken_, jetCorr);
     } else {
       evt.getByToken(jetCorrToken_, jetCorr);
     }
 
     typedef std::vector<T> JetCollection;
     edm::Handle<JetCollection> jets;
     evt.getByToken(token_, jets);
 
     int numJets = jets->size();
     for (int jetIndex = 0; jetIndex < numJets; ++jetIndex) {
       const T& jet = jets->at(jetIndex);
 
       const static PFJetMETcorrInputProducer_namespace::InputTypeCheckerT<T, Textractor> checkInputType{};
       checkInputType(jet);
 
       double emEnergyFraction = jet.chargedEmEnergyFraction() + jet.neutralEmEnergyFraction();
       if (skipEM_ && emEnergyFraction > skipEMfractionThreshold_)
         continue;
 
       const static PFJetMETcorrInputProducer_namespace::RawJetExtractorT<T> rawJetExtractor{};
       reco::Candidate::LorentzVector rawJetP4 = rawJetExtractor(jet);
 
       if (skipMuons_) {
         const std::vector<reco::CandidatePtr>& cands = jet.daughterPtrVector();
         for (std::vector<reco::CandidatePtr>::const_iterator cand = cands.begin(); cand != cands.end(); ++cand) {
           const reco::PFCandidate* pfcand = dynamic_cast<const reco::PFCandidate*>(cand->get());
           const reco::Candidate* mu =
               (pfcand != nullptr ? (pfcand->muonRef().isNonnull() ? pfcand->muonRef().get() : nullptr) : cand->get());
           if (mu != nullptr && (*skipMuonSelection_)(*mu)) {
             reco::Candidate::LorentzVector muonP4 = (*cand)->p4();
             rawJetP4 -= muonP4;
           }
         }
       }
 
       reco::Candidate::LorentzVector corrJetP4;
       if (checkInputType.isPatJet(jet))
         corrJetP4 = jetCorrExtractor_(jet, jetCorrLabel_.label(), jetCorrEtaMax_, &rawJetP4);
       else
         corrJetP4 = jetCorrExtractor_(jet, jetCorr.product(), jetCorrEtaMax_, &rawJetP4);
 
       if (corrJetP4.pt() > type1JetPtThreshold_) {
         reco::Candidate::LorentzVector rawJetP4offsetCorr = rawJetP4;
         if (!offsetCorrLabel_.label().empty()) {
           edm::Handle<reco::JetCorrector> offsetCorr;
           evt.getByToken(offsetCorrToken_, offsetCorr);
           if (checkInputType.isPatJet(jet))
             rawJetP4offsetCorr = jetCorrExtractor_(jet, offsetCorrLabel_.label(), jetCorrEtaMax_, &rawJetP4);
           else
             rawJetP4offsetCorr = jetCorrExtractor_(jet, offsetCorr.product(), jetCorrEtaMax_, &rawJetP4);
 
           for (typename std::vector<type2BinningEntryType*>::iterator type2BinningEntry = type2Binning_.begin();
                type2BinningEntry != type2Binning_.end();
                ++type2BinningEntry) {
             if (!(*type2BinningEntry)->binSelection_ || (*(*type2BinningEntry)->binSelection_)(corrJetP4)) {
               (*type2BinningEntry)->binOffsetEnergySum_.mex += (rawJetP4.px() - rawJetP4offsetCorr.px());
               (*type2BinningEntry)->binOffsetEnergySum_.mey += (rawJetP4.py() - rawJetP4offsetCorr.py());
               (*type2BinningEntry)->binOffsetEnergySum_.sumet += (rawJetP4.Et() - rawJetP4offsetCorr.Et());
             }
           }
         }
 
         //--- MET balances momentum of reconstructed particles,
         //    hence correction to jets and corresponding Type 1 MET correction are of opposite sign
         type1Correction->mex -= (corrJetP4.px() - rawJetP4offsetCorr.px());
         type1Correction->mey -= (corrJetP4.py() - rawJetP4offsetCorr.py());
         type1Correction->sumet += (corrJetP4.Et() - rawJetP4offsetCorr.Et());
       } else {
         for (typename std::vector<type2BinningEntryType*>::iterator type2BinningEntry = type2Binning_.begin();
              type2BinningEntry != type2Binning_.end();
              ++type2BinningEntry) {
           if (!(*type2BinningEntry)->binSelection_ || (*(*type2BinningEntry)->binSelection_)(corrJetP4)) {
             (*type2BinningEntry)->binUnclEnergySum_.mex += rawJetP4.px();
             (*type2BinningEntry)->binUnclEnergySum_.mey += rawJetP4.py();
             (*type2BinningEntry)->binUnclEnergySum_.sumet += rawJetP4.Et();
           }
         }
       }
     }
 
     //--- add
     //     o Type 1 MET correction                (difference corrected-uncorrected jet energy for jets of (corrected) Pt > 10 GeV)
     //     o momentum sum of "unclustered energy" (jets of (corrected) Pt < 10 GeV)
     //     o momentum sum of "offset energy"      (sum of energy attributed to pile-up/underlying event)
     //    to the event
     evt.put(std::move(type1Correction), "type1");
     for (typename std::vector<type2BinningEntryType*>::const_iterator type2BinningEntry = type2Binning_.begin();
          type2BinningEntry != type2Binning_.end();
          ++type2BinningEntry) {
       evt.put(std::unique_ptr<CorrMETData>(new CorrMETData((*type2BinningEntry)->binUnclEnergySum_)),
               (*type2BinningEntry)->getInstanceLabel_full("type2"));
       evt.put(std::unique_ptr<CorrMETData>(new CorrMETData((*type2BinningEntry)->binOffsetEnergySum_)),
               (*type2BinningEntry)->getInstanceLabel_full("offset"));
     }
   }