SubtractedDiscSurface.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
*/

// SubtractedDiscSurface.h, (c) ATLAS Detector software
 
#ifndef TRKGEOMETRYSURFACES_SUBTRACTEDDISCSURFACE_H
#define TRKGEOMETRYSURFACES_SUBTRACTEDDISCSURFACE_H
 
// Trk
#include "TrkDetDescrUtils/AreaExcluder.h"
#include "TrkEventPrimitives/ParamDefs.h"
#include "TrkSurfaces/DiscSurface.h"
 
// Geometry & Maths
#include "GeoPrimitives/GeoPrimitives.h"
 
class MsgStream;
 
namespace Trk {

 
class SubtractedDiscSurface final : public DiscSurface
{
public:
  SubtractedDiscSurface() = default;
  SubtractedDiscSurface(SubtractedDiscSurface&&) = default;
  SubtractedDiscSurface& operator=(SubtractedDiscSurface&&) = default;
  virtual ~SubtractedDiscSurface() = default;

  SubtractedDiscSurface(const DiscSurface& ps, std::shared_ptr<const AreaExcluder> vol, bool shared);

  SubtractedDiscSurface(const SubtractedDiscSurface& psf);

  SubtractedDiscSurface(const SubtractedDiscSurface& psf,
                        const Amg::Transform3D& shift);

  SubtractedDiscSurface& operator=(const SubtractedDiscSurface& psf);

  virtual bool operator==(const Surface& sf) const override;

  bool shared() const;

  virtual bool insideBounds(const Amg::Vector2D& locpos,
                            double tol1 = 0.,
                            double tol2 = 0.) const override final;

  const AreaExcluder* subtractedVolume() const;

  virtual std::string name() const override final
  {
    return "Trk::SubtractedDiscSurface";
  }
 
protected:
  std::shared_ptr<const AreaExcluder> m_subtrVol{nullptr};
  bool m_shared{true};
};
 
inline bool
SubtractedDiscSurface::insideBounds(const Amg::Vector2D& locpos,
                                    double tol1,
                                    double tol2) const
{
  // no subtracted Volume exists
  if (!m_subtrVol.get()) {
    return (this->bounds().inside(locpos, tol1, tol2));
  }
  // subtracted Volume exists, needs to be checked
  double rPos = locpos[Trk::locR];
  double phiPos = locpos[Trk::locPhi];
  Amg::Vector3D gp(rPos * cos(phiPos), rPos * sin(phiPos), 0.);
  if (m_shared) {
    return (this->bounds().inside(locpos, tol1, tol2) &&
            m_subtrVol.get()->inside(gp, 0.));
  }
  bool in(this->bounds().inside(locpos, tol1, tol2) &&
          !m_subtrVol.get()->inside(gp, 0.));
 
  return in;
}
 
inline bool
SubtractedDiscSurface::shared() const
{
  return m_shared;
}
 
inline const AreaExcluder*
SubtractedDiscSurface::subtractedVolume() const
{
  return m_subtrVol.get();
}
 
} // end of namespace
 
#endif // TRKGEOMETRYSURFACES_SUBTRACTEDDISCSURFACE_H