[Gate-users] Why do output file positron nuclides get deposited in the air

hao qiushi haoqiushi at outlook.com
Sat Sep 4 13:40:16 CEST 2021


I simulated the process of radiotherapy in mice by using mouse CT files. In the figure, the position of 14mm is the position where the beam enters the mouse body, and before that, it was air. What I don't understand is why some positron nuclides have higher activity values before they are injected into the mouse body.Can anyone help me.
[cid:8fdc6d6c-ad29-4e5a-a524-098d1613bd66]
My macro file:

/control/execute /home/bishe3/Desktop/gate/c200mouse/mac/verbose.mac
#/control/execute mac/visu.mac

/gate/geometry/setMaterialDatabase       /home/bishe3/Desktop/gate/c200mouse/data/GateMaterials.db

# World =======================

/gate/world/geometry/setXLength 2 m
/gate/world/geometry/setYLength 2 m
/gate/world/geometry/setZLength 2 m
/gate/world/setMaterial         Air

# geomerty ===================

/gate/HounsfieldMaterialGenerator/SetMaterialTable                  /home/bishe3/Desktop/gate/c200mouse/data/Schneider2000MaterialsTable.txt
/gate/HounsfieldMaterialGenerator/SetDensityTable                   /home/bishe3/Desktop/gate/c200mouse/data/Schneider2000DensitiesTable.txt
/gate/HounsfieldMaterialGenerator/SetDensityTolerance               0.1 g/cm3
/gate/HounsfieldMaterialGenerator/SetOutputMaterialDatabaseFilename /home/bishe3/Desktop/gate/c200mouse/data/patient-HUmaterials.db
/gate/HounsfieldMaterialGenerator/SetOutputHUMaterialFilename       /home/bishe3/Desktop/gate/c200mouse/data/patient-HU2mat.dat
/gate/HounsfieldMaterialGenerator/Generate

/gate/world/daughters/name                        mouse
/gate/world/daughters/insert                      ImageNestedParametrisedVolume
/gate/mouse/geometry/setImage                      /home/bishe3/Desktop/gate/c200mouse/data/ct_data/ct_380x992x208.hdr
/gate/geometry/setMaterialDatabase                /home/bishe3/Desktop/gate/c200mouse/data/patient-HUmaterials.db
/gate/mouse/geometry/setRangeToMaterialFile        /home/bishe3/Desktop/gate/c200mouse/data/patient-HU2mat.dat
/gate/mouse/geometry/setHUToMaterialFile           /home/bishe3/Desktop/gate/c200mouse/data/patient-HU2mat.dat

# PHYSICS=========

/gate/physics/addPhysicsList            QGSP_BIC_HP_EMZ

/gate/physics/Gamma/SetCutInRegion      world 10 mm
/gate/physics/Electron/SetCutInRegion   world 10 mm
/gate/physics/Positron/SetCutInRegion   world 10 mm

/gate/physics/Gamma/SetCutInRegion      mouse 1 mm
/gate/physics/Electron/SetCutInRegion   mouse 1 mm
/gate/physics/Positron/SetCutInRegion   mouse 1 mm

# output===================

/gate/actor/addActor                               DoseActor  doseDistribution
/gate/actor/doseDistribution/save                  /home/bishe3/Desktop/gate/c200mouse/output/output-C200.txt
/gate/actor/doseDistribution/attachTo         mouse
#/gate/actor/doseDistribution/setSize              38 99.2 20.8 mm
/gate/actor/doseDistribution/stepHitType           random
/gate/actor/doseDistribution/setPosition           0 0 0 cm
/gate/actor/doseDistribution/setResolution         380 992 208
#/gate/actor/doseDistribution/enableEdep            true
/gate/actor/doseDistribution/enableDose            true
/gate/actor/doseDistribution/enableNumberOfHits    true

/gate/actor/addActor                                   DoseActor  doseDistribution_txt
/gate/actor/doseDistribution_txt/save                  /home/bishe3/Desktop/gate/c200mouse/output/output-C200_y.txt
/gate/actor/doseDistribution_txt/attachTo              mouse
#/gate/actor/doseDistribution_txt/setSize              38 99.2 20.8 mm
/gate/actor/doseDistribution_txt/stepHitType           random
/gate/actor/doseDistribution_txt/setPosition           0 0 0 cm
/gate/actor/doseDistribution_txt/setResolution         1 992 1
#/gate/actor/doseDistribution_txt/enableEdep            true
/gate/actor/doseDistribution_txt/enableDose            true
/gate/actor/doseDistribution_txt/enableNumberOfHits    true

