geo::BoxBoundedGeo Class Reference

A base class aware of world box coordinatesAn object describing a simple shape can inherit from this one to gain access to a common boundary checking interface. More...

#include <BoxBoundedGeo.h>

Inheritance diagram for geo::BoxBoundedGeo:

Public Types
using	Coords_t = geo::Point_t
	Type of the coordinate triplet. More...
using	Coord_t = Coords_t::Scalar
	Type of the coordinate. More...

Public Member Functions
	BoxBoundedGeo ()=default
	Default constructor: sets an empty volume. More...
	BoxBoundedGeo (Coord_t x_min, Coord_t x_max, Coord_t y_min, Coord_t y_max, Coord_t z_min, Coord_t z_max)
	Constructor: sets the boundaries in world coordinates as specified. More...
	BoxBoundedGeo (Coords_t lower, Coords_t upper)
	Constructor: sets the boundaries in world coordinates as specified. More...
Dimension queries
double	MinX () const
	Returns the world x coordinate of the start of the box. More...
double	MaxX () const
	Returns the world x coordinate of the end of the box. More...
double	CenterX () const
	Returns the world x coordinate of the center of the box. More...
double	SizeX () const
	Returns the full size in the X dimension. More...
double	HalfSizeX () const
	Returns the size from the center to the border on X dimension. More...
double	MinY () const
	Returns the world y coordinate of the start of the box. More...
double	MaxY () const
	Returns the world y coordinate of the end of the box. More...
double	CenterY () const
	Returns the world y coordinate of the center of the box. More...
double	SizeY () const
	Returns the full size in the Y dimension. More...
double	HalfSizeY () const
	Returns the size from the center to the border on Y dimension. More...
double	MinZ () const
	Returns the world z coordinate of the start of the box. More...
double	MaxZ () const
	Returns the world z coordinate of the end of the box. More...
double	CenterZ () const
	Returns the world z coordinate of the center of the box. More...
double	SizeZ () const
	Returns the full size in the Z dimension. More...
double	HalfSizeZ () const
	Returns the size from the center to the border on Z dimension. More...
geo::Point_t	Min () const
	Returns the corner point with the smallest coordinates. More...
geo::Point_t	Max () const
	Returns the corner point with the largest coordinates. More...
geo::Point_t	Center () const
	Returns the center point of the box. More...
Containment in the full volume
bool	ContainsX (double x, double const wiggle=1) const
	Returns whether this TPC contains the specified world x coordinate. More...
bool	ContainsY (double y, double const wiggle=1) const
	Returns whether this TPC contains the specified world y coordinate. More...
bool	ContainsZ (double z, double const wiggle=1) const
	Returns whether this TPC contains the specified world z coordinate. More...
bool	ContainsYZ (double y, double z, double const wiggle=1) const
	Returns if TPC contains the specified world y and z coordinates. More...
bool	ContainsPosition (geo::Point_t const &point, double wiggle=1.0) const
	Returns whether this volume contains the specified point. More...
bool	ContainsPosition (TVector3 const &point, double wiggle=1.0) const
bool	ContainsPosition (double const *point, double wiggle=1.0) const
Containment in a fiducial volume
bool	InFiducialX (double x, double neg_margin, double pos_margin) const
	Returns whether TPC fiducial volume contains world x coordinate. More...
bool	InFiducialX (double x, double margin) const
	Returns whether TPC fiducial volume contains world x coordinate. More...
bool	InFiducialY (double y, double neg_margin, double pos_margin) const
	Returns whether TPC fiducial volume contains world y coordinate. More...
bool	InFiducialY (double y, double margin) const
	Returns whether TPC fiducial volume contains world y coordinate. More...
bool	InFiducialZ (double z, double neg_margin, double pos_margin) const
	Returns whether TPC fiducial volume contains world z coordinate. More...
bool	InFiducialZ (double z, double margin) const
	Returns whether TPC fiducial volume contains world z coordinate. More...
Overlaps
bool	OverlapsX (geo::BoxBoundedGeo const &other) const
	Returns if the x coordinates covered by this and other box overlap. More...
bool	OverlapsY (geo::BoxBoundedGeo const &other) const
	Returns if the y coordinates covered by this and other box overlap. More...
