#include <PFMETProducer.h>

Inheritance diagram for cms::PFMETProducer:

Public Member Functions
	PFMETProducer (const edm::ParameterSet &)

void	produce (edm::Event &, const edm::EventSetup &) override

	~PFMETProducer () 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

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
reco::METCovMatrix	getMETCovMatrix (const edm::Event &event, const edm::EventSetup &, const edm::Handle< edm::View< reco::Candidate >> &input) const

Private Attributes
bool	applyWeight_

bool	calculateSignificance_

double	globalThreshold_

edm::EDGetTokenT< edm::View< reco::Candidate > >	inputToken_

std::string	jetResPhiType_

std::string	jetResPtType_

std::string	jetSFType_

double	jetThreshold_

edm::EDGetTokenT< edm::View< reco::Jet > >	jetToken_

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

metsig::METSignificance *	metSigAlgo_

edm::EDGetTokenT< double >	rhoToken_

edm::InputTag	src_

edm::ValueMap< float > const *	weights_

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 57 of file PFMETProducer.h.

Constructor & Destructor Documentation

◆ PFMETProducer()

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

explicit

Definition at line 17 of file PFMETProducer.cc.

       : src_(iConfig.getParameter<edm::InputTag>("src")),
         inputToken_(consumes<edm::View<reco::Candidate>>(src_)),
         calculateSignificance_(iConfig.getParameter<bool>("calculateSignificance")),
         globalThreshold_(iConfig.getParameter<double>("globalThreshold")),
         applyWeight_(iConfig.getParameter<bool>("applyWeight")),
         weights_(nullptr) {
     if (applyWeight_) {
       edm::InputTag srcWeights = iConfig.getParameter<edm::InputTag>("srcWeights");
       if (srcWeights.label().empty())
         throw cms::Exception("InvalidInput") << "applyWeight set to True, but no weights given in PFMETProducer\n";
       if (srcWeights.label() == src_.label())
         edm::LogWarning("PFMETProducer")
             << "Particle and weights collection have the same label. You may be applying the same weights twice.\n";
       weightsToken_ = consumes<edm::ValueMap<float>>(srcWeights);
     }
     if (calculateSignificance_) {
       metSigAlgo_ = new metsig::METSignificance(iConfig);
  
       jetToken_ = mayConsume<edm::View<reco::Jet>>(iConfig.getParameter<edm::InputTag>("srcJets"));
       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(mayConsume<edm::View<reco::Candidate>>(*it));
       }
  
       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"));
     }
  
     std::string alias = iConfig.exists("alias") ? iConfig.getParameter<std::string>("alias") : "";
  
     produces<reco::PFMETCollection>().setBranchAlias(alias);
   }

References SiStripOfflineCRack_cfg::alias, applyWeight_, calculateSignificance_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), jetResPhiType_, jetResPtType_, jetSFType_, jetToken_, edm::InputTag::label(), lepTokens_, metSigAlgo_, METSignificance_cfi::METSignificance, rhoToken_, src_, HLT_FULL_cff::srcWeights, AlCaHLTBitMon_QueryRunRegistry::string, and weightsToken_.

◆ ~PFMETProducer()

cms::PFMETProducer::~PFMETProducer ( )

inlineoverride

Definition at line 60 of file PFMETProducer.h.

60 {}

Member Function Documentation

◆ fillDescriptions()

void PFMETProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 132 of file PFMETProducer.cc.

                                                                                  {
     edm::ParameterSetDescription desc;
     desc.add<edm::InputTag>("src", edm::InputTag("particleFlow"));
     desc.add<double>("globalThreshold", 0.);
     desc.add<std::string>("alias", "@module_label");
     desc.add<bool>("calculateSignificance", false);
     desc.addOptional<edm::InputTag>("srcJets");
     desc.addOptional<std::vector<edm::InputTag>>("srcLeptons");
     desc.addOptional<std::string>("srcJetSF");
     desc.addOptional<std::string>("srcJetResPt");
     desc.addOptional<std::string>("srcJetResPhi");
     desc.addOptional<edm::InputTag>("srcRho");
     edm::ParameterSetDescription params;
     params.setAllowAnything();  // FIXME: This still needs to be defined in METSignficance
     desc.addOptional<edm::ParameterSetDescription>("parameters", params);
     edm::ParameterSetDescription desc1 = desc;
     edm::ParameterSetDescription desc2 = desc;
     desc1.add<bool>("applyWeight", false);
     desc1.add<edm::InputTag>("srcWeights", edm::InputTag(""));
     descriptions.add("pfMet", desc1);
     desc2.add<bool>("applyWeight", true);
     desc2.add<edm::InputTag>("srcWeights", edm::InputTag("puppiNoLep"));
     descriptions.add("pfMetPuppi", desc2);
   }

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, HLT_FULL_cff::InputTag, CalibrationSummaryClient_cfi::params, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ getMETCovMatrix()

reco::METCovMatrix PFMETProducer::getMETCovMatrix	(	const edm::Event &	event,
		const edm::EventSetup &	setup,
		const edm::Handle< edm::View< reco::Candidate >> &	input
	)		const

private

Definition at line 86 of file PFMETProducer.cc.

                                                                                                               {
     // 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);
       /*
       for ( reco::CandidateView::const_iterator lepton = leptons_i->begin();
         lepton != leptons_i->end(); ++lepton ) {
         leptons.push_back(*lepton);
       }
      */
     }
  
     // jets
     edm::Handle<edm::View<reco::Jet>> inputJets;
     event.getByToken(jetToken_, inputJets);
  
     JME::JetResolution resPtObj = JME::JetResolution::get(setup, jetResPtType_);
     JME::JetResolution resPhiObj = JME::JetResolution::get(setup, jetResPhiType_);
     JME::JetResolutionScaleFactor resSFObj = JME::JetResolutionScaleFactor::get(setup, jetSFType_);
  
     edm::Handle<double> rho;
     event.getByToken(rhoToken_, rho);
  
     //Compute the covariance matrix and fill it
     double sumPtUnclustered = 0;
     reco::METCovMatrix cov = metSigAlgo_->getCovariance(*inputJets,
                                                         leptons,
                                                         candInput,
                                                         *rho,
                                                         resPtObj,
                                                         resPhiObj,
                                                         resSFObj,
                                                         event.isRealData(),
                                                         sumPtUnclustered,
                                                         weights_);
  
     return cov;
   }

References JME::JetResolution::get(), JME::JetResolutionScaleFactor::get(), metsig::METSignificance::getCovariance(), HLT_FULL_cff::inputJets, jetResPhiType_, jetResPtType_, jetSFType_, jetToken_, lepTokens_, HLT_FULL_cff::leptons, metSigAlgo_, rhoToken_, singleTopDQM_cfi::setup, met_cff::sumPtUnclustered, trackerHitRTTI::vector, and weights_.

Referenced by produce().

◆ produce()

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

override

Definition at line 54 of file PFMETProducer.cc.

                                                                          {
     edm::Handle<edm::View<reco::Candidate>> input;
     event.getByToken(inputToken_, input);
  
     if (applyWeight_)
       weights_ = &event.get(weightsToken_);
  
     METAlgo algo;
     CommonMETData commonMETdata;
     commonMETdata = algo.run(*input.product(), globalThreshold_, weights_);
  
     const math::XYZTLorentzVector p4(commonMETdata.mex, commonMETdata.mey, 0.0, commonMETdata.met);
     const math::XYZPoint vtx(0.0, 0.0, 0.0);
  
     PFSpecificAlgo pf;
     SpecificPFMETData specific;
     specific = pf.run(*input.product(), weights_);
  
     reco::PFMET pfmet(specific, commonMETdata.sumet, p4, vtx);
     pfmet.setIsWeighted(applyWeight_);
  
     if (calculateSignificance_) {
       reco::METCovMatrix sigcov = getMETCovMatrix(event, setup, input);
       pfmet.setSignificanceMatrix(sigcov);
     }
  
     auto pfmetcoll = std::make_unique<reco::PFMETCollection>();
  
     pfmetcoll->push_back(pfmet);
     event.put(std::move(pfmetcoll));
   }

References applyWeight_, calculateSignificance_, edm::ValueMap< T >::get(), getMETCovMatrix(), globalThreshold_, input, inputToken_, CommonMETData::met, CommonMETData::mex, CommonMETData::mey, eostools::move(), p4, packedPFCandidateRefMixer_cfi::pf, reco::MET::setIsWeighted(), reco::MET::setSignificanceMatrix(), singleTopDQM_cfi::setup, timingPdfMaker::specific, CommonMETData::sumet, extraflags_cff::vtx, weights_, and weightsToken_.

Member Data Documentation

◆ applyWeight_

bool cms::PFMETProducer::applyWeight_

private

Definition at line 83 of file PFMETProducer.h.

Referenced by PFMETProducer(), and produce().

◆ calculateSignificance_

bool cms::PFMETProducer::calculateSignificance_

private

Definition at line 71 of file PFMETProducer.h.

Referenced by PFMETProducer(), and produce().

◆ globalThreshold_

double cms::PFMETProducer::globalThreshold_

private

Definition at line 74 of file PFMETProducer.h.

Referenced by produce().

◆ inputToken_

edm::EDGetTokenT<edm::View<reco::Candidate> > cms::PFMETProducer::inputToken_

private

Definition at line 69 of file PFMETProducer.h.

Referenced by produce().

◆ jetResPhiType_

std::string cms::PFMETProducer::jetResPhiType_

private

Definition at line 81 of file PFMETProducer.h.

Referenced by getMETCovMatrix(), and PFMETProducer().

◆ jetResPtType_

std::string cms::PFMETProducer::jetResPtType_

private

Definition at line 80 of file PFMETProducer.h.

Referenced by getMETCovMatrix(), and PFMETProducer().

◆ jetSFType_

std::string cms::PFMETProducer::jetSFType_

private

Definition at line 79 of file PFMETProducer.h.

Referenced by getMETCovMatrix(), and PFMETProducer().

◆ jetThreshold_

double cms::PFMETProducer::jetThreshold_

private

Definition at line 75 of file PFMETProducer.h.

◆ jetToken_

edm::EDGetTokenT<edm::View<reco::Jet> > cms::PFMETProducer::jetToken_

private

Definition at line 77 of file PFMETProducer.h.

Referenced by getMETCovMatrix(), and PFMETProducer().

◆ lepTokens_

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

private

Definition at line 78 of file PFMETProducer.h.

Referenced by getMETCovMatrix(), and PFMETProducer().

◆ metSigAlgo_

metsig::METSignificance* cms::PFMETProducer::metSigAlgo_

private

Definition at line 72 of file PFMETProducer.h.

Referenced by getMETCovMatrix(), and PFMETProducer().

◆ rhoToken_

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

private

Definition at line 82 of file PFMETProducer.h.

Referenced by getMETCovMatrix(), and PFMETProducer().

◆ src_

edm::InputTag cms::PFMETProducer::src_

private

Definition at line 68 of file PFMETProducer.h.

Referenced by PFMETProducer().

◆ weights_

edm::ValueMap<float> const* cms::PFMETProducer::weights_

private

Definition at line 85 of file PFMETProducer.h.

Referenced by getMETCovMatrix(), and produce().

◆ weightsToken_

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

private

Definition at line 84 of file PFMETProducer.h.

Referenced by PFMETProducer(), and produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ PFMETProducer()

◆ ~PFMETProducer()

Member Function Documentation

◆ fillDescriptions()

◆ getMETCovMatrix()

◆ produce()

Member Data Documentation

◆ applyWeight_

◆ calculateSignificance_

◆ globalThreshold_

◆ inputToken_

◆ jetResPhiType_

◆ jetResPtType_

◆ jetSFType_

◆ jetThreshold_

◆ jetToken_

◆ lepTokens_

◆ metSigAlgo_

◆ rhoToken_

◆ src_

◆ weights_

◆ weightsToken_