#include <TrapezoidVolumeBounds.h>

Inheritance diagram for Trk::TrapezoidVolumeBounds:

Collaboration diagram for Trk::TrapezoidVolumeBounds:

Public Member Functions
	TrapezoidVolumeBounds ()
	Default Constructor. More...

	TrapezoidVolumeBounds (double minhlenghtx, double maxhlengthx, double hlenghty, double hlengthz)
	Constructor - the trapezoid boundaries (symmetric trapezoid) More...

	TrapezoidVolumeBounds (double minhlenghtx, double hlenghty, double hlengthz, double alpha, double beta)
	Constructor - the trapezoid boundaries (arbitrary trapezoid) More...

	TrapezoidVolumeBounds (const TrapezoidVolumeBounds &bobo)
	Copy Constructor. More...

virtual	~TrapezoidVolumeBounds ()
	Destructor. More...

TrapezoidVolumeBounds &	operator= (const TrapezoidVolumeBounds &bobo)
	Assignment operator. More...

TrapezoidVolumeBounds *	clone () const override final
	Virtual constructor. More...

bool	inside (const Amg::Vector3D &, double tol=0.) const override final
	This method checks if position in the 3D volume frame is inside the cylinder. More...

const std::vector< const Trk::Surface * > *	decomposeToSurfaces (const Amg::Transform3D &transform) override final
	Method to decompose the Bounds into Surfaces. More...

ObjectAccessor	boundarySurfaceAccessor (const Amg::Vector3D &gp, const Amg::Vector3D &dir, bool forceInside=false) const override
	Provide accessor for BoundarySurfaces. More...

double	minHalflengthX () const
	This method returns the minimal halflength in local x. More...

double	maxHalflengthX () const
	This method returns the maximal halflength in local x. More...

double	halflengthY () const
	This method returns the halflength in local y. More...

double	halflengthZ () const
	This method returns the halflength in local z. More...

double	alpha () const
	This method returns the opening angle in point A (negative local x) More...

double	beta () const
	This method returns the opening angle in point B (negative local x) More...

MsgStream &	dump (MsgStream &sl) const override final
	Output Method for MsgStream. More...

std::ostream &	dump (std::ostream &sl) const override final
	Output Method for std::ostream. More...

Private Member Functions
TrapezoidBounds *	faceXYTrapezoidBounds () const
	This method returns the associated TrapezoidBounds of the face PlaneSurface parallel to local xy plane. More...

RectangleBounds *	faceAlphaRectangleBounds () const
	This method returns the associated RecantleBounds of the face PlaneSurface attached to alpha (negative local x) More...

RectangleBounds *	faceBetaRectangleBounds () const
	This method returns the associated RecantleBounds of the face PlaneSurface attached to beta (positive local x) More...

RectangleBounds *	faceZXRectangleBoundsBottom () const
	This method returns the associated RecantleBounds of the face PlaneSurface parallel to local zx plane, negative local y. More...

RectangleBounds *	faceZXRectangleBoundsTop () const
	This method returns the associated RecantleBounds of the face PlaneSurface parallel to local zx plane, positive local y. More...

Private Attributes
double	m_minHalfX
	minimal halflength in x More...

double	m_maxHalfX
	maximal halflength in x More...

double	m_halfY
	halflength in y More...

double	m_halfZ
	halflength in z More...

double	m_alpha
	opening angle alpha (in point A) More...

double	m_beta
	opening angle beta (in point B) More...

SixObjectsAccessor	m_objectAccessor
	There's only one single object Acessor for the moment has to be implemented if Cuboids are used more widely. More...

Detailed Description

Bounds for a trapezoidal shaped Volume, the decomposeToSurfaces method creates a vector of 6 surfaces:

BoundarySurfaceFace [index]:

- negativeFaceXY     [0] : Trazpezoidal Trk::PlaneSurface,
                           parallel to \f$ xy \f$ plane at negative \f$ z

