# #
# 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"