#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

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::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<>
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 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 5 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_2018_cff::srcWeights, token_, and weightsToken_.

◆ ~PFCandMETcorrInputProducer()

PFCandMETcorrInputProducer::~PFCandMETcorrInputProducer ( )

override

Definition at line 28 of file PFCandMETcorrInputProducer.cc.

28 {}

Member Function Documentation

◆ produce()

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

overrideprivate

Definition at line 30 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::unique_ptr<CorrMETData>(new CorrMETData((*binningEntry)->binUnclEnergySum_)),
             (*binningEntry)->binLabel_);
   }
 }