#
#
##########################################################################
#
# neutronSource-cf252
#
# This macro file sets G4GPS parameters to produce an isotropic 
# point source of neutrons with a Cf252 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 Cf-252 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 	6.6448E-07	0.1561
/gps/hist/point 	3.9087E-06	0.3638
/gps/hist/point 	2.4853E-05	0.9366
/gps/hist/point 	7.2135E-05	1.6283
/gps/hist/point 	1.4427E-04	2.3068
/gps/hist/point 	2.8854E-04	3.2596
/gps/hist/point 	5.3608E-04	4.4543
/gps/hist/point 	8.4110E-04	5.5851
/gps/hist/point 	1.8205E-03	8.1268
/gps/hist/point 	4.3281E-03	12.5860
/gps/hist/point 	7.8305E-03	16.9617
/gps/hist/point 	1.4427E-02	22.8319
/gps/hist/point 	2.8854E-02	32.0059
/gps/hist/point 	4.9326E-02	41.2979
/gps/hist/point 	6.9521E-02	48.5251
/gps/hist/point 	8.9628E-02	54.2773
/gps/hist/point 	1.2332E-01	61.9469
/gps/hist/point 	1.7380E-01	71.3864
/gps/hist/point 	2.2407E-01	78.4661
/gps/hist/point 	2.7424E-01	83.7758
/gps/hist/point 	3.2436E-01	87.9056
/gps/hist/point 	3.7444E-01	91.4454
/gps/hist/point 	4.2451E-01	94.3953
/gps/hist/point 	4.7456E-01	96.1652
/gps/hist/point 	5.2460E-01	97.9351
/gps/hist/point 	5.7464E-01	99.1150
/gps/hist/point 	6.4872E-01	99.7050
/gps/hist/point 	7.4889E-01	100.0000
/gps/hist/point 	8.4902E-01	99.4100
/gps/hist/point 	9.4912E-01	98.2301
/gps/hist/point 	1.0970E+00	95.2802
/gps/hist/point 	1.2974E+00	90.2655
/gps/hist/point 	1.4978E+00	84.5133
/gps/hist/point 	1.6980E+00	78.3186
/gps/hist/point 	1.8982E+00	71.9764
/gps/hist/point 	2.1465E+00	64.4051
/gps/hist/point 	2.4469E+00	55.7522
/gps/hist/point 	2.7952E+00	46.6814
/gps/hist/point 	3.2436E+00	36.6372
/gps/hist/point 	3.7444E+00	27.6106
/gps/hist/point 	4.2451E+00	20.5900
/gps/hist/point 	4.7456E+00	15.2212
/gps/hist/point 	5.4848E+00	9.7345
/gps/hist/point 	6.4872E+00	5.1327
/gps/hist/point 	7.4889E+00	2.6578
/gps/hist/point 	8.4902E+00	1.3599
/gps/hist/point 	9.4912E+00	0.6873
/gps/hist/point 	1.0492E+01	0.3451
/gps/hist/point 	1.1493E+01	0.1720
/gps/hist/point 	1.2493E+01	0.0850
/gps/hist/point 	1.3494E+01	0.0419
#
# Linear interpolation between the data points is required to obtain
# the smooth energy spectrum
/gps/hist/inter Lin

/Solid/beamOn 100000