Trk::CylinderVolumeBounds Class Reference

#include <CylinderVolumeBounds.h>

Inheritance diagram for Trk::CylinderVolumeBounds:

Collaboration diagram for Trk::CylinderVolumeBounds:

Public Member Functions
	CylinderVolumeBounds ()
	Default Constructor. More...
	CylinderVolumeBounds (double radius, double halez)
	Constructor - full cylinder. More...
	CylinderVolumeBounds (double rinner, double router, double halez)
	Constructor - extruded cylinder. More...
	CylinderVolumeBounds (double rinner, double router, double halfPhiSector, double halez)
	Constructor - extruded cylinder segment. More...
	CylinderVolumeBounds (const CylinderVolumeBounds &cylbo)
	Copy Constructor. More...
virtual	~CylinderVolumeBounds ()
	Destructor. More...
CylinderVolumeBounds &	operator= (const CylinderVolumeBounds &cylbo)
	Assignment operator. More...
CylinderVolumeBounds *	clone () const override
	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...
virtual std::vector< std::unique_ptr< Trk::Surface > >	decomposeToSurfaces (const Amg::Transform3D &transform) override final
	Method to decompose the Bounds into boundarySurfaces. More...
ObjectAccessor	boundarySurfaceAccessor (const Amg::Vector3D &gp, const Amg::Vector3D &dir, bool forceInside=false) const override final
	Provide accessor for BoundarySurfaces. More...
double	innerRadius () const
	This method returns the inner radius. More...
double	outerRadius () const
	This method returns the outer radius. More...
double	mediumRadius () const
	This method returns the medium radius. More...
double	deltaRadius () const
	This method returns the delta radius. More...
double	halfPhiSector () const
	This method returns the halfPhiSector angle. More...
double	halflengthZ () const
	This method returns the halflengthZ. 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
std::shared_ptr< CylinderBounds >	innerCylinderBounds () const
	This method returns the associated CylinderBounds of the inner CylinderSurfaces. More...
std::shared_ptr< CylinderBounds >	outerCylinderBounds () const
	This method returns the associated CylinderBounds of the outer CylinderSurfaces. More...
std::shared_ptr< DiscBounds >	bottomDiscBounds () const
	This method returns the associated DiscBounds for the bottom/top DiscSurface. More...
std::shared_ptr< DiscBounds >	topDiscBounds () const
	This method returns the associated DiscBounds for the bottom/top DiscSurface. More...
std::shared_ptr< RectangleBounds >	sectorPlaneBounds () const
	This method returns the associated PlaneBounds limiting a sectoral CylinderVolume. More...
void	createBoundarySurfaceAccessors ()
	Private method to construct the accessors. More...

Private Attributes
double	m_innerRadius
double	m_outerRadius
double	m_halfPhiSector
double	m_halfZ
CylinderVolumeBoundaryAccessors	m_boundaryAccessors
	Accessors for Boundary surface access - static is not possible due to mismatched delete() / free () with TrkMagFieldUtils. More...

Static Private Attributes
static const double	s_numericalStable
	numerical stability More...

Detailed Description

Bounds for a cylindrical Volume, the decomposeToSurfaces method creates a vector of up to 6 surfaces:

case A) 3 Surfaces (full cylindrical tube): BoundarySurfaceFace [index]:

• negativeFaceXY [0] : Trk::DiscSurface with \( r_{inner}=0 \), parallel to \( xy \) plane at negative \( z \)
• positiveFaceXY [1] : Trk::DiscSurface with \( r_{inner}=0 \), parallel to \( xy \) plane at positive \( z \)
• cylinderCover [2] : Trk::CylinderSurface confining the Trk::Volume

case B) 4 Surfaces (tube with inner and outer radius): BoundarySurfaceFace [index]:

• negativeFaceXY [0] : Trk::DiscSurface with \( r_{inner}>0 \), parallel to \( xy \) plane at negative \( z \)
• positiveFaceXY [1] : Trk::DiscSurface with \( r_{inner}>0 \), parallel to \( xy \) plane at positive \( z \)
• tubeOuterCover [2] : Trk::CylinderSurface with \( r = r_{outer} \)
• tubeInnerCover [3] : Trk::CylinderSurface with \( r = r_{inner} \)

