HLTCaloTowerHtMhtProducer.cc

Go to the documentation of this file.

#include "HLTrigger/JetMET/interface/HLTCaloTowerHtMhtProducer.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "DataFormats/Common/interface/Handle.h"


// Constructor
HLTCaloTowerHtMhtProducer::HLTCaloTowerHtMhtProducer(const edm::ParameterSet & iConfig) :
  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>();
}

// Destructor
HLTCaloTowerHtMhtProducer::~HLTCaloTowerHtMhtProducer() = default;

// Fill descriptions
void HLTCaloTowerHtMhtProducer::fillDescriptions(edm::ConfigurationDescriptions & descriptions) {
    // 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);
}

// Produce the products
void HLTCaloTowerHtMhtProducer::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {

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