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

[maxime]

Dear Steven,

First, you are using correctly the stl to create a tesselated geometry of a heart.

But you are using a voxel source. A voxel source is intended for voxel geometries.

Instead here is what you should do:
- use a GPS source with, shape, volume, energy, etc… For example a simple cube
- attach it to the heart as you did (you can skip this if your heart will not move during the simulation)
- place it correctly with respect to the STL coordinates (check with visu like OGL)
- confine the source to heart: this final step is a trick to restrict the shape (cube) to the exact heart STL shape. So it is important that your simple shape (Cube, Sphere) is as close as possible but bigger than your heart.

For your second question, it could take a lot of time due to a side effect of your not correct source definition.

Best regards,
  Maxime Chauvin

> On 12 May 2020, at 06:24, Steven Marsh <steven.marsh at canterbury.ac.nz> wrote:
> 
> 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
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org <mailto:Gate-users at lists.opengatecollaboration.org>
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users <http://lists.opengatecollaboration.org/mailman/listinfo/gate-users>

> On 12 May 2020, at 06:24, Steven Marsh <steven.marsh at canterbury.ac.nz> wrote:
> 
> 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
> _______________________________________________
> Gate-users mailing list
> Gate-users at lists.opengatecollaboration.org <mailto:Gate-users at lists.opengatecollaboration.org>
> http://lists.opengatecollaboration.org/mailman/listinfo/gate-users <http://lists.opengatecollaboration.org/mailman/listinfo/gate-users>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200512/18e72dc7/attachment.html>