Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
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 ()
 
G4QuadrangularFacetoperator= (const G4QuadrangularFacet &right)
 
G4VFacetGetClone ()
 
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.

65 {
66  G4double delta = 1.0 * kCarTolerance; // dimension tolerance
67  G4double epsilon = 0.01 * kCarTolerance; // planarity tolerance
68 
69  fRadius = 0.0;
70 
71  G4ThreeVector e1, e2, e3;
72  SetVertex(0, vt0);
73  if (vertexType == ABSOLUTE)
74  {
75  SetVertex(1, vt1);
76  SetVertex(2, vt2);
77  SetVertex(3, vt3);
78 
79  e1 = vt1 - vt0;
80  e2 = vt2 - vt0;
81  e3 = vt3 - vt0;
82  }
83  else
84  {
85  SetVertex(1, vt0 + vt1);
86  SetVertex(2, vt0 + vt2);
87  SetVertex(3, vt0 + vt3);
88 
89  e1 = vt1;
90  e2 = vt2;
91  e3 = vt3;
92  }
93 
94  // Check length of sides and diagonals
95  //
96  G4double leng1 = e1.mag();
97  G4double leng2 = (e2-e1).mag();
98  G4double leng3 = (e3-e2).mag();
99  G4double leng4 = e3.mag();
100 
101  G4double diag1 = e2.mag();
102  G4double diag2 = (e3-e1).mag();
103 
104  if (leng1 <= delta || leng2 <= delta || leng3 <= delta || leng4 <= delta ||
105  diag1 <= delta || diag2 <= delta)
106  {
107  ostringstream message;
108  message << "Sides/diagonals of facet are too small." << G4endl
109  << "P0 = " << GetVertex(0) << G4endl
110  << "P1 = " << GetVertex(1) << G4endl
111  << "P2 = " << GetVertex(2) << G4endl
112  << "P3 = " << GetVertex(3) << G4endl
113  << "Side1 length (P0->P1) = " << leng1 << G4endl
114  << "Side2 length (P1->P2) = " << leng2 << G4endl
115  << "Side3 length (P2->P3) = " << leng3 << G4endl
116  << "Side4 length (P3->P0) = " << leng4 << G4endl
117  << "Diagonal1 length (P0->P2) = " << diag1 << G4endl
118  << "Diagonal2 length (P1->P3) = " << diag2;
119  G4Exception("G4QuadrangularFacet::G4QuadrangularFacet()",
120  "GeomSolids1001", JustWarning, message);
121  return;
122  }
123 
124  // Check that vertices are not collinear
125  //
126  G4double s1 = (e1.cross(e2)).mag()*0.5;
127  G4double s2 = ((e2-e1).cross(e3-e2)).mag()*0.5;
128  G4double s3 = (e2.cross(e3)).mag()*0.5;
129  G4double s4 = (e1.cross(e3)).mag()*0.5;
130 
131  G4double h1 = 2.*s1 / std::max(std::max(leng1,leng2),diag1);
132  G4double h2 = 2.*s2 / std::max(std::max(leng2,leng3),diag2);
133  G4double h3 = 2.*s3 / std::max(std::max(leng3,leng4),diag1);
134  G4double h4 = 2.*s4 / std::max(std::max(leng4,leng1),diag2);
135 
136  if (h1 <= delta || h2 <= delta || h3 <= delta || h4 <= delta )
137  {
138  ostringstream message;
139  message << "Facet has three or more collinear vertices." << G4endl
140  << "P0 = " << GetVertex(0) << G4endl
141  << "P1 = " << GetVertex(1) << G4endl
142  << "P2 = " << GetVertex(2) << G4endl
143  << "P3 = " << GetVertex(3) << G4endl
144  << "Height in P0-P1-P2 = " << h1 << G4endl
145  << "Height in P1-P2-P3 = " << h2 << G4endl
146  << "Height in P2-P3-P4 = " << h3 << G4endl
147  << "Height in P4-P0-P1 = " << h4;
148  G4Exception("G4QuadrangularFacet::G4QuadrangularFacet()",
149  "GeomSolids1001", JustWarning, message);
150  return;
151  }
152 
153  // Check that vertices are coplanar by computing minimal
154  // height of tetrahedron comprising of vertices
155  //
156  G4double smax = std::max( std::max(s1,s2), std::max(s3,s4) );
157  G4double hmin = 0.5 * std::fabs( e1.dot(e2.cross(e3)) ) / smax;
158  if (hmin >= epsilon)
159  {
160  ostringstream message;
161  message << "Facet is not planar." << G4endl
162  << "Disrepancy = " << hmin << G4endl
163  << "P0 = " << GetVertex(0) << G4endl
164  << "P1 = " << GetVertex(1) << G4endl
165  << "P2 = " << GetVertex(2) << G4endl
166  << "P3 = " << GetVertex(3);
167  G4Exception("G4QuadrangularFacet::G4QuadrangularFacet()",
168  "GeomSolids1001", JustWarning, message);
169  return;
170  }
171 
172  // Check that facet is convex by computing crosspoint
173  // of diagonals
174  //
175  G4ThreeVector normal = e2.cross(e3-e1);
176  G4double s = kInfinity, t = kInfinity, magnitude2 = normal.mag2();
177  if (magnitude2 > delta*delta) // check: magnitude2 != 0.
178  {
179  s = normal.dot(e1.cross(e3-e1)) / magnitude2;
180  t = normal.dot(e1.cross(e2)) / magnitude2;
181  }
182  if (s <= 0. || s >= 1. || t <= 0. || t >= 1.)
183  {
184  ostringstream message;
185  message << "Facet is not convex." << G4endl
186  << "Parameters of crosspoint of diagonals: "
187  << s << " and " << t << G4endl
188  << "should both be within (0,1) range" << G4endl
189  << "P0 = " << GetVertex(0) << G4endl
190  << "P1 = " << GetVertex(1) << G4endl
191  << "P2 = " << GetVertex(2) << G4endl
192  << "P3 = " << GetVertex(3);
193  G4Exception("G4QuadrangularFacet::G4QuadrangularFacet()",
194  "GeomSolids1001", JustWarning, message);
195  return;
196  }
197 
198  // Define facet
199  //
202 
203  normal = normal.unit();
204  fFacet1.SetSurfaceNormal(normal);
205  fFacet2.SetSurfaceNormal(normal);
206 
207  G4ThreeVector vtmp = 0.5 * (e1 + e2);
208  fCircumcentre = GetVertex(0) + vtmp;
209  G4double radiusSqr = vtmp.mag2();
210  fRadius = std::sqrt(radiusSqr);
211  // 29.02.2016 Remark by E.Tcherniaev: computation
212  // of fCircumcenter and fRadius is wrong, however
213  // it did not create any problem till now.
214  // Bizarre! Need to investigate!
215 }
static const G4double kInfinity
Definition: geomdefs.hh:42
static const G4double kCarTolerance
Definition: G4VFacet.hh:102
double dot(const Hep3Vector &) const
const G4int smax
static double normal(HepRandomEngine *eptr)
Definition: RandPoisson.cc:77
const XML_Char * s
Definition: expat.h:262
G4ThreeVector GetVertex(G4int i) const
void SetSurfaceNormal(G4ThreeVector normal)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
Definition: G4Exception.cc:41
void SetVertex(G4int i, const G4ThreeVector &val)
T max(const T t1, const T t2)
brief Return the largest of the two arguments
Hep3Vector unit() const
double mag2() const
#define G4endl
Definition: G4ios.hh:61
Hep3Vector cross(const Hep3Vector &) const
double G4double
Definition: G4Types.hh:76
double mag() const
double epsilon(double density, double temperature)

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.

