MultShiftMETcorrInputProducer Class Reference

#include <MultShiftMETcorrInputProducer.h>

Inheritance diagram for MultShiftMETcorrInputProducer:

Public Member Functions
	MultShiftMETcorrInputProducer (const edm::ParameterSet &)
	~MultShiftMETcorrInputProducer () 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 Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Static Private Member Functions
static int	translateTypeToAbsPdgId (reco::PFCandidate::ParticleType type)

Private Attributes
std::vector< edm::ParameterSet >	cfgCorrParameters_
std::vector< int >	counts_
std::vector< double >	etaMax_
std::vector< double >	etaMin_
std::vector< std::unique_ptr< TF1 > >	formula_x_
std::vector< std::unique_ptr< TF1 > >	formula_y_
std::string	moduleLabel_
edm::EDGetTokenT< edm::View< reco::Candidate > >	pflow_
std::vector< double >	sumPt_
std::vector< int >	type_
std::vector< int >	varType_
edm::EDGetTokenT< edm::View< reco::Vertex > >	vertices_
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

Compute MET correction to compensate systematic shift of MET in x/y-direction (cf. https://indico.cern.ch/getFile.py/access?contribId=1&resId=0&materialId=slides&confId=174318 )

Authors

Robert Schoefbeck, Vienna

Definition at line 32 of file MultShiftMETcorrInputProducer.h.

Constructor & Destructor Documentation

MultShiftMETcorrInputProducer()

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

explicit

Definition at line 34 of file MultShiftMETcorrInputProducer.cc.

    : moduleLabel_(cfg.getParameter<std::string>("@module_label")) {
  pflow_ = consumes<edm::View<reco::Candidate>>(cfg.getParameter<edm::InputTag>("srcPFlow"));
  vertices_ = consumes<edm::View<reco::Vertex>>(cfg.getParameter<edm::InputTag>("vertexCollection"));
 
  edm::InputTag srcWeights = cfg.getParameter<edm::InputTag>("srcWeights");
  if (!srcWeights.label().empty())
    weightsToken_ = consumes<edm::ValueMap<float>>(srcWeights);
 
  cfgCorrParameters_ = cfg.getParameter<std::vector<edm::ParameterSet>>("parameters");
  etaMin_.clear();
  etaMax_.clear();
  type_.clear();
  varType_.clear();
 
  produces<CorrMETData>();
 
  for (std::vector<edm::ParameterSet>::const_iterator v = cfgCorrParameters_.begin(); v != cfgCorrParameters_.end();
       v++) {
    TString corrPxFormula = v->getParameter<std::string>("fx");
    TString corrPyFormula = v->getParameter<std::string>("fy");
    std::vector<double> corrPxParams = v->getParameter<std::vector<double>>("px");
    std::vector<double> corrPyParams = v->getParameter<std::vector<double>>("py");
 
    formula_x_.push_back(std::unique_ptr<TF1>(new TF1(
        std::string(moduleLabel_).append("_").append(v->getParameter<std::string>("name")).append("_corrPx").c_str(),
        v->getParameter<std::string>("fx").c_str())));
    formula_y_.push_back(std::unique_ptr<TF1>(new TF1(
        std::string(moduleLabel_).append("_").append(v->getParameter<std::string>("name")).append("_corrPy").c_str(),
        v->getParameter<std::string>("fy").c_str())));
 
    for (unsigned i = 0; i < corrPxParams.size(); i++)
      formula_x_.back()->SetParameter(i, corrPxParams[i]);
    for (unsigned i = 0; i < corrPyParams.size(); i++)
      formula_y_.back()->SetParameter(i, corrPyParams[i]);
 
    counts_.push_back(0);
    sumPt_.push_back(0.);
    etaMin_.push_back(v->getParameter<double>("etaMin"));
    etaMax_.push_back(v->getParameter<double>("etaMax"));
    type_.push_back(v->getParameter<int>("type"));
    varType_.push_back(v->getParameter<int>("varType"));
  }
}

References mps_setup::append, looper::cfg, cfgCorrParameters_, counts_, etaMax_, etaMin_, formula_x_, formula_y_, mps_fire::i, moduleLabel_, pflow_, HLT_2018_cff::srcWeights, AlCaHLTBitMon_QueryRunRegistry::string, sumPt_, tools::TF1, type_, findQualityFiles::v, varType_, vertices_, and weightsToken_.

~MultShiftMETcorrInputProducer()

MultShiftMETcorrInputProducer::~MultShiftMETcorrInputProducer ( )

override

Definition at line 79 of file MultShiftMETcorrInputProducer.cc.

{}

Member Function Documentation

produce()

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

overrideprivate

Definition at line 81 of file MultShiftMETcorrInputProducer.cc.

                                                                                  {
  //get primary vertices
  edm::Handle<edm::View<reco::Vertex>> hpv;
  evt.getByToken(vertices_, hpv);
  if (!hpv.isValid()) {
    edm::LogError("MultShiftMETcorrInputProducer::produce") << "could not find vertex collection ";
  }
  std::vector<reco::Vertex> goodVertices;
  for (unsigned i = 0; i < hpv->size(); i++) {
    if ((*hpv)[i].ndof() > 4 && (fabs((*hpv)[i].z()) <= 24.) && (fabs((*hpv)[i].position().rho()) <= 2.0))
      goodVertices.push_back((*hpv)[i]);
  }
  int ngoodVertices = goodVertices.size();
 
  for (unsigned i = 0; i < counts_.size(); i++)
    counts_[i] = 0;
  for (unsigned i = 0; i < sumPt_.size(); i++)
    sumPt_[i] = 0.;
 
  edm::Handle<edm::View<reco::Candidate>> particleFlow;
  evt.getByToken(pflow_, particleFlow);
 
  edm::Handle<edm::ValueMap<float>> weights;
  if (!weightsToken_.isUninitialized())
    evt.getByToken(weightsToken_, weights);
  for (unsigned i = 0; i < particleFlow->size(); ++i) {
    const reco::Candidate& c = particleFlow->at(i);
    for (unsigned j = 0; j < type_.size(); j++) {
      if (abs(c.pdgId()) == translateTypeToAbsPdgId(reco::PFCandidate::ParticleType(type_[j]))) {
        if ((c.eta() > etaMin_[j]) and (c.eta() < etaMax_[j])) {
          float weight = (!weightsToken_.isUninitialized()) ? (*weights)[particleFlow->ptrAt(i)] : 1.0;
          counts_[j] += (weight > 0);
          sumPt_[j] += c.pt() * weight;
          continue;
        }
      }
    }
  }
 
  //MM: loop over all constituent types and sum each correction
  std::unique_ptr<CorrMETData> metCorr(new CorrMETData());
 
  double corx = 0.;
  double cory = 0.;
 
  for (std::vector<edm::ParameterSet>::const_iterator v = cfgCorrParameters_.begin(); v != cfgCorrParameters_.end();
       v++) {
    unsigned j = v - cfgCorrParameters_.begin();
 
    double val(0.);
    if (varType_[j] == 0) {
      val = counts_[j];
    }
    if (varType_[j] == 1) {
      val = ngoodVertices;
    }
    if (varType_[j] == 2) {
      val = sumPt_[j];
    }
 
    corx -= formula_x_[j]->Eval(val);
    cory -= formula_y_[j]->Eval(val);
 
  }  //end loop over corrections
 
  metCorr->mex = corx;
  metCorr->mey = cory;
  evt.put(std::move(metCorr), "");
}

References funct::abs(), HltBtagPostValidation_cff::c, cfgCorrParameters_, counts_, etaMax_, etaMin_, formula_x_, formula_y_, edm::Event::getByToken(), metFilters_cff::goodVertices, mps_fire::i, edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), dqmiolumiharvest::j, eostools::move(), pfLinker_cff::particleFlow, pflow_, position, edm::Event::put(), rho, sumPt_, translateTypeToAbsPdgId(), type_, findQualityFiles::v, heppy_batch::val, varType_, vertices_, mps_merge::weight, HLT_2018_cff::weights, weightsToken_, and z.

