[Gate-users] Filling an imported STL volume with a radionuclide

[Steven Marsh]

Hi

Perhaps one of you can assist me...

My aim is to import, into my simulation, anatomy (in stl format - though I can change this if necessary) and fill the heart with Tc-99m. I want to assign a total amount of activity to the heart volume - i.e. ideally GATE will evenly spread the activity within the heart.

I have imported the anatomy - via the following:

/gate/phantom/daughters/name heart
/gate/phantom/daughters/insert tessellated
/gate/heart/geometry/setPathToSTLFile data/Heart.stl
/gate/heart/setMaterial        Heart


And I have hopefully filled the heart with Tc-99m via the following:

/gate/source/addSource source voxel
/gate/source/source/gps/particle gamma
/gate/source/source/gps/type Volume
/gate/source/source/gps/angtype iso
/gate/source/source/gps/energytype UserSpectrum
/gate/source/source/gps/setSpectrumFile data/energy_spectrum_Tc99m.txt
/gate/source/source/setForcedHalfLife  21624.12 s  #6.0 hrs
/gate/source/source/setActivity 0.001 MBq
/gate/source/source/attachTo heart

This (in addition to other code not given above) runs.

My questions are:

Firstly - is the above doing what I am hoping it will do?  Will this spread 0.001MBq evenly throughout the heart volume?

Secondly - my code is running - though it is taking a very long time and is reporting it is up to run 4973100 with an elapsed time of 0.01 or so- does this relate to the release of primary photons - ie appx 5e6 photons have been released and tracked??

Any advice much appreciated

Best wishes
Steve
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200512/4dac180e/attachment.html>