G4Ray Class Reference

#include <G4Ray.hh>

Public Member Functions
	G4Ray ()
	G4Ray (const G4Point3D &start0, const G4Vector3D &dir0)
	~G4Ray ()
G4Point3D	GetPoint (G4double i) const
G4double	GetPPoint (const G4Point3D &p) const
const G4Vector3D &	GetDir () const
const G4Point3D &	GetStart () const
void	SetDir (const G4Vector3D &dir0)
void	SetStart (const G4Point3D &start0)
const G4Plane &	GetPlane (G4int number_of_plane) const
void	RayCheck ()
void	CreatePlanes ()
G4double	P2 (G4double x) const
G4int	NearZero (G4double val, G4double epsilon) const
void	MatVecOrtho (register G4Vector3D &out, register const G4Vector3D &in)
void	Vsetall (G4Vector3D &a, G4double s)
void	Init (const G4Point3D &start0, const G4Vector3D &dir0)

Static Public Member Functions
static G4int	CalcPlane3Pts (G4Plane &plane, const G4Point3D &a, const G4Point3D &b, const G4Point3D &c)
static void	Vcross (G4Plane &a, const G4Vector3D &b, const G4Vector3D &c)
static void	Vcross (G4Vector3D &a, const G4Vector3D &b, const G4Vector3D &c)
static void	Vmove (G4Point3D &a, const G4Point3D &b)
static void	Vadd2 (G4Point3D &a, const G4Point3D &b, const G4Vector3D &c)
static void	Vsub2 (G4Vector3D &a, const G4Point3D &b, const G4Point3D &c)
static void	Vscale (G4Plane &a, const G4Plane &b, G4double c)
static G4double	Vdot (const G4Plane &a, const G4Point3D &b)
static G4double	Magsq (const G4Plane &a)
static G4double	Magnitude (const G4Plane &a)

Detailed Description

Definition at line 48 of file G4Ray.hh.

Constructor & Destructor Documentation

G4Ray() [1/2]

G4Ray::G4Ray

(

)

Definition at line 39 of file G4Ray.cc.

  : r_min(0.), r_max(0.)
{
}

G4Ray() [2/2]

G4Ray::G4Ray	(	const G4Point3D &	start0,
		const G4Vector3D &	dir0
	)

Definition at line 44 of file G4Ray.cc.

{
  Init(start0, dir0);
}

~G4Ray()

G4Ray::~G4Ray

(

)

Definition at line 49 of file G4Ray.cc.

{
}

Member Function Documentation

CalcPlane3Pts()

G4int G4Ray::CalcPlane3Pts ( G4Plane &  plane, const G4Point3D &  a, const G4Point3D &  b, const G4Point3D &  c  )

static

Definition at line 212 of file G4Ray.cc.

{
  // Creates the two orthogonal planes which are needed in projecting the
  // surface into 2D.
 
  G4Vector3D    B_A;
  G4Vector3D    C_A;
  G4Vector3D    C_B;
  
  register G4double mag;
 
  Vsub2( B_A, b, a );
  Vsub2( C_A, c, a );
  Vsub2( C_B, c, b );
 
  Vcross( plane1, B_A, C_A );
 
  //    Ensure unit length Normal 
  mag = Magnitude(plane1);
  if( mag  <= SQRT_SMALL_FASTF )
  {
    return(-1);//    FAIL 
  }
  
  mag = 1/mag;
  
  G4Plane pl2(plane1);
  Vscale( plane1, pl2, mag );
 
  //     Find distance from the origin to the G4Plane 
  plane1.d = Vdot( plane1, a );
  
  return(0);    //ok
}

Referenced by CreatePlanes(), and G4FPlane::G4FPlane().

CreatePlanes()

void G4Ray::CreatePlanes

(

)

Definition at line 63 of file G4Ray.cc.

{
  // Creates two orthogonal planes(plane1,plane2) the ray (rray) 
  // situated in the intersection of the planes. The planes are 
  // used to project the surface (nurb) in two dimensions.
  
  G4Vector3D RayDir    = dir;
  G4Point3D  RayOrigin = start;
 
  G4Point3D  p1, p2, p3, p4;
  G4Vector3D dir1, dir2;
  G4Vector3D invdir = G4Vector3D( PINFINITY );
  
  if(!NearZero(RayDir.x(), SQRT_SMALL_FASTF)) 
    { invdir.setX(1.0 / RayDir.x()); }
    
  if(!NearZero(RayDir.y(), SQRT_SMALL_FASTF)) 
    { invdir.setY(1.0 / RayDir.y()); }
    
  if(!NearZero(RayDir.z(), SQRT_SMALL_FASTF)) 
    { invdir.setZ(1.0 / RayDir.z()); }
 
  MatVecOrtho(dir1, RayDir);
  
  Vcross( dir2, RayDir, dir1);
  Vmove(p1, RayOrigin);
  Vadd2(p2, RayOrigin, RayDir);
  Vadd2(p3, RayOrigin, dir1);
  Vadd2(p4, RayOrigin, dir2);
 
  CalcPlane3Pts( plane1, p1, p3, p2);
  CalcPlane3Pts( plane2, p1, p2, p4);
}

Referenced by G4FPlane::CalcNormal().

GetDir()

const G4Vector3D & G4Ray::GetDir ( ) const

inline

Referenced by G4FCylindricalSurface::Intersect(), G4SphericalSurface::Intersect(), G4ToroidalSurface::Intersect(), G4FPlane::Intersect(), G4ConicalSurface::Intersect(), G4CylindricalSurface::Intersect(), G4FConicalSurface::Intersect(), G4BREPSolid::Intersect(), G4BREPSolid::IsConvex(), G4BREPSolid::RemoveHiddenFaces(), and G4BoundingBox3D::Test().

GetPlane()

const G4Plane & G4Ray::GetPlane

(

G4int

number_of_plane

)

const

Definition at line 54 of file G4Ray.cc.

{
  if(number_of_plane==1)
    { return plane2; }
  else
    { return plane1; }
}

GetPoint()

G4Point3D G4Ray::GetPoint ( G4double  i ) const

inline

Referenced by G4SphericalSurface::Intersect(), G4ConicalSurface::Intersect(), and G4CylindricalSurface::Intersect().

GetPPoint()

G4double G4Ray::GetPPoint ( const G4Point3D &  p ) const

inline

GetStart()

const G4Point3D & G4Ray::GetStart ( ) const

inline

Referenced by G4FCylindricalSurface::Intersect(), G4SphericalSurface::Intersect(), G4ToroidalSurface::Intersect(), G4FPlane::Intersect(), G4ConicalSurface::Intersect(), G4CylindricalSurface::Intersect(), G4FConicalSurface::Intersect(), G4BREPSolid::Intersect(), and G4BoundingBox3D::Test().

Init()

void G4Ray::Init	(	const G4Point3D &	start0,
		const G4Vector3D &	dir0
	)

Referenced by G4Ray().

Magnitude()

G4double G4Ray::Magnitude ( const G4Plane &  a )

static

Definition at line 329 of file G4Ray.cc.

{
  return (std::sqrt( Magsq( a )) );
}

Referenced by CalcPlane3Pts().

Magsq()

G4double G4Ray::Magsq ( const G4Plane &  a )

static

Definition at line 323 of file G4Ray.cc.

{
  return ( a.a * a.a + a.b * a.b + a.c *a.c );
}

Referenced by Magnitude().

MatVecOrtho()

void G4Ray::MatVecOrtho	(	register G4Vector3D &	out,
		register const G4Vector3D &	in
	)

Definition at line 98 of file G4Ray.cc.

{
  register G4double f;
  G4int             i_Which;
 
  if( NearZero(in.x(), 0.0001)
   && NearZero(in.y(), 0.0001)
   && NearZero(in.z(), 0.0001) )  
  {
    Vsetall( out, 0 );
    return;
  }
 
  //      Find component closest to zero 
  f = std::fabs(in.x());
  i_Which=0;
  
  if( std::fabs(in.y()) < f )
  {
    f = std::fabs(in.y());
    i_Which=1;
  }
  
  if( std::fabs(in.z()) < f )
  {
    i_Which=2;
  }
  
  if(!i_Which)
  {
    f = std::sqrt((in.y())*(in.y())+(in.z())*(in.z()));    // hypot(in.y(),in.z())
  }
  else
  {
    if(i_Which==1)
    {
      f = std::sqrt((in.z())*(in.z())+(in.x())*(in.x()));  // hypot(in.z(),in.x())
    }
    else
    {
      f = std::sqrt((in.x())*(in.x())+(in.y())*(in.y()));  // hypot(in.x(),in.y())
    }
  }
  if( NearZero( f, SMALL ) )
  {
    Vsetall( out, 0 );
    return;
  }
    
  f = 1.0/f;
    
  if(!i_Which)
  {
    out.setX(0.0);
    out.setY(-in.z()*f);
    out.setZ( in.y()*f);
  }
  else
  {
    if(i_Which==1)
    {
      out.setY(0.0);
      out.setZ(-in.x()*f);
      out.setX( in.y()*f);
    }
    else
    {
      out.setZ(0.0);
      out.setX(-in.z()*f);
      out.setY( in.y()*f);
    }
  }
} 

