#include <gFEXJwoJAlgo.h>

Inheritance diagram for LVL1::gFEXJwoJAlgo:

Collaboration diagram for LVL1::gFEXJwoJAlgo:

Public Member Functions
	gFEXJwoJAlgo (const std::string &type, const std::string &name, const IInterface *parent)
	Constructors. More...

virtual StatusCode	initialize () override
	standard Athena-Algorithm method More...

virtual void	setAlgoConstant (int aFPGA_A, int bFPGA_A, int aFPGA_B, int bFPGA_B, int aFPGA_C, int bFPGA_C, int gXE_seedThrA, int gXE_seedThrB, int gXE_seedThrC) override

virtual std::vector< std::unique_ptr< gFEXJwoJTOB > >	jwojAlgo (const gTowersType &Atwr, const gTowersType &Btwr, const gTowersType &Ctwr, std::array< uint32_t, 4 > &outTOB) const override

ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard `StoreGateSvc` (event store) Returns (kind of) a pointer to the `StoreGateSvc`. More...

const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard `StoreGateSvc/DetectorStore` Returns (kind of) a pointer to the `StoreGateSvc`. More...

virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...

virtual StatusCode	sysStart () override
	Handle START transition. More...

virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...

virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...

Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...

Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...

void	updateVHKA (Gaudi::Details::PropertyBase &)

MsgStream &	msg () const

MsgStream &	msg (const MSG::Level lvl) 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 More...

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

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
void	gBlockAB (const gTowersType &twrs, gTowersType &gBlkSum, gTowersType &hasSeed, int seedThreshold) const

void	metFPGA (int FPGAnum, const gTowersType &twrs, const gTowersType &gBlkSum, int gBlockthreshold, int aFPGA, int bFPGA, int &MHT_x, int &MHT_y, int &MST_x, int &MST_y, int &MET_x, int &MET_y) const

void	etFPGA (int FPGAnum, const gTowersType &twrs, gTowersType &gBlkSum, int gBlockthreshold, int A, int B, int &eth, int &ets, int &etw) const

void	metTotal (int A_MET_x, int A_MET_y, int B_MET_x, int B_MET_y, int C_MET_x, int C_MET_y, int &MET_x, int &MET_y) const

void	etTotal (int A_ET, int B_ET, int C_ET, int &ET) const

float	sinLUT (unsigned int phiIDX, unsigned int aw) const

float	cosLUT (unsigned int phiIDX, unsigned int aw) const

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...

Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
float	m_aFPGA_A

float	m_bFPGA_A

float	m_aFPGA_B

float	m_bFPGA_B

float	m_aFPGA_C

float	m_bFPGA_C

float	m_gBlockthresholdA

float	m_gBlockthresholdB

float	m_gBlockthresholdC

StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...

StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...

std::vector< SG::VarHandleKeyArray * >	m_vhka

bool	m_varHandleArraysDeclared

Detailed Description

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

◆ gFEXJwoJAlgo()

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

Constructors.

Definition at line 22 of file gFEXJwoJAlgo.cxx.

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

Member Function Documentation

◆ cosLUT()

float LVL1::gFEXJwoJAlgo::cosLUT	(	unsigned int	phiIDX,
		unsigned int	aw
	)		const

private

Definition at line 541 of file gFEXJwoJAlgo.cxx.

 {
   float c = static_cast<float>(phiIDX)/std::pow(2,aw);
   float rad = (2*M_PI) *c;
   float rcos = std::cos(rad);
   int icos = std::round(rcos*(std::pow(2,aw) - 1));
   return icos;
 }

◆ declareGaudiProperty() [1/4]

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

inlineprivateinherited

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

Definition at line 170 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	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());
  
   }

◆ declareGaudiProperty() [3/4]

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

inlineprivateinherited

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

Definition at line 184 of file AthCommonDataStore.h.

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

◆ declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty	(	Gaudi::Property< T > &	t,
		const SG::NotHandleType &
	)

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(t);
   }

