#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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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_

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

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 46 of file METSignificanceProducer.h.

Constructor & Destructor Documentation

◆ METSignificanceProducer()

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

explicit

Definition at line 15 of file METSignificanceProducer.cc.

                                                                                  {
     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"));
  
     edm::InputTag srcWeights = iConfig.getParameter<edm::InputTag>("srcWeights");
     if (!srcWeights.label().empty())
       weightsToken_ = consumes<edm::ValueMap<float>>(srcWeights);
  
     metSigAlgo_ = new metsig::METSignificance(iConfig);
  
     produces<double>("METSignificance");
     produces<math::Error<2>::type>("METCovariance");
   }

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

◆ ~METSignificanceProducer()

METSignificanceProducer::~METSignificanceProducer ( )

override

Definition at line 41 of file METSignificanceProducer.cc.

41 { delete metSigAlgo_; }

References metSigAlgo_.

Member Function Documentation

◆ produce()

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

override

Definition at line 44 of file METSignificanceProducer.cc.

                                                                                    {
     //
     // 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);
  
     edm::Handle<edm::ValueMap<float>> weights;
     if (!weightsToken_.isUninitialized())
       event.getByToken(weightsToken_, weights);
  
     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 edm::ValueMap<float>* weightsPtr = nullptr;
     if (met.isWeighted()) {
       if (weightsToken_.isUninitialized())
         throw cms::Exception("InvalidInput")
             << "MET is weighted (e.g. PUPPI), but no weights given in METSignificanceProducer\n";
       weightsPtr = &*weights;
     }
     reco::METCovMatrix cov = metSigAlgo_->getCovariance(*jets,
                                                         leptons,
                                                         pfCandidates,
                                                         *rho,
                                                         resPtObj,
                                                         resPhiObj,
                                                         resSFObj,
                                                         event.isRealData(),
                                                         sumPtUnclustered,
                                                         weightsPtr);
     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");
   }