GeoIDSvc.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef ISF_SERVICES_GEOIDSVC_H
#define ISF_SERVICES_GEOIDSVC_H 1
 
// framework includes
#include "GaudiKernel/ServiceHandle.h"
#include "AthenaBaseComps/AthService.h"
 
// STL includes
#include <vector>
#include <list>
#include <set>
 
// DetectorDescription
#include "AtlasDetDescr/AtlasRegion.h"
 
// ISF includes
#include "ISF_Interfaces/IGeoIDSvc.h"
 
// EnvelopeDefinitionService
#include "SubDetectorEnvelopes/IEnvelopeDefSvc.h" //Required for RZPair definition.
 
namespace ISF {
 
  typedef std::pair<double, AtlasDetDescr::AtlasRegion>          RadiusGeoIDPair;
  typedef std::list<RZPair>                                      RZPairList;
 
  struct SortByRadius {
    bool operator() (const RadiusGeoIDPair& lhs, const RadiusGeoIDPair& rhs) const
    {return lhs.first<rhs.first;}
  };
  typedef std::set<RadiusGeoIDPair, SortByRadius> RadiusGeoIDPairSet;

  class GeoIDSvc : public extends<AthService, ISF::IGeoIDSvc> {
  public:
    GeoIDSvc(const std::string& name,ISvcLocator* svc);

    ~GeoIDSvc() = default;
 
    // Athena algtool's Hooks
    StatusCode  initialize();

    AtlasDetDescr::AtlasRegion    identifyGeoID(const Amg::Vector3D &pos) const;

    ISF::InsideType inside(const Amg::Vector3D &pos, AtlasDetDescr::AtlasRegion geoID) const;
 

    AtlasDetDescr::AtlasRegion identifyNextGeoID(const Amg::Vector3D &pos, const Amg::Vector3D &dir) const;
 
  private:
    std::unique_ptr<RZPairList> prepareRZPairs( AtlasDetDescr::AtlasRegion geoID);

    bool       checkSymmetric( RZPairList& positiveZ, RZPairList& negativeZ);
 
    template <typename T> inline int sign(T val);

    ServiceHandle<IEnvelopeDefSvc> m_envDefSvc{this, "EnvelopeDefSvc", "ISF_ISFEnvelopeDefSvc",
        "The EnvelopeDefinitionService describing the AtlasRegion boundaries."};

    DoubleProperty m_tolerance{this,"Tolerance", 1e-5,
        "Estimated tolerance within which coordinates are considered being equal."};

    double                *m_zBins{};
    int                    m_numZBins{0};
    RadiusGeoIDPair       *m_radiusBins{};
    int                    m_maxRBins{AtlasDetDescr::fNumAtlasRegions+1};
  };
 
  //
  // inline functions
  //

  template <typename T> inline int ISF::GeoIDSvc::sign(T val) {
    return (T(0) < val) - (val < T(0));
  }
 
} // ISF namespace
 
#endif //> !ISF_SERVICES_GEOIDSVC_H