Phase2L1CaloJetEmulator Class Reference

#include <L1Trigger/L1CaloTrigger/plugins/Phase2L1CaloJetEmulator.cc>

Inheritance diagram for Phase2L1CaloJetEmulator:

Public Member Functions
	Phase2L1CaloJetEmulator (const edm::ParameterSet &)
	~Phase2L1CaloJetEmulator () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
edm::EDGetTokenT< l1tp2::CaloTowerCollection >	caloTowerToken_
edm::ESGetToken< CaloTPGTranscoder, CaloTPGRecord >	decoderTag_
edm::EDGetTokenT< HcalTrigPrimDigiCollection >	hfToken_
edm::EDGetTokenT< l1t::HGCalTowerBxCollection >	hgcalTowerToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Description: Producing GCT calo jets using GCT barrel, HGCal and HF towers, based on firmware logic.

Implementation: Depends on producers for CaloTowerCollection, HGCalTowerBxCollection and HcalTrigPrimDigiCollection.

Definition at line 55 of file Phase2L1CaloJetEmulator.cc.

Constructor & Destructor Documentation

Phase2L1CaloJetEmulator()

Phase2L1CaloJetEmulator::Phase2L1CaloJetEmulator ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 75 of file Phase2L1CaloJetEmulator.cc.

    : caloTowerToken_(consumes<l1tp2::CaloTowerCollection>(iConfig.getParameter<edm::InputTag>("gctFullTowers"))),
      hgcalTowerToken_(consumes<l1t::HGCalTowerBxCollection>(iConfig.getParameter<edm::InputTag>("hgcalTowers"))),
      hfToken_(consumes<HcalTrigPrimDigiCollection>(iConfig.getParameter<edm::InputTag>("hcalDigis"))),
      decoderTag_(esConsumes<CaloTPGTranscoder, CaloTPGRecord>(edm::ESInputTag("", ""))) {
  produces<l1tp2::Phase2L1CaloJetCollection>("GCTJet");
}

~Phase2L1CaloJetEmulator()

Phase2L1CaloJetEmulator::~Phase2L1CaloJetEmulator ( )

override

Definition at line 83 of file Phase2L1CaloJetEmulator.cc.

{}

Member Function Documentation

fillDescriptions()

void Phase2L1CaloJetEmulator::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 457 of file Phase2L1CaloJetEmulator.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                         {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("gctFullTowers", edm::InputTag("l1tPhase2L1CaloEGammaEmulator", "GCTFullTowers"));
  desc.add<edm::InputTag>("hgcalTowers", edm::InputTag("l1tHGCalTowerProducer", "HGCalTowerProcessor"));
  desc.add<edm::InputTag>("hcalDigis", edm::InputTag("simHcalTriggerPrimitiveDigis"));
  descriptions.addWithDefaultLabel(desc);
}

produce()

void Phase2L1CaloJetEmulator::produce ( edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overrideprivate

Definition at line 90 of file Phase2L1CaloJetEmulator.cc.

References funct::abs(), reco::JetExtendedAssociation::allJets(), BXVector< T >::begin(), caloTowerToken_, gctobj::compareByEt(), decoderTag_, BXVector< T >::end(), hcalRecHitTable_cff::energy, EgHLTOffHistBins_cfi::et, PVValHelper::eta, edm::EventSetup::getData(), gctobj::getRegion(), hfToken_, l1tPhase2CaloJetEmulator_cfi::hgcalTowers, hgcalTowerToken_, mps_fire::i, hit::id, hcalRecHitTable_cff::ieta, iEvent, hcalRecHitTable_cff::iphi, dqmiolumiharvest::j, metsig::jet, l1tp2::Phase2L1CaloJet::jetEt(), l1tp2::Phase2L1CaloJet::jetEta(), l1tp2::Phase2L1CaloJet::jetPhi(), dqmdumpme::k, l1t::CaloTools::kHFBegin, l1t::CaloTools::kHFEnd, M_PI, gctobj::makeST(), gctobj::makeST_hf(), gctobj::makeST_hgcal(), eostools::move(), nBarrelEta, nBarrelPhi, nHfEta, nHfPhi, nHgcalEta, nHgcalPhi, nJets, nSTEta, nSTPhi, edm::Handle< T >::product(), l1tp2::Phase2L1CaloJet::setJetEt(), l1tp2::Phase2L1CaloJet::setJetEta(), l1tp2::Phase2L1CaloJet::setJetIEta(), l1tp2::Phase2L1CaloJet::setJetIPhi(), l1tp2::Phase2L1CaloJet::setJetPhi(), reco::LeafCandidate::setP4(), l1tp2::Phase2L1CaloJet::setTauEt(), l1tp2::Phase2L1CaloJet::setTowerEt(), l1tp2::Phase2L1CaloJet::setTowerEta(), l1tp2::Phase2L1CaloJet::setTowerIEta(), l1tp2::Phase2L1CaloJet::setTowerIPhi(), l1tp2::Phase2L1CaloJet::setTowerPhi(), jetUpdater_cfi::sort, groupFilesInBlocks::temp, and l1t::CaloTools::towerEta().

                                                                                   {
  using namespace edm;
  std::unique_ptr<l1tp2::Phase2L1CaloJetCollection> jetCands(make_unique<l1tp2::Phase2L1CaloJetCollection>());

  // Assign ETs to each eta-half of the barrel region (17x72 --> 18x72 to be able to make 3x3 super towers)
  edm::Handle<std::vector<l1tp2::CaloTower>> caloTowerCollection;
  if (!iEvent.getByToken(caloTowerToken_, caloTowerCollection))
    edm::LogError("Phase2L1CaloJetEmulator") << "Failed to get towers from caloTowerCollection!";

  iEvent.getByToken(caloTowerToken_, caloTowerCollection);
  float GCTintTowers[nBarrelEta][nBarrelPhi];
  float realEta[nBarrelEta][nBarrelPhi];
  float realPhi[nBarrelEta][nBarrelPhi];
  for (const l1tp2::CaloTower& i : *caloTowerCollection) {
    int ieta = i.towerIEta();
    int iphi = i.towerIPhi();
    if (i.ecalTowerEt() > 1.)
      GCTintTowers[ieta][iphi] = i.ecalTowerEt();  // suppress <= 1 GeV towers
    else
      GCTintTowers[ieta][iphi] = 0;
    realEta[ieta][iphi] = i.towerEta();
    realPhi[ieta][iphi] = i.towerPhi();
  }

  float temporary[nBarrelEta / 2][nBarrelPhi];

  // Assign ETs to each eta-half of the endcap region (18x72)
  edm::Handle<l1t::HGCalTowerBxCollection> hgcalTowerCollection;
  if (!iEvent.getByToken(hgcalTowerToken_, hgcalTowerCollection))
    edm::LogError("Phase2L1CaloJetEmulator") << "Failed to get towers from hgcalTowerCollection!";
  l1t::HGCalTowerBxCollection hgcalTowerColl;
  iEvent.getByToken(hgcalTowerToken_, hgcalTowerCollection);
  hgcalTowerColl = (*hgcalTowerCollection.product());
  float hgcalTowers[nHgcalEta][nHgcalPhi];
  float hgcalEta[nHgcalEta][nHgcalPhi];
  float hgcalPhi[nHgcalEta][nHgcalPhi];

  for (int iphi = 0; iphi < nHgcalPhi; iphi++) {
    for (int ieta = 0; ieta < nHgcalEta; ieta++) {
      hgcalTowers[ieta][iphi] = 0;
      if (ieta < nHgcalEta / 2)
        hgcalEta[ieta][iphi] = -3.045 + ieta * 0.087 + 0.0435;
      else
        hgcalEta[ieta][iphi] = 1.479 + (ieta - nHgcalEta / 2) * 0.087 + 0.0435;
      hgcalPhi[ieta][iphi] = -M_PI + (iphi * 2 * M_PI / nHgcalPhi) + (M_PI / nHgcalPhi);
    }
  }

  for (auto it = hgcalTowerColl.begin(0); it != hgcalTowerColl.end(0); it++) {
    float eta = it->eta();
    int ieta;
    if (eta < 0)
      ieta = 19 - it->id().iEta();
    else
      ieta = 20 + it->id().iEta();
    if (eta > 1.479)
      ieta = ieta - 4;
    int iphi = it->id().iPhi();
    if ((it->etEm() + it->etHad() > 1.) && abs(eta) > 1.479)
      hgcalTowers[ieta][iphi] = it->etEm() + it->etHad();  // suppress <= 1 GeV towers
  }

  float temporary_hgcal[nHgcalEta / 2][nHgcalPhi];

  // Assign ETs to each eta-half of the forward region (12x72)
  edm::Handle<HcalTrigPrimDigiCollection> hfHandle;
  if (!iEvent.getByToken(hfToken_, hfHandle))
    edm::LogError("Phase2L1CaloJetEmulator") << "Failed to get HcalTrigPrimDigi for HF!";
  iEvent.getByToken(hfToken_, hfHandle);
  float hfTowers[nHfEta][nHfPhi];
  float hfEta[nHfEta][nHfPhi];
  float hfPhi[nHfEta][nHfPhi];

  for (int iphi = 0; iphi < nHfPhi; iphi++) {
    for (int ieta = 0; ieta < nHfEta; ieta++) {
      hfTowers[ieta][iphi] = 0;
      int temp;
      if (ieta < nHfEta / 2)
        temp = ieta - l1t::CaloTools::kHFEnd;
      else
        temp = ieta - nHfEta / 2 + l1t::CaloTools::kHFBegin + 1;
      hfEta[ieta][iphi] = l1t::CaloTools::towerEta(temp);
      hfPhi[ieta][iphi] = -M_PI + (iphi * 2 * M_PI / nHfPhi) + (M_PI / nHfPhi);
    }
  }

  const auto& decoder = iSetup.getData(decoderTag_);
  for (const auto& hit : *hfHandle.product()) {
    double et = decoder.hcaletValue(hit.id(), hit.t0());
    int ieta = 0;
    if (abs(hit.id().ieta()) < l1t::CaloTools::kHFBegin)
      continue;
    if (abs(hit.id().ieta()) > l1t::CaloTools::kHFEnd)
      continue;
    if (hit.id().ieta() <= -(l1t::CaloTools::kHFBegin + 1)) {
      ieta = hit.id().ieta() + l1t::CaloTools::kHFEnd;
    } else if (hit.id().ieta() >= (l1t::CaloTools::kHFBegin + 1)) {
      ieta = nHfEta / 2 + (hit.id().ieta() - (l1t::CaloTools::kHFBegin + 1));
    }
    int iphi = 0;
    if (hit.id().iphi() <= nHfPhi / 2)
      iphi = hit.id().iphi() + (nHfPhi / 2 - 1);
    else if (hit.id().iphi() > nHfPhi / 2)
      iphi = hit.id().iphi() - (nHfPhi / 2 + 1);
    if (et > 1.)
      hfTowers[ieta][iphi] = et;  // suppress <= 1 GeV towers
  }

  float temporary_hf[nHfEta / 2][nHfPhi];

  // Begin creating jets
  // First create 3x3 super towers: 6x24 in barrel, endcap; 4x24 in forward
  // Then create up to 10 jets in each eta half of barrel, endcap, forward regions

  vector<l1tp2::Phase2L1CaloJet> halfBarrelJets, halfHgcalJets, halfHfJets;
  halfBarrelJets.clear();
  halfHgcalJets.clear();
  halfHfJets.clear();
  vector<l1tp2::Phase2L1CaloJet> allJets;
  allJets.clear();

  for (int k = 0; k < 2; k++) {
    halfBarrelJets.clear();
    halfHgcalJets.clear();
    halfHfJets.clear();
    gctobj::jetInfo jet[3 * nJets];

    // BARREL
    for (int iphi = 0; iphi < nBarrelPhi; iphi++) {
      for (int ieta = 0; ieta < nBarrelEta / 2; ieta++) {
        if (k == 0)
          temporary[ieta][iphi] = GCTintTowers[ieta][iphi];
        else
          temporary[ieta][iphi] = GCTintTowers[nBarrelEta / 2 + ieta][iphi];
      }
    }

    gctobj::GCTsupertower_t tempST[nSTEta][nSTPhi];
    gctobj::makeST(temporary, tempST);

    for (int i = 0; i < nJets; i++) {
      jet[i] = gctobj::getRegion(tempST);
      l1tp2::Phase2L1CaloJet tempJet;
      tempJet.setJetEt(jet[i].energy);
      tempJet.setTauEt(jet[i].tauEt);
      int gctjeteta = jet[i].etaCenter;
      int gctjetphi = jet[i].phiCenter;
      tempJet.setJetIEta(gctjeteta + k * nBarrelEta / 2);
      tempJet.setJetIPhi(gctjetphi);
      float jeteta = realEta[gctjeteta + k * nBarrelEta / 2][gctjetphi];
      float jetphi = realPhi[gctjeteta + k * nBarrelEta / 2][gctjetphi];
      tempJet.setJetEta(jeteta);
      tempJet.setJetPhi(jetphi);
      tempJet.setTowerEt(jet[i].energyMax);
      int gcttowereta = jet[i].etaMax;
      int gcttowerphi = jet[i].phiMax;
      tempJet.setTowerIEta(gcttowereta + k * nBarrelEta / 2);
      tempJet.setTowerIPhi(gcttowerphi);
      float towereta = realEta[gcttowereta + k * nBarrelEta / 2][gcttowerphi];
      float towerphi = realPhi[gcttowereta + k * nBarrelEta / 2][gcttowerphi];
      tempJet.setTowerEta(towereta);
      tempJet.setTowerPhi(towerphi);
      reco::Candidate::PolarLorentzVector tempJetp4;
      tempJetp4.SetPt(tempJet.jetEt());
      tempJetp4.SetEta(tempJet.jetEta());
      tempJetp4.SetPhi(tempJet.jetPhi());
      tempJetp4.SetM(0.);
      tempJet.setP4(tempJetp4);

      if (jet[i].energy > 0.)
        halfBarrelJets.push_back(tempJet);
    }

    // ENDCAP
    for (int iphi = 0; iphi < nHgcalPhi; iphi++) {
      for (int ieta = 0; ieta < nHgcalEta / 2; ieta++) {
        if (k == 0)
          temporary_hgcal[ieta][iphi] = hgcalTowers[ieta][iphi];
        else
          temporary_hgcal[ieta][iphi] = hgcalTowers[nHgcalEta / 2 + ieta][iphi];
      }
    }

    gctobj::GCTsupertower_t tempST_hgcal[nSTEta][nSTPhi];
    gctobj::makeST_hgcal(temporary_hgcal, tempST_hgcal);
    for (int i = nJets; i < 2 * nJets; i++) {
      jet[i] = gctobj::getRegion(tempST_hgcal);
      l1tp2::Phase2L1CaloJet tempJet;
      tempJet.setJetEt(jet[i].energy);
      tempJet.setTauEt(jet[i].tauEt);
      int hgcaljeteta = jet[i].etaCenter;
      int hgcaljetphi = jet[i].phiCenter;
      tempJet.setJetIEta(hgcaljeteta + k * nHgcalEta / 2);
      tempJet.setJetIPhi(hgcaljetphi);
      float jeteta = hgcalEta[hgcaljeteta + k * nHgcalEta / 2][hgcaljetphi];
      float jetphi = hgcalPhi[hgcaljeteta + k * nHgcalEta / 2][hgcaljetphi];
      tempJet.setJetEta(jeteta);
      tempJet.setJetPhi(jetphi);
      tempJet.setTowerEt(jet[i].energyMax);
      int hgcaltowereta = jet[i].etaMax;
      int hgcaltowerphi = jet[i].phiMax;
      tempJet.setTowerIEta(hgcaltowereta + k * nHgcalEta / 2);
      tempJet.setTowerIPhi(hgcaltowerphi);
      float towereta = hgcalEta[hgcaltowereta + k * nHgcalEta / 2][hgcaltowerphi];
      float towerphi = hgcalPhi[hgcaltowereta + k * nHgcalEta / 2][hgcaltowerphi];
      tempJet.setTowerEta(towereta);
      tempJet.setTowerPhi(towerphi);
      reco::Candidate::PolarLorentzVector tempJetp4;
      tempJetp4.SetPt(tempJet.jetEt());
      tempJetp4.SetEta(tempJet.jetEta());
      tempJetp4.SetPhi(tempJet.jetPhi());
      tempJetp4.SetM(0.);
      tempJet.setP4(tempJetp4);

      if (jet[i].energy > 0.)
        halfHgcalJets.push_back(tempJet);
    }

    // HF
    for (int iphi = 0; iphi < nHfPhi; iphi++) {
      for (int ieta = 0; ieta < nHfEta / 2; ieta++) {
        if (k == 0)
          temporary_hf[ieta][iphi] = hfTowers[ieta][iphi];
        else
          temporary_hf[ieta][iphi] = hfTowers[nHfEta / 2 + ieta][iphi];
      }
    }

    gctobj::GCTsupertower_t tempST_hf[nSTEta][nSTPhi];
    gctobj::makeST_hf(temporary_hf, tempST_hf);
    for (int i = 2 * nJets; i < 3 * nJets; i++) {
      jet[i] = gctobj::getRegion(tempST_hf);
      l1tp2::Phase2L1CaloJet tempJet;
      tempJet.setJetEt(jet[i].energy);
      tempJet.setTauEt(jet[i].tauEt);
      int hfjeteta = jet[i].etaCenter;
      int hfjetphi = jet[i].phiCenter;
      tempJet.setJetIEta(hfjeteta + k * nHfEta / 2);
      tempJet.setJetIPhi(hfjetphi);
      float jeteta = hfEta[hfjeteta + k * nHfEta / 2][hfjetphi];
      float jetphi = hfPhi[hfjeteta + k * nHfEta / 2][hfjetphi];
      tempJet.setJetEta(jeteta);
      tempJet.setJetPhi(jetphi);
      tempJet.setTowerEt(jet[i].energyMax);
      int hftowereta = jet[i].etaMax;
      int hftowerphi = jet[i].phiMax;
      tempJet.setTowerIEta(hftowereta + k * nHfEta / 2);
      tempJet.setTowerIPhi(hftowerphi);
      float towereta = hfEta[hftowereta + k * nHfEta / 2][hftowerphi];
      float towerphi = hfPhi[hftowereta + k * nHfEta / 2][hftowerphi];
      tempJet.setTowerEta(towereta);
      tempJet.setTowerPhi(towerphi);
      reco::Candidate::PolarLorentzVector tempJetp4;
      tempJetp4.SetPt(tempJet.jetEt());
      tempJetp4.SetEta(tempJet.jetEta());
      tempJetp4.SetPhi(tempJet.jetPhi());
      tempJetp4.SetM(0.);
      tempJet.setP4(tempJetp4);

      if (jet[i].energy > 0.)
        halfHfJets.push_back(tempJet);
    }

    // Stitching:
    // if the jet eta is at the boundary: for HB it should be within 0-1 in -ve eta, 32-33 in +ve eta; for HE it should be within 0-1/16-17 in -ve eta, 34-35/18-19 in +ve eta; for HF it should be within 10-11 in -ve eta, 12-13 in +ve eta
    // then get the phi of that jet and check if there is a neighbouring jet with the same phi, then merge to the jet that has higher ET
    // in both eta/phi allow a maximum of one tower between jet centers for stitching

    for (size_t i = 0; i < halfHgcalJets.size(); i++) {
      if (halfHgcalJets.at(i).jetIEta() >= (nHgcalEta / 2 - 2) && halfHgcalJets.at(i).jetIEta() < (nHgcalEta / 2 + 2)) {
        float hgcal_ieta = k * nBarrelEta + halfHgcalJets.at(i).jetIEta();
        for (size_t j = 0; j < halfBarrelJets.size(); j++) {
          float barrel_ieta = nHgcalEta / 2 + halfBarrelJets.at(j).jetIEta();
          if (abs(barrel_ieta - hgcal_ieta) <= 2 &&
              abs(halfBarrelJets.at(j).jetIPhi() - halfHgcalJets.at(i).jetIPhi()) <= 2) {
            float totalet = halfBarrelJets.at(j).jetEt() + halfHgcalJets.at(i).jetEt();
            float totalTauEt = halfBarrelJets.at(j).tauEt() + halfHgcalJets.at(i).tauEt();
            if (halfBarrelJets.at(j).jetEt() > halfHgcalJets.at(i).jetEt()) {
              halfHgcalJets.at(i).setJetEt(0.);
              halfHgcalJets.at(i).setTauEt(0.);
              halfBarrelJets.at(j).setJetEt(totalet);
              halfBarrelJets.at(j).setTauEt(totalTauEt);
              reco::Candidate::PolarLorentzVector tempJetp4;
              tempJetp4.SetPt(totalet);
              tempJetp4.SetEta(halfBarrelJets.at(j).jetEta());
              tempJetp4.SetPhi(halfBarrelJets.at(j).jetPhi());
              tempJetp4.SetM(0.);
              halfBarrelJets.at(j).setP4(tempJetp4);
            } else {
              halfHgcalJets.at(i).setJetEt(totalet);
              halfHgcalJets.at(i).setTauEt(totalTauEt);
              halfBarrelJets.at(j).setJetEt(0.);
              halfBarrelJets.at(j).setTauEt(0.);
              reco::Candidate::PolarLorentzVector tempJetp4;
              tempJetp4.SetPt(totalet);
              tempJetp4.SetEta(halfHgcalJets.at(i).jetEta());
              tempJetp4.SetPhi(halfHgcalJets.at(i).jetPhi());
              tempJetp4.SetM(0.);
              halfHgcalJets.at(i).setP4(tempJetp4);
            }
          }
        }
      } else if (halfHgcalJets.at(i).jetIEta() < 2 || halfHgcalJets.at(i).jetIEta() >= (nHgcalEta - 2)) {
        float hgcal_ieta = k * nBarrelEta + nHfEta / 2 + halfHgcalJets.at(i).jetIEta();
        for (size_t j = 0; j < halfHfJets.size(); j++) {
          float hf_ieta = k * nBarrelEta + k * nHgcalEta + halfHfJets.at(j).jetIEta();
          if (abs(hgcal_ieta - hf_ieta) < 3 && abs(halfHfJets.at(j).jetIPhi() - halfHgcalJets.at(i).jetIPhi()) < 3) {
            float totalet = halfHfJets.at(j).jetEt() + halfHgcalJets.at(i).jetEt();
            float totalTauEt = halfHfJets.at(j).tauEt() + halfHgcalJets.at(i).tauEt();
            if (halfHfJets.at(j).jetEt() > halfHgcalJets.at(i).jetEt()) {
              halfHgcalJets.at(i).setJetEt(0.);
              halfHgcalJets.at(i).setTauEt(0.);
              halfHfJets.at(j).setJetEt(totalet);
              halfHfJets.at(j).setTauEt(totalTauEt);
              reco::Candidate::PolarLorentzVector tempJetp4;
              tempJetp4.SetPt(totalet);
              tempJetp4.SetEta(halfHfJets.at(j).jetEta());
              tempJetp4.SetPhi(halfHfJets.at(j).jetPhi());
              tempJetp4.SetM(0.);
              halfHfJets.at(j).setP4(tempJetp4);
            } else {
              halfHgcalJets.at(i).setJetEt(totalet);
              halfHgcalJets.at(i).setTauEt(totalTauEt);
              halfHfJets.at(j).setJetEt(0.);
              halfHfJets.at(j).setTauEt(0.);
              reco::Candidate::PolarLorentzVector tempJetp4;
              tempJetp4.SetPt(totalet);
              tempJetp4.SetEta(halfHgcalJets.at(i).jetEta());
              tempJetp4.SetPhi(halfHgcalJets.at(i).jetPhi());
              tempJetp4.SetM(0.);
              halfHgcalJets.at(i).setP4(tempJetp4);
            }
          }
        }
      }
    }

    // Write up to 6 jets from each eta half of barrel, endcap, forward regions

    std::sort(halfBarrelJets.begin(), halfBarrelJets.end(), gctobj::compareByEt);
    for (size_t i = 0; i < halfBarrelJets.size(); i++) {
      if (halfBarrelJets.at(i).jetEt() > 0. && i < 6)
        allJets.push_back(halfBarrelJets.at(i));
    }

    std::sort(halfHgcalJets.begin(), halfHgcalJets.end(), gctobj::compareByEt);
    for (size_t i = 0; i < halfHgcalJets.size(); i++) {
      if (halfHgcalJets.at(i).jetEt() > 0. && i < 6)
        allJets.push_back(halfHgcalJets.at(i));
    }

    std::sort(halfHfJets.begin(), halfHfJets.end(), gctobj::compareByEt);
    for (size_t i = 0; i < halfHfJets.size(); i++) {
      if (halfHfJets.at(i).jetEt() > 0. && i < 6)
        allJets.push_back(halfHfJets.at(i));
    }
  }

  std::sort(allJets.begin(), allJets.end(), gctobj::compareByEt);
  for (size_t i = 0; i < allJets.size(); i++) {
    jetCands->push_back(allJets.at(i));
  }

  iEvent.put(std::move(jetCands), "GCTJet");
}

Member Data Documentation

caloTowerToken_

edm::EDGetTokenT<l1tp2::CaloTowerCollection> Phase2L1CaloJetEmulator::caloTowerToken_

private

Definition at line 66 of file Phase2L1CaloJetEmulator.cc.

Referenced by produce().

decoderTag_

edm::ESGetToken<CaloTPGTranscoder, CaloTPGRecord> Phase2L1CaloJetEmulator::decoderTag_

private

Definition at line 69 of file Phase2L1CaloJetEmulator.cc.

Referenced by produce().

hfToken_

edm::EDGetTokenT<HcalTrigPrimDigiCollection> Phase2L1CaloJetEmulator::hfToken_

private

Definition at line 68 of file Phase2L1CaloJetEmulator.cc.

Referenced by produce().

hgcalTowerToken_

edm::EDGetTokenT<l1t::HGCalTowerBxCollection> Phase2L1CaloJetEmulator::hgcalTowerToken_

private

Definition at line 67 of file Phase2L1CaloJetEmulator.cc.

Referenced by produce().