[Gate-users] DoseActor not recording any dose (empty 0's files)

G Tom gregthom992 at gmail.com
Wed Mar 3 18:46:43 CET 2021


I am using GATE version 9.0 (March 2020), Geant4 version Name:
geant4-10-06-patch-02    (29-May-2020)
My phase space files that produce the source particles are from IAEAphsp
files and the headers look like this:

$RECORD_CONTENTS:
  1   // X is stored ?
  1   // Y is stored ?
  0   // Z is stored ?
  1   // U is stored ?
  1   // V is stored ?
  1   // W is stored ?
  0   // Weight is stored ?
  0   // Extra floats stored ?
  0   // Extra longs stored ?

Thanks

GT



On Wed, Mar 3, 2021 at 12:43 PM G Tom <gregthom992 at gmail.com> wrote:

> Hi All,
>
> This has been puzzling me all week. I set up my simulation and I can
> visualize it working and interactions happening in the phantom to which I
> have attached a doseActor.
>
> [image: image.png]
>
> When I read the output after 10^8 histories, the dose actor recorded only
> 0s. What could I be missing ?
>
> Thanks
>
> GT
>
> My macro is below and placed just before the /gate/run/initialize line in
> my main macro and I start the simulations with
> #=====================================================
> # START
> #=====================================================
> /gate/random/setEngineName MersenneTwister
> /gate/random/setEngineSeed auto
>
> /gate/application/setTotalNumberOfPrimaries {NumParticles}
> /gate/application/start
>
>
> #=====================================================
> # ACTORS: stats, dose, fluence etc
> #=====================================================
> #<ACTORS>
> /control/alias BluePhantom2XLength 480.0
> /control/alias BluePhantom2YLength 480.0
> /control/alias BluePhantom2ZLength 310.0
> /control/alias BluePhantom2Material Water
>
>
> /control/divide HalfBluePhantom2XLength {BluePhantom2XLength} 2.0
> /control/divide HalfBluePhantom2YLength {BluePhantom2YLength} 2.0
> /control/divide HalfBluePhantom2ZLength {BluePhantom2ZLength} 2.0
>
> /control/alias BluePhantom2XPlacement 0.00
> /control/alias BluePhantom2YPlacement 0.00
>
> /control/subtract BluePhantom2ZPlacement 1000 {SSD}
> /control/subtract BluePhantom2ZPlacement {BluePhantom2ZPlacement}
> {HalfBluePhantom2ZLength}
>
> /gate/world/daughters/name BluePhantom2
> /gate/world/daughters/insert box
> /gate/BluePhantom2/setMaterial {BluePhantom2Material}
> /gate/BluePhantom2/geometry/setXLength {BluePhantom2XLength} mm
> /gate/BluePhantom2/geometry/setYLength {BluePhantom2YLength} mm
> /gate/BluePhantom2/geometry/setZLength {BluePhantom2ZLength} mm
> /gate/BluePhantom2/placement/setTranslation {BluePhantom2XPlacement}
> {BluePhantom2YPlacement} {BluePhantom2ZPlacement} mm
> /gate/BluePhantom2/vis/setColor cyan
> /gate/BluePhantom2/vis/forceWireframe 1
> /gate/BluePhantom2/vis/setVisible 1
>
> /gate/physics/Positron/SetCutInRegion BluePhantom2 0.01 mm
> /gate/physics/Electron/SetCutInRegion BluePhantom2 0.01 mm
> /gate/physics/Gamma/SetCutInRegion BluePhantom2 0.01 mm
>
>
>
> /control/alias DoseActorXVoxelSize 2 #mm
> /control/alias DoseActorYVoxelSize 2 #mm
> /control/alias DoseActorZVoxelSize 1 #mm
>
> /control/alias DoseActorXPlacement 0
> /control/alias DoseActorYPlacement 0
> /control/alias DoseActorZPlacement 0
>
> /control/alias DoseActorXSize {BluePhantom2XLength}
> /control/alias DoseActorYSize {BluePhantom2YLength}
> /control/alias DoseActorZSize {BluePhantom2ZLength}
>
> /control/divide DoseActorXResolution {DoseActorXSize} {DoseActorXVoxelSize}
> /control/divide DoseActorYResolution {DoseActorYSize} {DoseActorYVoxelSize}
> /control/divide DoseActorZResolution {DoseActorZSize} {DoseActorZVoxelSize}
> # Dose
> /gate/actor/addActor DoseActor doseMeasurement
> /gate/actor/doseMeasurement/attachTo BluePhantom2
> /gate/actor/doseMeasurement/save {OutputDirName}/3d.mhd
> /gate/actor/doseMeasurement/stepHitType random
> /gate/actor/doseMeasurement/setPosition {DoseActorXPlacement}
> {DoseActorYPlacement} {DoseActorZPlacement} mm
> /gate/actor/doseMeasurement/setVoxelSize   {DoseActorXVoxelSize}
> {DoseActorYVoxelSize} {DoseActorZVoxelSize} mm
> #/gate/actor/doseMeasurement/setResolution  {DoseActorXResolution}
> {DoseActorYResolution} {DoseActorZResolution}
> /gate/actor/doseMeasurement/setSize  {DoseActorXSize} {DoseActorYSize}
> {DoseActorZSize} mm
> /gate/actor/doseMeasurement/enableEdep  false
> /gate/actor/doseMeasurement/enableUncertaintyEdep  false
> /gate/actor/doseMeasurement/enableSquaredEdep  false
> /gate/actor/doseMeasurement/enableDose  true
> /gate/actor/doseMeasurement/enableSquaredDose  true
> /gate/actor/doseMeasurement/enableUncertaintyDose  false
> /gate/actor/doseMeasurement/enableNumberOfHits false
>
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210303/dc332f26/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 31916 bytes
Desc: not available
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210303/dc332f26/attachment-0001.png>


More information about the Gate-users mailing list