case C) 6 Surfaces (sectoral tube with inner and outer radius): BoundarySurfaceFace [index]:

• negativeFaceXY [0] : Trk::DiscSurface with \( r_{inner}>0 \) and \( \phi < \pi \), parallel to \( xy \) plane at negative \( z \)
• positiveFaceXY [1] : Trk::DiscSurface with \( r_{inner}>0 \) and \( \phi < \pi \), parallel to \( xy \) plane at positive \( z \)
• tubeSectorOuterCover [2] : Trk::CylinderSurface with \( r = r_{outer} \)
• tubeSectorInnerCover [3] : Trk::CylinderSurface with \( r = r_{inner} \)
• tubeSectorNegativePhi [4] : Rectangular Trk::PlaneSurface attached to [0] and [1] at negative \( \phi \)
• tubeSectorNegativePhi [5] : Rectangular Trk::PlaneSurface attached to [0] and [1] at positive \( \phi \)

Author

Andre.nosp@m.as.S.nosp@m.alzbu.nosp@m.rger.nosp@m.@cern.nosp@m..ch

Christos Anastopoulos (Athena MT modifications)

Definition at line 70 of file CylinderVolumeBounds.h.

Constructor & Destructor Documentation

CylinderVolumeBounds() [1/5]

Trk::CylinderVolumeBounds::CylinderVolumeBounds

(

)

Default Constructor.

Definition at line 30 of file CylinderVolumeBounds.cxx.

  : VolumeBounds()
  , m_innerRadius(0.)
  , m_outerRadius(0.)
  , m_halfPhiSector(0.)
  , m_halfZ(0.)
  , m_boundaryAccessors()
{}

CylinderVolumeBounds() [2/5]

Trk::CylinderVolumeBounds::CylinderVolumeBounds	(	double	radius,
		double	halez
	)

Constructor - full cylinder.

Definition at line 39 of file CylinderVolumeBounds.cxx.

  : VolumeBounds()
  , m_innerRadius(0.)
  , m_outerRadius(std::abs(radius))
  , m_halfPhiSector(M_PI)
  , m_halfZ(std::abs(halez))
  , m_boundaryAccessors()
{}

CylinderVolumeBounds() [3/5]

Trk::CylinderVolumeBounds::CylinderVolumeBounds	(	double	rinner,
		double	router,
		double	halez
	)

Constructor - extruded cylinder.

Definition at line 48 of file CylinderVolumeBounds.cxx.

  : VolumeBounds()
  , m_innerRadius(std::abs(rinner))
  , m_outerRadius(std::abs(router))
  , m_halfPhiSector(M_PI)
  , m_halfZ(std::abs(halez))
  , m_boundaryAccessors()
{}

CylinderVolumeBounds() [4/5]

Trk::CylinderVolumeBounds::CylinderVolumeBounds	(	double	rinner,
		double	router,
		double	halfPhiSector,
		double	halez
	)

Constructor - extruded cylinder segment.

Definition at line 60 of file CylinderVolumeBounds.cxx.

  : VolumeBounds()
  , m_innerRadius(std::abs(rinner))
  , m_outerRadius(std::abs(router))
  , m_halfPhiSector(haphi)
  , m_halfZ(std::abs(halez))
  , m_boundaryAccessors()
{}

CylinderVolumeBounds() [5/5]

Trk::CylinderVolumeBounds::CylinderVolumeBounds

(

const CylinderVolumeBounds &

cylbo

)

Copy Constructor.

Definition at line 73 of file CylinderVolumeBounds.cxx.

  : VolumeBounds()
  , m_innerRadius(cylbo.m_innerRadius)
  , m_outerRadius(cylbo.m_outerRadius)
  , m_halfPhiSector(cylbo.m_halfPhiSector)
  , m_halfZ(cylbo.m_halfZ)
  , m_boundaryAccessors()
{}

~CylinderVolumeBounds()

Trk::CylinderVolumeBounds::~CylinderVolumeBounds ( )

virtualdefault

Destructor.

Member Function Documentation

bottomDiscBounds()

std::shared_ptr< Trk::DiscBounds > Trk::CylinderVolumeBounds::bottomDiscBounds ( ) const

private

This method returns the associated DiscBounds for the bottom/top DiscSurface.

Definition at line 364 of file CylinderVolumeBounds.cxx.

{
  return std::make_shared<Trk::DiscBounds>(m_innerRadius, m_outerRadius, m_halfPhiSector);
}

boundarySurfaceAccessor()

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

finaloverridevirtual

Provide accessor for BoundarySurfaces.

Implements Trk::VolumeBounds.

Definition at line 163 of file CylinderVolumeBounds.cxx.

{
  // the tube case - most likely
  if (
    m_innerRadius > s_numericalStable &&
    std::abs(m_halfPhiSector - M_PI) < s_numericalStable) {
    // prepare the data
    double posZ = gp.z();
    double posR = gp.perp();
    // difference Vector
    Amg::Vector3D diff(gp + dir);
    // differences
    double deltaZ = diff.z() - posZ;
    double deltaR = diff.perp() - posR;
 
    // the isOnSurface cases + switchers (that change with the cases)
    bool isOnFace = false;
    bool intersectionIndicator = (deltaR > 0.);
    bool choiceIndicator = false;
 
    // on surface or look-a-likes
    // on zMin or slightly outside
    if (
      std::abs(posZ + m_halfZ) < s_numericalStable ||
      (posZ < (-m_halfZ) && deltaZ > 0.)) {
      isOnFace = true;
    } else if (
      std::abs(posZ - m_halfZ) < s_numericalStable ||
      (posZ > m_halfZ && deltaZ < 0.)) {
      // on zMax or slighly outside
      isOnFace = true;
    } else if (
      std::abs(posR - m_innerRadius) < s_numericalStable ||
      (posR < m_innerRadius && deltaR > 0.)) {
      // on innerRadius or slightly outside
      isOnFace = true;
      choiceIndicator =
        (deltaZ > 0.); // the side choice indicator <========== WRONG ==========
    } else if (
      std::abs(posR - m_outerRadius) < s_numericalStable ||
      (posR > m_outerRadius && deltaR < 0.)) {
      // on outRadius or slightly outside
      isOnFace = true;
      choiceIndicator =
        (deltaZ > 0.); // the side choice indicator <========== WRONG ==========
    }
 
    // the onSurface case
    // =================================================================================
    if (isOnFace) {
      if (intersectionIndicator) {
        // intersect the Rmax / with floatint point exception check
        Trk::CylinderIntersector intersectRmax(
          posR, posZ, deltaR != 0 ? deltaZ / deltaR : 0, m_outerRadius);
        double zOfIntersect = intersectRmax.yOfX;
        // now check if the intersect is inside m_halfZ
        if (std::abs(zOfIntersect) <= m_halfZ)
          return {(choiceIndicator || zOfIntersect > 0.)
                      ? m_boundaryAccessors.tubeAccessor(
                            Trk::TubeRincreaseZincrease)
                      : m_boundaryAccessors.tubeAccessor(
                            Trk::TubeRincreaseZdecrease)};
        // if the intersect is outside
        return {
            (choiceIndicator || zOfIntersect > 0.)
                ? m_boundaryAccessors.tubeAccessor(Trk::TubeZincreaseRincrease)
                : m_boundaryAccessors.tubeAccessor(
                      Trk::TubeZdecreaseRincrease)};
      }
      // intersect the Rmin
      Trk::CylinderIntersector intersectRmin(
        posR, posZ, deltaR != 0 ? deltaZ / deltaR : 0, m_innerRadius);
      double zOfIntersect = intersectRmin.yOfX;
      // now check if the intersect is inside m_halfZ
      if (std::abs(zOfIntersect) <= m_halfZ)
        return {
            (choiceIndicator || zOfIntersect > 0.)
                ? m_boundaryAccessors.tubeAccessor(Trk::TubeRdecreaseZincrease)
                : m_boundaryAccessors.tubeAccessor(
                      Trk::TubeRdecreaseZdecrease)};
      // if the intersect is outside
      return {
          (choiceIndicator || zOfIntersect > 0.)
              ? m_boundaryAccessors.tubeAccessor(Trk::TubeZincreaseRdecrease)
              : m_boundaryAccessors.tubeAccessor(Trk::TubeZdecreaseRdecrease)};
    }
    // =================================================================================================
 
    // ======================= the inside/outside part remains
    // =========================================
    //  (a) outside cases
    if (posR < m_innerRadius && deltaR < 0.)
      return {m_boundaryAccessors.tubeAccessor(Trk::TubeOutsideRminRdecrease)};
    if (posR > m_outerRadius && deltaR > 0.)
      return {m_boundaryAccessors.tubeAccessor(Trk::TubeOutsideRmaxRincrease)};
    if (posZ < -m_halfZ && deltaZ < 0.)
      return {m_boundaryAccessors.tubeAccessor(Trk::TubeOutsideZminZdecrease)};
    if (posZ > m_halfZ && deltaZ > 0.)
      return {m_boundaryAccessors.tubeAccessor(Trk::TubeOutsideZmaxZincrease)};
    // (b) inside cases
    // the increase R case
    if (deltaR > 0.) {
      // solve the linear equation for the outer Radius - with floating point
      // exception check
      Trk::CylinderIntersector intersectRmax(
        posR, posZ, deltaR != 0 ? deltaZ / deltaR : 0, m_outerRadius);
      double zOfIntersect = intersectRmax.yOfX;
 
      if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect > 0.)
        return {m_boundaryAccessors.tubeAccessor(Trk::TubeRincreaseZincrease)};
      if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect < 0.)
        return {m_boundaryAccessors.tubeAccessor(Trk::TubeRincreaseZdecrease)};
      if (std::abs(zOfIntersect) > m_halfZ && zOfIntersect < 0.)
        return {m_boundaryAccessors.tubeAccessor(Trk::TubeZdecreaseRincrease)};
      if (std::abs(zOfIntersect) > m_halfZ && zOfIntersect > 0.)
        return {m_boundaryAccessors.tubeAccessor(Trk::TubeZincreaseRincrease)};
 
    } else {
      // solve the linear equation for the inner Radius with floating point
      // exceptin check
      Trk::CylinderIntersector intersectRmin(
        posR, posZ, deltaR != 0 ? deltaZ / deltaR : 0, m_innerRadius);
      double zOfIntersect = intersectRmin.yOfX;
 
      if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect > 0.)
        return {m_boundaryAccessors.tubeAccessor(Trk::TubeRdecreaseZincrease)};
      if (std::abs(zOfIntersect) <= m_halfZ && zOfIntersect < 0.)
        return {m_boundaryAccessors.tubeAccessor(Trk::TubeRdecreaseZdecrease)};
      if (std::abs(zOfIntersect) > m_halfZ && zOfIntersect > 0.)
        return {m_boundaryAccessors.tubeAccessor(Trk::TubeZincreaseRdecrease)};
      return {m_boundaryAccessors.tubeAccessor(Trk::TubeZdecreaseRdecrease)};
    }
  }
  // the cylinder case
  if (
    m_innerRadius < s_numericalStable &&
    std::abs(m_halfPhiSector - M_PI) < s_numericalStable) {
 
    // ----------------------------- z / r information
    double posZ = gp.z();
    double posR = gp.perp();
 
    // is on R case
    bool isOnCylinder = (std::abs(posR - m_outerRadius) < s_numericalStable);
 
    // difference Vector
    Amg::Vector3D diff(gp + dir);
    // differences
    double deltaZ = diff.z() - posZ;
    double deltaR = diff.perp() - posR;
    // take the sign of the intersect
    int radiusSign = deltaR > 0. ? 1 : -1;
 
    // solve the linear equation for the cylinder Radius / with floating point
    // exceptin check
    Trk::CylinderIntersector intersectR(
      posR,
      posZ,
      deltaR != 0 ? deltaZ / deltaR : 0,
      radiusSign * m_outerRadius);
 
    //  the intersection point - a little trick for the OnFace cases
    double zOfIntersect = intersectR.yOfX;
    // check if intersection exists
    bool intersectsCylinder =
      !isOnCylinder ? (zOfIntersect * zOfIntersect <= m_halfZ * m_halfZ)
                    : false;
 
    // return the cases for going through the cylinder
    if (intersectsCylinder && zOfIntersect > 0.)
      return {m_boundaryAccessors.cylinderAccessor(Trk::CylinderZincrease)};
    if (intersectsCylinder && zOfIntersect <= 0.)
      return {m_boundaryAccessors.cylinderAccessor(Trk::CylinderZdecrease)};
    if (deltaZ > 0.)
      return {m_boundaryAccessors.cylinderAccessor(Trk::CylinderPositiveFace)};
    return {m_boundaryAccessors.cylinderAccessor(Trk::CylinderNegativeFace)};
  }
  // the sectoral cylinder case
  if (m_innerRadius != 0. && std::abs(m_halfPhiSector - M_PI) > 10e-3)
    return {m_boundaryAccessors.sectoralCylinderAccessor(
        Trk::StandardSectoralCylinder)};
  // it remains the sectoral tube case
  return {m_boundaryAccessors.sectoralTubeAccessor(Trk::StandardSectoralTube)};
}

clone()

CylinderVolumeBounds * Trk::CylinderVolumeBounds::clone ( ) const

inlineoverridevirtual

Virtual constructor.

Implements Trk::VolumeBounds.

Definition at line 171 of file CylinderVolumeBounds.h.

                                                               {
  return new CylinderVolumeBounds(*this);
}

createBoundarySurfaceAccessors()

void Trk::CylinderVolumeBounds::createBoundarySurfaceAccessors ( )

private

Private method to construct the accessors.

decomposeToSurfaces()

std::vector< std::unique_ptr< Trk::Surface > > Trk::CylinderVolumeBounds::decomposeToSurfaces ( const Amg::Transform3D &  transform )

finaloverridevirtual

Method to decompose the Bounds into boundarySurfaces.

Implements Trk::VolumeBounds.

Definition at line 98 of file CylinderVolumeBounds.cxx.

{
  auto retsf = std::vector<std::unique_ptr<Trk::Surface>>();
  // memory optimisation --- reserve the maximum
  retsf.reserve(6);
 
  // check if the transform is approximatively the identity
  bool isConcentric = transform.isApprox(Amg::Transform3D::Identity());
 
  Amg::RotationMatrix3D discRot = transform.linear();
  Amg::Vector3D cylCenter(transform.translation());
 
  // bottom Disc (negative z)
  Amg::RotationMatrix3D bottomDiscRot;
  bottomDiscRot.col(0) = discRot.col(1);
  bottomDiscRot.col(1) = discRot.col(0);
  bottomDiscRot.col(2) = -discRot.col(2);
  retsf.push_back(std::make_unique<Trk::DiscSurface>(
    Amg::Transform3D(
      transform * Amg::AngleAxis3D(M_PI, Amg::Vector3D(1., 0., 0.)) *
      Amg::Translation3D(Amg::Vector3D(0., 0., halflengthZ()))),
    bottomDiscBounds()));
  // top Disc (positive z)
  retsf.push_back(std::make_unique<Trk::DiscSurface>(
    Amg::Transform3D(transform* Amg::Translation3D( Amg::Vector3D(0, 0, halflengthZ()) )),
    topDiscBounds()));
  // outer Cylinder
  if (!isConcentric)
    retsf.push_back(std::make_unique<Trk::CylinderSurface>(
      Amg::Transform3D(transform), outerCylinderBounds()));
  else
    retsf.push_back(std::make_unique<Trk::CylinderSurface>(outerCylinderBounds()));
 
  // innermost Cylinder
  if (innerRadius() > s_numericalStable) {
    if (!isConcentric)
      retsf.push_back(std::make_unique<Trk::CylinderSurface>(
        Amg::Transform3D(transform), innerCylinderBounds()));
    else
      retsf.push_back(std::make_unique<Trk::CylinderSurface>(innerCylinderBounds()));
  }
 
  if (std::abs(halfPhiSector() - M_PI) > s_numericalStable) {
    // sectorPlane 1 (negative phi)
    retsf.push_back(std::make_unique<Trk::PlaneSurface>(
      Amg::Transform3D(
        transform *
        Amg::AngleAxis3D(-halfPhiSector(), Amg::Vector3D(0., 0., 1.)) *
        Amg::Translation3D(Amg::Vector3D(mediumRadius(), 0., 0.)) *
        Amg::AngleAxis3D(M_PI / 2, Amg::Vector3D(1., 0., 0.))),
      sectorPlaneBounds()));
    // sectorPlane 2 (positive phi)
    retsf.push_back(std::make_unique<Trk::PlaneSurface>(
      Amg::Transform3D(
        transform *
        Amg::AngleAxis3D(halfPhiSector(), Amg::Vector3D(0., 0., 1.)) *
        Amg::Translation3D(Amg::Vector3D(mediumRadius(), 0., 0.)) *
        Amg::AngleAxis3D(-M_PI / 2, Amg::Vector3D(1., 0., 0.))),
      sectorPlaneBounds()));
  }
  return retsf;
}

