#                                                          #
#                      I 125 - SL-125                      #
#                                                          #
#                                                          #
#                Latest version - May 2021                 #
#                  Author: Audran Poher                    #
#              Contact - fberumenm@gmail.com               #
#       https://doi.org/10.1016/j.brachy.2020.12.007       #
#                                                          #


##### MATERIAL DEFINITION #####
sv:Ma/MatSilverRod/Components = 1 "Silver"
uv:Ma/MatSilverRod/Fractions  = 1 1
d:Ma/MatSilverRod/Density     = 10.5 g/cm3

sv:Ma/MatRadioactiveLayer/Components = 1 "Iodine"
uv:Ma/MatRadioactiveLayer/Fractions  = 1 1
d:Ma/MatRadioactiveLayer/Density     = 4.93 g/cm3

##### SOURCE #####
s:So/ActiveSource/Type                       = "Volumetric"
s:So/ActiveSource/Component                  = "ActiveSourceGeometry"
sc:So/ActiveSource/ActiveMaterial            = "MatRadioactiveLayer"
s:So/ActiveSource/BeamParticle               = "gamma"
b:So/ActiveSource/RecursivelyIncludeChildren = "T"
i:So/ActiveSource/MaxNumberOfPointsToSample  = 1000000000
s:So/ActiveSource/BeamEnergySpectrumType     = "Discrete"
ic:So/ActiveSource/NumberOfHistoriesInRun    = 10

##### SPECTRUM ######
dv:So/ActiveSource/BeamEnergySpectrumValues = 7 3.77 27.202 27.472 30.944 30.995 31.704 35.3922 keV
uv:So/ActiveSource/BeamEnergySpectrumWeightsUnscaled = 7 0.149 0.401 0.740 0.0683 0.132 0.0380 0.0668
uv:So/ActiveSource/BeamEnergySpectrumWeights = 0.6269199 * So/ActiveSource/BeamEnergySpectrumWeightsUnscaled

##### SOURCE GEOMETRY #####

s:Ge/ActiveSourceGeometry/Type              = "Group"
s:Ge/ActiveSourceGeometry/Parent            = "World"
b:Ge/ActiveSourceGeometry/IsParallel        = "True"
s:Ge/ActiveSourceGeometry/ParallelWorldName = "ActiveSourceWorld"

========================
#125I sphere layer center
#========================

s:Ge/SphereLayerCenter/Type         = "TsSphere"
s:Ge/SphereLayerCenter/Parent       = "ActiveSourceGeometry"
s:Ge/SphereLayerCenter/Material     = "MatRadioactiveLayer"
d:Ge/SphereLayerCenter/RMin         = 0.249 mm
d:Ge/SphereLayerCenter/RMax         = 0.25 mm
d:Ge/SphereLayerCenter/SPhi         = 0. deg
d:Ge/SphereLayerCenter/DPhi         = 360. deg
d:Ge/SphereLayerCenter/STheta       = 0. deg
d:Ge/SphereLayerCenter/DTheta       = 180. deg
d:Ge/SphereLayerCenter/TransX       = 0. mm
d:Ge/SphereLayerCenter/TransY       = 0. mm
d:Ge/SphereLayerCenter/TransZ       = 0. mm
d:Ge/SphereLayerCenter/RotX         = 0. deg
d:Ge/SphereLayerCenter/RotY         = 0. deg
d:Ge/SphereLayerCenter/RotZ         = 0. deg
s:Ge/SphereLayerCenter/DrawingStyle = "FullWireFrame"
s:Ge/SphereLayerCenter/Color        = "yellow"
b:Ge/SphereLayerCenter/IsParallel   = "T"
s:Ge/SphereLayerCenter/ParallelWorldName = "ActiveSourceWorld"

#========================
#125I sphere layer mid right
#========================

s:Ge/SphereLayerMidRight/Type         = "TsSphere"
s:Ge/SphereLayerMidRight/Parent       = "ActiveSourceGeometry"
s:Ge/SphereLayerMidRight/Material     = "MatRadioactiveLayer"
d:Ge/SphereLayerMidRight/RMin         = 0.249 mm
d:Ge/SphereLayerMidRight/RMax         = 0.25 mm
d:Ge/SphereLayerMidRight/SPhi         = 0. deg
d:Ge/SphereLayerMidRight/DPhi         = 360. deg
d:Ge/SphereLayerMidRight/STheta       = 0. deg
d:Ge/SphereLayerMidRight/DTheta       = 180. deg
d:Ge/SphereLayerMidRight/TransX       = 0. mm
d:Ge/SphereLayerMidRight/TransY       = 0. mm
d:Ge/SphereLayerMidRight/TransZ       = -0.6 mm
d:Ge/SphereLayerMidRight/RotX         = 0. deg
d:Ge/SphereLayerMidRight/RotY         = 0. deg
d:Ge/SphereLayerMidRight/RotZ         = 0. deg
s:Ge/SphereLayerMidRight/DrawingStyle = "FullWireFrame"
s:Ge/SphereLayerMidRight/Color        = "yellow"
b:Ge/SphereLayerMidRight/IsParallel   = "True"
s:Ge/SphereLayerMidRight/ParallelWorldName = "ActiveSourceWorld"

#========================
#125I sphere layer far right
#========================

s:Ge/SphereLayerFarRight/Type         = "TsSphere"
s:Ge/SphereLayerFarRight/Parent       = "ActiveSourceGeometry"
s:Ge/SphereLayerFarRight/Material     = "MatRadioactiveLayer"
d:Ge/SphereLayerFarRight/RMin         = 0.249 mm
d:Ge/SphereLayerFarRight/RMax         = 0.25 mm
d:Ge/SphereLayerFarRight/SPhi         = 0. deg
d:Ge/SphereLayerFarRight/DPhi         = 360. deg
d:Ge/SphereLayerFarRight/STheta       = 0. deg
d:Ge/SphereLayerFarRight/DTheta       = 180. deg
d:Ge/SphereLayerFarRight/TransX       = 0. mm
d:Ge/SphereLayerFarRight/TransY       = 0. mm
d:Ge/SphereLayerFarRight/TransZ       = -1.2 mm
d:Ge/SphereLayerFarRight/RotX         = 0. deg
d:Ge/SphereLayerFarRight/RotY         = 0. deg
d:Ge/SphereLayerFarRight/RotZ         = 0. deg
s:Ge/SphereLayerFarRight/DrawingStyle = "FullWireFrame"
s:Ge/SphereLayerFarRight/Color        = "yellow"
b:Ge/SphereLayerFarRight/IsParallel   = "True"
s:Ge/SphereLayerFarRight/ParallelWorldName = "ActiveSourceWorld"

#========================
#125I sphere layer mid left
#========================

s:Ge/SphereLayerMidLeft/Type         = "TsSphere"
s:Ge/SphereLayerMidLeft/Parent       = "ActiveSourceGeometry"
s:Ge/SphereLayerMidLeft/Material     = "MatRadioactiveLayer"
d:Ge/SphereLayerMidLeft/RMin         = 0.249 mm
d:Ge/SphereLayerMidLeft/RMax         = 0.25 mm
d:Ge/SphereLayerMidLeft/SPhi         = 0. deg
d:Ge/SphereLayerMidLeft/DPhi         = 360. deg
d:Ge/SphereLayerMidLeft/STheta       = 0. deg
d:Ge/SphereLayerMidLeft/DTheta       = 180. deg
d:Ge/SphereLayerMidLeft/TransX       = 0. mm
d:Ge/SphereLayerMidLeft/TransY       = 0. mm
d:Ge/SphereLayerMidLeft/TransZ       = 0.6 mm
d:Ge/SphereLayerMidLeft/RotX         = 0. deg
d:Ge/SphereLayerMidLeft/RotY         = 0. deg
d:Ge/SphereLayerMidLeft/RotZ         = 0. deg
s:Ge/SphereLayerMidLeft/DrawingStyle = "FullWireFrame"
s:Ge/SphereLayerMidLeft/Color        = "yellow"
b:Ge/SphereLayerMidLeft/IsParallel   = "True"
s:Ge/SphereLayerMidLeft/ParallelWorldName = "ActiveSourceWorld"

#========================
#125I sphere layer far left
#========================

s:Ge/SphereLayerFarLeft/Type         = "TsSphere"
s:Ge/SphereLayerFarLeft/Parent       = "ActiveSourceGeometry"
s:Ge/SphereLayerFarLeft/Material     = "MatRadioactiveLayer"
d:Ge/SphereLayerFarLeft/RMin         = 0.249 mm
d:Ge/SphereLayerFarLeft/RMax         = 0.25 mm
d:Ge/SphereLayerFarLeft/SPhi         = 0. deg
d:Ge/SphereLayerFarLeft/DPhi         = 360. deg
d:Ge/SphereLayerFarLeft/STheta       = 0. deg
d:Ge/SphereLayerFarLeft/DTheta       = 180. deg
d:Ge/SphereLayerFarLeft/TransX       = 0. mm
d:Ge/SphereLayerFarLeft/TransY       = 0. mm
d:Ge/SphereLayerFarLeft/TransZ       = 1.2 mm
d:Ge/SphereLayerFarLeft/RotX         = 0. deg
d:Ge/SphereLayerFarLeft/RotY         = 0. deg
d:Ge/SphereLayerFarLeft/RotZ         = 0. deg
s:Ge/SphereLayerFarLeft/DrawingStyle = "FullWireFrame"
s:Ge/SphereLayerFarLeft/Color        = "yellow"
b:Ge/SphereLayerFarLeft/IsParallel   = "True"
s:Ge/SphereLayerFarLeft/ParallelWorldName = "ActiveSourceWorld"

#==================
#Silver bead center
#==================

s:Ge/SBC/Type         = "TsSphere"
s:Ge/SBC/Parent       = "ActiveSourceGeometry"
s:Ge/SBC/Material     = "MatSilverRod"
d:Ge/SBC/RMin         = 0. mm
d:Ge/SBC/RMax         = 0.249 mm
d:Ge/SBC/SPhi         = 0. deg
d:Ge/SBC/DPhi         = 360. deg
d:Ge/SBC/STheta       = 0. deg
d:Ge/SBC/DTheta       = 180. deg
d:Ge/SBC/TransX       = 0. mm
d:Ge/SBC/TransY       = 0. mm
d:Ge/SBC/TransZ       = 0. mm
d:Ge/SBC/RotX         = 0. deg
d:Ge/SBC/RotY         = 0. deg
d:Ge/SBC/RotZ         = 0. deg
s:Ge/SBC/DrawingStyle = "FullWireFrame"
s:Ge/SBC/Color        = "red"
b:Ge/SBC/IsParallel   = "True"
s:Ge/SBC/ParallelWorldName = "ActiveSourceWorld"

#=====================
#Silver bead mid right
#=====================

s:Ge/SBmR/Type         = "TsSphere"
s:Ge/SBmR/Parent       = "ActiveSourceGeometry"
s:Ge/SBmR/Material     = "MatSilverRod"
d:Ge/SBmR/RMin         = 0. mm
d:Ge/SBmR/RMax         = 0.249 mm
d:Ge/SBmR/SPhi         = 0. deg
d:Ge/SBmR/DPhi         = 360. deg
d:Ge/SBmR/STheta       = 0. deg
d:Ge/SBmR/DTheta       = 180. deg
d:Ge/SBmR/TransX       = 0. mm
d:Ge/SBmR/TransY       = 0. mm
d:Ge/SBmR/TransZ       = -0.6 mm
d:Ge/SBmR/RotX         = 0. deg
d:Ge/SBmR/RotY         = 0. deg
d:Ge/SBmR/RotZ         = 0. deg
s:Ge/SBmR/DrawingStyle = "FullWireFrame"
s:Ge/SBmR/Color        = "red"
b:Ge/SBmR/IsParallel   = "True"
s:Ge/SBmR/ParallelWorldName = "ActiveSourceWorld"

#=====================
#Silver bead far right
#=====================

