G4QuadrangularFacet Class Reference

#include <G4QuadrangularFacet.hh>

Inheritance diagram for G4QuadrangularFacet:

Collaboration diagram for G4QuadrangularFacet:

Public Member Functions
	G4QuadrangularFacet (const G4ThreeVector &Pt0, const G4ThreeVector &vt1, const G4ThreeVector &vt2, const G4ThreeVector &vt3, G4FacetVertexType)
	G4QuadrangularFacet (const G4QuadrangularFacet &right)
	~G4QuadrangularFacet ()
G4QuadrangularFacet &	operator= (const G4QuadrangularFacet &right)
G4VFacet *	GetClone ()
G4ThreeVector	Distance (const G4ThreeVector &p)
G4double	Distance (const G4ThreeVector &p, G4double minDist)
G4double	Distance (const G4ThreeVector &p, G4double minDist, const G4bool outgoing)
G4double	Extent (const G4ThreeVector axis)
G4bool	Intersect (const G4ThreeVector &p, const G4ThreeVector &v, const G4bool outgoing, G4double &distance, G4double &distFromSurface, G4ThreeVector &normal)
G4ThreeVector	GetSurfaceNormal () const
G4double	GetArea ()
G4ThreeVector	GetPointOnFace () const
G4GeometryType	GetEntityType () const
G4bool	IsDefined () const
G4int	GetNumberOfVertices () const
G4ThreeVector	GetVertex (G4int i) const
void	SetVertex (G4int i, const G4ThreeVector &val)
void	SetVertices (std::vector< G4ThreeVector > *v)
G4double	GetRadius () const
G4ThreeVector	GetCircumcentre () const
Public Member Functions inherited from G4VFacet
virtual	~G4VFacet ()
G4bool	operator== (const G4VFacet &right) const
void	ApplyTranslation (const G4ThreeVector v)
std::ostream &	StreamInfo (std::ostream &os) const
G4bool	IsInside (const G4ThreeVector &p) const

Private Member Functions
G4int	GetVertexIndex (G4int i) const
void	SetVertexIndex (G4int i, G4int val)
G4int	AllocatedMemory ()

Private Attributes
G4double	fRadius
G4ThreeVector	fCircumcentre
G4TriangularFacet	fFacet1
G4TriangularFacet	fFacet2

Additional Inherited Members
Static Protected Attributes inherited from G4VFacet
static const G4double	dirTolerance = 1.0E-14
static const G4double	kCarTolerance

Detailed Description

Definition at line 69 of file G4QuadrangularFacet.hh.

Constructor & Destructor Documentation

G4QuadrangularFacet() [1/2]

G4QuadrangularFacet::G4QuadrangularFacet	(	const G4ThreeVector &	Pt0,
		const G4ThreeVector &	vt1,
		const G4ThreeVector &	vt2,
		const G4ThreeVector &	vt3,
		G4FacetVertexType	vertexType
	)

Definition at line 54 of file G4QuadrangularFacet.cc.

{
  G4ThreeVector e1, e2, e3;

  SetVertex(0, vt0);
  if (vertexType == ABSOLUTE)
  {
    SetVertex(1, vt1);
    SetVertex(2, vt2);
    SetVertex(3, vt3);

    e1 = vt1 - vt0;
    e2 = vt2 - vt0;
    e3 = vt3 - vt0;
  }
  else
  {
    SetVertex(1, vt0 + vt1);
    SetVertex(2, vt0 + vt2);
    SetVertex(3, vt0 + vt3);

    e1 = vt1;
    e2 = vt2;
    e3 = vt3;
  }
  G4double length1 = e1.mag();
  G4double length2 = (GetVertex(2)-GetVertex(1)).mag();
  G4double length3 = (GetVertex(3)-GetVertex(2)).mag();
  G4double length4 = e3.mag();

  G4ThreeVector normal1 = e1.cross(e2).unit();
  G4ThreeVector normal2 = e2.cross(e3).unit(); 

  bool isDefined = (length1 > kCarTolerance && length2 > kCarTolerance &&
    length3 > kCarTolerance && length4 > kCarTolerance &&
    normal1.dot(normal2) >= 0.9999999999);

  if (isDefined)
  {
    fFacet1 = G4TriangularFacet (GetVertex(0),GetVertex(1),
                                 GetVertex(2),ABSOLUTE);
    fFacet2 = G4TriangularFacet (GetVertex(0),GetVertex(2),
                                 GetVertex(3),ABSOLUTE);

    G4TriangularFacet facet3 (GetVertex(0),GetVertex(1),GetVertex(3),ABSOLUTE);
    G4TriangularFacet facet4 (GetVertex(1),GetVertex(2),GetVertex(3),ABSOLUTE);

    G4ThreeVector normal12 = fFacet1.GetSurfaceNormal()
                           + fFacet2.GetSurfaceNormal();
    G4ThreeVector normal34 = facet3.GetSurfaceNormal()
                           + facet4.GetSurfaceNormal();
    G4ThreeVector normal = 0.25 * (normal12 + normal34);

    fFacet1.SetSurfaceNormal (normal);
    fFacet2.SetSurfaceNormal (normal);

    G4ThreeVector vtmp = 0.5 * (e1 + e2);
    fCircumcentre = GetVertex(0) + vtmp;
    G4double radiusSqr = vtmp.mag2();
    fRadius = std::sqrt(radiusSqr);
  }
  else
  {
    G4Exception("G4QuadrangularFacet::G4QuadrangularFacet()",
                "GeomSolids0002", JustWarning,
                "Length of sides of facet are too small or sides not planar.");
    G4cout << G4endl;
    G4cout << "P0 = " << GetVertex(0) << G4endl;
    G4cout << "P1 = " << GetVertex(1) << G4endl;
    G4cout << "P2 = " << GetVertex(2) << G4endl;
    G4cout << "P3 = " << GetVertex(3) << G4endl;
    G4cout << "Side lengths = P0->P1" << length1 << G4endl;    
    G4cout << "Side lengths = P1->P2" << length2 << G4endl;    
    G4cout << "Side lengths = P2->P3" << length3 << G4endl;    
    G4cout << "Side lengths = P3->P0" << length4 << G4endl;    
    G4cout << G4endl;
    fRadius = 0.0;
  }
}

Here is the call graph for this function:

Here is the caller graph for this function:

G4QuadrangularFacet() [2/2]

G4QuadrangularFacet::G4QuadrangularFacet

(

const G4QuadrangularFacet &

right

)

Definition at line 146 of file G4QuadrangularFacet.cc.

  : G4VFacet(rhs)
{
  fFacet1 = rhs.fFacet1;
  fFacet2 = rhs.fFacet2;
  fRadius = 0.0;
}

~G4QuadrangularFacet()

G4QuadrangularFacet::~G4QuadrangularFacet

(

)

Definition at line 140 of file G4QuadrangularFacet.cc.

{
}

Member Function Documentation

AllocatedMemory()

G4int G4QuadrangularFacet::AllocatedMemory ( )

inlineprivatevirtual

Implements G4VFacet.

Definition at line 206 of file G4QuadrangularFacet.hh.

{
  return sizeof(*this) + fFacet1.AllocatedMemory() + fFacet2.AllocatedMemory();
}

Here is the call graph for this function:

Distance() [1/3]

G4ThreeVector G4QuadrangularFacet::Distance

(

const G4ThreeVector &

p

)

Definition at line 181 of file G4QuadrangularFacet.cc.

{
  G4ThreeVector v1 = fFacet1.Distance(p);
  G4ThreeVector v2 = fFacet2.Distance(p);

  if (v1.mag2() < v2.mag2()) return v1;
  else return v2;
}

Here is the call graph for this function:

Here is the caller graph for this function:

Distance() [2/3]

G4double G4QuadrangularFacet::Distance ( const G4ThreeVector &  p, G4double  minDist  )

virtual

Implements G4VFacet.

Definition at line 192 of file G4QuadrangularFacet.cc.

{  
  G4double dist = Distance(p).mag();
  return dist;
}

Here is the call graph for this function:

Distance() [3/3]

G4double G4QuadrangularFacet::Distance ( const G4ThreeVector &  p, G4double  minDist, const G4bool  outgoing  )

virtual

Implements G4VFacet.

Definition at line 201 of file G4QuadrangularFacet.cc.

{
  G4double dist;

  G4ThreeVector v = Distance(p);
  G4double dir = v.dot(GetSurfaceNormal());
  if ( ((dir > dirTolerance) && (!outgoing))
    || ((dir < -dirTolerance) && outgoing))
    dist = kInfinity;
  else 
    dist = v.mag();
  return dist;
}

Here is the call graph for this function:

Extent()

G4double G4QuadrangularFacet::Extent ( const G4ThreeVector  axis )

virtual

Implements G4VFacet.

Definition at line 218 of file G4QuadrangularFacet.cc.

{
  G4double ss  = 0;

  for (G4int i = 0; i <= 3; ++i)
  {
    G4double sp = GetVertex(i).dot(axis);
    if (sp > ss) ss = sp;
  }
  return ss;
}

Here is the call graph for this function:

GetArea()

G4double G4QuadrangularFacet::GetArea ( )

virtual

Implements G4VFacet.

Definition at line 266 of file G4QuadrangularFacet.cc.

{
  G4double area = fFacet1.GetArea() + fFacet2.GetArea();
  return area;
}

Here is the call graph for this function:

GetCircumcentre()

G4ThreeVector G4QuadrangularFacet::GetCircumcentre ( ) const

inlinevirtual

Implements G4VFacet.

Definition at line 142 of file G4QuadrangularFacet.hh.

{
  return fCircumcentre;
}

GetClone()

G4VFacet * G4QuadrangularFacet::GetClone ( )

virtual

Implements G4VFacet.

Definition at line 171 of file G4QuadrangularFacet.cc.

{
  G4QuadrangularFacet *c = new G4QuadrangularFacet (GetVertex(0), GetVertex(1),
                                                    GetVertex(2), GetVertex(3),
                                                    ABSOLUTE);
  return c;
}

Here is the call graph for this function:

GetEntityType()

G4String G4QuadrangularFacet::GetEntityType ( ) const

virtual

Implements G4VFacet.

Definition at line 274 of file G4QuadrangularFacet.cc.

{
  return "G4QuadrangularFacet";
}

GetNumberOfVertices()

G4int G4QuadrangularFacet::GetNumberOfVertices ( ) const

inlinevirtual

Implements G4VFacet.

Definition at line 126 of file G4QuadrangularFacet.hh.

{
  return 4;
}

GetPointOnFace()

G4ThreeVector G4QuadrangularFacet::GetPointOnFace ( ) const

virtual

Implements G4VFacet.

Definition at line 255 of file G4QuadrangularFacet.cc.

{
  G4ThreeVector pr = (G4RandFlat::shoot(0.,1.) < 0.5)
                   ? fFacet1.GetPointOnFace() : fFacet2.GetPointOnFace();
  return pr;
}

Here is the call graph for this function:

GetRadius()

G4double G4QuadrangularFacet::GetRadius ( ) const

inlinevirtual

Implements G4VFacet.

Definition at line 137 of file G4QuadrangularFacet.hh.

{
  return fRadius;
}

GetSurfaceNormal()

G4ThreeVector G4QuadrangularFacet::GetSurfaceNormal ( ) const

virtual

Implements G4VFacet.

Definition at line 281 of file G4QuadrangularFacet.cc.

