#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 () const

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

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::G4QuadrangularFacet	(	const G4ThreeVector &	Pt0,
		const G4ThreeVector &	vt1,
		const G4ThreeVector &	vt2,
		const G4ThreeVector &	vt3,
		G4FacetVertexType	vertexType
	)

Definition at line 60 of file G4QuadrangularFacet.cc.

 {
   G4double delta   =  1.0 * kCarTolerance; // dimension tolerance
   G4double epsilon = 0.01 * kCarTolerance; // planarity tolerance
 
   fRadius = 0.0;
 
   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;
   }
 
   // Check length of sides and diagonals
   //
   G4double leng1 = e1.mag();
   G4double leng2 = (e2-e1).mag();
   G4double leng3 = (e3-e2).mag();
   G4double leng4 = e3.mag();
 
   G4double diag1 = e2.mag();
   G4double diag2 = (e3-e1).mag();
 
   if (leng1 <= delta || leng2 <= delta || leng3 <= delta || leng4 <= delta ||
       diag1 <= delta || diag2 <= delta)
   {
     ostringstream message;
     message << "Sides/diagonals of facet are too small." << G4endl
             << "P0 = " << GetVertex(0) << G4endl
             << "P1 = " << GetVertex(1) << G4endl
             << "P2 = " << GetVertex(2) << G4endl
             << "P3 = " << GetVertex(3) << G4endl
             << "Side1 length (P0->P1) = " << leng1 << G4endl
             << "Side2 length (P1->P2) = " << leng2 << G4endl
             << "Side3 length (P2->P3) = " << leng3 << G4endl
             << "Side4 length (P3->P0) = " << leng4 << G4endl
             << "Diagonal1 length (P0->P2) = " << diag1 << G4endl
             << "Diagonal2 length (P1->P3) = " << diag2;
     G4Exception("G4QuadrangularFacet::G4QuadrangularFacet()",
                 "GeomSolids1001", JustWarning, message);
     return;
   }
 
   // Check that vertices are not collinear
   //
   G4double s1 = (e1.cross(e2)).mag()*0.5;
   G4double s2 = ((e2-e1).cross(e3-e2)).mag()*0.5;
   G4double s3 = (e2.cross(e3)).mag()*0.5;
   G4double s4 = (e1.cross(e3)).mag()*0.5;
 
   G4double h1 = 2.*s1 / std::max(std::max(leng1,leng2),diag1);
   G4double h2 = 2.*s2 / std::max(std::max(leng2,leng3),diag2);
   G4double h3 = 2.*s3 / std::max(std::max(leng3,leng4),diag1);
   G4double h4 = 2.*s4 / std::max(std::max(leng4,leng1),diag2);
    
   if (h1 <= delta || h2 <= delta || h3 <= delta || h4 <= delta )
   {
     ostringstream message;
     message << "Facet has three or more collinear vertices." << G4endl
             << "P0 = " << GetVertex(0) << G4endl
             << "P1 = " << GetVertex(1) << G4endl
             << "P2 = " << GetVertex(2) << G4endl
             << "P3 = " << GetVertex(3) << G4endl
             << "Height in P0-P1-P2 = " << h1 << G4endl
             << "Height in P1-P2-P3 = " << h2 << G4endl
             << "Height in P2-P3-P4 = " << h3 << G4endl
             << "Height in P4-P0-P1 = " << h4;
     G4Exception("G4QuadrangularFacet::G4QuadrangularFacet()",
             "GeomSolids1001", JustWarning, message);
     return;
   }
 
   // Check that vertices are coplanar by computing minimal
   // height of tetrahedron comprising of vertices
   //
   G4double smax = std::max( std::max(s1,s2), std::max(s3,s4) ); 
   G4double hmin = 0.5 * std::fabs( e1.dot(e2.cross(e3)) ) / smax;
   if (hmin >= epsilon)
   {
     ostringstream message;
     message << "Facet is not planar." << G4endl
             << "Disrepancy = " << hmin << G4endl
             << "P0 = " << GetVertex(0) << G4endl
             << "P1 = " << GetVertex(1) << G4endl
             << "P2 = " << GetVertex(2) << G4endl
             << "P3 = " << GetVertex(3);
     G4Exception("G4QuadrangularFacet::G4QuadrangularFacet()",
             "GeomSolids1001", JustWarning, message);
     return;
   }
   
   // Check that facet is convex by computing crosspoint 
   // of diagonals
   //
   G4ThreeVector normal = e2.cross(e3-e1);
   G4double s = kInfinity, t = kInfinity, magnitude2 = normal.mag2();
   if (magnitude2 > delta*delta) // check: magnitude2 != 0.
   {
     s = normal.dot(e1.cross(e3-e1)) / magnitude2;
     t = normal.dot(e1.cross(e2))    / magnitude2;
   }
   if (s <= 0. || s >= 1. || t <= 0. || t >= 1.)
   {
     ostringstream message;
     message << "Facet is not convex." << G4endl
             << "Parameters of crosspoint of diagonals: "
             << s << " and " << t << G4endl
             << "should both be within (0,1) range" << G4endl
         << "P0 = " << GetVertex(0) << G4endl
         << "P1 = " << GetVertex(1) << G4endl
         << "P2 = " << GetVertex(2) << G4endl
         << "P3 = " << GetVertex(3);
     G4Exception("G4QuadrangularFacet::G4QuadrangularFacet()",
             "GeomSolids1001", JustWarning, message);
     return;
   }
 
   // Define facet
   //
   fFacet1 = G4TriangularFacet(GetVertex(0),GetVertex(1),GetVertex(2),ABSOLUTE);
   fFacet2 = G4TriangularFacet(GetVertex(0),GetVertex(2),GetVertex(3),ABSOLUTE);
 
   normal = normal.unit();
   fFacet1.SetSurfaceNormal(normal);
   fFacet2.SetSurfaceNormal(normal);
   
   G4ThreeVector vtmp = 0.5 * (e1 + e2);
   fCircumcentre = GetVertex(0) + vtmp;
   G4double radiusSqr = vtmp.mag2();
   fRadius = std::sqrt(radiusSqr);
   // 29.02.2016 Remark by E.Tcherniaev: computation
   // of fCircumcenter and fRadius is wrong, however
   // it did not create any problem till now.
   // Bizarre! Need to investigate!
 }

