#include <eFexTowerBuilder.h>

Inheritance diagram for LVL1::eFexTowerBuilder:

Collaboration diagram for LVL1::eFexTowerBuilder:

Public Member Functions
	eFexTowerBuilder (const std::string &name, ISvcLocator *pSvcLocator)

	~eFexTowerBuilder ()=default

virtual StatusCode	initialize ()

virtual StatusCode	execute (const EventContext &ctx) const

virtual StatusCode	sysInitialize () override
	Override sysInitialize. More...

virtual bool	isClonable () const override
	Specify if the algorithm is clonable. More...

virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant. More...

virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm. More...

virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies. More...

virtual bool	filterPassed (const EventContext &ctx) const

virtual void	setFilterPassed (bool state, const EventContext &ctx) const

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

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
StatusCode	fillTowers (const EventContext &ctx) const

StatusCode	fillMap (const EventContext &ctx) 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
std::mutex m_fillMapMutex	ATLAS_THREAD_SAFE

std::map< unsigned long long, std::pair< std::pair< int, int >, std::pair< int, int > > > m_scMap	ATLAS_THREAD_SAFE

SG::ReadHandleKey< xAOD::EventInfo >	m_eiKey {this,"EventInfoKey","EventInfo",""}

SG::ReadCondHandleKey< CaloSuperCellDetDescrManager >	m_ddmKey {this,"CaloSuperCellDetDescrManager","CaloSuperCellDetDescrManager",""}

SG::ReadHandleKey< CaloCellContainer >	m_scellKey { this, "CaloCellContainerReadKey", "SCell", "Read handle key for the supercells"}

SG::ReadHandleKey< xAOD::TriggerTowerContainer >	m_ttKey { this, "TriggerTowerContainerReadKey", "xAODTriggerTowers", "Read handle key for the triggerTowers"}

SG::WriteHandleKey< xAOD::eFexTowerContainer >	m_outKey {this, "eFexContainerWriteKey", "L1_eFexEmulatedTowers", "Name of the output container"}

Gaudi::Property< std::string >	m_mappingFile {this, "MappingFile", "L1CaloFEXByteStream/2023-02-13/scToEfexTowers.root", "PathResolver location to mapping file"}

ToolHandle< IeFEXSuperCellTowerIdProvider >	m_eFEXSuperCellTowerIdProviderTool {this, "eFEXSuperCellTowerIdProviderTool", "LVL1::eFEXSuperCellTowerIdProvider", "Tool that provides tower-FOGA mapping"}

Gaudi::Property< bool >	m_applyMasking {this,"ApplyMasking",true,"Apply masking of supercells based on provenance bits. Should be set to False for MC"}

DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks. More...

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 48 of file eFexTowerBuilder.h.

Member Typedef Documentation

◆ StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ eFexTowerBuilder()

LVL1::eFexTowerBuilder::eFexTowerBuilder	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 28 of file eFexTowerBuilder.cxx.

                                                                                     : AthReentrantAlgorithm( name, pSvcLocator ){
  
 }

◆ ~eFexTowerBuilder()

LVL1::eFexTowerBuilder::~eFexTowerBuilder ( )

default

Member Function Documentation

◆ cardinality()

unsigned int AthReentrantAlgorithm::cardinality ( ) const

overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 55 of file AthReentrantAlgorithm.cxx.

 {
   return 0;
 }

◆ declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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; }

◆ evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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; }

◆ execute()

StatusCode LVL1::eFexTowerBuilder::execute ( const EventContext & ctx ) const

virtual

Definition at line 341 of file eFexTowerBuilder.cxx.

                                                                   {
     ATH_MSG_DEBUG("Executing " << name() << "...");
     setFilterPassed(true, ctx);
  
  
     {
         std::lock_guard lock(m_fillMapMutex);
         if (m_scMap.empty()) CHECK( fillMap(ctx) );
     }
  
     return fillTowers(ctx);
  
 }

◆ extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

◆ extraOutputDeps()

const DataObjIDColl & AthReentrantAlgorithm::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 79 of file AthReentrantAlgorithm.cxx.

 {
   // If we didn't find any symlinks to add, just return the collection
   // from the base class.  Otherwise, return the extended collection.
   if (!m_extendedExtraObjects.empty()) {
     return m_extendedExtraObjects;
   }
   return Algorithm::extraOutputDeps();
 }

◆ fillMap()

StatusCode LVL1::eFexTowerBuilder::fillMap ( const EventContext & ctx ) const

private

Definition at line 197 of file eFexTowerBuilder.cxx.

                                                                   {
  
     ATH_MSG_INFO("Filling sc -> eFexTower map");
  
     SG::ReadCondHandle<CaloSuperCellDetDescrManager> ddm{m_ddmKey,ctx};
     SG::ReadHandle<CaloCellContainer> scells(m_scellKey,ctx); // 34048 is a full complement of scells
     if(!scells.isValid()){
         ATH_MSG_FATAL("Could not retrieve collection " << m_scellKey.key() );
         return StatusCode::FAILURE;
     }
     if (scells->size() != 34048) {
         ATH_MSG_FATAL("Cannot fill sc -> eFexTower mapping with an incomplete sc collection");
         return StatusCode::FAILURE;
     }
     struct TowerSCells {
         std::vector<unsigned long long> ps;
         std::vector<std::pair<float,unsigned long long>> l1;
         std::vector<std::pair<float,unsigned long long>> l2;
         std::vector<unsigned long long> l3;
         std::vector<unsigned long long> had;
         std::vector<unsigned long long> other;
     };
     static const auto etaIndex = [](float eta) { return int( eta*10 ) + ((eta<0) ? -1 : 1); }; // runs from -25 to 25, skipping over 0 (so gives outer edge eta)
     static const auto phiIndex = [](float phi) { return int( phi*32./ROOT::Math::Pi() ) + (phi<0 ? -1 : 1); }; // runs from -pi to pi, skipping over 0 (gives out edge phi)
     std::map<std::pair<int,int>,TowerSCells> towers;
     std::map<unsigned long long,int> eTowerSlots; // not used by this alg, but we produce the map for benefit of eFexTower->eTower alg
  
     for (auto digi: *scells) {
         Identifier id = digi->ID(); // this is if using supercells
  
         if (auto elem = ddm->get_element(id); elem && std::abs(elem->eta_raw())<2.5) {
             float eta = elem->eta_raw(); // this seems more symmetric
             int sampling = elem->getSampling();
             if(sampling==6 && ddm->getCaloCell_ID()->region(id)==0 && eta<0) eta-=0.01; // nudge this L2 endcap supercell into correct tower (right on boundary)
  
             unsigned long long val = id.get_compact();
  
             int towerid = -1;int slot = -1;bool issplit = false;
             CHECK(m_eFEXSuperCellTowerIdProviderTool->geteTowerIDandslot(id.get_compact(), towerid, slot, issplit));
             eTowerSlots[id.get_compact()] = slot;
  
             auto& sc = towers[std::pair(etaIndex(eta),phiIndex(elem->phi_raw()))];
             switch(sampling) {
                 case 0: case 4: //lar barrel/endcap presampler
                     sc.ps.push_back(val);
                     break;
                 case 1: case 5: //lar barrel/endcap l1
                     sc.l1.push_back({elem->eta(),val}); break;
                 case 2: case 6: //lar barrel/endcap l2
                     sc.l2.push_back({elem->eta(),val}); break;
                 case 3: case 7: //lar barrel/endcap l3
                     sc.l3.push_back(val); break;
                 case 8: case 9: case 10: case 11: //lar hec
                     sc.had.push_back(val); break;
                 default:
                     sc.other.push_back(val); break;
             }
         }
     }
  
  
     // sort (by increasing eta) l1/l2 sc and handle special cases
     // finally also output the eTower slot vector
     std::vector<size_t> slotVector(11);
     for(auto& [coord,sc] : towers) {
         std::sort(sc.l1.begin(),sc.l1.end());
         std::sort(sc.l2.begin(),sc.l2.end());
         // we have 5 l2 cells @ |eta|=1.45 ... put lowest |eta| one in l3 slot
         if (sc.l2.size()==5) {
             if (coord.first >= 0) {
                 sc.l3.push_back(sc.l2.front().second);
                 sc.l2.erase(sc.l2.begin()); // remove first
             } else {
                 sc.l3.push_back(sc.l2.back().second);
                 sc.l2.resize(sc.l2.size()-1); // remove last
             }
         }
         if (std::abs(coord.first)==15) { //|eta| = 1.45
             // in the overlap region it seems like the latome id with highest |eta| is swapped with next highest
             // so to compare we swap the first and second (3rd and 4th are fine) if eta < 0, or 3rd and 4th if eta > 0
             if (coord.first<0) {std::swap(sc.l1.at(0),sc.l1.at(1)); }
             else {std::swap(sc.l1.at(2),sc.l1.at(3));}
         }
         // handle case @ |eta|~1.8-2 with 6 L1 cells
         if (sc.l1.size()==6) {
             m_scMap[sc.l1.at(0).second] = std::pair(coord,std::pair(1,11));
             m_scMap[sc.l1.at(1).second] = std::pair(coord,std::pair(1,2));
             m_scMap[sc.l1.at(2).second] = std::pair(coord,std::pair(2,11));
             m_scMap[sc.l1.at(3).second] = std::pair(coord,std::pair(3,11));
             m_scMap[sc.l1.at(4).second] = std::pair(coord,std::pair(3,4));
             m_scMap[sc.l1.at(5).second] = std::pair(coord,std::pair(4,11));
             slotVector[1] = eTowerSlots[sc.l1.at(0).second];
             slotVector[2] = eTowerSlots[sc.l1.at(2).second];
             slotVector[3] = eTowerSlots[sc.l1.at(3).second];
             slotVector[4] = eTowerSlots[sc.l1.at(5).second];
         }
  
         // for |eta|>2.4 there's only 1 l1 sc, to match hardware this should be compared placed in the 'last' l1 input
         if (sc.l1.size()==1) {
             m_scMap[sc.l1.at(0).second] = std::pair(coord,std::pair(4,11));
             slotVector[1] = 1; slotVector[2] = 2; slotVector[3] = 3; slotVector[4] = eTowerSlots[sc.l1.at(0).second];
         }
  
         // fill the map with sc ids -> tower coord + slot
         if (!sc.ps.empty()) {m_scMap[sc.ps.at(0)] = std::pair(coord,std::pair(0,11)); slotVector[0] = eTowerSlots[sc.ps.at(0)]; }
         if(sc.l1.size()==4) for(size_t i=0;i<4;i++) if(sc.l1.size() > i) {m_scMap[sc.l1.at(i).second] = std::pair(coord,std::pair(i+1,11)); slotVector[i+1] = eTowerSlots[sc.l1.at(i).second]; }
         for(size_t i=0;i<4;i++) if(sc.l2.size() > i) { m_scMap[sc.l2.at(i).second] = std::pair(coord,std::pair(i+5,11)); slotVector[i+5] = eTowerSlots[sc.l2.at(i).second]; }
         if (!sc.l3.empty()) {m_scMap[sc.l3.at(0)] = std::pair(coord,std::pair(9,11)); slotVector[9] = eTowerSlots[sc.l3.at(0)]; }
         if (!sc.had.empty()) {m_scMap[sc.had.at(0)] = std::pair(coord,std::pair(10,11));slotVector[10] = eTowerSlots[sc.had.at(0)]; }
  
         // finally output the slotVector for this tower
         // do only for the slots that don't match
         // note to self: seems like everything is fine apart from the l1->ps remap for |eta|>2.4
         // so leaving this bit commented out for now ... useful to leave it here in case need to recheck in future
 //        for(size_t i=0;i<slotVector.size();i++) {
 //            if(slotVector[i] != i) {
 //                std::cout << coord.first << "," << coord.second << "," << i << "," << slotVector[i] << std::endl;
 //            }
 //        }
     }
  
     // save the map to disk
     TFile f("scToEfexTowers.root","RECREATE");
     TTree* t = new TTree("mapping","mapping");
     unsigned long long scid = 0;
     std::pair<int,int> coord = {0,0};
     std::pair<int,int> slot = {-1,-1};
     t->Branch("scid",&scid);
     t->Branch("etaIndex",&coord.first);
     t->Branch("phiIndex",&coord.second);
     t->Branch("slot1",&slot.first);
     t->Branch("slot2",&slot.second);
     for(auto& [id,val] : m_scMap) {
         scid = id; coord = val.first; slot = val.second;
         t->Fill();
     }
     t->Write();
     f.Close();
  
     return StatusCode::SUCCESS;
  
 }

◆ fillTowers()

StatusCode LVL1::eFexTowerBuilder::fillTowers ( const EventContext & ctx ) const

private

Definition at line 71 of file eFexTowerBuilder.cxx.

                                                                      {
  
  
     SG::ReadCondHandle<CaloSuperCellDetDescrManager> ddm{m_ddmKey,ctx};
     SG::ReadHandle<xAOD::TriggerTowerContainer> tTowers(m_ttKey,ctx);
     if(!tTowers.isValid()){
         ATH_MSG_FATAL("Could not retrieve collection " << m_ttKey.key() );
         return StatusCode::FAILURE;
     }
     SG::ReadHandle<CaloCellContainer> scells(m_scellKey,ctx); // n.b. 34048 is a full complement of scells
     if(!scells.isValid()){
         ATH_MSG_FATAL("Could not retrieve collection " << m_scellKey.key() );
         return StatusCode::FAILURE;
     }
  
     SG::ReadHandle<xAOD::EventInfo> ei(m_eiKey,ctx);
     if(!ei.isValid()) {
       ATH_MSG_FATAL("Cannot retrieve eventinfo");
       return StatusCode::FAILURE;
     }
     bool isMC = ei->eventType(xAOD::EventInfo::IS_SIMULATION); // currently only used to decide if should set a saturation code or not
     
  
     std::map<std::pair<int,int>,std::array<int,11>> towers;
  
     constexpr int INVALID_VALUE = -99999; // use this value to indicate invalid
     constexpr int SATURATED_VALUE = std::numeric_limits<int>::max()-1; // use this value to indicate saturation
  
     for (auto digi: *scells) {
         const auto itr = m_scMap.find(digi->ID().get_compact());
         if (itr == m_scMap.end()) { continue; } // not in map so not mapping to a tower
         int val =  std::round(digi->energy()/(12.5*std::cosh(digi->eta()))); // 12.5 is b.c. energy is in units of 12.5MeV per count
         // note: a val of < -99998 is what is produced if efex was sent an invalid code of 1022 (see LArRawtoSuperCell)
         bool isSaturated = (!isMC) ? (digi->quality()) : false; // not applying saturation codes in MC until the changes to trigger counts has been investigated
         bool isMasked = m_applyMasking ? ((digi)->provenance()&0x80) : false;
         bool isInvalid = m_applyMasking ? ((digi)->provenance()&0x40) : false;
         if(isInvalid) {
             val = INVALID_VALUE;
         }
         if(isSaturated) {
             val = SATURATED_VALUE;
         }
  
         auto& tower = towers[itr->second.first];
         if (itr->second.second.second<11) {
             // doing an energy split between slots ... don't include a masked channel (or invalid channel)
             if (!isMasked && val!=INVALID_VALUE) {
                 if(isSaturated) {
                     // mark both as saturated
                     tower.at(itr->second.second.first) = SATURATED_VALUE;
                     tower.at(itr->second.second.second) = SATURATED_VALUE;
                 }
                 if(tower.at(itr->second.second.first)!=(SATURATED_VALUE)) { // don't override saturation
                     // if the other contribution was masked or invalid, revert to 0 before adding this contribution
                     if (tower.at(itr->second.second.first)==std::numeric_limits<int>::max() || tower.at(itr->second.second.first)==INVALID_VALUE) {
                         tower.at(itr->second.second.first)=0;
                     }
                     tower.at(itr->second.second.first) += val >> 1;
                 }
                 if(tower.at(itr->second.second.second)!=(SATURATED_VALUE)) { // don't override saturation
                     // if the other contribution was masked or invalid, revert to 0 before adding this contribution
                     if (tower.at(itr->second.second.second)==std::numeric_limits<int>::max() || tower.at(itr->second.second.second)==INVALID_VALUE) {
                         tower.at(itr->second.second.second)=0;
                     }
                     tower.at(itr->second.second.second) += (val - (val >> 1)); // HW seems fixed now!
                 }
             }
             // hw is incorrectly ignoring masking on the second part
             // so always add the 2nd bit
             //tower.at(itr->second.second.second) += (val - (val >> 1)); // Removed b.c. of fix above - leaving this comment here until resolved!
         } else {
             auto& v = tower.at(itr->second.second.first);
             if (isMasked) {
                 // dont mark it masked if it already has a contribution
                 if(v==0) v = std::numeric_limits<int>::max();
             } else if(isSaturated) {
                 v = val;
             } else {
                 v += val;
             }
         }
  
     }
  
     // add tile energies from TriggerTowers
     static const auto etaIndex = [](float eta) { return int( eta*10 ) + ((eta<0) ? -1 : 1); };
     static const auto phiIndex = [](float phi) { return int( phi*32./M_PI ) + (phi<0 ? -1 : 1); };
     for(const xAOD::TriggerTower_v2* tTower : *tTowers) {
         if (std::abs(tTower->eta()) > 1.5) continue;
         if (tTower->sampling() != 1) continue;
         double phi = tTower->phi(); if(phi > M_PI) phi -= 2.*M_PI;
         towers[std::pair(etaIndex(tTower->eta()),phiIndex(phi))][10] = tTower->cpET();
     }
  
  
     SG::WriteHandle<xAOD::eFexTowerContainer> eTowers = SG::WriteHandle<xAOD::eFexTowerContainer>(m_outKey,ctx);
     ATH_CHECK( eTowers.record(std::make_unique<xAOD::eFexTowerContainer>(),std::make_unique<xAOD::eFexTowerAuxContainer>()) );
  
     static const auto calToFex = [](int calEt) {
         if(calEt == std::numeric_limits<int>::max()) return 0; // indicates masked channel
         if(calEt == SATURATED_VALUE) return 1023; // saturated channel
         if( calEt == INVALID_VALUE ) return 1022; // invalid channel value
         if(calEt<448) return std::max((calEt&~1)/2+32,1); // 25 MeV per eFexTower count
         if(calEt<1472) return (calEt-448)/4+256;          // 50 MeV per eFexTower count
         if(calEt<3520) return (calEt-1472)/8+512;         // 100 MeV ...
         if(calEt<11584) return (calEt-3520)/32+768;       // 400 MeV ...
         return 1020;
     };
  
     // now create the towers
     for(auto& [coord,counts] : towers) {
         size_t ni = (std::abs(coord.first)<=15) ? 10 : 11; // ensures we skip the tile towers for next line
         for(size_t i=0;i<ni;++i) counts[i] = (scells->empty() ? 1025 : calToFex(counts[i])); // do latome energy scaling to non-tile towers - if had no cells will use code "1025" to indicate
         eTowers->push_back( std::make_unique<xAOD::eFexTower>() );
         eTowers->back()->initialize( ( (coord.first<0 ? 0.5:-0.5) + coord.first)*0.1 ,
                                  ( (coord.second<0 ? 0.5:-0.5) + coord.second)*M_PI/32,
                                      std::vector<uint16_t>(counts.begin(), counts.end()),
                                  -1, /* module number */
                                  -1, /* fpga number */
                                  0,0 /* status flags ... could use to indicate which cells were actually present?? */);
     }
  
     return StatusCode::SUCCESS;
  
 }

◆ filterPassed()

virtual bool AthReentrantAlgorithm::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 135 of file AthReentrantAlgorithm.h.

                                                            {
     return execState( ctx ).filterPassed();
   }

◆ initialize()

StatusCode LVL1::eFexTowerBuilder::initialize ( )

virtual

Definition at line 33 of file eFexTowerBuilder.cxx.

                                         {
     ATH_MSG_INFO ("Initializing " << name() << "...");
  
     CHECK( m_ddmKey.initialize(true) );
     CHECK( m_ttKey.initialize(true) );
     CHECK( m_scellKey.initialize(true) );
     CHECK( m_outKey.initialize(true) );
     CHECK( m_eiKey.initialize(true) );
  
     if (auto fileName = PathResolverFindCalibFile( m_mappingFile ); !fileName.empty()) {
         std::unique_ptr<TFile> f( TFile::Open(fileName.c_str()) );
         if (f) {
             TTree* t = f->Get<TTree>("mapping");
             if(t) {
                 unsigned long long scid = 0;
                 std::pair<int,int> coord = {0,0};
                 std::pair<int,int> slot;
                 t->SetBranchAddress("scid",&scid);
                 t->SetBranchAddress("etaIndex",&coord.first);
                 t->SetBranchAddress("phiIndex",&coord.second);
                 t->SetBranchAddress("slot1",&slot.first);
                 t->SetBranchAddress("slot2",&slot.second);
                 for(Long64_t i=0;i<t->GetEntries();i++) {
                     t->GetEntry(i);
                     m_scMap[scid] = std::make_pair(coord,slot);
                 }
             }
         }
         if (m_scMap.empty()) {
             ATH_MSG_WARNING("Failed to load sc -> eFexTower map from " << fileName);
         } else {
             ATH_MSG_INFO("Loaded sc -> eFexTower map from " << fileName);
         }
     }
  
     return StatusCode::SUCCESS;
 }

◆ inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

◆ isClonable()

bool AthReentrantAlgorithm::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in Simulation::BeamEffectsAlg, InDet::SiTrackerSpacePointFinder, InDet::SCT_Clusterization, InDet::SiSPSeededTrackFinder, SCTRawDataProvider, InDet::GNNSeedingTrackMaker, SCT_PrepDataToxAOD, RoIBResultToxAOD, InDet::SiSPGNNTrackMaker, SCT_CablingCondAlgFromCoraCool, SCT_ReadCalibDataTestAlg, SCT_CablingCondAlgFromText, SCT_ReadCalibChipDataTestAlg, SCT_TestCablingAlg, SCT_ConfigurationConditionsTestAlg, ITkPixelCablingAlg, ITkStripCablingAlg, SCTEventFlagWriter, SCT_ConditionsSummaryTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_TdaqEnabledTestAlg, SCT_SiliconConditionsTestAlg, SCTSiLorentzAngleTestAlg, SCT_ByteStreamErrorsTestAlg, SCT_ConditionsParameterTestAlg, SCT_FlaggedConditionTestAlg, SCT_StripVetoTestAlg, SCT_RawDataToxAOD, and SCTSiPropertiesTestAlg.

Definition at line 44 of file AthReentrantAlgorithm.cxx.

 {
   // Reentrant algorithms are clonable.
   return true;
 }

◆ msg() [1/2]

MsgStream& AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                 {
     return this->msgStream();
   }

◆ msg() [2/2]

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

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                   {
     return this->msgStream(lvl);
   }

◆ msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::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< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

◆ renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                           {
     handlesArray.renounce();
   }

◆ setFilterPassed()

virtual void AthReentrantAlgorithm::setFilterPassed	(	bool	state,
		const EventContext &	ctx
	)		const

inlinevirtualinherited

Definition at line 139 of file AthReentrantAlgorithm.h.

                                                                             {
     execState( ctx ).setFilterPassed( state );
   }

◆ sysExecute()

StatusCode AthReentrantAlgorithm::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 67 of file AthReentrantAlgorithm.cxx.

 {
   return Gaudi::Algorithm::sysExecute (ctx);
 }

◆ sysInitialize()

StatusCode AthReentrantAlgorithm::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in InputMakerBase, and HypoBase.

Definition at line 96 of file AthReentrantAlgorithm.cxx.

                                                 {
   StatusCode sc=AthCommonDataStore<AthCommonMsg<Gaudi::Algorithm>>::sysInitialize();
  
   if (sc.isFailure()) {
     return sc;
   }
   
   ServiceHandle<ICondSvc> cs("CondSvc",name());
   for (auto h : outputHandles()) {
     if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
       // do this inside the loop so we don't create the CondSvc until needed
       if ( cs.retrieve().isFailure() ) {
         ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
         return StatusCode::SUCCESS;
       }      
       if (cs->regHandle(this,*h).isFailure()) {
         sc = StatusCode::FAILURE;
         ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                       << " with CondSvc");
       }
     }
   }
   return sc;  
 }

◆ sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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< Gaudi::Algorithm > >::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

◆ ATLAS_THREAD_SAFE [1/2]

std::mutex m_fillMapMutex LVL1::eFexTowerBuilder::ATLAS_THREAD_SAFE

mutableprivate

Definition at line 59 of file eFexTowerBuilder.h.

◆ ATLAS_THREAD_SAFE [2/2]

std::map<unsigned long long, std::pair<std::pair<int,int>,std::pair<int,int> > > m_scMap LVL1::eFexTowerBuilder::ATLAS_THREAD_SAFE

mutableprivate

Definition at line 60 of file eFexTowerBuilder.h.

◆ m_applyMasking

Gaudi::Property<bool> LVL1::eFexTowerBuilder::m_applyMasking {this,"ApplyMasking",true,"Apply masking of supercells based on provenance bits. Should be set to False for MC"}

private

Definition at line 76 of file eFexTowerBuilder.h.

◆ m_ddmKey

SG::ReadCondHandleKey<CaloSuperCellDetDescrManager> LVL1::eFexTowerBuilder::m_ddmKey {this,"CaloSuperCellDetDescrManager","CaloSuperCellDetDescrManager",""}

private

Definition at line 67 of file eFexTowerBuilder.h.

◆ m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

◆ m_eFEXSuperCellTowerIdProviderTool

ToolHandle<IeFEXSuperCellTowerIdProvider> LVL1::eFexTowerBuilder::m_eFEXSuperCellTowerIdProviderTool {this, "eFEXSuperCellTowerIdProviderTool", "LVL1::eFEXSuperCellTowerIdProvider", "Tool that provides tower-FOGA mapping"}

private

Definition at line 74 of file eFexTowerBuilder.h.

◆ m_eiKey

SG::ReadHandleKey<xAOD::EventInfo> LVL1::eFexTowerBuilder::m_eiKey {this,"EventInfoKey","EventInfo",""}

private

Definition at line 65 of file eFexTowerBuilder.h.

◆ m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

◆ m_extendedExtraObjects

DataObjIDColl AthReentrantAlgorithm::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 153 of file AthReentrantAlgorithm.h.

◆ m_mappingFile

Gaudi::Property<std::string> LVL1::eFexTowerBuilder::m_mappingFile {this, "MappingFile", "L1CaloFEXByteStream/2023-02-13/scToEfexTowers.root", "PathResolver location to mapping file"}

private

Definition at line 73 of file eFexTowerBuilder.h.

◆ m_outKey

SG::WriteHandleKey<xAOD::eFexTowerContainer> LVL1::eFexTowerBuilder::m_outKey {this, "eFexContainerWriteKey", "L1_eFexEmulatedTowers", "Name of the output container"}

private

Definition at line 71 of file eFexTowerBuilder.h.

◆ m_scellKey

SG::ReadHandleKey<CaloCellContainer> LVL1::eFexTowerBuilder::m_scellKey { this, "CaloCellContainerReadKey", "SCell", "Read handle key for the supercells"}

private

Definition at line 69 of file eFexTowerBuilder.h.

◆ m_ttKey

SG::ReadHandleKey<xAOD::TriggerTowerContainer> LVL1::eFexTowerBuilder::m_ttKey { this, "TriggerTowerContainerReadKey", "xAODTriggerTowers", "Read handle key for the triggerTowers"}

private

Definition at line 70 of file eFexTowerBuilder.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ eFexTowerBuilder()

◆ ~eFexTowerBuilder()

Member Function Documentation

◆ cardinality()

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

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ execute()

◆ extraDeps_update_handler()

◆ extraOutputDeps()

◆ fillMap()

◆ fillTowers()

◆ filterPassed()

◆ initialize()

◆ inputHandles()

◆ isClonable()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ setFilterPassed()

◆ sysExecute()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ ATLAS_THREAD_SAFE [1/2]

◆ ATLAS_THREAD_SAFE [2/2]

◆ m_applyMasking

◆ m_ddmKey

◆ m_detStore

◆ m_eFEXSuperCellTowerIdProviderTool

◆ m_eiKey

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_mappingFile

◆ m_outKey

◆ m_scellKey

◆ m_ttKey

◆ m_varHandleArraysDeclared

◆ m_vhka