bool	OverlapsZ (geo::BoxBoundedGeo const &other) const
	Returns if the z coordinates covered by this and other box overlap. More...
bool	Overlaps (geo::BoxBoundedGeo const &other) const
	Returns if this and other box overlap. More...
Setting dimensions
void	SetBoundaries (Coord_t x_min, Coord_t x_max, Coord_t y_min, Coord_t y_max, Coord_t z_min, Coord_t z_max)
	Sets the boundaries in world coordinates as specified. More...
void	SetBoundaries (Coords_t lower, Coords_t upper)
	Sets the boundaries in world coordinates as specified. More...
void	ExtendToInclude (Coord_t x, Coord_t y, Coord_t z)
	Extends the current box to also include the specified point. More...
void	ExtendToInclude (geo::Point_t const &point)
	Extends the current box to also include the specified point. More...
void	ExtendToInclude (BoxBoundedGeo const &box)
	Extends the current box to also include the specified one. More...
std::vector< TVector3 >	GetIntersections (TVector3 const &TrajectoryStart, TVector3 const &TrajectoryDirect) const
	Calculates the entry and exit points of a trajectory on the box surface. More...
std::vector< geo::Point_t >	GetIntersections (geo::Point_t const &TrajectoryStart, geo::Vector_t const &TrajectoryDirect) const

Static Public Member Functions
static bool	CoordinateContained (double c, double min, double max, double wiggle=1.)
	Returns whether the specified coordinate is in a range. More...
static bool	CoordinateContained (double c, double const *range, double wiggle=1.)
	Returns whether the specified coordinate is in a range. More...
static void	set_min (Coord_t &var, Coord_t value)
	Sets var to value if value is smaller than the current var value. More...
static void	set_max (Coord_t &var, Coord_t value)
	Sets var to value if value is larger than the current var value. More...
static void	set_min (Coords_t &var, geo::Point_t const &value)
	Sets each coordinate of var to the one in value if the latter is smaller. More...
static void	set_max (Coords_t &var, geo::Point_t const &value)
	Sets each coordinate of var to the one in value if the latter is larger. More...

Private Member Functions
void	SortCoordinates ()
	Makes sure each coordinate of the minimum point is smaller than maximum. More...

Private Attributes
Coords_t	c_min
	minimum coordinates (x, y, z) More...
Coords_t	c_max
	maximum coordinates (x, y, z) More...

Detailed Description

A base class aware of world box coordinates

An object describing a simple shape can inherit from this one to gain access to a common boundary checking interface.

The boundary is a simple box with axes parallel to the coordinate system.

Definition at line 33 of file BoxBoundedGeo.h.

Member Typedef Documentation

using geo::BoxBoundedGeo::Coord_t = Coords_t::Scalar

Type of the coordinate.

Definition at line 36 of file BoxBoundedGeo.h.

using geo::BoxBoundedGeo::Coords_t = geo::Point_t

Type of the coordinate triplet.

Definition at line 35 of file BoxBoundedGeo.h.

Constructor & Destructor Documentation

geo::BoxBoundedGeo::BoxBoundedGeo ( )

default

Default constructor: sets an empty volume.

See Also

SetBoundaries

We allow for a default construction since the derived object might need some specific action before being aware of its boundaries. In that case, SetBoundaries() will set the boundaries.

geo::BoxBoundedGeo::BoxBoundedGeo ( Coord_t  x_min, Coord_t  x_max, Coord_t  y_min, Coord_t  y_max, Coord_t  z_min, Coord_t  z_max  )

inline

Constructor: sets the boundaries in world coordinates as specified.

Parameters

x_min	lower x coordinate
x_max	upper x coordinate
y_min	lower y coordinate
y_max	upper y coordinate
z_min	lower z coordinate
z_max	upper z coordinate

See Also

SetBoundaries()

Note that is assumed that each minimum is larger than its maximum, and no check is performed.

Definition at line 62 of file BoxBoundedGeo.h.

       :
      c_min{ x_min, y_min, z_min }, c_max{ x_max, y_max, z_max }
      { SortCoordinates(); }

