#
#
##########################################################################
#
# neutronSource-AmBe
#
# This macro file sets G4GPS parameters to produce an isotropic 
# point source of neutrons with a AmBe energy spectrum.
# I. Pinera - May, 2015
##########################################################################
#
#
#
/IBD/on false
#/Solid/setConstructSM1Shielding false
/Solid/setConstructSM1FrontShielding false
#/Solid/setConstructSM1MuVeto false
#
/gps/particle neutron
/gps/pos/type Point
#for Pos1 (at 50 cm from the front HDPE neutron reflector)
/gps/pos/centre 0. 0. -54.725 cm
#for Pos2 (at 50 cm from the back HDPE neutron reflector)
#/gps/pos/centre 0. 0. 98.775 cm
/gps/ang/type iso
#/gps/ang/maxtheta 65 deg
#/gps/ang/rot2 0 -1 0

/gps/ene/type Arb
/gps/hist/type arb
#
# The AmBe neutron energy spectrum was taken from
# ISO 8529-1:2001
# Reference neutron radiations -- Part 1:
# Characteristics and methods of production
#
#                    Energy (MeV)    Rel. Yield
#
/gps/hist/point 	8.8069E-03	74.5014
/gps/hist/point 	2.0025E-01	86.4006
/gps/hist/point 	4.2642E-01	84.8238
/gps/hist/point 	6.3926E-01	76.1518
/gps/hist/point 	8.5529E-01	64.6711
/gps/hist/point 	1.0716E+00	57.9946
/gps/hist/point 	1.2869E+00	51.2195
/gps/hist/point 	1.5026E+00	47.4255
/gps/hist/point 	1.7129E+00	52.0325
/gps/hist/point 	1.9279E+00	57.6866
/gps/hist/point 	2.1433E+00	58.2656
/gps/hist/point 	2.3583E+00	58.2040
/gps/hist/point 	2.5736E+00	61.7886
/gps/hist/point 	2.7886E+00	76.3119
/gps/hist/point 	3.0038E+00	96.4770
/gps/hist/point 	3.2139E+00	100.0000
/gps/hist/point 	3.4288E+00	89.5048
/gps/hist/point 	3.6440E+00	83.1978
/gps/hist/point 	3.8590E+00	77.6053
/gps/hist/point 	4.0741E+00	72.8997
/gps/hist/point 	4.2841E+00	77.5068
/gps/hist/point 	4.4991E+00	82.2616
/gps/hist/point 	4.7142E+00	83.1978
/gps/hist/point 	4.9292E+00	86.1419
/gps/hist/point 	5.1443E+00	82.3848
/gps/hist/point 	5.3592E+00	70.8795
/gps/hist/point 	5.5743E+00	63.1436
/gps/hist/point 	5.7844E+00	55.8266
/gps/hist/point 	5.9993E+00	47.0806
/gps/hist/point 	6.2144E+00	47.9675
/gps/hist/point 	6.4294E+00	52.7716
/gps/hist/point 	6.6444E+00	49.5935
/gps/hist/point 	6.8545E+00	44.1734
/gps/hist/point 	7.0694E+00	43.4590
/gps/hist/point 	7.2845E+00	45.5285
/gps/hist/point 	7.4995E+00	48.6327
/gps/hist/point 	7.7145E+00	49.8645
/gps/hist/point 	7.9245E+00	45.7995
/gps/hist/point 	8.1395E+00	37.2506
/gps/hist/point 	8.3546E+00	26.2331
/gps/hist/point 	8.5695E+00	16.8662
/gps/hist/point 	8.7846E+00	11.5447
/gps/hist/point 	8.9996E+00	9.4937
/gps/hist/point 	9.2146E+00	10.3252
/gps/hist/point 	9.4246E+00	13.7127
/gps/hist/point 	9.6396E+00	16.1678
/gps/hist/point 	9.8546E+00	14.9593
/gps/hist/point 	1.0080E+01	11.0976
/gps/hist/point 	1.0300E+01	10.5285
/gps/hist/point 	1.0500E+01	7.9106
/gps/hist/point 	1.0700E+01	4.2967
#
# Linear interpolation between the data points is required to obtain
# the smooth energy spectrum
/gps/hist/inter Lin

/Solid/beamOn 100000