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	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () 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 33 of file HLTMETCleanerUsingJetID.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 14 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 32 of file HLTMETCleanerUsingJetID.cc.

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

                                                                                         {
    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 43 of file HLTMETCleanerUsingJetID.cc.

References funct::abs(), PVValHelper::eta, edm::Event::getByToken(), topObjectSelection_cff::goodJets, fwrapper::jets, m_theGoodJetToken, m_theJetToken, m_theMETToken, maxEta_, RazorAnalyzer::met, minPt_, eostools::move(), EnergyCorrector::pt, edm::Event::put(), 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 56 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

edm::InputTag HLTMETCleanerUsingJetID::jetsLabel_

private

Input tag for the 'all jets' collection.

Definition at line 53 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

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

private

Definition at line 60 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

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

private

Definition at line 59 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

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

private

Definition at line 58 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

double HLTMETCleanerUsingJetID::maxEta_

private

Maximum (abs) eta requirement for jets.

Definition at line 47 of file HLTMETCleanerUsingJetID.h.

Referenced by produce().

edm::InputTag HLTMETCleanerUsingJetID::metLabel_

private

Input tag for the MET collection.

Definition at line 50 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

double HLTMETCleanerUsingJetID::minPt_

private

Minimum pt requirement for jets.

Definition at line 44 of file HLTMETCleanerUsingJetID.h.

Referenced by produce().