ShiftedParticleProducer Class Reference

#include <ShiftedParticleProducer.h>

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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

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<>
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 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 33 of file ShiftedParticleProducer.h.

Member Typedef Documentation

CandidateView

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

private

Definition at line 34 of file ShiftedParticleProducer.h.

Constructor & Destructor Documentation

ShiftedParticleProducer()

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

explicit

Definition at line 4 of file ShiftedParticleProducer.cc.

                                                                           {
  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>();
}

References binning_, looper::cfg, moduleLabel_, shiftBy_, srcToken_, HLT_FULL_cff::srcWeights, AlCaHLTBitMon_QueryRunRegistry::string, and weightsToken_.

~ShiftedParticleProducer()

ShiftedParticleProducer::~ShiftedParticleProducer ( )

override

Definition at line 27 of file ShiftedParticleProducer.cc.

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

References binning_.

Member Function Documentation

getUncShift()

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

private

Definition at line 71 of file ShiftedParticleProducer.cc.

                                                                                 {
  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;
}

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

Referenced by produce().

produce()

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

overrideprivate

Definition at line 33 of file ShiftedParticleProducer.cc.

                                                                            {
  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));
}

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

Member Data Documentation

binning_

std::vector<binningEntryType*> ShiftedParticleProducer::binning_

private

Definition at line 73 of file ShiftedParticleProducer.h.

Referenced by getUncShift(), ShiftedParticleProducer(), and ~ShiftedParticleProducer().

moduleLabel_

std::string ShiftedParticleProducer::moduleLabel_

private

Definition at line 45 of file ShiftedParticleProducer.h.

Referenced by ExternalGeneratorFilter.ExternalGeneratorFilter::insertInto(), Modules.SwitchProducer::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), Modules.ESProducer::nameInProcessDesc_(), and ShiftedParticleProducer().

shiftBy_

double ShiftedParticleProducer::shiftBy_

private

Definition at line 75 of file ShiftedParticleProducer.h.

Referenced by produce(), and ShiftedParticleProducer().

srcToken_

edm::EDGetTokenT<CandidateView> ShiftedParticleProducer::srcToken_

private

Definition at line 47 of file ShiftedParticleProducer.h.

Referenced by produce(), and ShiftedParticleProducer().

weightsToken_

edm::EDGetTokenT<edm::ValueMap<float> > ShiftedParticleProducer::weightsToken_

private

Definition at line 48 of file ShiftedParticleProducer.h.

Referenced by produce(), and ShiftedParticleProducer().