| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer << 3 // * DISCLAIMER *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers *
6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which *
7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: *
8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license *
9 // * include a list of copyright holders. <<
10 // * 9 // * *
11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file << 14 // * use. *
16 // * for the full disclaimer and the limitatio <<
17 // * 15 // * *
18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the *
19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. *
20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work *
21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this *
22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. *
23 // * acceptance of all terms of the Geant4 Sof <<
24 // ******************************************* 21 // ********************************************************************
25 // 22 //
>> 23 //
>> 24 // $Id: G4VTwistSurface.hh,v 1.4 2005/12/09 13:38:07 link Exp $
>> 25 // GEANT4 tag $Name: geant4-08-00-patch-01 $
>> 26 //
>> 27 //
>> 28 // --------------------------------------------------------------------
>> 29 // GEANT 4 class header file
>> 30 //
>> 31 //
26 // G4VTwistSurface 32 // G4VTwistSurface
27 // 33 //
28 // Class description: 34 // Class description:
29 // 35 //
30 // Abstract base class for boundary surface of << 36 // Abstract base class for boundary surface of G4VSolid.
31 37
32 // 01-Aug-2002 - Kotoyo Hoshina (hoshina@hepbu << 38 // Author:
33 // 13-Nov-2003 - O.Link (Oliver.Link@cern.ch), << 39 // 01-Aug-2002 - Kotoyo Hoshina (hoshina@hepburn.s.chiba-u.ac.jp)
34 // from original version in Jupi << 40 //
>> 41 // History:
>> 42 // 13-Nov-2003 - O.Link (Oliver.Link@cern.ch), Integration in Geant4
>> 43 // from original version in Jupiter-2.5.02 application.
35 // ------------------------------------------- 44 // --------------------------------------------------------------------
36 #ifndef G4VTWISTSURFACE_HH << 45 #ifndef __G4VTWISTSURFACE__
37 #define G4VTWISTSURFACE_HH << 46 #define __G4VTWISTSURFACE__
38 <<
39 #include <CLHEP/Units/SystemOfUnits.h> <<
40 47
41 #include "G4VSolid.hh" 48 #include "G4VSolid.hh"
42 #include "geomdefs.hh" 49 #include "geomdefs.hh"
43 50
44 #include "G4RotationMatrix.hh" 51 #include "G4RotationMatrix.hh"
45 52
46 #define G4VSURFACENXX 10 53 #define G4VSURFACENXX 10
47 54
48 class G4VTwistSurface 55 class G4VTwistSurface
49 { 56 {
50 public: << 57 public: // without description
>> 58
>> 59 enum EValidate {kDontValidate = 0, kValidateWithTol = 1,
>> 60 kValidateWithoutTol = 2, kUninitialized = 3};
51 61
52 enum EValidate { kDontValidate = 0, kValida << 62 public: // with description
53 kValidateWithoutTol = 2, k <<
54 63
55 G4VTwistSurface (const G4String& name); << 64 G4VTwistSurface (const G4String &name);
56 G4VTwistSurface (const G4String& name, << 65 G4VTwistSurface (const G4String &name,
57 const G4RotationMatrix& ro << 66 const G4RotationMatrix &rot,
58 const G4ThreeVector& tl << 67 const G4ThreeVector &tlate,
59 G4int ha << 68 G4int handedness,
60 const EAxis ax << 69 const EAxis axis1,
61 const EAxis ax << 70 const EAxis axis2,
62 G4double ax << 71 G4double axis0min = -kInfinity,
63 G4double ax << 72 G4double axis1min = -kInfinity,
64 G4double ax << 73 G4double axis0max = kInfinity,
65 G4double ax << 74 G4double axis1max = kInfinity);
66 75
67 virtual ~G4VTwistSurface() = default; << 76 virtual ~G4VTwistSurface();
68 77
69 virtual G4int AmIOnLeftSide(const G4Thr << 78 virtual G4int AmIOnLeftSide(const G4ThreeVector &me,
70 const G4Thr << 79 const G4ThreeVector &vec,
71 G4boo 80 G4bool withTol = true);
72 81
73 virtual G4double DistanceToBoundary( 82 virtual G4double DistanceToBoundary( G4int areacode,
74 << 83 G4ThreeVector &xx,
75 const << 84 const G4ThreeVector &p) ;
76 85
77 virtual G4double DistanceToIn(const G4Thre << 86
78 const G4Thre << 87 virtual G4double DistanceToIn(const G4ThreeVector &gp,
79 G4Thre << 88 const G4ThreeVector &gv,
80 virtual G4double DistanceToOut(const G4Thr << 89 G4ThreeVector &gxxbest);
81 const G4Thr << 90 virtual G4double DistanceToOut(const G4ThreeVector &gp,
82 G4Thr << 91 const G4ThreeVector &gv,
83 virtual G4double DistanceTo(const G4ThreeV << 92 G4ThreeVector &gxxbest);
84 G4ThreeV << 93 virtual G4double DistanceTo(const G4ThreeVector &gp,
>> 94 G4ThreeVector &gxx);
85 95
86 virtual G4int DistanceToSurface(const G << 96 virtual G4int DistanceToSurface(const G4ThreeVector &gp,
87 const G << 97 const G4ThreeVector &gv,
88 G 98 G4ThreeVector gxx[],
89 G 99 G4double distance[],
90 G 100 G4int areacode[],
91 G 101 G4bool isvalid[],
92 EValida 102 EValidate validate=kValidateWithTol) = 0;
93 103
94 virtual G4int DistanceToSurface(const G << 104 virtual G4int DistanceToSurface(const G4ThreeVector &gp,
95 G 105 G4ThreeVector gxx[],
96 G 106 G4double distance[],
97 G 107 G4int areacode[]) = 0;
98 108
99 void DebugPrint() const; 109 void DebugPrint() const;
100 110
101 virtual G4ThreeVector GetNormal(const G4Thr << 111 // get methods
>> 112
>> 113 virtual G4ThreeVector GetNormal(const G4ThreeVector &xx,G4bool isGlobal) = 0;
102 114
103 virtual G4String GetName() const { ret 115 virtual G4String GetName() const { return fName; }
104 virtual void GetBoundaryParameters << 116 virtual void GetBoundaryParameters(const G4int &areacode,
105 << 117 G4ThreeVector &d,
106 << 118 G4ThreeVector &x0,
107 << 119 G4int &boundarytype) const;
108 virtual G4ThreeVector GetBoundaryAtPZ(G4int 120 virtual G4ThreeVector GetBoundaryAtPZ(G4int areacode,
109 const << 121 const G4ThreeVector &p) const;
110 122
111 inline G4double DistanceToPlaneWithV(const << 123 inline G4double DistanceToPlaneWithV(const G4ThreeVector &p,
112 const << 124 const G4ThreeVector &v,
113 const << 125 const G4ThreeVector &x0,
114 const << 126 const G4ThreeVector &n0,
115 << 127 G4ThreeVector &xx);
116 << 128
117 inline G4double DistanceToPlane(const G4Th << 129 inline G4double DistanceToPlane(const G4ThreeVector &p,
118 const G4Th << 130 const G4ThreeVector &x0,
119 const G4Th << 131 const G4ThreeVector &n0,
120 G4Th << 132 G4ThreeVector &xx);
121 133
122 inline G4double DistanceToPlane(const G4Th << 134 inline G4double DistanceToPlane(const G4ThreeVector &p,
123 const G4Th << 135 const G4ThreeVector &x0,
124 const G4Th << 136 const G4ThreeVector &t1,
125 const G4Th << 137 const G4ThreeVector &t2,
126 G4Th << 138 G4ThreeVector &xx,
127 G4Th << 139 G4ThreeVector &n);
128 << 140
129 inline G4double DistanceToLine (const G4Th << 141 inline G4double DistanceToLine (const G4ThreeVector &p,
130 const G4Th << 142 const G4ThreeVector &x0,
131 const G4Th << 143 const G4ThreeVector &d,
132 G4Th << 144 G4ThreeVector &xx);
133 << 145
134 inline G4bool IsAxis0 (G4int areacode) c << 146 inline G4bool IsAxis0 (G4int areacode) const;
135 inline G4bool IsAxis1 (G4int areacode) c << 147 inline G4bool IsAxis1 (G4int areacode) const;
136 inline G4bool IsOutside (G4int areacode) c << 148 inline G4bool IsOutside (G4int areacode) const;
137 inline G4bool IsInside (G4int areacode, G << 149 inline G4bool IsInside (G4int areacode, G4bool testbitmode = false) const;
138 inline G4bool IsBoundary (G4int areacode, G << 150 inline G4bool IsBoundary (G4int areacode, G4bool testbitmode = false) const;
139 inline G4bool IsCorner (G4int areacode, G << 151 inline G4bool IsCorner (G4int areacode, G4bool testbitmode = false) const;
140 inline G4bool IsValidNorm() const { return << 152 inline G4bool IsValidNorm() const { return fIsValidNorm; }
141 G4bool IsSameBoundary (G4VTwistSurface* sur << 153 G4bool IsSameBoundary (G4VTwistSurface *surface1, G4int areacode1,
142 G4VTwistSurface* sur << 154 G4VTwistSurface *surface2, G4int areacode2 ) const;
143 inline G4int GetAxisType(G4int areacode, G << 155 inline G4int GetAxisType(G4int areacode, G4int whichaxis) const;
144 << 156
145 inline G4ThreeVector ComputeGlobalPoint << 157 inline G4ThreeVector ComputeGlobalPoint (const G4ThreeVector &lp) const;
146 inline G4ThreeVector ComputeLocalPoint << 158 inline G4ThreeVector ComputeLocalPoint (const G4ThreeVector &gp) const;
147 inline G4ThreeVector ComputeGlobalDirection << 159 inline G4ThreeVector ComputeGlobalDirection (const G4ThreeVector &lp) const;
148 inline G4ThreeVector ComputeLocalDirection << 160 inline G4ThreeVector ComputeLocalDirection (const G4ThreeVector &gp) const;
149 161
150 // set methods 162 // set methods
151 163
152 inline void SetAxis(G4int i, const EAxis ax 164 inline void SetAxis(G4int i, const EAxis axis) { fAxis[i] = axis; }
153 inline void SetNeighbours(G4VTwistSurface* << 165 inline void SetNeighbours(G4VTwistSurface* axis0min, G4VTwistSurface* axis1min,
154 G4VTwistSurface* << 166 G4VTwistSurface* axis0max, G4VTwistSurface* axis1max);
155 167
156 virtual G4ThreeVector SurfacePoint(G4double << 168 virtual G4ThreeVector SurfacePoint(G4double , G4double,
157 G4bool i 169 G4bool isGlobal = false ) = 0 ;
158 virtual G4double GetBoundaryMin(G4double) = 170 virtual G4double GetBoundaryMin(G4double) = 0 ;
159 virtual G4double GetBoundaryMax(G4double) = 171 virtual G4double GetBoundaryMax(G4double) = 0 ;
160 virtual G4double GetSurfaceArea() = 0 ; 172 virtual G4double GetSurfaceArea() = 0 ;
161 virtual void GetFacets(G4int m, G4int n, G4 173 virtual void GetFacets(G4int m, G4int n, G4double xyz[][3],
162 G4int faces[][4], G4 174 G4int faces[][4], G4int iside) = 0 ;
163 G4int GetNode( G4int i, G4int j, G4int m, G << 175 G4int GetNode( G4int i, G4int j, G4int m, G4int n, G4int iside ) ;
164 G4int GetFace( G4int i, G4int j, G4int m, G << 176 G4int GetFace( G4int i, G4int j, G4int m, G4int n, G4int iside ) ;
165 G4int GetEdgeVisibility( G4int i, G4int j, << 177 G4int GetEdgeVisibility( G4int i, G4int j, G4int m, G4int n, G4int number, G4int orientation) ;
166 G4int number, G4in << 178
>> 179
>> 180 public: // without description
167 181
168 G4VTwistSurface(__void__&); 182 G4VTwistSurface(__void__&);
169 // Fake default constructor for usage res 183 // Fake default constructor for usage restricted to direct object
170 // persistency for clients requiring prea 184 // persistency for clients requiring preallocation of memory for
171 // persistifiable objects. 185 // persistifiable objects.
172 186
173 protected: << 187 protected: // with description
174 188
175 inline G4VTwistSurface** GetNeighbours() { << 189 // get methods
176 inline G4int GetNeighbours(G4int areacode, << 190
177 inline G4ThreeVector GetCorner(G4int areaco << 191 inline G4VTwistSurface** GetNeighbours() { return fNeighbours; }
178 void GetBoundaryAxis(G4int areacode, EAxis << 192 inline G4int GetNeighbours(G4int areacode, G4VTwistSurface* surfaces[]);
179 void GetBoundaryLimit(G4int areacode, G4dou << 193 inline G4ThreeVector GetCorner(G4int areacode) const;
180 virtual G4int GetAreaCode(const G4ThreeVect << 194 void GetBoundaryAxis(G4int areacode, EAxis axis[]) const;
181 << 195 void GetBoundaryLimit(G4int areacode, G4double limit[]) const;
182 virtual void SetBoundary(const G4int& << 196 virtual G4int GetAreaCode(const G4ThreeVector &xx, G4bool withtol=true) = 0;
183 const G4ThreeVecto << 197
184 const G4ThreeVecto << 198 // set methods
185 const G4int& << 199
>> 200 virtual void SetBoundary(const G4int &axiscode,
>> 201 const G4ThreeVector &direction,
>> 202 const G4ThreeVector &x0,
>> 203 const G4int &boundarytype);
186 // areacode must be one of them: 204 // areacode must be one of them:
187 // sAxis0 & sAxisMin, sAxis0 & sAxisMax, 205 // sAxis0 & sAxisMin, sAxis0 & sAxisMax,
188 // sAxis1 & sAxisMin, sAxis1 & sAxisMax. 206 // sAxis1 & sAxisMin, sAxis1 & sAxisMax.
189 // boundarytype represents the shape of l 207 // boundarytype represents the shape of locus
190 // from the start point to end point of b 208 // from the start point to end point of boundary.
191 // ex. 209 // ex.
192 // sAxisRho = linear line which start poi 210 // sAxisRho = linear line which start point is fixed at origin.
193 // sAxisPhi = part of circle which center 211 // sAxisPhi = part of circle which center placed at the origin.
194 212
195 void SetCorner(G4int areacode, G4double x, 213 void SetCorner(G4int areacode, G4double x, G4double y, G4double z);
196 214
197 private: 215 private:
198 216
199 virtual void SetBoundaries() = 0; 217 virtual void SetBoundaries() = 0;
200 virtual void SetCorners() = 0; 218 virtual void SetCorners() = 0;
201 219
202 // data members ----------------------------- 220 // data members ---------------------------------------------------------
203 221
204 public: 222 public:
205 223
206 static const G4int sOutside ; 224 static const G4int sOutside ;
207 static const G4int sInside ; 225 static const G4int sInside ;
208 static const G4int sBoundary; 226 static const G4int sBoundary;
209 static const G4int sCorner; 227 static const G4int sCorner;
210 static const G4int sC0Min1Min; 228 static const G4int sC0Min1Min;
211 static const G4int sC0Max1Min; 229 static const G4int sC0Max1Min;
212 static const G4int sC0Max1Max; 230 static const G4int sC0Max1Max;
213 static const G4int sC0Min1Max; 231 static const G4int sC0Min1Max;
214 static const G4int sAxisMin; 232 static const G4int sAxisMin;
215 static const G4int sAxisMax; 233 static const G4int sAxisMax;
216 static const G4int sAxisX; 234 static const G4int sAxisX;
217 static const G4int sAxisY; 235 static const G4int sAxisY;
218 static const G4int sAxisZ; 236 static const G4int sAxisZ;
219 static const G4int sAxisRho; 237 static const G4int sAxisRho;
220 static const G4int sAxisPhi; 238 static const G4int sAxisPhi;
221 static const G4int sAxis0; 239 static const G4int sAxis0;
222 static const G4int sAxis1; 240 static const G4int sAxis1;
223 static const G4int sSizeMask; 241 static const G4int sSizeMask;
224 static const G4int sAxisMask; 242 static const G4int sAxisMask;
225 static const G4int sAreaMask; 243 static const G4int sAreaMask;
226 244
227 protected: 245 protected:
228 246
229 class CurrentStatus 247 class CurrentStatus
230 { 248 {
231 public: << 249 public:
232 250
233 CurrentStatus(); 251 CurrentStatus();
234 virtual ~CurrentStatus(); 252 virtual ~CurrentStatus();
235 253
236 inline G4ThreeVector GetXX(G4int i) 254 inline G4ThreeVector GetXX(G4int i) const { return fXX[i]; }
237 inline G4double GetDistance(G4int i 255 inline G4double GetDistance(G4int i) const { return fDistance[i]; }
238 inline G4int GetAreacode(G4int i 256 inline G4int GetAreacode(G4int i) const { return fAreacode[i]; }
239 inline G4int GetNXX() 257 inline G4int GetNXX() const { return fNXX; }
240 inline G4bool IsDone() 258 inline G4bool IsDone() const { return fDone; }
241 inline G4bool IsValid(G4int i) 259 inline G4bool IsValid(G4int i) const { return fIsValid[i]; }
242 260
243 void SetCurrentStatus(G4int 261 void SetCurrentStatus(G4int i,
244 G4ThreeVector& << 262 G4ThreeVector &xx,
245 G4double& << 263 G4double &dist,
246 G4int& << 264 G4int &areacode,
247 G4bool& << 265 G4bool &isvalid,
248 G4int 266 G4int nxx,
249 EValidate 267 EValidate validate,
250 const G4ThreeVector* p, << 268 const G4ThreeVector *p,
251 const G4ThreeVector* v = << 269 const G4ThreeVector *v = 0);
252 270
253 void ResetfDone(EValidate validate, << 271 void ResetfDone(EValidate validate,
254 const G4ThreeVector* p, << 272 const G4ThreeVector *p,
255 const G4ThreeVector* v = nullp << 273 const G4ThreeVector *v = 0);
256 274
257 275
258 void DebugPrint() const; 276 void DebugPrint() const;
259 277
260 private: << 278 private:
261 279
262 G4double fDistance[G4VSURFAC 280 G4double fDistance[G4VSURFACENXX];
263 G4ThreeVector fXX[G4VSURFACENXX]; 281 G4ThreeVector fXX[G4VSURFACENXX];
264 G4int fAreacode[G4VSURFAC 282 G4int fAreacode[G4VSURFACENXX];
265 G4bool fIsValid[G4VSURFACE 283 G4bool fIsValid[G4VSURFACENXX];
266 G4int fNXX; 284 G4int fNXX;
267 G4ThreeVector fLastp; 285 G4ThreeVector fLastp;
268 G4ThreeVector fLastv; 286 G4ThreeVector fLastv;
269 EValidate fLastValidate; 287 EValidate fLastValidate;
270 G4bool fDone; 288 G4bool fDone;
271 }; 289 };
272 290
273 class Boundary 291 class Boundary
274 { 292 {
275 public: << 293 public:
276 << 294 Boundary();
277 Boundary() = default; << 295 virtual ~Boundary();
278 virtual ~Boundary() = default; <<
279 296
280 void SetFields(const G4int& areac << 297 void SetFields(const G4int &areacode,
281 const G4ThreeVector& d, << 298 const G4ThreeVector &d,
282 const G4ThreeVector& x0, << 299 const G4ThreeVector &x0,
283 const G4int& boun << 300 const G4int &boundarytype);
284 301
285 G4bool IsEmpty() const; 302 G4bool IsEmpty() const;
286 303
287 G4bool GetBoundaryParameters(const G4int << 304 G4bool GetBoundaryParameters(const G4int &areacode,
288 G4Thr << 305 G4ThreeVector &d,
289 G4Thr << 306 G4ThreeVector &x0,
290 G4int << 307 G4int &boundarytype) const;
291 308
292 private: << 309 private:
293 << 310 G4int fBoundaryAcode;
294 G4int fBoundaryAcode{-1}; <<
295 G4ThreeVector fBoundaryDirection; 311 G4ThreeVector fBoundaryDirection;
296 G4ThreeVector fBoundaryX0; 312 G4ThreeVector fBoundaryX0;
297 G4int fBoundaryType{0}; << 313 G4int fBoundaryType;
298 }; 314 };
299 315
300 EAxis fAxis[2]; 316 EAxis fAxis[2];
301 G4double fAxisMin[2]; 317 G4double fAxisMin[2];
302 G4double fAxisMax[2]; 318 G4double fAxisMax[2];
303 CurrentStatus fCurStatWithV; 319 CurrentStatus fCurStatWithV;
304 CurrentStatus fCurStat; 320 CurrentStatus fCurStat;
305 G4RotationMatrix fRot; 321 G4RotationMatrix fRot;
306 G4ThreeVector fTrans; 322 G4ThreeVector fTrans;
307 G4int fHandedness; 323 G4int fHandedness;
308 class G4SurfCurNormal 324 class G4SurfCurNormal
309 { 325 {
310 public: 326 public:
311 <<
312 G4ThreeVector p; 327 G4ThreeVector p;
313 G4ThreeVector normal; 328 G4ThreeVector normal;
314 }; 329 };
315 G4SurfCurNormal fCurrentNormal; 330 G4SurfCurNormal fCurrentNormal;
316 G4bool fIsValidNorm; 331 G4bool fIsValidNorm;
317 G4double kCarTolerance; <<
318 332
319 private: 333 private:
320 334
321 G4VTwistSurface* fNeighbours[4]; // {0,1,2, << 335 G4VTwistSurface *fNeighbours[4]; // {0,1,2,3} = sAxis0min, sAxis1min,
322 // << 336 // sAxis0max, sAxis1max
323 G4ThreeVector fCorners[4]; // corners o 337 G4ThreeVector fCorners[4]; // corners of the surface in local coordinate
324 Boundary fBoundaries[4]; // boundarie 338 Boundary fBoundaries[4]; // boundaries of the surface.
325 G4String fName; 339 G4String fName;
326 340
327 class G4SurfSideQuery 341 class G4SurfSideQuery
328 { 342 {
329 public: 343 public:
330 <<
331 G4ThreeVector me; 344 G4ThreeVector me;
332 G4ThreeVector vec; 345 G4ThreeVector vec;
333 G4bool withTol; 346 G4bool withTol;
334 G4int amIOnLeftSide; 347 G4int amIOnLeftSide;
335 }; 348 };
336 G4SurfSideQuery fAmIOnLeftSide; 349 G4SurfSideQuery fAmIOnLeftSide;
337 }; 350 };
338 351
339 //============================================ 352 //========================================================
340 // inline functions 353 // inline functions
341 //============================================ 354 //========================================================
342 355
343 struct Intersection 356 struct Intersection
344 { 357 {
345 G4double phi ; // parameter phi 358 G4double phi ; // parameter phi
346 G4double u ; // parameter u 359 G4double u ; // parameter u
347 G4ThreeVector xx ; // intersection point i 360 G4ThreeVector xx ; // intersection point in cartesian
348 G4double distance ; // distance to intersec 361 G4double distance ; // distance to intersection
349 G4int areacode ; // the areacode of the << 362 G4int areacode; // the areacode of the intersection
350 G4bool isvalid ; // valid intersection ? 363 G4bool isvalid ; // valid intersection ??
351 364
352 }; 365 };
353 366
354 inline 367 inline
355 G4bool DistanceSort( const Intersection& a, co << 368 G4bool DistanceSort( const Intersection &a, const Intersection &b)
356 { 369 {
357 return a.distance < b.distance ; 370 return a.distance < b.distance ;
358 } 371 }
359 372
360 inline 373 inline
361 G4bool EqualIntersection( const Intersection& << 374 G4bool EqualIntersection( const Intersection &a, const Intersection &b)
362 { 375 {
363 return ( ( a.xx - b.xx ).mag() < 1E-9*CLHEP: << 376 return ( ( a.xx - b.xx ).mag() < kCarTolerance ) ;
364 } 377 }
365 378
366 #include "G4VTwistSurface.icc" 379 #include "G4VTwistSurface.icc"
367 380
368 #endif 381 #endif
369 382