G4Ellipse Class Reference

#include <G4Ellipse.hh>

Inheritance diagram for G4Ellipse:

Public Member Functions
	G4Ellipse ()
virtual	~G4Ellipse ()
	G4Ellipse (const G4Ellipse &right)
G4Ellipse &	operator= (const G4Ellipse &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	GetSemiAxis1 () const
G4double	GetSemiAxis2 () const
void	Init (const G4Axis2Placement3D &position0, G4double semiAxis10, G4double semiAxis20)
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 46 of file G4Ellipse.hh.

Constructor & Destructor Documentation

G4Ellipse() [1/2]

G4Ellipse::G4Ellipse

(

)

Definition at line 42 of file G4Ellipse.cc.

  : semiAxis1(0.), semiAxis2(0.), ratioAxis2Axis1(0.), forTangent(.0)
{
}

Referenced by Project().

~G4Ellipse()

G4Ellipse::~G4Ellipse ( )

virtual

Definition at line 47 of file G4Ellipse.cc.

{
}

G4Ellipse() [2/2]

G4Ellipse::G4Ellipse

(

const G4Ellipse &

right

)

Definition at line 51 of file G4Ellipse.cc.

  : G4Conic(), semiAxis1(right.semiAxis1), semiAxis2(right.semiAxis2),
    ratioAxis2Axis1(right.ratioAxis2Axis1), toUnitCircle(right.toUnitCircle),
    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 G4Ellipse::GetPMax ( ) const

inlinevirtual

Implements G4Curve.

GetPoint()

G4Point3D G4Ellipse::GetPoint ( G4double  param ) const

inlinevirtual

Implements G4Curve.

Referenced by InitBounded().

GetPPoint()

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

inlinevirtual

Implements G4Curve.

GetSemiAxis1()

G4double G4Ellipse::GetSemiAxis1 ( ) const

inline

GetSemiAxis2()

G4double G4Ellipse::GetSemiAxis2 ( ) const

inline

Init()

void G4Ellipse::Init ( const G4Axis2Placement3D &  position0, G4double  semiAxis10, G4double  semiAxis20  )

inline

Referenced by G4CircularCurve::Project(), and Project().

InitBounded()

void G4Ellipse::InitBounded ( )

protectedvirtual

Implements G4Curve.

Definition at line 155 of file G4Ellipse.cc.

{
  // original implementation
  // const G4Point3D& center = position.GetLocation();
  // G4double maxEntent      = std::max(semiAxis1, semiAxis2);
  // G4Vector3D halfExtent(maxEntent, maxEntent, maxEntent);
  // bBox.Init(center+halfExtent, center-halfExtent);
 
  // 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 u= std::atan2(position.GetPY()(i)*semiAxis2,
              position.GetPX()(i)*semiAxis1);
    if (IsPOn(u)) 
      bBox.Extend(GetPoint(u));
    
    if (IsPOn(u+pi)) 
      bBox.Extend(GetPoint(u+pi));
  }
}

IntersectRay2D()

G4int G4Ellipse::IntersectRay2D ( const G4Ray &  ray )

inlinevirtual

Implements G4Curve.

operator=()

G4Ellipse & G4Ellipse::operator=

(

const G4Ellipse &

right

)

Definition at line 68 of file G4Ellipse.cc.

{
  if (&right == this) return *this;
 
  semiAxis1 = right.semiAxis1;
  semiAxis2 = right.semiAxis2;
  ratioAxis2Axis1 = right.ratioAxis2Axis1;
  toUnitCircle    = right.toUnitCircle;
  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 * G4Ellipse::Project ( const G4Transform3D &  tr = G4Transform3D::Identity )

virtual

Implements G4Curve.

Definition at line 91 of file G4Ellipse.cc.

{
  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( semiAxis1*xPrime );
  G4Vector3D b = G4Vector3D( semiAxis2*yPrime );
  
  G4double u;
  G4double abmag = a.mag2()-b.mag2();
  G4double prod = 2*a*b;
 
  if ((abmag > FLT_MAX) && (prod < -FLT_MAX))
    u = -pi/8;
  else if ((abmag < -FLT_MAX) && (prod > FLT_MAX))
    u = 3*pi/8;
  else if ((abmag < -FLT_MAX) && (prod < -FLT_MAX))
    u = -3*pi/8;
  else if ((std::abs(abmag) < perMillion) && (std::abs(prod) < perMillion))
    u = 0.;
  else
    u = std::atan2(prod,abmag) / 2;
 
  // get the coordinate axis directions and the semiaxis lengths
  G4Vector3D sAxis1          = G4Vector3D( a*std::cos(u)+b*std::sin(u) );
  G4Vector3D sAxis2          = G4Vector3D( a*std::cos(u+pi/2)+b*std::sin(u+pi/2) );
  G4double newSemiAxis1      = sAxis1.mag();
  G4double newSemiAxis2      = sAxis2.mag();
  G4Vector3D newRefDirection = sAxis1;
 
  // create the new ellipse
  G4Axis2Placement3D newPosition;
  newPosition.Init(newRefDirection, newAxis, newLocation);
  G4Ellipse* r= new G4Ellipse;
  r->Init(newPosition, newSemiAxis1, newSemiAxis2);
  
  // 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());
  
  // L. Broglia
  // copy sense of the curve
  r->SetSameSense(GetSameSense());
 
  return r;
}

Referenced by G4CircularCurve::Project().

Tangent()

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

virtual

Implements G4Curve.

Definition at line 182 of file G4Ellipse.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();
  if(GetSameSense())
    v = -v;
  
  return true;
}

---

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

• G4Ellipse.hh
• G4Ellipse.cc