Egamma1_OnlineMapNbhood.h

Go to the documentation of this file.

/*
 *   Copyright (C) 2002-2026 CERN for the benefit of the ATLAS collaboration
 */
 
#ifndef GLOBALSIM_EGAMMA1_ONLINEMAPNBHOOD_H
#define GLOBALSIM_EGAMMA1_ONLINEMAPNBHOOD_H
 
/*
  This Algorithm finds and outputs GlobalLArCells in the neighbourhood of eFEX
  RoIs. These neighhoods are used to run various Egamma1 Algorithms in GlobalSim.
  This method uses the Athena offline calo map to form the neighbourhood.
  It is seeded from the GlobalLArCells.
*/
 
#include "eFexRoIAlgTool.h"
#include "GlobalLArCell.h"
#include "GlobalLArCellContainer.h"
 
#include "../IO/LArStripNeighborhood.h"
#include "../IO/eEmNbhoodTOB.h"
 
#include "Identifier/Identifier.h"
 
#include "AthenaBaseComps/AthReentrantAlgorithm.h"
#include "xAODEventInfo/EventInfo.h"
 
#include <vector>
 
//forward declarations for calorimenter ID helpers
class CaloCell_ID;
class LArEM_ID;
 
namespace GlobalSim {
 
  class Egamma1_OnlineMapNbhood: public AthReentrantAlgorithm {
  public:
    
    Egamma1_OnlineMapNbhood(const std::string& name, ISvcLocator* pSvcLocator);

    virtual StatusCode  initialize() override;
    virtual StatusCode  execute(const EventContext& ) const override;    
 
  private:
 
    const CaloCell_ID*     m_calocell_id{};
    const LArEM_ID*        m_larem_id{};

    SG::ReadHandleKey<xAOD::EventInfo> m_eventInfoKey{
      this,
    "EventInfo",
    "EventInfo",
    "EventInfo name"};
    
    ToolHandle<eFexRoIAlgTool>
    m_roiAlgTool{
      this,
      "roiAlgTool",
      "GlobalSim::eFexRoIAlgTool",
      "AlgTool to provide a vector<const xAOD::eFexEMRoI*>"};
    
    SG::ReadHandleKey<GlobalSim::GlobalLArCellContainer>
    m_gblLArCellContainerKey {
      this,
      "GlobalLArCellsKey",
      "GlobalLArCells",
      "Key for the container of the LAr cells sent to Global"}; 

    SG::WriteHandleKey<std::vector<float>>
    m_eFEXetaKey {
      this,
      "eFEXetaKey",
      "eFEXeta"};

    SG::WriteHandleKey<std::vector<float>>
    m_eFEXphiKey {
      this,
      "eFEXphiKey",
      "eFEXphi"};

    SG::WriteHandleKey<std::vector<float>>
    m_FailedeFEXetaKey {
      this,
      "FailedeFEXetaKey",
      "FailedeFEXeta"};

    SG::WriteHandleKey<std::vector<float>>
    m_FailedeFEXphiKey {
      this,
      "FailedeFEXphiKey",
      "FailedeFEXphi"};

    SG::WriteHandleKey<IOBitwise::eEmNbhoodTOBContainer>
    m_neighKey {
      this,
      "stripNeighborhoodTOBKey",
      "stripNeighborhoodTOBContainer",
      "location to write strip neighborhoods of EFex RoIs, with the associated TOBs"};
 
    
    Gaudi::Property<bool> m_dump {
      this,
      "dump",
      false,
      "flag to enable dumps"};

    Gaudi::Property<bool> m_dumpTerse {
      this,
      "dumpTerse",
      false,
      "flag to enable terse dumps"};

    StatusCode
    findNeighborhoods_OnlineMap(const std::vector<const xAOD::eFexEMRoI*>&,
                const GlobalSim::GlobalLArCellContainer&,
                std::vector<bool>&,
                IOBitwise::eEmNbhoodTOBContainer&) const;

    StatusCode
    findNeighborhood_OnlineMap(const xAOD::eFexEMRoI*,
                   const GlobalSim::GlobalLArCellContainer&,
                   bool&,
                   IOBitwise::eEmNbhoodTOBContainer&) const;    

    bool
    findSeedCell(float,
         float,
         const GlobalSim::GlobalLArCellContainer&,
         Identifier&) const;

    bool
    findHalfStrips(int,
          float,
          float,
          const GlobalSim::GlobalLArCellContainer&,
          Identifier&) const;

    StatusCode
    findWindow(IdentifierHash,
           const GlobalSim::GlobalLArCellContainer&,
           std::vector<std::vector<std::shared_ptr<const GlobalSim::GlobalLArCell>> >&) const;    

    StatusCode
    findMaxima(IdentifierHash&,
           std::vector<std::vector<std::shared_ptr<const GlobalSim::GlobalLArCell>> >&) const;    

    std::shared_ptr<const GlobalLArCell> findMax(std::vector<std::shared_ptr<const GlobalLArCell>>&) const;
  };
 
}
#endif