Referenced by CreatePlanes().

NearZero()

G4int G4Ray::NearZero ( G4double  val, G4double  epsilon  ) const

inline

Referenced by CreatePlanes(), and MatVecOrtho().

P2()

G4double G4Ray::P2 ( G4double  x ) const

inline

RayCheck()

void G4Ray::RayCheck

(

)

Definition at line 251 of file G4Ray.cc.

{
  // Check that the ray has a G4Vector3D...
  if (dir==G4Vector3D(0, 0, 0)) 
  {
    G4Exception("G4Ray::RayCheck()", "GeomSolids0002", FatalException,
                "Invalid zero direction given !");
  }
 
  // Make sure that the vector is unit length
  dir= dir.unit();
  r_min = 0;
  r_max = 0;
}

Referenced by G4FPlane::CalcNormal().

SetDir()

void G4Ray::SetDir ( const G4Vector3D &  dir0 )

inline

SetStart()

void G4Ray::SetStart ( const G4Point3D &  start0 )

inline

Vadd2()

void G4Ray::Vadd2 ( G4Point3D &  a, const G4Point3D &  b, const G4Vector3D &  c  )

static

Definition at line 291 of file G4Ray.cc.

{
  a.setX(b.x() + c.x()) ;
  a.setY(b.y() + c.y()) ;
  a.setZ(b.z() + c.z()) ;
}   

Referenced by CreatePlanes().

Vcross() [1/2]

void G4Ray::Vcross ( G4Plane &  a, const G4Vector3D &  b, const G4Vector3D &  c  )

static

Definition at line 267 of file G4Ray.cc.

{ 
  a.a = b.y()  * c.z()  - b.z()  * c.y() ;
  a.b = b.z()  * c.x()  - b.x()  * c.z() ;
  a.c = b.x()  * c.y()  - b.y()  * c.x() ;
}

Referenced by CalcPlane3Pts(), and CreatePlanes().

Vcross() [2/2]

void G4Ray::Vcross ( G4Vector3D &  a, const G4Vector3D &  b, const G4Vector3D &  c  )

static

Definition at line 275 of file G4Ray.cc.

{ 
  a.setX(b.y()  * c.z()  - b.z()  * c.y()) ;
  a.setY(b.z()  * c.x()  - b.x()  * c.z()) ;
  a.setZ(b.x()  * c.y()  - b.y()  * c.x()) ;
}

Vdot()

G4double G4Ray::Vdot ( const G4Plane &  a, const G4Point3D &  b  )

static

Definition at line 315 of file G4Ray.cc.

{
  return (a.a * b.x() + 
      a.b * b.y() + 
      a.c * b.z());
}

Referenced by CalcPlane3Pts().

Vmove()

void G4Ray::Vmove ( G4Point3D &  a, const G4Point3D &  b  )

static

Definition at line 283 of file G4Ray.cc.

{ 
  a.setX(b.x());
  a.setY(b.y());
  a.setZ(b.z());
}

Referenced by CreatePlanes().

Vscale()

void G4Ray::Vscale ( G4Plane &  a, const G4Plane &  b, G4double  c  )

static

Definition at line 307 of file G4Ray.cc.

{ 
  a.a = b.a * c;
  a.b = b.b * c;
  a.c = b.c * c;
}

Referenced by CalcPlane3Pts().

Vsetall()

void G4Ray::Vsetall ( G4Vector3D &  a, G4double  s  )

inline

Referenced by MatVecOrtho().

Vsub2()

void G4Ray::Vsub2 ( G4Vector3D &  a, const G4Point3D &  b, const G4Point3D &  c  )

static

Definition at line 299 of file G4Ray.cc.

{
  a.setX(b.x() - c.x());
  a.setY(b.y() - c.y());
  a.setZ(b.z() - c.z());
}

Referenced by CalcPlane3Pts().

---

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

• G4Ray.hh
• G4Ray.cc