translateTypeToAbsPdgId()

int MultShiftMETcorrInputProducer::translateTypeToAbsPdgId ( reco::PFCandidate::ParticleType  type )

staticprivate

Definition at line 12 of file MultShiftMETcorrInputProducer.cc.

                                                                                           {
  switch (type) {
    case reco::PFCandidate::ParticleType::h:
      return 211;  // pi+
    case reco::PFCandidate::ParticleType::e:
      return 11;
    case reco::PFCandidate::ParticleType::mu:
      return 13;
    case reco::PFCandidate::ParticleType::gamma:
      return 22;
    case reco::PFCandidate::ParticleType::h0:
      return 130;  // K_L0
    case reco::PFCandidate::ParticleType::h_HF:
      return 1;  // dummy pdg code
    case reco::PFCandidate::ParticleType::egamma_HF:
      return 2;  // dummy pdg code
    case reco::PFCandidate::ParticleType::X:
    default:
      return 0;
  }
}

References MillePedeFileConverter_cfg::e, CustomPhysics_cfi::gamma, h, amptDefaultParameters_cff::mu, and X.

Referenced by produce().

Member Data Documentation

cfgCorrParameters_

std::vector<edm::ParameterSet> MultShiftMETcorrInputProducer::cfgCorrParameters_

private

