Heed::vec Class Reference

#include <vec.h>

Inheritance diagram for Heed::vec:

Public Member Functions
	vec (vfloat xx, vfloat yy, vfloat zz)
	Constructor.
	vec ()
	Default constructor.
virtual	~vec ()
	Destructor.
vfloat	length () const
vfloat	length2 () const
vec	down_new (const basis *fabas)
void	down (const basis *fabas)
vec	up_new (const basis *fabas_new)
void	up (const basis *fabas_new)
vec	down_new (const abssyscoor *fasc)
void	down (const abssyscoor *fasc)
	Convert numbering representation of object to basical system of fasc.
vec	up_new (const abssyscoor *fasc)
void	up (const abssyscoor *fasc)
	Convert numbering representation of objects to new system.
vec	turn_new (const vec &dir, vfloat angle)
void	turn (const vec &dir, vfloat angle)
	Turn this vector.
void	shift (const vec &dir)
void	random_round_vec ()
	Generate random unit vector in plane perpendicular to z-axis.
void	random_conic_vec (double theta)
void	random_sfer_vec ()
vec	operator/ (vfloat p) const
Public Member Functions inherited from Heed::absref
virtual	~absref ()
	Destructor.
virtual void	down (const abssyscoor *fasc)
	Convert numbering representation of object to basical system of fasc.
virtual void	up (const abssyscoor *fasc)
	Convert numbering representation of objects to new system.
virtual void	turn (const vec &dir, vfloat angle)
	Turn around axis doing via center of coordinate system along dir.
virtual void	shift (const vec &dir)

Public Attributes
vfloat	x
vfloat	y
vfloat	z

Friends
vec	operator* (const vec &v, vfloat p)
vec	operator*= (vec &v, vfloat p)
vec	operator* (vfloat p, const vec &v)
vec	operator/= (vec &v, vfloat p)
vec	operator+ (const vec &r1, const vec &r2)
vec &	operator+= (vec &r1, const vec &r2)
vec	operator- (const vec &r1, const vec &r2)
vec	operator-= (vec &r1, const vec &r2)
vec	operator- (const vec &r)
vfloat	operator* (const vec &r1, const vec &r2)
vec	operator|| (const vec &r1, const vec &r2)
	Vector product.
int	operator== (const vec &r1, const vec &r2)
	Return 1 if precisely the same vectors and 0 otherwise.
int	operator!= (const vec &r1, const vec &r2)
	Return 0 if precisely the same vectors and 1 otherwise.
bool	apeq (const vec &r1, const vec &r2, vfloat prec)
	Return true if two vectors are approximately the same.
vec	unit_vec (const vec &v)
vfloat	cos2vec (const vec &r1, const vec &r2)
vfloat	ang2vec (const vec &r1, const vec &r2)
vec	project_to_plane (const vec &r, const vec &normal)
vfloat	ang2projvec (const vec &r1, const vec &r2, const vec &normal)
vfloat	sin2vec (const vec &r1, const vec &r2)
int	check_par (const vec &r1, const vec &r2, vfloat prec)
int	check_perp (const vec &r1, const vec &r2, vfloat prec)
vec	switch_xyz (const vec &)

Detailed Description

Vector. Each vector is presented by three components and corresponds to some basis. // The components are a projection of the vector to unit basis vectors. If bas==NULL then it is the primary basis. So the concept of vector is more primary concept with comparison of basis, since one can not postulate basis while vectors do not exist.

Definition at line 186 of file vec.h.

Constructor & Destructor Documentation

vec() [1/2]

Heed::vec::vec ( vfloat  xx, vfloat  yy, vfloat  zz  )

inline

Constructor.

Definition at line 189 of file vec.h.

                                       {
    x = xx;
    y = yy;
    z = zz;
  }

vec() [2/2]

Heed::vec::vec ( )

inline

Default constructor.

Definition at line 195 of file vec.h.

        {
    x = 0;
    y = 0;
    z = 0;
  }

Referenced by operator/(), and turn_new().

~vec()

virtual Heed::vec::~vec ( )

inlinevirtual

Destructor.

Definition at line 201 of file vec.h.

{}

Member Function Documentation

down() [1/2]

void Heed::vec::down ( const abssyscoor *  fasc )

virtual

Convert numbering representation of object to basical system of fasc.

Reimplemented from Heed::absref.

Definition at line 225 of file vec.cpp.

{ down(fasc->Gabas()); }

down() [2/2]

void Heed::vec::down

(

const basis *

fabas

)

Definition at line 168 of file vec.cpp.

{ *this = this->down_new(fabas); }

Referenced by down(), Heed::point::down(), Heed::HeedParticle::physics(), Heed::HeedParticle_BGM::physics(), and Heed::HeedDeltaElectron::physics_after_new_speed().

down_new() [1/2]

vec Heed::vec::down_new

(

const abssyscoor *

fasc

)

Definition at line 224 of file vec.cpp.

{ return down_new(fasc->Gabas()); }

down_new() [2/2]

vec Heed::vec::down_new

(

const basis *

fabas

)

Definition at line 156 of file vec.cpp.

                                    {
  // pvecerror("vec vec::down_new(void)");
  vec r;
  vec ex = fabas->Gex();
  vec ey = fabas->Gey();
  vec ez = fabas->Gez();
  r.x = x * ex.x + y * ey.x + z * ez.x;
  r.y = x * ex.y + y * ey.y + z * ez.y;
  r.z = x * ex.z + y * ey.z + z * ez.z;
  return r;
}

Referenced by down(), and down_new().

length()

vfloat Heed::vec::length ( ) const

inline

Definition at line 205 of file vec.h.

{ return sqrt(x * x + y * y + z * z); }

Referenced by Heed::basis::basis(), Heed::circumf::check_point_in(), Heed::circumf::circumf(), Heed::plane::cross(), Heed::polyline::cross(), Heed::mparticle::curvature(), Heed::straight::distance(), Heed::trajestep::Gnextpoint(), Heed::trajestep::Gnextpoint1(), Heed::gparticle::gparticle(), Heed::HeedPhoton::HeedPhoton(), Heed::mparticle::new_speed(), Heed::HeedParticle::physics(), Heed::gparticle::print(), Heed::polygon::range(), Heed::splane::range(), Heed::trajestep::trajestep(), and turn_new().

length2()

vfloat Heed::vec::length2 ( ) const

inline

Definition at line 206 of file vec.h.

{ return x * x + y * y + z * z; }

operator/()

vec Heed::vec::operator/ ( vfloat  p ) const

inline

Definition at line 241 of file vec.h.

{ return vec(x / p, y / p, z / p); }

random_conic_vec()

void Heed::vec::random_conic_vec

(

double

theta

)

Generate random unit vector in any direction in conus with symmetry axis along z-axis and with angle theta (radian).

Definition at line 236 of file vec.cpp.

                                       {
  vfloat phi = M_PI * 2.0 * SRANLUX();
  double stheta = sin(theta);
  x = sin(phi) * stheta;
  y = cos(phi) * stheta;
  z = cos(theta);
}

Referenced by Heed::HeedParticle::physics(), and Heed::HeedParticle_BGM::physics().

random_round_vec()

void Heed::vec::random_round_vec

(

)

Generate random unit vector in plane perpendicular to z-axis.

Definition at line 229 of file vec.cpp.

                           {
  const vfloat phi = M_PI * 2.0 * SRANLUX();
  x = sin(phi);
  y = cos(phi);
  z = 0;
}

Referenced by Heed::HeedDeltaElectron::physics_after_new_speed(), and random_sfer_vec().

random_sfer_vec()

void Heed::vec::random_sfer_vec

(

)

Definition at line 244 of file vec.cpp.

                          {
  vfloat cteta = 2.0 * SRANLUX() - 1.0;
  random_round_vec();
  vfloat steta = sqrt(1.0 - cteta * cteta);
  *this = (*this) * steta;
  z = cteta;
}

Referenced by Heed::HeedPhoton::physics_after_new_speed().

shift()

void Heed::vec::shift ( const vec &  dir )

virtual

Reimplemented from Heed::absref.

Definition at line 220 of file vec.cpp.

                            {
  // Not defined for vectors
}

turn()

void Heed::vec::turn ( const vec &  dir, vfloat  angle  )

virtual

Turn this vector.

Reimplemented from Heed::absref.

Definition at line 216 of file vec.cpp.

                                           {
  *this = this->turn_new(dir, angle);
}

Referenced by Heed::trajestep::Gnextpoint(), Heed::trajestep::Gnextpoint1(), Heed::mparticle::new_speed(), Heed::spquadr::pt_angle_rad(), and Heed::splane::range().

turn_new()

vec Heed::vec::turn_new	(	const vec &	dir,
		vfloat	angle
	)

Definition at line 192 of file vec.cpp.

                                              {
  pvecerror("vec turn(vec& dir, vfloat& angle)");
  if ((*this).length() == 0) return vec(0, 0, 0);
  if (check_par(*this, dir, 0.0) != 0) {
    // parallel vectors are not changed
    return *this;
  }
  vfloat dirlen = dir.length();
  check_econd11a(dirlen, == 0, "cannot turn around zero vector", mcerr);
  vec u = dir / dirlen;  // unit vector
  vec constcomp = u * (*this) * u;
  vec ort1 = unit_vec(u || (*this));
  vec ort2 = ort1 || u;
  vec perpcomp = ort2 * (*this) * ort2;
  vfloat len = perpcomp.length();
  // mcout<<" constcomp="<<constcomp<<" ort1="<<ort1<<" ort2="<<ort2;
  ort1 = sin(angle) * len * ort1;
  ort2 = cos(angle) * len * ort2;
  // mcout<<" constcomp="<<constcomp<<" ort1="<<ort1<<" ort2="<<ort2
  //    <<" len="<<len<<" sin(angle)="<<sin(angle)<<" cos(angle)="<<cos(angle)
  //    <<" angle="<<angle<<'\n';
  return constcomp + ort1 + ort2;
}

Referenced by turn().

up() [1/2]

void Heed::vec::up ( const abssyscoor *  fasc )

virtual

Convert numbering representation of objects to new system.

Reimplemented from Heed::absref.

Definition at line 227 of file vec.cpp.

{ up(fasc->Gabas()); }

up() [2/2]

void Heed::vec::up

(

const basis *

fabas_new

)

Definition at line 190 of file vec.cpp.

{ *this = this->up_new(fabas_new); }

Referenced by Heed::manip_absvol::m_check_point_inside(), up(), and Heed::point::up().

up_new() [1/2]

vec Heed::vec::up_new

(

const abssyscoor *

fasc

)

Definition at line 226 of file vec.cpp.

{ return up_new(fasc->Gabas()); }

up_new() [2/2]

vec Heed::vec::up_new

(

const basis *

fabas_new

)

Definition at line 170 of file vec.cpp.

                                      {
  // it is assumed that fabas_new is derivative from old
  pvecerrorp("vec vec::up_new((const basis *pbas)");
  vec r;
  // check_econd11(fabas_new , ==NULL, mcerr);
  // not compiled in IRIX, reason is unkown
  if (fabas_new == NULL) {
    funnw.ehdr(mcerr);
    mcerr << "fabas_new==NULL\n";
    spexit(mcerr);
  }
  vec ex = fabas_new->Gex();
  vec ey = fabas_new->Gey();
  vec ez = fabas_new->Gez();
  r.x = x * ex.x + y * ex.y + z * ex.z;
  r.y = x * ey.x + y * ey.y + z * ey.z;
  r.z = x * ez.x + y * ez.y + z * ez.z;
  return r;
}

Referenced by up(), and up_new().

Friends And Related Function Documentation

ang2projvec

vfloat ang2projvec ( const vec &  r1, const vec &  r2, const vec &  normal  )

friend

Definition at line 136 of file vec.cpp.

                                                                    {
  pvecerror(
      "vfloat ang2projvec(const vec& r1, const vec& r2, const vec& normal)");
  vec rt1 = project_to_plane(r1, normal);
  vec rt2 = project_to_plane(r2, normal);
  if (rt1 == dv0 || rt2 == dv0) {
    vecerror = 1;
    return 0;
  }
  vfloat tang = ang2vec(rt1, rt2);
  if (tang == 0) return tang;  // projections are parallel
  vec at = rt1 || rt2;
  int i = check_par(at, normal, 0.0001);
  // mcout<<"r1="<<r1<<"r2="<<r2<<"normal="<<normal
  //     <<"rt1="<<rt1<<"rt2="<<rt2<<"\ntang="<<tang
  //     <<"\nat="<<at<<" i="<<i<<'\n';
  if (i == -1) return 2.0 * M_PI - tang;
  return tang;  // it works if angle <= PI
}

