[Gate-users] Confining heart source within tessalated STL body - loop 100000 errors.

[maxime]

Hi Steven,

can you tell us more about the errors you get ?

  Maxime

> On 19 May 2020, at 10:19, Steven Marsh <steven.marsh at canterbury.ac.nz> wrote:
> 
> Hi Gate users…
>  
> I have successfully confined a STL defined volume source (heart) containing Tc-99m and I can place that within a body (src_Heart in code below) so long as “body” is a simple shape( a box in macro1). When I attempt to place the same confined src_Heart  within a “body” described by a tessellated surface, I get loop 100000 errors and I can visualise dots indicating the source positions for decay are not located within the heart (macro2). Below I include the important parts of my macros which attempt to do this – note for any run I only attempt to run either macro2 OR macro2 – not both.
>  
> Is it possible to encapsulate a confined source volume within a tessalated surface (body in my case)?
>  
> Any ideas much appreciated
>  
> Best wishes
> Steve
> -----------------------------------------------------------------------------------------------------------------------
> Macro1)  This works well…body is a simple box
> gate/world/daughters/name body
> /gate/world/daughters/insert box
> /gate/body/geometry/setXLength 20 cm
> /gate/body/geometry/setYLength 20 cm
> /gate/body/geometry/setZLength 20 cm
> /gate/body/vis/setColor blue
> /gate/body/vis/forceWireframe
> /gate/body/setMaterial Air
>  
> ## ---IMPORT STL FOR HEART AND POSITION IN BODY
> /gate/body/daughters/name heart
> /gate/body/daughters/insert tessellated
> /gate/heart/geometry/setPathToSTLFile data/patientData/STLFiles/Heart_Only.stl
> /gate/heart/setMaterial        Heart
>  
> --------------------------------------------------------------------------------------------------------------------------- END Macro1
> Macro2)  This doesn’t work – body is a tessalated surface
> gate/world/daughters/name body
> /gate/world/daughters/insert tessellated
> /gate/body/geometry/setPathToSTLFile data/patientData/STLFiles/Body.stl
> /gate/body/setMaterial        Body
>  
> #-----IMPORT STL FOR HEART AND POSITION IN BODY 
> /gate/body/daughters/name heart
> /gate/body/daughters/insert tessellated
> /gate/heart/geometry/setPathToSTLFile data/patientData/STLFiles/Heart_Only.stl
> #/gate/heart/geometry/setPathToSTLFile data/patientData/STLFiles/Heart.stl
> /gate/heart/setMaterial       Heart
> ----------------------------------------------------------------------------------------------------------------------------- End Macro2
> Macro3) the following is used for either of the above to confine the src_Heart volume to the heart STL shape
> /gate/source/addSource src_Heart gps
> /gate/source/src_Heart/gps/particle gamma
> /gate/source/src_Heart/gps/pos/type Volume
> /gate/source/src_Heart/gps/pos/shape Para
> /gate/source/src_Heart/gps/pos/centre 0 0 0 cm 
> /gate/source/src_Heart/gps/pos/halfx 10 cm  
> /gate/source/src_Heart/gps/pos/halfy 10 cm  
> /gate/source/src_Heart/gps/pos/halfz 10 cm 
> /gate/source/src_Heart/gps/ang/type iso
> #Define the source as a Tc99m source
> /gate/source/src_Heart/gps/energytype UserSpectrum
> /gate/source/src_Heart/gps/setSpectrumFile data/energy_spectrum_Tc99m.txt
> /gate/source/src_Heart/setForcedHalfLife  21624.12 s  #6.0 hrs
> /gate/source/src_Heart/setForcedUnstableFlag true
> /gate/source/src_Heart/setForcedHalfLife 21624.12 s
> /gate/source/src_Heart/setActivity {ACTIVITY} MBq
> # Now confine src_Heart to the shape of the heart.  
> /gate/source/src_Heart/gps/pos/confine heart
> ----------------------------------------------------------------------------------------------------------------------------- End Macro3
>  
> _______________________________________________
> 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/20200519/a8314a77/attachment.html>