Definition at line 46 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

counts_

std::vector<int> MultShiftMETcorrInputProducer::counts_

private

Definition at line 49 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

etaMax_

std::vector<double> MultShiftMETcorrInputProducer::etaMax_

private

Definition at line 48 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

etaMin_

std::vector<double> MultShiftMETcorrInputProducer::etaMin_

private

Definition at line 48 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

formula_x_

std::vector<std::unique_ptr<TF1> > MultShiftMETcorrInputProducer::formula_x_

private

Definition at line 51 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

formula_y_

std::vector<std::unique_ptr<TF1> > MultShiftMETcorrInputProducer::formula_y_

private

Definition at line 52 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

moduleLabel_

std::string MultShiftMETcorrInputProducer::moduleLabel_

private

Definition at line 43 of file MultShiftMETcorrInputProducer.h.

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

pflow_

edm::EDGetTokenT<edm::View<reco::Candidate> > MultShiftMETcorrInputProducer::pflow_

private

Definition at line 41 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

sumPt_

std::vector<double> MultShiftMETcorrInputProducer::sumPt_

private

Definition at line 50 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

type_

std::vector<int> MultShiftMETcorrInputProducer::type_

private

Definition at line 49 of file MultShiftMETcorrInputProducer.h.

Referenced by Modules.Service::__str__(), Modules.Service::_placeImpl(), Modules.ESSource::_placeImpl(), Modules.ESProducer::_placeImpl(), Modules.ESPrefer::copy(), Mixins._TypedParameterizable::dumpPython(), Modules.ESPrefer::dumpPythonAs(), Modules.Service::dumpSequencePython(), RandomRunSource.RandomRunSource::insertInto(), Modules.Service::insertInto(), ExternalGeneratorFilter.ExternalGeneratorFilter::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::moduleLabel_(), Modules.ESProducer::moduleLabel_(), MultShiftMETcorrInputProducer(), Modules.ESSource::nameInProcessDesc_(), Modules.ESProducer::nameInProcessDesc_(), Modules.ESPrefer::nameInProcessDesc_(), and produce().

varType_

std::vector<int> MultShiftMETcorrInputProducer::varType_

private

Definition at line 49 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

vertices_

edm::EDGetTokenT<edm::View<reco::Vertex> > MultShiftMETcorrInputProducer::vertices_

private

Definition at line 42 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().

weightsToken_

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

private

Definition at line 44 of file MultShiftMETcorrInputProducer.h.

Referenced by MultShiftMETcorrInputProducer(), and produce().