HLTTrackMETProducer Class Reference

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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::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_, 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 ( )

overridedefault

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(), submitPVResolutionJobs::desc, 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  )

override

Definition at line 64 of file HLTTrackMETProducer.cc.

References funct::abs(), funct::cos(), PVValHelper::eta, excludePFMuons_, edm::Event::getByToken(), dqmiolumiharvest::j, fwrapper::jets, edm::InputTag::label(), m_theJetToken, m_thePFCandidateToken, m_theRecTrackToken, m_theTrackToken, maxEtaJet_, minNJet_, minPtJet_, eostools::move(), slimmedMuons_cfi::pfCandidates, pfCandidatesLabel_, phi, DiDispStaMuonMonitor_cfi::pt, edm::Event::put(), mps_fire::result, funct::sin(), mathSSE::sqrt(), 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().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().