◆ declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleBase &	hndl,
		const std::string &	doc,
		const SG::VarHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

   {
     this->declare(hndl.vhKey());
     hndl.vhKey().setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKey &	hndl,
		const std::string &	doc,
		const SG::VarHandleKeyType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

   {
     this->declare(hndl);
     hndl.setOwner(this);
  
     return PBASE::declareProperty(name,hndl,doc);
   }

◆ declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		SG::VarHandleKeyArray &	hndArr,
		const std::string &	doc,
		const SG::VarHandleKeyArrayType &
	)

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

   {
  
     // std::ostringstream ost;
     // ost << Algorithm::name() << " VHKA declareProp: " << name 
     //     << " size: " << hndArr.keys().size() 
     //     << " mode: " << hndArr.mode() 
     //     << "  vhka size: " << m_vhka.size()
     //     << "\n";
     // debug() << ost.str() << endmsg;
  
     hndArr.setOwner(this);
     m_vhka.push_back(&hndArr);
  
     Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
     if (p != 0) {
       p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
     } else {
       ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                     << name);
     }
  
     return p;
  
   }

◆ declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc,
		const SG::NotHandleType &
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

   {
     return PBASE::declareProperty(name, property, doc);
   }

◆ declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty	(	const std::string &	name,
		T &	property,
		const std::string &	doc = `"none"`
	)

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

   {
     typedef typename SG::HandleClassifier<T>::type htype;
     return declareProperty (name, property, doc, htype());
   }

◆ declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > & 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; }

◆ etFPGA()

void LVL1::gFEXJwoJAlgo::etFPGA	(	int	FPGAnum,
		const gTowersType &	twrs,
		gTowersType &	gBlkSum,
		int	gBlockthreshold,
		int	A,
		int	B,
		int &	eth,
		int &	ets,
		int &	etw
	)		const

private

Definition at line 423 of file gFEXJwoJAlgo.cxx.

                                                                                                  {
  
   gBlockthreshold = gBlockthreshold * 200 / 800; //gBlockthreshold is provided in counts with a resolution of 200 MeV, but here needs to be applied with a resolution of 800 GeV
  
   int64_t ethard_hi = 0;
   int64_t etsoft_hi = 0;
   int64_t ethard_lo = 0;
   int64_t etsoft_lo = 0;
  
   int64_t ethard = 0.0;
   int64_t etsoft = 0.0; 
  
   int multiplicitiveFactor = 0;
  
   if(FPGAnum < 2 ) {
    multiplicitiveFactor = cosLUT(0, 5);
   } else{
     multiplicitiveFactor = cosLUT(1, 5);
   }
  
 // firmware treats upper and lower columns differnetly 
  
   for( int irow = 0; irow < FEXAlgoSpaceDefs::ABCrows; irow++ ){
     for(int jcolumn = 0; jcolumn<6; jcolumn++){
         if(gBlkSum[irow][jcolumn] > gBlockthreshold){
       ethard_lo = ethard_lo + twrs[irow][jcolumn]*multiplicitiveFactor; 
     } else {
       etsoft_lo = etsoft_lo + twrs[irow][jcolumn]*multiplicitiveFactor; 
     }
     }
   }
  
   for( int irow = 0; irow < FEXAlgoSpaceDefs::ABCrows; irow++ ){
     for(int jcolumn = 6; jcolumn<12; jcolumn++){
         if(gBlkSum[irow][jcolumn] > gBlockthreshold){
       ethard_hi = ethard_hi + twrs[irow][jcolumn]*multiplicitiveFactor; 
     } else {
       etsoft_hi = etsoft_hi + twrs[irow][jcolumn]*multiplicitiveFactor; 
     }
     }
   }
  
   ethard = ethard_hi + ethard_lo;
   
   etsoft = etsoft_hi + etsoft_lo;
  
  
   int64_t etsum_hi = ethard_hi*A  + etsoft_hi*B  ;
   if ( etsum_hi < 0 ) etsum_hi = 0; 
  
   int64_t etsum_lo = ethard_lo*A  + etsoft_lo*B  ;
   if ( etsum_lo < 0 ) etsum_lo = 0; 
   
   int64_t etsum = etsum_hi + etsum_lo; 
  
  
   // convert 200 MeV LSB here 
   eth  = ethard>>3;
   ets  = etsoft>>3;
   etw  = (etsum  >>13 ) ;
  
   if(msgLvl(MSG::DEBUG)) { 
     std::cout << "DMS FPGA gTEJWOJ " << std::hex <<  FPGAnum << "et sum hard " << eth << "etsum soft" << ets << " A " << A << " B " << B << " weighted term " << etw << std::endl << std::dec; 
   }
 }

◆ etTotal()

void LVL1::gFEXJwoJAlgo::etTotal	(	int	A_ET,
		int	B_ET,
		int	C_ET,
		int &	ET
	)		const

private

Definition at line 512 of file gFEXJwoJAlgo.cxx.

                                              {
  
   //  leave at 200 MeV scale    
   ET = (A_ET + B_ET + C_ET) ; 
  
  
   // main vlaue of ET is always positive 
   if( ET > 0x0000FFF) ET =  0x0000FFF; 
  
 }

◆ evtStore() [1/2]

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

◆ evtStore() [2/2]

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

inlineinherited

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

Definition at line 90 of file AthCommonDataStore.h.

90 { return m_evtStore; }

◆ 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

◆ gBlockAB()

void LVL1::gFEXJwoJAlgo::gBlockAB	(	const gTowersType &	twrs,
		gTowersType &	gBlkSum,
		gTowersType &	hasSeed,
		int	seedThreshold
	)		const

private

Definition at line 280 of file gFEXJwoJAlgo.cxx.

                                                                                                                            {
  
   int rows = twrs.size();
   int cols = twrs[0].size();
   for( int irow = 0; irow < rows; irow++ ){
     for(int jcolumn = 0; jcolumn<cols; jcolumn++){
       // zero jet sum here
       gBlkSum[irow][jcolumn] = 0;
       int krowUp = (irow + 1)%32;
       int krowDn = (irow - 1 +32)%32;
       if( (jcolumn == 0) || (jcolumn == 6) ) {
         //left edge case
         gBlkSum[irow][jcolumn] =
         twrs[irow][jcolumn]   + twrs[krowUp][jcolumn]   + twrs[krowDn][jcolumn] +
         twrs[irow][jcolumn+1] + twrs[krowUp][jcolumn+1] + twrs[krowDn][jcolumn+1];
       } else if( (jcolumn == 5) || (jcolumn == 11)) {
         //  right edge case
         gBlkSum[irow][jcolumn] =
         twrs[irow][jcolumn]   + twrs[krowUp][jcolumn]   + twrs[krowDn][jcolumn] +
         twrs[irow][jcolumn-1] + twrs[krowUp][jcolumn-1] + twrs[krowDn][jcolumn-1];
       } else{
         // normal case
         gBlkSum[irow][jcolumn] =
         twrs[irow][jcolumn]   + twrs[krowUp][jcolumn]   + twrs[krowDn][jcolumn]   +
         twrs[irow][jcolumn-1] + twrs[krowUp][jcolumn-1] + twrs[krowDn][jcolumn-1] +
         twrs[irow][jcolumn+1] + twrs[krowUp][jcolumn+1] + twrs[krowDn][jcolumn+1];
       }
  
       if( gBlkSum[irow][jcolumn]  > seedThreshold) {
         hasSeed[irow][jcolumn] = 1;
       } else {
         hasSeed[irow][jcolumn] = 0;
       }
     
       if ( gBlkSum[irow][jcolumn] < 0 )       
         gBlkSum[irow][jcolumn] = 0;
  
       // was bits 11+3 downto 3, now is 11 downto 0 
       if ( gBlkSum[irow][jcolumn] > FEXAlgoSpaceDefs::gBlockMax )   {
         gBlkSum[irow][jcolumn] =  FEXAlgoSpaceDefs::gBlockMax;
       }
     }
   }
 }

◆ initialize()

StatusCode LVL1::gFEXJwoJAlgo::initialize ( )

overridevirtual

standard Athena-Algorithm method

