| Summary: | Wrong distance to In of G4Trap | ||
|---|---|---|---|
| Product: | Geant4 | Reporter: | John Apostolakis <John.Apostolakis> |
| Component: | geometry/solids | Assignee: | John Apostolakis <John.Apostolakis> |
| Status: | RESOLVED INVALID | ||
| Severity: | normal | ||
| Priority: | P5 | ||
| Version: | 9.5 | ||
| Hardware: | All | ||
| OS: | All | ||
|
Description
John Apostolakis
2012-11-23 12:24:19 CET
When running in check mode in a complex geometry, the following problem is reported
Solid gave DistanceToIn = 0 yet returns -kInside- for this point !
-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomNav1001
issued by : G4VoxelNavigation::ComputeStep()
Conflicting response from Solid at the projected *intersection* point with the solid:
Inaccurate solid DistanceToIn for solid E5AL
Original position (-7.671422800142409,-58.39097619339464,-56.72023453210477)
and direction (0.9482351604406044,0.2933922475180894,-0.1215362892327454)
Original distance to surface (Step): 0
Solid gave DistanceToIn = 0 yet returns -kInside- for this point !
Details of location and Solid:
Point = (-7.671422800142409,-58.39097619339464,-56.72023453210477)
Safety and distance values (from this intersection point):
Safety (from inside): DistanceToOut(p) = 15.48346644107644
Distance (from inside): DistanceToOut(p,v)= 32.43837579867777
Solid type: G4Trap.
Parameters of this solid are:
-----------------------------------------------------------
*** Dump for solid - E5AL ***
===================================================
Solid type: G4Trap
Parameters:
half length Z: 126.495 mm
half length Y of face -fDz: 60.049 mm
half length X of side -fDy1, face -fDz: 21.939 mm
half length X of side +fDy1, face -fDz: 21.939 mm
half length Y of face +fDz: 71.598 mm
half length X of side -fDy2, face +fDz: 26.355 mm
half length X of side +fDy2, face +fDz: 26.355 mm
std::tan(theta)*std::cos(phi): 0.0003145748193832724 degrees
std::tan(theta)*std::sin(phi): 20.18343727892806 degrees
std::tan(alpha), -fDz: 0 degrees
std::tan(alpha), +fDz: 0 degrees
trap side plane equations:
0 X + -0.956067407783281 Y + 0.2931469115896627 Z + -62.93170301622281 = 0
0 X + 0.9291424254719876 Y + -0.3697219944607462 Z + -61.15940644305538 = 0
-0.9998477879592079 X + 0 Y + -0.01744708895715516 Z + -24.14332453585099 = 0
0.9998475963754998 X + 0 Y + -0.01745806467324904 Z + -24.1433199096792 = 0
-----------------------------------------------------------
*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------
Dear John, By looking into the solid definition I found that given point is really Inside. Values of DistanceToOut(p)=15.48 and DistanceToOut(p,v)=32.44 confirm this. For the point, situated Inside, it is normal, that DistanceToIn replay 0. This problem seems to be more navigation problem, that solid problem. Best Regards, Tatiana. Dear John, The given point is really Inside of G4Trap, in this case DistanceToIn returns 0. This problem seems to be a problem of G4Naviagator. May be tested geometry has overlaps? Best Regards, Tatiana. Since it appears the problem is not in G4Trap I am closing it. |