Trk::TriangleBounds Class Reference

#include <TriangleBounds.h>

Inheritance diagram for Trk::TriangleBounds:

Collaboration diagram for Trk::TriangleBounds:

Public Types
enum	BoundValues {   bv_x1 = 0, bv_y1 = 1, bv_x2 = 2, bv_y2 = 3,   bv_x3 = 4, bv_y3 = 5, bv_length = 6 }
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 }

Public Member Functions
	TriangleBounds ()
	Default Constructor - needed for persistency. More...
	TriangleBounds (const TriangleBounds &tribo)=default
	Copy constructor. More...
TriangleBounds &	operator= (const TriangleBounds &recbo)=default
	Assignment Operator. More...
	TriangleBounds (TriangleBounds &&tribo) noexcept=default
	Move constructor. More...
TriangleBounds &	operator= (TriangleBounds &&recbo) noexcept=default
	Move Assignment Operator. More...
virtual	~TriangleBounds ()=default
	Destructor. More...
	TriangleBounds (const std::vector< std::pair< float, float >> &)
	Constructor with coordinates of vertices - floats. More...
	TriangleBounds (const std::vector< std::pair< double, double >> &)
	Constructor with coordinates of vertices - double. More...
	TriangleBounds (const Amg::Vector2D &p1, const Amg::Vector2D &p2, const Amg::Vector2D &p3)
	Constructor from three 2 Vectors. More...
virtual bool	operator== (const SurfaceBounds &sbo) const override
	Equality operator. More...
virtual TriangleBounds *	clone () const override
	Virtual constructor. More...
virtual BoundsType	type () const override final
	Return the type of the bounds for persistency. More...
virtual bool	inside (const Amg::Vector2D &locpo, double tol1=0., double tol2=0.) const override final
	This method checks if the provided local coordinates are inside the surface bounds. More...
virtual bool	inside (const Amg::Vector2D &locpo, const BoundaryCheck &bchk) const override final
virtual bool	insideLoc1 (const Amg::Vector2D &locpo, double tol1=0.) const override final
	This method checks inside bounds in loc1. More...
virtual bool	insideLoc2 (const Amg::Vector2D &locpo, double tol2=0.) const override final
	This method checks inside bounds in loc2. More...
virtual double	minDistance (const Amg::Vector2D &pos) const override final
	Minimal distance to boundary ( > 0 if outside and <=0 if inside) More...
std::vector< std::pair< TDD_real_t, TDD_real_t > >	vertices () const
	This method returns the coordinates of vertices. More...
virtual double	r () const override final
	This method returns the maximal extension on the local plane, i.e. More...
virtual MsgStream &	dump (MsgStream &sl) const override
	Output Method for MsgStream. More...
virtual std::ostream &	dump (std::ostream &sl) const override
	Output Method for std::ostream. More...
virtual bool	operator!= (const SurfaceBounds &sb) const
	Non-Equality operator. More...

Protected Member Functions
void	swap (double &b1, double &b2)
	Swap method to be called from DiscBounds or TrapezoidalBounds. More...
virtual void	initCache ()
	virtual initCache method for object persistency More...

Private Attributes
std::vector< TDD_real_t >	m_boundValues

Detailed Description

Bounds for a triangular, planar surface.

Author

sarka.nosp@m..tod.nosp@m.orova.nosp@m.@cer.nosp@m.n.ch

Definition at line 38 of file TriangleBounds.h.

Member Enumeration Documentation

BoundsType

enum Trk::SurfaceBounds::BoundsType

inherited

This enumerator simplifies the persistency, by saving a dynamic_cast to happen.

Other is reserved for the GeometrySurfaces implementation.

Enumerator
Cone
Cylinder
Diamond
Disc
Ellipse
Rectangle
RotatedTrapezoid
Trapezoid
Triangle
DiscTrapezoidal
Annulus
Other

Definition at line 58 of file SurfaceBounds.h.

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

BoundValues

enum Trk::TriangleBounds::BoundValues

for readability

Enumerator
bv_x1
bv_y1
bv_x2
bv_y2
bv_x3
bv_y3
bv_length

Definition at line 43 of file TriangleBounds.h.

  {
    bv_x1 = 0,
    bv_y1 = 1,
    bv_x2 = 2,
    bv_y2 = 3,
    bv_x3 = 4,
    bv_y3 = 5,
    bv_length = 6
  };

Constructor & Destructor Documentation

TriangleBounds() [1/6]

Trk::TriangleBounds::TriangleBounds

(

)

Default Constructor - needed for persistency.

Definition at line 19 of file TriangleBounds.cxx.

  : m_boundValues(TriangleBounds::bv_length, 0.)
{}

TriangleBounds() [2/6]

Trk::TriangleBounds::TriangleBounds ( const TriangleBounds &  tribo )

default

Copy constructor.

TriangleBounds() [3/6]

Trk::TriangleBounds::TriangleBounds ( TriangleBounds &&  tribo )

defaultnoexcept

Move constructor.

~TriangleBounds()

virtual Trk::TriangleBounds::~TriangleBounds ( )

virtualdefault

Destructor.

TriangleBounds() [4/6]

Trk::TriangleBounds::TriangleBounds

(

const std::vector< std::pair< float, float >> &

vertices

)

Constructor with coordinates of vertices - floats.

Definition at line 24 of file TriangleBounds.cxx.

  : m_boundValues(TriangleBounds::bv_length, 0.)
{
  size_t ib = 0;
  for (const std::pair<float, float>& p : vertices) {
    m_boundValues[2 * ib] = p.first;
    m_boundValues[2 * ib + 1] = p.second;
    if (ib == 2)
      break;
    ++ib;
  }
}

TriangleBounds() [5/6]

ATH_FLATTEN Trk::TriangleBounds::TriangleBounds

(

const std::vector< std::pair< double, double >> &

vertices

)

Constructor with coordinates of vertices - double.

Definition at line 40 of file TriangleBounds.cxx.

  : m_boundValues(TriangleBounds::bv_length, 0.)
{
  size_t ib = 0;
  for (const std::pair<double, double>& p : vertices) {
    m_boundValues[2 * ib] = p.first;
    m_boundValues[2 * ib + 1] = p.second;
    if (ib == 2)
      break;
    ++ib;
  }
}

TriangleBounds() [6/6]

Trk::TriangleBounds::TriangleBounds	(	const Amg::Vector2D &	p1,
		const Amg::Vector2D &	p2,
		const Amg::Vector2D &	p3
	)

Constructor from three 2 Vectors.

Definition at line 54 of file TriangleBounds.cxx.

  : m_boundValues(TriangleBounds::bv_length, 0.)
{
  m_boundValues[TriangleBounds::bv_x1] = p1.x();
  m_boundValues[TriangleBounds::bv_y1] = p1.y();
  m_boundValues[TriangleBounds::bv_x2] = p2.x();
  m_boundValues[TriangleBounds::bv_y2] = p2.y();
  m_boundValues[TriangleBounds::bv_x3] = p3.x();
  m_boundValues[TriangleBounds::bv_y3] = p3.y();
}

Member Function Documentation

clone()

virtual TriangleBounds* Trk::TriangleBounds::clone ( ) const

overridevirtual

Virtual constructor.

Implements Trk::SurfaceBounds.

dump() [1/2]

MsgStream & Trk::TriangleBounds::dump ( MsgStream &  sl ) const

overridevirtual

Output Method for MsgStream.

Implements Trk::SurfaceBounds.

Definition at line 229 of file TriangleBounds.cxx.

{
  sl << std::setiosflags(std::ios::fixed);
  sl << std::setprecision(7);
  sl << "Trk::TriangleBounds:  generating vertices (X, Y) " << '\n';
  sl << "(" << m_boundValues[TriangleBounds::bv_x1] << " , " << m_boundValues[TriangleBounds::bv_y1] << ") " << '\n';
  sl << "(" << m_boundValues[TriangleBounds::bv_x2] << " , " << m_boundValues[TriangleBounds::bv_y2] << ") " << '\n';
  sl << "(" << m_boundValues[TriangleBounds::bv_x3] << " , " << m_boundValues[TriangleBounds::bv_y3] << ") ";
  sl << std::setprecision(-1);
  return sl;
}

dump() [2/2]

std::ostream & Trk::TriangleBounds::dump ( std::ostream &  sl ) const

overridevirtual

Output Method for std::ostream.

Implements Trk::SurfaceBounds.

Definition at line 242 of file TriangleBounds.cxx.

{
  sl << std::setiosflags(std::ios::fixed);
  sl << std::setprecision(7);
  sl << "Trk::TriangleBounds:  generating vertices (X, Y)";
  sl << "(" << m_boundValues[TriangleBounds::bv_x1] << " , " << m_boundValues[TriangleBounds::bv_y1] << ") " << '\n';
  sl << "(" << m_boundValues[TriangleBounds::bv_x2] << " , " << m_boundValues[TriangleBounds::bv_y2] << ") " << '\n';
  sl << "(" << m_boundValues[TriangleBounds::bv_x3] << " , " << m_boundValues[TriangleBounds::bv_y3] << ") ";
  sl << std::setprecision(-1);
  return sl;
}

initCache()

virtual void Trk::SurfaceBounds::initCache ( )

inlineprotectedvirtualinherited

virtual initCache method for object persistency

Reimplemented in Trk::RotatedTrapezoidBounds, Trk::ConeBounds, Trk::RotatedDiamondBounds, and Trk::DiamondBounds.

Definition at line 129 of file SurfaceBounds.h.

{}

inside() [1/2]

bool Trk::TriangleBounds::inside ( const Amg::Vector2D &  locpo, const BoundaryCheck &  bchk  ) const

finaloverridevirtual

Implements Trk::SurfaceBounds.

Definition at line 122 of file TriangleBounds.cxx.

{
  if (bchk.bcType == 0)
    return TriangleBounds::inside(
      locpo, bchk.toleranceLoc1, bchk.toleranceLoc2);
 
  // a fast FALSE
  double fabsR = std::sqrt(locpo[Trk::locX] * locpo[Trk::locX] +
                           locpo[Trk::locY] * locpo[Trk::locY]);
  double max_ell = bchk.lCovariance(0, 0) > bchk.lCovariance(1, 1)
                     ? bchk.lCovariance(0, 0)
                     : bchk.lCovariance(1, 1);
  double limit = bchk.nSigmas * std::sqrt(max_ell);
  double r_max = TriangleBounds::r();
  if (fabsR > (r_max + limit))
    return false;
 
  // compute KDOP and axes for surface polygon
  std::vector<KDOP> elementKDOP(3);
  std::vector<Amg::Vector2D> elementP(3);
  float theta =
    (bchk.lCovariance(1, 0) != 0 &&
     (bchk.lCovariance(1, 1) - bchk.lCovariance(0, 0)) != 0)
      ? .5 * bchk.FastArcTan(2 * bchk.lCovariance(1, 0) /
                             (bchk.lCovariance(1, 1) - bchk.lCovariance(0, 0)))
      : 0.;
  sincosCache scResult = bchk.FastSinCos(theta);
  AmgMatrix(2, 2) rotMatrix;
  rotMatrix << scResult.cosC, scResult.sinC, -scResult.sinC, scResult.cosC;
  AmgMatrix(2, 2) normal;
  // cppcheck-suppress constStatement
  normal << 0, -1, 1, 0;
  // ellipse is always at (0,0), surface is moved to ellipse position and then
  // rotated
  Amg::Vector2D p;
  p << m_boundValues[TriangleBounds::bv_x1],
    m_boundValues[TriangleBounds::bv_y1];
  elementP[0] = (rotMatrix * (p - locpo));
  p << m_boundValues[TriangleBounds::bv_x2],
    m_boundValues[TriangleBounds::bv_y2];
  elementP[1] = (rotMatrix * (p - locpo));
  p << m_boundValues[TriangleBounds::bv_x3],
    m_boundValues[TriangleBounds::bv_y3];
  elementP[2] = (rotMatrix * (p - locpo));
  std::vector<Amg::Vector2D> axis = { normal * (elementP[1] - elementP[0]),
                                      normal * (elementP[2] - elementP[1]),
                                      normal * (elementP[2] - elementP[0]) };
  bchk.ComputeKDOP(elementP, axis, elementKDOP);
  // compute KDOP for error ellipse
  std::vector<KDOP> errelipseKDOP(3);
  bchk.ComputeKDOP(bchk.EllipseToPoly(3), axis, errelipseKDOP);
  // check if KDOPs overlap and return result
  return bchk.TestKDOPKDOP(elementKDOP, errelipseKDOP);
}

inside() [2/2]

bool Trk::TriangleBounds::inside ( const Amg::Vector2D &  locpo, double  tol1 = 0., double  tol2 = 0.  ) const

finaloverridevirtual

This method checks if the provided local coordinates are inside the surface bounds.

Implements Trk::SurfaceBounds.

Definition at line 76 of file TriangleBounds.cxx.

{
  std::pair<double, double> locB(m_boundValues[TriangleBounds::bv_x2] -
                                   m_boundValues[TriangleBounds::bv_x1],
                                 m_boundValues[TriangleBounds::bv_y2] -
                                   m_boundValues[TriangleBounds::bv_y1]);
  std::pair<double, double> locT(
    m_boundValues[TriangleBounds::bv_x3] - locpo[0],
    m_boundValues[TriangleBounds::bv_y3] - locpo[1]);
  std::pair<double, double> locV(
    m_boundValues[TriangleBounds::bv_x1] - locpo[0],
    m_boundValues[TriangleBounds::bv_y1] - locpo[1]);
 
  // special case :: third vertex ?
  if (locT.first * locT.first + locT.second * locT.second < tol1 * tol1)
    return true;
 
  // special case : lies on base ?
  double db = locB.first * locV.second - locB.second * locV.first;
  if (std::abs(db) < tol1) {
    double a =
      (locB.first != 0) ? -locV.first / locB.first : -locV.second / locB.second;
    return a > -tol2 && a - 1. < tol2;
  }
 
  double dn = locB.first * locT.second - locB.second * locT.first;
 
  if (std::abs(dn) > std::abs(tol1)) {
    double t = (locB.first * locV.second - locB.second * locV.first) / dn;
    if (t > 0.)
      return false;
 
    double a = (locB.first != 0.)
                 ? (t * locT.first - locV.first) / locB.first
                 : (t * locT.second - locV.second) / locB.second;
    if (a < -tol2 || a - 1. > tol2)
      return false;
  } else {
    return false;
  }
  return true;
}

insideLoc1()

virtual bool Trk::TriangleBounds::insideLoc1 ( const Amg::Vector2D &  locpo, double  tol1 = 0.  ) const

finaloverridevirtual

This method checks inside bounds in loc1.

• loc1/loc2 correspond to the natural coordinates of the surface

Implements Trk::SurfaceBounds.

insideLoc2()

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

finaloverridevirtual

This method checks inside bounds in loc2.

• loc1/loc2 correspond to the natural coordinates of the surface

Implements Trk::SurfaceBounds.

minDistance()

double Trk::TriangleBounds::minDistance ( const Amg::Vector2D &  pos ) const

finaloverridevirtual

Minimal distance to boundary ( > 0 if outside and <=0 if inside)

Implements Trk::SurfaceBounds.

Definition at line 179 of file TriangleBounds.cxx.

{
  const int Np = 3;
 
  double X[3] = { m_boundValues[TriangleBounds::bv_x1],
                  m_boundValues[TriangleBounds::bv_x2],
                  m_boundValues[TriangleBounds::bv_x3] };
  double Y[3] = { m_boundValues[TriangleBounds::bv_y1],
                  m_boundValues[TriangleBounds::bv_y2],
                  m_boundValues[TriangleBounds::bv_y3] };
 
  double dm = 1.e+20;
  double Ao = 0.;
  bool in = true;
 
  for (int i = 0; i != Np; ++i) {
 
    int j = (i == Np-1 ? 0 : i+1);
 
    double x = X[i] - pos[0];
    double y = Y[i] - pos[1];
    double dx = X[j] - X[i];
    double dy = Y[j] - Y[i];
    double A = x * dy - y * dx;
    double S = -(x * dx + y * dy);
 
    if (S <= 0.) {
      double d = x * x + y * y;
      if (d < dm)
        dm = d;
    } else {
      double a = dx * dx + dy * dy;
      if (S <= a) {
        double d = (A * A) / a;
        if (d < dm)
          dm = d;
      }
    }
    if (i && in && Ao * A < 0.)
      in = false;
    Ao = A;
  }
  if (in){
    return -sqrt(dm);
  }
  return sqrt(dm);
}

operator!=()

bool Trk::SurfaceBounds::operator!= ( const SurfaceBounds &  sb ) const

inlinevirtualinherited

Non-Equality operator.

Reimplemented in Trk::InvalidBounds.

Definition at line 141 of file SurfaceBounds.h.

{
  return !((*this) == sb);
}

operator=() [1/2]

TriangleBounds& Trk::TriangleBounds::operator= ( const TriangleBounds &  recbo )

default

Assignment Operator.

operator=() [2/2]

TriangleBounds& Trk::TriangleBounds::operator= ( TriangleBounds &&  recbo )

defaultnoexcept

Move Assignment Operator.

operator==()

bool Trk::TriangleBounds::operator== ( const SurfaceBounds &  sbo ) const

overridevirtual

Equality operator.

Implements Trk::SurfaceBounds.

Definition at line 66 of file TriangleBounds.cxx.

{
  // check the type first not to compare apples with oranges
  const Trk::TriangleBounds* tribo = dynamic_cast<const Trk::TriangleBounds*>(&sbo);
  if (!tribo)
    return false;
  return (m_boundValues == tribo->m_boundValues);
}

r()

virtual double Trk::TriangleBounds::r ( ) const

finaloverridevirtual

This method returns the maximal extension on the local plane, i.e.

$s\sqrt{h_{\phi}^2 + h_{\eta}^2}$

Implements Trk::SurfaceBounds.

swap()

void Trk::SurfaceBounds::swap ( double &  b1, double &  b2  )

inlineprotectedinherited

Swap method to be called from DiscBounds or TrapezoidalBounds.

Definition at line 133 of file SurfaceBounds.h.

{
  double tmp = b1;
  b1 = b2;
  b2 = tmp;
}

type()

virtual BoundsType Trk::TriangleBounds::type ( ) const

inlinefinaloverridevirtual

Return the type of the bounds for persistency.

Implements Trk::SurfaceBounds.

Definition at line 88 of file TriangleBounds.h.

{ return SurfaceBounds::Triangle; }

vertices()

std::vector<std::pair<TDD_real_t, TDD_real_t> > Trk::TriangleBounds::vertices

(

)

const

This method returns the coordinates of vertices.

Member Data Documentation

m_boundValues

std::vector<TDD_real_t> Trk::TriangleBounds::m_boundValues

private

Definition at line 118 of file TriangleBounds.h.

---

The documentation for this class was generated from the following files:

• TriangleBounds.h
• TriangleBounds.cxx