################################################################################
#
# Model of dose in spherical cavity for benchmarking bremsstrahlung
#
################################################################################


################################################################################
### Run
################################################################################
:start run control:
    ncase                       = 1e9
    geometry error limit        = 1e2
    statistical accuracy sought = 1e-5
    nbatch                      = 10
    nchunk                      = 10
:stop run control:

:start rng definition:
    initial seeds = 4523 22934
:stop rng definition:


################################################################################
### Geometry
################################################################################
:start geometry definition:

    ############################################################################
    ### world
    ############################################################################
    :start geometry:
        name     = box
        library  = egs_ndgeometry
        type     = EGS_XYZGeometry
        x-planes = -1.51 1.51
        y-planes = -1.51 1.51
        z-planes = -1.51 1.51
        :start media input:
            media = air
        :stop media input:
    :stop geometry:

    ############################################################################
    ### pmma sphere with air cavity
    ############################################################################
    :start geometry:
        name    = sphere
        library = egs_spheres
        radii   = 0.5 1.5
        :start media input:
            media = air pmma
            set medium = 1 1
        :stop media input:
    :stop geometry:

    ############################################################################
    ### inscribe sphere in air_box
    ############################################################################
    :start geometry:
        name    = world
        library = egs_genvelope
        base geometry = box
        inscribed geometries = sphere
    :stop geometry:

    ############################################################################
    ### simulation geometry
    ############################################################################
    simulation geometry = world

:stop geometry definition:


################################################################################
### Media
################################################################################
:start media definition:

    ### energy transport cutoffs
    ae =  0.521
    ap =  0.010
    ue = 30.511
    up = 30.000

    :start air:
        density correction file = air_dry_nearsealevel
    :stop air:

    :start pmma:
        density correction file = polymethylmethacrylate__lucite___perspex___plexiglas_
    :stop pmma:

:stop media definition:


################################################################################
### Source
################################################################################
:start source definition:

    ### parallel photon beam
    :start source:
        name    = parallel_beam
        library = egs_parallel_beam
        charge  = 0
        :start shape:
            library = egs_circle
            radius = 1.5
            :start transformation:
                translation = 0 0 -10
            :stop transformation:
        :stop shape:
        :start spectrum:
            type = monoenergetic
            energy = 23
        :stop spectrum:
    :stop source:

    ############################################################################
    ### use this source for the simulation
    ############################################################################
    simulation source = parallel_beam

:stop source definition:


################################################################################
### MC Transport Parameters
################################################################################
:start mc transport parameter:

    ### transport
    Global Ecut = 0.521
    Global Pcut = 0.010

    ### cross sections
    Photon cross sections           = epdl      # xcom | epdl | si
    Photoelectron angular sampling  = On        # On  | Off
    Rayleigh scattering             = On        # On  | Off | custom
    Triplet production              = On        # On  | Off
    Photonuclear attenuation        = On        # On  | Off
    Radiative Compton corrections   = Off       # On  | Off ("On" requires rad_compton1.mortran)
    Pair cross sections             = NRC       # NRC | BH
    Pair angular sampling           = KM        # KM  | Off | Simple
    Brems cross sections            = NRC       # NRC | NIST | BH
    Brems angular sampling          = KM        # KM  | Simple
    Spin effects                    = On        # On  | Off

    ### atomic binding
    Atomic relaxations              = Off        # On  | Off
    Bound Compton scattering        = Off        # On  | Off | norej
    Electron Impact Ionization      = Off        # On  | Off

    ### electron transport algorithm
    Electron-step algorithm         = EGSnrc    # EGSnrc | PRESTA-I | PRESTA-II | default
    Boundary crossing algorithm     = Exact     # Exact  | PRESTA-I
    Skin depth for BCA              = 3
    Ximax                           = 0.5
    ESTEPE                          = 0.25
    Global Smax                     = 1e10

:stop mc transport parameter:

################################################################################
### Scoring
################################################################################
:start scoring options:

    :start calculation geometry:
        geometry name   = world
        cavity geometry = cavity
        cavity regions  = 1
        cavity mass     = 0.0006308     # volume = 0.523599 cm^3
        enhance regions = 2
        enhancement     = 1024

    :stop calculation geometry:

:stop scoring options:


################################################################################
### Variance reduction
################################################################################
:start variance reduction:
    cs enhancement = 1
:stop variance reduction:

# OFF: cs enhancement = 0:   dose = 9.6485e-12 +/- 0.09%
# ON:  cs enhancement = 1:   dose = 5.5757e-12 +/- 0.33%