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
	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
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<>
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 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::HLTMETCleanerUsingJetID ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 13 of file HLTMETCleanerUsingJetID.cc.

References goodJetsLabel_, jetsLabel_, m_theGoodJetToken, m_theJetToken, m_theMETToken, and metLabel_.

    : 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::~HLTMETCleanerUsingJetID ( )

overridedefault

Member Function Documentation

fillDescriptions()

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

static

Definition at line 31 of file HLTMETCleanerUsingJetID.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and ProducerED_cfi::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);
}

produce()

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

override

Definition at line 42 of file HLTMETCleanerUsingJetID.cc.

References funct::abs(), PVValHelper::eta, topObjectSelection_cff::goodJets, iEvent, dqmiolumiharvest::j, PDWG_EXODelayedJetMET_cff::jets, m_theGoodJetToken, m_theJetToken, m_theMETToken, maxEta_, BTaggingMonitor_cfi::met, minPt_, eostools::move(), DiDispStaMuonMonitor_cfi::pt, 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

goodJetsLabel_

edm::InputTag HLTMETCleanerUsingJetID::goodJetsLabel_

private

Input tag for the 'good jets' collection.

Definition at line 54 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

jetsLabel_

edm::InputTag HLTMETCleanerUsingJetID::jetsLabel_

private

Input tag for the 'all jets' collection.

Definition at line 51 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

m_theGoodJetToken

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

private

Definition at line 58 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

m_theJetToken

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

private

Definition at line 57 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

m_theMETToken

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

private

Definition at line 56 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

maxEta_

double HLTMETCleanerUsingJetID::maxEta_

private

Maximum (abs) eta requirement for jets.

Definition at line 45 of file HLTMETCleanerUsingJetID.h.

Referenced by produce().

metLabel_

edm::InputTag HLTMETCleanerUsingJetID::metLabel_

private

Input tag for the MET collection.

Definition at line 48 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

minPt_

double HLTMETCleanerUsingJetID::minPt_

private

Minimum pt requirement for jets.

Definition at line 42 of file HLTMETCleanerUsingJetID.h.

Referenced by produce().