HLTMETCleanerUsingJetID Class Reference

This creates a MET object from the difference in MET between two input jet collections. More...

#include <HLTMETCleanerUsingJetID.h>

Inheritance diagram for HLTMETCleanerUsingJetID:

Public Member Functions
	HLTMETCleanerUsingJetID (const edm::ParameterSet &iConfig)
void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override
	~HLTMETCleanerUsingJetID () 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
edm::InputTag	goodJetsLabel_
	Input tag for the 'good jets' collection. More...
edm::InputTag	jetsLabel_
	Input tag for the 'all jets' collection. More...
edm::EDGetTokenT< reco::CaloJetCollection >	m_theGoodJetToken
edm::EDGetTokenT< reco::CaloJetCollection >	m_theJetToken
edm::EDGetTokenT< reco::CaloMETCollection >	m_theMETToken
double	maxEta_
	Maximum (abs) eta requirement for jets. More...
edm::InputTag	metLabel_
	Input tag for the MET collection. More...
double	minPt_
	Minimum pt requirement for jets. 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 creates a MET object from the difference in MET between two input jet collections.

Author

Jia Fu Low (Nov 2013)

This code creates a new MET vector defined as:

output MET = input MET + MET from 'good jets' - MET from 'all jets'

See header file for more information.

Author

a Jet/MET person

Definition at line 32 of file HLTMETCleanerUsingJetID.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file HLTMETCleanerUsingJetID.cc.

References goodJetsLabel_, jetsLabel_, m_theGoodJetToken, m_theJetToken, m_theMETToken, metLabel_, and ~HLTMETCleanerUsingJetID().

    : minPt_(iConfig.getParameter<double>("minPt")),
      maxEta_(iConfig.getParameter<double>("maxEta")),
      metLabel_(iConfig.getParameter<edm::InputTag>("metLabel")),
      jetsLabel_(iConfig.getParameter<edm::InputTag>("jetsLabel")),
      goodJetsLabel_(iConfig.getParameter<edm::InputTag>("goodJetsLabel")) {
  m_theMETToken = consumes<reco::CaloMETCollection>(metLabel_);
  m_theJetToken = consumes<reco::CaloJetCollection>(jetsLabel_);
  m_theGoodJetToken = consumes<reco::CaloJetCollection>(goodJetsLabel_);

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

HLTMETCleanerUsingJetID::~HLTMETCleanerUsingJetID ( )

overridedefault

Referenced by HLTMETCleanerUsingJetID().

Member Function Documentation

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

static

Definition at line 31 of file HLTMETCleanerUsingJetID.cc.

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

                                                                                         {
  edm::ParameterSetDescription desc;
  desc.add<double>("minPt", 20.);
  desc.add<double>("maxEta", 5.);
  desc.add<edm::InputTag>("metLabel", edm::InputTag("hltMet"));
  desc.add<edm::InputTag>("jetsLabel", edm::InputTag("hltAntiKT4CaloJets"));
  desc.add<edm::InputTag>("goodJetsLabel", edm::InputTag("hltCaloJetIDPassed"));
  descriptions.add("hltMETCleanerUsingJetID", desc);
}

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

override

Definition at line 42 of file HLTMETCleanerUsingJetID.cc.

References funct::abs(), PVValHelper::eta, edm::Event::getByToken(), topObjectSelection_cff::goodJets, dqmiolumiharvest::j, singleTopDQM_cfi::jets, m_theGoodJetToken, m_theJetToken, m_theMETToken, maxEta_, BTaggingMonitor_cfi::met, minPt_, eostools::move(), DiDispStaMuonMonitor_cfi::pt, edm::Event::put(), multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, mps_fire::result, and mathSSE::sqrt().

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

  edm::Handle<reco::CaloMETCollection> met;
  edm::Handle<reco::CaloJetCollection> jets;
  edm::Handle<reco::CaloJetCollection> goodJets;

  iEvent.getByToken(m_theMETToken, met);
  iEvent.getByToken(m_theJetToken, jets);
  iEvent.getByToken(m_theGoodJetToken, goodJets);

  double mex_jets = 0.;
  double mey_jets = 0.;
  double sumet_jets = 0.;
  if (!jets->empty()) {
    for (auto const& j : *jets) {
      double pt = j.pt();
      double eta = j.eta();
      double px = j.px();
      double py = j.py();

      if (pt > minPt_ && std::abs(eta) < maxEta_) {
        mex_jets -= px;
        mey_jets -= py;
        sumet_jets += pt;
      }
    }
  }

  double mex_goodJets = 0.;
  double mey_goodJets = 0.;
  double sumet_goodJets = 0.;
  if (!goodJets->empty()) {
    for (auto const& j : *goodJets) {
      double pt = j.pt();
      double eta = j.eta();
      double px = j.px();
      double py = j.py();

      if (pt > minPt_ && std::abs(eta) < maxEta_) {
        mex_goodJets -= px;
        mey_goodJets -= py;
        sumet_goodJets += pt;
      }
    }
  }

  if (!met->empty()) {
    double mex_diff = mex_goodJets - mex_jets;
    double mey_diff = mey_goodJets - mey_jets;
    //double sumet_diff = sumet_goodJets - sumet_jets;  // cannot set sumet...
    reco::Candidate::LorentzVector p4_clean(met->front().px() + mex_diff,
                                            mey_diff + met->front().py(),
                                            0,
                                            sqrt((met->front().px() + mex_diff) * (met->front().px() + mex_diff) +
                                                 (met->front().py() + mey_diff) * (met->front().py() + mey_diff)));

    reco::CaloMET cleanmet = met->front();
    cleanmet.setP4(p4_clean);
    result->push_back(cleanmet);
  }

  iEvent.put(std::move(result));
}

Member Data Documentation

edm::InputTag HLTMETCleanerUsingJetID::goodJetsLabel_

private

Input tag for the 'good jets' collection.

Definition at line 54 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

edm::InputTag HLTMETCleanerUsingJetID::jetsLabel_

private

Input tag for the 'all jets' collection.

Definition at line 51 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

edm::EDGetTokenT<reco::CaloJetCollection> HLTMETCleanerUsingJetID::m_theGoodJetToken

private

Definition at line 58 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

edm::EDGetTokenT<reco::CaloJetCollection> HLTMETCleanerUsingJetID::m_theJetToken

private

Definition at line 57 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

edm::EDGetTokenT<reco::CaloMETCollection> HLTMETCleanerUsingJetID::m_theMETToken

private

Definition at line 56 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

double HLTMETCleanerUsingJetID::maxEta_

private

Maximum (abs) eta requirement for jets.

Definition at line 45 of file HLTMETCleanerUsingJetID.h.

Referenced by produce().

edm::InputTag HLTMETCleanerUsingJetID::metLabel_

private

Input tag for the MET collection.

Definition at line 48 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

double HLTMETCleanerUsingJetID::minPt_

private

Minimum pt requirement for jets.

Definition at line 42 of file HLTMETCleanerUsingJetID.h.

Referenced by produce().