226  : G4VFacet(rhs)
227 {
228  fFacet1 = rhs.fFacet1;
229  fFacet2 = rhs.fFacet2;
230  fRadius = 0.0;
231 }
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.

261 {
262  G4ThreeVector v1 = fFacet1.Distance(p);
263  G4ThreeVector v2 = fFacet2.Distance(p);
264 
265  if (v1.mag2() < v2.mag2()) return v1;
266  else return v2;
267 }
double mag2() const
G4ThreeVector Distance(const G4ThreeVector &p)

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.

273 {
274  G4double dist = Distance(p).mag();
275  return dist;
276 }
double G4double
Definition: G4Types.hh:76
double mag() const
G4ThreeVector Distance(const G4ThreeVector &p)

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.

282 {
283  G4double dist;
284 
285  G4ThreeVector v = Distance(p);
286  G4double dir = v.dot(GetSurfaceNormal());
287  if ( ((dir > dirTolerance) && (!outgoing))
288  || ((dir < -dirTolerance) && outgoing))
289  dist = kInfinity;
290  else
291  dist = v.mag();
292  return dist;
293 }
static const G4double kInfinity
Definition: geomdefs.hh:42
double dot(const Hep3Vector &) const
G4ThreeVector GetSurfaceNormal() const
static const G4double dirTolerance
Definition: G4VFacet.hh:101
tuple v
Definition: test.py:18
double G4double
Definition: G4Types.hh:76
double mag() const
G4ThreeVector Distance(const G4ThreeVector &p)

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.

298 {
299  G4double ss = 0;
300 
301  for (G4int i = 0; i <= 3; ++i)
302  {
303  G4double sp = GetVertex(i).dot(axis);
304  if (sp > ss) ss = sp;
305  }
306  return ss;
307 }
double dot(const Hep3Vector &) const
int G4int
Definition: G4Types.hh:78
G4ThreeVector GetVertex(G4int i) const
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4double G4QuadrangularFacet::GetArea ( ) const
virtual

Implements G4VFacet.

Definition at line 346 of file G4QuadrangularFacet.cc.

347 {
348  G4double area = fFacet1.GetArea() + fFacet2.GetArea();
349  return area;
350 }
G4double GetArea() const
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4ThreeVector G4QuadrangularFacet::GetCircumcentre ( ) const
inlinevirtual

