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	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () 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 18 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 35 of file HLTCaloTowerHtMhtProducer.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), and HLT_2018_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(), edm::SortedCollection< T, SORT >::empty(), PVValHelper::eta, edm::Event::getByToken(), dqmiolumiharvest::j, m_theTowersToken, maxEtaTowerHt_, maxEtaTowerMht_, minPtTowerHt_, minPtTowerMht_, eostools::move(), p4, phi, DiDispStaMuonMonitor_cfi::pt, edm::Event::put(), multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, mps_fire::result, funct::sin(), mathSSE::sqrt(), HLT_2018_cff::towers, usePt_, and badGlobalMuonTaggersAOD_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));
 }