deltaRadius()

double Trk::CylinderVolumeBounds::deltaRadius ( ) const

inline

This method returns the delta radius.

Definition at line 199 of file CylinderVolumeBounds.h.

                                                      {
  return (m_outerRadius - m_innerRadius);
}

dump() [1/2]

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

finaloverridevirtual

Output Method for MsgStream.

Implements Trk::VolumeBounds.

Definition at line 379 of file CylinderVolumeBounds.cxx.

{
  std::stringstream temp_sl;
  temp_sl << std::setiosflags(std::ios::fixed);
  temp_sl << std::setprecision(2);
  temp_sl << "Trk::CylinderVolumeBounds: (rMin, rMax, halfPhi, halfZ) = ";
  temp_sl << m_innerRadius << ", " << m_outerRadius << ", " << m_halfPhiSector
          << ", " << m_halfZ;
  sl << temp_sl.str();
  return sl;
}

dump() [2/2]

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

finaloverridevirtual

Output Method for std::ostream.

Implements Trk::VolumeBounds.

Definition at line 392 of file CylinderVolumeBounds.cxx.

{
  std::stringstream temp_sl;
  temp_sl << std::setiosflags(std::ios::fixed);
  temp_sl << std::setprecision(2);
  temp_sl << "Trk::CylinderVolumeBounds: (rMin, rMax, halfPhi, halfZ) = ";
  temp_sl << m_innerRadius << ", " << m_outerRadius << ", " << m_halfPhiSector
          << ", " << m_halfZ;
  sl << temp_sl.str();
  return sl;
}

