Inheritance diagram for ShiftedParticleProducer:

Classes
struct	binningEntryType

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

	~ShiftedParticleProducer () 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 edm::View < reco::Candidate >	CandidateView

Private Member Functions
double	getUncShift (const reco::Candidate &originalParticle)

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

Private Attributes
std::vector< binningEntryType * >	binning_

std::string	moduleLabel_

double	shiftBy_

edm::EDGetTokenT< CandidateView >	srcToken_

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

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: Matthieu marionneau ETH

Definition at line 30 of file ShiftedParticleProducer.cc.

Member Typedef Documentation

typedef edm::View<reco::Candidate> ShiftedParticleProducer::CandidateView

private

Definition at line 31 of file ShiftedParticleProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 75 of file ShiftedParticleProducer.cc.

References binning_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), moduleLabel_, shiftBy_, srcToken_, HLT_FULL_cff::srcWeights, AlCaHLTBitMon_QueryRunRegistry::string, and weightsToken_.

                                                                            {
   moduleLabel_ = cfg.getParameter<std::string>("@module_label");
   srcToken_ = consumes<CandidateView>(cfg.getParameter<edm::InputTag>("src"));
   shiftBy_ = cfg.getParameter<double>("shiftBy");
   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_.push_back(new binningEntryType(*cfgBinningEntry, moduleLabel_));
     }
   } else {
     std::string uncertainty = cfg.getParameter<std::string>("uncertainty");
     binning_.push_back(new binningEntryType(uncertainty, moduleLabel_));
   }
 
   produces<reco::CandidateCollection>();
 }

ShiftedParticleProducer::~ShiftedParticleProducer ( )

override

Definition at line 98 of file ShiftedParticleProducer.cc.

References binning_.

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

Member Function Documentation

double ShiftedParticleProducer::getUncShift ( const reco::Candidate & originalParticle )

private

Definition at line 142 of file ShiftedParticleProducer.cc.

References binning_, reco::Candidate::energy(), reco::Candidate::eta(), and reco::Candidate::pt().

Referenced by produce().

                                                                                  {
   double valx = 0;
   double valy = 0;
   for (std::vector<binningEntryType*>::iterator binningEntry = binning_.begin(); binningEntry != binning_.end();
        ++binningEntry) {
     if ((!(*binningEntry)->binSelection_) || (*(*binningEntry)->binSelection_)(originalParticle)) {
       if ((*binningEntry)->energyDep_)
         valx = originalParticle.energy();
       else
         valx = originalParticle.pt();
 
       valy = originalParticle.eta();
       return (*binningEntry)->binUncFormula_->Eval(valx, valy);
     }
   }
   return 0;
 }

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

overrideprivate

Definition at line 104 of file ShiftedParticleProducer.cc.

References edm::Event::getByToken(), getUncShift(), mps_fire::i, edm::Ptr< T >::id(), edm::isNotFinite(), edm::EDGetTokenT< T >::isUninitialized(), eostools::move(), reco::Candidate::p4(), edm::Event::put(), reco::LeafCandidate::setP4(), edm::shift, shiftBy_, srcToken_, histoStyle::weight, HLT_FULL_cff::weights, and weightsToken_.

                                                                             {
   edm::Handle<CandidateView> originalParticles;
   evt.getByToken(srcToken_, originalParticles);
 
   edm::Handle<edm::ValueMap<float>> weights;
   if (!weightsToken_.isUninitialized())
     evt.getByToken(weightsToken_, weights);
 
   auto shiftedParticles = std::make_unique<reco::CandidateCollection>();
 
   for (unsigned i = 0; i < originalParticles->size(); ++i) {
     float weight = 1.0;
     if (!weightsToken_.isUninitialized()) {
       edm::Ptr<reco::Candidate> particlePtr = originalParticles->ptrAt(i);
       while (!weights->contains(particlePtr.id()) && (particlePtr->numberOfSourceCandidatePtrs() > 0))
         particlePtr = particlePtr->sourceCandidatePtr(0);
       weight = (*weights)[particlePtr];
     }
     const reco::Candidate& originalParticle = originalParticles->at(i);
     reco::LeafCandidate weightedParticle(originalParticle);
     weightedParticle.setP4(originalParticle.p4() * weight);
     double uncertainty = getUncShift(weightedParticle);
     double shift = shiftBy_ * uncertainty;
 
     reco::Candidate::LorentzVector shiftedParticleP4 = originalParticle.p4();
     //leave 0*nan = 0
     if ((weight > 0) && (!(edm::isNotFinite(shift) && shiftedParticleP4.mag2() == 0)))
       shiftedParticleP4 *= (1. + shift);
 
     std::unique_ptr<reco::Candidate> shiftedParticle = std::make_unique<reco::LeafCandidate>(originalParticle);
     shiftedParticle->setP4(shiftedParticleP4);
 
     shiftedParticles->push_back(shiftedParticle.release());
   }
 
   evt.put(std::move(shiftedParticles));
 }