Project CMSSW displayed by LXR CMSSW_14_1_X_2024-07-25-2300/​L1Trigger/​GlobalCaloTrigger/​src/​L1GlobalCaloTrigger.cc	Source navigation Diff markup Identifier search General search
	Version: CMSSW_14_1_X_2024-07-25-2300 CMSSW_14_1_X_2024-07-24-2300 CMSSW_14_1_X_2024-07-23-2300 CMSSW_14_1_X_2024-07-22-2300 CMSSW_14_1_X_2024-07-21-2300 CMSSW_14_1_X_2024-07-19-2300 Revert   Change

File indexing completed on 2024-04-06 12:19:53

 #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);
 }
 
 /// GCT Destructor
 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();
 }
 
 ///=================================================================================================
 ///
 /// Methods to reset all processors and process an event (consisting of multiple bunch crossings)
 ///
 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()) {
     /// Sort the input data by bunch crossing number
     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++;
     }
   }
 }
 
 /// Sort the input data by bunch crossing number
 void L1GlobalCaloTrigger::sortInputData() {
   std::sort(m_allInputEmCands.begin(), m_allInputEmCands.end(), emcBxComparator);
   std::sort(m_allInputRegions.begin(), m_allInputRegions.end(), rgnBxComparator);
 }
 
 /// Setup bunch crossing range (depending on input data)
 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);
 }
 
 /// Partial reset for a new bunch crossing
 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);
 }
 
 /// Process a new bunch crossing
 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();
 }
 
 ///=================================================================================================
 /// Configuration options for the GCT
 ///
 /// setup the Jet Finder parameters
 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);
   }
 }
 
 /// setup the Jet Calibration Lut
 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);
   }
 }
 
 /// Setup the tau algorithm parameters
 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);
   }
 }
 
 /// setup scale for missing Ht
 void L1GlobalCaloTrigger::setHtMissScale(const L1CaloEtScale* const scale) {
   if (theEnergyFinalStage != nullptr) {
     theEnergyFinalStage->setHtMissScale(scale);
   }
 }
 
 /// setup Hf sum LUTs
 void L1GlobalCaloTrigger::setupHfSumLuts(const L1CaloEtScale* const scale) {
   if (getHfSumProcessor() != nullptr) {
     getHfSumProcessor()->setupLuts(scale);
   }
 }
 
 /// setup the input channel mask
 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);
   }
 }
 
 /// check we have done all the setup
 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;
 }
 
 /// provide access to hf sum processor
 L1GctGlobalHfSumAlgos* L1GlobalCaloTrigger::getHfSumProcessor() const {
   L1GctGlobalHfSumAlgos* result = nullptr;
   if (theEnergyFinalStage != nullptr) {
     result = theEnergyFinalStage->getHfSumProcessor();
   }
   return result;
 }
 
 /// setup the bunch crossing range to be processed
 /// process crossings from (firstBx) to (lastBx)
 void L1GlobalCaloTrigger::setBxRange(const int firstBx, const int lastBx) {
   m_bxStart = firstBx;
   m_numOfBx = lastBx - firstBx + 1;
   m_bxRangeAuto = false;
 }
 /// process crossings from (-numOfBx) to (numOfBx)
 void L1GlobalCaloTrigger::setBxRangeSymmetric(const int numOfBx) {
   m_bxStart = -numOfBx;
   m_numOfBx = 2 * numOfBx + 1;
   m_bxRangeAuto = false;
 }
 /// process all crossings present in the input (and only those crossings)
 void L1GlobalCaloTrigger::setBxRangeAutomatic() {
   m_bxStart = 0;
   m_numOfBx = 1;
   m_bxRangeAuto = true;
 }
 
 ///=================================================================================================
 /// Input data set methods
 ///
 /// Use the following two methods for full emulator operation
 /// set jet regions from the RCT at the input to be processed
 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());
 }
 
 /// set electrons from the RCT at the input to be processed
 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());
 }
 
 /// Private method to send one bunch crossing's worth of regions to the processors
 void L1GlobalCaloTrigger::fillRegions(vector<L1CaloRegion>::iterator& rgn, const int bx) {
   while (rgn != m_allInputRegions.end() && rgn->bx() == bx) {
     setRegion(*rgn++);
   }
 }
 
 /// Private method to send one bunch crossing's worth of electrons to the processors
 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++;
   }
 }
 
 /// Set a jet region at the input to be processed
 /// Called from fillRegions() above - also available to be called directly
 /// (but the "user" has to take care of any multiple bunch crossing issues)
 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);
   }
 }
 
 /// Construct a jet region and set it at the input to be processed.
 /// For testing/debugging only.
 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);
 }
 
 /// Set an isolated EM candidate to be processed
 /// Called from fillEmCands() above - also available to be called directly
 /// (but the "user" has to take care of any multiple bunch crossing issues)
 void L1GlobalCaloTrigger::setIsoEm(const L1CaloEmCand& em) {
   if (!m_inputChannelMask->emCrateMask(em.rctCrate()))
     theIsoElectronSorters.at(sorterNo(em))->setInputEmCand(em);
 }
 
 /// Set a non-isolated EM candidate to be processed
 /// Called from fillEmCands() above - also available to be called directly
 /// (but the "user" has to take care of any multiple bunch crossing issues)
 void L1GlobalCaloTrigger::setNonIsoEm(const L1CaloEmCand& em) {
   if (!m_inputChannelMask->emCrateMask(em.rctCrate()))
     theNonIsoElectronSorters.at(sorterNo(em))->setInputEmCand(em);
 }
 
 ///=================================================================================================
 /// Print method
 ///
 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;
 }
 
 ///=================================================================================================
 /// Output data get methods
 ///
 // 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(); }
 
 /// all jets from jetfinders in raw format
 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;
 }
 
 /// control output messages
 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);
 }
 
 /// ordering of the electron sorters to give the correct
 /// priority to the candidates in the final sort
 /// The priority ordering is:
 ///    crates  4 - 8 : priority 0 (highest)
 ///    crates  0 - 3 : priority 1
 ///    crates 13 -17 : priority 2
 ///    crates  9 -12 : priority 3 (lowest)
 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;
 }

This page was automatically generated by the 2.2.1 LXR engine. The LXR team