halflengthZ()

double Trk::CylinderVolumeBounds::halflengthZ ( ) const

inline

This method returns the halflengthZ.

Definition at line 207 of file CylinderVolumeBounds.h.

{ return m_halfZ; }

halfPhiSector()

double Trk::CylinderVolumeBounds::halfPhiSector ( ) const

inline

This method returns the halfPhiSector angle.

Definition at line 203 of file CylinderVolumeBounds.h.

                                                        {
  return m_halfPhiSector;
}

innerCylinderBounds()

std::shared_ptr< Trk::CylinderBounds > Trk::CylinderVolumeBounds::innerCylinderBounds ( ) const

private

This method returns the associated CylinderBounds of the inner CylinderSurfaces.

Definition at line 352 of file CylinderVolumeBounds.cxx.

{
  return std::make_shared<Trk::CylinderBounds>(m_innerRadius, m_halfPhiSector, m_halfZ);
}

innerRadius()

double Trk::CylinderVolumeBounds::innerRadius ( ) const

inline

This method returns the inner radius.

Definition at line 187 of file CylinderVolumeBounds.h.

                                                      {
  return m_innerRadius;
}

inside()

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

inlinefinaloverridevirtual

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

Implements Trk::VolumeBounds.

Definition at line 175 of file CylinderVolumeBounds.h.

                                                           {
  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);
}

