Message: Re: Track stuck or not moving with too small voxel sizes Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Track stuck or not moving with too small voxel sizes 

Forum: Geometry
Re: None Track stuck or not moving with too small voxel sizes (shirin Abbasinejad Enger)
Re: None Re: Track stuck or not moving with too small voxel sizes (shirin Abbasinejad Enger)
Re: None Re: Track stuck or not moving with too small voxel sizes (John Apostolakis)
Date: 27 Nov, 2012
From: Joseph Perl <Joseph Perl>

I've reproduced the problem myself. It does not have to do with the number or size of the voxels. I think it's a fairly basic problem that we just didn't realize because we are generally focused on external beams hitting our detectors rather than point sources within our detectors.

Here's what I do to reproduce the problem. Place an isotropic point source at the center of the world. Place a nested parameterization box at the center of the world.

If the box is undivided, no problem.

If the box is divided into an even number of voxels, I see many track stuck warnings.

If the box is divided into an odd number of voxels, I do not see track stuck warnings.

The behavior is consistent whether I divided the box by a small even number, such as 2 x 2 x 2 (gives problems) or 3 x 3 x 3 (no problems), or by a large even number, such as 200 x 200 x 200 (gives problems) or 201 x 201 x 201 (no problems).

If the box is divided into an even number of voxels, but I offset my box from the center of the world, even by a tiny bit, the warnings go away.

In my case, I am using 200 MeV protons with a modular physics list consisting only of g4em-standard_opt0. But I do not believe the particle type or physics list is the issue.

 Add Message Add Message
to: "Re: Track stuck or not moving with too small voxel sizes"

 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 ]