Definition at line 29 of file gFEXJwoJAlgo.cxx.

                                    {
  
   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::IgFEXJwoJAlgo::interfaceID ( )

inlinestaticinherited

Definition at line 39 of file IgFEXJwoJAlgo.h.

   {
     return IID_IgFEXJwoJAlgo;
   }

◆ jwojAlgo()

std::vector< std::unique_ptr< gFEXJwoJTOB > > LVL1::gFEXJwoJAlgo::jwojAlgo	(	const gTowersType &	Atwr,
		const gTowersType &	Btwr,
		const gTowersType &	Ctwr,
		std::array< uint32_t, 4 > &	outTOB
	)		const

overridevirtual

Implements LVL1::IgFEXJwoJAlgo.

Definition at line 54 of file gFEXJwoJAlgo.cxx.

                                                                                                        {
  
  
   // input towers have 200 MeV LSB
  
   // find gBlocks
   gTowersType AgBlk;
   gTowersType Ascaled;
  
   gTowersType BgBlk;
   gTowersType Bscaled;
  
   gTowersType CgBlk;
   gTowersType Cscaled;
  
   gTowersType hasSeed;
  
   gBlockAB(Atwr, AgBlk, hasSeed, m_gBlockthresholdA);
   gBlockAB(Btwr, BgBlk, hasSeed, m_gBlockthresholdB);
   gBlockAB(Ctwr, CgBlk, hasSeed, m_gBlockthresholdC);
  
  
   // switch to 10 bit number
   // DMS -- do we eventaully need to check for overflows here? 
   
   for(int irow = 0; irow<FEXAlgoSpaceDefs::ABCrows; irow++){
     for(int jcolumn = 0; jcolumn<FEXAlgoSpaceDefs::ABcolumns; jcolumn++){
       Ascaled[irow][jcolumn] = Atwr[irow][jcolumn] >> 2;
       AgBlk[irow][jcolumn] = AgBlk[irow][jcolumn] >> 2;
  
       Bscaled[irow][jcolumn] = Btwr[irow][jcolumn] >> 2;
       BgBlk[irow][jcolumn] = BgBlk[irow][jcolumn] >> 2;
       
       Cscaled[irow][jcolumn] = Ctwr[irow][jcolumn] >> 2;
       CgBlk[irow][jcolumn] = CgBlk[irow][jcolumn] >> 2;
  
     }
   }
  
  
  //FPGA A observables
   int A_MHT_x = 0x0;
   int A_MHT_y = 0x0;
   int A_MST_x = 0x0;
   int A_MST_y = 0x0;
   int A_MET_x = 0x0;
   int A_MET_y = 0x0;
  
   int A_eth = 0x0;
   int A_ets = 0x0;
   int A_etw = 0x0;
  
   //FPGA B observables
   int B_MHT_x = 0x0;
   int B_MHT_y = 0x0;
   int B_MST_x = 0x0;
   int B_MST_y = 0x0;
   int B_MET_x = 0x0;
   int B_MET_y = 0x0;
  
   int B_eth = 0x0;
   int B_ets = 0x0;
   int B_etw = 0x0;
  
  
   //FPGA C observables
   int C_MHT_x = 0x0;
   int C_MHT_y = 0x0;
   int C_MST_x = 0x0;
   int C_MST_y = 0x0;
   int C_MET_x = 0x0;
   int C_MET_y = 0x0;
  
   int C_eth = 0x0;
   int C_ets = 0x0;
   int C_etw = 0x0;
  
  
   //Global observables
   int MHT_x = 0x0;
   int MHT_y = 0x0;
   int MST_x = 0x0;
   int MST_y = 0x0;
   int MET_x = 0x0;
   int MET_y = 0x0;
  
   int ETH =  0x0; 
   int ETS =  0x0;
   int ETW =  0x0;
  
   int total_sumEt = 0x0;
   int MET = 0x0; 
  
  
   // will need to hard code etFPGA ,a's and b's 
   int etBprime =0;
  
  
   metFPGA(0, Ascaled, AgBlk, m_gBlockthresholdA, m_aFPGA_A, m_bFPGA_A, A_MHT_x, A_MHT_y, A_MST_x, A_MST_y, A_MET_x, A_MET_y);
   etFPGA (0, Ascaled, AgBlk, m_gBlockthresholdA, m_aFPGA_A, etBprime, A_eth, A_ets, A_etw); 
  
   metFPGA(1, Bscaled, BgBlk, m_gBlockthresholdB, m_aFPGA_B, m_bFPGA_B, B_MHT_x, B_MHT_y, B_MST_x, B_MST_y, B_MET_x, B_MET_y);
   etFPGA (1, Bscaled, BgBlk, m_gBlockthresholdB, m_aFPGA_B, etBprime, B_eth, B_ets, B_etw); 
  
   metFPGA(2, Cscaled, CgBlk, m_gBlockthresholdC, m_aFPGA_C, m_bFPGA_C, C_MHT_x, C_MHT_y, C_MST_x, C_MST_y, C_MET_x, C_MET_y);
   etFPGA (2, Cscaled, CgBlk, m_gBlockthresholdC, m_aFPGA_C, etBprime, C_eth, C_ets, C_etw); 
  
   metTotal(A_MHT_x, A_MHT_y, B_MHT_x, B_MHT_y, C_MHT_x, C_MHT_y, MHT_x, MHT_y);
   metTotal(A_MST_x, A_MST_y, B_MST_x, B_MST_y, C_MST_x, C_MST_y, MST_x, MST_y);
   metTotal(A_MET_x, A_MET_y, B_MET_x, B_MET_y, C_MET_x, C_MET_y, MET_x, MET_y);
  
   etTotal(A_eth, B_eth, C_eth, ETH);
   etTotal(A_ets, B_ets, C_ets, ETS);
   etTotal(A_etw, B_etw, C_etw, ETW);
  
   total_sumEt = ETW;      
  
   // components should all be less than 12 bits at this point with 200 MeV LSB
   int MET2 = MET_x * MET_x + MET_y * MET_y;
  
   if (MET2 > 0x0FFFFFF) {
       MET = 0x000FFF;
   } else {
       // repeat the byte stream converter calculation here -- not what the hardware actually does
       MET = std::sqrt(MET2); 
  
       
           // best guess at current hardware.  Note that this is computed in the bytestream converter      
       // take most signficant 12 bits 
       //int MET12 = MET2 >> 12; 
       // simulate the look up -- only 6 most signficant bits currently set -- to be checked 
       //MET = ( (int)(std::sqrt(MET12)) << 6) &  0x00000FC0   ;
   }
  
  
   //Define a vector to be filled with all the TOBs of one event
   std::vector<std::unique_ptr<gFEXJwoJTOB>> tobs_v;
   tobs_v.resize(4);
  
  
   // fill in TOBs
   // The order of the TOBs is given according to the TOB ID (TODO: check how it's done in fw)
  
   // First TOB is (MET, SumEt)
   outTOB[0] = (total_sumEt &  0x00000FFF) << 0; //set the Quantity2 to the corresponding slot (LSB)
   outTOB[0] = outTOB[0] | (MET  &  0x00000FFF) << 12;//Quantity 1 (in bit number 12)
   if (total_sumEt != 0) outTOB[0] = outTOB[0] | 0x00000001 << 24;//Status bit for Quantity 2 (0 if quantity is null)
   if (MET != 0) outTOB[0] = outTOB[0] | 0x00000001 << 25;//Status bit for Quantity 1 (0 if quantity is null)
   outTOB[0] = outTOB[0] | (1  &  0x0000001F) << 26;//TOB ID is 1 for scalar values (5 bits starting at 26)
  
   // std::cout << "DMS MET " << std::hex << MET  << " total_sumEt "  << total_sumEt <<  std::endl << std::dec;
  
 // Second TOB is (MET_x, MET_y)
   outTOB[1] = (MET_y &  0x00000FFF) << 0; //set the Quantity2 to the corresponding slot (LSB)
   outTOB[1] = outTOB[1] | (MET_x  &  0x00000FFF) << 12;//Quantity 1 (in bit number 12)
   if (MET_y != 0) outTOB[1] = outTOB[1] | 0x00000001 << 24;//Status bit for Quantity 2 (0 if quantity is null)
   if (MET_x != 0) outTOB[1] = outTOB[1] | 0x00000001 << 25;//Status bit for Quantity 1 (0 if quantity is null)
   outTOB[1] = outTOB[1] | (2  &  0x0000001F) << 26;//TOB ID is 2 for MET_x, MET_y (5 bits starting at 26)
  
 // Third TOB is hard components (MHT_x, MHT_y)
   outTOB[2] = (MHT_y &  0x00000FFF) << 0; //set the Quantity2 to the corresponding slot (LSB)
   outTOB[2] = outTOB[2] | (MHT_x  &  0x00000FFF) << 12;//Quantity 1 (in bit number 12)
   if (MHT_y != 0) outTOB[2] = outTOB[2] | 0x00000001 << 24;//Status bit for Quantity 2 (0 if quantity is null)
   if (MHT_x != 0) outTOB[2] = outTOB[2] | 0x00000001 << 25;//Status bit for Quantity 1 (0 if quantity is null)
   outTOB[2] = outTOB[2] | (3  &  0x0000001F) << 26;//TOB ID is 3 for hard components (5 bits starting at 26)
  
   // Fourth TOB is hard components (MST_x, MST_y)
   outTOB[3] = (MST_y &  0x00000FFF) << 0; //set the Quantity2 to the corresponding slot (LSB)
   outTOB[3] = outTOB[3] | (MST_x  &  0x00000FFF) << 12;//Quantity 1 (in bit number 12)
   if (MST_y != 0) outTOB[3] = outTOB[3] | 0x00000001 << 24;//Status bit for Quantity 2 (0 if quantity is null)
   if (MST_x != 0) outTOB[3] = outTOB[3] | 0x00000001 << 25;//Status bit for Quantity 1 (0 if quantity is null)
   outTOB[3] = outTOB[3] | (4  &  0x0000001F) << 26;//TOB ID is 4 for soft components (5 bits starting at 26)
  
  
   tobs_v[0] = std::make_unique<gFEXJwoJTOB>();
   tobs_v[0]->setWord(outTOB[0]);
   tobs_v[0]->setQuantity1(MET);
   tobs_v[0]->setQuantity2(total_sumEt);
   tobs_v[0]->setSaturation(0); //Always 0 for now, need a threshold?
   tobs_v[0]->setTobID(1);
   if( MET != 0 ) tobs_v[0]->setStatus1(1);
   else tobs_v[0]->setStatus1(0);
   if(total_sumEt!= 0) tobs_v[0]->setStatus2(1);
   else tobs_v[0]->setStatus2(0);
  
   tobs_v[1] = std::make_unique<gFEXJwoJTOB>();
   tobs_v[1]->setWord(outTOB[1]);
   tobs_v[1]->setQuantity1(MET_x);
   tobs_v[1]->setQuantity2(MET_y);
   tobs_v[1]->setSaturation(0); //Always 0 for now, need a threshold?
   tobs_v[1]->setTobID(2);
   if( MET_x != 0 ) tobs_v[1]->setStatus1(1);
   else tobs_v[1]->setStatus1(0);
   if(MET_y!= 0) tobs_v[1]->setStatus2(1);
   else tobs_v[1]->setStatus2(0);
  
   tobs_v[2] = std::make_unique<gFEXJwoJTOB>();
   tobs_v[2]->setWord(outTOB[2]);
   tobs_v[2]->setQuantity1(MHT_x);
   tobs_v[2]->setQuantity2(MHT_y);
   tobs_v[2]->setSaturation(0); //Always 0 for now, need a threshold?
   tobs_v[2]->setTobID(3);
   if( MHT_x != 0 ) tobs_v[2]->setStatus1(1);
   else tobs_v[2]->setStatus1(0);
   if(MHT_y!= 0) tobs_v[2]->setStatus2(1);
   else tobs_v[2]->setStatus2(0);
  
   tobs_v[3] = std::make_unique<gFEXJwoJTOB>();
   tobs_v[3]->setWord(outTOB[3]);
   tobs_v[3]->setQuantity1(MST_x);
   tobs_v[3]->setQuantity2(MST_y);
   tobs_v[3]->setSaturation(0); //Always 0 for now, need a threshold?
   tobs_v[3]->setTobID(4);
   if( MST_x != 0 ) tobs_v[3]->setStatus1(1);
   else tobs_v[2]->setStatus1(0);
   if(MST_y!= 0) tobs_v[3]->setStatus2(1);
   else tobs_v[3]->setStatus2(0);
  
  
   return tobs_v;
  
 }

◆ metFPGA()

void LVL1::gFEXJwoJAlgo::metFPGA	(	int	FPGAnum,
		const gTowersType &	twrs,
		const gTowersType &	gBlkSum,
		int	gBlockthreshold,
		int	aFPGA,
		int	bFPGA,
		int &	MHT_x,
		int &	MHT_y,
		int &	MST_x,
		int &	MST_y,
		int &	MET_x,
		int &	MET_y
	)		const

private

Definition at line 326 of file gFEXJwoJAlgo.cxx.

                                                            {
   
   gBlockthreshold = gBlockthreshold * 200 / 800; //gBlockthreshold is provided in counts with a resolution of 200 MeV, but here needs to be applied with a resolution of 800 GeV
   // in the RTL  code these are 19+ 5 = 24 bits 
   int64_t h_tx_hi = 0;
   int64_t h_ty_hi = 0;
   int64_t h_tx_lw = 0;
   int64_t h_ty_lw = 0;
   
   int64_t e_tx_hi = 0;
   int64_t e_ty_hi = 0;
   int64_t e_tx_lw = 0;
   int64_t e_ty_lw = 0;   
   
   for( int irow = 0; irow < FEXAlgoSpaceDefs::ABCrows; irow++ ){
     for(int jcolumn = 6; jcolumn<12; jcolumn++){
       if( FPGAnum == 2){
         int frow = 2*(irow/2)  + 1; 
         
         if(gBlkSum[irow][jcolumn] > gBlockthreshold){
           h_tx_hi += (twrs[irow][jcolumn])*(cosLUT(frow, 5));
           h_ty_hi += (twrs[irow][jcolumn])*(sinLUT(frow, 5));
         } else {
           e_tx_hi += (twrs[irow][jcolumn])*(cosLUT(frow, 5));
           e_ty_hi += (twrs[irow][jcolumn])*(sinLUT(frow, 5));
         }
  
       } else {
   
         if(gBlkSum[irow][jcolumn] > gBlockthreshold){
           h_tx_hi += (twrs[irow][jcolumn])*(cosLUT(irow, 5));
           h_ty_hi += (twrs[irow][jcolumn])*(sinLUT(irow, 5));
         } else {
           e_tx_hi += (twrs[irow][jcolumn])*(cosLUT(irow, 5));
           e_ty_hi += (twrs[irow][jcolumn])*(sinLUT(irow, 5));
         }
       }
     }
      
     for(int jcolumn = 0; jcolumn<6; jcolumn++){
       if( FPGAnum == 2){
         int frow = 2*(irow/2)  + 1;
         
         if(gBlkSum[irow][jcolumn] > gBlockthreshold){
           h_tx_lw += (twrs[irow][jcolumn])*(cosLUT(frow, 5));
           h_ty_lw += (twrs[irow][jcolumn])*(sinLUT(frow, 5));
         } else{
           e_tx_lw += (twrs[irow][jcolumn])*(cosLUT(frow, 5));
           e_ty_lw += (twrs[irow][jcolumn])*(sinLUT(frow, 5));
         }
       } else {
   
         if(gBlkSum[irow][jcolumn] > gBlockthreshold){
           h_tx_lw += (twrs[irow][jcolumn])*(cosLUT(irow, 5));
           h_ty_lw += (twrs[irow][jcolumn])*(sinLUT(irow, 5));
         } else {
           e_tx_lw += (twrs[irow][jcolumn])*(cosLUT(irow, 5));
           e_ty_lw += (twrs[irow][jcolumn])*(sinLUT(irow, 5));
         }
       }
     }
   }
  
   // According to https://gitlab.cern.ch/atlas-l1calo/gfex/firmware/-/issues/406#note_5662344
   // there is no division -- so LSB is indeed 25 MeV
  
   //but later changed to 200 MeV so divide by 8
  
   long int fMHT_x =  (h_tx_hi + h_tx_lw) ;
   long int fMHT_y =  (h_ty_hi + h_ty_lw) ;
   long int fMST_x =  (e_tx_hi + e_tx_lw) ;
   long int fMST_y =  (e_ty_hi + e_ty_lw) ;
     
   MHT_x =  (h_tx_hi + h_tx_lw) >> 3;
   MHT_y =  (h_ty_hi + h_ty_lw) >> 3;
   MST_x =  (e_tx_hi + e_tx_lw) >> 3;
   MST_y =  (e_ty_hi + e_ty_lw) >> 3;
  
   //  a and b coffecients are 10 bits
   // multiplication has an addtional 2^10
   // constant JWJ_OW        : integer := 35;--Out width
   // values are 35 bits long and top 16 bits are taken -- so divide by 2^19
   //  2^10/2^19 = 1/2^9 = 1/512
    
   long int fMET_x = ( aFPGA * (fMHT_x) + bFPGA * (fMST_x) ) >> 13 ;
   long int fMET_y = ( aFPGA * (fMHT_y) + bFPGA * (fMST_y) ) >> 13 ;
  
   MET_x  = fMET_x;
   MET_y  = fMET_y;
     
 }

◆ metTotal()

void LVL1::gFEXJwoJAlgo::metTotal	(	int	A_MET_x,
		int	A_MET_y,
		int	B_MET_x,
		int	B_MET_y,
		int	C_MET_x,
		int	C_MET_y,
		int &	MET_x,
		int &	MET_y
	)		const

private

Definition at line 491 of file gFEXJwoJAlgo.cxx.

                                                             {
  
   
   MET_x = A_MET_x + B_MET_x + C_MET_x;
   MET_y = A_MET_y + B_MET_y+ C_MET_y;
  
   // Truncation of the result, as the individual quantities are 16 bits, while the TOB field is 12 bits
   // MET_x = MET_x >> 4;
   // MET_y = MET_y >> 4;
  
   if (MET_x < -0x008000) MET_x = -0x08000;
   if (MET_y < -0x008000) MET_y = -0x00800;
  
   if (MET_x > 0x0007FF) MET_x  = 0x0007FF;
   if (MET_y > 0x0007FF) MET_y  = 0x0007FF;
  
 }

◆ msg() [1/2]

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

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

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

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

◆ setAlgoConstant()

void LVL1::gFEXJwoJAlgo::setAlgoConstant	(	int	aFPGA_A,
		int	bFPGA_A,
		int	aFPGA_B,
		int	bFPGA_B,
		int	aFPGA_C,
		int	bFPGA_C,
		int	gXE_seedThrA,
		int	gXE_seedThrB,
		int	gXE_seedThrC
	)

overridevirtual

Implements LVL1::IgFEXJwoJAlgo.

Definition at line 36 of file gFEXJwoJAlgo.cxx.

                                                                                          {
   m_aFPGA_A = aFPGA_A;
   m_bFPGA_A = bFPGA_A;
   m_aFPGA_B = aFPGA_B;
   m_bFPGA_B = bFPGA_B;
   m_aFPGA_C = aFPGA_C;
   m_bFPGA_C = bFPGA_C;
   m_gBlockthresholdA = gXE_seedThrA;
   m_gBlockthresholdB = gXE_seedThrB;
   m_gBlockthresholdC = gXE_seedThrC;
  
 }

◆ sinLUT()

float LVL1::gFEXJwoJAlgo::sinLUT	(	unsigned int	phiIDX,
		unsigned int	aw
	)		const

private

Definition at line 528 of file gFEXJwoJAlgo.cxx.

 {
   float c = static_cast<float>(phiIDX)/std::pow(2,aw);
   float rad = (2*M_PI) *c;
   float rsin = std::sin(rad);
   int isin = std::round(rsin*(std::pow(2,aw) - 1));
   return isin;
  
 }

◆ 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 DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

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

float LVL1::gFEXJwoJAlgo::m_aFPGA_A

private

Definition at line 47 of file gFEXJwoJAlgo.h.

◆ m_aFPGA_B

float LVL1::gFEXJwoJAlgo::m_aFPGA_B

private

Definition at line 49 of file gFEXJwoJAlgo.h.

◆ m_aFPGA_C

float LVL1::gFEXJwoJAlgo::m_aFPGA_C

private

Definition at line 51 of file gFEXJwoJAlgo.h.

◆ m_bFPGA_A

float LVL1::gFEXJwoJAlgo::m_bFPGA_A

private

Definition at line 48 of file gFEXJwoJAlgo.h.

◆ m_bFPGA_B

float LVL1::gFEXJwoJAlgo::m_bFPGA_B

private

Definition at line 50 of file gFEXJwoJAlgo.h.

◆ m_bFPGA_C

float LVL1::gFEXJwoJAlgo::m_bFPGA_C

private

Definition at line 52 of file gFEXJwoJAlgo.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_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_gBlockthresholdA

float LVL1::gFEXJwoJAlgo::m_gBlockthresholdA

private

Definition at line 53 of file gFEXJwoJAlgo.h.

◆ m_gBlockthresholdB

float LVL1::gFEXJwoJAlgo::m_gBlockthresholdB

private

Definition at line 54 of file gFEXJwoJAlgo.h.

◆ m_gBlockthresholdC

float LVL1::gFEXJwoJAlgo::m_gBlockthresholdC

private

Definition at line 55 of file gFEXJwoJAlgo.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

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ gFEXJwoJAlgo()

Member Function Documentation

◆ cosLUT()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ etFPGA()

◆ etTotal()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ gBlockAB()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ jwojAlgo()

◆ metFPGA()

◆ metTotal()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ setAlgoConstant()

◆ sinLUT()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_aFPGA_A

◆ m_aFPGA_B

◆ m_aFPGA_C

◆ m_bFPGA_A

◆ m_bFPGA_B

◆ m_bFPGA_C

◆ m_detStore

◆ m_evtStore

◆ m_gBlockthresholdA

◆ m_gBlockthresholdB

◆ m_gBlockthresholdC

◆ m_varHandleArraysDeclared

◆ m_vhka