SurfaceBounds.h

Go to the documentation of this file.

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

// SurfaceBounds.h, (c) ATLAS Detector software
 
#ifndef TRKSURFACES_SURFACEBOUNDS_H
#define TRKSURFACES_SURFACEBOUNDS_H
 
// STD
#include <iomanip>
#include <iostream>
 
// GaudiKernel
#include "GaudiKernel/MsgStream.h"
 
// Geo & Math library
#include "GeoPrimitives/GeoPrimitives.h"
 
// Trk included
#include "TrkSurfaces/BoundaryCheck.h"
 
#ifdef TRKDETDESCR_USEFLOATPRECISON
typedef float TDD_real_t;
#else
typedef double TDD_real_t;
#endif
 
namespace Trk {

 
class SurfaceBounds
{
 
public:
  enum BoundsType
  {
    Cone = 0,
    Cylinder = 1,
    Diamond = 2,
    Disc = 3,
    Ellipse = 5,
    Rectangle = 6,
    RotatedTrapezoid = 7,
    Trapezoid = 8,
    Triangle = 9,
    DiscTrapezoidal = 10,
    Annulus = 11,
    Other = 12
 
  };

  SurfaceBounds() = default;
 
  SurfaceBounds(const SurfaceBounds&) = default;
  SurfaceBounds& operator=(const SurfaceBounds&) = default;
  SurfaceBounds(SurfaceBounds&&) noexcept = default;
  SurfaceBounds& operator=(SurfaceBounds&&) noexcept = default;
  virtual ~SurfaceBounds() = default;

  virtual SurfaceBounds* clone() const = 0;

  virtual bool operator==(const SurfaceBounds& sb) const = 0;

  virtual bool operator!=(const SurfaceBounds& sb) const;

  virtual BoundsType type() const = 0;

  virtual bool inside(const Amg::Vector2D& locpo, double tol1 = 0., double tol2 = 0.) const = 0;
  virtual bool inside(const Amg::Vector2D& locpo, const BoundaryCheck& bchk) const = 0;
  virtual bool insideLoc1(const Amg::Vector2D& locpo, double tol1 = 0.) const = 0;

  virtual bool insideLoc2(const Amg::Vector2D& locpo, double tol2 = 0.) const = 0;

  virtual double minDistance(const Amg::Vector2D& pos) const = 0;

  virtual double r() const = 0;

  virtual MsgStream& dump(MsgStream& sl) const = 0;

  virtual std::ostream& dump(std::ostream& sl) const = 0;
 
protected:
  void swap(double& b1, double& b2);

  virtual void initCache() {}
};
 
inline void
SurfaceBounds::swap(double& b1, double& b2)
{
  double tmp = b1;
  b1 = b2;
  b2 = tmp;
}
 
inline bool
SurfaceBounds::operator!=(const SurfaceBounds& sb) const
{
  return !((*this) == sb);
}

MsgStream&
operator<<(MsgStream& sl, const SurfaceBounds& sb);
std::ostream&
operator<<(std::ostream& sl, const SurfaceBounds& sb);
 
} // end of namespace
 
#endif // TRKSURFACES_SURFACEBOUNDS_H