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)

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

	~HLTTrackMETProducer () override

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

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::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

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 41 of file HLTTrackMETProducer.h.

Constructor & Destructor Documentation

◆ HLTTrackMETProducer()

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

explicit

Definition at line 17 of file HLTTrackMETProducer.cc.

     : 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>();
 }

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

◆ ~HLTTrackMETProducer()

HLTTrackMETProducer::~HLTTrackMETProducer ( )

overridedefault

Member Function Documentation

◆ fillDescriptions()

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

static

Definition at line 44 of file HLTTrackMETProducer.cc.

                                                                                      {
   // 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);
 }

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

◆ produce()

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

override

Definition at line 64 of file HLTTrackMETProducer.cc.

                                                                                {
   // Create a pointer to the products
   std::unique_ptr<reco::METCollection> result(new reco::METCollection());
  
   if (pfCandidatesLabel_.label().empty())
     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->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;
       }
     }
  
   } else if (useTracks_ && !tracks->empty()) {
     for (auto const& j : *tracks) {
       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->empty()) {
     for (auto const& j : *pfRecTracks) {
       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->empty()) {
     for (auto const& j : *pfCandidates) {
       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 (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));
 }

References funct::abs(), funct::cos(), PVValHelper::eta, excludePFMuons_, iEvent, dqmiolumiharvest::j, singleTopDQM_cfi::jets, edm::InputTag::label(), m_theJetToken, m_thePFCandidateToken, m_theRecTrackToken, m_theTrackToken, maxEtaJet_, minNJet_, minPtJet_, eostools::move(), p4, zmumugammaAnalyzer_cfi::pfCandidates, pfCandidatesLabel_, phi, DiDispStaMuonMonitor_cfi::pt, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, mps_fire::result, funct::sin(), mathSSE::sqrt(), PDWG_EXOHSCP_cff::tracks, usePFCandidates_, usePFCandidatesCharged_, usePFRecTracks_, usePt_, useTracks_, and badGlobalMuonTaggersAOD_cff::vtx.

Member Data Documentation

◆ excludePFMuons_

bool HLTTrackMETProducer::excludePFMuons_

private

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

Definition at line 70 of file HLTTrackMETProducer.h.

Referenced by produce().

◆ jetsLabel_

edm::InputTag HLTTrackMETProducer::jetsLabel_

private

Input jet, track, PFRecTrack, PFCandidate collections.

Definition at line 83 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer().

◆ m_theJetToken

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

private

Definition at line 88 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

◆ m_thePFCandidateToken

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

private

Definition at line 91 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

◆ m_theRecTrackToken

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

private

Definition at line 90 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

◆ m_theTrackToken

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

private

Definition at line 89 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

◆ maxEtaJet_

double HLTTrackMETProducer::maxEtaJet_

private

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

Definition at line 80 of file HLTTrackMETProducer.h.

Referenced by produce().

◆ minNJet_

int HLTTrackMETProducer::minNJet_

private

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

Definition at line 74 of file HLTTrackMETProducer.h.

Referenced by produce().

◆ minPtJet_

double HLTTrackMETProducer::minPtJet_

private

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

Definition at line 77 of file HLTTrackMETProducer.h.

Referenced by produce().

◆ pfCandidatesLabel_

edm::InputTag HLTTrackMETProducer::pfCandidatesLabel_

private

Definition at line 86 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

◆ pfRecTracksLabel_

edm::InputTag HLTTrackMETProducer::pfRecTracksLabel_

private

Definition at line 85 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer().

◆ tracksLabel_

edm::InputTag HLTTrackMETProducer::tracksLabel_

private

Definition at line 84 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer().

◆ usePFCandidates_

bool HLTTrackMETProducer::usePFCandidates_

private

Use PF candidates as input instead of jets.

Definition at line 66 of file HLTTrackMETProducer.h.

Referenced by produce().

◆ usePFCandidatesCharged_

bool HLTTrackMETProducer::usePFCandidatesCharged_

private

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

Definition at line 63 of file HLTTrackMETProducer.h.

Referenced by produce().

◆ usePFRecTracks_

bool HLTTrackMETProducer::usePFRecTracks_

private

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

Definition at line 59 of file HLTTrackMETProducer.h.

Referenced by produce().

◆ usePt_

bool HLTTrackMETProducer::usePt_

private

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

Definition at line 51 of file HLTTrackMETProducer.h.

Referenced by produce().

◆ useTracks_

bool HLTTrackMETProducer::useTracks_

private

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

Definition at line 55 of file HLTTrackMETProducer.h.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HLTTrackMETProducer()

◆ ~HLTTrackMETProducer()

Member Function Documentation

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ excludePFMuons_

◆ jetsLabel_

◆ m_theJetToken

◆ m_thePFCandidateToken

◆ m_theRecTrackToken

◆ m_theTrackToken

◆ maxEtaJet_

◆ minNJet_

◆ minPtJet_

◆ pfCandidatesLabel_

◆ pfRecTracksLabel_

◆ tracksLabel_

◆ usePFCandidates_

◆ usePFCandidatesCharged_

◆ usePFRecTracks_

◆ usePt_

◆ useTracks_