This produces a reco::MET object that stores MHT (or MET) More...

#include <HLTMhtProducer.h>

Inheritance diagram for HLTMhtProducer:

Public Member Functions
	HLTMhtProducer (const edm::ParameterSet &iConfig)

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

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

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Attributes
bool	excludePFMuons_

edm::InputTag	jetsLabel_
	Input jet, PFCandidate collections. More...

edm::EDGetTokenT< reco::JetView >	m_theJetToken

edm::EDGetTokenT< reco::PFCandidateCollection >	m_thePFCandidateToken

double	maxEtaJet_
	Maximum (abs) eta requirement for jets. More...

int	minNJet_
	Minimum number of jets passing pt and eta requirements. More...

double	minPtJet_
	Minimum pt requirement for jets. More...

edm::InputTag	pfCandidatesLabel_

bool	usePt_
	Use pt; otherwise, use et. More...

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

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

This produces a reco::MET object that stores MHT (or MET)

Author: Steven Lowette; Michele de Gruttola, Jia Fu Low (Nov 2013)

MHT (or MET) is calculated using input CaloJet or PFJet collection. MHT can include or exclude the contribution from muons.

See header file for documentation

Author: Steven Lowette

Definition at line 34 of file HLTMhtProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file HLTMhtProducer.cc.

References excludePFMuons_, jetsLabel_, edm::InputTag::label(), m_theJetToken, m_thePFCandidateToken, pfCandidatesLabel_, and ~HLTMhtProducer().

                                                               :
   usePt_                  ( iConfig.getParameter<bool>("usePt") ),
   excludePFMuons_         ( iConfig.getParameter<bool>("excludePFMuons") ),
   minNJet_                ( iConfig.getParameter<int>("minNJet") ),
   minPtJet_               ( iConfig.getParameter<double>("minPtJet") ),
   maxEtaJet_              ( iConfig.getParameter<double>("maxEtaJet") ),
   jetsLabel_              ( iConfig.getParameter<edm::InputTag>("jetsLabel") ),
   pfCandidatesLabel_      ( iConfig.getParameter<edm::InputTag>("pfCandidatesLabel") ) {
     m_theJetToken = consumes<edm::View<reco::Jet>>(jetsLabel_);
     if (pfCandidatesLabel_.label() == "") excludePFMuons_ = false;
     if (excludePFMuons_) m_thePFCandidateToken = consumes<reco::PFCandidateCollection>(pfCandidatesLabel_);
 
     // Register the products
     produces<reco::METCollection>();
 }

HLTMhtProducer::~HLTMhtProducer ( )

overridedefault

Referenced by HLTMhtProducer().

Member Function Documentation

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

static

Definition at line 38 of file HLTMhtProducer.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

                                                                                  {
     // Current default is for hltPFMET
     edm::ParameterSetDescription desc;
     desc.add<bool>("usePt", true);
     desc.add<bool>("excludePFMuons", false);
     desc.add<int>("minNJet",0);
     desc.add<double>("minPtJet", 0.);
     desc.add<double>("maxEtaJet", 999.);
     desc.add<edm::InputTag>("jetsLabel", edm::InputTag("hltAntiKT4PFJets"));
     desc.add<edm::InputTag>("pfCandidatesLabel",  edm::InputTag("hltParticleFlow"));
     descriptions.add("hltMhtProducer", desc);
 }

void HLTMhtProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 52 of file HLTMhtProducer.cc.

References funct::abs(), edm::View< T >::begin(), funct::cos(), edm::View< T >::empty(), edm::View< T >::end(), PVValHelper::eta, excludePFMuons_, edm::Event::getByToken(), fwrapper::jets, m_theJetToken, m_thePFCandidateToken, maxEtaJet_, minNJet_, minPtJet_, eostools::move(), p4, slimmedMuons_cfi::pfCandidates, phi, EnergyCorrector::pt, edm::Event::put(), mps_fire::result, funct::sin(), mathSSE::sqrt(), usePt_, and badGlobalMuonTaggersAOD_cff::vtx.

                                                                           {
 
     // Create a pointer to the products
     std::unique_ptr<reco::METCollection> result(new reco::METCollection());
 
     edm::Handle<reco::JetView> jets;
     iEvent.getByToken(m_theJetToken, jets);
 
     edm::Handle<reco::PFCandidateCollection> pfCandidates;
     if (excludePFMuons_)
         iEvent.getByToken(m_thePFCandidateToken, pfCandidates);
 
     int nj = 0;
     double sumet = 0., mhx = 0., mhy = 0.;
 
     if (!jets->empty()) {
         for(reco::JetView::const_iterator j = jets->begin(); j != jets->end(); ++j) {
             double pt = usePt_ ? j->pt() : j->et();
             double eta = j->eta();
             double phi = j->phi();
             double px = usePt_ ? j->px() : j->et() * cos(phi);
             double py = usePt_ ? j->py() : j->et() * sin(phi);
 
             if (pt > minPtJet_ && std::abs(eta) < maxEtaJet_) {
                 mhx -= px;
                 mhy -= py;
                 sumet += pt;
                 ++nj;
             }
         }
     }
 
     if (excludePFMuons_) {
         for (auto const & j : *pfCandidates) {
             if (std::abs(j.pdgId()) == 13) {
                 mhx += j.px();
                 mhy += j.py();
             }
         }
     }
 
     if (nj < minNJet_) { sumet = 0; mhx = 0; mhy = 0; }
 
     reco::MET::LorentzVector p4(mhx, mhy, 0, sqrt(mhx*mhx + mhy*mhy));
     reco::MET::Point vtx(0, 0, 0);
     reco::MET mht(sumet, p4, vtx);
     result->push_back(mht);
 
     // Put the products into the Event
     iEvent.put(std::move(result));
 }

Member Data Documentation

bool HLTMhtProducer::excludePFMuons_

private

Exclude PF muons in the MHT calculation (but not HT) Ignored if pfCandidatesLabel_ is empty.

Definition at line 47 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer(), and produce().

edm::InputTag HLTMhtProducer::jetsLabel_

private

Input jet, PFCandidate collections.

Definition at line 59 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer().

edm::EDGetTokenT<reco::JetView> HLTMhtProducer::m_theJetToken

private

Definition at line 62 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer(), and produce().

edm::EDGetTokenT<reco::PFCandidateCollection> HLTMhtProducer::m_thePFCandidateToken

private

Definition at line 63 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer(), and produce().

double HLTMhtProducer::maxEtaJet_

private

Maximum (abs) eta requirement for jets.

Definition at line 56 of file HLTMhtProducer.h.

Referenced by produce().

int HLTMhtProducer::minNJet_

private

Minimum number of jets passing pt and eta requirements.

Definition at line 50 of file HLTMhtProducer.h.

Referenced by produce().

double HLTMhtProducer::minPtJet_

private

Minimum pt requirement for jets.

Definition at line 53 of file HLTMhtProducer.h.

Referenced by produce().

edm::InputTag HLTMhtProducer::pfCandidatesLabel_

private

Definition at line 60 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer().

bool HLTMhtProducer::usePt_

private

Use pt; otherwise, use et.

Definition at line 43 of file HLTMhtProducer.h.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation