[Gate-users] G4sphere not available on nanometer scale

[John Apostolakis]

Dear Wonmo Sung,

Could you confirm the version of Geant4 you are using?

I have two suggestions. The first is to try using the version of G4Sphere from the newest patch of Geant4. You can copy the source file into your installation, and, if needed only the header file and inline icc file from the include directory - and of course then the rebuild Geant4 libraries.

The second concerns the tolerance. The default value is 1.e-09 mm, which is not small enough for your geometry. You can see in the Users Guide how to adapt it to the size of your setup.

If your problem remains, I suggest to report it in the Geant4 Problem reporting system.

Regards,
John Apostolakis

On 10 Nov 2015, at 08:15, 성원모 <wsung at snu.ac.kr<mailto:wsung at snu.ac.kr>> wrote:


Dear Gate users


Hello

I'm currently using Gate v7.1.


I tried to use setRmin for hollow sphere on nanometer scale.

e.g.,

/gate/sphere/geometry/setRmin 1 nm

/gate/sphere/geometry/setRmax 20 nm


However, I got G4Exception as below.

This error does not occur if I replace nm with mm OR if I use 0 nm for setRmin (=full sphere).

Could you provide any comments or help to solve this problem? I appreciate any help in advance.


-----------------------------------------------------------
    *** Dump for solid - sphere_solid ***
    ===================================================
 Solid type: G4Sphere
 Parameters:
    inner radius: 1e-06 mm
    outer radius: 2.000000000000001e-05 mm
    starting phi of segment  : 0 degrees
    delta phi of segment     : 360 degrees
    starting theta of segment: 0 degrees
    delta theta of segment   : 180 degrees
-----------------------------------------------------------

-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : GeomSolids1002
      issued by : G4Sphere::DistanceToOut(p,v,..)
Logic error: snxt = kInfinity  ???
Position:

p.x() = 5.106433541911882e-07 mm
p.y() = -7.11253374326104e-07 mm
p.z() = 9.224550730873363e-08 mm

Rp = 8.80424460864498e-07 mm

Direction:

v.x() = 0.9926097739516453
v.y() = -0.02063914065141685
v.z() = 0.1195820326338147

Proposed distance :

snxt = 9e+99 mm

*** This is just a warning message. ***
-------- WWWW -------- G4Exception-END --------- WWWW -------


Regards,

Wonmo Sung

_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org<mailto:Gate-users at lists.opengatecollaboration.org>
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/mailman/private/gate-users/attachments/20151110/a922823c/attachment.html>