TopoSteeringStructure.cxx File Reference

#include "L1TopoCoreSim/TopoSteeringStructure.h"
#include "L1TopoCommon/StringUtils.h"
#include "L1TopoCommon/Exception.h"
#include "L1TopoCommon/Types.h"
#include "L1TopoConfig/LayoutConstraints.h"
#include "L1TopoConfig/L1TopoMenu.h"
#include "L1TopoInterfaces/DecisionAlg.h"
#include "L1TopoInterfaces/SortingAlg.h"
#include "L1TopoInterfaces/CountingAlg.h"
#include "L1TopoCoreSim/Connector.h"
#include "L1TopoCoreSim/InputConnector.h"
#include "L1TopoCoreSim/SortingConnector.h"
#include "L1TopoCoreSim/DecisionConnector.h"
#include "L1TopoCoreSim/CountingConnector.h"
#include "L1TopoHardware/L1TopoHardware.h"
#include "TrigConfData/L1Menu.h"
#include <set>
#include <iostream>
#include <iomanip>
#include <memory>
#include <boost/lexical_cast.hpp>

Go to the source code of this file.

Functions
TCS::StatusCode TCS::TopoSteeringStructure::setupFromMenu	ATLAS_NOT_THREAD_SAFE (const TrigConf::L1Menu &l1menu, bool legacy, bool debug)
TCS::StatusCode TCS::TopoSteeringStructure::instantiateAlgorithms	ATLAS_NOT_THREAD_SAFE (bool debug)

Function Documentation

ATLAS_NOT_THREAD_SAFE() [1/2]

TCS::StatusCode TCS::TopoSteeringStructure::instantiateAlgorithms ATLAS_NOT_THREAD_SAFE

(

bool

debug

)

Definition at line 443 of file TopoSteeringStructure.cxx.

                                                                                 {
 
   for(TCS::Connector* conn: m_connectors) {
 
      if(conn->isInputConnector()) continue;
 
      // for each connector instantiate the algorithm and add to connector
      const std::string & alg = conn->algorithmName();
       
      // split into name and type
      std::string algType(alg, 0, alg.find('/'));
      std::string algName(alg, alg.find('/')+1);
 
      ConfigurableAlg * algInstance = TCS::AlgFactory::mutable_instance().algorithm(algName);
      if(algInstance==0) {
         if(debug)
            cout << "Instantiating " << alg << endl;
         algInstance = TCS::AlgFactory::mutable_instance().create(algType, algName);
      } else {
         if(algInstance->className() != algType) {
            TCS_EXCEPTION("L1 TopoSteering: duplicate algorithm names:  algorithm " << algName << " has already been instantiated but with different type");
         }
      }
      conn->setAlgorithm(algInstance);
   }
   return TCS::StatusCode::SUCCESS;
}

ATLAS_NOT_THREAD_SAFE() [2/2]

TCS::StatusCode TCS::TopoSteeringStructure::setupFromMenu ATLAS_NOT_THREAD_SAFE	(	const TrigConf::L1Menu &	l1menu,
		bool	legacy,
		bool	debug
	)

Definition at line 125 of file TopoSteeringStructure.cxx.

                                                                                                                    {
 
   if(debug)
      cout << "/***************************************************************************/" << endl
           << "           L1Topo steering structure: configuring from L1 Topo Menu          " << endl
           << "/***************************************************************************/" << endl;
 
   //   set<TCS::inputTOBType_t> neededInputs; // Stores inputs for DecisionConnectors
   vector<string> storedConn; // Stores decision connectors in order to avoid double counting
   vector<vector<string>> confAlgorithms; // Stores algorithm name/category that have been configured in L1Menu to be used for setting the parameters
   vector<string> confMultAlgorithms; // Stores algorithm names that have been configured in L1Menu to be used for setting the multiplicity thresholds
   // Loop over boards in L1Menu and skip the ones that are not TOPO. Use Run-2 boards if legacy flag is on
 
   string AvailableMultAlgs[] = { "eEmMultiplicity", 
                  "jEmMultiplicity", 
                  "eTauMultiplicity", 
                  "jTauMultiplicity",
                  "cTauMultiplicity",
                  "jJetMultiplicity",
                  "jLJetMultiplicity",
                  "gJetMultiplicity",
                  "gLJetMultiplicity",
                  "EnergyThreshold",
                                  "LArSaturation" };
 
   for (const string & boardName : l1menu.boardNames() ){
     
     auto & l1board = l1menu.board(boardName);
 
     if (l1board.type() != "TOPO") continue;
     if (l1board.legacy() != legacy) continue;
 
     // Access the connectors in the boards
 
     for (const string & connName : l1board.connectorNames() ){
 
       auto & l1conn = l1menu.connector(connName);
 
        // First configure decision algorithms
        // Look at all Topo algorithms in each connector, and get inputs from each algorithm to configure SortingConnectors
         if ( l1conn.connectorType() == TrigConf::L1Connector::ConnectorType::ELECTRICAL ) {
 
          for( size_t fpga : { 0, 1} ) {
              for( size_t clock : { 0, 1} ) {
                for( auto & tl : l1conn.triggerLines(fpga, clock) ) {
         
                  const string & tlName = tl.name();
                  auto & algo = l1menu.algorithmFromTriggerline(tlName);
              const string algoName = (legacy?"R2_"+algo.name():algo.name());
 
                  // One algorithm can have multiple trigger lines. Check the connector/algorithm has not been stored already
                  auto it = find(storedConn.begin(), storedConn.end(), algoName);
                  if (it == storedConn.end()) { // Algorithm/Connector does not exist: create and store it
 
                    storedConn.push_back(algoName);
                    vector<string> inputNames;
                    for( auto & input : algo.inputs() ) {
                        if( sortingConnector(input) == 0 ) { // if connector does not exist, create it
                          if(debug) 
                            cout << "L1TopoSteering: Decision algo( " << algo.name() << " ) input is not defined: " << input << ". Creating sortingConnector" << endl;
           
                          auto & sortAlgo = l1menu.algorithm(input, algo.category());
                          if(!(sortAlgo.type() == TrigConf::L1TopoAlgorithm::AlgorithmType::SORTING ) )
                            TCS_EXCEPTION("L1TopoSteering: Decision algo " << algo.name() << ") has as input " << input << " which is not associated to a sorting algorithm");
 
                          // Create connector
                          SortingConnector * sortConn = new SortingConnector(sortAlgo.inputs().at(0), sortAlgo.klass()+"/"+input, sortAlgo.outputs().at(0));
                          if(debug)
                            cout << "Adding sorting connector " << "[" << *sortConn << "]" << endl;
                          addSortingConnector( sortConn );
                          confAlgorithms.push_back({sortAlgo.name(), sortAlgo.category()});
 
                          } // if connector does not exist
         
                        inputNames.push_back(input);
 
                    } // loop over inputs
 
                    DecisionConnector * conn = new DecisionConnector(algoName, inputNames, algo.klass()+"/"+algoName, algo.outputs());
                    conn->m_decision.setNBits( algo.outputs().size() );
 
                    if(tl.name() != "UNDEF")
                        conn->m_triggers.push_back(tl);
 
                    if(debug)
                        cout << "Adding decision connector " << "[" << *conn << "]" << endl;
                    addDecisionConnector( conn );
                    confAlgorithms.push_back({algo.name(), algo.category()});
 
                  } else { // Connector already exists: look for it and add the trigger line
                    for(const auto & out : algo.outputs()){
                        auto c = m_outputLookup.find(out);
                        if (c != m_outputLookup.end()){
                          auto conn = c->second;
                          if(tl.name() != "UNDEF")
                          conn->m_triggers.push_back(tl);
                          break;
                        }
                    }
                  }
 
                } // Trigger Line
 
              } // Clock
 
          } // FPGA
 
         } else { // Configure optical connectors - multiplicity algorithms
 
          // In multiplicity boards all trigger lines are copied into the L1Menu::Connector 4 times (fpga 0,1 and clock 0,1)
          // Take (fpga, clock) = (0, 0)
 
            for( auto & tl : l1conn.triggerLines(0, 0) ) {
 
          const string & tlName = tl.name();
          auto & algo = l1menu.algorithmFromTriggerline(tlName);
 
              string algo_klass = algo.klass();
              if(algo_klass=="eEmVarMultiplicity") algo_klass="eEmMultiplicity"; // in sim, use the same multiplicity algo for fixed and variable thresholds
 
          //Zero Bias events can't be simulated.
          //The following line is set to simply ignore it in the simulation
          //ZeroBias Multiplicity algorithms are not added to confMultAlgorithms
          if ( algo_klass == "ZeroBias" ) continue;
 
              auto it = find(storedConn.begin(), storedConn.end(), algo.name());
          if (it == storedConn.end()) { // Algorithm/Connector does not exist: create and store it
 
        storedConn.push_back(algo.name());
        if(debug)
          cout << "L1TopoSteering: Multiplicity algo( " << algo.name() << " ) has as input " << algo.inputs().at(0) << endl;
             
        CountingConnector * conn = new CountingConnector(algo.name(), algo.inputs().at(0), algo_klass+"/"+algo.name(), algo.outputs().at(0));
        conn->m_count.setNBits( tl.nbits() );
        conn->m_count.setFirstBit( tl.startbit() );
        
        if(tl.name() != "UNDEF")
          conn->m_triggers.push_back(tl);
        
        if(debug)
          cout << "Adding count connector " << "[" << *conn << "]" << endl;
        addCountingConnector( conn );
        confMultAlgorithms.push_back( algo.name() );
          }
          
        } // Trigger Line
        
     } // Optical connectors - multiplicities
     
     } // Connector in l1board
     
   } // Board in l1menu
 
  
   if(debug)
     cout << "... building input connectors" << endl;
   for(const auto & sortConn : m_sortedLookup) {
     const string & in = sortConn.second->inputNames()[0]; // name of input
     
     if( m_inputLookup.count(in) > 0 ) continue; // InputConnector already exists
 
     InputConnector * conn = new InputConnector(in);
     m_connectors.push_back(conn);
     m_inputLookup[in] = conn;
     if(debug)
       cout << "Adding input connector " << "[" << *conn << "]" << endl;
   }
   for(const auto & countConn : m_countLookup) {
     const string & in = countConn.second->inputNames()[0]; // name of input
 
     if( m_inputLookup.count(in) > 0 ) continue; // InputConnector already exists
 
     InputConnector * conn = new InputConnector(in);
     m_connectors.push_back(conn);
     m_inputLookup[in] = conn;
     if(debug)
       cout << "Adding input connector " << "[" << *conn << "]" << endl;
   }
   
   // link the connector objects together
   TCS::StatusCode sc = linkConnectors();
 
   // instantiate the algorithms from the algorithm names in the connectors
   if(debug)
     cout << "... instantiating algorithms" << endl;
   sc &= instantiateAlgorithms(debug);
 
 
   // set algorithm parameters
   if(debug)
     cout << "... setting algorithm parameters" << endl;
 
   for ( auto & confAlgo : confAlgorithms ) {
     
      auto & l1algo = l1menu.algorithm(confAlgo.at(0), confAlgo.at(1));
      auto l1algoName = confAlgo.at(0);
      if (l1algo.type() == TrigConf::L1TopoAlgorithm::AlgorithmType::DECISION && legacy)
          l1algoName ="R2_"+confAlgo.at(0);
 
      if(debug)
         cout << "TopoSteeringStructure: Parameters for algorithm with name " << l1algoName << " going to be configured." << endl;
      ConfigurableAlg * alg = AlgFactory::mutable_instance().algorithm(l1algoName);
     
      if(alg->isDecisionAlg()){
         ( static_cast<DecisionAlg *>(alg) )->setNumberOutputBits(l1algo.outputs().size());
      }
      // create ParameterSpace for this algorithm
      auto ps = std::make_unique<ParameterSpace>(alg->name());
 
      for(auto & pe : l1algo.parameters()) {
     
         auto & pname = pe.name();
         uint32_t val = pe.value();
         uint32_t sel = pe.selection();
 
         if(debug)
       cout << "Algo Name: " << l1algoName << " parameter " << ": " << setw(20) << left << pname << " value = " << setw(3) << left << val << " (selection " << sel << ")" << endl;
         ps->addParameter( pname, val, sel);
     
      }
       
      for(auto & gen : l1algo.generics().getKeys()) {
     
      auto pe = l1algo.generics().getObject(gen);
         string pname = gen;
         uint32_t val = interpretGenericParam(pe.getAttribute("value"));
         if (pname == "NumResultBits"){
             if(val != l1algo.outputs().size()) {
                TCS_EXCEPTION("Algorithm " << pname << " parameter OutputBits (" << val << ") is different from output size (" << l1algo.outputs().size() << ")");
             }
          continue; // ignore this, because it is defined through the output list
         }
         if(debug)
             cout << " fixed parameter : " << setw(20) << left << pname << " value = " << setw(3) << left << val << endl;
         ps->addParameter( pname, val );
     
      }
       
       
      if(debug)
         cout << " (setting parameters)";
      alg->setParameters( *ps );
       
      if(debug)
         cout << " --> (parameters set)";
       
     if(debug)
       cout << " --> (parameters stored)" << endl;
   } // Set parameters for Sorting/Decision algorithms
 
   // // set thresholds for multiplicity algorithms
   for ( auto & multAlgo : confMultAlgorithms ) {
 
      auto & l1algo = l1menu.algorithm(multAlgo, "MULTTOPO");
 
      //Zero Bias events can't be simulated.
      //The following line is set to simply ignore it in the simulation
      if ( l1algo.klass() == "ZeroBias" ) continue;
 
      ConfigurableAlg * alg = AlgFactory::mutable_instance().algorithm(l1algo.name());
 
      // Get L1Threshold object and pass it to CountingAlg, from where it will be propagated to and decoded in each algorithm
      // The output of each algorithm and the threshold name is the same - use output name to retrieve L1Threshold object
      auto & l1thr = l1menu.threshold( l1algo.outputs().at(0) );
      auto pCountAlg = dynamic_cast<CountingAlg *>(alg);
      if (not pCountAlg) continue;
      pCountAlg->setThreshold(l1thr);
 
   } // Set thresholds for multiplicity algorithms
   
   
   m_isConfigured = true;
   
   if(debug)
     cout << "... L1TopoSteering successfully configured" << endl;
   
   return sc;
}