#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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 20 of file MuonMETcorrInputProducer.cc.

References edm::Event::getByToken(), eostools::move(), HLT_FULL_cff::muon, muonCorrectionMapToken_, patZpeak::muons, reco::MuonMETCorrectionData::NotUsed, patZpeak::numMuons, 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));
 }