Here is the call graph for this function:

Here is the caller graph for this function:

G4QuadrangularFacet::G4QuadrangularFacet ( const G4QuadrangularFacet & right )

Definition at line 225 of file G4QuadrangularFacet.cc.

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

G4QuadrangularFacet::~G4QuadrangularFacet ( )

Definition at line 219 of file G4QuadrangularFacet.cc.

220 {

221 }

Member Function Documentation

G4ThreeVector G4QuadrangularFacet::Distance ( const G4ThreeVector & p )

Definition at line 260 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:

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

virtual

Implements G4VFacet.

Definition at line 271 of file G4QuadrangularFacet.cc.

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

Here is the call graph for this function:

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

virtual

Implements G4VFacet.

Definition at line 280 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:

G4double G4QuadrangularFacet::Extent ( const G4ThreeVector axis )

virtual

Implements G4VFacet.

Definition at line 297 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:

G4double G4QuadrangularFacet::GetArea ( ) const

virtual

Implements G4VFacet.

Definition at line 346 of file G4QuadrangularFacet.cc.

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

Here is the call graph for this function:

G4ThreeVector G4QuadrangularFacet::GetCircumcentre ( ) const

inlinevirtual

Implements G4VFacet.

Definition at line 142 of file G4QuadrangularFacet.hh.

 {
   return fCircumcentre;
 }

G4VFacet * G4QuadrangularFacet::GetClone ( )

virtual

Implements G4VFacet.

Definition at line 250 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:

G4String G4QuadrangularFacet::GetEntityType ( ) const

virtual

Implements G4VFacet.

Definition at line 354 of file G4QuadrangularFacet.cc.

 {
   return "G4QuadrangularFacet";
 }

G4int G4QuadrangularFacet::GetNumberOfVertices ( ) const

inlinevirtual

Implements G4VFacet.

Definition at line 126 of file G4QuadrangularFacet.hh.

 {
   return 4;
 }

G4ThreeVector G4QuadrangularFacet::GetPointOnFace ( ) const

virtual

Implements G4VFacet.

Definition at line 334 of file G4QuadrangularFacet.cc.

 {
   G4double s1 = fFacet1.GetArea();
   G4double s2 = fFacet2.GetArea();
   return ((s1+s2)*G4UniformRand() < s1) ?
     fFacet1.GetPointOnFace() : fFacet2.GetPointOnFace();
 }

Here is the call graph for this function:

G4double G4QuadrangularFacet::GetRadius ( ) const

inlinevirtual

Implements G4VFacet.

Definition at line 137 of file G4QuadrangularFacet.hh.

 {
   return fRadius;
 }

G4ThreeVector G4QuadrangularFacet::GetSurfaceNormal ( ) const

virtual

Implements G4VFacet.

Definition at line 361 of file G4QuadrangularFacet.cc.

 {
   return fFacet1.GetSurfaceNormal();
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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:

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

virtual

Implements G4VFacet.

Definition at line 311 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:

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:

G4QuadrangularFacet & G4QuadrangularFacet::operator= ( const G4QuadrangularFacet & right )

Definition at line 236 of file G4QuadrangularFacet.cc.

 {
   if (this == &rhs)
     return *this;
 
   fFacet1 = rhs.fFacet1;
   fFacet2 = rhs.fFacet2;
   fRadius = 0.0;
 
   return *this;
 }

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:

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:

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

geant4.10.03.p01/source/geometry/solids/specific/include/G4QuadrangularFacet.hh
geant4.10.03.p01/source/geometry/solids/specific/src/G4QuadrangularFacet.cc

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation