HLTMhtProducer Class Reference

This produces a reco::MET object that stores MHT (or MET) More...

#include <HLTMhtProducer.h>

Inheritance diagram for HLTMhtProducer:

Public Member Functions
	HLTMhtProducer (const edm::ParameterSet &iConfig)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	~HLTMhtProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () 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	maxEtaJet_
	Maximum (abs) eta requirement for jets. More...
int	minNJet_
	Minimum number of jets passing pt and eta requirements. More...
double	minPtJet_
	Minimum pt requirement for jets. More...
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 MHT (or MET)

Author

Steven Lowette

Michele de Gruttola, Jia Fu Low (Nov 2013)

MHT (or MET) is calculated using input CaloJet or PFJet collection. MHT can include or exclude the contribution from muons.

See header file for documentation

Author

Steven Lowette

Definition at line 33 of file HLTMhtProducer.h.

Constructor & Destructor Documentation

HLTMhtProducer()

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

explicit

Definition at line 17 of file HLTMhtProducer.cc.

    : usePt_(iConfig.getParameter<bool>("usePt")),
      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")),
      pfCandidatesLabel_(iConfig.getParameter<edm::InputTag>("pfCandidatesLabel")) {
  m_theJetToken = consumes<edm::View<reco::Jet>>(jetsLabel_);
  if (pfCandidatesLabel_.label().empty())
    excludePFMuons_ = false;
  if (excludePFMuons_)
    m_thePFCandidateToken = consumes<reco::PFCandidateCollection>(pfCandidatesLabel_);
 
  // Register the products
  produces<reco::METCollection>();
}

References excludePFMuons_, jetsLabel_, edm::InputTag::label(), m_theJetToken, m_thePFCandidateToken, and pfCandidatesLabel_.

~HLTMhtProducer()

HLTMhtProducer::~HLTMhtProducer ( )

overridedefault

Member Function Documentation

fillDescriptions()

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

static

Definition at line 39 of file HLTMhtProducer.cc.

                                                                                {
  // Current default is for hltPFMET
  edm::ParameterSetDescription desc;
  desc.add<bool>("usePt", true);
  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>("pfCandidatesLabel", edm::InputTag("hltParticleFlow"));
  descriptions.add("hltMhtProducer", desc);
}

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and HLT_FULL_cff::InputTag.

produce()

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

override

Definition at line 53 of file HLTMhtProducer.cc.

                                                                          {
  // Create a pointer to the products
  std::unique_ptr<reco::METCollection> result(new reco::METCollection());
 
  edm::Handle<reco::JetView> jets;
  iEvent.getByToken(m_theJetToken, jets);
 
  edm::Handle<reco::PFCandidateCollection> pfCandidates;
  if (excludePFMuons_)
    iEvent.getByToken(m_thePFCandidateToken, pfCandidates);
 
  int nj = 0;
  double sumet = 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 > 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, m_theJetToken, m_thePFCandidateToken, maxEtaJet_, minNJet_, minPtJet_, eostools::move(), p4, zmumugammaAnalyzer_cfi::pfCandidates, phi, DiDispStaMuonMonitor_cfi::pt, multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, mps_fire::result, funct::sin(), mathSSE::sqrt(), usePt_, and extraflags_cff::vtx.

Member Data Documentation

excludePFMuons_

bool HLTMhtProducer::excludePFMuons_

private

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

Definition at line 46 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer(), and produce().

jetsLabel_

edm::InputTag HLTMhtProducer::jetsLabel_

private

Input jet, PFCandidate collections.

Definition at line 58 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer().

m_theJetToken

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

private

Definition at line 61 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer(), and produce().

m_thePFCandidateToken

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

private

Definition at line 62 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer(), and produce().

maxEtaJet_

double HLTMhtProducer::maxEtaJet_

private

Maximum (abs) eta requirement for jets.

Definition at line 55 of file HLTMhtProducer.h.

Referenced by produce().

minNJet_

int HLTMhtProducer::minNJet_

private

Minimum number of jets passing pt and eta requirements.

Definition at line 49 of file HLTMhtProducer.h.

Referenced by produce().

minPtJet_

double HLTMhtProducer::minPtJet_

private

Minimum pt requirement for jets.

Definition at line 52 of file HLTMhtProducer.h.

Referenced by produce().

pfCandidatesLabel_

edm::InputTag HLTMhtProducer::pfCandidatesLabel_

private

Definition at line 59 of file HLTMhtProducer.h.

Referenced by HLTMhtProducer().

usePt_

bool HLTMhtProducer::usePt_

private

Use pt; otherwise, use et.

Definition at line 42 of file HLTMhtProducer.h.

Referenced by produce().