# ProductionAndStoppingActor==================

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_11
/gate/actor/prdandstp_11/save               /home/bishe3/Desktop/gate/c200mouse/output/product_stop11_3D.txt
/gate/actor/prdandstp_11/attachTo           mouse
/gate/actor/prdandstp_11/setResolution      380 992 208
/gate/actor/prdandstp_11/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_10
/gate/actor/prdandstp_10/save               /home/bishe3/Desktop/gate/c200mouse/output/product_stop10_3D.txt
/gate/actor/prdandstp_10/attachTo           mouse
/gate/actor/prdandstp_10/setResolution      380 992 208
/gate/actor/prdandstp_10/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_15
/gate/actor/prdandstp_15/save               /home/bishe3/Desktop/gate/c200mouse/output/product_stop15_3D.txt
/gate/actor/prdandstp_15/attachTo           mouse
/gate/actor/prdandstp_15/setResolution      380 992 208
/gate/actor/prdandstp_15/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_11_txt
/gate/actor/prdandstp_11_txt/save               /home/bishe3/Desktop/gate/c200mouse/output/product_stop11.txt
/gate/actor/prdandstp_11_txt/attachTo           mouse
/gate/actor/prdandstp_11_txt/setResolution      1 992 1
/gate/actor/prdandstp_11_txt/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_10_txt
/gate/actor/prdandstp_10_txt/save               /home/bishe3/Desktop/gate/c200mouse/output/product_stop10.txt
/gate/actor/prdandstp_10_txt/attachTo           mouse
/gate/actor/prdandstp_10_txt/setResolution      1 992 1
/gate/actor/prdandstp_10_txt/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_15_txt
/gate/actor/prdandstp_15_txt/save               /home/bishe3/Desktop/gate/c200mouse/output/product_stop15.txt
/gate/actor/prdandstp_15_txt/attachTo           mouse
/gate/actor/prdandstp_15_txt/setResolution      1 992 1
/gate/actor/prdandstp_15_txt/stepHitType        random

/gate/actor/addActor                      SimulationStatisticActor stat
/gate/actor/stat/save                     /home/bishe3/Desktop/gate/c200mouse/output/stat-C100.txt
/gate/actor/stat/saveEveryNSeconds        10

# filtter====================

/gate/actor/prdandstp_11/addFilter                     particleFilter
/gate/actor/prdandstp_11/particleFilter/addParticleZ      6
/gate/actor/prdandstp_11/particleFilter/addParticleA      11

/gate/actor/prdandstp_10/addFilter                     particleFilter
/gate/actor/prdandstp_10/particleFilter/addParticleZ      6
/gate/actor/prdandstp_10/particleFilter/addParticleA      10

/gate/actor/prdandstp_15/addFilter                     particleFilter
/gate/actor/prdandstp_15/particleFilter/addParticleZ      8
/gate/actor/prdandstp_15/particleFilter/addParticleA      15

/gate/actor/prdandstp_11_txt/addFilter                     particleFilter
/gate/actor/prdandstp_11_txt/particleFilter/addParticleZ      6
/gate/actor/prdandstp_11_txt/particleFilter/addParticleA      11

/gate/actor/prdandstp_10_txt/addFilter                     particleFilter
/gate/actor/prdandstp_10_txt/particleFilter/addParticleZ      6
/gate/actor/prdandstp_10_txt/particleFilter/addParticleA      10

/gate/actor/prdandstp_15_txt/addFilter                     particleFilter
/gate/actor/prdandstp_15_txt/particleFilter/addParticleZ      8
/gate/actor/prdandstp_15_txt/particleFilter/addParticleA      15
# INITIALISATION============

/gate/run/initialize
#/gate/physics/processList Available
#/gate/physics/processList Enabled

# BEAMS==================

/gate/source/addSource                         mybeam  PencilBeam
#/gate/source/mybeam/setParticleType            proton
/gate/source/mybeam/setParticleType            GenericIon
/gate/source/mybeam/setIonProperties           6 12 6 0
/gate/source/mybeam/setEnergy                  1080 MeV
/gate/source/mybeam/setSigmaEnergy             2 MeV
/gate/source/mybeam/setPosition                0 -70 1.6 mm
/gate/source/mybeam/setRotationAxis            1 0 0
/gate/source/mybeam/setRotationAngle           -90 deg
/gate/source/mybeam/setSigmaX                  2 mm
/gate/source/mybeam/setSigmaY                  2 mm
/gate/source/mybeam/setSigmaTheta              0.01 mrad
/gate/source/mybeam/setSigmaPhi                0.01 mrad
/gate/source/mybeam/setEllipseXThetaEmittance  0.01 mm*mrad
/gate/source/mybeam/setEllipseYPhiEmittance    0.01 mm*mrad

# JamesRandom Ranlux64 MersenneTwister
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto

/gate/application/noGlobalOutput
#/gate/application/readTimeSlicesIn               data/timeslice.dat
#/gate/application/readNumberOfPrimariesInAFile   data/primary.dat
/gate/application/setTotalNumberOfPrimaries       1000000
/gate/application/start
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210904/ad473db9/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 111463 bytes
Desc: image.png
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210904/ad473db9/attachment-0001.png>


More information about the Gate-users mailing list