#include <JTowerMappingDataCondAlg.h>

Inheritance diagram for LVL1::JTowerMappingDataCondAlg:

Collaboration diagram for LVL1::JTowerMappingDataCondAlg:

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

virtual	~JTowerMappingDataCondAlg () override

virtual StatusCode	initialize () override

virtual StatusCode	execute (const EventContext &ctx) const final

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

Protected Attributes
const CaloCell_SuperCell_ID *	m_scid {nullptr}

const CaloCell_ID *	m_ccid {nullptr}

Private Types
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
virtual const JGTowerBase_ID *	getJGTowerID () const override

virtual StatusCode	buildForwardMapping (JGTowerMappingData &data, const CaloSuperCellDetDescrManager *mgr) const override

virtual StatusCode	loadTowerAreas (JGTowerMappingData &data) const override

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
Gaudi::Property< std::string >	m_towerAreasFile

const JTower_ID *	m_jtowerID {nullptr}

std::vector< float >	m_towerAreas

SG::WriteCondHandleKey< JGTowerMappingData >	m_outputKey

SG::ReadCondHandleKey< CaloSuperCellDetDescrManager >	m_caloSuperCellMgrKey

Gaudi::Property< bool >	m_mapTileCells {this, "MapTileCells", false}

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 13 of file JTowerMappingDataCondAlg.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

◆ JTowerMappingDataCondAlg()

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

Definition at line 14 of file JTowerMappingDataCondAlg.cxx.

       : JGTowerMappingDataCondAlgBase(name, pSvcLocator)
   {
   }

◆ ~JTowerMappingDataCondAlg()

LVL1::JTowerMappingDataCondAlg::~JTowerMappingDataCondAlg ( )

overridevirtual

Definition at line 19 of file JTowerMappingDataCondAlg.cxx.

19 {}

Member Function Documentation

◆ buildForwardMapping()

StatusCode LVL1::JTowerMappingDataCondAlg::buildForwardMapping	(	JGTowerMappingData &	data,
		const CaloSuperCellDetDescrManager *	mgr
	)		const

overrideprivatevirtual

Implements LVL1::JGTowerMappingDataCondAlgBase.

Definition at line 49 of file JTowerMappingDataCondAlg.cxx.

   {
     ATH_MSG_INFO("# central JTowers " << data.size());
     // Create a single JTower per forward supercell
     for (std::size_t scHashIdx = 0; scHashIdx < m_scid->calo_cell_hash_max(); ++scHashIdx)
     {
       Identifier scID = m_scid->cell_id(scHashIdx);
       if (m_scid->is_tile(scID) && m_scid->sampling(scID) != 2)
         continue; // skipping tile SCs for now
       const CaloDetDescrElement *dde = scDetMgr->get_element(scID);
       if (!dde)
       {
         ATH_MSG_ERROR("Error loading CaloDetDescrElement");
         return StatusCode::FAILURE;
       }
       float superCellEta = dde->eta_raw();
       float superCellPhi = dde->phi_raw();
       if (std::abs(superCellEta) < 3.1)
         continue;
       float superCellDEta = dde->deta();
       float superCellDPhi = dde->dphi();
       JGTowerHelper helper(superCellEta, superCellDEta, superCellPhi, superCellDPhi);
       helper.SetSCIndices(scHashIdx);
  
       std::string detName = CaloSampling::getSamplingName(m_ccid->calo_sample(scID));
       if (detName.find("FCAL0") != std::string::npos)
         helper.SetSampling(2);
       else if (detName.find("FCAL1") != std::string::npos)
         helper.SetSampling(3);
       else if (detName.find("FCAL2") != std::string::npos)
         helper.SetSampling(4);
       else
         helper.SetSampling(0);
       data.push_back(std::move(helper));
     }
     ATH_MSG_INFO("# towers (including forward)");
     return StatusCode::SUCCESS;
   }

◆ cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality

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 75 of file AthCommonReentrantAlgorithm.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::JGTowerMappingDataCondAlgBase::execute ( const EventContext & ctx ) const

finalvirtualinherited

Definition at line 34 of file JGTowerMappingDataCondAlgBase.cxx.

   {
     SG::WriteCondHandle<JGTowerMappingData> outputHandle(m_outputKey, ctx);
     if (outputHandle.isValid())
     {
       // In serial athena conditions algorithms are run on every event.
       // Therefore return early if the output object already exists.
       return StatusCode::SUCCESS;
     }
     const JGTowerBase_ID *jgTowerID = getJGTowerID();
     SG::ReadCondHandle<CaloSuperCellDetDescrManager> caloSuperCellMgrHandle{m_caloSuperCellMgrKey,ctx};
     const CaloSuperCellDetDescrManager *scDetMgr = *caloSuperCellMgrHandle;
     // JGTowerMappingData is typedef of std::vector<JGTowerHelper>
     auto data = std::make_unique<JGTowerMappingData>();
  
     // Store the towers in a map by their location to identify ones which are in front/behind each other
     // tuple elements: pos_neg, region, eta, phi
     std::map<std::tuple<int, int, int, int>, std::vector<std::size_t>> towerGrid;
     // For speed later, also store the max phi index and min/max eta indices per (int) region code
     std::map<int, std::tuple<int, int, int>> regionParams;
  
     for (std::size_t towerIdx = 0; towerIdx < jgTowerID->tower_hash_max(); ++towerIdx)
     {
       Identifier towerID = jgTowerID->tower_id(towerIdx);
       Identifier regionID = jgTowerID->region_id(towerID);
  
       int detSide = jgTowerID->pos_neg(towerID) < 0 ? -1 : +1;
  
       // Get the phi and eta boundaries of the towers (only advisable for barrel)
       float jDEta = jgTowerID->etaGranularity(regionID);
       float jDPhi = jgTowerID->phiGranularity(regionID);
       int nTowers = (int)std::ceil(TMath::Pi() / jDPhi);
       jDPhi = TMath::Pi() / nTowers;
  
       float jEta = detSide * ((jgTowerID->eta(towerID) + 0.5) * jDEta + jgTowerID->eta0(regionID));
       float jPhi = TVector2::Phi_mpi_pi((jgTowerID->phi(towerID) + 0.5) * jDPhi + jgTowerID->phi0(regionID));
  
       data->emplace_back(jEta, jDEta, jPhi, jDPhi);
       JGTowerHelper &helper = data->back();
       helper.SetSampling(jgTowerID->sampling(towerID));
  
       std::vector<std::size_t> &gridLoc = towerGrid[std::make_tuple(
           jgTowerID->pos_neg(towerID),
           jgTowerID->region(towerID),
           jgTowerID->eta(towerID),
           jgTowerID->phi(towerID))];
       regionParams[jgTowerID->region(towerID)] = std::make_tuple(
           jgTowerID->phi_max(regionID),
           jgTowerID->eta_min(regionID),
           jgTowerID->eta_max(regionID));
  
       // Now insert our index into this such that the vector remains sorted in ascending order of sampling
       gridLoc.insert(std::lower_bound(
                          gridLoc.begin(), gridLoc.end(), towerIdx,
                          [&data](std::size_t lhs, std::size_t rhs) { return data->at(lhs).sampling() < data->at(rhs).sampling(); }),
                      towerIdx);
  
       for (std::size_t scIdx = 0; scIdx < m_scid->calo_cell_hash_max(); ++scIdx)
       {
         Identifier scID = m_scid->cell_id(scIdx);
         const CaloDetDescrElement *dde = scDetMgr->get_element(scID);
         if (!dde)
         {
           ATH_MSG_ERROR("Error loading CaloDetDescrElement");
           return StatusCode::FAILURE;
         }
         if (m_scid->is_tile(scID) && m_scid->sampling(scID) != 2)
         {
           if (!m_mapTileCells)
             continue;
           if (jgTowerID->sampling(towerID) != 0 && helper.inBox(dde->eta_raw(), dde->phi_raw()))
             helper.SetTileIndices(scIdx);
         }
         else
         {
           float superCellEta = dde->eta_raw();
           // Treat forward super cells separately
           if (std::abs(superCellEta) > 3.2)
             continue;
           // Very special case being handle here. At the end of the barrel there is a constant-eta ring that has
           // eta_raw=1.4 (which is the midpoint of the scell). These will be put into the g/jTower that starts at eta=1.4
           if (std::abs(std::abs(superCellEta) - 1.4) < 0.001 &&
               m_scid->region(scID) == 0 &&
               m_scid->sampling(scID) == 2)
             superCellEta += (superCellEta > 0 ? 0.05 : -0.05);
           if (helper.inBox(superCellEta, dde->phi_raw()) &&
               jgTowerID->pos_neg(towerID) * m_scid->pos_neg(scID) > 0 &&
               ((m_scid->is_em(scID) && jgTowerID->sampling(towerID) == 0) ||
                (!m_scid->is_em(scID) && jgTowerID->sampling(towerID) == 1)))
             helper.SetSCIndices(scIdx);
         }
       } //> end loop over SC indices
     }   //> end loop over tower IDs
  
     for (auto &gridLocPair : towerGrid)
     {
       std::vector<std::size_t> &towerIndices = gridLocPair.second;
       // Figure out the relations to this square in the grid
       int pos_neg = std::get<0>(gridLocPair.first);
       int region = std::get<1>(gridLocPair.first);
       int iEta = std::get<2>(gridLocPair.first);
       int iPhi = std::get<3>(gridLocPair.first);
       int phiMax = std::get<0>(regionParams[region]);
       int etaMin = std::get<1>(regionParams[region]);
       int etaMax = std::get<2>(regionParams[region]);
       auto nextPhiItr = towerGrid.find(std::make_tuple(pos_neg, region, iEta, iPhi == phiMax ? 0 : iPhi + 1));
       auto prevPhiItr = towerGrid.find(std::make_tuple(pos_neg, region, iEta, iPhi == 0 ? phiMax : iPhi - 1));
       auto nextEtaItr = towerGrid.end();
       auto prevEtaItr = towerGrid.end();
       // We want 'next' and 'prev' to refer to increasing eta, not increasing |eta|
       // Increasing 'region' increases |eta| so in what follows we correct for this by
       // choosing whether we set nextEtaItr or prevEtaItr based on the sign of pos_neg
       if (iEta == etaMax)
       {
         // Check if the next region a) exists and b) has the right phi granularity
         auto itr = regionParams.find(region + 1);
         if (itr != regionParams.end() && std::get<0>(itr->second) == phiMax)
           // Set the next eta index to the next region and the eta_min index of that region (probably 0)
           (pos_neg > 0 ? nextEtaItr : prevEtaItr) = towerGrid.find(std::make_tuple(pos_neg, region + 1, std::get<1>(itr->second), iPhi));
       }
       else
         (pos_neg > 0 ? nextEtaItr : prevEtaItr) = towerGrid.find(std::make_tuple(pos_neg, region, iEta + 1, iPhi));
       if (iEta == etaMin)
       {
         if (region == 0)
           (pos_neg > 0 ? prevEtaItr : nextEtaItr) = towerGrid.find(std::make_tuple(pos_neg * -1, region, iEta, iPhi));
         else
         {
           auto itr = regionParams.find(region - 1);
           if (itr != regionParams.end() && std::get<0>(itr->second) == phiMax)
             // Set the previous eta index to the previous region and the eta_max index of that region
             (pos_neg > 0 ? prevEtaItr : nextEtaItr) = towerGrid.find(std::make_tuple(pos_neg, region - 1, std::get<2>(itr->second), iPhi));
         }
       }
       else
         (pos_neg > 0 ? prevEtaItr : nextEtaItr) = towerGrid.find(std::make_tuple(pos_neg, region, iEta - 1, iPhi));
  
       std::vector<std::size_t> nextEtaTowers;
       std::vector<std::size_t> prevEtaTowers;
       std::vector<std::size_t> nextPhiTowers;
       std::vector<std::size_t> prevPhiTowers;
       if (nextEtaItr != towerGrid.end())
         nextEtaTowers = nextEtaItr->second;
       if (prevEtaItr != towerGrid.end())
         prevEtaTowers = prevEtaItr->second;
       if (nextPhiItr != towerGrid.end())
         nextPhiTowers = nextPhiItr->second;
       if (prevPhiItr != towerGrid.end())
         prevPhiTowers = prevPhiItr->second;
       nextEtaTowers.resize(towerIndices.size(), SIZE_MAX);
       prevEtaTowers.resize(towerIndices.size(), SIZE_MAX);
       nextPhiTowers.resize(towerIndices.size(), SIZE_MAX);
       prevPhiTowers.resize(towerIndices.size(), SIZE_MAX);
  
       std::size_t previousIdx = SIZE_MAX;
       for (std::size_t iIdx = 0; iIdx < towerIndices.size(); ++iIdx)
       {
         // NB: iIdx here is an index of _indices_, i.e. the position in the towerIndex
         // vector as well as the next/prev eta/phi towers vectors
         std::size_t currentIdx = towerIndices.at(iIdx);
         JGTowerHelper &helper = data->at(currentIdx);
  
         helper.setNextEtaIndex(nextEtaTowers.at(iIdx));
         helper.setPreviousEtaIndex(prevEtaTowers.at(iIdx));
         helper.setNextPhiIndex(nextPhiTowers.at(iIdx));
         helper.setPreviousPhiIndex(prevPhiTowers.at(iIdx));
  
         helper.setIndexInFront(previousIdx);
         if (previousIdx != SIZE_MAX)
           data->at(previousIdx).setIndexBehind(currentIdx);
         previousIdx = currentIdx;
       }
     }
  
     ATH_CHECK(buildForwardMapping(*data, scDetMgr));
     ATH_CHECK(loadTowerAreas(*data));
  
     // Set this object to be valid for every event
     outputHandle.addDependency(IOVInfiniteRange::infiniteRunLB());
     ATH_CHECK(outputHandle.record(std::move(data)));
  
     return StatusCode::SUCCESS;
   }

◆ 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 & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps

overridevirtualinherited

Return the list of extra output dependencies.

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

Definition at line 94 of file AthCommonReentrantAlgorithm.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 BaseAlg::extraOutputDeps();
 }

◆ filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

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

◆ getJGTowerID()

virtual const JGTowerBase_ID* LVL1::JTowerMappingDataCondAlg::getJGTowerID ( ) const

inlineoverrideprivatevirtual

Implements LVL1::JGTowerMappingDataCondAlgBase.

Definition at line 25 of file JTowerMappingDataCondAlg.h.

25 { return m_jtowerID; }

◆ initialize()

StatusCode LVL1::JTowerMappingDataCondAlg::initialize ( )

overridevirtual

Reimplemented from LVL1::JGTowerMappingDataCondAlgBase.

Definition at line 21 of file JTowerMappingDataCondAlg.cxx.

   {
     ATH_CHECK(JGTowerMappingDataCondAlgBase::initialize());
     ATH_CHECK(detStore()->retrieve(m_jtowerID));
  
     std::unique_ptr<TFile> fIn(TFile::Open(
         PathResolver::find_calib_file(m_towerAreasFile.value()).c_str()));
     if (!fIn || fIn->IsZombie())
     {
       ATH_MSG_ERROR("Failed to open tower areas file " << m_towerAreasFile);
       return StatusCode::FAILURE;
     }
     TH1F *areaHist = (TH1F *)fIn->Get("jTowerArea_final_hist");
     if (!areaHist)
     {
       ATH_MSG_ERROR("Failed to read area histogram!");
       return StatusCode::FAILURE;
     }
     m_towerAreas.clear();
     m_towerAreas.reserve(areaHist->GetNbinsX());
     // Normalise the areas so that the first tower has an area of 1
     float scale = 1.0 / areaHist->GetBinContent(1);
     for (int towerIdx = 0; towerIdx < areaHist->GetNbinsX(); ++towerIdx)
       m_towerAreas.push_back(areaHist->GetBinContent(towerIdx + 1) * scale);
  
     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 AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

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

◆ loadTowerAreas()

StatusCode LVL1::JTowerMappingDataCondAlg::loadTowerAreas ( JGTowerMappingData & data ) const

overrideprivatevirtual

Implements LVL1::JGTowerMappingDataCondAlgBase.

Definition at line 90 of file JTowerMappingDataCondAlg.cxx.

   {
     if (data.size() != m_towerAreas.size())
     {
       ATH_MSG_ERROR("Number of towers (" << data.size() << ") and number of tower areas (" << m_towerAreas.size() << ") do not match!");
       return StatusCode::FAILURE;
     }
  
     for (std::size_t idx = 0; idx < m_towerAreas.size(); ++idx)
       data.at(idx).setArea(m_towerAreas.at(idx));
     return StatusCode::SUCCESS;
   }

◆ 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 AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed	(	bool	state,
		const EventContext &	ctx
	)		const

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

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

◆ sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::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 85 of file AthCommonReentrantAlgorithm.cxx.

 {
   return BaseAlg::sysExecute (ctx);
 }

◆ sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::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 61 of file AthCommonReentrantAlgorithm.cxx.

                                                                {
   StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::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

◆ m_caloSuperCellMgrKey

SG::ReadCondHandleKey<CaloSuperCellDetDescrManager> LVL1::JGTowerMappingDataCondAlgBase::m_caloSuperCellMgrKey

privateinherited

Initial value:

{

this,"CaloSuperCellDetDescrManager","CaloSuperCellDetDescrManager","SG key of the resulting CaloSuperCellDetDescrManager"}

Definition at line 31 of file JGTowerMappingDataCondAlgBase.h.

◆ m_ccid

const CaloCell_ID* LVL1::JGTowerMappingDataCondAlgBase::m_ccid {nullptr}

protectedinherited

Definition at line 42 of file JGTowerMappingDataCondAlgBase.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_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 AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

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

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

◆ m_jtowerID

const JTower_ID* LVL1::JTowerMappingDataCondAlg::m_jtowerID {nullptr}

private

Definition at line 24 of file JTowerMappingDataCondAlg.h.

◆ m_mapTileCells

Gaudi::Property<bool> LVL1::JGTowerMappingDataCondAlgBase::m_mapTileCells {this, "MapTileCells", false}

privateinherited

Definition at line 33 of file JGTowerMappingDataCondAlgBase.h.

◆ m_outputKey

SG::WriteCondHandleKey<JGTowerMappingData> LVL1::JGTowerMappingDataCondAlgBase::m_outputKey

privateinherited

Initial value:

{

this, "MappingData", ""}

Definition at line 29 of file JGTowerMappingDataCondAlgBase.h.

◆ m_scid

const CaloCell_SuperCell_ID* LVL1::JGTowerMappingDataCondAlgBase::m_scid {nullptr}

protectedinherited

Definition at line 41 of file JGTowerMappingDataCondAlgBase.h.

◆ m_towerAreas

std::vector<float> LVL1::JTowerMappingDataCondAlg::m_towerAreas

private

Definition at line 30 of file JTowerMappingDataCondAlg.h.

◆ m_towerAreasFile

Gaudi::Property<std::string> LVL1::JTowerMappingDataCondAlg::m_towerAreasFile

private

Initial value:

{

this, "TowerAreasFile", "Run3L1CaloSimulation/Noise/jTowerCorrection.20210209.r11881.root"}

Definition at line 22 of file JTowerMappingDataCondAlg.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

Protected Attributes

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ JTowerMappingDataCondAlg()

◆ ~JTowerMappingDataCondAlg()

Member Function Documentation

◆ buildForwardMapping()

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

◆ filterPassed()

◆ getJGTowerID()

◆ initialize()

◆ inputHandles()

◆ isClonable()

◆ loadTowerAreas()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ renounce()

◆ renounceArray()

◆ setFilterPassed()

◆ sysExecute()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_caloSuperCellMgrKey

◆ m_ccid

◆ m_detStore

◆ m_evtStore

◆ m_extendedExtraObjects

◆ m_jtowerID

◆ m_mapTileCells

◆ m_outputKey

◆ m_scid

◆ m_towerAreas

◆ m_towerAreasFile

◆ m_varHandleArraysDeclared

◆ m_vhka