SeedingRegionByL1.cc

Go to the documentation of this file.

/*
   Author: Swagata Mukherjee

   Date: Feb 2021

   TICL is currently seeded by tracks, or just globally.
   Here, adding option to seed TICL by L1 e/gamma objects (L1 TkEm).
   This is expected to be useful for CPU timing at the HLT.
*/

#include "SeedingRegionByL1.h"

#include <algorithm>
#include <set>
#include <vector>

#include "FWCore/MessageLogger/interface/MessageLogger.h"

ticl::SeedingRegionByL1::SeedingRegionByL1(const edm::ParameterSet &conf, edm::ConsumesCollector &sumes)
    : SeedingRegionAlgoBase(conf, sumes),
      l1TkEmsToken_(sumes.consumes<std::vector<l1t::TkEm>>(conf.getParameter<edm::InputTag>("l1TkEmColl"))),
      algoVerbosity_(conf.getParameter<int>("algo_verbosity")),
      minPt_(conf.getParameter<double>("minPt")),
      minAbsEta_(conf.getParameter<double>("minAbsEta")),
      maxAbsEta_(conf.getParameter<double>("maxAbsEta")),
      endcapScalings_(conf.getParameter<std::vector<double>>("endcapScalings")),
      quality_(conf.getParameter<int>("quality")) {}

void ticl::SeedingRegionByL1::makeRegions(const edm::Event &ev,
                                          const edm::EventSetup &es,
                                          std::vector<TICLSeedingRegion> &result) {
  auto l1TrkEms = ev.getHandle(l1TkEmsToken_);
  edm::ProductID l1tkemsId = l1TrkEms.id();

  for (size_t indx = 0; indx < (*l1TrkEms).size(); indx++) {
    const auto &l1TrkEm = (*l1TrkEms)[indx];
    double offlinePt = this->tkEmOfflineEt(l1TrkEm.pt());
    if ((offlinePt < minPt_) || (std::abs(l1TrkEm.eta()) < minAbsEta_) || (std::abs(l1TrkEm.eta()) > maxAbsEta_) ||
        (l1TrkEm.EGRef()->hwQual() != quality_)) {
      continue;
    }

    int iSide = int(l1TrkEm.eta() > 0);
    result.emplace_back(GlobalPoint(l1TrkEm.p4().X(), l1TrkEm.p4().Y(), l1TrkEm.p4().Z()),
                        GlobalVector(l1TrkEm.px(), l1TrkEm.py(), l1TrkEm.pz()),
                        iSide,
                        indx,
                        l1tkemsId);
  }

  std::sort(result.begin(), result.end(), [](const TICLSeedingRegion &a, const TICLSeedingRegion &b) {
    return a.directionAtOrigin.perp2() > b.directionAtOrigin.perp2();
  });
}

double ticl::SeedingRegionByL1::tkEmOfflineEt(double et) const {
  return (endcapScalings_.at(0) + et * endcapScalings_.at(1) + et * et * endcapScalings_.at(2));
}

void ticl::SeedingRegionByL1::fillPSetDescription(edm::ParameterSetDescription &desc) {
  desc.add<edm::InputTag>("l1TkEmColl", edm::InputTag("L1TkPhotonsHGC", "EG"));
  desc.add<double>("minPt", 10);
  desc.add<double>("minAbsEta", 1.479);
  desc.add<double>("maxAbsEta", 4.0);
  desc.add<std::vector<double>>("endcapScalings", {3.17445, 1.13219, 0.0});
  desc.add<int>("quality", 5);
  SeedingRegionAlgoBase::fillPSetDescription(desc);
}