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

[Steven Marsh]

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

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200519/a52b47dc/attachment.html>