Implements G4VFacet.

Definition at line 142 of file G4QuadrangularFacet.hh.

143 {
144  return fCircumcentre;
145 }
G4VFacet * G4QuadrangularFacet::GetClone ( )
virtual

Implements G4VFacet.

Definition at line 250 of file G4QuadrangularFacet.cc.

251 {
253  GetVertex(2), GetVertex(3),
254  ABSOLUTE);
255  return c;
256 }
G4QuadrangularFacet(const G4ThreeVector &Pt0, const G4ThreeVector &vt1, const G4ThreeVector &vt2, const G4ThreeVector &vt3, G4FacetVertexType)
G4ThreeVector GetVertex(G4int i) const
tuple c
Definition: test.py:13

Here is the call graph for this function:

G4String G4QuadrangularFacet::GetEntityType ( ) const
virtual

Implements G4VFacet.

Definition at line 354 of file G4QuadrangularFacet.cc.

355 {
356  return "G4QuadrangularFacet";
357 }
G4int G4QuadrangularFacet::GetNumberOfVertices ( ) const
inlinevirtual

Implements G4VFacet.

Definition at line 126 of file G4QuadrangularFacet.hh.

127 {
128  return 4;
129 }
G4ThreeVector G4QuadrangularFacet::GetPointOnFace ( ) const
virtual

Implements G4VFacet.

Definition at line 334 of file G4QuadrangularFacet.cc.

335 {
336  G4double s1 = fFacet1.GetArea();
337  G4double s2 = fFacet2.GetArea();
338  return ((s1+s2)*G4UniformRand() < s1) ?
339  fFacet1.GetPointOnFace() : fFacet2.GetPointOnFace();
340 }
G4double GetArea() const
#define G4UniformRand()
Definition: Randomize.hh:97
double G4double
Definition: G4Types.hh:76
G4ThreeVector GetPointOnFace() const

Here is the call graph for this function:

G4double G4QuadrangularFacet::GetRadius ( ) const
inlinevirtual

Implements G4VFacet.

Definition at line 137 of file G4QuadrangularFacet.hh.

138 {
139  return fRadius;
140 }
G4ThreeVector G4QuadrangularFacet::GetSurfaceNormal ( ) const
virtual

Implements G4VFacet.

Definition at line 361 of file G4QuadrangularFacet.cc.

362 {
363  return fFacet1.GetSurfaceNormal();
364 }
G4ThreeVector GetSurfaceNormal() const

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.

132 {
133  return i == 3 ? fFacet2.GetVertex(2) : fFacet1.GetVertex(i);
134 }
G4ThreeVector GetVertex(G4int i) const

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.

317 {
318  G4bool intersect =
319  fFacet1.Intersect(p,v,outgoing,distance,distFromSurface,normal);
320  if (!intersect) intersect =
321  fFacet2.Intersect(p,v,outgoing,distance,distFromSurface,normal);
322  if (!intersect)
323  {
324  distance = distFromSurface = kInfinity;
325  normal.set(0,0,0);
326  }
327  return intersect;
328 }
void set(double x, double y, double z)
G4bool Intersect(const G4ThreeVector &p, const G4ThreeVector &v, const G4bool outgoing, G4double &distance, G4double &distFromSurface, G4ThreeVector &normal)
static const G4double kInfinity
Definition: geomdefs.hh:42
bool G4bool
Definition: G4Types.hh:79

Here is the call graph for this function:

G4bool G4QuadrangularFacet::IsDefined ( ) const
inlinevirtual

Implements G4VFacet.

Definition at line 174 of file G4QuadrangularFacet.hh.

175 {
176  return fFacet1.IsDefined();
177 }
G4bool IsDefined() const

Here is the call graph for this function:

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

Definition at line 236 of file G4QuadrangularFacet.cc.

237 {
238  if (this == &rhs)
239  return *this;
240 
241  fFacet1 = rhs.fFacet1;
242  fFacet2 = rhs.fFacet2;
243  fRadius = 0.0;
244 
245  return *this;
246 }
void G4QuadrangularFacet::SetVertex ( G4int  i,
const G4ThreeVector val 
)
inlinevirtual

Implements G4VFacet.

Definition at line 147 of file G4QuadrangularFacet.hh.

148 {
149  switch (i)
150  {
151  case 0:
152  fFacet1.SetVertex(0, val);
153  fFacet2.SetVertex(0, val);
154  break;
155  case 1:
156  fFacet1.SetVertex(1, val);
157  break;
158  case 2:
159  fFacet1.SetVertex(2, val);
160  fFacet2.SetVertex(1, val);
161  break;
162  case 3:
163  fFacet2.SetVertex(2, val);
164  break;
165  }
166 }
void SetVertex(G4int i, const G4ThreeVector &val)

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.

169 {
170  fFacet1.SetVertices(v);
171  fFacet2.SetVertices(v);
172 }
void SetVertices(std::vector< G4ThreeVector > *v)
tuple v
Definition: test.py:18

Here is the call graph for this function:


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