# Script file: example.s
#
#    Example experiment script file for the ABC simulator.  First, a universe
# is defined as a box with 20m sides, and given a color 99 (white) and material 
# 1 (air).  A new particle, the ABCino, is created and given several possible
# decay modes.  The experiment consists of a cylindrical detector with a small
# magnetic field.  Surrounding the cylinder are vertical slabs of aluminum,
# and a thin slab of iron in front.  Ten generators are defined, one for each
# decay mode of the ABCino, along with decays for kaons, muons, and pions.
# Finally, protons are generated in the backwards direction to test their 
# effects in the iron.

# Define universe:
#   col mat xmin xmax ymin ymax zmin zmax
  u 99  01  -10  +10  -10  +10  -10  +10

# Define new particle:
#   num/name     col  chrg  mass(GeV)  life(m)  decay
  p 99 ABCino    01   +0.0  4.00       0.5      4
       070809 0.25                              # pi- pi+ pi0
       121314 0.25                              # K- K+ KL
       121315 0.25                              # K- K+ KS
       171819 0.25                              # S- S+ S0

# Define objects:
#   num/name   shp col mat det bty/mg  pos-(x,y,z)   siz-(1,2,3)   ang-(a,b,c)
  o 01 chamber cyl 01  01  1   d  .20  +0.0 +0.0  0  4.0 10. 10.     0  0  0
  o 11 block1  box 02  03  1   0  0.0  +3.5 +6.06 0  1.0 3.0 10.    +60 0  0
  o 12 block1  box 02  03  1   0  0.0  +6.06 +3.5 0  1.0 3.0 10.    +30 0  0
  o 13 block1  box 02  03  1   0  0.0  +7.00  0   0  1.0 3.0 10.     0  0  0
  o 14 block1  box 02  03  1   0  0.0  +6.06 -3.5 0  1.0 3.0 10.    -30 0  0
  o 15 block1  box 02  03  1   0  0.0  +3.5 -6.06 0  1.0 3.0 10.    -60 0  0
  o 20 plate   box 03  06  1   0  0.0  -7.0 +0.0  0  0.2 10. 10.     0  0  0

# Define generator:
#   num prt flg       pos-(x,y,z)        mom-(x,y,z)        delmom
  g 01  99  1      +0.0  +0.0  +0.0   0.50   0.00   0.00     0.00  # ABCino
  g 02  99  2      +0.0  +0.0  +0.0   0.50   0.00   0.00     0.00  #   dec2
  g 03  99  3      +0.0  +0.0  +0.0   0.50   0.00   0.00     0.00  #   dec3
  g 04  99  4      +0.0  +0.0  +0.0   0.50   0.00   0.00     0.00  #   dec4
  g 05  12  2      +0.0  +0.0  +0.0   1.00   0.00   0.00     0.00  # K- dec2
  g 06  13  3      +0.0  +0.0  +0.0   1.00   0.00   0.00     0.00  # K+ dec3
  g 07  14  6      +0.0  +0.0  +0.0   1.00   0.00   0.00     0.00  # KL dec6
  g 08  05  1      +0.0  +0.0  +0.0   0.00   0.00   0.00     0.00  # mu- decay
  g 09  07  1      +0.0  +0.0  +0.0   0.50   0.00   0.00     0.00  # pi- decay
  g 10  10  0      +0.0  +0.0  +0.0  -1.00   0.00   0.00     0.50  # proton