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


:start geometry definition:

    :start geometry:
        library = egs_cones
        type = EGS_ConeStack
        name = cyl_phantom
        axis = 0 0 0.0 0 0 1
        :start layer:
           thickness = 10
           top radii = 20
           bottom radii = 20
           media = AIR521ICRU
        :stop layer:
        :start layer:
           thickness = 10
           top radii = 20
           bottom radii = 20
           media = AIR521ICRU
        :stop layer:
    :stop geometry:

    simulation geometry = cyl_phantom
: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_dynamic_shape
           :start shape:
                library = egs_rectangle
                rectangle = -1 -1 1 1
           :stop shape:
           :start motion:
                control point 1 = 0.0 0 0 0 0 0 0
                control point 2 = 1.0 0 0 0 0 0 10
                control point 3 = 0.5 0 0 0 0 0 20
                control point 4 = 0.75 0 0 0 0 0 30
                control point 5 = 1.0 0 0 0 0 0 40
           :stop motion:
        :stop shape:
        charge = 0
        direction = 0 0 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= .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= On
 Compton cross sections= default
 Pair angular sampling= Simple
 Pair cross sections= BH
 Photoelectron angular sampling= On
 Rayleigh scattering= Off
 Atomic relaxations= On
 Electron impact ionization= On
 Photon cross sections= xcom
 Photon cross-sections output= Off

:stop MC transport parameter:

:start ausgab object definition:

:start ausgab object:
      library = egs_track_scoring
      name = sample_tracks
      score electrons = yes
      score positrons = no
      include time index = yes
      stop scoring = 10000000
      buffer size = 100024
:stop ausgab object:

:start ausgab object:
      library = egs_phsp_scoring
      name = dyn_shape_test
      output format = EGSnrc
      particle type = all
      from regions = 0
      to regions = 1
:stop ausgab object:

:stop ausgab object definition:

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

######################################### Scoring options.
#        The entire section that follows can be missing. If missing, just
#        energy deposition fractions will be calculated. If present,
#        the user can specify in which regions to score a pulse height
#        distribution.
#