[Gate-users] dosimetry MCNP vs GATE

[Ashok Tiwari]

 Hi Mohamad, 
You may check our paper for the comparison of GATE vs MCNP results:
Ashok Tiwari, Stephen A. Graves, and John Sunderland (2020) The Impact of Tissue Type and Density on Dose Point Kernels for Patient-Specific Voxel-Wise Dosimetry: A Monte Carlo Investigation. Radiation Research In-Press.

Thanks, Ashok

    On Monday, 18 May 2020, 03:25:58 pm GMT-4, Eng. M Pharaon <pharaon_m at hotmail.com> wrote:  
 
  #yiv2861119228 P {margin-top:0;margin-bottom:0;}Dear Gate Users,Did any one try to compare between gate and mcnp in energy deposited in  region  and dose by regions ?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.I am using this code please take a look and share with me your experience 
#/control/execute verbose.mac
#/control/execute ../mac/visu.mac



#==============================
# GEOMETRY
#==============================

/gate/geometry/setMaterialDatabase ../data/GateMaterials.db

#----------------------
# World
#----------------------

/gate/world/geometry/setXLength 300 cm
/gate/world/geometry/setYLength 300 cm
/gate/world/geometry/setZLength 300 cm
/gate/world/setMaterial Air
/gate/world/vis/setVisible 0

# Body
/gate/world/daughters/name phantom
/gate/world/daughters/insert ImageNestedParametrisedVolume
/gate/phantom/geometry/setImage ../data/Segm.mhd

/gate/phantom/geometry/setRangeToMaterialFile ../data/range_geo.dat

#==============================
# PHYSICS
#==============================

#----------------------
# Physics List
#----------------------

/gate/physics/addPhysicsList emlivermore


#----------------------
# Cuts and Step limiters
#----------------------
/gate/physics/Electron/SetCutInRegion world 0.01 mm
/gate/physics/Gamma/SetCutInRegion world 0.01 mm

/gate/physics/Electron/SetCutInRegion phantom 0.001 mm
/gate/physics/Gamma/SetCutInRegion phantom 0.001 mm

/gate/physics/SetMaxStepSizeInRegion phantom 0.001 mm
/gate/physics/ActivateStepLimiter e-

/gate/physics/displayCuts


#########/gate/physics/addPhysicsList                                      emlivermore

#########/gate/physics/Gamma/SetCutInRegion                                world 10.0 mm
##########/gate/physics/Electron/SetCutInRegion                             world 10.0 mm
##########/gate/physics/Positron/SetCutInRegion                             world 10.0 mm

##########/gate/physics/Gamma/SetCutInRegion                                phantom 0.1 mm
##########/gate/physics/Electron/SetCutInRegion                             phantom 0.1 mm
##########/gate/physics/Positron/SetCutInRegion                             phantom 0.1 mm


#==============================
# DETECTORS
#==============================

# It is possible to save a DoseActor as a 3D matrix.

/gate/actor/addActor                                              DoseActor  dose3D
/gate/actor/dose3D/attachTo                                       phantom
/gate/actor/dose3D/stepHitType                                    random
/gate/actor/dose3D/setResolution                                  1346 693 147
/gate/actor/dose3D/enableDose                                     true
/gate/actor/dose3D/save                                           ../output/output.mhd
/gate/actor/dose3D/inputDoseByRegions                             ../data/Segm.mhd
/gate/actor/dose3D/outputDoseByRegions                            ../output/DoseByRegions.txt

#==============================
# INITIALISATION
#==============================

/gate/run/initialize

#==============================
# BEAM
#==============================
# Source

/gate/source/addSource                                            phantom voxel
# DECLARATION THAT THE VOXELIZED SOURCE WILL BE ENTERED USING IMAGE DATA
/gate/source/phantom/reader/insert                              image

/gate/source/phantom/imageReader/translator/insert              range
/gate/source/phantom/imageReader/rangeTranslator/readTable      ../data/ActivityRange.dat
/gate/source/phantom/imageReader/rangeTranslator/describe       1

# THE FOLLOWING LINE ALLOWS YOU TO INSERT THE IMAGE DATA USING THE APPROPRIATE EXTENSION FILE
/gate/source/phantom/imageReader/readFile                       ../data/Segm.mhd

# THE DEFAULT POSITION OF THE VOXELIZED SOURCE IS IN THE 1ST QUARTER
# SO THE VOXELIZED SOURCE HAS TO BE SHIFTED OVER HALF ITS DIMENSION IN THE NEGATIVE DIRECTION ON EACH AXIS
/gate/source/phantom/setPosition                                 -673  -346.5  -73.5 mm


# THE FOLLOWING LINES CHARACTERIZE THE SIZE (NO DIFFERENCE WITH AN ANALYTICAL SOURCE)
/gate/source/phantom/gps/particle gamma
/gate/source/phantom/gps/energytype Mono
/gate/source/phantom/gps/monoenergy 0.1 MeV
/gate/source/phantom/gps/angtype iso

#==============================
#==============================
# START SIMULATION
#==============================

/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto


/gate/application/setTotalNumberOfPrimaries 1000
/gate/application/start
==============================================================#activity range file1
135   135   1
===================
Best Regards
                     Mohamad Pharaon
                           Nuclear Engineer 
        
_______________________________________________
Gate-users mailing list
Gate-users at lists.opengatecollaboration.org
http://lists.opengatecollaboration.org/mailman/listinfo/gate-users  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200518/54722f40/attachment-0001.html>