s:Ge/SBfR/Type         = "TsSphere"
s:Ge/SBfR/Parent       = "ActiveSourceGeometry"
s:Ge/SBfR/Material     = "MatSilverRod"
d:Ge/SBfR/RMin         = 0. mm
d:Ge/SBfR/RMax         = 0.249 mm
d:Ge/SBfR/SPhi         = 0. deg
d:Ge/SBfR/DPhi         = 360. deg
d:Ge/SBfR/STheta       = 0. deg
d:Ge/SBfR/DTheta       = 180. deg
d:Ge/SBfR/TransX       = 0. mm
d:Ge/SBfR/TransY       = 0. mm
d:Ge/SBfR/TransZ       = -1.2 mm
d:Ge/SBfR/RotX         = 0. deg
d:Ge/SBfR/RotY         = 0. deg
d:Ge/SBfR/RotZ         = 0. deg
s:Ge/SBfR/DrawingStyle = "FullWireFrame"
s:Ge/SBfR/Color        = "red"
b:Ge/SBfR/IsParallel   = "True"
s:Ge/SBfR/ParallelWorldName = "ActiveSourceWorld"

#====================
#Silver bead mid left
#====================

s:Ge/SBmL/Type         = "TsSphere"
s:Ge/SBmL/Parent       = "ActiveSourceGeometry"
s:Ge/SBmL/Material     = "MatSilverRod"
d:Ge/SBmL/RMin         = 0. mm
d:Ge/SBmL/RMax         = 0.249 mm
d:Ge/SBmL/SPhi         = 0. deg
d:Ge/SBmL/DPhi         = 360. deg
d:Ge/SBmL/STheta       = 0. deg
d:Ge/SBmL/DTheta       = 180. deg
d:Ge/SBmL/TransX       = 0. mm
d:Ge/SBmL/TransY       = 0. mm
d:Ge/SBmL/TransZ       = 0.6 mm
d:Ge/SBmL/RotX         = 0. deg
d:Ge/SBmL/RotY         = 0. deg
d:Ge/SBmL/RotZ         = 0. deg
s:Ge/SBmL/DrawingStyle = "FullWireFrame"
s:Ge/SBmL/Color        = "red"
b:Ge/SBmL/IsParallel   = "True"
s:Ge/SBmL/ParallelWorldName = "ActiveSourceWorld"

#====================
#Silver bead far left
#====================

s:Ge/SBfL/Type         = "TsSphere"
s:Ge/SBfL/Parent       = "ActiveSourceGeometry"
s:Ge/SBfL/Material     = "MatSilverRod"
d:Ge/SBfL/RMin         = 0. mm
d:Ge/SBfL/RMax         = 0.249 mm
d:Ge/SBfL/SPhi         = 0. deg
d:Ge/SBfL/DPhi         = 360. deg
d:Ge/SBfL/STheta       = 0. deg
d:Ge/SBfL/DTheta       = 180. deg
d:Ge/SBfL/TransX       = 0. mm
d:Ge/SBfL/TransY       = 0. mm
d:Ge/SBfL/TransZ       = 1.2 mm
d:Ge/SBfL/RotX         = 0. deg
d:Ge/SBfL/RotY         = 0. deg
d:Ge/SBfL/RotZ         = 0. deg
s:Ge/SBfL/DrawingStyle = "FullWireFrame"
s:Ge/SBfL/Color        = "red"
b:Ge/SBfL/IsParallel   = "True"
s:Ge/SBfL/ParallelWorldName = "ActiveSourceWorld"

##### FILLING GEOMETRY #####

s:Ge/FillingGeometry/Type              = "Group"
s:Ge/FillingGeometry/Parent            = "World"
b:Ge/FillingGeometry/IsParallel        = "True"
s:Ge/FillingGeometry/ParallelWorldName = "FillingWorld"

#====================
#Dry Air cylinder
#====================

s:Ge/DryAirCylinder/Type         = "TsCylinder"
s:Ge/DryAirCylinder/Parent       = "FillingGeometry"
s:Ge/DryAirCylinder/Material     = "G4_AIR"
d:Ge/DryAirCylinder/RMin         = 0. mm
d:Ge/DryAirCylinder/RMax         = 0.35 mm
d:Ge/DryAirCylinder/SPhi         = 0. deg
d:Ge/DryAirCylinder/DPhi         = 360. deg
d:Ge/DryAirCylinder/HL           = 1.6 mm
d:Ge/DryAirCylinder/TransX       = 0. mm
d:Ge/DryAirCylinder/TransY       = 0. mm
d:Ge/DryAirCylinder/TransZ       = 0. mm
d:Ge/DryAirCylinder/RotX         = 0. deg
d:Ge/DryAirCylinder/RotY         = 0. deg
d:Ge/DryAirCylinder/RotZ         = 0. deg
s:Ge/DryAirCylinder/DrawingStyle = "FullWireFrame"
s:Ge/DryAirCylinder/Color        = "white"
b:Ge/DryAirCylinder/IsParallel   = "True"
s:Ge/DryAirCylinder/ParallelWorldName = "FillingWorld"

#===============================
#Left cap of air 
#===============================

s:Ge/AirLeftCap/Type         = "TsSphere"
s:Ge/AirLeftCap/Parent       = "FillingGeometry"
s:Ge/AirLeftCap/Material     = "G4_AIR"
d:Ge/AirLeftCap/RMin         = 0. mm
d:Ge/AirLeftCap/RMax         = 0.35 mm
d:Ge/AirLeftCap/SPhi         = 0. deg
d:Ge/AirLeftCap/DPhi         = 360. deg
d:Ge/AirLeftCap/STheta       = 0. deg
d:Ge/AirLeftCap/DTheta       = 90. deg
d:Ge/AirLeftCap/TransX       = 0. mm
d:Ge/AirLeftCap/TransY       = 0. mm
d:Ge/AirLeftCap/TransZ       = 1.6 mm
d:Ge/AirLeftCap/RotX         = 0. deg
d:Ge/AirLeftCap/RotY         = 0. deg
d:Ge/AirLeftCap/RotZ         = 0. deg
s:Ge/AirLeftCap/DrawingStyle = "FullWireFrame"
s:Ge/AirLeftCap/Color        = "white"
b:Ge/AirLeftCap/IsParallel   = "True"
s:Ge/AirLeftCap/ParallelWorldName = "FillingWorld"

#================================
#Right cap of air
#================================

s:Ge/AirRightCap/Type         = "TsSphere"
s:Ge/AirRightCap/Parent       = "FillingGeometry"
s:Ge/AirRightCap/Material     = "G4_AIR"
d:Ge/AirRightCap/RMin         = 0. mm
d:Ge/AirRightCap/RMax         = 0.35 mm
d:Ge/AirRightCap/SPhi         = 0. deg
d:Ge/AirRightCap/DPhi         = 360. deg
d:Ge/AirRightCap/STheta       = 0. deg
d:Ge/AirRightCap/DTheta       = 90. deg
d:Ge/AirRightCap/TransX       = 0. mm
d:Ge/AirRightCap/TransY       = 0. mm
d:Ge/AirRightCap/TransZ       = -1.6 mm
d:Ge/AirRightCap/RotX         = 180. deg
d:Ge/AirRightCap/RotY         = 0. deg
d:Ge/AirRightCap/RotZ         = 0. deg
s:Ge/AirRightCap/DrawingStyle = "FullWireFrame"
s:Ge/AirRightCap/Color        = "white"
b:Ge/AirRightCap/IsParallel   = "True"
s:Ge/AirRightCap/ParallelWorldName = "FillingWorld"

##### SHIELDING #####

s:Ge/ShieldingGeo/Type              = "Group"
s:Ge/ShieldingGeo/Parent            = "World"
b:Ge/ShieldingGeo/IsParallel        = "True"
s:Ge/ShieldingGeo/ParallelWorldName = "ShieldingWorld"

#=====================
#Cylinder of Titanium
#=====================