geo::BoxBoundedGeo::BoxBoundedGeo ( Coords_t  lower, Coords_t  upper  )

inline

Constructor: sets the boundaries in world coordinates as specified.

Parameters

lower	lower coordinates (x, y, z)
upper	upper coordinates (x, y, z)

See Also

SetBoundaries()

Note that is assumed that each minimum is larger than its maximum, and no check is performed.

Definition at line 79 of file BoxBoundedGeo.h.

                                                 :
      c_min(lower), c_max(upper)
      { SortCoordinates(); }

Member Function Documentation

geo::Point_t geo::BoxBoundedGeo::Center ( ) const

inline

Returns the center point of the box.

Definition at line 139 of file BoxBoundedGeo.h.

      { return { CenterX(), CenterY(), CenterZ() }; }

double geo::BoxBoundedGeo::CenterX ( ) const

inline

Returns the world x coordinate of the center of the box.

Definition at line 94 of file BoxBoundedGeo.h.

{ return (MinX() + MaxX()) / 2.; }

double geo::BoxBoundedGeo::CenterY ( ) const

inline

Returns the world y coordinate of the center of the box.

Definition at line 109 of file BoxBoundedGeo.h.

{ return (MinY() + MaxY()) / 2.; }

double geo::BoxBoundedGeo::CenterZ ( ) const

inline

Returns the world z coordinate of the center of the box.

Definition at line 124 of file BoxBoundedGeo.h.

{ return (MinZ() + MaxZ()) / 2.; }

bool geo::BoxBoundedGeo::ContainsPosition ( geo::Point_t const &  point, double  wiggle = 1.0  ) const

inline

Returns whether this volume contains the specified point.

Parameters

point	the point [cm]
wiggle	expansion factor for the range

Returns

whether the specified coordinate is in this volume

If the wiggle is larger than 1, each size of the volume is expanded by the wiggle factor. If the wiggle is less than 1, each size is shrinked.

Definition at line 205 of file BoxBoundedGeo.h.

      {
        return ContainsX(point.X(), wiggle)
          && ContainsYZ(point.Y(), point.Z(), wiggle);
      } // ContainsPosition()

bool geo::BoxBoundedGeo::ContainsPosition	(	TVector3 const &	point,
		double	wiggle = 1.0
	)		const

See Also

ContainsPosition(geo::Point_t const&, double) const.

Definition at line 26 of file BoxBoundedGeo.cxx.

    { return ContainsPosition(geo::vect::toPoint(point), wiggle); }

bool geo::BoxBoundedGeo::ContainsPosition	(	double const *	point,
		double	wiggle = 1.0
	)		const

See Also

ContainsPosition(geo::Point_t const&, double) const.

Definition at line 31 of file BoxBoundedGeo.cxx.

    { return ContainsPosition(geo::vect::makePointFromCoords(point), wiggle); }

bool geo::BoxBoundedGeo::ContainsX ( double  x, double const  wiggle = 1  ) const

inline

Returns whether this TPC contains the specified world x coordinate.

Parameters

x	the absolute ("world") coordinate x
wiggle	expansion factor for the range (see ContainsPosition())

Returns

whether the specified coordinate is in this TPC

See Also

ContainsPosition()

Note that x is by definition the drift direction, and a reconstructed x typically depends on an assumption respect to the event time.

Definition at line 159 of file BoxBoundedGeo.h.

      { return CoordinateContained(x, MinX(), MaxX(), wiggle); }

bool geo::BoxBoundedGeo::ContainsY ( double  y, double const  wiggle = 1  ) const

inline

Returns whether this TPC contains the specified world y coordinate.

Parameters

y	the absolute ("world") coordinate y
wiggle	expansion factor for the range (see ContainsPosition())

Returns

whether the specified coordinate is in this TPC

See Also

ContainsPosition()

Definition at line 169 of file BoxBoundedGeo.h.

      { return CoordinateContained(y, MinY(), MaxY(), wiggle); }

bool geo::BoxBoundedGeo::ContainsYZ ( double  y, double  z, double const  wiggle = 1  ) const

inline

Returns if TPC contains the specified world y and z coordinates.

Parameters

y	the absolute ("world") coordinate y
z	the absolute ("world") coordinate z
wiggle	expansion factor for the range (see ContainsPosition())

Returns

whether the specified coordinate is in this TPC

See Also

ContainsPosition()

Definition at line 190 of file BoxBoundedGeo.h.

      { return ContainsY(y, wiggle) && ContainsZ(z, wiggle); }

bool geo::BoxBoundedGeo::ContainsZ ( double  z, double const  wiggle = 1  ) const

inline

Returns whether this TPC contains the specified world z coordinate.

Parameters

z	the absolute ("world") coordinate z
wiggle	expansion factor for the range (see ContainsPosition())

Returns

whether the specified coordinate is in this TPC

See Also

ContainsPosition()

Definition at line 179 of file BoxBoundedGeo.h.

      { return CoordinateContained(z, MinZ(), MaxZ(), wiggle); }

static bool geo::BoxBoundedGeo::CoordinateContained ( double  c, double  min, double  max, double  wiggle = 1.  )

inlinestatic

Returns whether the specified coordinate is in a range.

Parameters

c	the coordinate
min	lower boundary of the range
max	upper boundary of the range
wiggle	expansion factor for the range

Returns

whether the specified coordinate is in a range

If the wiggle is larger than 1, the range is expanded by the wiggle factor. If the wiggle is less than 1, the range is shrinked.

Definition at line 355 of file BoxBoundedGeo.h.

      {
        return (c >= (min > 0? min / wiggle: min * wiggle))
          && (c <= (max < 0? max / wiggle: max * wiggle));
      } // CoordinateContained()

static bool geo::BoxBoundedGeo::CoordinateContained ( double  c, double const *  range, double  wiggle = 1.  )

inlinestatic

Returns whether the specified coordinate is in a range.

Parameters

c	the coordinate
range	pointer to [ lower boundary, upper boundary ] of the range
wiggle	expansion factor for the range

Returns

whether the specified coordinate is in a range

See Also

CoordinateContained(double, double, double, double)

If the wiggle is larger than 1, the range is expanded by the wiggle factor. If the wiggle is less than 1, the range is shrinked.

Definition at line 373 of file BoxBoundedGeo.h.

      { return CoordinateContained(c, range[0], range[1], wiggle); }

void geo::BoxBoundedGeo::ExtendToInclude ( Coord_t  x, Coord_t  y, Coord_t  z  )

inline

Extends the current box to also include the specified point.

Parameters

x	x coordinate of the point to include
y	y coordinate of the point to include
z	z coordinate of the point to include

Definition at line 414 of file BoxBoundedGeo.h.

      { ExtendToInclude(geo::Point_t(x, y, z)); }

void geo::BoxBoundedGeo::ExtendToInclude ( geo::Point_t const &  point )

inline

Extends the current box to also include the specified point.

Parameters

point

coordinates of the point to include

Definition at line 421 of file BoxBoundedGeo.h.

      {
        set_min(c_min, point);
        set_max(c_max, point);
      }

void geo::BoxBoundedGeo::ExtendToInclude ( BoxBoundedGeo const &  box )

inline

Extends the current box to also include the specified one.

Parameters

box	the box to include

It is assumed that the box has its boundaries properly sorted.

Definition at line 433 of file BoxBoundedGeo.h.

      {
        set_min(c_min, box.Min());
        set_max(c_max, box.Max());
      } // ExtendToInclude()

std::vector< TVector3 > geo::BoxBoundedGeo::GetIntersections	(	TVector3 const &	TrajectoryStart,
		TVector3 const &	TrajectoryDirect
	)		const

Calculates the entry and exit points of a trajectory on the box surface.

Author

Christoph Rudolf von Rohr (crohr.nosp@m.@fna.nosp@m.l.gov)

Date

July 27th, 2015

Parameters

TrajectoryStart	position of the trajectory source
TrajectoryDirect	direction vector of the trajectory

This member is public since it just gives an output and does not change any member. The algorithm works only for a box shaped active volume with facing walls parallel to axis. If the return std::vector is empty the trajectory does not intersect with the box. Normally the return value should have one (if the trajectory originates in the box) or two (else) entries. If the return value has two entries the first represents the entry point and the second the exit point

