###############################################################################
#
#  EGSnrc egs++ tutor7pp application sample input file
#  Copyright (C) 2015 National Research Council Canada
#
#  This file is part of EGSnrc.
#
#  EGSnrc 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.
#
#  EGSnrc 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.
#
#  You should have received a copy of the GNU Affero General Public License
#  along with EGSnrc. If not, see <http://www.gnu.org/licenses/>.
#
###############################################################################
#
#  Author:          Iwan Kawrakow, 2005
#
#  Contributors:
#
###############################################################################
#
#  An example input file for tutor7pp.
#
###############################################################################


######################################## Define the geometry:
#                                        a 0.1mm plate of W
#
:start geometry definition:

    :start geometry:
        library = egs_ndgeometry
        type = EGS_XYZGeometry
        name = target
        x-planes = -100 100
        y-planes = -100 100
        z-planes = 0 0.01 10.0 10.001
         :start media input:
            media = VACUUM W521ICRU
            set medium = 0 1
        :stop media input:
    :stop geometry:

    :start geometry:
        library = egs_box
        name = surroundings
        box size = 40.0
        :start media input:
           media = VACUUM
        :stop media input:
    :stop geometry:

    :start geometry:
        library = egs_genvelope
        name = xtube
        base geometry = surroundings
        inscribed geometries = target
    :stop geometry:

    simulation geometry = target
:stop geometry definition:

######################################## Define the source:
#                                       shared library source
#
:start source definition:
    :start source:
        library = egs_parallel_beam
        name = e_beam
        :start shape:
           library = egs_circle
           radius = 0.01
        :stop shape:
        direction = 0 0 1
        charge = -1
        :start spectrum:
            type = monoenergetic
            energy = 0.45
        :stop spectrum:
    :stop source:

    simulation source = e_beam

:stop source definition:

######################################### Transport parameters
#                                         you can add any other option
#                                         supported by EGSnrc
#
:start MC transport parameter:
 Global ECUT= 0.521
 Global PCUT= 0.01
 Global SMAX= 5
 ESTEPE= 0.25
 XIMAX= 0.5
 Boundary crossing algorithm= exact
 Skin depth for BCA= 0
 Electron-step algorithm= PRESTA-II
 Spin effects= On
 Brems angular sampling= Simple
 Brems cross sections= NIST
 Bound Compton scattering= Off
 Compton cross sections= default
 Pair angular sampling= Simple
 Pair cross sections= BH
 Photoelectron angular sampling= On
 Rayleigh scattering= On
 Atomic relaxations= Off
 Electron impact ionization= Off
 Photon cross sections= xcom
 Photon cross-sections output= Off

:stop MC transport parameter:

:start ausgab object definition:
    :start ausgab object:
      library = egs_phsp_scoring
      name = 450kVp_tutor7pp_split_test_relax_off
      output format = EGSnrc
      particle type = all
      from regions = 1
      to regions = 2
    :stop ausgab object:

    :start ausgab object:
      library = egs_radiative_splitting
      name = brem_splitting
      splitting type = directional
      splitting = 1000
      field size = 10
      ssd = 10.0
    :stop ausgab object:

:stop ausgab object definition:

######################################## Run control
:start run control:
    ncase = 8000000
    #statistical accuracy sought = 1  # this input is optional
    max cpu hours allowed = 15       # this input is optional
    #calculation = combine   # could als be first or restart or analyze
                             # the input is optional and if missing, a
                             # fresh calculation will be done

:stop run control: