cms::METSignificanceProducer Class Reference

#include <METSignificanceProducer.h>

Inheritance diagram for cms::METSignificanceProducer:

Public Member Functions
	METSignificanceProducer (const edm::ParameterSet &)
void	produce (edm::Event &, const edm::EventSetup &) override
	~METSignificanceProducer () 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 Attributes
std::string	jetResPhiType_
std::string	jetResPtType_
std::string	jetSFType_
std::vector< edm::EDGetTokenT< edm::View< reco::Candidate > > >	lepTokens_
metsig::METSignificance *	metSigAlgo_
edm::EDGetTokenT< edm::View< reco::MET > >	metToken_
edm::EDGetTokenT< edm::View< reco::Candidate > >	pfCandidatesToken_
edm::EDGetTokenT< edm::View< reco::Jet > >	pfjetsToken_
edm::EDGetTokenT< double >	rhoToken_

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

Definition at line 45 of file METSignificanceProducer.h.

Constructor & Destructor Documentation

METSignificanceProducer::METSignificanceProducer ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 15 of file METSignificanceProducer.cc.

References edm::ParameterSet::getParameter(), jetResPhiType_, jetResPtType_, jetSFType_, lepTokens_, metSigAlgo_, METSignificance_cfi::METSignificance, metToken_, pfCandidatesToken_, pfjetsToken_, rhoToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                 {
    std::vector<edm::InputTag> srcLeptonsTags = iConfig.getParameter<std::vector<edm::InputTag> >("srcLeptons");
    for (std::vector<edm::InputTag>::const_iterator it = srcLeptonsTags.begin(); it != srcLeptonsTags.end(); it++) {
      lepTokens_.push_back(consumes<edm::View<reco::Candidate> >(*it));
    }

    pfjetsToken_ = consumes<edm::View<reco::Jet> >(iConfig.getParameter<edm::InputTag>("srcPfJets"));

    metToken_ = consumes<edm::View<reco::MET> >(iConfig.getParameter<edm::InputTag>("srcMet"));
    pfCandidatesToken_ = consumes<edm::View<reco::Candidate> >(iConfig.getParameter<edm::InputTag>("srcPFCandidates"));

    jetSFType_ = iConfig.getParameter<std::string>("srcJetSF");
    jetResPtType_ = iConfig.getParameter<std::string>("srcJetResPt");
    jetResPhiType_ = iConfig.getParameter<std::string>("srcJetResPhi");
    rhoToken_ = consumes<double>(iConfig.getParameter<edm::InputTag>("srcRho"));

    metSigAlgo_ = new metsig::METSignificance(iConfig);

    produces<double>("METSignificance");
    produces<math::Error<2>::type>("METCovariance");
  }

METSignificanceProducer::~METSignificanceProducer ( )

override

Definition at line 37 of file METSignificanceProducer.cc.

References metSigAlgo_.

{ delete metSigAlgo_; }

Member Function Documentation

void METSignificanceProducer::produce ( edm::Event &  event, const edm::EventSetup &  setup  )

override

Definition at line 40 of file METSignificanceProducer.cc.

References DEFINE_FWK_MODULE, JME::JetResolution::get(), JME::JetResolutionScaleFactor::get(), metsig::METSignificance::getCovariance(), metsig::METSignificance::getSignificance(), edm::EventBase::isRealData(), jetResPhiType_, jetResPtType_, singleTopDQM_cfi::jets, jetSFType_, lepTokens_, HLT_2018_cff::leptons, BTaggingMonitor_cfi::met, metSigAlgo_, metToken_, eostools::move(), zmumugammaAnalyzer_cfi::pfCandidates, pfCandidatesToken_, pfjetsToken_, rhoToken_, and met_cff::significance.

                                                                                   {
    //
    // met
    //
    edm::Handle<edm::View<reco::MET> > metHandle;
    event.getByToken(metToken_, metHandle);
    const reco::MET& met = (*metHandle)[0];

    //
    // candidates
    //
    edm::Handle<reco::CandidateView> pfCandidates;
    event.getByToken(pfCandidatesToken_, pfCandidates);

    //
    // leptons
    //
    std::vector<edm::Handle<reco::CandidateView> > leptons;
    for (std::vector<edm::EDGetTokenT<edm::View<reco::Candidate> > >::const_iterator srcLeptons_i = lepTokens_.begin();
         srcLeptons_i != lepTokens_.end();
         ++srcLeptons_i) {
      edm::Handle<reco::CandidateView> leptons_i;
      event.getByToken(*srcLeptons_i, leptons_i);
      leptons.push_back(leptons_i);
    }

    //
    // jets
    //
    edm::Handle<edm::View<reco::Jet> > jets;
    event.getByToken(pfjetsToken_, jets);

    edm::Handle<double> rho;
    event.getByToken(rhoToken_, rho);

    JME::JetResolution resPtObj = JME::JetResolution::get(setup, jetResPtType_);
    JME::JetResolution resPhiObj = JME::JetResolution::get(setup, jetResPhiType_);
    JME::JetResolutionScaleFactor resSFObj = JME::JetResolutionScaleFactor::get(setup, jetSFType_);

    //
    // compute the significance
    //
    double sumPtUnclustered = 0;
    const reco::METCovMatrix cov = metSigAlgo_->getCovariance(
        *jets, leptons, pfCandidates, *rho, resPtObj, resPhiObj, resSFObj, event.isRealData(), sumPtUnclustered);
    double sig = metSigAlgo_->getSignificance(cov, met);

    auto significance = std::make_unique<double>();
    (*significance) = sig;

    auto covPtr = std::make_unique<math::Error<2>::type>();
    (*covPtr)(0, 0) = cov(0, 0);
    (*covPtr)(1, 0) = cov(1, 0);
    (*covPtr)(1, 1) = cov(1, 1);

    event.put(std::move(covPtr), "METCovariance");
    event.put(std::move(significance), "METSignificance");
  }

Member Data Documentation

std::string cms::METSignificanceProducer::jetResPhiType_

private

Definition at line 61 of file METSignificanceProducer.h.

Referenced by METSignificanceProducer(), and produce().

std::string cms::METSignificanceProducer::jetResPtType_

private

Definition at line 60 of file METSignificanceProducer.h.

Referenced by METSignificanceProducer(), and produce().

std::string cms::METSignificanceProducer::jetSFType_

private

Definition at line 59 of file METSignificanceProducer.h.

Referenced by METSignificanceProducer(), and produce().

std::vector<edm::EDGetTokenT<edm::View<reco::Candidate> > > cms::METSignificanceProducer::lepTokens_

private

Definition at line 57 of file METSignificanceProducer.h.

Referenced by METSignificanceProducer(), and produce().

metsig::METSignificance* cms::METSignificanceProducer::metSigAlgo_

private

Definition at line 63 of file METSignificanceProducer.h.

Referenced by METSignificanceProducer(), produce(), and ~METSignificanceProducer().

edm::EDGetTokenT<edm::View<reco::MET> > cms::METSignificanceProducer::metToken_

private

Definition at line 55 of file METSignificanceProducer.h.

Referenced by METSignificanceProducer(), and produce().

edm::EDGetTokenT<edm::View<reco::Candidate> > cms::METSignificanceProducer::pfCandidatesToken_

private

Definition at line 56 of file METSignificanceProducer.h.

Referenced by METSignificanceProducer(), and produce().

edm::EDGetTokenT<edm::View<reco::Jet> > cms::METSignificanceProducer::pfjetsToken_

private

Definition at line 54 of file METSignificanceProducer.h.

Referenced by METSignificanceProducer(), and produce().

edm::EDGetTokenT<double> cms::METSignificanceProducer::rhoToken_

private

Definition at line 58 of file METSignificanceProducer.h.

Referenced by METSignificanceProducer(), and produce().