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

bool	hasAbilityToProduceInRuns () 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 36 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 45 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 53 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()));
     }
   }
 }