geoalgo::Vector Class Reference

#include "GeoVector.h"

Inheritance diagram for geoalgo::Vector:

Public Member Functions
	Vector ()
	Default ctor. More...
	Vector (size_t n)
	Ctor to instantiate with invalid value. More...
	Vector (const std::vector< double > &obj)
	Default ctor w/ a bare std::vector<double> More...
	Vector (double const x, double const y)
	ctor w/ x & y More...
	Vector (double const x, double const y, double const z)
	ctor w/ x, y & z More...
	Vector (TVector3 const &pt)
	ctor w/ TVector3 More...
	Vector (TLorentzVector const &pt)
	ctor w/ TLorentzVector More...
void	Normalize ()
	Normalize itself. More...
bool	IsValid () const
	Check if point is valid. More...
double	SqLength () const
	Compute the squared length of the vector. More...
double	Length () const
	Compute the length of the vector. More...
Vector	Dir () const
	Return a direction unit vector. More...
double	Phi () const
	Compute the angle Phi. More...
double	Theta () const
	Compute the angle theta. More...
double	SqDist (Vector const &obj) const
	Compute the squared distance to another vector. More...
double	Dist (Vector const &obj) const
	Compute the distance to another vector. More...
double	Dot (Vector const &obj) const
Vector	Cross (Vector const &obj) const
	Compute a dot product of two vectors. More...
double	Angle (Vector const &obj) const
	Compute a cross product of two vectors. More...
TLorentzVector	ToTLorentzVector () const
	Compute an opening angle w.r.t. the given vector. More...
void	compat (Vector const &obj) const
	Dimensional check for a compatibility. More...
void	RotateX (double const &theta)
	rotation operations More...
void	RotateY (double const &theta)
void	RotateZ (double const &theta)
Vector &	operator+= (Vector const &rhs)
Vector &	operator-= (Vector const &rhs)
Vector &	operator*= (double const rhs)
Vector &	operator/= (double const rhs)
Vector &	operator= (Vector const &rhs)
Vector	operator+ (Vector const &rhs) const
Vector	operator- (Vector const &rhs) const
double	operator* (Vector const &rhs) const
Vector	operator* (double const &rhs) const
Vector	operator/ (double const &rhs) const
bool	operator< (Vector const &rhs) const
bool	operator< (double const &rhs) const
bool	operator== (Vector const &rhs) const
bool	operator!= (Vector const &rhs) const

Public Attributes
T	elements
	STL member. More...

Protected Member Functions
double	_SqDist_ (Vector const &obj) const
	Compute the squared-distance to another vector w/o dimension check. More...
double	_Dist_ (Vector const &obj) const
	Compute the distance to another vector w/o dimension check. More...
double	_Dot_ (Vector const &obj) const
	Compute a dot product w/o dimention check. More...
Vector	_Cross_ (Vector const &obj) const
	Compute a cross product w/o dimension check. More...
double	_Angle_ (Vector const &obj) const
	Compute the angle in degrees between 2 vectors w/o dimension check. More...

Friends
class	Trajectory
class	HalfLine
class	LineSegment
class	Sphere
class	GeoAlgo
std::ostream &	operator<< (std::ostream &o,::geoalgo::Vector const &a)
	Streamer. More...

Detailed Description

This class represents an n-dimensional vector

Definition at line 33 of file GeoVector.h.

Constructor & Destructor Documentation

geoalgo::Vector::Vector ( )

inline

Default ctor.

Definition at line 42 of file GeoVector.h.

Referenced by Vector().

: std::vector<double>() {}

geoalgo::Vector::Vector ( size_t  n )

inline

Ctor to instantiate with invalid value.

Definition at line 45 of file GeoVector.h.

: std::vector<double>(n, kINVALID_DOUBLE) {}

geoalgo::Vector::Vector ( const std::vector< double > &  obj )

inline

Default ctor w/ a bare std::vector<double>

Definition at line 48 of file GeoVector.h.

References _Angle_(), _Cross_(), _Dist_(), _Dot_(), _SqDist_(), Angle(), compat(), Cross(), Dir(), Dist(), Dot(), IsValid(), Length(), Normalize(), Phi(), pt, RotateX(), RotateY(), RotateZ(), SqDist(), SqLength(), Theta(), ToTLorentzVector(), Vector(), x, y, and z.

: std::vector<double>(obj) {}

geoalgo::Vector::Vector	(	double const	x,
		double const	y
	)

ctor w/ x & y

Definition at line 10 of file GeoVector.cxx.

References x, and y.

                                               : Vector(2)
  {
    (*this)[0] = x;
    (*this)[1] = y;
  }

geoalgo::Vector::Vector	(	double const	x,
		double const	y,
		double const	z
	)

ctor w/ x, y & z

