#include <ShiftedMETcorrInputProducer.h>

Inheritance diagram for ShiftedMETcorrInputProducer:

Classes
struct	binningEntryType

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

	~ShiftedMETcorrInputProducer () 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 Types
typedef std::vector< edm::InputTag >	vInputTag

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

Private Attributes
std::vector< binningEntryType * >	binning_

double	shiftBy_

vInputTag	src_

std::vector< edm::EDGetTokenT< CorrMETData > >	srcTokens_

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

Vary px, py and sumEt of "unclustered energy" (PFJets of Pt < 10 GeV plus PFCandidates not within jets) by +/- 1 standard deviation, in order to estimate resulting uncertainty on MET

Author: Christian Veelken, LLR

Definition at line 25 of file ShiftedMETcorrInputProducer.h.

Member Typedef Documentation

typedef std::vector<edm::InputTag> ShiftedMETcorrInputProducer::vInputTag

private

Definition at line 33 of file ShiftedMETcorrInputProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 5 of file ShiftedMETcorrInputProducer.cc.

References binning_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), shiftBy_, src_, and srcTokens_.

                                                                                    {
   src_ = cfg.getParameter<vInputTag>("src");
 
   //--- check that all InputTags refer to the same module label
   //   (i.e. differ by instance label only)
   for (vInputTag::const_iterator src_ref = src_.begin(); src_ref != src_.end(); ++src_ref) {
     for (vInputTag::const_iterator src_test = src_ref; src_test != src_.end(); ++src_test) {
       if (src_test->label() != src_ref->label())
         throw cms::Exception("ShiftedMETcorrInputProducer")
             << "InputTags specified by 'src' Configuration parameter must not refer to different module labels !!\n";
     }
   }
 
   shiftBy_ = cfg.getParameter<double>("shiftBy");
 
   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_.push_back(new binningEntryType(*cfgBinningEntry));
     }
   } else {
     double uncertainty = cfg.getParameter<double>("uncertainty");
     binning_.push_back(new binningEntryType(uncertainty));
   }
 
   for (vInputTag::const_iterator src_i = src_.begin(); src_i != src_.end(); ++src_i) {
     for (std::vector<binningEntryType*>::const_iterator binningEntry = binning_.begin(); binningEntry != binning_.end();
          ++binningEntry) {
       srcTokens_.push_back(consumes<CorrMETData>(
           edm::InputTag(src_i->label(), (*binningEntry)->getInstanceLabel_full(src_i->instance()))));
       produces<CorrMETData>((*binningEntry)->getInstanceLabel_full(src_i->instance()));
     }
   }
 }

ShiftedMETcorrInputProducer::~ShiftedMETcorrInputProducer ( )

override

Definition at line 42 of file ShiftedMETcorrInputProducer.cc.

References binning_.

                                                           {
   for (std::vector<binningEntryType*>::const_iterator it = binning_.begin(); it != binning_.end(); ++it) {
     delete (*it);
   }
 }

Member Function Documentation

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

overrideprivate

Definition at line 48 of file ShiftedMETcorrInputProducer.cc.

References binning_, DEFINE_FWK_MODULE, edm::Event::getByToken(), eostools::move(), edm::Event::put(), edm::shift, shiftBy_, src_, and srcTokens_.

                                                                                 {
   unsigned countToken(0);
   for (vInputTag::const_iterator src_i = src_.begin(); src_i != src_.end(); ++src_i) {
     for (std::vector<binningEntryType*>::iterator binningEntry = binning_.begin(); binningEntry != binning_.end();
          ++binningEntry) {
       edm::Handle<CorrMETData> originalObject;
       evt.getByToken(srcTokens_.at(countToken), originalObject);
       ++countToken;
 
       double shift = shiftBy_ * (*binningEntry)->binUncertainty_;
 
       auto shiftedObject = std::make_unique<CorrMETData>(*originalObject);
       //--- MET balances momentum of reconstructed particles,
       //    hence variations of "unclustered energy" and MET are opposite in sign
       shiftedObject->mex = -shift * originalObject->mex;
       shiftedObject->mey = -shift * originalObject->mey;
       shiftedObject->sumet = shift * originalObject->sumet;
 
       evt.put(std::move(shiftedObject), (*binningEntry)->getInstanceLabel_full(src_i->instance()));
     }
   }
 }