<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;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
<p>Dears Gate users,</p>
<p><br>
</p>
<p>I am running an example as exercise for basic dose calculations. It consists of a cylindrical phantom (patient) contains a target treated by
<b>ONE</b> proton PBS. Then I summed all the resulted dose in phantom and target separately (phantom 18.81 cGy, target 48.3 cGy). I run the example again with
<b>TWO</b> beams and I got the same numbers for both dose actors. The number of primaries was the same in both simulations. </p>
<p>Is the primaries are split on the two beams by default? Am I on the correct way of dose calc.?</p>
<p>Here a copy of the main macro.</p>
<p><br>
</p>
<p>Thanks in advance,</p>
<p>Awad</p>
<p><br>
</p>
<p></p>
<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># WATERCYLINDER</div>
<div>/gate/world/daughters/name water_cylinder</div>
<div>/gate/world/daughters/insert cylinder</div>
<div>/gate/water_cylinder/setMaterial Water</div>
<div>/gate/water_cylinder/geometry/setRmin 0 cm</div>
<div>/gate/water_cylinder/geometry/setRmax 15 cm</div>
<div>/gate/water_cylinder/geometry/setHeight 40 cm</div>
<div>/gate/water_cylinder/placement/setRotationAxis 1 0 0</div>
<div>/gate/water_cylinder/placement/setRotationAngle 90 deg</div>
<div>/gate/water_cylinder/placement/setTranslation 0.0 0.0 0.0 cm</div>
<div>/gate/water_cylinder/vis/setColor blue</div>
<div>/gate/water_cylinder/vis/forceWireframe </div>
<div><br>
</div>
<div>#Tumor</div>
<div>/gate/water_cylinder/daughters/name target</div>
<div>/gate/water_cylinder/daughters/insert sphere</div>
<div>/gate/target/setMaterial Water</div>
<div>/gate/target/geometry/setRmin 0 cm</div>
<div>/gate/target/geometry/setRmax 3 cm</div>
<div>/gate/target/placement/setTranslation 0.0 0.0 0.0 cm</div>
<div>/gate/target/vis/setColor red</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_cylinder 0.1 mm</div>
<div>/gate/physics/Electron/SetCutInRegion water_cylinder 0.1 mm</div>
<div>/gate/physics/Positron/SetCutInRegion water_cylinder 0.1 mm</div>
<div><br>
</div>
<div>/gate/physics/Gamma/SetCutInRegion target 0.1 mm</div>
<div>/gate/physics/Electron/SetCutInRegion target 0.1 mm</div>
<div>/gate/physics/Positron/SetCutInRegion target 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/phantom2Fields.txt</div>
<div>/gate/actor/dose/attachTo <span class="Apple-tab-span" style="white-space:pre">
</span> water_cylinder</div>
<div>/gate/actor/dose/stepHitType random</div>
<div>/gate/actor/dose/setPosition 0 0 0 cm</div>
<div>/gate/actor/dose/setResolution 100 100 300</div>
<div>/gate/actor/dose/setSize 8 8 30 cm</div>
<div>/gate/actor/dose/saveEveryNSeconds 60</div>
<div> /gate/actor/dose/enableEdep false</div>
<div>/gate/actor/dose/enableUncertaintyEdep false</div>
<div>/gate/actor/dose/enableDose true</div>
<div>/gate/actor/dose/enableUncertaintyDose true</div>
<div>/gate/actor/dose/enableNumberOfHits false</div>
<div><br>
</div>
<div>/gate/actor/addActor DoseActor dose1</div>
<div>/gate/actor/dose1/save output/target2Fields.txt</div>
<div>/gate/actor/dose1/attachTo <span class="Apple-tab-span" style="white-space:pre">
</span> target</div>
<div>/gate/actor/dose1/stepHitType random</div>
<div>/gate/actor/dose1/setPosition 0 0 0 cm</div>
<div>/gate/actor/dose1/setResolution 100 100 100</div>
<div>/gate/actor/dose1/setSize 6 6 6 cm</div>
<div>/gate/actor/dose1/saveEveryNSeconds 60 </div>
<div><span style="font-size: 12pt;">/gate/actor/dose1/enableEdep false</span></div>
<div>/gate/actor/dose1/enableUncertaintyEdep false</div>
<div>/gate/actor/dose1/enableDose true</div>
<div>/gate/actor/dose1/enableUncertaintyDose true</div>
<div>/gate/actor/dose1/enableNumberOfHits false</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 -50 cm</div>
<div>/gate/source/PBS/setEnergy 150 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 5 mm*mrad</div>
<div>/gate/source/PBS/setEllipseXThetaRotationNorm positive</div>
<div>/gate/source/PBS/setEllipseYPhiEmittance 5 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 0 deg</div>
<div><br>
</div>
<div>/gate/source/addSource PBS2 PencilBeam</div>
<div>/gate/source/PBS2/setPosition 0 0 50 cm</div>
<div>/gate/source/PBS2/setEnergy 150 MeV</div>
<div>/gate/source/PBS2/setSigmaEnergy 0.1 MeV</div>
<div>/gate/source/PBS2/setSigmaX 3 mm</div>
<div>/gate/source/PBS2/setSigmaY 3 mm</div>
<div>/gate/source/PBS2/setSigmaTheta 1 mrad</div>
<div>/gate/source/PBS2/setSigmaPhi 1 mrad</div>
<div>/gate/source/PBS2/setEllipseXThetaEmittance 5 mm*mrad</div>
<div>/gate/source/PBS2/setEllipseXThetaRotationNorm positive</div>
<div>/gate/source/PBS2/setEllipseYPhiEmittance 5 mm*mrad</div>
<div>/gate/source/PBS2/setEllipseYPhiRotationNorm positive</div>
<div>/gate/source/PBS2/setParticleType proton</div>
<div>/gate/source/PBS2/setRotationAxis 1 0 0</div>
<div>/gate/source/PBS2/setRotationAngle 180 deg</div>
<div><br>
</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">Outlook</a><br>
</p>
</div>
</div>
</body>
</html>