<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Dears Gate users,</p>
<p><br>
</p>
<p>I tried to run simulations for one proton PBS at different source to surface distances (nozzle-patient-distances)  to investigate the effect of SSD variation on the dose distribution in water phantom. SSD were set to be  (100,75,50,25 and 5 cm) at every
 run I changed the source position in the z-direction. Dose was normalised to the maximum dose. But, I got the same distribution for the first four values which I believe is not correct. Also, I did a similar simulation for proton beam with a fix energy and
 different diameters (changing sigma_x and sigma_y) and I got the same dose distribution with a slight change. Until now I am unable to detect the error in my set-up.</p>
<p><br>
</p>
<p>I attached a copy of the main macro file and I would be grateful if you could help me to find out the error source. </p>
<p><br>
</p>
<p>Thanks in advance,</p>
<p>Awad</p>
<p><br>
</p>
<p></p>
<div><br>
</div>
<div>#-------------------oooooOOOOO00000OOOOOooooo---------------------#</div>
<div><br>
</div>
<div># Example of electron beam in  phantom (water)</div>
<div>#-------------------oooooOOOOO00000OOOOOooooo---------------------#</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># VERBOSE and VISUALISATION</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/control/execute mac/verbose.mac</div>
<div>#/control/execute mac/visu.mac</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># GEOMETRY</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/geometry/setMaterialDatabase data/GateMaterials.db</div>
<div>/gate/geometry/setMaterialDatabase data/patient-HUmaterials.db</div>
<div><br>
</div>
<div># WORLD</div>
<div>/gate/world/setMaterial                                          Air</div>
<div>/gate/world/geometry/setXLength                       5.0 m</div>
<div>/gate/world/geometry/setYLength                       5.0 m</div>
<div>/gate/world/geometry/setZLength                       5.0 m</div>
<div>#/gate/world/placement/setTranslation             0 0 0 mm</div>
<div><br>
</div>
<div># WATERBOX</div>
<div>/gate/world/daughters/name                               water_box</div>
<div>/gate/world/daughters/insert                               box</div>
<div>/gate/water_box/setMaterial                                Water</div>
<div>/gate/water_box/geometry/setXLength              30 cm</div>
<div>/gate/water_box/geometry/setYLength              30 cm</div>
<div>/gate/water_box/geometry/setZLength              50 cm</div>
<div>/gate/water_box/placement/setTranslation       0.0 0.0 0.0 cm</div>
<div>/gate/water_box/vis/setColor blue</div>
<div><br>
</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># PHYSICS</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/physics/addPhysicsList emstandard_opt3</div>
<div><br>
</div>
<div>/gate/physics/Gamma/SetCutInRegion      world 1 mm</div>
<div>/gate/physics/Electron/SetCutInRegion   world 1 mm</div>
<div>/gate/physics/Positron/SetCutInRegion   world 1 mm</div>
<div><br>
</div>
<div>/gate/physics/Gamma/SetCutInRegion      water_box 0.1 mm</div>
<div>/gate/physics/Electron/SetCutInRegion   water_box 0.1 mm</div>
<div>/gate/physics/Positron/SetCutInRegion   water_box 0.1 mm</div>
<div><br>
</div>
<div>/gate/physics/SetMaxStepSizeInRegion world 1 mm</div>
<div>/gate/physics/ActivateStepLimiter proton</div>
<div>/gate/physics/ActivateStepLimiter deuteron</div>
<div>/gate/physics/ActivateStepLimiter triton</div>
<div>/gate/physics/ActivateStepLimiter alpha</div>
<div>/gate/physics/ActivateStepLimiter GenericIon</div>
<div><br>
</div>
<div>/gate/physics/ActivateStepLimiter proton</div>
<div>/gate/physics/ActivateStepLimiter e-</div>
<div><br>
</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># Statistics actor</div>
<div>#=====================================================</div>
<div>/gate/actor/addActor                    SimulationStatisticActor stat</div>
<div>/gate/actor/stat/saveEveryNSeconds      60</div>
<div>/gate/actor/stat/save                   output/statistics.txt</div>
<div><br>
</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># Dose Actor</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/actor/addActor                                             DoseActor  dose</div>
<div>/gate/actor/dose/save                                           output/120.txt</div>
<div>/gate/actor/dose/attachTo    <span class="Apple-tab-span" style="white-space:pre">
</span>                            water_box</div>
<div>/gate/actor/dose/stepHitType                              random</div>
<div>/gate/actor/dose/setPosition                                0 0 0 cm</div>
<div>/gate/actor/dose/setResolution                           50 50 500</div>
<div>/gate/actor/dose/setSize                                        5  5 50 cm</div>
<div>/gate/actor/dose/saveEveryNSeconds                 60</div>
<div>/gate/actor/dose/enableEdep                               true </div>
<div>/gate/actor/dose/enableUncertaintyEdep          true</div>
<div>/gate/actor/dose/enableDose                               true</div>
<div>/gate/actor/dose/enableUncertaintyDose          true</div>
<div>/gate/actor/dose/enableNumberOfHits              true</div>
<div><br>
</div>
<div><br>
</div>
<div>#=====================================</div>
<div># SOURCE</div>
<div>#=====================================</div>
<div>/gate/source/addSource                                              PBS PencilBeam</div>
<div>/gate/source/PBS/setPosition                                      0  0 -100 cm</div>
<div>/gate/source/PBS/setEnergy                                        120 MeV</div>
<div>/gate/source/PBS/setSigmaEnergy                               0.1 MeV</div>
<div>/gate/source/PBS/setSigmaX                                        3 mm</div>
<div>/gate/source/PBS/setSigmaY                                        3 mm</div>
<div>/gate/source/PBS/setSigmaTheta                                1 mrad</div>
<div>/gate/source/PBS/setSigmaPhi                                    1 mrad</div>
<div>/gate/source/PBS/setEllipseXThetaEmittance          3 mm*mrad</div>
<div>/gate/source/PBS/setEllipseXThetaRotationNorm   positive</div>
<div>/gate/source/PBS/setEllipseYPhiEmittance               3 mm*mrad</div>
<div>/gate/source/PBS/setEllipseYPhiRotationNorm       positive</div>
<div>/gate/source/PBS/setParticleType                             proton</div>
<div>#/gate/source/PBS/setRotationAxis                          1 0 0</div>
<div>#/gate/source/PBS/setRotationAngle                       90 deg</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># INITIALISATION</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/run/initialize</div>
<div><br>
</div>
<div><br>
</div>
<div>#=====================================================</div>
<div># START BEAMS</div>
<div>#=====================================================</div>
<div><br>
</div>
<div>/gate/random/setEngineName MersenneTwister</div>
<div>/gate/random/setEngineSeed auto</div>
<div><br>
</div>
<div># When only using the EmCalculator actor for various configurations</div>
<div># (without caring of the dose deposition)</div>
<div>#/gate/application/setTotalNumberOfPrimaries 1</div>
<div># Otherwise</div>
<div>/gate/application/setTotalNumberOfPrimaries 1000000</div>
<div><br>
</div>
<div>/gate/application/start</div>
<br>
<p></p>
<p><br>
</p>
<div id="Signature">
<p>Sent from <a href="http://aka.ms/weboutlook" id="LPNoLP" tabindex="-1" disabled="true">
Outlook</a><br>
</p>
</div>
</div>
</body>
</html>