Definition at line 16 of file GeoVector.cxx.

References x, y, and z.

                                                               : Vector(3)
  {
    (*this)[0] = x;
    (*this)[1] = y;
    (*this)[2] = z;
  }

geoalgo::Vector::Vector

(

TVector3 const &

pt

)

ctor w/ TVector3

Definition at line 23 of file GeoVector.cxx.

                                   : Vector(3)
  {
    (*this)[0] = pt[0];
    (*this)[1] = pt[1];
    (*this)[2] = pt[2];
  }

geoalgo::Vector::Vector

(

TLorentzVector const &

pt

)

ctor w/ TLorentzVector

Definition at line 30 of file GeoVector.cxx.

                                         : Vector(3)
  {
    (*this)[0] = pt[0];
    (*this)[1] = pt[1];
    (*this)[2] = pt[2];
  }

Member Function Documentation

double geoalgo::Vector::_Angle_ ( Vector const &  obj ) const

protected

Compute the angle in degrees between 2 vectors w/o dimension check.

Definition at line 165 of file GeoVector.cxx.

References _Dot_(), and Length().

Referenced by Angle(), and Vector().

  {
    return acos(_Dot_(obj) / Length() / obj.Length());
  }

Vector geoalgo::Vector::_Cross_ ( Vector const &  obj ) const

protected

Compute a cross product w/o dimension check.

Definition at line 156 of file GeoVector.cxx.

Referenced by Cross(), and Vector().

  {
    Vector res(3);
    res[0] = (*this)[1] * obj[2] - obj[1] * (*this)[2];
    res[1] = (*this)[2] * obj[0] - obj[2] * (*this)[0];
    res[2] = (*this)[0] * obj[1] - obj[0] * (*this)[1];
    return res;
  }

double geoalgo::Vector::_Dist_ ( Vector const &  obj ) const

protected

Compute the distance to another vector w/o dimension check.

Definition at line 146 of file GeoVector.cxx.

References _SqDist_().

Referenced by geoalgo::Sphere::Contain(), and Vector().

  {
    return sqrt(_SqDist_(obj));
  }

double geoalgo::Vector::_Dot_ ( Vector const &  obj ) const

protected

Compute a dot product w/o dimention check.

Definition at line 151 of file GeoVector.cxx.

Referenced by _Angle_(), Dot(), and Vector().

  {
    return (*this) * obj;
  }

double geoalgo::Vector::_SqDist_ ( Vector const &  obj ) const

protected

Compute the squared-distance to another vector w/o dimension check.

Definition at line 138 of file GeoVector.cxx.

References larg4::dist(), and util::size().

Referenced by _Dist_(), geoalgo::GeoAlgo::_SqDist_(), geoalgo::GeoAlgo::Intersection(), SqDist(), and Vector().

  {
    double dist = 0;
    for (size_t i = 0; i < size(); ++i)
      dist += ((*this)[i] - obj[i]) * ((*this)[i] - obj[i]);
    return dist;
  }

double geoalgo::Vector::Angle

(

Vector const &

obj

)

const

Compute a cross product of two vectors.

Definition at line 101 of file GeoVector.cxx.

References _Angle_(), compat(), and util::size().

Referenced by geoalgo::Cylinder::Contain(), and Vector().

  {
    compat(obj);
    if (size() != 2 && size() != 3)
      throw GeoAlgoException("<<Angle>> only possible for 2 or 3-dimensional vectors!");
    return _Angle_(obj);
  }

void geoalgo::Vector::compat

(

Vector const &

obj

)

const

Dimensional check for a compatibility.

Definition at line 128 of file GeoVector.cxx.

References util::size().

Referenced by Angle(), geoalgo::GeoAlgo::ClosestPt(), geoalgo::GeoAlgo::commonOrigin(), geoalgo::Sphere::Contain(), Dot(), operator<(), operator==(), SqDist(), geoalgo::GeoAlgo::SqDist(), and Vector().

  {
    if (size() != obj.size()) {
      std::ostringstream msg;
      msg << "<<" << __FUNCTION__ << ">>"
          << " size mismatch: " << size() << " != " << obj.size() << std::endl;
      throw GeoAlgoException(msg.str());
    }
  }

Vector geoalgo::Vector::Cross

(

Vector const &

obj

)

const

Compute a dot product of two vectors.

Definition at line 79 of file GeoVector.cxx.

References _Cross_(), and util::size().

Referenced by Vector().

  {

    if (size() != 3 || obj.size() != 3)

      throw GeoAlgoException("<<Cross>> only possible for 3-dimensional vectors!");

    return _Cross_(obj);
  }

Vector geoalgo::Vector::Dir

(

)

const

Return a direction unit vector.

Definition at line 121 of file GeoVector.cxx.

References Length().

