L1HLTTauMatching.cc

Go to the documentation of this file.

#include "RecoTauTag/HLTProducers/interface/L1HLTTauMatching.h"
#include "Math/GenVector/VectorUtil.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/L1Trigger/interface/L1JetParticle.h"
#include "DataFormats/L1Trigger/interface/L1JetParticleFwd.h"
#include "DataFormats/HLTReco/interface/TriggerTypeDefs.h"
#include "FWCore/Utilities/interface/InputTag.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "DataFormats/TauReco/interface/PFTau.h"

//
// class decleration
//
using namespace reco;
using namespace std;
using namespace edm;
using namespace l1extra;

L1HLTTauMatching::L1HLTTauMatching(const edm::ParameterSet& iConfig)
    : jetSrc(consumes<PFTauCollection>(iConfig.getParameter<InputTag>("JetSrc"))),
      tauTrigger(consumes<trigger::TriggerFilterObjectWithRefs>(iConfig.getParameter<InputTag>("L1TauTrigger"))),
      mEt_Min(iConfig.getParameter<double>("EtMin")) {
  produces<PFTauCollection>();
}

void L1HLTTauMatching::produce(edm::StreamID iSId, edm::Event& iEvent, const edm::EventSetup& iES) const {
  using namespace edm;
  using namespace std;
  using namespace reco;
  using namespace trigger;
  using namespace l1extra;

  unique_ptr<PFTauCollection> tauL2jets(new PFTauCollection);

  constexpr double matchingR2 = 0.5 * 0.5;

  //Getting HLT jets to be matched
  edm::Handle<PFTauCollection> tauJets;
  iEvent.getByToken(jetSrc, tauJets);

  edm::Handle<trigger::TriggerFilterObjectWithRefs> l1TriggeredTaus;
  iEvent.getByToken(tauTrigger, l1TriggeredTaus);

  vector<l1extra::L1JetParticleRef> tauCandRefVec;
  vector<l1extra::L1JetParticleRef> jetCandRefVec;
  l1TriggeredTaus->getObjects(trigger::TriggerL1TauJet, tauCandRefVec);
  l1TriggeredTaus->getObjects(trigger::TriggerL1CenJet, jetCandRefVec);

  math::XYZPoint a(0., 0., 0.);

  for (unsigned int iL1Tau = 0; iL1Tau < tauCandRefVec.size(); iL1Tau++) {
    for (unsigned int iJet = 0; iJet < tauJets->size(); iJet++) {
      //Find the relative L2TauJets, to see if it has been reconstructed
      const PFTau& myJet = (*tauJets)[iJet];
      double deltaR2 = ROOT::Math::VectorUtil::DeltaR2(myJet.p4().Vect(), (tauCandRefVec[iL1Tau]->p4()).Vect());
      if (deltaR2 < matchingR2) {
        //       LeafCandidate myLC(myJet);
        if (myJet.leadChargedHadrCand().isNonnull()) {
          a = myJet.leadChargedHadrCand()->vertex();
        }
        PFTau myPFTau(std::numeric_limits<int>::quiet_NaN(), myJet.p4(), a);
        if (myJet.pt() > mEt_Min) {
          //          tauL2LC->push_back(myLC);
          tauL2jets->push_back(myPFTau);
        }
        break;
      }
    }
  }

  for (unsigned int iL1Tau = 0; iL1Tau < jetCandRefVec.size(); iL1Tau++) {
    for (unsigned int iJet = 0; iJet < tauJets->size(); iJet++) {
      const PFTau& myJet = (*tauJets)[iJet];
      //Find the relative L2TauJets, to see if it has been reconstructed
      double deltaR2 = ROOT::Math::VectorUtil::DeltaR2(myJet.p4().Vect(), (jetCandRefVec[iL1Tau]->p4()).Vect());
      if (deltaR2 < matchingR2) {
        //       LeafCandidate myLC(myJet);
        if (myJet.leadChargedHadrCand().isNonnull()) {
          a = myJet.leadChargedHadrCand()->vertex();
        }

        PFTau myPFTau(std::numeric_limits<int>::quiet_NaN(), myJet.p4(), a);
        if (myJet.pt() > mEt_Min) {
          //          tauL2LC->push_back(myLC);
          tauL2jets->push_back(myPFTau);
        }
        break;
      }
    }
  }

  iEvent.put(std::move(tauL2jets));
  // iEvent.put(std::move(tauL2LC));
}

void L1HLTTauMatching::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("L1TauTrigger", edm::InputTag("hltL1sDoubleIsoTau40er"))
      ->setComment("Name of trigger filter");
  desc.add<edm::InputTag>("JetSrc", edm::InputTag("hltSelectedPFTausTrackPt1MediumIsolationReg"))
      ->setComment("Input collection of PFTaus");
  desc.add<double>("EtMin", 0.0)->setComment("Minimal pT of PFTau to match");
  descriptions.setComment(
      "This module produces collection of PFTaus matched to L1ExtraTaus/Jets passing a HLT filter (Only p4 and vertex "
      "of returned PFTaus are set).");
  descriptions.add("L1HLTJetsMatching", desc);
}

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(L1HLTTauMatching);