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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T...>

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T...>

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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 18 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 ( )

overridedefault

Member Function Documentation

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

static

Definition at line 35 of file HLTCaloTowerHtMhtProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), submitPVResolutionJobs::desc, 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
	)

override

Definition at line 48 of file HLTCaloTowerHtMhtProducer.cc.

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

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