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)

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

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

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, towersLabel_, and ~HLTCaloTowerHtMhtProducer().

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

overridedefault

Referenced by HLTCaloTowerHtMhtProducer().

Member Function Documentation

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

static

Definition at line 36 of file HLTCaloTowerHtMhtProducer.cc.

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

                                                                                             {
     // 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
	)

override

Definition at line 49 of file HLTCaloTowerHtMhtProducer.cc.

References funct::abs(), funct::cos(), edm::SortedCollection< T, SORT >::empty(), PVValHelper::eta, edm::Event::getByToken(), m_theTowersToken, maxEtaTowerHt_, maxEtaTowerMht_, minPtTowerHt_, minPtTowerMht_, eostools::move(), p4, phi, EnergyCorrector::pt, edm::Event::put(), mps_fire::result, funct::sin(), mathSSE::sqrt(), usePt_, and extraflags_cff::vtx.

                                                                                      {
 
     // Create a pointer to the products
     std::unique_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->empty()) {
         for(auto const & j : *towers) {
             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(std::move(result));
 }