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)

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

	~HLTMETCleanerUsingJetID ()

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (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

Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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 ( )

default

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
	)

virtual

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->size() > 0 ) {
         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->size() > 0) {
         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->size() > 0) {
         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().

Public Member Functions

Static Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation