L1GlobalCaloTrigger.cc

Go to the documentation of this file.

#include "L1Trigger/GlobalCaloTrigger/interface/L1GlobalCaloTrigger.h"
 
#include "CondFormats/L1TObjects/interface/L1GctJetFinderParams.h"
#include "CondFormats/L1TObjects/interface/L1GctChannelMask.h"
 
#include "L1Trigger/GlobalCaloTrigger/interface/L1GctJetEtCalibrationLut.h"
#include "L1Trigger/GlobalCaloTrigger/interface/L1GctEmLeafCard.h"
#include "L1Trigger/GlobalCaloTrigger/interface/L1GctWheelJetFpga.h"
#include "L1Trigger/GlobalCaloTrigger/interface/L1GctWheelEnergyFpga.h"
#include "L1Trigger/GlobalCaloTrigger/interface/L1GctJetFinalStage.h"
#include "L1Trigger/GlobalCaloTrigger/interface/L1GctGlobalEnergyAlgos.h"
#include "L1Trigger/GlobalCaloTrigger/interface/L1GctGlobalHfSumAlgos.h"
#include "L1Trigger/GlobalCaloTrigger/interface/L1GctElectronFinalSort.h"
 
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
using std::vector;
 
//DEFINE STATICS
const int L1GlobalCaloTrigger::N_JET_LEAF_CARDS = 6;
const int L1GlobalCaloTrigger::N_EM_LEAF_CARDS = 2;
const int L1GlobalCaloTrigger::N_WHEEL_CARDS = 2;
 
// constructor
L1GlobalCaloTrigger::L1GlobalCaloTrigger(const L1GctJetLeafCard::jetFinderType jfType, unsigned jetLeafMask)
    : theJetLeafCards(N_JET_LEAF_CARDS),
      theJetFinders(N_JET_LEAF_CARDS * 3),
      theEmLeafCards(N_EM_LEAF_CARDS),
      theIsoElectronSorters(N_EM_LEAF_CARDS * 2),
      theNonIsoElectronSorters(N_EM_LEAF_CARDS * 2),
      theWheelJetFpgas(N_WHEEL_CARDS),
      theWheelEnergyFpgas(N_WHEEL_CARDS),
      m_jetFinderParams(nullptr),
      m_jetEtCalLuts(),
      m_inputChannelMask(nullptr),
      m_bxRangeAuto(true),
      m_bxStart(0),
      m_numOfBx(1),
      m_allInputEmCands(),
      m_allInputRegions() {
  // construct hardware
  build(jfType, jetLeafMask);
}
 
L1GlobalCaloTrigger::~L1GlobalCaloTrigger() {
  // Delete the components of the GCT that we made in build()
  // (But not the LUTs, since these don't belong to us)
 
  if (theNonIsoEmFinalStage != nullptr)
    delete theNonIsoEmFinalStage;
 
  if (theIsoEmFinalStage != nullptr)
    delete theIsoEmFinalStage;
 
  if (theEnergyFinalStage != nullptr)
    delete theEnergyFinalStage;
 
  if (theJetFinalStage != nullptr)
    delete theJetFinalStage;
 
  for (unsigned i = 0; i < theWheelEnergyFpgas.size(); ++i) {
    if (theWheelEnergyFpgas.at(i) != nullptr)
      delete theWheelEnergyFpgas.at(i);
  }
  theWheelEnergyFpgas.clear();
 
  for (unsigned i = 0; i < theWheelJetFpgas.size(); ++i) {
    if (theWheelJetFpgas.at(i) != nullptr)
      delete theWheelJetFpgas.at(i);
  }
  theWheelJetFpgas.clear();
 
  for (unsigned i = 0; i < theEmLeafCards.size(); ++i) {
    if (theEmLeafCards.at(i) != nullptr)
      delete theEmLeafCards.at(i);
  }
  theEmLeafCards.clear();
 
  for (unsigned i = 0; i < theJetLeafCards.size(); ++i) {
    if (theJetLeafCards.at(i) != nullptr)
      delete theJetLeafCards.at(i);
  }
  theJetLeafCards.clear();
}
 
void L1GlobalCaloTrigger::reset() {
  // Input data
  m_allInputEmCands.clear();
  m_allInputRegions.clear();
 
  if (m_bxRangeAuto) {
    m_bxStart = 0;
    m_numOfBx = 1;
  }
 
  // EM Leaf Card
  for (int i = 0; i < N_EM_LEAF_CARDS; i++) {
    theEmLeafCards.at(i)->reset();
  }
 
  // Jet Leaf cards
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->reset();
  }
 
  // Jet Finders
  for (int i = 0; i < N_JET_LEAF_CARDS * 3; i++) {
    theJetFinders.at(i)->reset();
  }
 
  // Wheel Cards
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelJetFpgas.at(i)->reset();
  }
 
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelEnergyFpgas.at(i)->reset();
  }
 
  // Electron Final Stage
  theIsoEmFinalStage->reset();
  theNonIsoEmFinalStage->reset();
 
  // Jet Final Stage
  theJetFinalStage->reset();
 
  // Energy Final Stage
  theEnergyFinalStage->reset();
}
 
void L1GlobalCaloTrigger::process() {
  // Shouldn't get here unless the setup has been completed
  if (setupOk()) {
    sortInputData();
    // Extract the earliest and latest bunch crossing
    // in the input if required, and forward to the processors
    // to determine the size of the output vectors
    bxSetup();
 
    vector<L1CaloEmCand>::iterator emc = m_allInputEmCands.begin();
    vector<L1CaloRegion>::iterator rgn = m_allInputRegions.begin();
    int bx = m_bxStart;
 
    // Loop over bunch crossings
    for (int i = 0; i < m_numOfBx; i++) {
      // Perform partial reset (reset processing logic but preserve pipeline contents)
      bxReset(bx);
      // Fill input data into processors for this bunch crossing
      fillEmCands(emc, bx);
      fillRegions(rgn, bx);
      // Process this bunch crossing
      bxProcess(bx);
      bx++;
    }
  }
}
 
void L1GlobalCaloTrigger::sortInputData() {
  std::sort(m_allInputEmCands.begin(), m_allInputEmCands.end(), emcBxComparator);
  std::sort(m_allInputRegions.begin(), m_allInputRegions.end(), rgnBxComparator);
}
 
void L1GlobalCaloTrigger::bxSetup() {
  // Assume input data have been sorted by bunch crossing number
  if (m_bxRangeAuto) {
    // Find parameters defining the range of bunch crossings to be processed
    int16_t firstBxEmCand = (m_allInputEmCands.empty() ? 0 : m_allInputEmCands.front().bx());
    int16_t firstBxRegion = (m_allInputRegions.empty() ? 0 : m_allInputRegions.front().bx());
    int16_t lastBxEmCand = (m_allInputEmCands.empty() ? 0 : m_allInputEmCands.back().bx());
    int16_t lastBxRegion = (m_allInputRegions.empty() ? 0 : m_allInputRegions.back().bx());
    m_bxStart = std::min(firstBxEmCand, firstBxRegion);
    m_numOfBx = std::max(lastBxEmCand, lastBxRegion) - m_bxStart + 1;
  } else {
    // Remove any input from before the start of the requested range
    for (vector<L1CaloEmCand>::iterator emc = m_allInputEmCands.begin(); emc != m_allInputEmCands.end(); emc++) {
      if (emc->bx() >= m_bxStart)
        break;
      m_allInputEmCands.erase(emc);
    }
 
    for (vector<L1CaloRegion>::iterator rgn = m_allInputRegions.begin(); rgn != m_allInputRegions.end(); rgn++) {
      if (rgn->bx() >= m_bxStart)
        break;
      m_allInputRegions.erase(rgn);
    }
  }
 
  // Setup pipeline lengths
  // EM Leaf Card
  for (int i = 0; i < N_EM_LEAF_CARDS; i++) {
    theEmLeafCards.at(i)->setBxRange(m_bxStart, m_numOfBx);
  }
 
  // Jet Leaf cards
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->setBxRange(m_bxStart, m_numOfBx);
  }
 
  // Wheel Cards
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelJetFpgas.at(i)->setBxRange(m_bxStart, m_numOfBx);
  }
 
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelEnergyFpgas.at(i)->setBxRange(m_bxStart, m_numOfBx);
  }
 
  // Electron Final Stage
  theIsoEmFinalStage->setBxRange(m_bxStart, m_numOfBx);
  theNonIsoEmFinalStage->setBxRange(m_bxStart, m_numOfBx);
 
  // Jet Final Stage
  theJetFinalStage->setBxRange(m_bxStart, m_numOfBx);
 
  // Energy Final Stage
  theEnergyFinalStage->setBxRange(m_bxStart, m_numOfBx);
}
 
void L1GlobalCaloTrigger::bxReset(const int bx) {
  // EM Leaf Card
  for (int i = 0; i < N_EM_LEAF_CARDS; i++) {
    theEmLeafCards.at(i)->setNextBx(bx);
  }
 
  // Jet Leaf cards
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->setNextBx(bx);
  }
 
  // Wheel Cards
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelJetFpgas.at(i)->setNextBx(bx);
  }
 
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelEnergyFpgas.at(i)->setNextBx(bx);
  }
 
  // Electron Final Stage
  theIsoEmFinalStage->setNextBx(bx);
  theNonIsoEmFinalStage->setNextBx(bx);
 
  // Jet Final Stage
  theJetFinalStage->setNextBx(bx);
 
  // Energy Final Stage
  theEnergyFinalStage->setNextBx(bx);
}
 
void L1GlobalCaloTrigger::bxProcess(const int bx) {
  // EM Leaf Card
  for (int i = 0; i < N_EM_LEAF_CARDS; i++) {
    theEmLeafCards.at(i)->fetchInput();
    theEmLeafCards.at(i)->process();
  }
 
  // Jet Leaf cards - first stage processing
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->fetchInput();
  }
 
  // Jet Leaf cards - second stage processing
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->process();
  }
 
  // Wheel Cards
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelJetFpgas.at(i)->fetchInput();
    theWheelJetFpgas.at(i)->process();
  }
 
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelEnergyFpgas.at(i)->fetchInput();
    theWheelEnergyFpgas.at(i)->process();
  }
 
  // Electron Final Stage
  theIsoEmFinalStage->fetchInput();
  theIsoEmFinalStage->process();
 
  theNonIsoEmFinalStage->fetchInput();
  theNonIsoEmFinalStage->process();
 
  // Jet Final Stage
  theJetFinalStage->fetchInput();
  theJetFinalStage->process();
 
  // Energy Final Stage
  theEnergyFinalStage->fetchInput();
  theEnergyFinalStage->process();
}
 
void L1GlobalCaloTrigger::setJetFinderParams(const L1GctJetFinderParams* const jfpars) {
  // Some parameters not (yet?) implemented
  if ((jfpars->getCenForJetEtaBoundary() == 7) && (jfpars->getCenJetEtSeedGct() == jfpars->getTauJetEtSeedGct())) {
    m_jetFinderParams = jfpars;
    // Need to propagate the new parameters to all the JetFinders
    for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
      theJetLeafCards.at(i)->getJetFinderA()->setJetFinderParams(jfpars);
      theJetLeafCards.at(i)->getJetFinderB()->setJetFinderParams(jfpars);
      theJetLeafCards.at(i)->getJetFinderC()->setJetFinderParams(jfpars);
    }
    // Also send to the final energy calculation (for missing Ht)
    theEnergyFinalStage->setJetFinderParams(jfpars);
  }
}
 
void L1GlobalCaloTrigger::setJetEtCalibrationLuts(const L1GlobalCaloTrigger::lutPtrVector& jfluts) {
  m_jetEtCalLuts = jfluts;
  // Need to propagate the new lut to all the JetFinders
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->getJetFinderA()->setJetEtCalibrationLuts(jfluts);
    theJetLeafCards.at(i)->getJetFinderB()->setJetEtCalibrationLuts(jfluts);
    theJetLeafCards.at(i)->getJetFinderC()->setJetEtCalibrationLuts(jfluts);
  }
}
 
void L1GlobalCaloTrigger::setupTauAlgo(const bool useImprovedAlgo, const bool ignoreVetoBitsForIsolation) {
  // Need to propagate the new parameters to all the JetFinders
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->getJetFinderA()->setupTauAlgo(useImprovedAlgo, ignoreVetoBitsForIsolation);
    theJetLeafCards.at(i)->getJetFinderB()->setupTauAlgo(useImprovedAlgo, ignoreVetoBitsForIsolation);
    theJetLeafCards.at(i)->getJetFinderC()->setupTauAlgo(useImprovedAlgo, ignoreVetoBitsForIsolation);
  }
}
 
void L1GlobalCaloTrigger::setHtMissScale(const L1CaloEtScale* const scale) {
  if (theEnergyFinalStage != nullptr) {
    theEnergyFinalStage->setHtMissScale(scale);
  }
}
 
void L1GlobalCaloTrigger::setupHfSumLuts(const L1CaloEtScale* const scale) {
  if (getHfSumProcessor() != nullptr) {
    getHfSumProcessor()->setupLuts(scale);
  }
}
 
void L1GlobalCaloTrigger::setChannelMask(const L1GctChannelMask* const mask) {
  m_inputChannelMask = mask;
  // Need to propagate the new mask to all the JetFinders
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->getJetFinderA()->setEnergySumMasks(mask);
    theJetLeafCards.at(i)->getJetFinderB()->setEnergySumMasks(mask);
    theJetLeafCards.at(i)->getJetFinderC()->setEnergySumMasks(mask);
  }
}
 
bool L1GlobalCaloTrigger::setupOk() const {
  bool result = true;
  result &= (m_inputChannelMask != nullptr);
  // EM Leaf Card
  for (int i = 0; i < N_EM_LEAF_CARDS; i++) {
    result &= theEmLeafCards.at(i)->setupOk();
  }
 
  // Jet Leaf cards
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    result &= theJetLeafCards.at(i)->setupOk();
  }
 
  // Jet Finders
  for (int i = 0; i < N_JET_LEAF_CARDS * 3; i++) {
    result &= theJetFinders.at(i)->setupOk();
  }
 
  // Wheel Cards
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    result &= theWheelJetFpgas.at(i)->setupOk();
  }
 
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    result &= theWheelEnergyFpgas.at(i)->setupOk();
  }
 
  // Electron Final Stage
  result &= theIsoEmFinalStage->setupOk();
  result &= theNonIsoEmFinalStage->setupOk();
 
  // Jet Final Stage
  result &= theJetFinalStage->setupOk();
 
  // Energy Final Stage
  result &= theEnergyFinalStage->setupOk();
 
  // All done.
  return result;
}
 
L1GctGlobalHfSumAlgos* L1GlobalCaloTrigger::getHfSumProcessor() const {
  L1GctGlobalHfSumAlgos* result = nullptr;
  if (theEnergyFinalStage != nullptr) {
    result = theEnergyFinalStage->getHfSumProcessor();
  }
  return result;
}
 
void L1GlobalCaloTrigger::setBxRange(const int firstBx, const int lastBx) {
  m_bxStart = firstBx;
  m_numOfBx = lastBx - firstBx + 1;
  m_bxRangeAuto = false;
}
void L1GlobalCaloTrigger::setBxRangeSymmetric(const int numOfBx) {
  m_bxStart = -numOfBx;
  m_numOfBx = 2 * numOfBx + 1;
  m_bxRangeAuto = false;
}
void L1GlobalCaloTrigger::setBxRangeAutomatic() {
  m_bxStart = 0;
  m_numOfBx = 1;
  m_bxRangeAuto = true;
}
 
void L1GlobalCaloTrigger::fillRegions(const vector<L1CaloRegion>& rgn) {
  // To enable multiple bunch crossing operation, we copy the input regions into a vector,
  // from which they will be extracted one bunch crossing at a time and sent to the processors
  vector<L1CaloRegion>::iterator itr = m_allInputRegions.end();
  m_allInputRegions.insert(itr, rgn.begin(), rgn.end());
}
 
void L1GlobalCaloTrigger::fillEmCands(const vector<L1CaloEmCand>& em) {
  // To enable multiple bunch crossing operation, we copy the input electrons into a vector,
  // from which they will be extracted one bunch crossing at a time and sent to the processors
  vector<L1CaloEmCand>::iterator itr = m_allInputEmCands.end();
  m_allInputEmCands.insert(itr, em.begin(), em.end());
}
 
void L1GlobalCaloTrigger::fillRegions(vector<L1CaloRegion>::iterator& rgn, const int bx) {
  while (rgn != m_allInputRegions.end() && rgn->bx() == bx) {
    setRegion(*rgn++);
  }
}
 
void L1GlobalCaloTrigger::fillEmCands(vector<L1CaloEmCand>::iterator& emc, const int bx) {
  while (emc != m_allInputEmCands.end() && emc->bx() == bx) {
    if (emc->isolated()) {
      setIsoEm(*emc);
    } else {
      setNonIsoEm(*emc);
    }
    emc++;
  }
}
 
void L1GlobalCaloTrigger::setRegion(const L1CaloRegion& region) {
  if (!m_inputChannelMask->regionMask(region.gctEta(), region.gctPhi())) {
    unsigned crate = region.rctCrate();
    // Find the relevant jetFinders
    static const unsigned NPHI = L1CaloRegionDetId::N_PHI / 2;
    unsigned prevphi = crate % NPHI;
    unsigned thisphi = (crate + 1) % NPHI;
    unsigned nextphi = (crate + 2) % NPHI;
 
    // Send the region to six jetFinders.
    theJetFinders.at(thisphi)->setInputRegion(region);
    theJetFinders.at(nextphi)->setInputRegion(region);
    theJetFinders.at(prevphi)->setInputRegion(region);
    theJetFinders.at(thisphi + NPHI)->setInputRegion(region);
    theJetFinders.at(nextphi + NPHI)->setInputRegion(region);
    theJetFinders.at(prevphi + NPHI)->setInputRegion(region);
  }
}
 
void L1GlobalCaloTrigger::setRegion(
    const unsigned et, const unsigned ieta, const unsigned iphi, const bool overFlow, const bool fineGrain) {
  //  L1CaloRegion temp = L1CaloRegion::makeRegionFromGctIndices(et, overFlow, fineGrain, false, false, ieta, iphi, 0);
  L1CaloRegion temp(et, overFlow, fineGrain, false, false, ieta, iphi, 0);
  setRegion(temp);
}
 
void L1GlobalCaloTrigger::setIsoEm(const L1CaloEmCand& em) {
  if (!m_inputChannelMask->emCrateMask(em.rctCrate()))
    theIsoElectronSorters.at(sorterNo(em))->setInputEmCand(em);
}
 
void L1GlobalCaloTrigger::setNonIsoEm(const L1CaloEmCand& em) {
  if (!m_inputChannelMask->emCrateMask(em.rctCrate()))
    theNonIsoElectronSorters.at(sorterNo(em))->setInputEmCand(em);
}
 
void L1GlobalCaloTrigger::print() {
  using edm::LogInfo;
  using std::endl;
 
  LogInfo("L1GlobalCaloTrigger") << "=== Global Calo Trigger ===" << endl;
  LogInfo("L1GlobalCaloTrigger") << "=== START DEBUG OUTPUT  ===" << endl;
 
  LogInfo("L1GlobalCaloTrigger") << endl;
  LogInfo("L1GlobalCaloTrigger") << "N Jet Leaf Cards " << theJetLeafCards.size() << endl;
  LogInfo("L1GlobalCaloTrigger") << "N Wheel Jet Fpgas " << theWheelJetFpgas.size() << endl;
  LogInfo("L1GlobalCaloTrigger") << "N Wheel Energy Fpgas " << theWheelEnergyFpgas.size() << endl;
  LogInfo("L1GlobalCaloTrigger") << "N Em Leaf Cards " << theEmLeafCards.size() << endl;
  LogInfo("L1GlobalCaloTrigger") << endl;
 
  for (unsigned i = 0; i < theJetLeafCards.size(); i++) {
    LogInfo("L1GlobalCaloTrigger") << "Jet Leaf Card " << i << " : " << theJetLeafCards.at(i) << endl;
    LogInfo("L1GlobalCaloTrigger") << (*theJetLeafCards.at(i));
  }
  LogInfo("L1GlobalCaloTrigger") << endl;
 
  for (unsigned i = 0; i < theWheelJetFpgas.size(); i++) {
    LogInfo("L1GlobalCaloTrigger") << "Wheel Jet FPGA " << i << " : " << theWheelJetFpgas.at(i) << endl;
    LogInfo("L1GlobalCaloTrigger") << (*theWheelJetFpgas.at(i));
  }
  LogInfo("L1GlobalCaloTrigger") << endl;
 
  for (unsigned i = 0; i < theWheelEnergyFpgas.size(); i++) {
    LogInfo("L1GlobalCaloTrigger") << "Wheel Energy FPGA " << i << " : " << theWheelEnergyFpgas.at(i) << endl;
    LogInfo("L1GlobalCaloTrigger") << (*theWheelEnergyFpgas.at(i));
  }
  LogInfo("L1GlobalCaloTrigger") << endl;
 
  LogInfo("L1GlobalCaloTrigger") << (*theJetFinalStage);
  LogInfo("L1GlobalCaloTrigger") << endl;
 
  LogInfo("L1GlobalCaloTrigger") << (*theEnergyFinalStage);
  LogInfo("L1GlobalCaloTrigger") << endl;
 
  for (unsigned i = 0; i < theEmLeafCards.size(); i++) {
    LogInfo("L1GlobalCaloTrigger") << ((i == 0) ? "Positive eta " : "Negative eta ");
    LogInfo("L1GlobalCaloTrigger") << "EM Leaf Card " << i << " : " << theEmLeafCards.at(i) << endl;
    LogInfo("L1GlobalCaloTrigger") << (*theEmLeafCards.at(i));
  }
  LogInfo("L1GlobalCaloTrigger") << endl;
 
  LogInfo("L1GlobalCaloTrigger") << (*theIsoEmFinalStage);
  LogInfo("L1GlobalCaloTrigger") << endl;
 
  LogInfo("L1GlobalCaloTrigger") << (*theNonIsoEmFinalStage);
 
  LogInfo("L1GlobalCaloTrigger") << "=== Global Calo Trigger ===" << endl;
  LogInfo("L1GlobalCaloTrigger") << "===  END DEBUG OUTPUT   ===" << endl;
}
 
// isolated EM outputs
L1GctEmCandCollection L1GlobalCaloTrigger::getIsoElectrons() const { return theIsoEmFinalStage->getOutputCands(); }
 
// non isolated EM outputs
L1GctEmCandCollection L1GlobalCaloTrigger::getNonIsoElectrons() const {
  return theNonIsoEmFinalStage->getOutputCands();
}
 
// central jet outputs to GT
L1GctJetCandCollection L1GlobalCaloTrigger::getCentralJets() const { return theJetFinalStage->getCentralJets(); }
 
// forward jet outputs to GT
L1GctJetCandCollection L1GlobalCaloTrigger::getForwardJets() const { return theJetFinalStage->getForwardJets(); }
 
// tau jet outputs to GT
L1GctJetCandCollection L1GlobalCaloTrigger::getTauJets() const { return theJetFinalStage->getTauJets(); }
 
L1GctInternJetDataCollection L1GlobalCaloTrigger::getInternalJets() const {
  L1GctInternJetDataCollection allJets, jfJets;
 
  // Loop over jetfinders, find the internal jets and add them to the list
  for (unsigned jf = 0; jf < theJetFinders.size(); jf++) {
    jfJets = theJetFinders.at(jf)->getInternalJets();
    allJets.insert(allJets.end(), jfJets.begin(), jfJets.end());
  }
 
  return allJets;
}
 
// total Et output
L1GctEtTotalCollection L1GlobalCaloTrigger::getEtSumCollection() const {
  L1GctEtTotalCollection result(m_numOfBx);
  int bx = m_bxStart;
  for (int i = 0; i < m_numOfBx; i++) {
    L1GctEtTotal temp(
        theEnergyFinalStage->getEtSumColl().at(i).value(), theEnergyFinalStage->getEtSumColl().at(i).overFlow(), bx++);
    result.at(i) = temp;
  }
  return result;
}
 
L1GctEtHadCollection L1GlobalCaloTrigger::getEtHadCollection() const {
  L1GctEtHadCollection result(m_numOfBx);
  int bx = m_bxStart;
  for (int i = 0; i < m_numOfBx; i++) {
    L1GctEtHad temp(
        theEnergyFinalStage->getEtHadColl().at(i).value(), theEnergyFinalStage->getEtHadColl().at(i).overFlow(), bx++);
    result.at(i) = temp;
  }
  return result;
}
 
L1GctEtMissCollection L1GlobalCaloTrigger::getEtMissCollection() const {
  L1GctEtMissCollection result(m_numOfBx);
  int bx = m_bxStart;
  for (int i = 0; i < m_numOfBx; i++) {
    L1GctEtMiss temp(theEnergyFinalStage->getEtMissColl().at(i).value(),
                     theEnergyFinalStage->getEtMissPhiColl().at(i).value(),
                     theEnergyFinalStage->getEtMissColl().at(i).overFlow(),
                     bx++);
    result.at(i) = temp;
  }
  return result;
}
 
L1GctHtMissCollection L1GlobalCaloTrigger::getHtMissCollection() const {
  L1GctHtMissCollection result(m_numOfBx);
  int bx = m_bxStart;
  for (int i = 0; i < m_numOfBx; i++) {
    L1GctHtMiss temp(theEnergyFinalStage->getHtMissColl().at(i).value(),
                     theEnergyFinalStage->getHtMissPhiColl().at(i).value(),
                     theEnergyFinalStage->getHtMissColl().at(i).overFlow(),
                     bx++);
    result.at(i) = temp;
  }
  return result;
}
 
L1GctInternEtSumCollection L1GlobalCaloTrigger::getInternalEtSums() const {
  L1GctInternEtSumCollection allSums, procSums;
 
  // Go through all the processor types that process et sums
  // JetFinders
  for (unsigned jf = 0; jf < theJetFinders.size(); jf++) {
    procSums = theJetFinders.at(jf)->getInternalEtSums();
    allSums.insert(allSums.end(), procSums.begin(), procSums.end());
  }
 
  // Jet Leaf cards
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    procSums = theJetLeafCards.at(i)->getInternalEtSums();
    allSums.insert(allSums.end(), procSums.begin(), procSums.end());
  }
 
  // Wheel Cards
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    procSums = theWheelEnergyFpgas.at(i)->getInternalEtSums();
    allSums.insert(allSums.end(), procSums.begin(), procSums.end());
  }
 
  return allSums;
}
 
L1GctInternHtMissCollection L1GlobalCaloTrigger::getInternalHtMiss() const {
  L1GctInternHtMissCollection allSums, procSums;
 
  // Go through all the processor types that process et sums
  // JetFinders
  for (unsigned jf = 0; jf < theJetFinders.size(); jf++) {
    procSums = theJetFinders.at(jf)->getInternalHtMiss();
    allSums.insert(allSums.end(), procSums.begin(), procSums.end());
  }
 
  // Jet Leaf cards
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    procSums = theJetLeafCards.at(i)->getInternalHtMiss();
    allSums.insert(allSums.end(), procSums.begin(), procSums.end());
  }
 
  // Wheel Cards
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    procSums = theWheelJetFpgas.at(i)->getInternalHtMiss();
    allSums.insert(allSums.end(), procSums.begin(), procSums.end());
  }
 
  return allSums;
}
 
L1GctHFBitCountsCollection L1GlobalCaloTrigger::getHFBitCountsCollection() const {
  L1GctHFBitCountsCollection result(m_numOfBx);
  if (getHfSumProcessor() != nullptr) {
    int bx = m_bxStart;
    for (int i = 0; i < m_numOfBx; i++) {
      L1GctHFBitCounts temp = L1GctHFBitCounts::fromGctEmulator(
          static_cast<int16_t>(bx),
          getHfSumProcessor()->hfSumsOutput(L1GctHfEtSumsLut::bitCountPosEtaRing1).at(i),
          getHfSumProcessor()->hfSumsOutput(L1GctHfEtSumsLut::bitCountNegEtaRing1).at(i),
          getHfSumProcessor()->hfSumsOutput(L1GctHfEtSumsLut::bitCountPosEtaRing2).at(i),
          getHfSumProcessor()->hfSumsOutput(L1GctHfEtSumsLut::bitCountNegEtaRing2).at(i));
      result.at(i) = temp;
      bx++;
    }
  }
  return result;
}
 
L1GctHFRingEtSumsCollection L1GlobalCaloTrigger::getHFRingEtSumsCollection() const {
  L1GctHFRingEtSumsCollection result(m_numOfBx);
  if (getHfSumProcessor() != nullptr) {
    int bx = m_bxStart;
    for (int i = 0; i < m_numOfBx; i++) {
      L1GctHFRingEtSums temp = L1GctHFRingEtSums::fromGctEmulator(
          static_cast<int16_t>(bx),
          getHfSumProcessor()->hfSumsOutput(L1GctHfEtSumsLut::etSumPosEtaRing1).at(i),
          getHfSumProcessor()->hfSumsOutput(L1GctHfEtSumsLut::etSumNegEtaRing1).at(i),
          getHfSumProcessor()->hfSumsOutput(L1GctHfEtSumsLut::etSumPosEtaRing2).at(i),
          getHfSumProcessor()->hfSumsOutput(L1GctHfEtSumsLut::etSumNegEtaRing2).at(i));
      result.at(i) = temp;
      bx++;
    }
  }
  return result;
}
 
void L1GlobalCaloTrigger::setVerbose() {
  // EM Leaf Card
  for (int i = 0; i < N_EM_LEAF_CARDS; i++) {
    theEmLeafCards.at(i)->setVerbose();
  }
 
  // Jet Leaf cards
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->setVerbose();
  }
 
  // Jet Finders
  for (int i = 0; i < N_JET_LEAF_CARDS * 3; i++) {
    theJetFinders.at(i)->setVerbose();
  }
 
  // Wheel Cards
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelJetFpgas.at(i)->setVerbose();
  }
 
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelEnergyFpgas.at(i)->setVerbose();
  }
 
  // Electron Final Stage
  theIsoEmFinalStage->setVerbose();
  theNonIsoEmFinalStage->setVerbose();
 
  // Jet Final Stage
  theJetFinalStage->setVerbose();
 
  // Energy Final Stage
  theEnergyFinalStage->setVerbose();
}
 
void L1GlobalCaloTrigger::setTerse() {
  // EM Leaf Card
  for (int i = 0; i < N_EM_LEAF_CARDS; i++) {
    theEmLeafCards.at(i)->setTerse();
  }
 
  // Jet Leaf cards
  for (int i = 0; i < N_JET_LEAF_CARDS; i++) {
    theJetLeafCards.at(i)->setTerse();
  }
 
  // Jet Finders
  for (int i = 0; i < N_JET_LEAF_CARDS * 3; i++) {
    theJetFinders.at(i)->setTerse();
  }
 
  // Wheel Cards
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelJetFpgas.at(i)->setTerse();
  }
 
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    theWheelEnergyFpgas.at(i)->setTerse();
  }
 
  // Electron Final Stage
  theIsoEmFinalStage->setTerse();
  theNonIsoEmFinalStage->setTerse();
 
  // Jet Final Stage
  theJetFinalStage->setTerse();
 
  // Energy Final Stage
  theEnergyFinalStage->setTerse();
}
 
/* PRIVATE METHODS */
 
