MultShiftMETcorrDBInputProducer Class Reference

#include <MultShiftMETcorrDBInputProducer.h>

Inheritance diagram for MultShiftMETcorrDBInputProducer:

Public Member Functions
	MultShiftMETcorrDBInputProducer (const edm::ParameterSet &)
	~MultShiftMETcorrDBInputProducer () 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 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_
int	counts_
std::vector< double >	etaMax_
std::vector< double >	etaMin_
std::unique_ptr< TF1 >	formula_x_
std::unique_ptr< TF1 >	formula_y_
bool	mIsData
edm::ESGetToken< MEtXYcorrectParametersCollection, MEtXYcorrectRecord >	mMEtXYcorParaColl_
std::string	moduleLabel_
std::string	mPayloadName
std::string	mSampleType
edm::EDGetTokenT< edm::View< reco::Candidate > >	pflow_
double	sumPt_
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

SangEun Lee,

Robert Schoefbeck, Vienna

Definition at line 36 of file MultShiftMETcorrDBInputProducer.h.

Constructor & Destructor Documentation

MultShiftMETcorrDBInputProducer()

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

explicit

Definition at line 38 of file MultShiftMETcorrDBInputProducer.cc.

    : moduleLabel_(cfg.getParameter<std::string>("@module_label")) {
  mPayloadName = cfg.getParameter<std::string>("payloadName");
  mSampleType = (cfg.exists("sampleType")) ? cfg.getParameter<std::string>("sampleType") : "MC";
  mIsData = cfg.getParameter<bool>("isData");
 
  pflow_ = consumes<edm::View<reco::Candidate>>(cfg.getParameter<edm::InputTag>("srcPFlow"));
  vertices_ = consumes<edm::View<reco::Vertex>>(cfg.getParameter<edm::InputTag>("vertexCollection"));
  mMEtXYcorParaColl_ =
      esConsumes<MEtXYcorrectParametersCollection, MEtXYcorrectRecord>(edm::ESInputTag("", mPayloadName));
 
  edm::InputTag srcWeights = cfg.getParameter<edm::InputTag>("srcWeights");
  if (!srcWeights.label().empty())
    weightsToken_ = consumes<edm::ValueMap<float>>(srcWeights);
 
  etaMin_.clear();
  etaMax_.clear();
 
  produces<CorrMETData>();
}

References looper::cfg, etaMax_, etaMin_, mIsData, mMEtXYcorParaColl_, mPayloadName, mSampleType, pflow_, HLT_FULL_cff::srcWeights, AlCaHLTBitMon_QueryRunRegistry::string, vertices_, and weightsToken_.

~MultShiftMETcorrDBInputProducer()

MultShiftMETcorrDBInputProducer::~MultShiftMETcorrDBInputProducer ( )

override

Definition at line 59 of file MultShiftMETcorrDBInputProducer.cc.

{}

Member Function Documentation

produce()

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

overrideprivate

