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 ( const edm::ParameterSet & iConfig )

explicit

Definition at line 17 of file HLTTrackMETProducer.cc.

References jetsLabel_, m_theJetToken, m_thePFCandidateToken, m_theRecTrackToken, m_theTrackToken, pfCandidatesLabel_, pfRecTracksLabel_, tracksLabel_, and ~HLTTrackMETProducer().

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

overridedefault

Referenced by HLTTrackMETProducer().

Member Function Documentation

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

static

Definition at line 44 of file HLTTrackMETProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLT_2018_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
	)

override

Definition at line 64 of file HLTTrackMETProducer.cc.

References funct::abs(), edm::View< T >::begin(), funct::cos(), edm::View< T >::empty(), edm::View< T >::end(), PVValHelper::eta, excludePFMuons_, edm::Event::getByToken(), 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, edm::Event::put(), 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.

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

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

Referenced by produce().

edm::InputTag HLTTrackMETProducer::jetsLabel_

private

Input jet, track, PFRecTrack, PFCandidate collections.

Definition at line 83 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer().

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

private

Definition at line 88 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

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

private

Definition at line 91 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

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

private

Definition at line 90 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

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

private

Definition at line 89 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 80 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 74 of file HLTTrackMETProducer.h.

Referenced by produce().

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

edm::InputTag HLTTrackMETProducer::pfCandidatesLabel_

private

Definition at line 86 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer(), and produce().

edm::InputTag HLTTrackMETProducer::pfRecTracksLabel_

private

Definition at line 85 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer().

edm::InputTag HLTTrackMETProducer::tracksLabel_

private

Definition at line 84 of file HLTTrackMETProducer.h.

Referenced by HLTTrackMETProducer().

bool HLTTrackMETProducer::usePFCandidates_

private

Use PF candidates as input instead of jets.

Definition at line 66 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 63 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 59 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 51 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 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

Member Function Documentation

Member Data Documentation