84 : ncorners(0), corners(0)
97 r[0] = tail->
r;
z[0] = tail->
z;
98 r[1] = head->
r;
z[1] = head->
z;
164 zNormEdge[0] /= lAdj;
176 zNormEdge[1] /= lAdj;
185 : startPhi(0.), deltaPhi(0.), phiIsOpen(false), allBehind(false),
186 cone(0), rNorm(0.), zNorm(0.), rS(0.), zS(0.), length(0.),
187 prevRS(0.), prevZS(0.), nextRS(0.), nextZS(0.), ncorners(0), corners(0),
188 kCarTolerance(0.), fSurfaceArea(0.), instanceID(0)
224 if (
this == &source) {
return *
this; }
296 G4double normSign = outgoing ? +1 : -1;
304 if (nside == 0)
return false;
316 if (normSign*v.dot(normal) > 0)
333 if (normSign*v.dot(pNormal) > 0)
339 distFromSurface = -normSign*
DistanceAway( p,
false, distOutside2 );
340 if (distOutside2 < surfTolerance*surfTolerance)
342 if (distFromSurface > -surfTolerance)
354 distFromSurface = s1;
367 if (nside==1)
return false;
379 if (normSign*v.dot(normal) > 0)
384 if (normSign*v.dot(pNormal) > 0)
387 distFromSurface = -normSign*
DistanceAway( p,
false, distOutside2 );
388 if (distOutside2 < surfTolerance*surfTolerance)
390 if (distFromSurface > -surfTolerance)
398 distFromSurface = s2;
420 G4double normSign = outgoing ? -1 : +1;
433 return std::sqrt( distFrom*distFrom + distOut2 );
435 return std::fabs(distFrom);
446 return std::sqrt( distFrom*distFrom + distOut2 );
448 return std::fabs(distFrom);
466 dist2 = distFrom*distFrom + distOut2;
468 *bestDistance = std::sqrt( dist2);
472 if ( (std::fabs(edgeRZnorm) < tolerance)
473 && (distOut2< tolerance*tolerance) )
475 else if (edgeRZnorm < 0)
494 *bestDistance = std::sqrt( dFrom*dFrom + dOut2 );
542 if (bd > ad) ad = bd;
543 if (
cd > ad) ad =
cd;
544 if (dd > ad) ad = dd;
599 G4double rFudge = 1.0/std::cos(0.5*sigPhi);
729 cosPhi = std::cos(phi),
730 sinPhi = std::sin(phi);
733 v1( r1*cosPhi, r1*sinPhi, z1 ),
748 cosPhi = std::cos(phi),
749 sinPhi = std::sin(phi);
777 polygon.
SetNormal( deltaV.cross(v1-v0).unit() );
799 polygon.
SetNormal( deltaV.cross(v0-v2).unit() );
812 }
while( --numPhi > 0 );
826 a1(
r[1]*cosPhi,
r[1]*sinPhi,
z[1] ),
827 b0( r0*cosPhi, r0*sinPhi,
z[0] ),
828 b1( r1*cosPhi, r1*sinPhi,
z[1] );
938 if (opposite) rx = -rx;
943 G4double deltaR = rx -
r[0], deltaZ = zx -
z[0];
968 if (edgeRZnorm) *edgeRZnorm = answer;
990 if (d2 < d1) d1 =
d2;
997 distOutside2 += dist*dist;
1000 *edgeRZnorm =
std::max(std::fabs(*edgeRZnorm),std::fabs(dist));
1022 G4double rx = p.perp(), zx = p.z();
1028 if (rx < 0) part = -1;
1033 G4double deltaR = rx -
r[0]*part, deltaZ = zx -
z[0];
1053 deltaR = rx -
r[1]*part;
1061 if (edgeRZnorm) *edgeRZnorm = answer;
1083 if (d2 < d1) d1 =
d2;
1090 distOutside2 += dist*dist;
1093 *edgeRZnorm =
std::max(std::fabs(*edgeRZnorm),std::fabs(dist));
1117 if (!
cone->
HitOn( rx, hit.z() ))
return false;
1128 while( phi <
startPhi-phiTolerant )
1140 qb = qx - corners[3];
1143 if (normSign*qacb.dot(v) < 0)
return false;
1145 else if (phi < phiTolerant)
1149 qb = qx - corners[0];
1152 if (normSign*qacb.dot(v) < 0)
return false;
1186 G4double deter = tx1*ty2 - tx2*ty1;
1188 G4double s1 = ((x2-x1)*ty2 - tx2*(y2-y1))/deter;
1189 G4double s2 = ((x2-x1)*ty1 - tx1*(y2-y1))/deter;
1195 x = 0.5*( x1+s1*tx1 + x2+s2*tx2 );
1196 y = 0.5*( y1+s1*ty1 + y2+s2*ty2 );
1221 dr=
r[1]-
r[0];dz=
z[1]-
z[0];
1242 zz = z[0]+(rr-r[0])*dz/dr;
G4ThreeVector GetPointOnFace()
static G4GEOM_DLL G4PlSideManager subInstanceManager
static const G4double kInfinity
CLHEP::Hep3Vector G4ThreeVector
static const G4float tolerance
G4bool Intersect(const G4ThreeVector &p, const G4ThreeVector &v, G4bool outgoing, G4double surfTolerance, G4double &distance, G4double &distFromSurface, G4ThreeVector &normal, G4bool &isAllBehind)
G4double GetSurfaceTolerance() const
void SetNormal(const G4ThreeVector &newNormal)
virtual G4bool PartialClip(const G4VoxelLimits &voxelLimit, const EAxis IgnoreMe)
virtual void AddVertexInOrder(const G4ThreeVector vertex)
EInside Inside(const G4ThreeVector &p, G4double tolerance, G4double *bestDistance)
virtual void ClearAllVertices()
G4bool HitOn(const G4double r, const G4double z)
G4double Extent(const G4ThreeVector axis)
G4ThreeVector Normal(const G4ThreeVector &p, G4double *bestDistance)
G4int LineHitsCone(const G4ThreeVector &p, const G4ThreeVector &v, G4double *s1, G4double *s2)
const G4int kMinMeshSections
G4int CreateSubInstance()
void CopyStuff(const G4PolyconeSide &source)
static void FindLineIntersect(G4double x1, G4double y1, G4double tx1, G4double ty1, G4double x2, G4double y2, G4double tx2, G4double ty2, G4double &x, G4double &y)
virtual ~G4PolyconeSide()
static double normal(HepRandomEngine *eptr)
void CalculateExtent(const EAxis axis, const G4VoxelLimits &voxelLimit, const G4AffineTransform &tranform, G4SolidExtentList &extentList)
G4PolyconeSide(const G4PolyconeSideRZ *prevRZ, const G4PolyconeSideRZ *tail, const G4PolyconeSideRZ *head, const G4PolyconeSideRZ *nextRZ, G4double phiStart, G4double deltaPhi, G4bool phiIsOpen, G4bool isAllBehind=false)
static const double twopi
void AddSurface(const G4ClippablePolygon &surface)
static const G4PlSideManager & GetSubInstanceManager()
G4IntersectingCone * cone
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double GetPhi(const G4ThreeVector &p)
const G4double x[NPOINTSGL]
G4double DistanceAway(const G4ThreeVector &p, G4bool opposite, G4double &distOutside2, G4double *rzNorm=0)
G4bool PointOnCone(const G4ThreeVector &hit, G4double normSign, const G4ThreeVector &p, const G4ThreeVector &v, G4ThreeVector &normal)
G4double Distance(const G4ThreeVector &p, G4bool outgoing)
G4PolyconeSide & operator=(const G4PolyconeSide &source)
const G4double kMeshAngleDefault
static G4GeometryTolerance * GetInstance()
const G4int kMaxMeshSections