#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

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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	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_

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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 60 of file PFMETProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 16 of file PFMETProducer.cc.

References met_cff::alias, calculateSignificance_, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), jetResPhiType_, jetResPtType_, jetSFType_, jetToken_, lepTokens_, edm::EDConsumerBase::mayConsume(), metSigAlgo_, METSignificance_cfi::METSignificance, rhoToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

     : inputToken_(consumes<edm::View<reco::Candidate> >(iConfig.getParameter<edm::InputTag>("src")))
     , calculateSignificance_(iConfig.getParameter<bool>("calculateSignificance"))
     , globalThreshold_(iConfig.getParameter<double>("globalThreshold"))
   {
     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);
   }

cms::PFMETProducer::~PFMETProducer ( )

inlineoverride

Definition at line 64 of file PFMETProducer.h.

References input.

64 { }

Member Function Documentation

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 73 of file PFMETProducer.cc.

References DEFINE_FWK_MODULE, JME::JetResolution::get(), JME::JetResolutionScaleFactor::get(), metsig::METSignificance::getCovariance(), JetEnergyShift_cfi::inputJets, edm::EventBase::isRealData(), jetResPhiType_, jetResPtType_, jetSFType_, jetToken_, lepTokens_, AK4GenJetFlavourInfos_cfi::leptons, metSigAlgo_, and rhoToken_.

Referenced by produce().

                                                                                                                                                               {
 
     // 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
     reco::METCovMatrix cov = metSigAlgo_->getCovariance( *inputJets, leptons, candInput, *rho, resPtObj, resPhiObj, resSFObj, event.isRealData() );
 
     return cov;
   }

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

override

Definition at line 43 of file PFMETProducer.cc.

References patPFMETCorrections_cff::algo, calculateSignificance_, getMETCovMatrix(), globalThreshold_, input, inputToken_, CommonMETData::met, CommonMETData::mex, CommonMETData::mey, eostools::move(), p4, packedPFCandidateRefMixer_cfi::pf, edm::Handle< T >::product(), PFSpecificAlgo::run(), METAlgo::run(), reco::MET::setSignificanceMatrix(), CommonMETData::sumet, and badGlobalMuonTaggersAOD_cff::vtx.

   {
     edm::Handle<edm::View<reco::Candidate> > input;
     event.getByToken(inputToken_, input);
 
     METAlgo algo;
     CommonMETData commonMETdata = algo.run(*input.product(), globalThreshold_);
 
     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 = pf.run(*input.product());
 
     reco::PFMET pfmet(specific, commonMETdata.sumet, p4, vtx);
 
     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));
   }