# seed name: OncoSeed_6711
# air kerma: 3.7666E-14

#  Copyright (C) 2020, Rowan M Thomson and D.W.O. Rogers
#
# This file is part of the egs_brachy package to be used with EGSnrc
#
#  egs_brachy is free software: you can redistribute it and/or modify it
#  under the terms of the GNU Affero General Public License as published
#  by the Free Software Foundation, either version 3 of the License, or
#  (at your option) any later version.
#
#  egs_brachy is distributed in the hope that it will be useful, but
#  WITHOUT ANY WARRANTY; without even the implied warranty of
#  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
#  Affero General Public License for more details:
#  <http://www.gnu.org/licenses/>.
#
# If you use this file, please cite
#       M Chamberland, R E P Taylor, D W O Rogers and R M Thomson
#       egs_brachy: a versatile and fast Monte Carlo code for brachytherapy
#       Phys Med Biol 61(2016) 8214 -- 8231
# and
#       Habib Safigholi, Marc J.P. Chamberland, Randle E.P. Taylor,
#       Christian H. Allen, Martin P. Martinov, D. W. O. Rogers and Rowan
#       M. Thomson,  Update of the CLRP TG-43 parameter database for
#       low-energy brachytherapy sources, Med. Phys. 47 (2020) 4656-4669
#

:start geometry definition:
    # the set of planes needed for the CD geometry
    :start geometry:
        library = egs_planes
        name = planes
        type = EGS_Zplanes
        positions = -0.2275 -0.1875 0.1875 0.2275
    :stop geometry:

    :start geometry:
        library = egs_spheres
        name = sph_end_cap_1
        midpoint = 0 0  -0.1875
        type = EGS_cSpheres
        radii = 0.04

        :start media input:
            media = Ti
            set medium = 0 0
        :stop media input:
    :stop geometry:

    :start geometry:
        library = egs_cylinders
        name = cladding
        midpoint = 0 0 0
        type = EGS_ZCylinders
        radii = 0.033 0.04

        :start media input:
            media = AIR_TG43 Ti
            set medium = 0 0
            set medium = 1 1
        :stop media input:
    :stop geometry:

    :start geometry:

        library = egs_cones
       type    = EGS_ConeStack
        axis    = 0 0 -0.140   0 0 1
        name    = source
 
      #   source end cylinder
        :start layer:
            thickness    = 0.000175
            top radii    = 0.0175
            bottom radii = 0.017675
            media        = AgBrAgI_6.20
        :stop layer:

       #  source end cone
        :start layer:
           thickness    = 0.007325
           top radii    = 0.0175 0.017675
           bottom radii = 0.024825 0.025
           media        = Ag AgBrAgI_6.20
        :stop layer:
 
      #  mid source
        :start layer:
           thickness    = 0.265
           bottom radii = 0.024825 0.025
           media        =  Ag AgBrAgI_6.20
        :stop layer:
 
     #   source end cone
        :start layer:
           thickness    = 0.007325
           bottom radii = 0.0175 0.017675
           media        = Ag AgBrAgI_6.20
        :stop layer:

       # source end cylinder
        :start layer:
           thickness    = 0.000175
           top radii    = 0.017675
           bottom radii = 0.0175
           media        = AgBrAgI_6.20
        :stop layer:

    :stop geometry:

    :start geometry:
        library = egs_genvelope
        name = seed_center
        base geometry = cladding
        inscribed geometries = source
    :stop geometry:

    :start geometry:
        library = egs_spheres
        name = sph_end_cap_2
        midpoint = 0 0  0.1875
        type = EGS_cSpheres
        radii = 0.04

        :start media input:
            media = Ti
            set medium = 0 0
        :stop media input:
    :stop geometry:

    :start geometry:
        library = egs_cdgeometry
        name = seed_unwrapped
        base geometry = planes
        set geometry = 0 sph_end_cap_1
        set geometry = 1 seed_center
        set geometry = 2 sph_end_cap_2
        boundary tolerance = 1e-6
    :stop geometry:
	
    :start geometry:
		library = egs_rz
		name = wrapper
		radii = 0.041
		z-planes = -0.2285 0.2285
		:start media input:
			media = WATER_0.998
		:stop media input:
    :stop geometry:

    :start geometry:
        library = egs_genvelope
        name = seed
        base geometry = wrapper
        inscribed geometries = seed_unwrapped
    :stop geometry:

    simulation geometry = seed
	
:stop geometry definition: