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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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...

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

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

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

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

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, and metLabel_.

       : usePt_         (iConfig.getParameter<bool>("usePt")),
         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 ( )

Definition at line 30 of file HLTMETCleanerUsingJetID.cc.

30 {}

Member Function Documentation

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

static

Definition at line 33 of file HLTMETCleanerUsingJetID.cc.

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

                                                                                          {
     edm::ParameterSetDescription desc;
     desc.add<bool>("usePt", false);
     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

Implements edm::stream::EDProducerBase.

Definition at line 45 of file HLTMETCleanerUsingJetID.cc.

References funct::abs(), funct::cos(), eta(), edm::Event::getByToken(), topObjectSelection_cff::goodJets, j, fwrapper::jets, m_theGoodJetToken, m_theJetToken, m_theMETToken, maxEta_, caloMETBenchmarkGeneric_cfi::met, minPt_, reco::LeafCandidate::p4(), phi, EnergyCorrector::pt, edm::Event::put(), query::result, reco::LeafCandidate::setP4(), funct::sin(), mathSSE::sqrt(), and usePt_.

                                                                                    {
 
     // Create a pointer to the products
     std::auto_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(reco::CaloJetCollection::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 > 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(reco::CaloJetCollection::const_iterator j = goodJets->begin(); j != goodJets->end(); ++j) {
             double pt = usePt_ ? j->pt() : j->et();
             double eta = j->eta();
             double phi = j->phi();
             double px = usePt_ ? j->px() : j->pt() * cos(phi);
             double py = usePt_ ? j->py() : j->pt() * sin(phi);
 
             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_diff(mex_diff, mey_diff, 0, sqrt(mex_diff*mex_diff + mey_diff*mey_diff));
 
         reco::CaloMET cleanmet = met->front();
         cleanmet.setP4(cleanmet.p4() + p4_diff);
         result->push_back(cleanmet);
     }
 
     iEvent.put( result );
 }

Member Data Documentation

edm::InputTag HLTMETCleanerUsingJetID::goodJetsLabel_

private

Input tag for the 'good jets' collection.

Definition at line 58 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

edm::InputTag HLTMETCleanerUsingJetID::jetsLabel_

private

Input tag for the 'all jets' collection.

Definition at line 55 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

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

private

Definition at line 62 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

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

private

Definition at line 61 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

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

private

Definition at line 60 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID(), and produce().

double HLTMETCleanerUsingJetID::maxEta_

private

Maximum (abs) eta requirement for jets.

Definition at line 49 of file HLTMETCleanerUsingJetID.h.

Referenced by produce().

edm::InputTag HLTMETCleanerUsingJetID::metLabel_

private

Input tag for the MET collection.

Definition at line 52 of file HLTMETCleanerUsingJetID.h.

Referenced by HLTMETCleanerUsingJetID().

double HLTMETCleanerUsingJetID::minPt_

private

Minimum pt requirement for jets.

Definition at line 46 of file HLTMETCleanerUsingJetID.h.

Referenced by produce().

bool HLTMETCleanerUsingJetID::usePt_

private

Use pt; otherwise, use et.

Definition at line 43 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