LVL1::eFEXSysSim Class Reference

#include <eFEXSysSim.h>

Inheritance diagram for LVL1::eFEXSysSim:

Collaboration diagram for LVL1::eFEXSysSim:

Public Member Functions
	eFEXSysSim (const std::string &type, const std::string &name, const IInterface *parent)
	Constructors.
eFEXSysSim &&	operator= (const eFEXSysSim &)=delete
	Destructor.
virtual StatusCode	initialize ()
	standard Athena-Algorithm method
virtual StatusCode	execute (eFEXOutputCollection *inputOutputCollection)
virtual int	calcTowerID (int eta, int phi, int mod) const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Static Public Member Functions
static const InterfaceID &	interfaceID ()

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed.

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
StatusCode	StoreTauTOBs (std::map< int, std::vector< std::unique_ptr< eFEXtauTOB > > > &allTauTobObjects, SG::WriteHandleKey< xAOD::eFexTauRoIContainer > &eFexTauxTOBOutKey, SG::WriteHandleKey< xAOD::eFexTauRoIContainer > &eFexTauOutKey)
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Static Private Member Functions
template<class TOBObjectClass>
static bool	TOBetSort (const TOBObjectClass &i, const TOBObjectClass &j, bool isTau)
	Internal data.

Private Attributes
ToolHandle< eFEXSim >	m_eFEXSimTool {this, "eFEXSimTool", "LVL1::eFEXSim", "Tool that creates the eFEX Simulation"}
ToolHandle< eFEXFillEDM >	m_eFEXFillEDMTool {this, "eFEXFillEDMTool", "LVL1::eFEXFillEDM", "Tool to fille eFEX EDMs"}
SG::ReadHandleKey< LVL1::eTowerContainer >	m_eTowerContainerSGKey {this, "MyETowers", "eTowerContainer", "Input container for eTowers"}
SG::ReadHandleKey< TrigConf::L1Menu >	m_l1MenuKey {this, "L1TriggerMenu", "DetectorStore+L1TriggerMenu","Name of the L1Menu object to read configuration from"}
SG::WriteHandleKey< xAOD::eFexEMRoIContainer >	m_eFexOutKey {this,"Key_eFexEMOutputContainer","L1_eEMRoI","Output eFexEM TOB container"}
SG::WriteHandleKey< xAOD::eFexEMRoIContainer >	m_eFexEMxTOBOutKey {this,"Key_eFexEMxTOBOutputContainer","L1_eEMxRoI","Output eFexEM xTOB container"}
SG::WriteHandleKey< xAOD::eFexTauRoIContainer >	m_eFexTauActiveOutKey {this,"Key_eFexTauOutputContainer","L1_eTauRoI","Output eFexTau active (BDT/heuristic) algorithm TOB container"}
SG::WriteHandleKey< xAOD::eFexTauRoIContainer >	m_eFexTauActivexTOBOutKey {this,"Key_eFexTauxTOBOutputContainer","L1_eTauxRoI","Output eFexTau active (BDT/heuristic) algorithm xTOB container"}
SG::WriteHandleKey< xAOD::eFexTauRoIContainer >	m_eFexTauAltOutKey {this,"Key_eFexAltTauOutputContainer","","Output eFexTau alternative (BDT/heuristic) algorithm TOB container"}
SG::WriteHandleKey< xAOD::eFexTauRoIContainer >	m_eFexTauAltxTOBOutKey {this,"Key_eFexAltTauxTOBOutputContainer","","Output eFexTau alternative (BDT/heuristic) algorithm xTOB container"}
ToolHandle< eFEXFPGATowerIdProvider >	m_eFEXFPGATowerIdProviderTool {this, "eFEXFPGATowerIdProviderTool", "LVL1::eFEXFPGATowerIdProvider", "Tool that provides tower-FPGA mapping"}
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 42 of file eFEXSysSim.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

eFEXSysSim()

LVL1::eFEXSysSim::eFEXSysSim	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent )

Constructors.

Definition at line 37 of file eFEXSysSim.cxx.

                                                                                            :
    AthAlgTool(type,name,parent)
  {
    declareInterface<eFEXSysSim>(this);
 
  }

Member Function Documentation

calcTowerID()

int LVL1::eFEXSysSim::calcTowerID ( int eta, int phi, int mod ) const

virtual

Definition at line 74 of file eFEXSysSim.cxx.

                                                             {
 
    return ((64*eta) + phi + mod);
  }

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & hndl, const SG::VarHandleKeyType &  )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T, V, H > & t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                     {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode LVL1::eFEXSysSim::execute ( eFEXOutputCollection * inputOutputCollection )

virtual

Definition at line 79 of file eFEXSysSim.cxx.

                                                                            {
 
    SG::ReadHandle<LVL1::eTowerContainer> this_eTowerContainer(m_eTowerContainerSGKey/*,ctx*/);
    if(!this_eTowerContainer.isValid()){
      ATH_MSG_FATAL("Could not retrieve eTowerContainer " << m_eTowerContainerSGKey.key());
      return StatusCode::FAILURE;
    }
 
    // remove TOBs of the previous events from the map
    std::map<int, std::vector<std::unique_ptr<eFEXegTOB>> >  allEmTobObjects;
    std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> > allTauHeuristicTobObjects;
    std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> > allTauBDTTobObjects;
 
    // do mapping with preloaded csv file if it is available
    if (m_eFEXFPGATowerIdProviderTool->ifhaveinputfile()) {
      int tmp_eTowersIDs_subset_eFEX[10][18];
      for (int i_efex{ 0 }; i_efex < 24; i_efex++) {
          ATH_CHECK(m_eFEXFPGATowerIdProviderTool->getRankedTowerIDineFEX(i_efex, tmp_eTowersIDs_subset_eFEX));
          m_eFEXSimTool->init(i_efex);
          ATH_CHECK(m_eFEXSimTool->NewExecute(tmp_eTowersIDs_subset_eFEX, inputOutputCollection));
          // Get TOBs from this eFEX
          allEmTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXegTOB>> >::value_type(i_efex, m_eFEXSimTool->getEmTOBs() ));
          allTauHeuristicTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> >::value_type(i_efex, m_eFEXSimTool->getTauHeuristicTOBs() ));
          allTauBDTTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> >::value_type(i_efex, m_eFEXSimTool->getTauBDTTOBs() ));
          m_eFEXSimTool->reset();
      }
    } else {
    // We need to split the towers into 3 blocks in eta and 8 blocks in phi.
 
    // boundaries in eta: -2.5, -0.8, 0.8, 2.5
    // REAL boundaries in eta (overlaps must occur for sliding window algorithms!): -2.5, -0.7, -0.9, 0.9, 0.7, 2.5
    // Written explicitly:
    // -2.5 -> -0.7
    // -0.9 -> 0.9
    // 0.7 -> 2.5
 
    // boundaries in phi: 0.2, 1.0, 1.8, 2.6, 3.4, 4.2, 5.0, 5.8
    // Written explicitly with REAL boundaries in phi (overlaps must occur for sliding window algorithms!)
    // 0.1 -> 1.1
    // 0.9 -> 1.9
    // 1.7 -> 2.7
    // 2.5 -> 3.5
    // 3.3 -> 4.3
    // 4.1 -> 5.1
    // 4.9 -> 5.9
    // 5.7 -> 0.3
 
    // C-SIDE NEGATIVE EFEXs
    // DO THE LEFT-MOST (NEGATIVE ETA) EFEXs FIRST
    int fexcounter = 0;
    // Example values for eFEX 0
    //id_modifier + phi + (64 * eta)
    int emecEta = 24; int emecPhi = 1; int emecMod = 500000;
    int initialEMEC = calcTowerID(emecEta,emecPhi,emecMod); //501537;
    int transEta = 14; int transPhi = 1; int transMod = 300000;
    int initialTRANS = calcTowerID(transEta,transPhi,transMod); //300897;
    int embEta = 13; int embPhi = 1; int embMod = 100000;
    int initialEMB = calcTowerID(embEta,embPhi,embMod); //100833;
    int eFEXa = 0;
 
    for (int thisEFEX=eFEXa; thisEFEX<=21; thisEFEX+=3){
      
      if(fexcounter > 0){ initialEMEC += 8; initialTRANS += 8; initialEMB += 8; } // TODO // SOMEHOW REMOVE HARD-CODING?
 
      // decide which subset of towers (and therefore supercells) should go to the eFEX
      std::map<int,eTower> tmp_eTowersColl_subset;
      
      // let's try doing this with an array initially just containing tower IDs.
      int tmp_eTowersIDs_subset [10][18];
      
      int rows = sizeof tmp_eTowersIDs_subset / sizeof tmp_eTowersIDs_subset[0];
      int cols = sizeof tmp_eTowersIDs_subset[0] / sizeof tmp_eTowersIDs_subset[0][0];
 
      // set the EMEC part
      for(int thisCol=0; thisCol<10; thisCol++){
    for(int thisRow=0; thisRow<rows; thisRow++){
 
      int towerid = initialEMEC - (thisCol * 64) + thisRow;
 
      if( (thisEFEX == 21) && (thisRow >= 7)){ towerid -= 64; };
      
      tmp_eTowersIDs_subset[thisRow][thisCol] = towerid;
      tmp_eTowersColl_subset.insert( std::map<int, eTower>::value_type(towerid,  *(this_eTowerContainer->findTower(towerid))));
 
    }
      }
 
      // set the TRANS part
      for(int thisRow = 0; thisRow < rows; thisRow++){
 
    int towerid = initialTRANS + thisRow;
 
    if( (thisEFEX == 21) && (thisRow >= 7)){ towerid -= 64; };
 
    tmp_eTowersIDs_subset[thisRow][10] = towerid;
    tmp_eTowersColl_subset.insert( std::map<int, eTower>::value_type(towerid,  *(this_eTowerContainer->findTower(towerid))));
 
      }
 
      // set the EMB part
      for(int thisCol = 11; thisCol < cols; thisCol++){
    for(int thisRow=0; thisRow<rows; thisRow++){
 
      int towerid = initialEMB - ( (thisCol-11) * 64) + thisRow;
 
      if( (thisEFEX == 21) && (thisRow >= 7)){ towerid -= 64; };
 
      tmp_eTowersIDs_subset[thisRow][thisCol] = towerid;
      tmp_eTowersColl_subset.insert( std::map<int, eTower>::value_type(towerid,  *(this_eTowerContainer->findTower(towerid))));
 
    }
      }
      
 
      if(false){
    ATH_MSG_DEBUG("CONTENTS OF eFEX " << thisEFEX << " :");
        for (int thisRow=rows-1; thisRow>=0; thisRow--){
          for (int thisCol=0; thisCol<cols; thisCol++){
            int tmptowerid = tmp_eTowersIDs_subset[thisRow][thisCol];
            const float tmptowereta = this_eTowerContainer->findTower(tmptowerid)->eta();
            const float tmptowerphi = this_eTowerContainer->findTower(tmptowerid)->phi();
            if(thisCol != cols-1){ ATH_MSG_DEBUG("|  " << tmptowerid << "([" << tmptowerphi << "][" << tmptowereta << "])  "); }
            else { ATH_MSG_DEBUG("|  " << tmptowerid << "([" << tmptowereta << "][" << tmptowerphi << "])  |"); }
          }
        }
      }
 
 
      m_eFEXSimTool->init(thisEFEX);
      ATH_CHECK(m_eFEXSimTool->NewExecute(tmp_eTowersIDs_subset, inputOutputCollection));
      allEmTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXegTOB>> >::value_type(thisEFEX, (m_eFEXSimTool->getEmTOBs() ) ));
      allTauHeuristicTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> >::value_type(thisEFEX, (m_eFEXSimTool->getTauHeuristicTOBs() ) ));
      allTauBDTTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> >::value_type(thisEFEX, (m_eFEXSimTool->getTauBDTTOBs() ) ));
      m_eFEXSimTool->reset();
 
      fexcounter++;
    }
 
    // CENTRAL EFEXs
    // DO THE CENTRAL EFEXs SECOND
    fexcounter = 0;
    int embnegEta = 8; int embnegPhi = 1; int embnegMod = 100000;
    int initialEMB_neg = calcTowerID(embnegEta,embnegPhi,embnegMod); //100513;
    int embposEta = 0; int embposPhi = 1; int embposMod = 200000;
    int initialEMB_pos = calcTowerID(embposEta,embposPhi,embposMod); //200001;
    int eFEXb = 1;
 
    for (int thisEFEX=eFEXb; thisEFEX<=22; thisEFEX+=3){
 
      if(fexcounter > 0){  initialEMB_neg += 8; initialEMB_pos += 8; }
      
      // decide which subset of towers (and therefore supercells) should go to the eFEX
      std::map<int,eTower> tmp_eTowersColl_subset;
      
      // doing this with an array initially just containing tower IDs.
      int tmp_eTowersIDs_subset [10][18];
 
      int rows = sizeof tmp_eTowersIDs_subset / sizeof tmp_eTowersIDs_subset[0];
      int cols = sizeof tmp_eTowersIDs_subset[0] / sizeof tmp_eTowersIDs_subset[0][0];
 
      // set the EMB part
      for(int thisCol = 0; thisCol < cols; thisCol++){
        for(int thisRow=0; thisRow<rows; thisRow++){
      int towerid = -1;
 
      int tmp_initEMB = initialEMB_neg;
 
      if(thisCol < 9){
        towerid = tmp_initEMB - ( (thisCol) * 64) + thisRow;
      }
      else{
        tmp_initEMB = initialEMB_pos;
        towerid = tmp_initEMB + ( (thisCol-9) * 64) + thisRow;
      }
 
      if( (thisEFEX == 22) && (thisRow >= 7)){ towerid -= 64; };
 
          tmp_eTowersIDs_subset[thisRow][thisCol] = towerid;
 
          tmp_eTowersColl_subset.insert( std::map<int, eTower>::value_type(towerid,  *(this_eTowerContainer->findTower(towerid))));
        }
      }
      
 
      if(false){
          ATH_MSG_DEBUG("CONTENTS OF eFEX " << thisEFEX << " :");
        for (int thisRow=rows-1; thisRow>=0; thisRow--){
          for (int thisCol=0; thisCol<cols; thisCol++){
            int tmptowerid = tmp_eTowersIDs_subset[thisRow][thisCol];
            const float tmptowereta = this_eTowerContainer->findTower(tmptowerid)->eta();
            const float tmptowerphi = this_eTowerContainer->findTower(tmptowerid)->phi();
            if(thisCol != cols-1){ ATH_MSG_DEBUG("|  " << tmptowerid << "([" << tmptowereta << "][" << tmptowerphi << "])  "); }
            else { ATH_MSG_DEBUG("|  " << tmptowerid << "([" << tmptowereta << "][" << tmptowerphi << "])  |"); }
          }
        }
      }
 
      //tool use instead
      m_eFEXSimTool->init(thisEFEX);
      ATH_CHECK(m_eFEXSimTool->NewExecute(tmp_eTowersIDs_subset, inputOutputCollection));
      allEmTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXegTOB>> >::value_type(thisEFEX, (m_eFEXSimTool->getEmTOBs() ) ));
      allTauHeuristicTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> >::value_type(thisEFEX, (m_eFEXSimTool->getTauHeuristicTOBs() ) ));
      allTauBDTTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> >::value_type(thisEFEX, (m_eFEXSimTool->getTauBDTTOBs() ) ));
      m_eFEXSimTool->reset();
 
      fexcounter++;
    }
 
    // POSITIVE EFEXs
    // LET'S DO THE RIGHT-MOST (POSTITIVE ETA) EFEXs THIRD
    fexcounter = 0;
    // Example values for eFEX 0
    emecEta = 15; emecPhi = 1; emecMod = 600000;
    initialEMEC = calcTowerID(emecEta,emecPhi,emecMod); //600961;
    transEta = 14; transPhi = 1; transMod = 400000;
    initialTRANS = calcTowerID(transEta,transPhi,transMod); //400897;
    embEta = 7; embPhi = 1; embMod = 200000;
    initialEMB = calcTowerID(embEta,embPhi,embMod); //200449;
    int eFEXc = 2;
 
    for (int thisEFEX=eFEXc; thisEFEX<=23; thisEFEX+=3){
 
      if(fexcounter > 0){ initialEMEC += 8; initialTRANS += 8; initialEMB += 8; }
 
      // decide which subset of towers (and therefore supercells) should go to the eFEX
      std::map<int,eTower> tmp_eTowersColl_subset;
 
      // doing this with an array initially just containing tower IDs.
      int tmp_eTowersIDs_subset [10][18];
 
      int rows = sizeof tmp_eTowersIDs_subset / sizeof tmp_eTowersIDs_subset[0];
      int cols = sizeof tmp_eTowersIDs_subset[0] / sizeof tmp_eTowersIDs_subset[0][0];
 
      // set the EMB part
      for(int thisCol = 0; thisCol < 7; thisCol++){
        for(int thisRow=0; thisRow<rows; thisRow++){
          int towerid = initialEMB + ( (thisCol) * 64) + thisRow;
 
      if( (thisEFEX == 23) && (thisRow >= 7)){ towerid -= 64; };
 
          tmp_eTowersIDs_subset[thisRow][thisCol] = towerid;
          tmp_eTowersColl_subset.insert( std::map<int, eTower>::value_type(towerid,  *(this_eTowerContainer->findTower(towerid))));
        }
      }
      // set the TRANS part
      for(int thisRow = 0; thisRow < rows; thisRow++){
        int towerid = initialTRANS + thisRow;
 
    if( (thisEFEX == 23) && (thisRow >= 7)){ towerid -= 64; };
 
        tmp_eTowersIDs_subset[thisRow][7] = towerid;
        tmp_eTowersColl_subset.insert( std::map<int, eTower>::value_type(towerid,  *(this_eTowerContainer->findTower(towerid))));
      }
      // set the EMEC part
      for(int thisCol=8; thisCol<cols; thisCol++){
        for(int thisRow=0; thisRow<rows; thisRow++){
          int towerid = initialEMEC + ( (thisCol-8) * 64) + thisRow;
 
      if( (thisEFEX == 23) && (thisRow >= 7)){ towerid -= 64; };
 
          tmp_eTowersIDs_subset[thisRow][thisCol] = towerid;
          tmp_eTowersColl_subset.insert( std::map<int, eTower>::value_type(towerid,  *(this_eTowerContainer->findTower(towerid))));
        }
      }
 
      // Debug printout
      if(false){
    ATH_MSG_DEBUG("CONTENTS OF eFEX " << thisEFEX << " :");
    for (int thisRow=rows-1; thisRow>=0; thisRow--){
      for (int thisCol=0; thisCol<cols; thisCol++){
        int tmptowerid = tmp_eTowersIDs_subset[thisRow][thisCol];
        const float tmptowereta = this_eTowerContainer->findTower(tmptowerid)->eta();
        const float tmptowerphi = this_eTowerContainer->findTower(tmptowerid)->phi();
        if(thisCol != cols-1){ ATH_MSG_DEBUG("|  " << tmptowerid << "([" << tmptowereta << "][" << tmptowerphi << "])  "); }
        else { ATH_MSG_DEBUG("|  " << tmptowerid << "([" << tmptowereta << "][" << tmptowerphi << "])  |"); }
      }
    }
      }
 
      //tool use instead
      m_eFEXSimTool->init(thisEFEX);
      ATH_CHECK(m_eFEXSimTool->NewExecute(tmp_eTowersIDs_subset, inputOutputCollection));
      allEmTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXegTOB>> >::value_type(thisEFEX, (m_eFEXSimTool->getEmTOBs() ) ));
      allTauHeuristicTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> >::value_type(thisEFEX, (m_eFEXSimTool->getTauHeuristicTOBs() ) ));
      allTauBDTTobObjects.insert( std::map<int, std::vector<std::unique_ptr<eFEXtauTOB>> >::value_type(thisEFEX, (m_eFEXSimTool->getTauBDTTOBs() ) ));
      m_eFEXSimTool->reset();
 
      fexcounter++;
    }
    
    }//close the non-csv loop over eFEXes
 
    // EM TOBs and xTOBs
 
    // TOB Containers
    auto eContainer = std::make_unique<xAOD::eFexEMRoIContainer> ();
    auto eAuxContainer = std::make_unique<xAOD::eFexEMRoIAuxContainer> ();
    eContainer->setStore(eAuxContainer.get());
 
    // xTOB Containers
    auto xeContainer = std::make_unique<xAOD::eFexEMRoIContainer> ();
    auto xeAuxContainer = std::make_unique<xAOD::eFexEMRoIAuxContainer> ();
    xeContainer->setStore(xeAuxContainer.get());
 
    // iterate over all Em Tobs and fill xTOB EDM with them
    for( auto const& [efex, tobObjects] : allEmTobObjects ){
      for(auto &tobObject : tobObjects){
        m_eFEXFillEDMTool->fillEmEDM(xeContainer, efex, tobObject, true);
      }
    }
 
    // Form list of TOBs, sorted and truncated from each eFEX
    // Vector to store sorted TOBs from all eFEXes
    // Note that this step means moving TOBs from the all objects vector, which is why we do it last
    std::vector<std::unique_ptr<eFEXegTOB>> emTOBs;
    // Loop through eFEXes and sort TOBs from each
    auto iter = allEmTobObjects.begin();
    while (iter != allEmTobObjects.end()) {
        std::vector<std::unique_ptr<eFEXegTOB>> tobsSort = std::move(iter->second);
        // sort tobs by their et (last 12 bits of the 32 bit tob word)
        std::sort (tobsSort.begin(), tobsSort.end(), std::bind(TOBetSort<std::unique_ptr<eFEXegTOB>>, std::placeholders::_1, std::placeholders::_2, false));
        // Truncate at 6 TOBs per eFEX
        if (tobsSort.size() > 6) tobsSort.resize(6);
        // Append to system TOB list
        for (unsigned int t = 0; t < tobsSort.size(); ++t) emTOBs.push_back(std::move(tobsSort[t]));
        // Next eFEX
        ++iter;
    }
 
    // iterate over sorted eFEX EM TOBs and fill TOB EDM with them
    for(auto &tobObject : emTOBs){
        int efex = tobObject->geteFEXID();
        m_eFEXFillEDMTool->fillEmEDM(eContainer, efex, tobObject);
    }
 
    // Match xTOBs to TOBs and set isTOB flags if matched
    matchTOBs(eContainer, xeContainer);
   
    // Record EDMs in StoreGate   
    SG::WriteHandle<xAOD::eFexEMRoIContainer> outputeFexEMxTOBHandle(m_eFexEMxTOBOutKey/*, ctx*/);
    ATH_MSG_DEBUG("  write: " << outputeFexEMxTOBHandle.key() << " = " << "..." );
    ATH_CHECK(outputeFexEMxTOBHandle.record(std::move(xeContainer),std::move(xeAuxContainer)));
 
    SG::WriteHandle<xAOD::eFexEMRoIContainer> outputeFexHandle(m_eFexOutKey/*, ctx*/);
    ATH_MSG_DEBUG("  write: " << outputeFexHandle.key() << " = " << "..." );
    ATH_CHECK(outputeFexHandle.record(std::move(eContainer),std::move(eAuxContainer)));
 
    SG::ReadHandle<TrigConf::L1Menu> l1Menu (m_l1MenuKey/*, ctx*/);
    ATH_CHECK(l1Menu.isValid());
 
    auto & thr_eTAU = l1Menu->thrExtraInfo().eTAU();
    int activeAlgo = (thr_eTAU.algoVersion() == 0 || thr_eTAU.algoVersion() == 3) ? xAOD::eFexTauRoI_v1::Heuristic : xAOD::eFexTauRoI_v1::BDT;
    bool omitAltTauContainer = m_eFexTauAltxTOBOutKey.empty() || m_eFexTauAltOutKey.empty();
 
    // Repeat for Tau TOBs and xTOBs
    if (activeAlgo == xAOD::eFexTauRoI_v1::Heuristic) {
        ATH_CHECK(StoreTauTOBs(allTauHeuristicTobObjects, m_eFexTauActivexTOBOutKey, m_eFexTauActiveOutKey));
        if (!omitAltTauContainer) {
            ATH_CHECK(StoreTauTOBs(allTauBDTTobObjects, m_eFexTauAltxTOBOutKey, m_eFexTauAltOutKey));
        }
    } else if (activeAlgo == xAOD::eFexTauRoI_v1::BDT) {
        ATH_CHECK(StoreTauTOBs(allTauBDTTobObjects, m_eFexTauActivexTOBOutKey, m_eFexTauActiveOutKey));
        if (!omitAltTauContainer) {
            ATH_CHECK(StoreTauTOBs(allTauHeuristicTobObjects, m_eFexTauAltxTOBOutKey, m_eFexTauAltOutKey));
        }
    }
 
    //Send TOBs to bytestream?
    // ToDo
    // To implement
    // {--Implement--}
 
    return StatusCode::SUCCESS;
 
  }

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase & ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

initialize()

StatusCode LVL1::eFEXSysSim::initialize ( )

virtual

standard Athena-Algorithm method

Destructor.

Definition at line 52 of file eFEXSysSim.cxx.

  {
    
    ATH_CHECK( m_eTowerContainerSGKey.initialize() );
 
    ATH_CHECK( m_eFEXSimTool.retrieve() );
 
    ATH_CHECK( m_eFexOutKey.initialize() );
    ATH_CHECK( m_eFexEMxTOBOutKey.initialize() );
 
    ATH_CHECK( m_eFexTauActiveOutKey.initialize() );
    ATH_CHECK( m_eFexTauActivexTOBOutKey.initialize() );
    ATH_CHECK( m_eFexTauAltOutKey.initialize(SG::AllowEmpty) );
    ATH_CHECK( m_eFexTauAltxTOBOutKey.initialize(SG::AllowEmpty) );
 
    ATH_CHECK( m_eFEXFPGATowerIdProviderTool.retrieve() );
 
    ATH_CHECK(m_l1MenuKey.initialize());
 
    return StatusCode::SUCCESS;
  }

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

interfaceID()

const InterfaceID & LVL1::eFEXSysSim::interfaceID ( )

inlinestatic

Definition at line 45 of file eFEXSysSim.h.

{ return IID_IeFEXSysSim; };

msg()

MsgStream & AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

operator=()

eFEXSysSim && LVL1::eFEXSysSim::operator= ( const eFEXSysSim & )

delete

Destructor.

outputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

renounce()

std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void > AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray & handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

StoreTauTOBs()

StatusCode LVL1::eFEXSysSim::StoreTauTOBs ( std::map< int, std::vector< std::unique_ptr< eFEXtauTOB > > > & allTauTobObjects, SG::WriteHandleKey< xAOD::eFexTauRoIContainer > & eFexTauxTOBOutKey, SG::WriteHandleKey< xAOD::eFexTauRoIContainer > & eFexTauOutKey )

private

Definition at line 455 of file eFEXSysSim.cxx.

 {
    std::unique_ptr< xAOD::eFexTauRoIContainer > tauContainer;
    std::unique_ptr< xAOD::eFexTauRoIAuxContainer > tauAuxContainer;
    std::unique_ptr< xAOD::eFexTauRoIContainer > xtauContainer;
    std::unique_ptr< xAOD::eFexTauRoIAuxContainer > xtauAuxContainer;
 
    tauContainer = std::make_unique<xAOD::eFexTauRoIContainer> ();
    tauAuxContainer = std::make_unique<xAOD::eFexTauRoIAuxContainer> ();
    tauContainer->setStore(tauAuxContainer.get());
 
    xtauContainer = std::make_unique<xAOD::eFexTauRoIContainer> ();
    xtauAuxContainer = std::make_unique<xAOD::eFexTauRoIAuxContainer> ();
    xtauContainer->setStore(xtauAuxContainer.get());
 
    // iterate over all tau TOBs and fill xTOB EDM with them
    for( auto const& [efex, tobObjects] : allTauTobObjects ){
      for( auto &tobObject: tobObjects ){
        m_eFEXFillEDMTool->fillTauEDM(xtauContainer, efex, tobObject, true);
      }
    }
 
    // Form list of TOBs, sorted and truncated from each eFEX
    // Vector to store sorted TOBs from all eFEXes
    // Note that this step means moving TOBs from the all objects vector, which is why we do it last
    std::vector<std::unique_ptr<eFEXtauTOB>> tauTOBs;
    // Loop through eFEXes and sort TOBs from each
    auto iterTau = allTauTobObjects.begin();
    while (iterTau != allTauTobObjects.end()) {
        std::vector<std::unique_ptr<eFEXtauTOB>> tobsSort = std::move(iterTau->second);
        // sort tobs by their et (last 12 bits of the 32 bit tob word)
        std::sort (tobsSort.begin(), tobsSort.end(), std::bind(TOBetSort<std::unique_ptr<eFEXtauTOB>>, std::placeholders::_1, std::placeholders::_2, true));
        // Truncate at 6 TOBs per eFEX
        if (tobsSort.size() > 6) tobsSort.resize(6);
        // Append to system TOB list
        for (unsigned int t = 0; t < tobsSort.size(); ++t) tauTOBs.push_back(std::move(tobsSort[t]));
        // Next eFEX
        ++iterTau;
    }
 
    // iterate over sorted eFEX Tau TOBs and fill TOB EDM with them
    for(auto &tobObject : tauTOBs){
        int efex = tobObject->geteFEXID();
        m_eFEXFillEDMTool->fillTauEDM(tauContainer, efex, tobObject);
    }
 
    // Match xTOBs and TOBs and set isTOB flags if matched
    matchTOBs(tauContainer, xtauContainer);
 
    // Record containers in StoreGate
    SG::WriteHandle<xAOD::eFexTauRoIContainer> outputeFexTauxTOBHandle(eFexTauxTOBOutKey/*, ctx*/);
    ATH_MSG_DEBUG(" write: " << outputeFexTauxTOBHandle.key() << " = " << "..." );
    ATH_CHECK(outputeFexTauxTOBHandle.record(std::move(xtauContainer), std::move(xtauAuxContainer)));
 
    SG::WriteHandle<xAOD::eFexTauRoIContainer> outputeFexTauHandle(eFexTauOutKey/*, ctx*/);
    ATH_MSG_DEBUG(" write: " << outputeFexTauHandle.key() << " = " << "..." );
    ATH_CHECK(outputeFexTauHandle.record(std::move(tauContainer), std::move(tauAuxContainer)));
 
    return StatusCode::SUCCESS;
 
  }

sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in asg::AsgMetadataTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and DerivationFramework::CfAthAlgTool.

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

TOBetSort()

template<class TOBObjectClass>

bool LVL1::eFEXSysSim::TOBetSort ( const TOBObjectClass & i, const TOBObjectClass & j, bool isTau )

inlinestaticprivate

Internal data.

Definition at line 62 of file eFEXSysSim.h.

                                                                                                                         {
      auto et_i = (i->getTobword() >> 0 ) & 0xfff;
      auto et_j = (j->getTobword() >> 0 ) & 0xfff;
      if(et_i > et_j) return true;
      if(et_i==et_j) {
        // resolve ties with procNumber (2,1,3,0 for em, 2,1,0,3 for tau), then phi, then eta
        auto procOrder = (isTau) ? std::map<unsigned int,unsigned int>{{2,3},{1,2},{0,1},{3,0}} :
                                   std::map<unsigned int,unsigned int>{{2,3},{1,2},{3,1},{0,0}};
        auto proc_i = procOrder.at((i->getTobword()) >> 30);
        auto proc_j = procOrder.at((j->getTobword()) >> 30);
        if(proc_i > proc_j) return true;
        if (proc_i == proc_j) {
          auto phi_i = (i->getTobword()) & 0x07000000;
          auto phi_j = (j->getTobword()) & 0x07000000;
          if(phi_i > phi_j) return true;
          if(phi_i == phi_j) {
            auto eta_i = (i->getTobword()) & 0x38000000;
            auto eta_j = (j->getTobword()) & 0x38000000;
            if(eta_i > eta_j) return true;
          }
 
        }
      }
      return false;
    }

updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase & )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_eFexEMxTOBOutKey

SG::WriteHandleKey< xAOD::eFexEMRoIContainer > LVL1::eFEXSysSim::m_eFexEMxTOBOutKey {this,"Key_eFexEMxTOBOutputContainer","L1_eEMxRoI","Output eFexEM xTOB container"}

private

Definition at line 102 of file eFEXSysSim.h.

{this,"Key_eFexEMxTOBOutputContainer","L1_eEMxRoI","Output eFexEM xTOB container"};

m_eFEXFillEDMTool

ToolHandle<eFEXFillEDM> LVL1::eFEXSysSim::m_eFEXFillEDMTool {this, "eFEXFillEDMTool", "LVL1::eFEXFillEDM", "Tool to fille eFEX EDMs"}

private

Definition at line 96 of file eFEXSysSim.h.

{this, "eFEXFillEDMTool", "LVL1::eFEXFillEDM", "Tool to fille eFEX EDMs"};

m_eFEXFPGATowerIdProviderTool

ToolHandle<eFEXFPGATowerIdProvider> LVL1::eFEXSysSim::m_eFEXFPGATowerIdProviderTool {this, "eFEXFPGATowerIdProviderTool", "LVL1::eFEXFPGATowerIdProvider", "Tool that provides tower-FPGA mapping"}

