CaloTowerFromL1TSeededCreatorForTauHLT.cc

Go to the documentation of this file.

// makes CaloTowerCandidates from CaloTowers
// original author: L.Lista INFN, modifyed by: F.Ratnikov UMd
// Author for regionality A. Nikitenko
// Modified by S. Gennai + T. Strebler

#include "DataFormats/RecoCandidate/interface/RecoCaloTowerCandidate.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Math/interface/deltaR.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "RecoTauTag/HLTProducers/interface/CaloTowerFromL1TSeededCreatorForTauHLT.h"
// Math
#include "Math/GenVector/VectorUtil.h"
#include <cmath>

using namespace edm;
using namespace reco;
using namespace std;

CaloTowerFromL1TSeededCreatorForTauHLT::CaloTowerFromL1TSeededCreatorForTauHLT(const ParameterSet& p)
    : m_verbose(p.getUntrackedParameter<int>("verbose", 0)),
      m_towers_token(consumes<CaloTowerCollection>(p.getParameter<InputTag>("towers"))),
      m_cone(p.getParameter<double>("UseTowersInCone")),
      m_tauTrigger_token(consumes<trigger::TriggerFilterObjectWithRefs>(p.getParameter<InputTag>("TauTrigger"))),
      m_EtThreshold(p.getParameter<double>("minimumEt")),
      m_EThreshold(p.getParameter<double>("minimumE")) {
  produces<CaloTowerCollection>();
}

CaloTowerFromL1TSeededCreatorForTauHLT::~CaloTowerFromL1TSeededCreatorForTauHLT() = default;

void CaloTowerFromL1TSeededCreatorForTauHLT::produce(StreamID sid, Event& evt, const EventSetup& stp) const {
  edm::Handle<CaloTowerCollection> caloTowers;
  evt.getByToken(m_towers_token, caloTowers);

  double m_cone2 = m_cone * m_cone;

  // L1 seeds
  edm::Handle<trigger::TriggerFilterObjectWithRefs> l1TriggeredTaus;
  evt.getByToken(m_tauTrigger_token, l1TriggeredTaus);

  auto cands = std::make_unique<CaloTowerCollection>();
  cands->reserve(caloTowers->size());

  l1t::TauVectorRef tauCandRefVec;
  l1TriggeredTaus->getObjects(trigger::TriggerL1Tau, tauCandRefVec);

  for (auto const& tauCandRef : tauCandRefVec) {
    for (auto const& cal : *caloTowers) {
      bool isAccepted = false;
      if (m_verbose == 2) {
        edm::LogInfo("JetDebugInfo") << "CaloTowerFromL1TSeededCreatorForTauHLT::produce->  tower et/eta/phi/e: "
                                     << cal.et() << '/' << cal.eta() << '/' << cal.phi() << '/' << cal.energy()
                                     << " is...";
      }
      if (cal.et() >= m_EtThreshold && cal.energy() >= m_EThreshold) {
        math::PtEtaPhiELorentzVector p(cal.et(), cal.eta(), cal.phi(), cal.energy());
        double delta2 = deltaR2((tauCandRef->p4()).Vect(), p);
        if (delta2 < m_cone2) {
          isAccepted = true;
          cands->push_back(cal);
        }
      }

      if (m_verbose == 2) {
        if (isAccepted)
          edm::LogInfo("JetDebugInfo") << "accepted \n";
        else
          edm::LogInfo("JetDebugInfo") << "rejected \n";
      }
    }
  }

  evt.put(std::move(cands));
}

void CaloTowerFromL1TSeededCreatorForTauHLT::fillDescriptions(edm::ConfigurationDescriptions& desc) {
  edm::ParameterSetDescription aDesc;

  aDesc.add<edm::InputTag>("TauTrigger", edm::InputTag("hltL1sDoubleIsoTau40er"))
      ->setComment("Name of trigger filter for L1 seeds");
  aDesc.add<edm::InputTag>("towers", edm::InputTag("towerMaker"))->setComment("Input tower collection");
  aDesc.add<double>("UseTowersInCone", 0.8)->setComment("Radius of cone around seed");
  aDesc.add<double>("minimumE", 0.8)->setComment("Minimum tower energy");
  aDesc.add<double>("minimumEt", 0.5)->setComment("Minimum tower ET");
  aDesc.addUntracked<int>("verbose", 0)->setComment("Verbosity level; 0=silent");

  desc.add("CaloTowerFromL1TSeededCreatorForTauHLT", aDesc);
  desc.setComment("Produce tower collection around L1 particle seed.");
}