This produces a reco::MET object that stores HT and MHT. More...

#include <HLTCaloTowerHtMhtProducer.h>

Inheritance diagram for HLTCaloTowerHtMhtProducer:

Public Member Functions
	HLTCaloTowerHtMhtProducer (const edm::ParameterSet &iConfig)

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

	~HLTCaloTowerHtMhtProducer ()

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)

	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 (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &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::EDGetTokenT < CaloTowerCollection >	m_theTowersToken

double	maxEtaTowerHt_
	Maximum (abs) eta requirement for jets. More...

double	maxEtaTowerMht_

double	minPtTowerHt_
	Minimum pt requirement for jets. More...

double	minPtTowerMht_

edm::InputTag	towersLabel_
	Input CaloTower collection. 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

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 produces a reco::MET object that stores HT and MHT.

Author: Steven Lowette; Michele de Gruttola, Jia Fu Low (Nov 2013); Thiago Tomei (added the "HT from CaloTower code")

HT & MHT are calculated using input CaloTower collection. HT is stored as sumet_, MHT is stored as p4_ in the output.

See header file for documentation

Author: Steven Lowette; Thiago Tomei

Definition at line 31 of file HLTCaloTowerHtMhtProducer.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 19 of file HLTCaloTowerHtMhtProducer.cc.

References m_theTowersToken, and towersLabel_.

                                                                                     :
   usePt_                  ( iConfig.getParameter<bool>("usePt") ),
   minPtTowerHt_             ( iConfig.getParameter<double>("minPtTowerHt") ),
   minPtTowerMht_            ( iConfig.getParameter<double>("minPtTowerMht") ),
   maxEtaTowerHt_            ( iConfig.getParameter<double>("maxEtaTowerHt") ),
   maxEtaTowerMht_           ( iConfig.getParameter<double>("maxEtaTowerMht") ),
   towersLabel_              ( iConfig.getParameter<edm::InputTag>("towersLabel") ) {
     m_theTowersToken = consumes<CaloTowerCollection>(towersLabel_);
     
     // Register the products
     produces<reco::METCollection>();
 }

HLTCaloTowerHtMhtProducer::~HLTCaloTowerHtMhtProducer ( )

Definition at line 33 of file HLTCaloTowerHtMhtProducer.cc.

33 {}

Member Function Documentation

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

static

Definition at line 36 of file HLTCaloTowerHtMhtProducer.cc.

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

                                                                                             {
     // Current default is for hltHtMht
     edm::ParameterSetDescription desc;
     desc.add<bool>("usePt", false);
     desc.add<double>("minPtTowerHt", 1.);
     desc.add<double>("minPtTowerMht", 1.);
     desc.add<double>("maxEtaTowerHt", 5.);
     desc.add<double>("maxEtaTowerMht", 5.);
     desc.add<edm::InputTag>("towersLabel", edm::InputTag("hltTowerMakerForAll"));
     descriptions.add("hltCaloTowerHtMhtProducer", desc);
 }

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

virtual

Implements edm::stream::EDProducerBase.

Definition at line 49 of file HLTCaloTowerHtMhtProducer.cc.

References funct::abs(), funct::cos(), eta, edm::Event::getByToken(), j, m_theTowersToken, maxEtaTowerHt_, maxEtaTowerMht_, minPtTowerHt_, minPtTowerMht_, p4, phi, EnergyCorrector::pt, edm::Event::put(), mps_fire::result, funct::sin(), mathSSE::sqrt(), HLT_FULL_cff::towers, and usePt_.

                                                                                      {
 
     // Create a pointer to the products
     std::auto_ptr<reco::METCollection> result(new reco::METCollection());
 
     edm::Handle<CaloTowerCollection> towers;
     iEvent.getByToken(m_theTowersToken, towers);
 
     double ht = 0., mhx = 0., mhy = 0.;
 
     if (towers->size() > 0) {
         for(CaloTowerCollection::const_iterator j = towers->begin(); j != towers->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 > minPtTowerHt_ && std::abs(eta) < maxEtaTowerHt_) {
                 ht += pt;
             }
 
             if (pt > minPtTowerMht_ && std::abs(eta) < maxEtaTowerMht_) {
                 mhx -= px;
                 mhy -= py;
             }
         }
     }
 
     reco::MET::LorentzVector p4(mhx, mhy, 0, sqrt(mhx*mhx + mhy*mhy));
     reco::MET::Point vtx(0, 0, 0);
     reco::MET htmht(ht, p4, vtx);
     result->push_back(htmht);
     
     // Put the products into the Event
     iEvent.put(result);
 }