Message: G4Sphere::DistanceToOut problem Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question G4Sphere::DistanceToOut problem 

Forum: Event and Track Management
Date: 15 Sep, 2016
From: Anastasia Makarova <Anastasia Makarova>

Dear Geant4 community,

I am modelling a 10 nm gold sphere in water phantom. To score the dose around it I used G4VPVParameterisation class and created concentric spheres. All in all: 990 concentric spheres with 1 nm step starting from the surface of the particle.
The particle and detector spheres are placed in origin.

I have the following warning messages:
-----------------------------------------------------------
    *** Dump for solid - CocentrSphere ***
    ===================================================
 Solid type: G4Sphere
 Parameters:
    inner radius: 1e-05 mm
    outer radius: 1.1e-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() = -3.957318597258563e-06 mm
p.y() = -1.048985396323587e-06 mm
p.z() = 9.123555174573557e-06 mm

Rp = 9.999999993269828e-06 mm

Direction:

v.x() = 0
v.y() = 0
v.z() = 1

Proposed distance :

snxt = 9e+99 mm

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

In my parallel world Construct() function I used the following code:

G4VPVParameterisation* sphereParam = new GoldSphereParameterisation;
new G4PVParameterised("Spheres", CocentrSphereLog, outerSphereLog, kUndefined, 990, sphereParam);

My parametrisation class ComputeDimentions() function is as follows:

void ComputeDimensions (G4Sphere& sphere, const G4int copyNo, const G4VPhysicalVolume* physVol) const
    {
        G4double step = 1.0*nm;
        sphere.SetInnerRadius(10*nm + step*copyNo);
        sphere.SetOuterRadius(10*nm + step*(copyNo+1));
    }

How can I tackle this problem? The warning messages consume a lot of computing time.

Thank you for your help,
Anastasia





   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/09/15/01.06-26250-GoldParallelWorld.cc
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2016/09/15/01.06-78235-SphereParameterisation.hh

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: G4Sphere::DistanceToOut problem   (Makoto Asai - 15 Sep, 2016)
 Add Message Add Message
to: "G4Sphere::DistanceToOut problem"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]