#include <jFEXForwardElecAlgo.h>

Inheritance diagram for LVL1::jFEXForwardElecAlgo:

Collaboration diagram for LVL1::jFEXForwardElecAlgo:

Public Member Functions
	jFEXForwardElecAlgo (const std::string &type, const std::string &name, const IInterface *parent)
	Constructors.
virtual StatusCode	initialize () override
	standard Athena-Algorithm method
virtual	~jFEXForwardElecAlgo ()
	Destructor.
virtual StatusCode	safetyTest () override
	Standard methods.
virtual StatusCode	reset () override
virtual void	setup (int inputTable[FEXAlgoSpaceDefs::jFEX_algoSpace_height][FEXAlgoSpaceDefs::jFEX_wide_algoSpace_width], int jfex, int fpga) override
virtual std::unordered_map< uint, jFEXForwardElecInfo >	calculateEDM () override
virtual void	setFPGAEnergy (std::unordered_map< int, std::vector< int > > etmapEM, std::unordered_map< int, std::vector< int > > etmapHAD) override
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
virtual std::array< float, 2 >	getEtaPhi (uint) override
virtual std::array< int, 2 >	getEtEmHad (uint) const override
bool	getEMSat (unsigned int ttID)
bool	isValidSeed (uint seedTTID) const
void	findAndFillNextTT (jFEXForwardElecInfo &elCluster, int neta, int nphi)
StatusCode	ReadfromFile (const std::string &, std::unordered_map< unsigned int, std::vector< unsigned int > > &) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
SG::ReadHandleKey< LVL1::jTowerContainer >	m_jTowerContainerKey {this, "MyjTowers", "jTowerContainer", "jTower input container"}
SG::ReadHandle< jTowerContainer >	m_jTowerContainer
std::unordered_map< int, std::vector< int > >	m_map_Etvalues_EM
std::unordered_map< int, std::vector< int > >	m_map_Etvalues_HAD
int	m_jFEXalgoTowerID [FEXAlgoSpaceDefs::jFEX_algoSpace_height][FEXAlgoSpaceDefs::jFEX_wide_algoSpace_width]
int	m_jfex
int	m_fpga
const int	m_Edge_dR2 = std::round( (std::pow(2M_PI/32,2)) 1e5 )
const int	m_Edge_dR3 = std::round( (std::pow(3M_PI/32,2)) 1e5 )
const int	m_Edge_dR4 = std::round( (std::pow(4M_PI/32,2)) 1e5 )
Gaudi::Property< std::string >	m_IsoMapStr {this, "IsoMap", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_iso.dat", "Contains Trigger towers in (forward) EM layer used for isolation"}
Gaudi::Property< std::string >	m_Frac1MapStr {this, "Frac1Map", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_frac.dat", "Contains Trigger towers in FCal layer2 used for hadronic fraction 1 discriminant"}
Gaudi::Property< std::string >	m_Frac2MapStr {this, "Frac2Map", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_frac2.dat", "Contains Trigger towers in FCal layer3 used for hadronic fraction 2 discriminant"}
Gaudi::Property< std::string >	m_SearchGTauStr {this, "SearchGTauMap", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_searchGTau.dat" , "Contains Trigger tower to find local max (greater than)"}
Gaudi::Property< std::string >	m_SearchGeTauStr {this, "SearchGeTauMap", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_searchGeTau.dat", "Contains Trigger tower to find local max (greater or equal than)"}
std::unordered_map< unsigned int, std::vector< unsigned int > >	m_SeedRingMap
std::unordered_map< unsigned int, std::vector< unsigned int > >	m_1stRingMap
std::unordered_map< unsigned int, std::vector< unsigned int > >	m_IsoMap
std::unordered_map< unsigned int, std::vector< unsigned int > >	m_Frac1Map
std::unordered_map< unsigned int, std::vector< unsigned int > >	m_Frac2Map
std::unordered_map< unsigned int, std::vector< unsigned int > >	m_SearchGTauMap
std::unordered_map< unsigned int, std::vector< unsigned int > >	m_SearchGeTauMap
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

Static Private Attributes
static constexpr float	m_2PI = 2*M_PI
static constexpr float	m_TT_Size_phi = M_PI/32

Detailed Description

Definition at line 25 of file jFEXForwardElecAlgo.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

◆ jFEXForwardElecAlgo()

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

Constructors.

Definition at line 28 of file jFEXForwardElecAlgo.cxx.

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

◆ ~jFEXForwardElecAlgo()

LVL1::jFEXForwardElecAlgo::~jFEXForwardElecAlgo ( )

virtual

Destructor.

Definition at line 33 of file jFEXForwardElecAlgo.cxx.

33 {

34 }

Member Function Documentation

◆ calculateEDM()

std::unordered_map< uint, LVL1::jFEXForwardElecInfo > LVL1::jFEXForwardElecAlgo::calculateEDM ( void )

overridevirtual

Implements LVL1::IjFEXForwardElecAlgo.

Definition at line 231 of file jFEXForwardElecAlgo.cxx.

                                                                                            {
    std::unordered_map<uint, LVL1::jFEXForwardElecInfo> clusterList;
    std::vector<int> lower_centre_neta;
    std::vector<int> upper_centre_neta;
   
    //check if we are in module 0 or 5 and assign corrrect eta FEXAlgoSpace parameters
    if(m_jfex == 0) {
      //Module 0 
      lower_centre_neta.assign({FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMB_start_eta, // 28
                                FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMIE_start_eta, // 24
                                FEXAlgoSpaceDefs::jFEX_algoSpace_C_FCAL_start_eta}); // 12
      
      upper_centre_neta.assign({FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMB_end_eta, // 37
                                FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMIE_end_eta, // 28
                                FEXAlgoSpaceDefs::jFEX_algoSpace_C_FCAL_end_eta }); // 24
    }
    else {
      //Module 5
      lower_centre_neta.assign({FEXAlgoSpaceDefs::jFEX_algoSpace_A_EMB_eta, // 8
                                FEXAlgoSpaceDefs::jFEX_algoSpace_A_EMIE_eta, // 17
                                FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL_start_eta}); // 21
                                
      upper_centre_neta.assign({FEXAlgoSpaceDefs::jFEX_algoSpace_A_EMIE_eta, // 17
                                FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL_start_eta, // 21 
                                FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL_end_eta}); // 33
    }
 
    //define phi FEXAlgoSpace parameters
    std::vector<int> lower_centre_nphi{FEXAlgoSpaceDefs::jFEX_algoSpace_EMB_start_phi, FEXAlgoSpaceDefs::jFEX_algoSpace_EMIE_start_phi,  FEXAlgoSpaceDefs::jFEX_algoSpace_FCAL_start_phi};
    std::vector<int> upper_centre_nphi{FEXAlgoSpaceDefs::jFEX_algoSpace_EMB_end_phi, FEXAlgoSpaceDefs::jFEX_algoSpace_EMIE_end_phi,  FEXAlgoSpaceDefs::jFEX_algoSpace_FCAL_end_phi};
 
    //loop over different EM/FCAL1 eta phi core fpga regions with different granularities. These are potential seed  towers for electron clusters 
    for(uint region = 0; region<3; region++) {
      for(int nphi = lower_centre_nphi[region]; nphi < upper_centre_nphi[region]; nphi++) {
          for(int neta = lower_centre_neta[region]; neta < upper_centre_neta[region]; neta++) {
          
          // ignore  seeds for |eta| < 2.3 or from the first FCAL eta bin                
          if (m_jfex == 0 && neta >= FEXAlgoSpaceDefs::jFEX_algoSpace_C_FwdEl_start) continue;
          if (m_jfex == 5 && neta <= FEXAlgoSpaceDefs::jFEX_algoSpace_A_FwdEl_start) continue;
          if (m_jfex == 0 && neta == FEXAlgoSpaceDefs::jFEX_algoSpace_C_FCAL1_1st)   continue;
          if (m_jfex == 5 && neta == FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL1_1st)   continue;
             
          // define ttID (only FCAL1 in the third region) which will be the key for class in map, ignore tower ID = 0
          uint ttID = m_jFEXalgoTowerID[nphi][neta];
          if(ttID == 0) continue;
          //check if seed candidate is actual seed (passes local maximum criteria)
          if(!isValidSeed(ttID)) continue;
          
          //gather some first, basic information for resulting cluster/TOB
          jFEXForwardElecInfo elCluster;
          elCluster.setup(m_jfex, ttID, neta, nphi);
          const auto [centreTT_eta, centreTT_phi] = getEtaPhi(ttID);
          const auto [centreTT_EtEM, centreTT_EtHad] = getEtEmHad(ttID);
          elCluster.setCoreTTfPhi(centreTT_phi);
          elCluster.setCoreTTfEta(centreTT_eta);
          elCluster.setCoreTTEtEM(centreTT_EtEM);  
          elCluster.setCoreTTSatEM(getEMSat(ttID));                
          elCluster.setNextTTEtEM(0);
          elCluster.setNextTTID(0);
          elCluster.setTTEtEMiso(0);
 
          //find "NextTT", i.e., highest ET neighbour
          findAndFillNextTT(elCluster, neta, nphi);
 
          // sum up EM isolation using the isolation map and remove cluster ET
          {
            int sumEtEM = 0;
            auto it_iso_map = m_IsoMap.find(ttID);
            if(it_iso_map != m_IsoMap.end()) {
              for(const auto& gtt : it_iso_map->second){
                  auto [tmp_EtEM,tmp_EtHad] = getEtEmHad(gtt);
                  sumEtEM += tmp_EtEM;  
              }
              elCluster.setTTEtEMiso(sumEtEM-elCluster.getNextTTEtEM());
            } else {
                ATH_MSG_ERROR("Could not find TT" << ttID << " in jEM isolation map file.");
            }
          }
          
          if(fabs(centreTT_eta) < 3.2) {
            // for non-FCal positions only Frac1 is meaningful and has a "trivial" mapping
            elCluster.setTTEtHad1(centreTT_EtHad);
            elCluster.setTTEtHad2(0);
          } else {
            // sum up Et for hadronic fraction 1
            {
              int sumEtHad1 = 0;
              auto it_frac1_map = m_Frac1Map.find(ttID);
              if(it_frac1_map != m_Frac1Map.end()) {
                for(const auto& gtt : it_frac1_map->second){
                    auto [tmp_EtEM,tmp_EtHad] = getEtEmHad(gtt);
                    sumEtHad1 += tmp_EtHad;  
                }
                elCluster.setTTEtHad1(sumEtHad1);
              } else { 
                ATH_MSG_ERROR("Could not find TT" << ttID << " in jEM frac1 map file.");
              }
            }
            
            // sum up Et for hadronic fraction 2 (only FCal!)
            {
              int sumEtHad2 = 0;
              auto it_frac2_map = m_Frac2Map.find(ttID);
              if(it_frac2_map != m_Frac2Map.end()) {
                for(const auto& gtt : it_frac2_map->second) {
                    auto [tmp_EtEM,tmp_EtHad] = getEtEmHad(gtt);
                    sumEtHad2 += tmp_EtHad;  
                }
                elCluster.setTTEtHad2(sumEtHad2);
              } else { 
                ATH_MSG_ERROR("Could not find TT" << ttID << " in jEM frac2 map file.");
              }
            }
          }
 
          // save this cluster in the list
          clusterList[ttID] = std::move(elCluster);
        }//eta
      }//phi
    }// 3 regions
    return clusterList;
  } 

◆ 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.

95{ return m_detStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

StoreGateSvc_t m_detStore

Definition AthCommonDataStore.h:393

◆ 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.

85{ return m_evtStore; }

AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

StoreGateSvc_t m_evtStore

Definition AthCommonDataStore.h:390

◆ 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

◆ findAndFillNextTT()

void LVL1::jFEXForwardElecAlgo::findAndFillNextTT	(	jFEXForwardElecInfo &	elCluster,
		int	neta,
		int	nphi )

private

Definition at line 146 of file jFEXForwardElecAlgo.cxx.

                                                                                                    {
    //determine direction for offsets (-1 : C-side, +1: A-side)
    int direction = m_jfex < 3 ? -1 : +1;
    
    std::vector<std::pair<int,int>> neighbours;
    //eta/phi index offsets depend on the position of the seed the cluster is created from
    //differentiate the (many) different cases:
    if ( direction>0 ? neta < FEXAlgoSpaceDefs::jFEX_algoSpace_A_EMIE_eta - 1    // A-side condition
                     : neta > FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMB_start_eta)  // C-side condition
    { //Barrel
      neighbours = { {neta  , nphi+1},
                     {neta  , nphi-1},
                     {neta+1, nphi  },
                     {neta-1, nphi  }
                   };
    } else if ( direction>0 ? neta == FEXAlgoSpaceDefs::jFEX_algoSpace_A_EMIE_eta - 1
                            : neta == FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMB_start_eta ) 
    { //Barrel next to Endcap
      neighbours = { {neta, nphi+1},
                     {neta, nphi-1},
                     {neta-direction, nphi},
                     {neta+direction, nphi/2} //boundary crossing into endcap -> reduced phi granularity
                   };
    } else if ( direction>0 ? neta == FEXAlgoSpaceDefs::jFEX_algoSpace_A_EMIE_eta 
                            : neta == FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMB_start_eta - 1 ) 
    { //Endcap next to Barrel
      neighbours = { {neta, nphi+1},
                     {neta, nphi-1},
                     {neta+direction, nphi},
                     {neta-direction, 2*nphi+0}, //crossing into barrel region, higher phi granularity
                     {neta-direction, 2*nphi+1}, // -> consider both "touching" towers 
                   };
    } else if ( direction>0 ? neta > FEXAlgoSpaceDefs::jFEX_algoSpace_A_EMIE_eta         && neta < FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL_start_eta - 1
                            : neta < FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMIE_end_eta - 1 && neta > FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMIE_start_eta ) 
    { //Endcap
      neighbours = { {neta  , nphi+1},
                     {neta  , nphi-1},
                     {neta+1, nphi  },
                     {neta-1, nphi  }
                   };
    } else if ( direction>0 ? neta == FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL_start_eta - 1
                            : neta == FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMIE_start_eta ) 
    { //Endcap next to FCal
      neighbours = { {neta, nphi+1},
                     {neta, nphi-1},
                     {neta-direction, nphi},
                     {neta+2*direction, nphi/2} //boundary crossing into FCal -> reduced phi granularity and skip first FCal bin
                   };
    } else if ( direction>0 ? neta == FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL_start_eta + 1   //first FCal bin must be skipped!
                            : neta == FEXAlgoSpaceDefs::jFEX_algoSpace_C_EMIE_start_eta - 2 ) //first FCal bin must be skipped!
    { //FCal next to Endcap
      //phi spacing in FCal is very wide, no longer consider adding towers in phi direction
      neighbours = { {neta-2*direction, 2*nphi+0}, //boundary crossing into endcap, higher phi granularity
                     {neta-2*direction, 2*nphi+1}, // -> consider both "touching" towers
                     {neta+direction, nphi}
                   };
    } else if ( direction>0 ? neta > FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL_start_eta + 1 && neta < FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL_end_eta - 1
                            : neta < FEXAlgoSpaceDefs::jFEX_algoSpace_C_FCAL_end_eta - 2   && neta > FEXAlgoSpaceDefs::jFEX_algoSpace_C_FCAL_start_eta ) 
    { //FCal
      neighbours = { {neta+1, nphi},
                     {neta-1, nphi}
                   };
    } else if ( direction>0 ? neta == FEXAlgoSpaceDefs::jFEX_algoSpace_A_FCAL_end_eta - 1 
                            : neta == FEXAlgoSpaceDefs::jFEX_algoSpace_C_FCAL_start_eta  ) 
    { //FCal, last eta bin
      neighbours = { {neta-direction, nphi}
                   };    
    } else {
      ATH_MSG_ERROR("Eta index " << neta << " (side: "<< (direction>0?"A":"C") << ") does not seem to belong to any valid seed region");
    }
    
    //iterate over neighbours, find most energetic one
    for (const auto& candEtaPhi: neighbours) {
      uint candID = m_jFEXalgoTowerID[candEtaPhi.second][candEtaPhi.first];
      const auto [candTT_EtEM, candTT_EtHad] = getEtEmHad(candID);
      if (candTT_EtEM > elCluster.getNextTTEtEM()) {
        elCluster.setNextTTEtEM(candTT_EtEM);
        elCluster.setNextTTID(candID);
        elCluster.setNextTTSatEM(getEMSat(candID));
      }
    }
    
  }

◆ getEMSat()

bool LVL1::jFEXForwardElecAlgo::getEMSat ( unsigned int ttID )

private

Definition at line 93 of file jFEXForwardElecAlgo.cxx.

                                                           {
    if(ttID == 0) {
        return false;
    } 
    
    const LVL1::jTower * tmpTower = m_jTowerContainer->findTower(ttID);
    return tmpTower->getEMSat();
  }

◆ getEtaPhi()

std::array< float, 2 > LVL1::jFEXForwardElecAlgo::getEtaPhi ( uint ttID )

overrideprivatevirtual

Implements LVL1::IjFEXForwardElecAlgo.

Definition at line 69 of file jFEXForwardElecAlgo.cxx.

                                                                {
    if(ttID == 0) {
      return {999,999};
    }
    const LVL1::jTower *tmpTower = m_jTowerContainer->findTower(ttID);
    return {tmpTower->centreEta(),tmpTower->centrePhi()};
  }

◆ getEtEmHad()

std::array< int, 2 > LVL1::jFEXForwardElecAlgo::getEtEmHad ( uint ttID ) const

overrideprivatevirtual

Implements LVL1::IjFEXForwardElecAlgo.

Definition at line 77 of file jFEXForwardElecAlgo.cxx.

                                                                     {
    if(ttID == 0) {
      return {0,0};
    }
    int TT_EtEM = 0;
    if(m_map_Etvalues_EM.find(ttID) != m_map_Etvalues_EM.end()) {
      TT_EtEM = m_map_Etvalues_EM.at(ttID)[0];
    }
    int TT_EtHad = 0;
    if(m_map_Etvalues_HAD.find(ttID) != m_map_Etvalues_HAD.end()) {
      TT_EtHad = m_map_Etvalues_HAD.at(ttID)[0];
    }
 
    return {TT_EtEM, TT_EtHad};
  }

◆ initialize()

StatusCode LVL1::jFEXForwardElecAlgo::initialize ( )

overridevirtual

standard Athena-Algorithm method

Definition at line 36 of file jFEXForwardElecAlgo.cxx.

                                                 {
    ATH_CHECK(m_jTowerContainerKey.initialize());
 
    ATH_CHECK(ReadfromFile(PathResolver::find_calib_file(m_IsoMapStr), m_IsoMap ));
    ATH_CHECK(ReadfromFile(PathResolver::find_calib_file(m_Frac1MapStr), m_Frac1Map ));
    ATH_CHECK(ReadfromFile(PathResolver::find_calib_file(m_Frac2MapStr), m_Frac2Map ));
    ATH_CHECK(ReadfromFile(PathResolver::find_calib_file(m_SearchGTauStr), m_SearchGTauMap));
    ATH_CHECK(ReadfromFile(PathResolver::find_calib_file(m_SearchGeTauStr), m_SearchGeTauMap));
 
    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::IjFEXForwardElecAlgo::interfaceID ( )

inlinestaticinherited

Definition at line 40 of file IjFEXForwardElecAlgo.h.

{
    return IID_IjFEXForwardElecAlgo;
}

◆ isValidSeed()

bool LVL1::jFEXForwardElecAlgo::isValidSeed ( uint seedTTID ) const

private

Definition at line 110 of file jFEXForwardElecAlgo.cxx.

                                                               {
    auto [centreTT_EtEM,centreTT_EtHad] = getEtEmHad(seedTTID);
    // check if seed has strictly more energy than its neighbours
    {
      auto it_seed_map = m_SearchGTauMap.find(seedTTID);
      if(it_seed_map == m_SearchGTauMap.end()) {
          ATH_MSG_ERROR("Could not find TT" << seedTTID << " in the seach (>) local maxima for tau/em file.");
          return false;
      }
      for (const auto& gtt : it_seed_map->second ){
        auto [tmp_EtEM,tmp_EtHad] = getEtEmHad(gtt);
        if(tmp_EtEM>=centreTT_EtEM) {
          return false;
        }
      }
    }
 
    // check if seed has equal or more energy than its neighbours
    {
      auto it_seed_map = m_SearchGeTauMap.find(seedTTID);
      if(it_seed_map == m_SearchGeTauMap.end()) {
          ATH_MSG_ERROR("Could not find TT" << seedTTID << " in the seach (>=) local maxima for tau/em file.");
          return false;
      }
      for (const auto& gtt : it_seed_map->second ){
        auto [tmp_EtEM,tmp_EtHad] = getEtEmHad(gtt);
        if( tmp_EtEM>centreTT_EtEM) {
          return false;
        }
      }
    }
 
    return true;  
  }

◆ 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);
  }

◆ 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.

◆ ReadfromFile()

StatusCode LVL1::jFEXForwardElecAlgo::ReadfromFile	(	const std::string &	fileName,
		std::unordered_map< unsigned int, std::vector< unsigned int > > &	fillingMap ) const

private

Definition at line 354 of file jFEXForwardElecAlgo.cxx.

                                                                                                                                                         {
    std::string myline;
    //opening file with ifstream
    std::ifstream myfile(fileName);
    if ( !myfile.is_open() ){
        ATH_MSG_ERROR("Could not open file:" << fileName);
        return StatusCode::FAILURE;
    }
    
    //loading the mapping information
    while ( std::getline (myfile, myline) ) {
        //removing the header of the file (it is just information!)
        if(myline[0] == '#') continue;
        
        //Splitting myline in different substrings
        std::stringstream oneLine(myline);
        
        //reading elements
        std::vector<unsigned int> elements;
        std::string element;
        while(std::getline(oneLine, element, ' '))
        {
            elements.push_back(std::stoi(element));
        }
        
        // We should have at least two elements! Central TT and (at least) itself
        if(elements.size() < 1){
            ATH_MSG_ERROR("Unexpected number of elemennts (<1 expected) in file: "<< fileName);
            return StatusCode::FAILURE;
        }
        
        //Central TiggerTower
        unsigned int TTID = elements.at(0);
        
        // rest of TTs that need to be checked
        elements.erase(elements.begin());
        fillingMap[TTID] = std::move(elements);        
    }
    myfile.close();
 
    return StatusCode::SUCCESS;
  }   

◆ 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();
  }

◆ reset()

StatusCode LVL1::jFEXForwardElecAlgo::reset ( )

overridevirtual

Implements LVL1::IjFEXForwardElecAlgo.

Definition at line 58 of file jFEXForwardElecAlgo.cxx.

                                            {
    return StatusCode::SUCCESS;
  }

◆ safetyTest()

StatusCode LVL1::jFEXForwardElecAlgo::safetyTest ( )

overridevirtual

Standard methods.

Implements LVL1::IjFEXForwardElecAlgo.

Definition at line 49 of file jFEXForwardElecAlgo.cxx.

                                                 {
    m_jTowerContainer = SG::ReadHandle<jTowerContainer>(m_jTowerContainerKey);
    if(! m_jTowerContainer.isValid()) {
      ATH_MSG_ERROR("Could not retrieve jTowerContainer " << m_jTowerContainerKey.key());
      return StatusCode::FAILURE;
    }
    return StatusCode::SUCCESS;
  }

◆ setFPGAEnergy()

void LVL1::jFEXForwardElecAlgo::setFPGAEnergy	(	std::unordered_map< int, std::vector< int > >	etmapEM,
		std::unordered_map< int, std::vector< int > >	etmapHAD )

overridevirtual

Implements LVL1::IjFEXForwardElecAlgo.

Definition at line 102 of file jFEXForwardElecAlgo.cxx.

                                                    {
    m_map_Etvalues_EM=std::move(etmapEM); 
    m_map_Etvalues_HAD=std::move(etmapHAD);
  }

◆ setup()

void LVL1::jFEXForwardElecAlgo::setup	(	int	inputTable[FEXAlgoSpaceDefs::jFEX_algoSpace_height][FEXAlgoSpaceDefs::jFEX_wide_algoSpace_width],
		int	jfex,
		int	fpga )

overridevirtual

Implements LVL1::IjFEXForwardElecAlgo.

Definition at line 62 of file jFEXForwardElecAlgo.cxx.

                                                                                                                                                          {
    std::copy(&inputTable[0][0], &inputTable[0][0] + (FEXAlgoSpaceDefs::jFEX_algoSpace_height*FEXAlgoSpaceDefs::jFEX_wide_algoSpace_width), &m_jFEXalgoTowerID[0][0]);
    m_jfex=jfex;
    m_fpga=fpga;
  }

◆ 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.

◆ 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_1stRingMap

std::unordered_map<unsigned int, std::vector<unsigned int> > LVL1::jFEXForwardElecAlgo::m_1stRingMap

private

Definition at line 71 of file jFEXForwardElecAlgo.h.

◆ m_2PI

float LVL1::jFEXForwardElecAlgo::m_2PI = 2*M_PI

staticconstexprprivate

Definition at line 58 of file jFEXForwardElecAlgo.h.

◆ 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_Edge_dR2

const int LVL1::jFEXForwardElecAlgo::m_Edge_dR2 = std::round( (std::pow(2*M_PI/32,2)) * 1e5 )

private

Definition at line 60 of file jFEXForwardElecAlgo.h.

◆ m_Edge_dR3

const int LVL1::jFEXForwardElecAlgo::m_Edge_dR3 = std::round( (std::pow(3*M_PI/32,2)) * 1e5 )

private

Definition at line 61 of file jFEXForwardElecAlgo.h.

◆ m_Edge_dR4

const int LVL1::jFEXForwardElecAlgo::m_Edge_dR4 = std::round( (std::pow(4*M_PI/32,2)) * 1e5 )

private

Definition at line 62 of file jFEXForwardElecAlgo.h.

◆ 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_fpga

int LVL1::jFEXForwardElecAlgo::m_fpga

private

Definition at line 57 of file jFEXForwardElecAlgo.h.

◆ m_Frac1Map

std::unordered_map<unsigned int, std::vector<unsigned int> > LVL1::jFEXForwardElecAlgo::m_Frac1Map

private

Definition at line 73 of file jFEXForwardElecAlgo.h.

◆ m_Frac1MapStr

Gaudi::Property<std::string> LVL1::jFEXForwardElecAlgo::m_Frac1MapStr {this, "Frac1Map", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_frac.dat", "Contains Trigger towers in FCal layer2 used for hadronic fraction 1 discriminant"}

private

Definition at line 65 of file jFEXForwardElecAlgo.h.

65{this, "Frac1Map", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_frac.dat", "Contains Trigger towers in FCal layer2 used for hadronic fraction 1 discriminant"};

◆ m_Frac2Map

std::unordered_map<unsigned int, std::vector<unsigned int> > LVL1::jFEXForwardElecAlgo::m_Frac2Map

private

Definition at line 74 of file jFEXForwardElecAlgo.h.

◆ m_Frac2MapStr

Gaudi::Property<std::string> LVL1::jFEXForwardElecAlgo::m_Frac2MapStr {this, "Frac2Map", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_frac2.dat", "Contains Trigger towers in FCal layer3 used for hadronic fraction 2 discriminant"}

private

Definition at line 66 of file jFEXForwardElecAlgo.h.

66{this, "Frac2Map", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_frac2.dat", "Contains Trigger towers in FCal layer3 used for hadronic fraction 2 discriminant"};

◆ m_IsoMap

std::unordered_map<unsigned int, std::vector<unsigned int> > LVL1::jFEXForwardElecAlgo::m_IsoMap

private

Definition at line 72 of file jFEXForwardElecAlgo.h.

◆ m_IsoMapStr

Gaudi::Property<std::string> LVL1::jFEXForwardElecAlgo::m_IsoMapStr {this, "IsoMap", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_iso.dat", "Contains Trigger towers in (forward) EM layer used for isolation"}

private

Definition at line 64 of file jFEXForwardElecAlgo.h.

64{this, "IsoMap", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_iso.dat", "Contains Trigger towers in (forward) EM layer used for isolation"};

◆ m_jfex

int LVL1::jFEXForwardElecAlgo::m_jfex

private

Definition at line 56 of file jFEXForwardElecAlgo.h.

◆ m_jFEXalgoTowerID

int LVL1::jFEXForwardElecAlgo::m_jFEXalgoTowerID[FEXAlgoSpaceDefs::jFEX_algoSpace_height][FEXAlgoSpaceDefs::jFEX_wide_algoSpace_width]

private

Definition at line 55 of file jFEXForwardElecAlgo.h.

◆ m_jTowerContainer

SG::ReadHandle<jTowerContainer> LVL1::jFEXForwardElecAlgo::m_jTowerContainer

private

Definition at line 52 of file jFEXForwardElecAlgo.h.

◆ m_jTowerContainerKey

SG::ReadHandleKey<LVL1::jTowerContainer> LVL1::jFEXForwardElecAlgo::m_jTowerContainerKey {this, "MyjTowers", "jTowerContainer", "jTower input container"}

private

Definition at line 51 of file jFEXForwardElecAlgo.h.

51{this, "MyjTowers", "jTowerContainer", "jTower input container"};

◆ m_map_Etvalues_EM

std::unordered_map<int,std::vector<int> > LVL1::jFEXForwardElecAlgo::m_map_Etvalues_EM

private

Definition at line 53 of file jFEXForwardElecAlgo.h.

◆ m_map_Etvalues_HAD

std::unordered_map<int,std::vector<int> > LVL1::jFEXForwardElecAlgo::m_map_Etvalues_HAD

private

Definition at line 54 of file jFEXForwardElecAlgo.h.

◆ m_SearchGeTauMap

std::unordered_map<unsigned int, std::vector<unsigned int> > LVL1::jFEXForwardElecAlgo::m_SearchGeTauMap

private

Definition at line 76 of file jFEXForwardElecAlgo.h.

◆ m_SearchGeTauStr

Gaudi::Property<std::string> LVL1::jFEXForwardElecAlgo::m_SearchGeTauStr {this, "SearchGeTauMap", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_searchGeTau.dat", "Contains Trigger tower to find local max (greater or equal than)"}

private

Definition at line 68 of file jFEXForwardElecAlgo.h.

68{this, "SearchGeTauMap", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_searchGeTau.dat", "Contains Trigger tower to find local max (greater or equal than)"};

◆ m_SearchGTauMap

std::unordered_map<unsigned int, std::vector<unsigned int> > LVL1::jFEXForwardElecAlgo::m_SearchGTauMap

private

Definition at line 75 of file jFEXForwardElecAlgo.h.

◆ m_SearchGTauStr

Gaudi::Property<std::string> LVL1::jFEXForwardElecAlgo::m_SearchGTauStr {this, "SearchGTauMap", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_searchGTau.dat" , "Contains Trigger tower to find local max (greater than)"}

private

Definition at line 67 of file jFEXForwardElecAlgo.h.

67{this, "SearchGTauMap", "Run3L1CaloSimulation/JetMaps/2024_04_09/jFEX_FWD_searchGTau.dat" , "Contains Trigger tower to find local max (greater than)"};

◆ m_SeedRingMap

std::unordered_map<unsigned int, std::vector<unsigned int> > LVL1::jFEXForwardElecAlgo::m_SeedRingMap

private

Definition at line 70 of file jFEXForwardElecAlgo.h.

◆ m_TT_Size_phi

float LVL1::jFEXForwardElecAlgo::m_TT_Size_phi = M_PI/32

staticconstexprprivate

Definition at line 59 of file jFEXForwardElecAlgo.h.

◆ 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:

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ jFEXForwardElecAlgo()

◆ ~jFEXForwardElecAlgo()

Member Function Documentation

◆ calculateEDM()

◆ declareGaudiProperty()

◆ declareProperty()

◆ detStore()

◆ evtStore()

◆ extraDeps_update_handler()

◆ findAndFillNextTT()

◆ getEMSat()

◆ getEtaPhi()

◆ getEtEmHad()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ isValidSeed()

◆ msg()

◆ msgLvl()

◆ outputHandles()

◆ ReadfromFile()

◆ renounce()

◆ renounceArray()

◆ reset()

◆ safetyTest()

◆ setFPGAEnergy()

◆ setup()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_1stRingMap

◆ m_2PI

◆ m_detStore

◆ m_Edge_dR2

◆ m_Edge_dR3

◆ m_Edge_dR4

◆ m_evtStore

◆ m_fpga

◆ m_Frac1Map

◆ m_Frac1MapStr

◆ m_Frac2Map

◆ m_Frac2MapStr

◆ m_IsoMap

◆ m_IsoMapStr

◆ m_jfex

◆ m_jFEXalgoTowerID

◆ m_jTowerContainer

◆ m_jTowerContainerKey

◆ m_map_Etvalues_EM

◆ m_map_Etvalues_HAD

◆ m_SearchGeTauMap

◆ m_SearchGeTauStr

◆ m_SearchGTauMap

◆ m_SearchGTauStr

◆ m_SeedRingMap

◆ m_TT_Size_phi

◆ m_varHandleArraysDeclared

◆ m_vhka