#include <PFJetMETcorrInputProducerT.h>

Inheritance diagram for PFJetMETcorrInputProducerT< T, Textractor >:

Classes
struct	type2BinningEntryType

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

	~PFJetMETcorrInputProducerT ()

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

Private Attributes
double	jetCorrEtaMax_

Textractor	jetCorrExtractor_

edm::InputTag	jetCorrLabel_

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::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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 70 of file PFJetMETcorrInputProducerT.h.

Constructor & Destructor Documentation

template<typename T , typename Textractor >

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

inlineexplicit

Definition at line 74 of file PFJetMETcorrInputProducerT.h.

     : moduleLabel_(cfg.getParameter<std::string>("@module_label")),
       offsetCorrLabel_(""),
       skipMuonSelection_(0)
   {
     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");
     jetCorrToken_ = consumes<reco::JetCorrector>(jetCorrLabel_);
 
     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 ( )

inline

Definition at line 122 of file PFJetMETcorrInputProducerT.h.

References PFJetMETcorrInputProducerT< T, Textractor >::skipMuonSelection_, and PFJetMETcorrInputProducerT< T, Textractor >::type2Binning_.

   {
     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 >

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

inlineprivatevirtual

Implements edm::EDProducer.

Definition at line 134 of file PFJetMETcorrInputProducerT.h.

   {
     std::auto_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;
     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& rawJet = jets->at(jetIndex);
 
       static PFJetMETcorrInputProducer_namespace::InputTypeCheckerT<T, Textractor> checkInputType;
       checkInputType(rawJet);
 
       double emEnergyFraction = rawJet.chargedEmEnergyFraction() + rawJet.neutralEmEnergyFraction();
       if ( skipEM_ && emEnergyFraction > skipEMfractionThreshold_ ) continue;
 
       static PFJetMETcorrInputProducer_namespace::RawJetExtractorT<T> rawJetExtractor;
       reco::Candidate::LorentzVector rawJetP4 = rawJetExtractor(rawJet);
       if ( skipMuons_ ) {
     const std::vector<reco::CandidatePtr> & cands = rawJet.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 != 0 ? ( pfcand->muonRef().isNonnull() ? pfcand->muonRef().get() : 0) : cand->get());
       if ( mu != 0 && (*skipMuonSelection_)(*mu) ) {
         reco::Candidate::LorentzVector muonP4 = (*cand)->p4();
         rawJetP4 -= muonP4;
       }
     }
       }
 
       reco::Candidate::LorentzVector corrJetP4;
       if ( checkInputType.isPatJet(rawJet) )
         corrJetP4 = jetCorrExtractor_(rawJet, jetCorrLabel_.label(), jetCorrEtaMax_, &rawJetP4);
       else
         corrJetP4 = jetCorrExtractor_(rawJet, 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(rawJet) )
             rawJetP4offsetCorr = jetCorrExtractor_(rawJet, offsetCorrLabel_.label(), jetCorrEtaMax_, &rawJetP4);
           else
         rawJetP4offsetCorr = jetCorrExtractor_(rawJet, 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(type1Correction, "type1");
     for ( typename std::vector<type2BinningEntryType*>::const_iterator type2BinningEntry = type2Binning_.begin();
       type2BinningEntry != type2Binning_.end(); ++type2BinningEntry ) {
       evt.put(std::auto_ptr<CorrMETData>(new CorrMETData((*type2BinningEntry)->binUnclEnergySum_)), (*type2BinningEntry)->getInstanceLabel_full("type2"));
       evt.put(std::auto_ptr<CorrMETData>(new CorrMETData((*type2BinningEntry)->binOffsetEnergySum_)), (*type2BinningEntry)->getInstanceLabel_full("offset"));
     }
   }