Referenced by Vector().

  {
    Vector res(*this);
    res /= res.Length();
    return res;
  }

double geoalgo::Vector::Dist

(

Vector const &

obj

)

const

Compute the distance to another vector.

Definition at line 68 of file GeoVector.cxx.

References SqDist().

Referenced by geoalgo::Sphere::Sphere(), and Vector().

  {
    return sqrt(SqDist(obj));
  }

double geoalgo::Vector::Dot

(

Vector const &

obj

)

const

Definition at line 73 of file GeoVector.cxx.

References _Dot_(), and compat().

Referenced by geoalgo::GeoAlgo::_commonOrigin_(), geoalgo::Sphere::Sphere(), and Vector().

  {
    compat(obj);
    return _Dot_(obj);
  }

bool geoalgo::Vector::IsValid

(

)

const

Check if point is valid.

Definition at line 37 of file GeoVector.cxx.

References geoalgo::kINVALID_DOUBLE.

Referenced by geoalgo::GeoAlgo::Intersection(), and Vector().

  {

    for (auto const& v : (*this)) {
      // if any point is different from kINVALID_DOUBLE
      // then the point is valid
      if (v != kINVALID_DOUBLE) return true;
    }

    return false;
  }

double geoalgo::Vector::Length

(

)

const

Compute the length of the vector.

Definition at line 57 of file GeoVector.cxx.

References SqLength().

Referenced by _Angle_(), geoalgo::GeoAlgo::_ClosestPt_(), Dir(), Normalize(), geoalgo::HalfLine::Normalize(), operator<(), geoalgo::Sphere::Sphere(), and Vector().

  {
    return sqrt(SqLength());
  }

void geoalgo::Vector::Normalize

(

)

Normalize itself.

Definition at line 116 of file GeoVector.cxx.

References Length().

Referenced by geoalgo::GeoAlgo::_commonOrigin_(), geoalgo::GeoAlgo::_RemainingPoints_(), and Vector().

  {
    (*this) /= this->Length();
  }

bool geoalgo::Vector::operator!= ( Vector const &  rhs ) const

inline

Definition at line 187 of file GeoVector.h.

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

double geoalgo::Vector::operator* ( Vector const &  rhs ) const

inline

Definition at line 147 of file GeoVector.h.

References util::size().

    {
      double res = 0;
      for (size_t i = 0; i < size(); ++i)
        res += (*this)[i] * rhs[i];
      return res;
    }

Vector geoalgo::Vector::operator* ( double const &  rhs ) const

inline

Definition at line 155 of file GeoVector.h.

    {
      Vector res((*this));
      res *= rhs;
      return res;
    }

Vector& geoalgo::Vector::operator*= ( double const  rhs )

inline

Definition at line 111 of file GeoVector.h.

    {
      for (auto& v : *this)
        v *= rhs;
      return *this;
    }

Vector geoalgo::Vector::operator+ ( Vector const &  rhs ) const

inline

Definition at line 133 of file GeoVector.h.

    {
      Vector res((*this));
      res += rhs;
      return res;
    }

Vector& geoalgo::Vector::operator+= ( Vector const &  rhs )

inline

Definition at line 97 of file GeoVector.h.

References util::size().

    {
      for (size_t i = 0; i < size(); ++i)
        (*this)[i] += rhs[i];
      return *this;
    }

Vector geoalgo::Vector::operator- ( Vector const &  rhs ) const

inline

Definition at line 140 of file GeoVector.h.

    {
      Vector res((*this));
      res -= rhs;
      return res;
    }

Vector& geoalgo::Vector::operator-= ( Vector const &  rhs )

inline

Definition at line 104 of file GeoVector.h.

References util::size().

    {
      for (size_t i = 0; i < size(); ++i)
        (*this)[i] -= rhs[i];
      return *this;
    }

Vector geoalgo::Vector::operator/ ( double const &  rhs ) const

inline

Definition at line 162 of file GeoVector.h.

    {
      Vector res((*this));
      res /= rhs;
      return res;
    }

Vector& geoalgo::Vector::operator/= ( double const  rhs )

inline

Definition at line 118 of file GeoVector.h.

    {
      for (auto& v : *this)
        v /= rhs;
      return *this;
    }

bool geoalgo::Vector::operator< ( Vector const &  rhs ) const

inline

Definition at line 169 of file GeoVector.h.

References compat(), and util::size().

    {
      compat(rhs);
      for (size_t i = 0; i < size(); ++i)
        if ((*this)[i] < rhs[i]) return true;
      return false;
    }

bool geoalgo::Vector::operator< ( double const &  rhs ) const

inline

Definition at line 177 of file GeoVector.h.

References Length().

{ return Length() < rhs; }

Vector& geoalgo::Vector::operator= ( Vector const &  rhs )

inline

