Problem 1396

Summary: Wrong distance to In of G4Trap
Product: Geant4 Reporter: John Apostolakis <John.Apostolakis>
Component: geometry/solidsAssignee: 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

    
Comment 1 John Apostolakis 2012-11-23 12:25:27 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 -------
Comment 2 tatiana.nikitina 2012-12-11 10:26:23 CET
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.
Comment 3 tatiana.nikitina 2013-02-18 10:19:38 CET
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.
Comment 4 John Apostolakis 2013-04-03 16:43:40 CEST
Since it appears the problem is not in G4Trap I am closing it.