#include <PFCandMETcorrInputProducer.h>

Inheritance diagram for PFCandMETcorrInputProducer:

Classes
struct	binningEntryType

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

	~PFCandMETcorrInputProducer () 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::vector< std::unique_ptr< binningEntryType > >	binning_

std::string	moduleLabel_

edm::EDGetTokenT< edm::View< reco::Candidate > >	token_

edm::EDGetTokenT< edm::ValueMap< float > >	weightsToken_

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 PFCandidates not within jets ("unclustered energy"), needed as input for Type 2 MET corrections

Authors: Michael Schmitt, Richard Cavanaugh, The University of Florida Florent Lacroix, University of Illinois at Chicago Christian Veelken, LLR

Definition at line 31 of file PFCandMETcorrInputProducer.h.

Constructor & Destructor Documentation

◆ PFCandMETcorrInputProducer()

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

explicit

Definition at line 7 of file PFCandMETcorrInputProducer.cc.

     : moduleLabel_(cfg.getParameter<std::string>("@module_label")) {
   token_ = consumes<edm::View<reco::Candidate>>(cfg.getParameter<edm::InputTag>("src"));
   edm::InputTag srcWeights = cfg.getParameter<edm::InputTag>("srcWeights");
   if (!srcWeights.label().empty())
     weightsToken_ = consumes<edm::ValueMap<float>>(srcWeights);
  
   if (cfg.exists("binning")) {
     typedef std::vector<edm::ParameterSet> vParameterSet;
     vParameterSet cfgBinning = cfg.getParameter<vParameterSet>("binning");
     for (vParameterSet::const_iterator cfgBinningEntry = cfgBinning.begin(); cfgBinningEntry != cfgBinning.end();
          ++cfgBinningEntry) {
       binning_.emplace_back(new binningEntryType(*cfgBinningEntry));
     }
   } else {
     binning_.emplace_back(new binningEntryType());
   }
  
   for (auto binningEntry = binning_.begin(); binningEntry != binning_.end(); ++binningEntry) {
     produces<CorrMETData>((*binningEntry)->binLabel_);
   }
 }

References binning_, looper::cfg, HLT_FULL_cff::srcWeights, token_, and weightsToken_.

◆ ~PFCandMETcorrInputProducer()

PFCandMETcorrInputProducer::~PFCandMETcorrInputProducer ( )

override

Definition at line 30 of file PFCandMETcorrInputProducer.cc.

30 {}

Member Function Documentation

◆ produce()

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

overrideprivate

Definition at line 32 of file PFCandMETcorrInputProducer.cc.

                                                                                {
   for (auto binningEntry = binning_.begin(); binningEntry != binning_.end(); ++binningEntry) {
     (*binningEntry)->binUnclEnergySum_ = CorrMETData();
   }
  
   typedef edm::View<reco::Candidate> CandidateView;
   edm::Handle<CandidateView> cands;
   evt.getByToken(token_, cands);
  
   edm::Handle<edm::ValueMap<float>> weights;
   if (!weightsToken_.isUninitialized())
     evt.getByToken(weightsToken_, weights);
  
   for (edm::View<reco::Candidate>::const_iterator cand = cands->begin(); cand != cands->end(); ++cand) {
     float weight = (!weightsToken_.isUninitialized()) ? (*weights)[cands->ptrAt(cand - cands->begin())] : 1.0;
     for (auto binningEntry = binning_.begin(); binningEntry != binning_.end(); ++binningEntry) {
       if (!(*binningEntry)->binSelection_ || (*(*binningEntry)->binSelection_)(cand->p4() * weight)) {
         (*binningEntry)->binUnclEnergySum_.mex += cand->px() * weight;
         (*binningEntry)->binUnclEnergySum_.mey += cand->py() * weight;
         (*binningEntry)->binUnclEnergySum_.sumet += cand->et() * weight;
       }
     }
   }
  
   //--- add momentum sum of PFCandidates not within jets ("unclustered energy") to the event
   for (auto binningEntry = binning_.cbegin(); binningEntry != binning_.cend(); ++binningEntry) {
     evt.put(std::make_unique<CorrMETData>((*binningEntry)->binUnclEnergySum_), (*binningEntry)->binLabel_);
   }
 }