SimplePolygonBrepVolumeBounds.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// SimplePolygonBrepVolumeBounds.h, (c) ATLAS Detector software
 
#ifndef TRKVOLUMES_SIMPLEPOLYGONBREPVOLUMEBOUNDS_H
#define TRKVOLUMES_SIMPLEPOLYGONBREPVOLUMEBOUNDS_H
 
// Trk
#include "GeoPrimitives/GeoPrimitives.h"
#include "TrkDetDescrUtils/ObjectsAccessor.h"
#include "TrkVolumes/VolumeBounds.h"
 
class MsgStream;
 
namespace Trk {
 
class Surface;
class PlaneSurface;
class RectangleBounds;
 
class SimplePolygonBrepVolumeBounds final: public VolumeBounds {
 public:
  SimplePolygonBrepVolumeBounds();
 
  SimplePolygonBrepVolumeBounds(std::vector<std::pair<float, float> > xyvtx,
                                float hlengthz);
 
  SimplePolygonBrepVolumeBounds(std::vector<std::pair<double, double> > xyvtx,
                                double hlengthz);
 
  SimplePolygonBrepVolumeBounds(const SimplePolygonBrepVolumeBounds& bobo);
 
  virtual ~SimplePolygonBrepVolumeBounds();
 
  SimplePolygonBrepVolumeBounds& operator=(
      const SimplePolygonBrepVolumeBounds& bobo);
 
  SimplePolygonBrepVolumeBounds* clone() const override final;
 
  bool inside(const Amg::Vector3D&, double tol = 0.) const override final; 
 
  const std::vector<const Trk::Surface*>* decomposeToSurfaces
  (const Amg::Transform3D& transform) override final;
 
  ObjectAccessor boundarySurfaceAccessor(
      const Amg::Vector3D& gp, const Amg::Vector3D& dir,
      bool forceInside = false) const override final;
  std::vector<std::pair<double, double> > xyVertices() const;
  double halflengthZ() const;
 
  const Trk::Volume* combinedVolume() const;
 
  const Trk::Volume* envelope() const;
 
  MsgStream& dump(MsgStream& sl) const override;
 
  std::ostream& dump(std::ostream& sl) const override;
 
 private:
  void processSubVols();
  Trk::PlaneSurface* sideSurf(const Amg::Transform3D&, unsigned int,
                              unsigned int) const;
  static bool Xor(bool x, bool y) ;
 
  bool Left(std::pair<double, double> a, std::pair<double, double> b,
            std::pair<double, double> c) const;
 
  bool Intersect(std::pair<double, double> a, std::pair<double, double> b,
                 std::pair<double, double> c,
                 std::pair<double, double> d) const;
 
  bool InCone(int i, int j,
              const std::vector<std::pair<double, double> >& inputVertices) const;
 
  bool Diagonalie(int i, int j,
                  const std::vector<std::pair<double, double> >& inputVertices) const;
 
  bool Diagonal(
      int i, int j,
      const std::vector<std::pair<double, double> >& inputVertices) const;
 
  std::vector<std::pair<double, double> > TriangulatePolygon(
      const std::vector<std::pair<double, double> >& Vertices) const;
 
  std::vector<std::pair<double, double> > TriangulatePolygonCheck(
      const std::vector<std::pair<double, double> >& Vertices);
 
  std::vector<std::pair<double, double> > m_xyVtx;  
  double m_halfX;  
  double m_halfY;  
  double m_halfZ;  
 
  int m_ordering;  
  const Trk::Volume* m_combinedVolume;  
  const Trk::Volume* m_envelope;        
 
  Trk::EightObjectsAccessor m_objectAccessor;  
};
 
inline SimplePolygonBrepVolumeBounds* SimplePolygonBrepVolumeBounds::clone()
    const {
  return new SimplePolygonBrepVolumeBounds(*this);
}
inline std::vector<std::pair<double, double> >
SimplePolygonBrepVolumeBounds::xyVertices() const {
  return m_xyVtx;
}
 
inline double SimplePolygonBrepVolumeBounds::halflengthZ() const {
  return m_halfZ;
}
 
inline const Trk::Volume* SimplePolygonBrepVolumeBounds::combinedVolume()
    const {
  return m_combinedVolume;
}
 
inline const Trk::Volume* SimplePolygonBrepVolumeBounds::envelope() const {
  return m_envelope;
}
 
inline ObjectAccessor SimplePolygonBrepVolumeBounds::boundarySurfaceAccessor(
    const Amg::Vector3D&, const Amg::Vector3D&, bool) const {
  return Trk::ObjectAccessor(m_objectAccessor);
}
 
}  // namespace Trk
 
#endif  // TRKVOLUMES_SIMPLEPOLYGONBREPVOLUMEBOUNDS_H