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	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () 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_2018_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_2018_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().