<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 style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Dear Gate Users,</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
Did any one try to compare between gate and mcnp in energy deposited in region and dose by regions ?</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
i am finding a lot of differences in the deposited energy value - that calculated in MCNP compered with Gate - when i am using energy higher than 0.3 Mev.</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
I am using this code please take a look and share with me your experience </div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>#/control/execute verbose.mac<br>
</span>
<div>#/control/execute ../mac/visu.mac<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>#==============================<br>
</div>
<div># GEOMETRY<br>
</div>
<div>#==============================<br>
</div>
<div><br>
</div>
<div>/gate/geometry/setMaterialDatabase ../data/GateMaterials.db<br>
</div>
<div><br>
</div>
<div>#----------------------<br>
</div>
<div># World<br>
</div>
<div>#----------------------<br>
</div>
<div><br>
</div>
<div>/gate/world/geometry/setXLength 300 cm<br>
</div>
<div>/gate/world/geometry/setYLength 300 cm<br>
</div>
<div>/gate/world/geometry/setZLength 300 cm<br>
</div>
<div>/gate/world/setMaterial Air<br>
</div>
<div>/gate/world/vis/setVisible 0<br>
</div>
<div><br>
</div>
<div># Body<br>
</div>
<div>/gate/world/daughters/name phantom<br>
</div>
<div>/gate/world/daughters/insert ImageNestedParametrisedVolume<br>
</div>
<div>/gate/phantom/geometry/setImage ../data/Segm.mhd<br>
</div>
<div><br>
</div>
<div>/gate/phantom/geometry/setRangeToMaterialFile ../data/range_geo.dat<br>
</div>
<div><br>
</div>
<div>#==============================<br>
</div>
<div># PHYSICS<br>
</div>
<div>#==============================<br>
</div>
<div><br>
</div>
<div>#----------------------<br>
</div>
<div># Physics List<br>
</div>
<div>#----------------------<br>
</div>
<div><br>
</div>
<div>/gate/physics/addPhysicsList emlivermore<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>#----------------------<br>
</div>
<div># Cuts and Step limiters<br>
</div>
<div>#----------------------<br>
</div>
<div>/gate/physics/Electron/SetCutInRegion world 0.01 mm<br>
</div>
<div>/gate/physics/Gamma/SetCutInRegion world 0.01 mm<br>
</div>
<div><br>
</div>
<div>/gate/physics/Electron/SetCutInRegion phantom 0.001 mm<br>
</div>
<div>/gate/physics/Gamma/SetCutInRegion phantom 0.001 mm<br>
</div>
<div><br>
</div>
<div>/gate/physics/SetMaxStepSizeInRegion phantom 0.001 mm<br>
</div>
<div>/gate/physics/ActivateStepLimiter e-<br>
</div>
<div><br>
</div>
<div>/gate/physics/displayCuts<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>#########/gate/physics/addPhysicsList emlivermore<br>
</div>
<div><br>
</div>
<div>#########/gate/physics/Gamma/SetCutInRegion world 10.0 mm<br>
</div>
<div>##########/gate/physics/Electron/SetCutInRegion world 10.0 mm<br>
</div>
<div>##########/gate/physics/Positron/SetCutInRegion world 10.0 mm<br>
</div>
<div><br>
</div>
<div>##########/gate/physics/Gamma/SetCutInRegion phantom 0.1 mm<br>
</div>
<div>##########/gate/physics/Electron/SetCutInRegion phantom 0.1 mm<br>
</div>
<div>##########/gate/physics/Positron/SetCutInRegion phantom 0.1 mm<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>#==============================<br>
</div>
<div># DETECTORS<br>
</div>
<div>#==============================<br>
</div>
<div><br>
</div>
<div># It is possible to save a DoseActor as a 3D matrix.<br>
</div>
<div><br>
</div>
<div>/gate/actor/addActor DoseActor dose3D<br>
</div>
<div>/gate/actor/dose3D/attachTo phantom<br>
</div>
<div>/gate/actor/dose3D/stepHitType random<br>
</div>
<div>/gate/actor/dose3D/setResolution 1346 693 147<br>
</div>
<div>/gate/actor/dose3D/enableDose true<br>
</div>
<div>/gate/actor/dose3D/save ../output/output.mhd<br>
</div>
<div>/gate/actor/dose3D/inputDoseByRegions ../data/Segm.mhd<br>
</div>
<div>/gate/actor/dose3D/outputDoseByRegions ../output/DoseByRegions.txt<br>
</div>
<div><br>
</div>
<div>#==============================<br>
</div>
<div># INITIALISATION<br>
</div>
<div>#==============================<br>
</div>
<div><br>
</div>
<div>/gate/run/initialize<br>
</div>
<div><br>
</div>
<div>#==============================<br>
</div>
<div># BEAM<br>
</div>
<div>#==============================<br>
</div>
<div># Source<br>
</div>
<div><br>
</div>
<div>/gate/source/addSource phantom voxel<br>
</div>
<div># DECLARATION THAT THE VOXELIZED SOURCE WILL BE ENTERED USING IMAGE DATA<br>
</div>
<div>/gate/source/phantom/reader/insert image<br>
</div>
<div><br>
</div>
<div>/gate/source/phantom/imageReader/translator/insert range<br>
</div>
<div>/gate/source/phantom/imageReader/rangeTranslator/readTable ../data/ActivityRange.dat<br>
</div>
<div>/gate/source/phantom/imageReader/rangeTranslator/describe 1<br>
</div>
<div><br>
</div>
<div># THE FOLLOWING LINE ALLOWS YOU TO INSERT THE IMAGE DATA USING THE APPROPRIATE EXTENSION FILE<br>
</div>
<div>/gate/source/phantom/imageReader/readFile ../data/Segm.mhd<br>
</div>
<div><br>
</div>
<div># THE DEFAULT POSITION OF THE VOXELIZED SOURCE IS IN THE 1ST QUARTER<br>
</div>
<div># SO THE VOXELIZED SOURCE HAS TO BE SHIFTED OVER HALF ITS DIMENSION IN THE NEGATIVE DIRECTION ON EACH AXIS<br>
</div>
<div>/gate/source/phantom/setPosition -673 -346.5 -73.5 mm<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div># THE FOLLOWING LINES CHARACTERIZE THE SIZE (NO DIFFERENCE WITH AN ANALYTICAL SOURCE)<br>
</div>
<div>/gate/source/phantom/gps/particle gamma<br>
</div>
<div>/gate/source/phantom/gps/energytype Mono<br>
</div>
<div>/gate/source/phantom/gps/monoenergy 0.1 MeV<br>
</div>
<div>/gate/source/phantom/gps/angtype iso<br>
</div>
<div><br>
</div>
<div>#==============================<br>
</div>
<div>#==============================<br>
</div>
<div># START SIMULATION<br>
</div>
<div>#==============================<br>
</div>
<div><br>
</div>
<div>/gate/random/setEngineName MersenneTwister<br>
</div>
<div>/gate/random/setEngineSeed auto<br>
</div>
<div><br>
</div>
<div><br>
</div>
<div>/gate/application/setTotalNumberOfPrimaries 1000<br>
</div>
<span>/gate/application/start</span><br>
</div>
<div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
==============================================================</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
#activity range file</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>1<br>
</span><span>135 135 1</span><br>
</div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>===================</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><br>
</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span>Best Regards</span></div>
<div style="font-family: Calibri, Helvetica, sans-serif; font-size: 12pt; color: rgb(0, 0, 0);">
<span><br>
</span></div>
<div id="Signature"> <font size="2" style="font-size:10pt"> <font color="#666666"><strong>Mohamad Pharaon</strong></font><br>
<font color="#666666"><strong> Nuclear Engineer </strong>
</font><br>
<font color="#666666"><strong> </strong></font></font><br>
</div>
</div>
</body>
</html>