This produces a reco::MET object that stores HT and MHT. More...

#include <HLTHtMhtProducer.h>

Inheritance diagram for HLTHtMhtProducer:

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

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

	~HLTHtMhtProducer ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

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

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, std::unordered_multimap< std::string, edm::ProductResolverIndex > 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	maxEtaJetHt_
	Maximum (abs) eta requirement for jets. More...

double	maxEtaJetMht_

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

int	minNJetMht_

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

double	minPtJetMht_

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
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 HT and MHT.

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

HT & MHT are calculated using input CaloJet or PFJet collection. MHT can include or exclude the contribution from muons. HT is stored as sumet_, MHT is stored as p4_ in the output.

See header file for documentation

Author: Steven Lowette

Definition at line 35 of file HLTHtMhtProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file HLTHtMhtProducer.cc.

References jetsLabel_, m_theJetToken, m_thePFCandidateToken, pfCandidatesLabel_, and ~HLTHtMhtProducer().

                                                                   :
   usePt_                  ( iConfig.getParameter<bool>("usePt") ),
   excludePFMuons_         ( iConfig.getParameter<bool>("excludePFMuons") ),
   minNJetHt_              ( iConfig.getParameter<int>("minNJetHt") ),
   minNJetMht_             ( iConfig.getParameter<int>("minNJetMht") ),
   minPtJetHt_             ( iConfig.getParameter<double>("minPtJetHt") ),
   minPtJetMht_            ( iConfig.getParameter<double>("minPtJetMht") ),
   maxEtaJetHt_            ( iConfig.getParameter<double>("maxEtaJetHt") ),
   maxEtaJetMht_           ( iConfig.getParameter<double>("maxEtaJetMht") ),
   jetsLabel_              ( iConfig.getParameter<edm::InputTag>("jetsLabel") ),
   pfCandidatesLabel_      ( iConfig.getParameter<edm::InputTag>("pfCandidatesLabel") ) {
     m_theJetToken = consumes<edm::View<reco::Jet>>(jetsLabel_);
     m_thePFCandidateToken = consumes<reco::PFCandidateCollection>(pfCandidatesLabel_);
 
     // Register the products
     produces<reco::METCollection>();
 }

HLTHtMhtProducer::~HLTHtMhtProducer ( )

default

Referenced by HLTHtMhtProducer().

Member Function Documentation

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

static

Definition at line 40 of file HLTHtMhtProducer.cc.

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

                                                                                    {
     // Current default is for hltHtMht
     edm::ParameterSetDescription desc;
     desc.add<bool>("usePt", false);
     desc.add<bool>("excludePFMuons", false);
     desc.add<int>("minNJetHt", 0);
     desc.add<int>("minNJetMht", 0);
     desc.add<double>("minPtJetHt", 40.);
     desc.add<double>("minPtJetMht", 30.);
     desc.add<double>("maxEtaJetHt", 3.);
     desc.add<double>("maxEtaJetMht", 5.);
     desc.add<edm::InputTag>("jetsLabel", edm::InputTag("hltCaloJetL1FastJetCorrected"));
     desc.add<edm::InputTag>("pfCandidatesLabel",  edm::InputTag("hltParticleFlow"));
     descriptions.add("hltHtMhtProducer", desc);
 }

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

virtual

Definition at line 57 of file HLTHtMhtProducer.cc.

References funct::abs(), edm::View< T >::begin(), funct::cos(), edm::View< T >::end(), PVValHelper::eta, excludePFMuons_, edm::Event::getByToken(), fwrapper::jets, edm::InputTag::label(), m_theJetToken, m_thePFCandidateToken, maxEtaJetHt_, maxEtaJetMht_, minNJetHt_, minNJetMht_, minPtJetHt_, minPtJetMht_, eostools::move(), p4, slimmedMuons_cfi::pfCandidates, pfCandidatesLabel_, phi, EnergyCorrector::pt, edm::Event::put(), mps_fire::result, funct::sin(), edm::View< T >::size(), mathSSE::sqrt(), usePt_, and badGlobalMuonTaggersAOD_cff::vtx.

                                                                             {
 
     // Create a pointer to the products
     std::unique_ptr<reco::METCollection> result(new reco::METCollection());
 
     if (pfCandidatesLabel_.label() == "")
         excludePFMuons_ = false;
 
     edm::Handle<reco::JetView> jets;
     iEvent.getByToken(m_theJetToken, jets);
 
     edm::Handle<reco::PFCandidateCollection> pfCandidates;
     if (excludePFMuons_)
         iEvent.getByToken(m_thePFCandidateToken, pfCandidates);
 
     int nj_ht = 0, nj_mht = 0;
     double ht = 0., mhx = 0., mhy = 0.;
 
     if (jets->size() > 0) {
         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 > minPtJetHt_ && std::abs(eta) < maxEtaJetHt_) {
                 ht += pt;
                 ++nj_ht;
             }
 
             if (pt > minPtJetMht_ && std::abs(eta) < maxEtaJetMht_) {
                 mhx -= px;
                 mhy -= py;
                 ++nj_mht;
             }
         }
     }
 
     if (excludePFMuons_) {
         for (auto const & j : *pfCandidates) {
             if (std::abs(j.pdgId()) == 13) {
                 mhx += j.px();
                 mhy += j.py();
             }
         }
     }
 
     if (nj_ht  < minNJetHt_ ) { ht = 0; }
     if (nj_mht < minNJetMht_) { 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 htmht(ht, p4, vtx);
     result->push_back(htmht);
 
     // Put the products into the Event
     iEvent.put(std::move(result));
 }

Member Data Documentation

bool HLTHtMhtProducer::excludePFMuons_

private

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

Definition at line 48 of file HLTHtMhtProducer.h.

Referenced by produce().

edm::InputTag HLTHtMhtProducer::jetsLabel_

private

Input jet, PFCandidate collections.

Definition at line 63 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer().

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

private

Definition at line 66 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer(), and produce().

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

private

Definition at line 67 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer(), and produce().

double HLTHtMhtProducer::maxEtaJetHt_

private

Maximum (abs) eta requirement for jets.

Definition at line 59 of file HLTHtMhtProducer.h.

Referenced by produce().

double HLTHtMhtProducer::maxEtaJetMht_

private

Definition at line 60 of file HLTHtMhtProducer.h.

Referenced by produce().

int HLTHtMhtProducer::minNJetHt_

private

Minimum number of jets passing pt and eta requirements.

Definition at line 51 of file HLTHtMhtProducer.h.

Referenced by produce().

int HLTHtMhtProducer::minNJetMht_

private

Definition at line 52 of file HLTHtMhtProducer.h.

Referenced by produce().

double HLTHtMhtProducer::minPtJetHt_

private

Minimum pt requirement for jets.

Definition at line 55 of file HLTHtMhtProducer.h.

Referenced by produce().

double HLTHtMhtProducer::minPtJetMht_

private

Definition at line 56 of file HLTHtMhtProducer.h.

Referenced by produce().

edm::InputTag HLTHtMhtProducer::pfCandidatesLabel_

private

Definition at line 64 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer(), and produce().

bool HLTHtMhtProducer::usePt_

private

Use pt; otherwise, use et.

Definition at line 44 of file HLTHtMhtProducer.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