// instantiate hardware/algorithms
void L1GlobalCaloTrigger::build(L1GctJetLeafCard::jetFinderType jfType, unsigned jetLeafMask) {
  // The first half of the jet leaf cards are at negative eta,
  // followed by positive eta
  // Jet Leaf cards
  if (jetLeafMask == 0) {
    for (int jlc = 0; jlc < N_JET_LEAF_CARDS; jlc++) {
      theJetLeafCards.at(jlc) = new L1GctJetLeafCard(jlc, jlc % 3, jfType);
      theJetFinders.at(3 * jlc) = theJetLeafCards.at(jlc)->getJetFinderA();
      theJetFinders.at(3 * jlc + 1) = theJetLeafCards.at(jlc)->getJetFinderB();
      theJetFinders.at(3 * jlc + 2) = theJetLeafCards.at(jlc)->getJetFinderC();
    }
  } else {
    // Setup for hardware testing with reduced number of leaf cards
    unsigned mask = jetLeafMask;
    for (int jlc = 0; jlc < N_JET_LEAF_CARDS; jlc++) {
      if ((mask & 1) == 0) {
        theJetLeafCards.at(jlc) = new L1GctJetLeafCard(jlc, jlc % 3, jfType);
      } else {
        theJetLeafCards.at(jlc) = new L1GctJetLeafCard(jlc, jlc % 3, L1GctJetLeafCard::nullJetFinder);
      }
      theJetFinders.at(3 * jlc) = theJetLeafCards.at(jlc)->getJetFinderA();
      theJetFinders.at(3 * jlc + 1) = theJetLeafCards.at(jlc)->getJetFinderB();
      theJetFinders.at(3 * jlc + 2) = theJetLeafCards.at(jlc)->getJetFinderC();
      mask = mask >> 1;
    }
  }
 
  //Link jet leaf cards together
  vector<L1GctJetLeafCard*> neighbours(2);
  for (int jlc = 0; jlc < N_JET_LEAF_CARDS / 2; jlc++) {
    // Define local constant for ease of typing
    static const int NL = N_JET_LEAF_CARDS / 2;
    int nlc = (jlc + 1) % NL;
    int mlc = (jlc + (NL - 1)) % NL;
    neighbours.at(0) = theJetLeafCards.at(mlc);
    neighbours.at(1) = theJetLeafCards.at(nlc);
    theJetLeafCards.at(jlc)->setNeighbourLeafCards(neighbours);
    neighbours.at(0) = theJetLeafCards.at(NL + mlc);
    neighbours.at(1) = theJetLeafCards.at(NL + nlc);
    theJetLeafCards.at(NL + jlc)->setNeighbourLeafCards(neighbours);
  }
 
  // EM leaf cards
  // Card 0 is positive eta, card 1 is negative eta
  for (int i = 0; i < N_EM_LEAF_CARDS; i++) {
    theEmLeafCards.at(i) = new L1GctEmLeafCard(i);
    theIsoElectronSorters.at(2 * i) = theEmLeafCards.at(i)->getIsoElectronSorterU1();
    theIsoElectronSorters.at(2 * i + 1) = theEmLeafCards.at(i)->getIsoElectronSorterU2();
    theNonIsoElectronSorters.at(2 * i) = theEmLeafCards.at(i)->getNonIsoElectronSorterU1();
    theNonIsoElectronSorters.at(2 * i + 1) = theEmLeafCards.at(i)->getNonIsoElectronSorterU2();
  }
 
  // Wheel Fpgas
  vector<L1GctJetLeafCard*> wheelJetLeafCards(3);
  vector<L1GctJetLeafCard*> wheelEnergyLeafCards(3);
 
  // The first wheel card is at negative eta,
  // the second one is at positive eta
  for (int i = 0; i < N_WHEEL_CARDS; i++) {
    for (int j = 0; j < 3; j++) {
      wheelJetLeafCards.at(j) = theJetLeafCards.at(i * 3 + j);
      wheelEnergyLeafCards.at(j) = theJetLeafCards.at(i * 3 + j);
    }
    theWheelJetFpgas.at(i) = new L1GctWheelJetFpga(i, wheelJetLeafCards);
    theWheelEnergyFpgas.at(i) = new L1GctWheelEnergyFpga(i, wheelEnergyLeafCards);
  }
 
  // Jet Final Stage
  theJetFinalStage = new L1GctJetFinalStage(theWheelJetFpgas);
 
  // Electron Final Sort
  theIsoEmFinalStage = new L1GctElectronFinalSort(true, theEmLeafCards.at(0), theEmLeafCards.at(1));
  theNonIsoEmFinalStage = new L1GctElectronFinalSort(false, theEmLeafCards.at(0), theEmLeafCards.at(1));
 
  // Global Energy Algos
  theEnergyFinalStage = new L1GctGlobalEnergyAlgos(theWheelEnergyFpgas, theWheelJetFpgas);
}
 
unsigned L1GlobalCaloTrigger::sorterNo(const L1CaloEmCand& em) const {
  unsigned crate = em.rctCrate();
  unsigned result = (((crate % 9) < 4) ? 1 : 0);
  if (crate >= 9)
    result += 2;
  if (crate >= 18)
    result = 0;
  return result;
}