mediumRadius()

double Trk::CylinderVolumeBounds::mediumRadius ( ) const

inline

This method returns the medium radius.

Definition at line 195 of file CylinderVolumeBounds.h.

                                                       {
  return 0.5 * (m_innerRadius + m_outerRadius);
}

operator=()

Trk::CylinderVolumeBounds & Trk::CylinderVolumeBounds::operator=

(

const CylinderVolumeBounds &

cylbo

)

Assignment operator.

Definition at line 86 of file CylinderVolumeBounds.cxx.

{
  if (this != &cylbo) {
    m_innerRadius = cylbo.m_innerRadius;
    m_outerRadius = cylbo.m_outerRadius;
    m_halfPhiSector = cylbo.m_halfPhiSector;
    m_halfZ = cylbo.m_halfZ;
  }
  return *this;
}

outerCylinderBounds()

std::shared_ptr< Trk::CylinderBounds > Trk::CylinderVolumeBounds::outerCylinderBounds ( ) const

private

This method returns the associated CylinderBounds of the outer CylinderSurfaces.

Definition at line 358 of file CylinderVolumeBounds.cxx.

{
  return std::make_shared<Trk::CylinderBounds>(m_outerRadius, m_halfPhiSector, m_halfZ);
}

outerRadius()

double Trk::CylinderVolumeBounds::outerRadius ( ) const

inline

This method returns the outer radius.

Definition at line 191 of file CylinderVolumeBounds.h.

                                                      {
  return m_outerRadius;
}

sectorPlaneBounds()

std::shared_ptr< Trk::RectangleBounds > Trk::CylinderVolumeBounds::sectorPlaneBounds ( ) const

private

This method returns the associated PlaneBounds limiting a sectoral CylinderVolume.

Definition at line 370 of file CylinderVolumeBounds.cxx.

{
  return std::make_shared<Trk::RectangleBounds>(
    0.5 * (m_outerRadius - m_innerRadius), m_halfZ);
}

topDiscBounds()

std::shared_ptr< DiscBounds > Trk::CylinderVolumeBounds::topDiscBounds ( ) const

inlineprivate

This method returns the associated DiscBounds for the bottom/top DiscSurface.

Definition at line 209 of file CylinderVolumeBounds.h.

                                                                           {
  return this->bottomDiscBounds();
}

Member Data Documentation

m_boundaryAccessors

CylinderVolumeBoundaryAccessors Trk::CylinderVolumeBounds::m_boundaryAccessors

private

Accessors for Boundary surface access - static is not possible due to mismatched delete() / free () with TrkMagFieldUtils.

Definition at line 165 of file CylinderVolumeBounds.h.

m_halfPhiSector

double Trk::CylinderVolumeBounds::m_halfPhiSector

private

Definition at line 159 of file CylinderVolumeBounds.h.

m_halfZ

double Trk::CylinderVolumeBounds::m_halfZ

private

Definition at line 160 of file CylinderVolumeBounds.h.

m_innerRadius

double Trk::CylinderVolumeBounds::m_innerRadius

private

Definition at line 157 of file CylinderVolumeBounds.h.

m_outerRadius

double Trk::CylinderVolumeBounds::m_outerRadius

private

Definition at line 158 of file CylinderVolumeBounds.h.

s_numericalStable

const double Trk::CylinderVolumeBounds::s_numericalStable

staticprivate

Initial value:

=
  10e-2 * Gaudi::Units::mm

numerical stability

Definition at line 168 of file CylinderVolumeBounds.h.

---

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

• CylinderVolumeBounds.h
• CylinderVolumeBounds.cxx