{
  return fFacet1.GetSurfaceNormal();
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetVertex()

G4ThreeVector G4QuadrangularFacet::GetVertex ( G4int  i ) const

inlinevirtual

Implements G4VFacet.

Definition at line 131 of file G4QuadrangularFacet.hh.

{
  return i == 3 ? fFacet2.GetVertex(2) : fFacet1.GetVertex(i);
}

Here is the call graph for this function:

Here is the caller graph for this function:

GetVertexIndex()

G4int G4QuadrangularFacet::GetVertexIndex ( G4int  i ) const

inlineprivatevirtual

Implements G4VFacet.

Definition at line 179 of file G4QuadrangularFacet.hh.

{
  return i == 3 ? fFacet2.GetVertexIndex(2) : fFacet1.GetVertexIndex(i);
}

Here is the call graph for this function:

Intersect()

G4bool G4QuadrangularFacet::Intersect ( const G4ThreeVector &  p, const G4ThreeVector &  v, const G4bool  outgoing, G4double &  distance, G4double &  distFromSurface, G4ThreeVector &  normal  )

virtual

Implements G4VFacet.

Definition at line 232 of file G4QuadrangularFacet.cc.

{
  G4bool intersect =
    fFacet1.Intersect(p,v,outgoing,distance,distFromSurface,normal);
  if (!intersect) intersect =
    fFacet2.Intersect(p,v,outgoing,distance,distFromSurface,normal);
  if (!intersect)
  {
    distance = distFromSurface = kInfinity;
    normal.set(0,0,0);
  }
  return intersect;
}

Here is the call graph for this function:

IsDefined()

G4bool G4QuadrangularFacet::IsDefined ( ) const

inlinevirtual

Implements G4VFacet.

Definition at line 174 of file G4QuadrangularFacet.hh.

{
  return fFacet1.IsDefined();
}

Here is the call graph for this function:

operator=()

G4QuadrangularFacet & G4QuadrangularFacet::operator=

(

const G4QuadrangularFacet &

right

)

Definition at line 157 of file G4QuadrangularFacet.cc.

{
  if (this == &rhs)
    return *this;

  fFacet1 = rhs.fFacet1;
  fFacet2 = rhs.fFacet2;
  fRadius = 0.0;

  return *this;
}

SetVertex()

void G4QuadrangularFacet::SetVertex ( G4int  i, const G4ThreeVector &  val  )

inlinevirtual

Implements G4VFacet.

Definition at line 147 of file G4QuadrangularFacet.hh.

{
  switch (i)
  {
    case 0:
      fFacet1.SetVertex(0, val);
      fFacet2.SetVertex(0, val);
      break;
    case 1:
      fFacet1.SetVertex(1, val);
      break;
    case 2:
      fFacet1.SetVertex(2, val);
      fFacet2.SetVertex(1, val);
      break;
    case 3:
      fFacet2.SetVertex(2, val);
      break;
    }
}

Here is the call graph for this function:

SetVertexIndex()

void G4QuadrangularFacet::SetVertexIndex ( G4int  i, G4int  val  )

inlineprivatevirtual

Implements G4VFacet.

Definition at line 185 of file G4QuadrangularFacet.hh.

{
  switch (i)
  {
    case 0:
      fFacet1.SetVertexIndex(0, val);
      fFacet2.SetVertexIndex(0, val);
      break;
    case 1:
      fFacet1.SetVertexIndex(1, val);
      break;
    case 2:
      fFacet1.SetVertexIndex(2, val);
      fFacet2.SetVertexIndex(1, val);
      break;
    case 3:
      fFacet2.SetVertexIndex(2, val);
      break;
  }
}

Here is the call graph for this function:

SetVertices()

void G4QuadrangularFacet::SetVertices ( std::vector< G4ThreeVector > *  v )

inlinevirtual

Implements G4VFacet.

Definition at line 168 of file G4QuadrangularFacet.hh.

{
  fFacet1.SetVertices(v);
  fFacet2.SetVertices(v);
}

Here is the call graph for this function:

Member Data Documentation

fCircumcentre

G4ThreeVector G4QuadrangularFacet::fCircumcentre

private

Definition at line 117 of file G4QuadrangularFacet.hh.

fFacet1

G4TriangularFacet G4QuadrangularFacet::fFacet1

private

Definition at line 119 of file G4QuadrangularFacet.hh.

fFacet2

G4TriangularFacet G4QuadrangularFacet::fFacet2

private

Definition at line 119 of file G4QuadrangularFacet.hh.

fRadius

G4double G4QuadrangularFacet::fRadius

private

Definition at line 116 of file G4QuadrangularFacet.hh.

---

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

• G4QuadrangularFacet.hh
• G4QuadrangularFacet.cc