private

Definition at line 107 of file eFEXSysSim.h.

{this, "eFEXFPGATowerIdProviderTool", "LVL1::eFEXFPGATowerIdProvider", "Tool that provides tower-FPGA mapping"};

m_eFexOutKey

SG::WriteHandleKey< xAOD::eFexEMRoIContainer > LVL1::eFEXSysSim::m_eFexOutKey {this,"Key_eFexEMOutputContainer","L1_eEMRoI","Output eFexEM TOB container"}

private

Definition at line 101 of file eFEXSysSim.h.

{this,"Key_eFexEMOutputContainer","L1_eEMRoI","Output eFexEM TOB container"};

m_eFEXSimTool

ToolHandle<eFEXSim> LVL1::eFEXSysSim::m_eFEXSimTool {this, "eFEXSimTool", "LVL1::eFEXSim", "Tool that creates the eFEX Simulation"}

private

Definition at line 94 of file eFEXSysSim.h.

{this, "eFEXSimTool",    "LVL1::eFEXSim",    "Tool that creates the eFEX Simulation"};

m_eFexTauActiveOutKey

SG::WriteHandleKey< xAOD::eFexTauRoIContainer > LVL1::eFEXSysSim::m_eFexTauActiveOutKey {this,"Key_eFexTauOutputContainer","L1_eTauRoI","Output eFexTau active (BDT/heuristic) algorithm TOB container"}