s:Ge/TitaniumCylinder/Type         = "TsCylinder"
s:Ge/TitaniumCylinder/Parent       = "ShieldingGeo"
s:Ge/TitaniumCylinder/Material     = "Titanium"
d:Ge/TitaniumCylinder/RMin         = 0. mm
d:Ge/TitaniumCylinder/RMax         = 0.40 mm
d:Ge/TitaniumCylinder/SPhi         = 180. deg
d:Ge/TitaniumCylinder/DPhi         = 360. deg
d:Ge/TitaniumCylinder/HL           = 1.85 mm
d:Ge/TitaniumCylinder/TransX       = 0. mm
d:Ge/TitaniumCylinder/TransY       = 0. mm
d:Ge/TitaniumCylinder/TransZ       = 0. mm
d:Ge/TitaniumCylinder/RotX         = 0. deg
d:Ge/TitaniumCylinder/RotY         = 0. deg
d:Ge/TitaniumCylinder/RotZ         = 0. deg
s:Ge/TitaniumCylinder/Color        = "blue"
s:Ge/TitaniumCylinder/DrawingStyle = "FullWireFrame"
b:Ge/TitaniumCylinder/IsParallel   = "True"
s:Ge/TitaniumCylinder/ParallelWorldName = "ShieldingWorld"

#===============================
#Left half-sphere of Titanium 
#===============================

s:Ge/TitaniumLeftCap/Type         = "TsSphere"
s:Ge/TitaniumLeftCap/Parent       = "ShieldingGeo"
s:Ge/TitaniumLeftCap/Material     = "Titanium"
d:Ge/TitaniumLeftCap/RMin         = 0 mm
d:Ge/TitaniumLeftCap/RMax         = 0.4 mm
d:Ge/TitaniumLeftCap/SPhi         = 0. deg
d:Ge/TitaniumLeftCap/DPhi         = 360. deg
d:Ge/TitaniumLeftCap/STheta       = 0. deg
d:Ge/TitaniumLeftCap/DTheta       = 90. deg
d:Ge/TitaniumLeftCap/TransX       = 0. mm
d:Ge/TitaniumLeftCap/TransY       = 0. mm
d:Ge/TitaniumLeftCap/TransZ       = 1.85 mm
d:Ge/TitaniumLeftCap/RotX         = 0. deg
d:Ge/TitaniumLeftCap/RotY         = 0. deg
d:Ge/TitaniumLeftCap/RotZ         = 0. deg
s:Ge/TitaniumLeftCap/DrawingStyle = "FullWireFrame"
s:Ge/TitaniumLeftCap/Color        = "blue"
b:Ge/TitaniumLeftCap/IsParallel   = "True"
s:Ge/TitaniumLeftCap/ParallelWorldName = "ShieldingWorld"

#================================
#Right half-sphere of titanium
#================================

s:Ge/TitaniumRightCap/Type         = "TsSphere"
s:Ge/TitaniumRightCap/Parent       = "ShieldingGeo"
s:Ge/TitaniumRightCap/Material     = "Titanium"
d:Ge/TitaniumRightCap/RMin         = 0 mm
d:Ge/TitaniumRightCap/RMax         = 0.4 mm
d:Ge/TitaniumRightCap/SPhi         = 0. deg
d:Ge/TitaniumRightCap/DPhi         = 360. deg
d:Ge/TitaniumRightCap/STheta       = 0. deg
d:Ge/TitaniumRightCap/DTheta       = 90. deg
d:Ge/TitaniumRightCap/TransX       = 0. mm
d:Ge/TitaniumRightCap/TransY       = 0. mm
d:Ge/TitaniumRightCap/TransZ       = -1.85 mm
d:Ge/TitaniumRightCap/RotX         = 180. deg
d:Ge/TitaniumRightCap/RotY         = 0. deg
d:Ge/TitaniumRightCap/RotZ         = 0. deg
s:Ge/TitaniumRightCap/DrawingStyle = "FullWireFrame"
s:Ge/TitaniumRightCap/Color        = "blue"
b:Ge/TitaniumRightCap/IsParallel   = "True"
s:Ge/TitaniumRightCap/ParallelWorldName = "ShieldingWorld"

sv:Ph/Default/LayeredMassGeometryWorlds =  3 "ShieldingWorld" "FillingWorld" "ActiveSourceWorld"