Definition at line 125 of file GeoVector.h.

    {
      this->resize(rhs.size());
      for (size_t i = 0; i < rhs.size(); ++i)
        (*this)[i] = rhs[i];
      return (*this);
    }

bool geoalgo::Vector::operator== ( Vector const &  rhs ) const

inline

Definition at line 179 of file GeoVector.h.

References compat(), and util::size().

    {
      compat(rhs);
      for (size_t i = 0; i < size(); ++i)
        if ((*this)[i] != rhs[i]) return false;
      return true;
    }

double geoalgo::Vector::Phi

(

)

const

Compute the angle Phi.

Definition at line 89 of file GeoVector.cxx.

Referenced by Vector().

  {
    return (*this)[0] == 0.0 && (*this)[1] == 0.0 ? 0.0 : atan2((*this)[1], (*this)[0]);
  }

void geoalgo::Vector::RotateX

(

double const &

theta

)

rotation operations

Definition at line 170 of file GeoVector.cxx.

Referenced by Vector().

  {

    double c = cos(theta);
    double s = sin(theta);

    double ynew = (*this)[1] * c - (*this)[2] * s;
    double znew = (*this)[1] * s + (*this)[2] * c;

    (*this)[1] = ynew;
    (*this)[2] = znew;

    return;
  }

void geoalgo::Vector::RotateY

(

double const &

theta

)

Definition at line 185 of file GeoVector.cxx.

Referenced by Vector().

  {

    double c = cos(theta);
    double s = sin(theta);

    double xnew = (*this)[0] * c + (*this)[2] * s;
    double znew = -(*this)[0] * s + (*this)[2] * c;

    (*this)[0] = xnew;
    (*this)[2] = znew;

    return;
  }

void geoalgo::Vector::RotateZ

(

double const &

theta

)

Definition at line 200 of file GeoVector.cxx.

Referenced by Vector().

  {

    double c = cos(theta);
    double s = sin(theta);

    double xnew = (*this)[0] * c - (*this)[1] * s;
    double ynew = (*this)[0] * s + (*this)[1] * c;

    (*this)[0] = xnew;
    (*this)[1] = ynew;

    return;
  }

double geoalgo::Vector::SqDist

(

Vector const &

obj

)

const

Compute the squared distance to another vector.

Definition at line 62 of file GeoVector.cxx.

References _SqDist_(), and compat().

Referenced by Dist(), and Vector().

  {
    compat(obj);
    return _SqDist_(obj);
  }

double geoalgo::Vector::SqLength

(

)

const

Compute the squared length of the vector.

Definition at line 49 of file GeoVector.cxx.

Referenced by geoalgo::GeoAlgo::_ClosestPt_(), geoalgo::GeoAlgo::_SqDist_(), Length(), and Vector().

  {
    double res = 0;
    for (auto const& v : (*this))
      res += v * v;
    return res;
  }

double geoalgo::Vector::Theta

(

)

const

Compute the angle theta.

Definition at line 94 of file GeoVector.cxx.

References util::size().

Referenced by Vector().

  {
    if (size() != 3) throw GeoAlgoException("<<Theta>> Only possible for 3-dimensional vectors!");

    return (*this).Length() == 0.0 ? 0.0 : acos((*this)[2] / (*this).Length());
  }

TLorentzVector geoalgo::Vector::ToTLorentzVector

(

)

const

Compute an opening angle w.r.t. the given vector.

Convert geovector to TLorentzVector (with 4th element set equal to 0)

Definition at line 109 of file GeoVector.cxx.

References util::size().

Referenced by Vector().

  {
    if (size() != 3)
      throw GeoAlgoException("<<ToTLorentsVector>> only possible for 3-dimensional vectors!");
    return TLorentzVector((*this)[0], (*this)[1], (*this)[2], 0.);
  }

Friends And Related Function Documentation

friend class GeoAlgo

friend

Definition at line 38 of file GeoVector.h.

friend class HalfLine

friend

Definition at line 35 of file GeoVector.h.

friend class LineSegment

friend

Definition at line 36 of file GeoVector.h.

std::ostream& operator<< ( std::ostream &  o, ::geoalgo::Vector const &  a  )

friend

Streamer.

Definition at line 191 of file GeoVector.h.

    {
      o << "Vector (";
      for (auto const& v : a)
        o << v << " ";
      o << ")";
      return o;
    }

friend class Sphere

friend

Definition at line 37 of file GeoVector.h.

friend class Trajectory

friend

Definition at line 34 of file GeoVector.h.

Member Data Documentation

T std::vector< T >::elements

inherited

STL member.

---

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

• larcorealg/v10_00_02/source/larcorealg/GeoAlgo/GeoVector.h
• larcorealg/v10_00_02/source/larcorealg/GeoAlgo/GeoVector.cxx