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

#include <HLTTrackMETProducer.h>

Inheritance diagram for HLTTrackMETProducer:

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

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

	~HLTTrackMETProducer ()

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

	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 &	itemsToGetFromEvent () 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, track, PFRecTrack, PFCandidate collections. More...

edm::EDGetTokenT< reco::JetView >	m_theJetToken

edm::EDGetTokenT < reco::PFCandidateCollection >	m_thePFCandidateToken

edm::EDGetTokenT < reco::PFRecTrackCollection >	m_theRecTrackToken

edm::EDGetTokenT < reco::TrackCollection >	m_theTrackToken

double	maxEtaJet_
	Maximum (abs) eta requirement for jets (or objects used as input) More...

int	minNJet_

double	minPtJet_
	Minimum pt requirement for jets (or objects used as input) More...

edm::InputTag	pfCandidatesLabel_

edm::InputTag	pfRecTracksLabel_

edm::InputTag	tracksLabel_

bool	usePFCandidates_
	Use PF candidates as input instead of jets. More...

bool	usePFCandidatesCharged_

bool	usePFRecTracks_

bool	usePt_

bool	useTracks_

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 MHT (or MET)

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

MHT (or MET) is calculated using one of the input collections:

Calo or PF jets
tracks
PF tracks
PF charged candidates
PF candidates MHT can include or exclude the contribution from muons.

See header file for documentation

Author: Steven Lowette

Definition at line 42 of file HLTTrackMETProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 18 of file HLTTrackMETProducer.cc.

References jetsLabel_, m_theJetToken, m_thePFCandidateToken, m_theRecTrackToken, m_theTrackToken, pfCandidatesLabel_, pfRecTracksLabel_, and tracksLabel_.

                                                                         :
   usePt_                  ( iConfig.getParameter<bool>("usePt") ),
   useTracks_              ( iConfig.getParameter<bool>("useTracks") ),
   usePFRecTracks_         ( iConfig.getParameter<bool>("usePFRecTracks") ),
   usePFCandidatesCharged_ ( iConfig.getParameter<bool>("usePFCandidatesCharged") ),
   usePFCandidates_        ( iConfig.getParameter<bool>("usePFCandidates") ),
   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") ),
   tracksLabel_            ( iConfig.getParameter<edm::InputTag>("tracksLabel") ),
   pfRecTracksLabel_       ( iConfig.getParameter<edm::InputTag>("pfRecTracksLabel") ),
   pfCandidatesLabel_      ( iConfig.getParameter<edm::InputTag>("pfCandidatesLabel") ) {
     m_theJetToken = consumes<edm::View<reco::Jet>>(jetsLabel_);
     m_theTrackToken = consumes<reco::TrackCollection>(tracksLabel_);
     m_theRecTrackToken = consumes<reco::PFRecTrackCollection>(pfRecTracksLabel_);
     m_thePFCandidateToken = consumes<reco::PFCandidateCollection>(pfCandidatesLabel_);
 
     // Register the products
     produces<reco::METCollection>();
 }

HLTTrackMETProducer::~HLTTrackMETProducer ( )

Definition at line 42 of file HLTTrackMETProducer.cc.

42 {}

Member Function Documentation

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

static

Definition at line 45 of file HLTTrackMETProducer.cc.

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

                                                                                       {
     // Current default is for hltPFMET
     edm::ParameterSetDescription desc;
     desc.add<bool>("usePt", true);
     desc.add<bool>("useTracks", false);
     desc.add<bool>("usePFRecTracks", false);
     desc.add<bool>("usePFCandidatesCharged", true);
     desc.add<bool>("usePFCandidates", false);
     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>("tracksLabel",  edm::InputTag("hltL3Muons"));
     desc.add<edm::InputTag>("pfRecTracksLabel",  edm::InputTag("hltLightPFTracks"));
     desc.add<edm::InputTag>("pfCandidatesLabel",  edm::InputTag("hltParticleFlow"));
     descriptions.add("hltTrackMETProducer", desc);
 }

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

virtual

Implements edm::stream::EDProducerBase.

Definition at line 65 of file HLTTrackMETProducer.cc.

