Computes the Type-1 corrections for pfMET. A specific version of the Type1MET class from the JetMETCorrections/Type1MET package. More...

#include <Type1PFMET.h>

Inheritance diagram for Type1PFMET:

Public Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

	Type1PFMET (const edm::ParameterSet &)

	Type1PFMET ()

	~Type1PFMET () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	run (const reco::METCollection &uncorMET, const reco::JetCorrector &corrector, const reco::PFJetCollection &uncorJet, double jetPTthreshold, double jetEMfracLimit, double jetMufracLimit, reco::METCollection *corMET)

Private Attributes
edm::EDGetTokenT< reco::JetCorrector >	correctorToken

double	jetEMfracLimit

double	jetMufracLimit

double	jetPTthreshold

edm::EDGetTokenT< reco::PFJetCollection >	tokenUncorJets

edm::EDGetTokenT< reco::METCollection >	tokenUncorMet

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

Computes the Type-1 corrections for pfMET. A specific version of the Type1MET class from the JetMETCorrections/Type1MET package.

Author: Michal Bluj

Date: February 2009

Definition at line 31 of file Type1PFMET.h.

Constructor & Destructor Documentation

Type1PFMET::Type1PFMET ( const edm::ParameterSet & iConfig )

explicit

Definition at line 14 of file Type1PFMET.cc.

References edm::ParameterSet::getParameter().

 {
   tokenUncorMet  = consumes<METCollection>(iConfig.getParameter<edm::InputTag>("inputUncorMetLabel"));
   tokenUncorJets = consumes<PFJetCollection>(iConfig.getParameter<edm::InputTag>("inputUncorJetsTag"));
   correctorToken = consumes<JetCorrector>(iConfig.getParameter<edm::InputTag>("corrector"));
   jetPTthreshold      = iConfig.getParameter<double>("jetPTthreshold");
   jetEMfracLimit      = iConfig.getParameter<double>("jetEMfracLimit");
   jetMufracLimit      = iConfig.getParameter<double>("jetMufracLimit");
   produces<METCollection>();
 }

Type1PFMET::Type1PFMET ( )

explicit

Definition at line 24 of file Type1PFMET.cc.

24 {}

Type1PFMET::~Type1PFMET ( )

override

Definition at line 27 of file Type1PFMET.cc.

27 {}

Member Function Documentation

void Type1PFMET::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 30 of file Type1PFMET.cc.

References mitigatedMETSequence_cff::corrector, edm::Event::getByToken(), eostools::move(), convertSQLitetoXML_cfg::output, edm::Handle< T >::product(), edm::Event::put(), and writedatasetfile::run.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

 {
   using namespace edm;
   Handle<PFJetCollection> inputUncorJets;
   iEvent.getByToken( tokenUncorJets, inputUncorJets );
   Handle<JetCorrector> corrector;
   iEvent.getByToken( correctorToken, corrector );
   Handle<METCollection> inputUncorMet;                     //Define Inputs
   iEvent.getByToken( tokenUncorMet,  inputUncorMet );     //Get Inputs
   std::unique_ptr<METCollection> output( new METCollection() );  //Create empty output
   run( *(inputUncorMet.product()), *(corrector.product()), *(inputUncorJets.product()),
        jetPTthreshold, jetEMfracLimit, jetMufracLimit,
        &*output );                                         //Invoke the algorithm
   iEvent.put(std::move(output));                                        //Put output into Event
 }

void Type1PFMET::run	(	const reco::METCollection &	uncorMET,
		const reco::JetCorrector &	corrector,
		const reco::PFJetCollection &	uncorJet,
		double	jetPTthreshold,
		double	jetEMfracLimit,
		double	jetMufracLimit,
		reco::METCollection *	corMET
	)

private

Definition at line 46 of file Type1PFMET.cc.

References MessageLogger_cfi::cerr, corr, reco::JetCorrector::correction(), delta, metsig::jet, nanoDQM_cfi::MET, reco::MET::mEtCorr(), CorrMETData::mex, CorrMETData::mey, reco::LeafCandidate::px(), reco::LeafCandidate::py(), mps_fire::result, mathSSE::sqrt(), CorrMETData::sumet, reco::MET::sumEt(), and reco::LeafCandidate::vertex().

Referenced by Types.EventID::cppID(), Types.LuminosityBlockID::cppID(), and o2olib.O2OTool::execute().

 {
   if (!corMET) {
     std::cerr << "Type1METAlgo_run-> undefined output MET collection. Stop. " << std::endl;
     return;
   }
 
   double DeltaPx = 0.0;
   double DeltaPy = 0.0;
   double DeltaSumET = 0.0;
   // ---------------- Calculate jet corrections, but only for those uncorrected jets
   // ---------------- which are above the given threshold.  This requires that the
   // ---------------- uncorrected jets be matched with the corrected jets.
   for( PFJetCollection::const_iterator jet = uncorJet.begin(); jet != uncorJet.end(); ++jet) {
     if( jet->pt() > jetPTthreshold ) {
       double emEFrac =
     jet->chargedEmEnergyFraction() + jet->neutralEmEnergyFraction();
       double muEFrac = jet->chargedMuEnergyFraction();
       if( emEFrac < jetEMfracLimit
       && muEFrac < jetMufracLimit ) {
     double corr = corrector.correction (*jet) - 1.; // correction itself
     DeltaPx +=  jet->px() * corr;
     DeltaPy +=  jet->py() * corr;
     DeltaSumET += jet->et() * corr;
       }
     }
   }
   //----------------- Calculate and set deltas for new MET correction
   CorrMETData delta;
   delta.mex   =  - DeltaPx;    //correction to MET (from Jets) is negative,
   delta.mey   =  - DeltaPy;    //since MET points in direction opposite of jets
   delta.sumet =  DeltaSumET;
   //----------------- Fill holder with corrected MET (= uncorrected + delta) values
   const MET* u = &(uncorMET.front());
   double corrMetPx = u->px()+delta.mex;
   double corrMetPy = u->py()+delta.mey;
   MET::LorentzVector correctedMET4vector( corrMetPx, corrMetPy, 0.,
                       sqrt (corrMetPx*corrMetPx + corrMetPy*corrMetPy)
                       );
   //----------------- get previous corrections and push into new corrections
   std::vector<CorrMETData> corrections = u->mEtCorr();
   corrections.push_back( delta );
   //----------------- Push onto MET Collection
   MET result = MET(u->sumEt()+delta.sumet,
            corrections,
            correctedMET4vector,
            u->vertex() );
   corMET->push_back(result);
 
   return;
 }

Member Data Documentation

edm::EDGetTokenT<reco::JetCorrector> Type1PFMET::correctorToken

private

Definition at line 41 of file Type1PFMET.h.

double Type1PFMET::jetEMfracLimit

private

Definition at line 43 of file Type1PFMET.h.

double Type1PFMET::jetMufracLimit

private

Definition at line 44 of file Type1PFMET.h.

double Type1PFMET::jetPTthreshold

private

Definition at line 42 of file Type1PFMET.h.

edm::EDGetTokenT<reco::PFJetCollection> Type1PFMET::tokenUncorJets

private

Definition at line 40 of file Type1PFMET.h.

edm::EDGetTokenT<reco::METCollection> Type1PFMET::tokenUncorMet

private

Definition at line 39 of file Type1PFMET.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation