#
#
##########################################################################
# neutroSource-cosmic.mac
#
# This macro file simulates neutrons induced by
# cosmic rays. The spectrum is based on 
# Gordon et al. 2004  "Measurement of the flux and energy spectrum of 
# cosmic-ray induced neutrons on the ground" IEEE Transactions on 
#  Nuclear Science 51(6)
#
##########################################################################
#
#
#
/IBD/on false
#
/gps/particle neutron
#/gps/pos/type Point
#/gps/pos/centre -100. 0. 0. cm
/gps/pos/type Surface
/gps/pos/shape Sphere
/gps/pos/centre 0. 0. 0. cm
/gps/pos/radius 150. cm
/gps/ang/type iso
# /gps/ang/maxtheta 15 deg
# /gps/ang/rot2 0 0 1

/gps/ene/type Arb
/gps/hist/type arb
#
# The neutron energy spectrum was taken from
# arXiv:hep-ex/0101049v1
#
#                Energy (MeV)  Rel. Yield
#

/gps/hist/point 1.054   6.83E-04
/gps/hist/point 1.165   8.19E-04
/gps/hist/point 1.287   7.61E-04
/gps/hist/point 1.423   7.02E-04
/gps/hist/point 1.572   6.00E-04
/gps/hist/point 1.738   5.72E-04
/gps/hist/point 1.920   5.06E-04
/gps/hist/point 2.122   5.02E-04
/gps/hist/point 2.346   5.44E-04
/gps/hist/point 2.592   4.30E-04
/gps/hist/point 2.865   3.34E-04
/gps/hist/point 3.166   2.65E-04
/gps/hist/point 3.499   1.86E-04
/gps/hist/point 3.867   1.64E-04
/gps/hist/point 4.274   1.73E-04
/gps/hist/point 4.724   1.88E-04
/gps/hist/point 5.220   1.53E-04
/gps/hist/point 5.769   1.25E-04
/gps/hist/point 6.376   1.16E-04
/gps/hist/point 7.047   8.90E-05
/gps/hist/point 7.788   7.16E-05
/gps/hist/point 8.607   6.73E-05
/gps/hist/point 9.512   5.53E-05
/gps/hist/point 10.51   4.58E-05
/gps/hist/point 11.62   4.09E-05
/gps/hist/point 12.84   3.80E-05
/gps/hist/point 14.19   3.44E-05
/gps/hist/point 16.16   3.02E-05
/gps/hist/point 18.52   3.22E-05
/gps/hist/point 25.70   2.59E-05
/gps/hist/point 44.19   2.09E-05
/gps/hist/point 75.98   1.53E-05
/gps/hist/point 130.7   9.64E-06
/gps/hist/point 224.6   4.30E-06
/gps/hist/point 386.3   1.33E-06
/gps/hist/point 664.2   3.99E-07
/gps/hist/point 1142    1.02E-07
/gps/hist/point 1964    2.24E-08
/gps/hist/point 3376    3.36E-09
/gps/hist/point 5805    4.71E-10
/gps/hist/point 9982    9.87E-11
/gps/hist/point 17160   3.83E-11

#
# Linear interpolation between the data points is required to obtain
# the smooth energy spectrum
/gps/hist/inter Lin

/Solid/beamOn 1000