Definition at line 126 of file BoxBoundedGeo.cxx.

  {
    std::vector<TVector3> intersections;
    for (auto const& point: GetIntersections(geo::Point_t(TrajectoryStart), geo::Vector_t(TrajectoryDirect)))
      intersections.emplace_back(point.X(), point.Y(), point.Z());
    return intersections;
  } // GetIntersections(TVector3)

std::vector< geo::Point_t > geo::BoxBoundedGeo::GetIntersections	(	geo::Point_t const &	TrajectoryStart,
		geo::Vector_t const &	TrajectoryDirect
	)		const

Definition at line 35 of file BoxBoundedGeo.cxx.

          {

    std::vector<geo::Point_t> IntersectionPoints;
    std::vector<double> LineParameters;

    // Generate normal vectors and offsets for every plane of the box
    // All normal vectors are headed outwards
    // BUG the double brace syntax is required to work around clang bug 21629
    // (https://bugs.llvm.org/show_bug.cgi?id=21629)
    static std::array<geo::Vector_t, 6U> const NormalVectors = {{
      -geo::Xaxis<geo::Vector_t>(), geo::Xaxis<geo::Vector_t>(), // anode, cathode,
      -geo::Yaxis<geo::Vector_t>(), geo::Yaxis<geo::Vector_t>(), // bottom, top,
      -geo::Zaxis<geo::Vector_t>(), geo::Zaxis<geo::Vector_t>()  // upstream, downstream
    }};
    // BUG the double brace syntax is required to work around clang bug 21629
    // (https://bugs.llvm.org/show_bug.cgi?id=21629)
    std::array<geo::Point_t, 6U> const NormalVectorOffsets = {{
      geo::Point_t{ Min().X(), Min().Y(), Min().Z() }, // Anode offset
      geo::Point_t{ Max().X(), Min().Y(), Min().Z() }, // Cathode
      geo::Point_t{ Min().X(), Min().Y(), Min().Z() }, // Bottom
      geo::Point_t{ Min().X(), Max().Y(), Min().Z() }, // Top
      geo::Point_t{ Min().X(), Min().Y(), Min().Z() }, // upstream
      geo::Point_t{ Min().X(), Min().Y(), Max().Z() }  // downstream
    }};

    // Loop over all surfaces of the box
    for(unsigned int face_no = 0; face_no < NormalVectors.size(); face_no++)
    {
      // Check if trajectory and surface are not parallel
      if(NormalVectors[face_no].Dot(TrajectoryDirect))
      {
        // Calculate the line parameter for the intersection points
        LineParameters.push_back( NormalVectors[face_no].Dot(NormalVectorOffsets.at(face_no) - TrajectoryStart)
                                / NormalVectors[face_no].Dot(TrajectoryDirect) );
      }
      else continue;

      // Calculate intersection point using the line parameter
      IntersectionPoints.push_back( TrajectoryStart + LineParameters.back()*TrajectoryDirect );

      // Coordinate which should be ignored when checking for limits added by Christoph Rudolf von Rohr 05/21/2016
      unsigned int NoCheckCoord;

      // Calculate NoCheckCoord out of the face_no
      if(face_no % 2)
      {
          // Convert odd face number to coordinate
          NoCheckCoord = (face_no - 1)/2;
      }
      else
      {
          // Convert even face number to coordinate
          NoCheckCoord = face_no/2;
      }

      // Loop over all three space coordinates
      unsigned int coord = 0;
      for(auto extractCoord: geo::vect::coordReaders<geo::Point_t>())
      {
        auto const lastPointCoord = geo::vect::bindCoord(IntersectionPoints.back(), extractCoord);
        auto const minCoord = geo::vect::bindCoord(c_min, extractCoord);
        auto const maxCoord = geo::vect::bindCoord(c_max, extractCoord);

        // Changed by Christoph Rudolf von Rohr 05/21/2016
        // Then check if point is not within the surface limits at this coordinate, without looking
        // at the plane normal vector coordinate. We can assume, that our algorithm already found this coordinate correctily.
        // In rare cases, looking for boundaries in this coordinate causes unexpected behavior due to floating point inaccuracies.
        if( coord++ != NoCheckCoord && ((lastPointCoord() > maxCoord()) || (lastPointCoord() < minCoord)) )
        {
          // if off limits, get rid of the useless data and break the coordinate loop
          LineParameters.pop_back();
          IntersectionPoints.pop_back();
          break;
        }
      }// coordinate loop
    }// Surcaces loop

    // sort points according to their parameter value (first is entry, second is exit)
    if(LineParameters.size() == 2 && LineParameters.front() > LineParameters.back())
    {
      std::swap(IntersectionPoints.front(),IntersectionPoints.back());
    }

    return IntersectionPoints;
  } // GetIntersections()

double geo::BoxBoundedGeo::HalfSizeX ( ) const

inline

Returns the size from the center to the border on X dimension.

Definition at line 100 of file BoxBoundedGeo.h.

{ return SizeX() / 2.0; }

double geo::BoxBoundedGeo::HalfSizeY ( ) const

inline

Returns the size from the center to the border on Y dimension.

Definition at line 115 of file BoxBoundedGeo.h.

{ return SizeY() / 2.0; }

double geo::BoxBoundedGeo::HalfSizeZ ( ) const

inline

Returns the size from the center to the border on Z dimension.

Definition at line 130 of file BoxBoundedGeo.h.

{ return SizeZ() / 2.0; }

bool geo::BoxBoundedGeo::InFiducialX ( double  x, double  neg_margin, double  pos_margin  ) const

inline

Returns whether TPC fiducial volume contains world x coordinate.

Parameters

x	the absolute ("world") coordinate x
neg_margin	how far within the TPC the fiducial region starts
pos_margin	how far before the TPC the fiducial region ends

Returns

whether the specified coordinate is in this TPC fiducial volume

The fiducial volume is defined by the specified margins, that denote how much of the coordinate is not fiducial. For example,

bool bWithin = tpc->InFiducialX(x, 5., 8.);

on a TPC with x from 0 to 100 (cm, presumably) will return true for all x between 5 and 92. If positive margin is not specified, it's assumed to be the same as the negative one. Specifying a negative mergin effectively extends the TPC volume. Note that x is by definition the drift direction, and a reconstructed x typically depends on an assumption respect to the event time.

Definition at line 241 of file BoxBoundedGeo.h.

      {
        return CoordinateContained(x, MinX() + neg_margin, MaxX() - pos_margin);
      }

bool geo::BoxBoundedGeo::InFiducialX ( double  x, double  margin  ) const

inline

Returns whether TPC fiducial volume contains world x coordinate.

See Also

InFiducialX(double, double, double) const

Margins are symmetric.

Definition at line 251 of file BoxBoundedGeo.h.

      { return InFiducialX(x, margin, margin); }

bool geo::BoxBoundedGeo::InFiducialY ( double  y, double  neg_margin, double  pos_margin  ) const

inline

Returns whether TPC fiducial volume contains world y coordinate.

Parameters

y	the absolute ("world") coordinate y
neg_margin	how far within the TPC the fiducial region starts
pos_margin	how far before the TPC the fiducial region ends

Returns

whether the specified coordinate is in this TPC fiducial volume

The fiducial volume is defined by the specified margins, that denote how much of the coordinate is not fiducial. For example,

bool bWithin = tpc->InFiducialY(y, 5., 8.);

on a TPC with y from 0 to 100 (cm, presumably) will return true for all y between 5 and 92. If positive margin is not specified, it's assumed to be the same as the negative one. Specifying a negative mergin effectively extends the TPC volume.

Definition at line 272 of file BoxBoundedGeo.h.

      {
        return CoordinateContained(y, MinY() + neg_margin, MaxY() - pos_margin);
      }

bool geo::BoxBoundedGeo::InFiducialY ( double  y, double  margin  ) const

inline

Returns whether TPC fiducial volume contains world y coordinate.

See Also

InFiducialY(double, double, double) const

Margins are symmetric.

Definition at line 282 of file BoxBoundedGeo.h.

      { return InFiducialY(y, margin, margin); }

bool geo::BoxBoundedGeo::InFiducialZ ( double  z, double  neg_margin, double  pos_margin  ) const

inline

Returns whether TPC fiducial volume contains world z coordinate.

Parameters

z	the absolute ("world") coordinate y
neg_margin	how far within the TPC the fiducial region starts
pos_margin	how far before the TPC the fiducial region ends

Returns

whether the specified coordinate is in this TPC fiducial volume

The fiducial volume is defined by the specified margins, that denote how much of the coordinate is not fiducial. For example,

bool bWithin = tpc->InFiducialZ(z, 5., 8.);

on a TPC with z from 0 to 100 (cm, presumably) will return true for all z between 5 and 92. If positive margin is not specified, it's assumed to be the same as the negative one. Specifying a negative mergin effectively extends the TPC volume.

Definition at line 303 of file BoxBoundedGeo.h.

      {
        return CoordinateContained(z, MinZ() + neg_margin, MaxZ() - pos_margin);
      }

bool geo::BoxBoundedGeo::InFiducialZ ( double  z, double  margin  ) const

inline

Returns whether TPC fiducial volume contains world z coordinate.

See Also

InFiducialZ(double, double, double) const

Margins are symmetric.

Definition at line 313 of file BoxBoundedGeo.h.

      { return InFiducialZ(z, margin, margin); }

geo::Point_t geo::BoxBoundedGeo::Max ( ) const

inline

Returns the corner point with the largest coordinates.

Definition at line 136 of file BoxBoundedGeo.h.

{ return c_max; }

double geo::BoxBoundedGeo::MaxX ( ) const

inline

Returns the world x coordinate of the end of the box.

Definition at line 91 of file BoxBoundedGeo.h.

{ return c_max.X(); }

double geo::BoxBoundedGeo::MaxY ( ) const

inline

Returns the world y coordinate of the end of the box.

Definition at line 106 of file BoxBoundedGeo.h.

{ return c_max.Y(); }

double geo::BoxBoundedGeo::MaxZ ( ) const

inline

Returns the world z coordinate of the end of the box.

Definition at line 121 of file BoxBoundedGeo.h.

{ return c_max.Z(); }

geo::Point_t geo::BoxBoundedGeo::Min ( ) const

inline

Returns the corner point with the smallest coordinates.

Definition at line 133 of file BoxBoundedGeo.h.

{ return c_min; }

double geo::BoxBoundedGeo::MinX ( ) const

inline

Returns the world x coordinate of the start of the box.

Definition at line 88 of file BoxBoundedGeo.h.

{ return c_min.X(); }

double geo::BoxBoundedGeo::MinY ( ) const

inline

Returns the world y coordinate of the start of the box.

Definition at line 103 of file BoxBoundedGeo.h.

{ return c_min.Y(); }

double geo::BoxBoundedGeo::MinZ ( ) const

inline

Returns the world z coordinate of the start of the box.

Definition at line 118 of file BoxBoundedGeo.h.

{ return c_min.Z(); }

bool geo::BoxBoundedGeo::Overlaps ( geo::BoxBoundedGeo const &  other ) const

inline

Returns if this and other box overlap.

Definition at line 336 of file BoxBoundedGeo.h.

      { return OverlapsX(other) && OverlapsY(other) && OverlapsZ(other); }

bool geo::BoxBoundedGeo::OverlapsX ( geo::BoxBoundedGeo const &  other ) const

inline

Returns if the x coordinates covered by this and other box overlap.

Definition at line 324 of file BoxBoundedGeo.h.

      { return std::min(MaxX(), other.MaxX()) > std::max(MinX(), other.MinX()); }

bool geo::BoxBoundedGeo::OverlapsY ( geo::BoxBoundedGeo const &  other ) const

inline

Returns if the y coordinates covered by this and other box overlap.

Definition at line 328 of file BoxBoundedGeo.h.

      { return std::min(MaxY(), other.MaxY()) > std::max(MinY(), other.MinY()); }

bool geo::BoxBoundedGeo::OverlapsZ ( geo::BoxBoundedGeo const &  other ) const

inline

Returns if the z coordinates covered by this and other box overlap.

Definition at line 332 of file BoxBoundedGeo.h.

      { return std::min(MaxZ(), other.MaxZ()) > std::max(MinZ(), other.MinZ()); }

static void geo::BoxBoundedGeo::set_max ( Coord_t &  var, Coord_t  value  )

inlinestatic

Sets var to value if value is larger than the current var value.

Definition at line 468 of file BoxBoundedGeo.h.

      { if (value > var) var = value; }

static void geo::BoxBoundedGeo::set_max ( Coords_t &  var, geo::Point_t const &  value  )

inlinestatic

Sets each coordinate of var to the one in value if the latter is larger.

Definition at line 480 of file BoxBoundedGeo.h.

      {
        if (value.X() > var.X()) var.SetX(value.X());
        if (value.Y() > var.Y()) var.SetY(value.Y());
        if (value.Z() > var.Z()) var.SetZ(value.Z());
      }

static void geo::BoxBoundedGeo::set_min ( Coord_t &  var, Coord_t  value  )

inlinestatic

Sets var to value if value is smaller than the current var value.

Definition at line 464 of file BoxBoundedGeo.h.

      { if (value < var) var = value; }

static void geo::BoxBoundedGeo::set_min ( Coords_t &  var, geo::Point_t const &  value  )

inlinestatic

Sets each coordinate of var to the one in value if the latter is smaller.

Definition at line 472 of file BoxBoundedGeo.h.

      {
        if (value.X() < var.X()) var.SetX(value.X());
        if (value.Y() < var.Y()) var.SetY(value.Y());
        if (value.Z() < var.Z()) var.SetZ(value.Z());
      }

void geo::BoxBoundedGeo::SetBoundaries ( Coord_t  x_min, Coord_t  x_max, Coord_t  y_min, Coord_t  y_max, Coord_t  z_min, Coord_t  z_max  )

inline

Sets the boundaries in world coordinates as specified.

Parameters

x_min	lower x coordinate
x_max	upper x coordinate
y_min	lower y coordinate
y_max	upper y coordinate
z_min	lower z coordinate
z_max	upper z coordinate

Definition at line 389 of file BoxBoundedGeo.h.

      {
        c_min.SetXYZ(x_min, y_min, z_min);
        c_max.SetXYZ(x_max, y_max, z_max);
        SortCoordinates();
      }

void geo::BoxBoundedGeo::SetBoundaries ( Coords_t  lower, Coords_t  upper  )

inline

Sets the boundaries in world coordinates as specified.

Parameters

lower	lower coordinates (x, y, z)
upper	upper coordinates (x, y, z)

Definition at line 405 of file BoxBoundedGeo.h.

      { c_min = lower; c_max = upper; SortCoordinates(); }

double geo::BoxBoundedGeo::SizeX ( ) const

inline

Returns the full size in the X dimension.

Definition at line 97 of file BoxBoundedGeo.h.

{ return MaxX() - MinX(); }

double geo::BoxBoundedGeo::SizeY ( ) const

inline

Returns the full size in the Y dimension.

Definition at line 112 of file BoxBoundedGeo.h.

{ return MaxY() - MinY(); }

double geo::BoxBoundedGeo::SizeZ ( ) const

inline

Returns the full size in the Z dimension.

Definition at line 127 of file BoxBoundedGeo.h.

{ return MaxZ() - MinZ(); }

void geo::BoxBoundedGeo::SortCoordinates ( )

private

Makes sure each coordinate of the minimum point is smaller than maximum.

Definition at line 135 of file BoxBoundedGeo.cxx.

                                      {

    for (auto coordMan: geo::vect::coordManagers<geo::Point_t>()) {
      auto min = geo::vect::bindCoord(c_min, coordMan);
      auto max = geo::vect::bindCoord(c_max, coordMan);
      if (min() > max()) {
        auto temp = min();
        min = max();
        max = temp;
      }
    } // for

  } // BoxBoundedGeo::SortCoordinates()

Member Data Documentation

Coords_t geo::BoxBoundedGeo::c_max

private

maximum coordinates (x, y, z)

Definition at line 491 of file BoxBoundedGeo.h.

Coords_t geo::BoxBoundedGeo::c_min

private

minimum coordinates (x, y, z)

Definition at line 490 of file BoxBoundedGeo.h.

---

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

• BoxBoundedGeo.h
• BoxBoundedGeo.cxx