#include <ShiftedParticleProducerT.h>

Inheritance diagram for ShiftedParticleProducerT< T >:

Classes
struct	binningEntryType

Public Member Functions
	ShiftedParticleProducerT (const edm::ParameterSet &cfg)

	~ShiftedParticleProducerT () 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 Types
typedef std::vector< T >	ParticleCollection

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

Private Attributes
std::vector< binningEntryType * >	binning_

std::string	moduleLabel_

double	shiftBy_

edm::EDGetTokenT < ParticleCollection >	srcToken_

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

template<typename T>
class ShiftedParticleProducerT< T >

Vary energy of electrons/muons/tau-jets by +/- 1 standard deviation, in order to estimate resulting uncertainty on MET

NOTE: energy scale uncertainties need to be specified in python config

Author: Christian Veelken, LLR

Definition at line 31 of file ShiftedParticleProducerT.h.

Member Typedef Documentation

template<typename T >

typedef std::vector<T> ShiftedParticleProducerT< T >::ParticleCollection

private

Definition at line 32 of file ShiftedParticleProducerT.h.

Constructor & Destructor Documentation

template<typename T >

ShiftedParticleProducerT< T >::ShiftedParticleProducerT ( const edm::ParameterSet & cfg )

inlineexplicit

Definition at line 35 of file ShiftedParticleProducerT.h.

References ShiftedParticleProducerT< T >::binning_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), ShiftedParticleProducerT< T >::shiftBy_, and ShiftedParticleProducerT< T >::srcToken_.

       : moduleLabel_(cfg.getParameter<std::string>("@module_label")) {
     srcToken_ = consumes<ParticleCollection>(cfg.getParameter<edm::InputTag>("src"));
 
     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));
     }
 
     produces<ParticleCollection>();
   }

template<typename T >

ShiftedParticleProducerT< T >::~ShiftedParticleProducerT ( )

inlineoverride

Definition at line 55 of file ShiftedParticleProducerT.h.

References ShiftedParticleProducerT< T >::binning_.

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

Member Function Documentation

template<typename T >

void ShiftedParticleProducerT< T >::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

inlineoverrideprivate

Definition at line 62 of file ShiftedParticleProducerT.h.

References ShiftedParticleProducerT< T >::binning_, edm::Event::getByToken(), edm::isNotFinite(), eostools::move(), edm::Event::put(), edm::shift, ShiftedParticleProducerT< T >::shiftBy_, and ShiftedParticleProducerT< T >::srcToken_.

                                                                 {
     edm::Handle<ParticleCollection> originalParticles;
     evt.getByToken(srcToken_, originalParticles);
 
     auto shiftedParticles = std::make_unique<ParticleCollection>();
 
     for (typename ParticleCollection::const_iterator originalParticle = originalParticles->begin();
          originalParticle != originalParticles->end();
          ++originalParticle) {
       double uncertainty = 0.;
       for (typename std::vector<binningEntryType*>::iterator binningEntry = binning_.begin();
            binningEntry != binning_.end();
            ++binningEntry) {
         if ((!(*binningEntry)->binSelection_) || (*(*binningEntry)->binSelection_)(*originalParticle)) {
           uncertainty = (*binningEntry)->binUncertainty_;
           break;
         }
       }
 
       double shift = shiftBy_ * uncertainty;
 
       reco::Candidate::LorentzVector shiftedParticleP4 = originalParticle->p4();
       //leave 0*nan = 0
       if (!(edm::isNotFinite(shift) && shiftedParticleP4.mag2() == 0))
         shiftedParticleP4 *= (1. + shift);
 
       T shiftedParticle(*originalParticle);
       shiftedParticle.setP4(shiftedParticleP4);
 
       shiftedParticles->push_back(shiftedParticle);
     }
 
     evt.put(std::move(shiftedParticles));
   }