#include <MuonMETcorrInputProducer.h>

Inheritance diagram for MuonMETcorrInputProducer:

Public Member Functions
	MuonMETcorrInputProducer (const edm::ParameterSet &)

	~MuonMETcorrInputProducer () override

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

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

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

Private Attributes
std::string	moduleLabel_

edm::EDGetTokenT< edm::ValueMap< reco::MuonMETCorrectionData > >	muonCorrectionMapToken_

edm::EDGetTokenT< reco::MuonCollection >	token_

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

Sum CaloMET muon corrections, needed to compute Type 2 MET corrections in case muon corrected CaloMET is used as input to CorrectedCaloMETProducer

Author: Christian Veelken, LLR

Definition at line 27 of file MuonMETcorrInputProducer.h.

Constructor & Destructor Documentation

MuonMETcorrInputProducer::MuonMETcorrInputProducer ( const edm::ParameterSet & cfg )

explicit

Definition at line 7 of file MuonMETcorrInputProducer.cc.

References edm::ParameterSet::getParameter(), muonCorrectionMapToken_, and token_.

   : moduleLabel_(cfg.getParameter<std::string>("@module_label"))
 {
   token_ = consumes<reco::MuonCollection>(cfg.getParameter<edm::InputTag>("src"));
   muonCorrectionMapToken_ = consumes<edm::ValueMap<reco::MuonMETCorrectionData> >(cfg.getParameter<edm::InputTag>("srcMuonCorrections"));
 
   produces<CorrMETData>();
 }

MuonMETcorrInputProducer::~MuonMETcorrInputProducer ( )

override

Definition at line 16 of file MuonMETcorrInputProducer.cc.

 {
 // nothing to be done yet...
 }

Member Function Documentation

void MuonMETcorrInputProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

overrideprivate

Definition at line 21 of file MuonMETcorrInputProducer.cc.

References DEFINE_FWK_MODULE, edm::Event::getByToken(), eostools::move(), muonCorrectionMapToken_, extraflags_cff::muons, reco::MuonMETCorrectionData::NotUsed, reco::LeafCandidate::pt(), edm::Event::put(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), token_, and reco::MuonMETCorrectionData::type().

 {
   std::unique_ptr<CorrMETData> unclEnergySum(new CorrMETData());
 
   edm::Handle<reco::MuonCollection> muons;
   evt.getByToken(token_, muons);
 
   typedef edm::ValueMap<reco::MuonMETCorrectionData> MuonMETCorrectionMap;
   edm::Handle<MuonMETCorrectionMap> muonCorrections;
   evt.getByToken(muonCorrectionMapToken_, muonCorrections);
 
 //--- sum muon corrections. 
 //
 //    NOTE: MET = -(jets + muon corrections + "unclustered energy"),
 //          so "unclustered energy" = -(MET + jets + muons),
 //          i.e. muon corrections enter the sum of "unclustered energy" with negative sign.
 //
   int numMuons = muons->size();
   for ( int muonIndex = 0; muonIndex < numMuons; ++muonIndex ) {
     const reco::Muon& muon = muons->at(muonIndex);
     
     reco::MuonRef muonRef(muons, muonIndex);
 
     reco::MuonMETCorrectionData muonCorrection = (*muonCorrections)[muonRef];
     if ( muonCorrection.type() != reco::MuonMETCorrectionData::NotUsed ) {
       unclEnergySum->mex   -= muon.px();
       unclEnergySum->mey   -= muon.py();
       unclEnergySum->sumet -= muon.pt();
     }
   }
 
 //--- add sum of muon corrections to the event
   evt.put(std::move(unclEnergySum));
 }