References funct::abs(), funct::cos(), eta, excludePFMuons_, edm::Event::getByToken(), j, fwrapper::jets, edm::InputTag::label(), m_theJetToken, m_thePFCandidateToken, m_theRecTrackToken, m_theTrackToken, maxEtaJet_, minNJet_, minPtJet_, eostools::move(), p4, reco::tau::pfCandidates(), pfCandidatesLabel_, phi, EnergyCorrector::pt, edm::Event::put(), mps_fire::result, funct::sin(), mathSSE::sqrt(), testEve_cfg::tracks, usePFCandidates_, usePFCandidatesCharged_, usePFRecTracks_, usePt_, and useTracks_.

                                                                                {
 
     // Create a pointer to the products
     std::unique_ptr<reco::METCollection> result(new reco::METCollection());
 
     if (pfCandidatesLabel_.label() == "")
         excludePFMuons_ = false;
 
     bool useJets = !useTracks_ && !usePFRecTracks_ && !usePFCandidatesCharged_ && !usePFCandidates_;
     if (!useJets) {
         minNJet_ = 0;
     }
 
     edm::Handle<reco::JetView> jets;
     if (useJets) iEvent.getByToken(m_theJetToken, jets);
 
     edm::Handle<reco::TrackCollection> tracks;
     if (useTracks_) iEvent.getByToken(m_theTrackToken, tracks);
 
     edm::Handle<reco::PFRecTrackCollection> pfRecTracks;
     if (usePFRecTracks_) iEvent.getByToken(m_theRecTrackToken, pfRecTracks);
 
     edm::Handle<reco::PFCandidateCollection> pfCandidates;
     if (excludePFMuons_ || usePFCandidatesCharged_ || usePFCandidates_)
         iEvent.getByToken(m_thePFCandidateToken, pfCandidates);
 
     int nj = 0;
     double sumet = 0., mhx = 0., mhy = 0.;
 
     if (useJets && 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 > minPtJet_ && std::abs(eta) < maxEtaJet_) {
                 mhx -= px;
                 mhy -= py;
                 sumet += pt;
                 ++nj;
             }
         }
 
     } else if (useTracks_ && tracks->size() > 0) {
         for (reco::TrackCollection::const_iterator j = tracks->begin(); j != tracks->end(); ++j) {
             double pt = j->pt();
             double px = j->px();
             double py = j->py();
             double eta = j->eta();
 
             if (pt > minPtJet_ && std::abs(eta) < maxEtaJet_) {
                 mhx -= px;
                 mhy -= py;
                 sumet += pt;
                 ++nj;
             }
         }
 
     } else if (usePFRecTracks_ && pfRecTracks->size() > 0) {
         for (reco::PFRecTrackCollection::const_iterator j = pfRecTracks->begin(); j != pfRecTracks->end(); ++j) {
             double pt = j->trackRef()->pt();
             double px = j->trackRef()->px();
             double py = j->trackRef()->py();
             double eta = j->trackRef()->eta();
 
             if (pt > minPtJet_ && std::abs(eta) < maxEtaJet_) {
                 mhx -= px;
                 mhy -= py;
                 sumet += pt;
                 ++nj;
             }
         }
 
     } else if ((usePFCandidatesCharged_ || usePFCandidates_) && pfCandidates->size() > 0) {
         for (reco::PFCandidateCollection::const_iterator j = pfCandidates->begin(); j != pfCandidates->end(); ++j) {
             if (usePFCandidatesCharged_ && j->charge() == 0)  continue;
             double pt = j->pt();
             double px = j->px();
             double py = j->py();
             double eta = j->eta();
 
             if (pt > minPtJet_ && std::abs(eta) < maxEtaJet_) {
                 mhx -= px;
                 mhy -= py;
                 sumet += pt;
                 ++nj;
             }
         }
     }
 
     if (excludePFMuons_) {
         for (reco::PFCandidateCollection::const_iterator j = pfCandidates->begin(); j != pfCandidates->end(); ++j) {
             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 HLTTrackMETProducer::excludePFMuons_

private

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

Definition at line 71 of file HLTTrackMETProducer.h.

Referenced by produce().

edm::InputTag HLTTrackMETProducer::jetsLabel_

private

Input jet, track, PFRecTrack, PFCandidate collections.

Definition at line 84 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer().

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

private

Definition at line 89 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

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

private

Definition at line 92 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

edm::EDGetTokenT<reco::PFRecTrackCollection> HLTTrackMETProducer::m_theRecTrackToken

private

Definition at line 91 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

edm::EDGetTokenT<reco::TrackCollection> HLTTrackMETProducer::m_theTrackToken

private

Definition at line 90 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

double HLTTrackMETProducer::maxEtaJet_

private

Maximum (abs) eta requirement for jets (or objects used as input)

Definition at line 81 of file HLTTrackMETProducer.h.

Referenced by produce().

int HLTTrackMETProducer::minNJet_

private

Minimum number of jets passing pt and eta requirements Ignored if jets are not used as input

Definition at line 75 of file HLTTrackMETProducer.h.

Referenced by produce().

double HLTTrackMETProducer::minPtJet_

private

Minimum pt requirement for jets (or objects used as input)

Definition at line 78 of file HLTTrackMETProducer.h.

Referenced by produce().

edm::InputTag HLTTrackMETProducer::pfCandidatesLabel_

private

Definition at line 87 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

edm::InputTag HLTTrackMETProducer::pfRecTracksLabel_

private

Definition at line 86 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer().

edm::InputTag HLTTrackMETProducer::tracksLabel_

private

Definition at line 85 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer().

bool HLTTrackMETProducer::usePFCandidates_

private

Use PF candidates as input instead of jets.

Definition at line 67 of file HLTTrackMETProducer.h.

Referenced by produce().

bool HLTTrackMETProducer::usePFCandidatesCharged_

private

Use PF charged candidates as input instead of jets. If true, it overrides usePFCandidates_.

Definition at line 64 of file HLTTrackMETProducer.h.

Referenced by produce().

bool HLTTrackMETProducer::usePFRecTracks_

private

Use PF tracks as input instead of jets. If true, it overrides usePFCandidatesCharged_ & usePFCandidates_.

Definition at line 60 of file HLTTrackMETProducer.h.

Referenced by produce().

bool HLTTrackMETProducer::usePt_

private

Use pt; otherwise, use et. Ignored if jets are not used as input.

Definition at line 52 of file HLTTrackMETProducer.h.

Referenced by produce().

bool HLTTrackMETProducer::useTracks_

private

Use reco tracks as input instead of jets. If true, it overrides usePFRecTracks, usePFCandidatesCharged_ & usePFCandidates_.

Definition at line 56 of file HLTTrackMETProducer.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