<div dir="ltr"><div>I am using GATE version 9.0 (March 2020), Geant4 version Name: geant4-10-06-patch-02    (29-May-2020)</div><div>My phase space files that produce the source particles are from IAEAphsp files and the headers look like this:</div><div><br></div><div>$RECORD_CONTENTS:<br>  1   // X is stored ?<br>  1   // Y is stored ?<br>  0   // Z is stored ?<br>  1   // U is stored ?<br>  1   // V is stored ?<br>  1   // W is stored ?<br>  0   // Weight is stored ?<br>  0   // Extra floats stored ?<br>  0   // Extra longs stored ?</div><div><br></div><div>Thanks</div><div><br></div><div>GT<br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Mar 3, 2021 at 12:43 PM G Tom <<a href="mailto:gregthom992@gmail.com">gregthom992@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi All, <br></div><div><br></div><div>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.</div><div><br></div><div>
<img src="cid:ii_kltpwvrp0" alt="image.png" width="541" height="326">

</div><div><br></div><div>When I read the output after 10^8 histories, the dose actor recorded only 0s. What could I be missing ? <br></div><div><br></div><div>Thanks <br></div><div><br></div><div>GT<br></div><div><br></div><div>My macro is below and placed just before the /gate/run/initialize line in my main macro and I start the simulations with <br></div><div>#=====================================================<br># START<br>#=====================================================<br>/gate/random/setEngineName MersenneTwister<br>/gate/random/setEngineSeed auto<br><br>/gate/application/setTotalNumberOfPrimaries {NumParticles}<br>/gate/application/start</div><div><br></div><div><br></div><div>#=====================================================<br># ACTORS: stats, dose, fluence etc<br>#=====================================================<br>#<ACTORS><br>/control/alias BluePhantom2XLength 480.0<br>/control/alias BluePhantom2YLength 480.0<br>/control/alias BluePhantom2ZLength 310.0<br>/control/alias BluePhantom2Material Water<br><br><br>/control/divide HalfBluePhantom2XLength {BluePhantom2XLength} 2.0<br>/control/divide HalfBluePhantom2YLength {BluePhantom2YLength} 2.0<br>/control/divide HalfBluePhantom2ZLength {BluePhantom2ZLength} 2.0<br><br>/control/alias BluePhantom2XPlacement 0.00<br>/control/alias BluePhantom2YPlacement 0.00<br><br>/control/subtract BluePhantom2ZPlacement 1000 {SSD}<br>/control/subtract BluePhantom2ZPlacement {BluePhantom2ZPlacement} {HalfBluePhantom2ZLength}<br><br>/gate/world/daughters/name BluePhantom2<br>/gate/world/daughters/insert box<br>/gate/BluePhantom2/setMaterial {BluePhantom2Material}<br>/gate/BluePhantom2/geometry/setXLength {BluePhantom2XLength} mm<br>/gate/BluePhantom2/geometry/setYLength {BluePhantom2YLength} mm<br>/gate/BluePhantom2/geometry/setZLength {BluePhantom2ZLength} mm<br>/gate/BluePhantom2/placement/setTranslation {BluePhantom2XPlacement} {BluePhantom2YPlacement} {BluePhantom2ZPlacement} mm<br>/gate/BluePhantom2/vis/setColor cyan<br>/gate/BluePhantom2/vis/forceWireframe 1<br>/gate/BluePhantom2/vis/setVisible 1<br><br>/gate/physics/Positron/SetCutInRegion BluePhantom2 0.01 mm<br>/gate/physics/Electron/SetCutInRegion BluePhantom2 0.01 mm<br>/gate/physics/Gamma/SetCutInRegion BluePhantom2 0.01 mm<br><br><br><br>/control/alias DoseActorXVoxelSize 2 #mm<br>/control/alias DoseActorYVoxelSize       2 #mm<br>/control/alias DoseActorZVoxelSize       1 #mm<br><br>/control/alias DoseActorXPlacement     0<br>/control/alias DoseActorYPlacement   0<br>/control/alias DoseActorZPlacement   0<br><br>/control/alias DoseActorXSize      {BluePhantom2XLength}<br>/control/alias DoseActorYSize    {BluePhantom2YLength}<br>/control/alias DoseActorZSize    {BluePhantom2ZLength}<br><br>/control/divide DoseActorXResolution   {DoseActorXSize}        {DoseActorXVoxelSize}<br>/control/divide DoseActorYResolution     {DoseActorYSize}        {DoseActorYVoxelSize}<br>/control/divide DoseActorZResolution     {DoseActorZSize}        {DoseActorZVoxelSize}<br># Dose<br>/gate/actor/addActor DoseActor doseMeasurement<br>/gate/actor/doseMeasurement/attachTo BluePhantom2<br>/gate/actor/doseMeasurement/save {OutputDirName}/3d.mhd<br>/gate/actor/doseMeasurement/stepHitType random<br>/gate/actor/doseMeasurement/setPosition {DoseActorXPlacement} {DoseActorYPlacement} {DoseActorZPlacement} mm<br>/gate/actor/doseMeasurement/setVoxelSize   {DoseActorXVoxelSize} {DoseActorYVoxelSize} {DoseActorZVoxelSize} mm<br>#/gate/actor/doseMeasurement/setResolution  {DoseActorXResolution} {DoseActorYResolution} {DoseActorZResolution}<br>/gate/actor/doseMeasurement/setSize  {DoseActorXSize} {DoseActorYSize} {DoseActorZSize} mm<br>/gate/actor/doseMeasurement/enableEdep  false<br>/gate/actor/doseMeasurement/enableUncertaintyEdep  false<br>/gate/actor/doseMeasurement/enableSquaredEdep  false<br>/gate/actor/doseMeasurement/enableDose  true<br>/gate/actor/doseMeasurement/enableSquaredDose  true<br>/gate/actor/doseMeasurement/enableUncertaintyDose  false<br>/gate/actor/doseMeasurement/enableNumberOfHits false</div><div><br></div><div><br></div><div><br></div><div><br></div></div>
</blockquote></div>