Definition at line 61 of file MultShiftMETcorrDBInputProducer.cc.

                                                                                    {
  // Get para.s from DB
  edm::ESHandle<MEtXYcorrectParametersCollection> MEtXYcorParaColl = es.getHandle(mMEtXYcorParaColl_);
 
  // get the sections from Collection (pair of section and METCorr.Par class)
  std::vector<MEtXYcorrectParametersCollection::key_type> keys;
  // save level to keys for each METParameter in METParameter collection
  MEtXYcorParaColl->validKeys(keys);
 
  //get primary vertices
  edm::Handle<edm::View<reco::Vertex>> hpv;
  evt.getByToken(vertices_, hpv);
  if (!hpv.isValid()) {
    edm::LogError("MultShiftMETcorrDBInputProducer::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();
 
  edm::Handle<edm::View<reco::Candidate>> particleFlow;
  evt.getByToken(pflow_, particleFlow);
 
  edm::Handle<edm::ValueMap<float>> weights;
  if (!weightsToken_.isUninitialized())
    evt.getByToken(weightsToken_, weights);
 
  //loop over all constituent types and sum each correction
  //std::unique_ptr<CorrMETData> metCorr(new CorrMETData());
  std::unique_ptr<CorrMETData> metCorr(new CorrMETData());
 
  double corx = 0.;
  double cory = 0.;
 
  // check DB
  for (std::vector<MEtXYcorrectParametersCollection::key_type>::const_iterator ikey = keys.begin(); ikey != keys.end();
       ++ikey) {
    if (mIsData) {
      if (!MEtXYcorParaColl->isShiftData(*ikey))
        throw cms::Exception("MultShiftMETcorrDBInputProducer::produce")
            << "DB is not for Data. Set proper option: \"corrPfMetXYMultDB.isData\" !!\n";
    } else {
      if (MEtXYcorParaColl->isShiftData(*ikey))
        throw cms::Exception("MultShiftMETcorrDBInputProducer::produce")
            << "DB is for Data. Set proper option: \"corrPfMetXYMultDB.isData\" !!\n";
    }
  }
 
  for (std::vector<MEtXYcorrectParametersCollection::key_type>::const_iterator ikey = keys.begin(); ikey != keys.end();
       ++ikey) {
    if (!mIsData) {
      if (mSampleType == "MC") {
        if (!MEtXYcorParaColl->isShiftMC(*ikey))
          continue;
      } else if (mSampleType == "DY") {
        if (!MEtXYcorParaColl->isShiftDY(*ikey))
          continue;
      } else if (mSampleType == "TTJets") {
        if (!MEtXYcorParaColl->isShiftTTJets(*ikey))
          continue;
      } else if (mSampleType == "WJets") {
        if (!MEtXYcorParaColl->isShiftWJets(*ikey))
          continue;
      } else
        throw cms::Exception("MultShiftMETcorrDBInputProducer::produce")
            << "SampleType: " << mSampleType << " is not reserved !!!\n";
    }
 
    std::string sectionName = MEtXYcorParaColl->findLabel(*ikey);
    MEtXYcorrectParameters const& MEtXYcorParams = (*MEtXYcorParaColl)[*ikey];
 
    counts_ = 0;
    sumPt_ = 0;
 
    for (unsigned i = 0; i < particleFlow->size(); ++i) {
      const reco::Candidate& c = particleFlow->at(i);
      if (abs(c.pdgId()) ==
          translateTypeToAbsPdgId(reco::PFCandidate::ParticleType(MEtXYcorParams.definitions().PtclType()))) {
        if ((c.eta() > MEtXYcorParams.record(0).xMin(0)) and (c.eta() < MEtXYcorParams.record(0).xMax(0))) {
          float weight = (!weightsToken_.isUninitialized()) ? (*weights)[particleFlow->ptrAt(i)] : 1.0;
          counts_ += (weight > 0);
          sumPt_ += c.pt() * weight;
          continue;
        }
      }
    }
    double val(0.);
    unsigned parVar = MEtXYcorParams.definitions().parVar(0);
 
    if (parVar == 0) {
      val = counts_;
 
    } else if (parVar == 1) {
      val = ngoodVertices;
 
    } else if (parVar == 2) {
      val = sumPt_;
 
    } else {
      throw cms::Exception("MultShiftMETcorrDBInputProducer::produce")
          << "parVar: " << parVar << " is not reserved !!!\n";
    }
 
    formula_x_ = std::make_unique<TF1>("corrPx", MEtXYcorParams.definitions().formula().c_str());
    formula_y_ = std::make_unique<TF1>("corrPy", MEtXYcorParams.definitions().formula().c_str());
 
    for (unsigned i(0); i < MEtXYcorParams.record(0).nParameters(); i++) {
      formula_x_->SetParameter(i, MEtXYcorParams.record(0).parameter(i));
    }
    for (unsigned i(0); i < MEtXYcorParams.record(1).nParameters(); i++) {
      formula_y_->SetParameter(i, MEtXYcorParams.record(1).parameter(i));
    }
 
    corx -= formula_x_->Eval(val);
    cory -= formula_y_->Eval(val);
 
  }  //end loop over corrections
 
  metCorr->mex = corx;
  metCorr->mey = cory;
  evt.put(std::move(metCorr), "");
}

References funct::abs(), HltBtagPostValidation_cff::c, counts_, MEtXYcorrectParameters::definitions(), Exception, MEtXYcorrectParametersCollection::findLabel(), MEtXYcorrectParameters::Definitions::formula(), formula_x_, formula_y_, edm::Event::getByToken(), edm::EventSetup::getHandle(), metFilters_cff::goodVertices, mps_fire::i, MEtXYcorrectParametersCollection::isShiftData(), MEtXYcorrectParametersCollection::isShiftDY(), MEtXYcorrectParametersCollection::isShiftMC(), MEtXYcorrectParametersCollection::isShiftTTJets(), MEtXYcorrectParametersCollection::isShiftWJets(), edm::EDGetTokenT< T >::isUninitialized(), edm::HandleBase::isValid(), relativeConstraints::keys, mIsData, mMEtXYcorParaColl_, eostools::move(), mSampleType, MEtXYcorrectParameters::Record::nParameters(), MEtXYcorrectParameters::Record::parameter(), pfLinker_cff::particleFlow, MEtXYcorrectParameters::Definitions::parVar(), pflow_, position, MEtXYcorrectParameters::Definitions::PtclType(), edm::Event::put(), MEtXYcorrectParameters::record(), rho, AlCaHLTBitMon_QueryRunRegistry::string, sumPt_, translateTypeToAbsPdgId(), heppy_batch::val, MEtXYcorrectParametersCollection::validKeys(), vertices_, mps_merge::weight, HLT_FULL_cff::weights, weightsToken_, MEtXYcorrectParameters::Record::xMax(), MEtXYcorrectParameters::Record::xMin(), and z.

translateTypeToAbsPdgId()

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

staticprivate

Definition at line 16 of file MultShiftMETcorrDBInputProducer.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> MultShiftMETcorrDBInputProducer::cfgCorrParameters_

private

Definition at line 54 of file MultShiftMETcorrDBInputProducer.h.

counts_

int MultShiftMETcorrDBInputProducer::counts_

private

Definition at line 57 of file MultShiftMETcorrDBInputProducer.h.

Referenced by produce().

etaMax_

std::vector<double> MultShiftMETcorrDBInputProducer::etaMax_

private

Definition at line 56 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer().

etaMin_

std::vector<double> MultShiftMETcorrDBInputProducer::etaMin_

private

Definition at line 56 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer().

formula_x_

std::unique_ptr<TF1> MultShiftMETcorrDBInputProducer::formula_x_

private

Definition at line 59 of file MultShiftMETcorrDBInputProducer.h.

Referenced by produce().

formula_y_

std::unique_ptr<TF1> MultShiftMETcorrDBInputProducer::formula_y_

private

Definition at line 60 of file MultShiftMETcorrDBInputProducer.h.

Referenced by produce().

mIsData

bool MultShiftMETcorrDBInputProducer::mIsData

private

Definition at line 51 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().

mMEtXYcorParaColl_

edm::ESGetToken<MEtXYcorrectParametersCollection, MEtXYcorrectRecord> MultShiftMETcorrDBInputProducer::mMEtXYcorParaColl_

private

Definition at line 47 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().

moduleLabel_

std::string MultShiftMETcorrDBInputProducer::moduleLabel_

private

Definition at line 48 of file MultShiftMETcorrDBInputProducer.h.

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

mPayloadName

std::string MultShiftMETcorrDBInputProducer::mPayloadName

private

Definition at line 49 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer().

mSampleType

std::string MultShiftMETcorrDBInputProducer::mSampleType

private

Definition at line 50 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().

pflow_

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

private

Definition at line 45 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().

sumPt_

double MultShiftMETcorrDBInputProducer::sumPt_

private

Definition at line 58 of file MultShiftMETcorrDBInputProducer.h.

Referenced by produce().

vertices_

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

private

Definition at line 46 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().

weightsToken_

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

private

Definition at line 52 of file MultShiftMETcorrDBInputProducer.h.

Referenced by MultShiftMETcorrDBInputProducer(), and produce().