<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>