[Gate-users] Dose in Air

[Maikol Salas Ramirez]

Dear Sara

It looks like the cross section for air in Geant4 and EGSnrc are diferent.
EGSnrc has a good perform at low energy, then it is normal that for high
energy both software match.

Another different in the stepfuntion, I mean the algorithm is different.
Now as Tony said, you can change the step size, it is a tricky way to
chance the cross section for an specific material.

If you want to go in deep, can be good check the physic manual of each code
and ask directly in the Geant4 and EGSnrc blogs.

Best regards
Maikol


El mié., 31 oct. 2018 a las 11:19, Sara Neira Castro (<
sara.neira.castro at hotmail.com>) escribió:

> Dear Gate users,
>
> We have recently started to use Gate to simulate dose deposition in
> phantoms/CTs, and we have found some inconsistencies in doses deposited in
> low  density regions (mainly air) when crosschecking with EGSnrc: doses
> calculated with Gate are much lower than doses calculated with EGSnrc. In
> the  figure we report PDDs for different energies (photons) in a phantom
> formed by several layers of material (Water-Air-Water-Lung-Water). Material
> definitions are identical for Gate and EGS. Doses in water are identical,
> and doses in lung are very similar, but doses in air differ for low
> energies.  Interestingly, this effect seems to vanish with increasing
> energy (>1 MeV).
>
> We have checked the input thousands of times, and yet we cannot find any
> explanation. We have tried different Gate physics packages
> (emstandard_opt3,  emstandard_opt4,penelope ...) and cutoffs (typically 10
> keV or below), but the systematic difference is still there.  Is there any
> known limitation of Gate for simulating dose deposition of low energy
> photons in low density media?, or is it some kind of error in our
> simulation?. We would be grateful for any suggestion. Below are the main
> features of the Gate code:
>
> # GEOMETRY
> /gate/world/daughters/name     cajaAgua
> /gate/world/daughters/insert   box
> /gate/cajaAgua/geometry/setXLength 29 cm
> /gate/cajaAgua/geometry/setYLength 29 cm
> /gate/cajaAgua/geometry/setZLength 30. cm
> /gate/cajaAgua/setMaterial Water
>
> /gate/cajaAgua/daughters/name     aire1
> /gate/cajaAgua/daughters/insert   box
> /gate/aire1/geometry/setXLength 29 cm
> /gate/aire1/geometry/setYLength 29 cm
> /gate/aire1/geometry/setZLength 3. cm
> /gate/aire1/setMaterial Air
> /gate/aire1/placement/setTranslation 0. 0. -10.5 cm
>
> /gate/cajaAgua/daughters/name     aire2
> /gate/cajaAgua/daughters/insert   box
> /gate/aire2/geometry/setXLength 29 cm
> /gate/aire2/geometry/setYLength 29 cm
> /gate/aire2/geometry/setZLength 3. cm
> /gate/aire2/setMaterial Lung
> /gate/aire2/placement/setTranslation 0 0. -4.5 cm
>
>
> The source consists in a plane of monoenergetic photons which  is located
> on the front face of the geometry.  The code was run multiple times for
> several energies (200, 511, 1000 and 3000 keV).
>
> # SOURCE
> /gate/source/addSource fuente gps
> /gate/source/fuente/gps/particle gamma
> /gate/source/fuente/gps/energytype Mono
> /gate/source/fuente/gps/monoenergy 200 keV
>
> /gate/source/fuente/gps/type Plane
> /gate/source/fuente/gps/shape Rectangle
> /gate/source/fuente/gps/halfx 2.5 cm
> /gate/source/fuente/gps/halfy 2.5 cm
> /gate/source/fuente/gps/centre 0. 0. -15 cm
>
> /gate/source/fuente/gps/angtype iso
> /gate/source/fuente/gps/mintheta 0. deg
> /gate/source/fuente/gps/maxtheta 0. deg
> /gate/source/fuente/gps/minphi 0. deg
> /gate/source/fuente/gps/maxphi 360. deg
>
>
> Thank you very much in advance,
> Sara
>
> _______________________________________________
> 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/20181105/5f38e4d9/attachment.html>