CylinderVolumeBounds.h

Go to the documentation of this file.

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

// CylinderVolumeBounds.h, (c) ATLAS Detector software
 
#ifndef TRKVOLUMES_CYLINDERVOLUMESBOUNDS_H
#define TRKVOLUMES_CYLINDERVOLUMESBOUNDS_H
 
// Trk
#include "GeoPrimitives/GeoPrimitives.h"
#include "TrkVolumes/CylinderVolumeBoundaryAccessors.h"
#include "TrkVolumes/VolumeBounds.h"
 
class MsgStream;
 
namespace Trk {
 
class Surface;
class DiscBounds;
class CylinderBounds;
class RectangleBounds;

 
class CylinderVolumeBounds final: public VolumeBounds {
 public:
  CylinderVolumeBounds();

  CylinderVolumeBounds(double radius, double halez);

  CylinderVolumeBounds(double rinner, double router, double halez);

  CylinderVolumeBounds(double rinner, double router, double halfPhiSector,
                       double halez);

  CylinderVolumeBounds(const CylinderVolumeBounds& cylbo);

  virtual ~CylinderVolumeBounds();

  CylinderVolumeBounds& operator=(const CylinderVolumeBounds& cylbo);

  CylinderVolumeBounds* clone() const override;

  bool inside(const Amg::Vector3D&, double tol = 0.) const override final;

  virtual std::vector<std::unique_ptr<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;

  double innerRadius() const;

  double outerRadius() const;

  double mediumRadius() const;

  double deltaRadius() const;

  double halfPhiSector() const;

  double halflengthZ() const;

  MsgStream& dump(MsgStream& sl) const override final;
  std::ostream& dump(std::ostream& sl) const override final;
 
 private:
  std::shared_ptr<CylinderBounds> innerCylinderBounds() const;

  std::shared_ptr<CylinderBounds> outerCylinderBounds() const;

  std::shared_ptr<DiscBounds> bottomDiscBounds() const;

  std::shared_ptr<DiscBounds> topDiscBounds() const;

  std::shared_ptr<RectangleBounds> sectorPlaneBounds() const;

  void createBoundarySurfaceAccessors();
 
  double m_innerRadius;
  double m_outerRadius;
  double m_halfPhiSector;
  double m_halfZ;

  CylinderVolumeBoundaryAccessors m_boundaryAccessors;

  static const double s_numericalStable;
};
 
inline CylinderVolumeBounds* CylinderVolumeBounds::clone() const {
  return new CylinderVolumeBounds(*this);
}
 
inline bool CylinderVolumeBounds::inside(const Amg::Vector3D& pos,
                                         double tol) const {
  double ros = pos.perp();
  // bool insidePhi =  fabs(pos.phi()) <= m_halfPhiSector + tol;
  bool insidePhi = cos(pos.phi()) >= cos(m_halfPhiSector) - tol;
  bool insideR =
      insidePhi ? ((ros >= m_innerRadius - tol) && (ros <= m_outerRadius + tol))
                : false;
  bool insideZ = insideR ? (fabs(pos.z()) <= m_halfZ + tol) : false;
  return (insideZ && insideR && insidePhi);
}
 
inline double CylinderVolumeBounds::innerRadius() const {
  return m_innerRadius;
}
 
inline double CylinderVolumeBounds::outerRadius() const {
  return m_outerRadius;
}
 
inline double CylinderVolumeBounds::mediumRadius() const {
  return 0.5 * (m_innerRadius + m_outerRadius);
}
 
inline double CylinderVolumeBounds::deltaRadius() const {
  return (m_outerRadius - m_innerRadius);
}
 
inline double CylinderVolumeBounds::halfPhiSector() const {
  return m_halfPhiSector;
}
 
inline double CylinderVolumeBounds::halflengthZ() const { return m_halfZ; }
 
inline std::shared_ptr<DiscBounds> CylinderVolumeBounds::topDiscBounds() const {
  return this->bottomDiscBounds();
}
 
}  // namespace Trk
 
#endif  // TRKVOLUMES_CYLINDERVOLUMESBOUNDS_H