ang2vec

vfloat ang2vec ( const vec &  r1, const vec &  r2  )

friend

Definition at line 89 of file vec.cpp.

                                             {
  // angle between vectors
  // instead of return acos(cos2vec(r1,r2)); which produces NaN on linux at
  // parallel vectors
  vfloat cs = cos2vec(r1, r2);
  if (vecerror != 0) return 0;
  if (cs > 0.707106781187 || cs < -0.707106781187) {  // 1.0/sqrt(2)
    // pass to sin, it will be more exactly
    vfloat sn = sin2vec(r1, r2);
    if (vecerror != 0) return 0;
    if (cs > 0.0)
      return asin(sn);
    else
      return M_PI - asin(sn);
  }
  return acos(cs);
}

apeq

bool apeq ( const vec &  r1, const vec &  r2, vfloat  prec  )

friend

Return true if two vectors are approximately the same.

check_par

int check_par ( const vec &  r1, const vec &  r2, vfloat  prec  )

friend

Check whether two vectors are parallel, or anti-parallel. Returns: 1 - parallel, -1 - antiparallel, 0 not parallel. Also returns 0 if one or both vectors have zero length. Thus, if angle between vectors < prec, they are parallel.

Referenced by turn_new().

check_perp

int check_perp ( const vec &  r1, const vec &  r2, vfloat  prec  )

friend

Check whether two vectors are perpendicular. Returns: 1 perpendicular, 0 not perpendicular. Also returns 0 if one or both vectors have zero length. Thus, if angle between vectors a > 0.5 * M_PI - max(prec, vprecision) and a < 0.5 * M_PI + max(prec, vprecision), they are perpendicular.

cos2vec

vfloat cos2vec ( const vec &  r1, const vec &  r2  )

friend

Definition at line 66 of file vec.cpp.

                                             {
  // cosinus of angle between vectors
  // If one of vectors has zero length, it returns 2.
  pvecerror("vfloat cos2vec(const vec& r1, const vec& r2)");
  vfloat lr1 = r1.length2();
  vfloat lr2 = r2.length2();
  // mcout<<"cos2vec:\n";
  // Iprintn(mcout, lr1);
  // Iprintn(mcout, lr2);
  if (lr1 == 0 || lr2 == 0) {
    vecerror = 1;
    return 0;
  }
  vfloat cs = r1 * r2;
  int sign = 1;
  if (cs < 0) sign = -1;
  cs = cs * cs;
  cs = sign * sqrt(cs / (lr1 * lr2));
  // mcout<<"r1="<<r1<<"r2="<<r2<<"cos="<<cs<<'\n';
  return cs;
  // return r1*r2/(lr1*lr2);
}

operator!=

int operator!= ( const vec &  r1, const vec &  r2  )

friend

Return 0 if precisely the same vectors and 1 otherwise.

operator* [1/3]

vfloat operator* ( const vec &  r1, const vec &  r2  )

friend

Definition at line 261 of file vec.h.

                                                        {
    return r1.x * r2.x + r1.y * r2.y + r1.z * r2.z;
  }

operator* [2/3]

vec operator* ( const vec &  v, vfloat  p  )

friend

Definition at line 231 of file vec.h.

                                               {
    return vec(v.x * p, v.y * p, v.z * p);
  }

operator* [3/3]

vec operator* ( vfloat  p, const vec &  v  )

friend

Definition at line 238 of file vec.h.

                                               {
    return vec(v.x * p, v.y * p, v.z * p);
  }

operator*=

vec operator*= ( vec &  v, vfloat  p  )

friend

Definition at line 234 of file vec.h.

                                          {
    v = v * p;
    return v;
  }

operator+

vec operator+ ( const vec &  r1, const vec &  r2  )

friend

Definition at line 246 of file vec.h.

                                                     {
    return vec(r1.x + r2.x, r1.y + r2.y, r1.z + r2.z);
  }

operator+=

vec & operator+= ( vec &  r1, const vec &  r2  )

friend

Definition at line 249 of file vec.h.

                                                 {
    r1 = r1 + r2;
    return r1;
  }

operator- [1/2]

vec operator- ( const vec &  r )

friend

Definition at line 260 of file vec.h.

{ return vec(-r.x, -r.y, -r.z); }

operator- [2/2]

vec operator- ( const vec &  r1, const vec &  r2  )

friend

Definition at line 253 of file vec.h.

                                                     {
    return vec(r1.x - r2.x, r1.y - r2.y, r1.z - r2.z);
  }

operator-=

vec operator-= ( vec &  r1, const vec &  r2  )

friend

Definition at line 256 of file vec.h.

                                                {
    r1 = r1 - r2;
    return r1;
  }

operator/=

vec operator/= ( vec &  v, vfloat  p  )

friend

Definition at line 242 of file vec.h.

                                          {
    v = v / p;
    return v;
  }

operator==

int operator== ( const vec &  r1, const vec &  r2  )

friend

Return 1 if precisely the same vectors and 0 otherwise.

operator||

vec operator|| ( const vec &  r1, const vec &  r2  )

friend

Vector product.

Definition at line 265 of file vec.h.

                                                      {
    return vec(r1.y * r2.z - r1.z * r2.y, r1.z * r2.x - r1.x * r2.z,
               r1.x * r2.y - r1.y * r2.x);
  }

project_to_plane

vec project_to_plane ( const vec &  r, const vec &  normal  )

friend

Definition at line 124 of file vec.cpp.

                                                      {
  pvecerror("vec project_to_plane(const vec& r, const vec& normal)");
  vec per(normal || r);
  if (per == dv0) {
    // either one of vectors is 0 or they are parallel
    return dv0;
  }
  vec ax = unit_vec(per || normal);
  vfloat v = ax * r;
  return v * ax;
}

sin2vec

vfloat sin2vec ( const vec &  r1, const vec &  r2  )

friend

Definition at line 107 of file vec.cpp.

                                             {
  // sinus of angle between vectors
  pvecerror("vfloat sin2vec(const vec& r1, const vec& r2)");
  vfloat lr1 = r1.length2();
  vfloat lr2 = r2.length2();
  if (lr1 == 0 || lr2 == 0) {
    vecerror = 1;
    return 0;
  }
  vfloat sn = (r1 || r2).length();
  sn = sn * sn;
  sn = sqrt(sn / (lr1 * lr2));
  // mcout<<"r1="<<r1<<"r2="<<r2<<"sin="<<sn<<'\n';
  return sn;
  // return sin(ang2vec(r1,r2));
}

switch_xyz

vec switch_xyz ( const vec &  )

friend

unit_vec

vec unit_vec ( const vec &  v )

friend

Referenced by turn_new().

Member Data Documentation

x

vfloat Heed::vec::x

Definition at line 203 of file vec.h.

Referenced by Heed::box::check_point_inside(), down_new(), Heed::HeedFieldMap::field_map(), Garfield::TrackHeed::GetCluster(), Heed::HeedFieldMap::inside(), length(), length2(), operator/(), Heed::operator<<(), Heed::HeedDeltaElectron::physics_after_new_speed(), random_conic_vec(), random_round_vec(), Heed::box::range_ext(), Heed::straight::straight(), Garfield::TrackHeed::TransportPhoton(), up_new(), and vec().

y

vfloat Heed::vec::y

Definition at line 203 of file vec.h.

Referenced by Heed::box::check_point_inside(), down_new(), Heed::HeedFieldMap::field_map(), Garfield::TrackHeed::GetCluster(), Heed::HeedFieldMap::inside(), length(), length2(), operator/(), Heed::operator<<(), Heed::HeedDeltaElectron::physics_after_new_speed(), random_conic_vec(), random_round_vec(), Heed::box::range_ext(), Heed::straight::straight(), Garfield::TrackHeed::TransportPhoton(), up_new(), vec(), and Heed::straight::vecdistance().

z

vfloat Heed::vec::z

Definition at line 203 of file vec.h.

Referenced by Heed::box::check_point_inside(), down_new(), Heed::HeedFieldMap::field_map(), Garfield::TrackHeed::GetCluster(), Heed::HeedFieldMap::inside(), length(), length2(), operator/(), Heed::operator<<(), random_conic_vec(), random_round_vec(), random_sfer_vec(), Heed::box::range_ext(), Heed::straight::straight(), Garfield::TrackHeed::TransportPhoton(), up_new(), and vec().

---

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

• vec.h
• vec.cpp