\( - positiveFaceXY [1] : Trazpezoidal Trk::PlaneSurface, parallel to \) xy \( plane at positive \) z \( - trapezoidFaceAlpha [2] : Rectangular Trk::PlaneSurface, attached to [0] and [1] at negative \) x \( (associated to alpha) - trapezoidFaceBeta [3] : Rectangular Trk::PlaneSurface, attached to [0] and [1] at positive \) x \( (associated to beta) - negativeFaceZX [4] : Rectangular Trk::PlaneSurface, parallel to \) zx \( plane at negative \) y \( - positiveFaceZX [5] : Rectangular Trk::PlaneSurface, parallel to \) zx \( plane at positive \) y

Definition at line 57 of file TrapezoidVolumeBounds.h.

Constructor & Destructor Documentation

◆ TrapezoidVolumeBounds() [1/4]

Trk::TrapezoidVolumeBounds::TrapezoidVolumeBounds ( )

Default Constructor.

Definition at line 26 of file TrapezoidVolumeBounds.cxx.

   : VolumeBounds()
   , m_minHalfX(0.)
   , m_maxHalfX(0.)
   , m_halfY(0.)
   , m_halfZ(0.)
   , m_alpha(0.)
   , m_beta(0.)
   , m_objectAccessor()
 {}

◆ TrapezoidVolumeBounds() [2/4]

Trk::TrapezoidVolumeBounds::TrapezoidVolumeBounds	(	double	minhlenghtx,
		double	maxhlengthx,
		double	hlenghty,
		double	hlengthz
	)

Constructor - the trapezoid boundaries (symmetric trapezoid)

Definition at line 37 of file TrapezoidVolumeBounds.cxx.

   : VolumeBounds()
   , m_minHalfX(minhalex)
   , m_maxHalfX(maxhalex)
   , m_halfY(haley)
   , m_halfZ(halez)
   , m_alpha(0.)
   , m_beta(0.)
   , m_objectAccessor()
 {
   m_alpha = atan((m_maxHalfX - m_minHalfX) / 2 / m_halfY) + 0.5 * M_PI;
   m_beta = m_alpha;
 }

◆ TrapezoidVolumeBounds() [3/4]

Trk::TrapezoidVolumeBounds::TrapezoidVolumeBounds	(	double	minhlenghtx,
		double	hlenghty,
		double	hlengthz,
		double	alpha,
		double	beta
	)

Constructor - the trapezoid boundaries (arbitrary trapezoid)

Definition at line 55 of file TrapezoidVolumeBounds.cxx.

   : VolumeBounds()
   , m_minHalfX(minhalex)
   , m_maxHalfX(0.)
   , m_halfY(haley)
   , m_halfZ(halez)
   , m_alpha(alpha)
   , m_beta(beta)
   , m_objectAccessor()
 {
   double gamma = (alpha > beta) ? (alpha - 0.5 * M_PI) : (beta - 0.5 * M_PI);
   m_maxHalfX = m_minHalfX + (2. * m_halfY) * tan(gamma);
 }

◆ TrapezoidVolumeBounds() [4/4]

Trk::TrapezoidVolumeBounds::TrapezoidVolumeBounds ( const TrapezoidVolumeBounds & bobo )

Copy Constructor.

Definition at line 74 of file TrapezoidVolumeBounds.cxx.

   : VolumeBounds()
   , m_minHalfX(trabo.m_minHalfX)
   , m_maxHalfX(trabo.m_maxHalfX)
   , m_halfY(trabo.m_halfY)
   , m_halfZ(trabo.m_halfZ)
   , m_alpha(trabo.m_alpha)
   , m_beta(trabo.m_beta)
   , m_objectAccessor(trabo.m_objectAccessor)
 {}

◆ ~TrapezoidVolumeBounds()

Trk::TrapezoidVolumeBounds::~TrapezoidVolumeBounds ( )

virtualdefault

Destructor.

Member Function Documentation

◆ alpha()

double Trk::TrapezoidVolumeBounds::alpha ( ) const

inline

This method returns the opening angle in point A (negative local x)

Definition at line 165 of file TrapezoidVolumeBounds.h.

165 { return m_alpha; }

◆ beta()

double Trk::TrapezoidVolumeBounds::beta ( ) const

inline

This method returns the opening angle in point B (negative local x)

Definition at line 166 of file TrapezoidVolumeBounds.h.

166 { return m_beta; }

◆ boundarySurfaceAccessor()

ObjectAccessor Trk::TrapezoidVolumeBounds::boundarySurfaceAccessor	(	const Amg::Vector3D &	gp,
		const Amg::Vector3D &	dir,
		bool	forceInside = `false`
	)		const

inlineoverridevirtual

Provide accessor for BoundarySurfaces.

Implements Trk::VolumeBounds.

Definition at line 168 of file TrapezoidVolumeBounds.h.

                                                         {
   return Trk::ObjectAccessor(m_objectAccessor);
 }

◆ clone()

TrapezoidVolumeBounds * Trk::TrapezoidVolumeBounds::clone ( ) const

inlinefinaloverridevirtual

Virtual constructor.

Implements Trk::VolumeBounds.

Definition at line 153 of file TrapezoidVolumeBounds.h.

                                                                  {
   return new TrapezoidVolumeBounds(*this);
 }

◆ decomposeToSurfaces()

const std::vector< const Trk::Surface * > * Trk::TrapezoidVolumeBounds::decomposeToSurfaces ( const Amg::Transform3D & transform )

finaloverridevirtual

Method to decompose the Bounds into Surfaces.

Implements Trk::VolumeBounds.

Definition at line 104 of file TrapezoidVolumeBounds.cxx.

 {
   std::vector<const Trk::Surface*>* retsf =
     new std::vector<const Trk::Surface*>;
  
   // face surfaces xy
   Amg::RotationMatrix3D trapezoidRotation(transform.rotation());
   Amg::Vector3D trapezoidX(trapezoidRotation.col(0));
   Amg::Vector3D trapezoidY(trapezoidRotation.col(1));
   Amg::Vector3D trapezoidZ(trapezoidRotation.col(2));
   Amg::Vector3D trapezoidCenter(transform.translation());
  
   //   (1) - at negative local z
   retsf->push_back(new Trk::PlaneSurface(
     Amg::Transform3D(
       transform *
       Amg::AngleAxis3D(180 * Gaudi::Units::deg, Amg::Vector3D(0., 1., 0.)) *
       Amg::Translation3D(Amg::Vector3D(0., 0., this->halflengthZ()))),
     this->faceXYTrapezoidBounds()));
   //   (2) - at positive local z
   retsf->push_back(new Trk::PlaneSurface(
     Amg::Transform3D(
       transform *
       Amg::Translation3D(Amg::Vector3D(0., 0., this->halflengthZ()))),
     this->faceXYTrapezoidBounds()));
   // face surfaces yz
   // transmute cyclical
   //   (3) - at point A, attached to alpha opening angle
   Amg::Vector3D A(
     this->minHalflengthX(), this->halflengthY(), trapezoidCenter.z());
   Amg::RotationMatrix3D alphaZRotation =
     (s_idRotation *
      Amg::AngleAxis3D(this->alpha() - 0.5 * M_PI, Amg::Vector3D(0., 0., 1.)))
       .toRotationMatrix();
   // CLHEP::HepRotation  alphaRotation(alphaZRotation*trapezoidRotation);
   Amg::RotationMatrix3D faceAlphaRotation;
   faceAlphaRotation.col(0) = alphaZRotation.col(1);
   faceAlphaRotation.col(1) = -alphaZRotation.col(2);
   faceAlphaRotation.col(2) = -alphaZRotation.col(0);
   RectangleBounds* faceAlphaBounds = this->faceAlphaRectangleBounds();
   // Amg::Vector3D
   // faceAlphaPosition(A+faceAlphaRotation.colX()*faceAlphaBounds->halflengthX());
   Amg::Vector3D faceAlphaPosition0(
     -0.5 * (this->minHalflengthX() + this->maxHalflengthX()), 0., 0.);
   Amg::Vector3D faceAlphaPosition = transform * faceAlphaPosition0;
   retsf->push_back(new Trk::PlaneSurface(
     Amg::Transform3D(
       (trapezoidRotation * faceAlphaRotation) *
       Amg::Translation3D(faceAlphaPosition)),
     faceAlphaBounds));
   //   (4) - at point B, attached to beta opening angle
   Amg::Vector3D B(
     this->minHalflengthX(), -this->halflengthY(), trapezoidCenter.z());
   Amg::RotationMatrix3D betaZRotation =
     (s_idRotation *
      Amg::AngleAxis3D(-(this->beta() - 0.5 * M_PI), Amg::Vector3D(0., 0., 1.)))
       .toRotationMatrix();
   // CLHEP::HepRotation  betaRotation(betaZRotation*trapezoidRotation);
   Amg::RotationMatrix3D faceBetaRotation;
   faceBetaRotation.col(0) = betaZRotation.col(1);
   faceBetaRotation.col(1) = betaZRotation.col(2);
   faceBetaRotation.col(2) = betaZRotation.col(0);
   RectangleBounds* faceBetaBounds = this->faceBetaRectangleBounds();
   // Amg::Vector3D
   // faceBetaPosition(B+faceBetaRotation.colX()*faceBetaBounds->halflengthX());
   Amg::Vector3D faceBetaPosition0(
     0.5 * (this->minHalflengthX() + this->maxHalflengthX()), 0., 0.);
   Amg::Vector3D faceBetaPosition = transform * faceBetaPosition0;
   retsf->push_back(new Trk::PlaneSurface(
     Amg::Transform3D(
       trapezoidRotation * faceBetaRotation *
       Amg::Translation3D(faceBetaPosition)),
     faceBetaBounds));
   // face surfaces zx
   //   (5) - at negative local x
   retsf->push_back(new Trk::PlaneSurface(
     Amg::Transform3D(
       transform *
       Amg::AngleAxis3D(180. * Gaudi::Units::deg, Amg::Vector3D(1., 0., 0.)) *
       Amg::Translation3D(Amg::Vector3D(0., this->halflengthY(), 0.)) *
       Amg::AngleAxis3D(-90 * Gaudi::Units::deg, Amg::Vector3D(0., 1., 0.)) *
       Amg::AngleAxis3D(-90. * Gaudi::Units::deg, Amg::Vector3D(1., 0., 0.))),
     this->faceZXRectangleBoundsBottom()));
   //   (6) - at positive local x
   retsf->push_back(new Trk::PlaneSurface(
     Amg::Transform3D(
       transform *
       Amg::Translation3D(Amg::Vector3D(0., this->halflengthY(), 0.)) *
       Amg::AngleAxis3D(-90 * Gaudi::Units::deg, Amg::Vector3D(0., 1., 0.)) *
       Amg::AngleAxis3D(-90. * Gaudi::Units::deg, Amg::Vector3D(1., 0., 0.))),
     this->faceZXRectangleBoundsTop()));
  
   return retsf;
 }

◆ dump() [1/2]

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

finaloverridevirtual

Output Method for MsgStream.

Implements Trk::VolumeBounds.

Definition at line 251 of file TrapezoidVolumeBounds.cxx.

 {
   std::stringstream temp_sl;
   temp_sl << std::setiosflags(std::ios::fixed);
   temp_sl << std::setprecision(7);
   temp_sl
     << "Trk::TrapezoidVolumeBounds: (minhalfX, halfY, halfZ, alpha, beta) = ";
   temp_sl << "(" << m_minHalfX << ", " << m_halfY << ", " << m_halfZ;
   temp_sl << ", " << m_alpha << ", " << m_beta << ")";
   sl << temp_sl.str();
   return sl;
 }

◆ dump() [2/2]

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

finaloverridevirtual

Output Method for std::ostream.

Implements Trk::VolumeBounds.

Definition at line 265 of file TrapezoidVolumeBounds.cxx.

 {
   std::stringstream temp_sl;
   temp_sl << std::setiosflags(std::ios::fixed);
   temp_sl << std::setprecision(7);
   temp_sl
     << "Trk::TrapezoidVolumeBounds: (minhalfX, halfY, halfZ, alpha, beta) = ";
   temp_sl << "(" << m_minHalfX << ", " << m_halfY << ", " << m_halfZ;
   temp_sl << ", " << m_alpha << ", " << m_beta << ")";
   sl << temp_sl.str();
   return sl;
 }

◆ faceAlphaRectangleBounds()

Trk::RectangleBounds * Trk::TrapezoidVolumeBounds::faceAlphaRectangleBounds ( ) const

private

This method returns the associated RecantleBounds of the face PlaneSurface attached to alpha (negative local x)

Definition at line 209 of file TrapezoidVolumeBounds.cxx.

 {
   return new Trk::RectangleBounds(m_halfY / cos(m_alpha - 0.5 * M_PI), m_halfZ);
 }

◆ faceBetaRectangleBounds()

Trk::RectangleBounds * Trk::TrapezoidVolumeBounds::faceBetaRectangleBounds ( ) const

private

This method returns the associated RecantleBounds of the face PlaneSurface attached to beta (positive local x)

Definition at line 215 of file TrapezoidVolumeBounds.cxx.

 {
   return new Trk::RectangleBounds(m_halfY / cos(m_beta - 0.5 * M_PI), m_halfZ);
 }

◆ faceXYTrapezoidBounds()

Trk::TrapezoidBounds * Trk::TrapezoidVolumeBounds::faceXYTrapezoidBounds ( ) const

private

This method returns the associated TrapezoidBounds of the face PlaneSurface parallel to local xy plane.

Definition at line 202 of file TrapezoidVolumeBounds.cxx.

 {
   // return new Trk::TrapezoidBounds(m_minHalfX,m_halfY, m_alpha, m_beta);
   return new Trk::TrapezoidBounds(m_minHalfX, m_maxHalfX, m_halfY);
 }

◆ faceZXRectangleBoundsBottom()

Trk::RectangleBounds * Trk::TrapezoidVolumeBounds::faceZXRectangleBoundsBottom ( ) const

private

This method returns the associated RecantleBounds of the face PlaneSurface parallel to local zx plane, negative local y.

Definition at line 221 of file TrapezoidVolumeBounds.cxx.

 {
   return new Trk::RectangleBounds(m_halfZ, m_minHalfX);
 }

◆ faceZXRectangleBoundsTop()

Trk::RectangleBounds * Trk::TrapezoidVolumeBounds::faceZXRectangleBoundsTop ( ) const

private

This method returns the associated RecantleBounds of the face PlaneSurface parallel to local zx plane, positive local y.

Definition at line 227 of file TrapezoidVolumeBounds.cxx.

 {
   // double delta = (m_alpha < m_beta) ? m_alpha - M_PI/2. : m_beta - M_PI/2.;
   // return new Trk::RectangleBounds(m_halfZ,
   // 0.5*(m_minHalfX+m_minHalfX+2.*m_halfY/cos(delta)));
   return new Trk::RectangleBounds(m_halfZ, m_maxHalfX);
 }

◆ halflengthY()

double Trk::TrapezoidVolumeBounds::halflengthY ( ) const

inline

This method returns the halflength in local y.

Definition at line 163 of file TrapezoidVolumeBounds.h.

163 { return m_halfY; }

◆ halflengthZ()

double Trk::TrapezoidVolumeBounds::halflengthZ ( ) const

inline

This method returns the halflength in local z.

Definition at line 164 of file TrapezoidVolumeBounds.h.

164 { return m_halfZ; }

◆ inside()

bool Trk::TrapezoidVolumeBounds::inside	(	const Amg::Vector3D &	pos,
		double	tol = `0.`
	)		const

finaloverridevirtual

This method checks if position in the 3D volume frame is inside the cylinder.

Implements Trk::VolumeBounds.

Definition at line 236 of file TrapezoidVolumeBounds.cxx.

 {
   if (std::abs(pos.z()) > m_halfZ + tol)
     return false;
   if (std::abs(pos.y()) > m_halfY + tol)
     return false;
   Trk::TrapezoidBounds* faceXYBounds = this->faceXYTrapezoidBounds();
   Amg::Vector2D locp(pos.x(), pos.y());
   bool inside(faceXYBounds->inside(locp, tol, tol));
   delete faceXYBounds;
   return inside;
 }

◆ maxHalflengthX()

double Trk::TrapezoidVolumeBounds::maxHalflengthX ( ) const

inline

This method returns the maximal halflength in local x.

Definition at line 160 of file TrapezoidVolumeBounds.h.

                                                           {
   return m_maxHalfX;
 }

◆ minHalflengthX()

double Trk::TrapezoidVolumeBounds::minHalflengthX ( ) const

inline

This method returns the minimal halflength in local x.

Definition at line 157 of file TrapezoidVolumeBounds.h.

                                                           {
   return m_minHalfX;
 }

◆ operator=()

Trk::TrapezoidVolumeBounds & Trk::TrapezoidVolumeBounds::operator= ( const TrapezoidVolumeBounds & bobo )

Assignment operator.

Definition at line 89 of file TrapezoidVolumeBounds.cxx.

 {
   if (this != &trabo) {
     m_minHalfX = trabo.m_minHalfX;
     m_maxHalfX = trabo.m_maxHalfX;
     m_halfY = trabo.m_halfY;
     m_halfZ = trabo.m_halfZ;
     m_alpha = trabo.m_alpha;
     m_beta = trabo.m_beta;
     m_objectAccessor = trabo.m_objectAccessor;
   }
   return *this;
 }

Member Data Documentation

◆ m_alpha

double Trk::TrapezoidVolumeBounds::m_alpha

private

opening angle alpha (in point A)

Definition at line 145 of file TrapezoidVolumeBounds.h.

◆ m_beta

double Trk::TrapezoidVolumeBounds::m_beta

private

opening angle beta (in point B)

Definition at line 146 of file TrapezoidVolumeBounds.h.

◆ m_halfY

double Trk::TrapezoidVolumeBounds::m_halfY

private

halflength in y

Definition at line 142 of file TrapezoidVolumeBounds.h.

◆ m_halfZ

double Trk::TrapezoidVolumeBounds::m_halfZ

private

halflength in z

Definition at line 143 of file TrapezoidVolumeBounds.h.

◆ m_maxHalfX

double Trk::TrapezoidVolumeBounds::m_maxHalfX

private

maximal halflength in x

Definition at line 141 of file TrapezoidVolumeBounds.h.

◆ m_minHalfX

double Trk::TrapezoidVolumeBounds::m_minHalfX

private

minimal halflength in x

Definition at line 140 of file TrapezoidVolumeBounds.h.

◆ m_objectAccessor

SixObjectsAccessor Trk::TrapezoidVolumeBounds::m_objectAccessor

private

There's only one single object Acessor for the moment has to be implemented if Cuboids are used more widely.

Definition at line 150 of file TrapezoidVolumeBounds.h.

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TrapezoidVolumeBounds() [1/4]

◆ TrapezoidVolumeBounds() [2/4]

◆ TrapezoidVolumeBounds() [3/4]

◆ TrapezoidVolumeBounds() [4/4]

◆ ~TrapezoidVolumeBounds()

Member Function Documentation

◆ alpha()

◆ beta()

◆ boundarySurfaceAccessor()

◆ clone()

◆ decomposeToSurfaces()

◆ dump() [1/2]

◆ dump() [2/2]

◆ faceAlphaRectangleBounds()

◆ faceBetaRectangleBounds()

◆ faceXYTrapezoidBounds()

◆ faceZXRectangleBoundsBottom()

◆ faceZXRectangleBoundsTop()

◆ halflengthY()

◆ halflengthZ()

◆ inside()

◆ maxHalflengthX()

◆ minHalflengthX()

◆ operator=()

Member Data Documentation

◆ m_alpha

◆ m_beta

◆ m_halfY

◆ m_halfZ

◆ m_maxHalfX

◆ m_minHalfX

◆ m_objectAccessor