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
	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, PFCandidate collections. More...
edm::EDGetTokenT < reco::CandidateView >	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<>
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 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 31 of file HLTHtMhtProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 20 of file HLTHtMhtProducer.cc.

References jetsLabel_, m_theJetToken, m_thePFCandidateToken, and pfCandidatesLabel_.

    : 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<reco::CandidateView>(jetsLabel_);
  m_thePFCandidateToken = consumes<reco::PFCandidateCollection>(pfCandidatesLabel_);

  // Register the products
  produces<reco::METCollection>();
}

HLTHtMhtProducer::~HLTHtMhtProducer ( )

overridedefault

Member Function Documentation

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

static

Definition at line 42 of file HLTHtMhtProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, and HLT_FULL_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 59 of file HLTHtMhtProducer.cc.

References funct::abs(), funct::cos(), PVValHelper::eta, excludePFMuons_, edm::Event::getByToken(), fwrapper::jets, edm::InputTag::label(), m_theJetToken, m_thePFCandidateToken, maxEtaJetHt_, maxEtaJetMht_, minNJetHt_, minNJetMht_, minPtJetHt_, minPtJetMht_, eostools::move(), slimmedMuons_cfi::pfCandidates, pfCandidatesLabel_, phi, DiDispStaMuonMonitor_cfi::pt, edm::Event::put(), mps_fire::result, funct::sin(), mathSSE::sqrt(), and usePt_.

                                                                            {
  // Create a pointer to the products
  std::unique_ptr<reco::METCollection> result(new reco::METCollection());

  if (pfCandidatesLabel_.label().empty())
    excludePFMuons_ = false;

  edm::Handle<reco::CandidateView> 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.;

  for (auto const& aJet : *jets) {
    double const pt = usePt_ ? aJet.pt() : aJet.et();
    double const eta = aJet.eta();
    double const phi = aJet.phi();
    double const px = usePt_ ? aJet.px() : aJet.et() * cos(phi);
    double const py = usePt_ ? aJet.py() : aJet.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& aCand : *pfCandidates) {
      if (std::abs(aCand.pdgId()) == 13) {
        mhx += aCand.px();
        mhy += aCand.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 44 of file HLTHtMhtProducer.h.

Referenced by produce().

edm::InputTag HLTHtMhtProducer::jetsLabel_

private

Input jet, PFCandidate collections.

Definition at line 59 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer().

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

private

Definition at line 62 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer(), and produce().

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

private

Definition at line 63 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer(), and produce().

double HLTHtMhtProducer::maxEtaJetHt_

private

Maximum (abs) eta requirement for jets.

Definition at line 55 of file HLTHtMhtProducer.h.

Referenced by produce().

double HLTHtMhtProducer::maxEtaJetMht_

private

Definition at line 56 of file HLTHtMhtProducer.h.

Referenced by produce().

int HLTHtMhtProducer::minNJetHt_

private

Minimum number of jets passing pt and eta requirements.

Definition at line 47 of file HLTHtMhtProducer.h.

Referenced by produce().

int HLTHtMhtProducer::minNJetMht_

private

Definition at line 48 of file HLTHtMhtProducer.h.

Referenced by produce().

double HLTHtMhtProducer::minPtJetHt_

private

Minimum pt requirement for jets.

Definition at line 51 of file HLTHtMhtProducer.h.

Referenced by produce().

double HLTHtMhtProducer::minPtJetMht_

private

Definition at line 52 of file HLTHtMhtProducer.h.

Referenced by produce().

edm::InputTag HLTHtMhtProducer::pfCandidatesLabel_

private

Definition at line 60 of file HLTHtMhtProducer.h.

Referenced by HLTHtMhtProducer(), and produce().

bool HLTHtMhtProducer::usePt_

private

Use pt; otherwise, use et.

Definition at line 40 of file HLTHtMhtProducer.h.

Referenced by produce().