G4Hyperbola Class Reference

#include <G4Hyperbola.hh>

Inheritance diagram for G4Hyperbola:

Public Member Functions
	G4Hyperbola ()
virtual	~G4Hyperbola ()
	G4Hyperbola (const G4Hyperbola &right)
G4Hyperbola &	operator= (const G4Hyperbola &right)
G4Curve *	Project (const G4Transform3D &tr=G4Transform3D::Identity)
G4bool	Tangent (G4CurvePoint &cp, G4Vector3D &v)
G4double	GetPMax () const
G4Point3D	GetPoint (G4double param) const
G4double	GetPPoint (const G4Point3D &p) const
G4double	GetSemiAxis () const
G4double	GetSemiImagAxis () const
void	Init (const G4Axis2Placement3D &position0, G4double semiAxis0, G4double semiImagAxis0)
G4int	IntersectRay2D (const G4Ray &ray)
Public Member Functions inherited from G4Conic
	G4Conic ()
virtual	~G4Conic ()
	G4Conic (const G4Conic &right)
G4Conic &	operator= (const G4Conic &right)
const G4Axis2Placement3D *	GetPosition () const
G4double	GetPShift () const
void	SetPShift (G4double pShift0)
Public Member Functions inherited from G4Curve
	G4Curve ()
virtual	~G4Curve ()
	G4Curve (const G4Curve &c)
G4Curve &	operator= (const G4Curve &c)
G4bool	operator== (const G4Curve &right) const
virtual G4String	GetEntityType () const
virtual G4Curve *	Project (const G4Transform3D &tr=G4Transform3D::Identity)=0
virtual G4bool	Tangent (G4CurvePoint &cp, G4Vector3D &v)=0
virtual G4int	IntersectRay2D (const G4Ray &ray)=0
const G4Point3D &	GetStart () const
const G4Point3D &	GetEnd () const
G4double	GetPStart () const
G4double	GetPEnd () const
void	SetBounds (G4double p1, G4double p2)
void	SetBounds (G4double p1, const G4Point3D &p2)
void	SetBounds (const G4Point3D &p1, G4double p2)
void	SetBounds (const G4Point3D &p1, const G4Point3D &p2)
G4bool	IsBounded () const
G4bool	IsPOn (G4double param) const
void	SetSameSense (G4int sameSense0)
G4int	GetSameSense () const
virtual G4double	GetPMax () const =0
virtual G4Point3D	GetPoint (G4double param) const =0
virtual G4double	GetPPoint (const G4Point3D &p) const =0
const G4BoundingBox3D *	BBox () const
virtual const char *	Name () const
virtual void	SetParentSrfPtr (const G4Surface *)

Protected Member Functions
void	InitBounded ()
virtual void	InitBounded ()=0

Additional Inherited Members
Protected Attributes inherited from G4Conic
G4Axis2Placement3D	position
G4double	pShift
Protected Attributes inherited from G4Curve
G4BoundingBox3D	bBox
G4Point3D	start
G4Point3D	end
G4double	pStart
G4double	pEnd
G4double	pRange
G4bool	bounded
G4int	sameSense
G4double	kCarTolerance

Detailed Description

Definition at line 44 of file G4Hyperbola.hh.

Constructor & Destructor Documentation

G4Hyperbola() [1/2]

G4Hyperbola::G4Hyperbola

(

)

Definition at line 41 of file G4Hyperbola.cc.

  : semiAxis(0.), semiImagAxis(0.), ratioAxisImagAxis(0.), forTangent(0.)
{
}

Referenced by Project().

~G4Hyperbola()

G4Hyperbola::~G4Hyperbola ( )

virtual

Definition at line 46 of file G4Hyperbola.cc.

{
}

G4Hyperbola() [2/2]

G4Hyperbola::G4Hyperbola

(

const G4Hyperbola &

right

)

Definition at line 50 of file G4Hyperbola.cc.

  : G4Conic(), Focus1(right.Focus1), Focus2(right.Focus2),
    ProjFocus1(right.ProjFocus1), ProjFocus2(right.ProjFocus2),
    semiAxis(right.semiAxis), semiImagAxis(right.semiImagAxis),
    ratioAxisImagAxis(right.ratioAxisImagAxis),
    toUnitHyperbola(right.toUnitHyperbola), forTangent(right.forTangent)
{
  pShift    = right.pShift;
  position  = right.position;
  bBox      = right.bBox;
  start     = right.start;
  end       = right.end;
  pStart    = right.pStart;
  pEnd      = right.pEnd;
  pRange    = right.pRange;
  bounded   = right.bounded;
  sameSense = right.sameSense;
}

Member Function Documentation

GetPMax()

G4double G4Hyperbola::GetPMax ( ) const

inlinevirtual

Implements G4Curve.

GetPoint()

G4Point3D G4Hyperbola::GetPoint ( G4double  param ) const

inlinevirtual

Implements G4Curve.

Referenced by InitBounded().

GetPPoint()

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

inlinevirtual

Implements G4Curve.

GetSemiAxis()

G4double G4Hyperbola::GetSemiAxis ( ) const

inline

GetSemiImagAxis()

G4double G4Hyperbola::GetSemiImagAxis ( ) const

inline

Init()

void G4Hyperbola::Init ( const G4Axis2Placement3D &  position0, G4double  semiAxis0, G4double  semiImagAxis0  )

inline

Referenced by Project().

InitBounded()

