HLTJetL1TMatchProducer.cc

Go to the documentation of this file.

#include <string>

#include "HLTrigger/JetMET/interface/HLTJetL1TMatchProducer.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/Math/interface/deltaPhi.h"

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "HLTrigger/HLTcore/interface/defaultModuleLabel.h"

template <typename T>
HLTJetL1TMatchProducer<T>::HLTJetL1TMatchProducer(const edm::ParameterSet& iConfig) {
  jetsInput_ = iConfig.template getParameter<edm::InputTag>("jetsInput");
  L1Jets_ = iConfig.template getParameter<edm::InputTag>("L1Jets");
  DeltaR_ = iConfig.template getParameter<double>("DeltaR");

  typedef std::vector<T> TCollection;
  m_theJetToken = consumes<TCollection>(jetsInput_);
  m_theL1JetToken = consumes<l1t::JetBxCollection>(L1Jets_);
  produces<TCollection>();
}

template <typename T>
void HLTJetL1TMatchProducer<T>::beginJob() {}

template <typename T>
HLTJetL1TMatchProducer<T>::~HLTJetL1TMatchProducer() = default;

template <typename T>
void HLTJetL1TMatchProducer<T>::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("jetsInput", edm::InputTag("hltAntiKT5PFJets"));
  desc.add<edm::InputTag>("L1Jets", edm::InputTag("hltCaloStage2Digis"));
  desc.add<double>("DeltaR", 0.5);
  descriptions.add(defaultModuleLabel<HLTJetL1TMatchProducer<T>>(), desc);
}

template <typename T>
void HLTJetL1TMatchProducer<T>::produce(edm::Event& iEvent, const edm::EventSetup& iSetup) {
  typedef std::vector<T> TCollection;

  edm::Handle<TCollection> jets;
  iEvent.getByToken(m_theJetToken, jets);

  std::unique_ptr<TCollection> result(new TCollection);

  edm::Handle<l1t::JetBxCollection> l1Jets;
  iEvent.getByToken(m_theL1JetToken, l1Jets);
  bool trigger_bx_only = true;  // selection of BX not implemented

  if (l1Jets.isValid()) {
    typename TCollection::const_iterator jet_iter;
    for (jet_iter = jets->begin(); jet_iter != jets->end(); ++jet_iter) {
      bool isMatched = false;
      for (int ibx = l1Jets->getFirstBX(); ibx <= l1Jets->getLastBX(); ++ibx) {
        if (trigger_bx_only && (ibx != 0))
          continue;
        for (auto it = l1Jets->begin(ibx); it != l1Jets->end(ibx); it++) {
          if (it->et() == 0)
            continue;  // if you don't care about L1T candidates with zero ET.
          const double deltaeta = jet_iter->eta() - it->eta();
          const double deltaphi = deltaPhi(jet_iter->phi(), it->phi());
          if (sqrt(deltaeta * deltaeta + deltaphi * deltaphi) < DeltaR_)
            isMatched = true;
          //cout << "bx:  " << ibx << "  et:  "  << it->et() << "  eta:  "  << it->eta() << "  phi:  "  << it->phi() << "\n";
        }
      }
      if (isMatched == true)
        result->push_back(*jet_iter);
    }  // jet_iter
  } else {
    edm::LogWarning("MissingProduct") << "L1Upgrade l1Jets bx collection not found." << std::endl;
  }

  iEvent.put(std::move(result));
}