Changeset 1347 for trunk/source/geometry/solids/CSG/test/testG4Sphere.cc
- Timestamp:
- Dec 22, 2010, 3:52:27 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/geometry/solids/CSG/test/testG4Sphere.cc
r1316 r1347 25 25 // 26 26 27 // $Id: testG4Sphere.cc,v 1. 29 2009/05/14 13:22:44 tnikitinExp $28 // GEANT4 tag $Name: geant4-09-04- beta-cand-01$27 // $Id: testG4Sphere.cc,v 1.30 2010/03/24 09:50:03 grichine Exp $ 28 // GEANT4 tag $Name: geant4-09-04-ref-00 $ 29 29 // 30 30 // G4Sphere Test File … … 96 96 { 97 97 G4double kAngTolerance = G4GeometryTolerance::GetInstance()->GetAngularTolerance(); 98 98 EInside inside; 99 99 G4int i, iMax; 100 100 G4double Dist, vol, volCheck; … … 292 292 G4ThreeVector vb830(0.1644697995,0.18507236,0.9688642354); 293 293 294 295 G4Sphere s18_03_10("s18_03_10", 47*mm, 48*mm, 0*deg, 200*deg, 80*deg, 100*deg); 296 G4ThreeVector p18_03_10(43.37539710867407*mm, 16.12036885157033*mm, -8.224548565698871*mm); 297 G4ThreeVector v18_03_10(-0.4161958548132512, 0.6603942936714858, -0.6250283092167956); 298 299 inside = s18_03_10.Inside(p18_03_10) ; 300 G4cout<<"s18_03_10.Inside(p18_03_10 ... = "<<OutputInside(inside)<<G4endl ; 301 302 303 294 304 #ifdef NDEBUG 295 305 G4Exception("FAIL: *** Assertions must be compiled in! ***"); … … 347 357 // Check G4Sphere::Inside 348 358 349 EInsideinside = s7.Inside(G4ThreeVector(1399.984667238032,359 inside = s7.Inside(G4ThreeVector(1399.984667238032, 350 360 5.9396696802500299, 351 361 -2.7661927818688308) ) ; … … 501 511 Dist=sn22.DistanceToOut(G4ThreeVector(-45.,0.,0.),vx,calcNorm,pgoodNorm,pNorm); 502 512 assert(ApproxEqual(Dist,45.)); 503 504 513 G4cout<<"Dist from Center ="<<sn22.DistanceToOut(G4ThreeVector(0.,49.,0),G4ThreeVector(0,-1,0))<<G4endl; 514 G4cout<<"Dist from Center ="<<sn22.DistanceToOut(G4ThreeVector(-45.,0.,0),G4ThreeVector(1,0,0))<<G4endl; 505 515 506 516 // … … 1049 1059 G4double d, G4ThreeVector dir, EInside exp); 1050 1060 1051 1061 G4Sphere SpAroundX("SpAroundX", 10.*mm, 1000.*mm, -1.0*degree, 2.0*degree, 0.*degree, 180.0*degree ); 1052 1062 1053 1063 G4double sinOneDeg = std::sin( 1.0 * degree ); … … 1091 1101 1092 1102 // Try the edges 1093 1103 G4ThreeVector NegEdgePt( radOne * std::cos(-270.0*degree), radOne * std::sin(-270.0*degree), 0.0); 1094 1104 G4ThreeVector PosEdgePt( radOne * std::cos(10.0*degree), radOne * std::sin(10.0*degree), 0.0); 1095 1105
Note: See TracChangeset
for help on using the changeset viewer.