#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	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () 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 20 of file MuonMETcorrInputProducer.cc.

References DEFINE_FWK_MODULE, edm::Event::getByToken(), eostools::move(), muonCorrectionMapToken_, PDWG_BPHSkim_cff::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));
 }