void G4Hyperbola::InitBounded ( )

protectedvirtual

Implements G4Curve.

Definition at line 157 of file G4Hyperbola.cc.

{
  // the bbox must include the start and endpoints as well as the
  // extreme points if they lie on the curve
  bBox.Init(GetStart(), GetEnd());
 
  // the parameter values
  // belonging to the points with an extreme x, y and z coordinate
  for (G4int i=0; i<3; i++) 
  {
    G4double x_i= position.GetPX()(i);
    
    if (std::abs(x_i) <= 
        G4GeometryTolerance::GetInstance()->GetAngularTolerance()) 
    {
      G4double tanhu= - (semiImagAxis*position.GetPY()(i)) / (semiAxis*x_i);
      
      if (std::abs(tanhu)<=1) 
      {
        G4double u= 0.5*std::log((1+tanhu)/(1-tanhu));  // atanh(tanhu)
    if (IsPOn(u))
      bBox.Extend(GetPoint(u));
      }
    }
  }
}

IntersectRay2D()

G4int G4Hyperbola::IntersectRay2D ( const G4Ray &  ray )

inlinevirtual

Implements G4Curve.

operator=()

G4Hyperbola & G4Hyperbola::operator=

(

const G4Hyperbola &

right

)

Definition at line 69 of file G4Hyperbola.cc.

{
  if (&right == this) return *this;
 
  Focus1 = right.Focus1;
  Focus2 = right.Focus2;
  ProjFocus1   = right.ProjFocus1;
  ProjFocus2   = right.ProjFocus2;
  semiAxis     = right.semiAxis;
  semiImagAxis = right.semiImagAxis;
  ratioAxisImagAxis = right.ratioAxisImagAxis;
  toUnitHyperbola   = right.toUnitHyperbola;
  forTangent        = right.forTangent;
  pShift    = right.pShift;
  position  = right.position;
  bBox      = right.bBox;
  start     = right.start;
  end       = right.end;
  pStart    = right.pStart;
  pEnd      = right.pEnd;
  pRange    = right.pRange;
  bounded   = right.bounded;
  sameSense = right.sameSense;
 
  return *this;
}

Project()

G4Curve * G4Hyperbola::Project ( const G4Transform3D &  tr = G4Transform3D::Identity )

virtual

!!!!!!!!!!!

!!!!!!!!!!!

Implements G4Curve.

Definition at line 96 of file G4Hyperbola.cc.

{
  G4Exception("G4Hyperbola::Project()", "GeomSolids0001",
              FatalException, "Sorry, not yet implemented.");
 
  G4Point3D newLocation= tr*position.GetLocation();
  newLocation.setZ(0);
  G4double axisZ= (tr*position.GetPZ()).unit().z();
 
  if (std::abs(axisZ)<G4GeometryTolerance::GetInstance()->GetAngularTolerance()) 
  {
    return 0;
  }
 
  G4Vector3D newAxis(0, 0, axisZ>0? +1: -1);
 
  // get the parameter of an endpoint of an axis
  // (this is a point the distance of which from the center is extreme)
  G4Vector3D xPrime = tr*position.GetPX();
  xPrime.setZ(0);
  
  G4Vector3D yPrime = tr*position.GetPY();
  yPrime.setZ(0);
  
  G4Vector3D a = G4Vector3D( semiAxis*xPrime );
  G4Vector3D b = G4Vector3D( semiImagAxis*yPrime );
 
  G4double xval = -2*a*b/(a.mag2()+b.mag2());
  
  G4double u= (0.5*std::log((1+xval)/(1-xval)))/2;  // atanh(xval)/2
 
  // get the coordinate axis directions and the semiaxis lengths
  G4Vector3D sAxis= G4Vector3D( a*std::cosh(u)+b*std::sinh(u) );
  
  //!!!!!!!!!!!!
  G4Vector3D sImagAxis= G4Vector3D( a*std::cosh(u+pi/2)+b*std::sinh(u+pi/2) );
  
  //!!!!!!!!!!!!
  G4double   newSemiAxis     = sAxis.mag();
  G4double   newSemiImagAxis = sImagAxis.mag();
  G4Vector3D newRefDirection = sAxis;
 
  // create the new hyperbola
  G4Axis2Placement3D newPosition;
  newPosition.Init(newRefDirection, newAxis, newLocation);
 
  G4Hyperbola* r= new G4Hyperbola;
  r->Init(newPosition, newSemiAxis, newSemiImagAxis);
 
  // introduce the shift in the parametrization
  // maybe the Sign must be changed?
  r->SetPShift(u);
  
  // set the bounds when necessary
  if (IsBounded()) 
    r->SetBounds(GetPStart(), GetPEnd());
  
  return r;
}

Tangent()

G4bool G4Hyperbola::Tangent ( G4CurvePoint &  cp, G4Vector3D &  v  )

virtual

Implements G4Curve.

Definition at line 184 of file G4Hyperbola.cc.

{
  // The tangent is computed from the 3D point representation
  // for all conics. An alternaive implementation (based on
  // the parametric point) might be worthwhile adding
  // for efficiency.
  
  const G4Axis2Placement3D& pos= *(GetPosition());
  G4Point3D p= pos.GetToPlacementCoordinates() * cp.GetPoint();
  
  v= forTangent*p.y()*pos.GetPX() + p.x()*pos.GetPY();
  
  return true;
}

---

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

• G4Hyperbola.hh
• G4Hyperbola.cc