HLTHtMhtProducer Class Reference

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)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	~HLTHtMhtProducer () 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, 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

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 34 of file HLTHtMhtProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 17 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 ( )

overridedefault

Referenced by HLTHtMhtProducer().

Member Function Documentation

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

static

Definition at line 39 of file HLTHtMhtProducer.cc.

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

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

override

Definition at line 56 of file HLTHtMhtProducer.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, maxEtaJetHt_, maxEtaJetMht_, minNJetHt_, minNJetMht_, minPtJetHt_, minPtJetMht_, 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(), usePt_, 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;

  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->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 > 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 47 of file HLTHtMhtProducer.h.

Referenced by produce().

edm::InputTag HLTHtMhtProducer::jetsLabel_

private

Input jet, PFCandidate collections.

Definition at line 62 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer().

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

private

Definition at line 65 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer(), and produce().

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

private

Definition at line 66 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer(), and produce().

double HLTHtMhtProducer::maxEtaJetHt_

private

Maximum (abs) eta requirement for jets.

Definition at line 58 of file HLTHtMhtProducer.h.

Referenced by produce().

double HLTHtMhtProducer::maxEtaJetMht_

private

Definition at line 59 of file HLTHtMhtProducer.h.

Referenced by produce().

int HLTHtMhtProducer::minNJetHt_

private

Minimum number of jets passing pt and eta requirements.

Definition at line 50 of file HLTHtMhtProducer.h.

Referenced by produce().

int HLTHtMhtProducer::minNJetMht_

private

Definition at line 51 of file HLTHtMhtProducer.h.

Referenced by produce().

double HLTHtMhtProducer::minPtJetHt_

private

Minimum pt requirement for jets.

Definition at line 54 of file HLTHtMhtProducer.h.

Referenced by produce().

double HLTHtMhtProducer::minPtJetMht_

private

Definition at line 55 of file HLTHtMhtProducer.h.

Referenced by produce().

edm::InputTag HLTHtMhtProducer::pfCandidatesLabel_

private

Definition at line 63 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer(), and produce().

bool HLTHtMhtProducer::usePt_

private

Use pt; otherwise, use et.

Definition at line 43 of file HLTHtMhtProducer.h.

Referenced by produce().