private

Definition at line 103 of file eFEXSysSim.h.

{this,"Key_eFexTauOutputContainer","L1_eTauRoI","Output eFexTau active (BDT/heuristic) algorithm TOB container"};

m_eFexTauActivexTOBOutKey

SG::WriteHandleKey< xAOD::eFexTauRoIContainer > LVL1::eFEXSysSim::m_eFexTauActivexTOBOutKey {this,"Key_eFexTauxTOBOutputContainer","L1_eTauxRoI","Output eFexTau active (BDT/heuristic) algorithm xTOB container"}

private

Definition at line 104 of file eFEXSysSim.h.

{this,"Key_eFexTauxTOBOutputContainer","L1_eTauxRoI","Output eFexTau active (BDT/heuristic) algorithm xTOB container"};

m_eFexTauAltOutKey

SG::WriteHandleKey< xAOD::eFexTauRoIContainer > LVL1::eFEXSysSim::m_eFexTauAltOutKey {this,"Key_eFexAltTauOutputContainer","","Output eFexTau alternative (BDT/heuristic) algorithm TOB container"}

private

Definition at line 105 of file eFEXSysSim.h.

{this,"Key_eFexAltTauOutputContainer","","Output eFexTau alternative (BDT/heuristic) algorithm TOB container"};

m_eFexTauAltxTOBOutKey

SG::WriteHandleKey< xAOD::eFexTauRoIContainer > LVL1::eFEXSysSim::m_eFexTauAltxTOBOutKey {this,"Key_eFexAltTauxTOBOutputContainer","","Output eFexTau alternative (BDT/heuristic) algorithm xTOB container"}

private

Definition at line 106 of file eFEXSysSim.h.

{this,"Key_eFexAltTauxTOBOutputContainer","","Output eFexTau alternative (BDT/heuristic) algorithm xTOB container"};

m_eTowerContainerSGKey

SG::ReadHandleKey<LVL1::eTowerContainer> LVL1::eFEXSysSim::m_eTowerContainerSGKey {this, "MyETowers", "eTowerContainer", "Input container for eTowers"}

private

Definition at line 98 of file eFEXSysSim.h.

{this, "MyETowers", "eTowerContainer", "Input container for eTowers"};

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_l1MenuKey

SG::ReadHandleKey<TrigConf::L1Menu> LVL1::eFEXSysSim::m_l1MenuKey {this, "L1TriggerMenu", "DetectorStore+L1TriggerMenu","Name of the L1Menu object to read configuration from"}

private

Definition at line 99 of file eFEXSysSim.h.

{this, "L1TriggerMenu", "DetectorStore+L1TriggerMenu","Name of the L1Menu object to read configuration from"}; 

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

The documentation for this class was generated from the following files:

• eFEXSysSim.h
• eFEXSysSim.cxx