def test_time_coincidence_window(self):
"""
test whether the timing window is correct (-1 when all timing is included)
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
#ht = fh.Histories('92235_18MeV_timing.cgmf',nevents=nread)
h = fh.Histories(
'../../utils/cgmf/tests/u235nf-th-events/histories.cgmf.parallel.0.reference',
nevents=nread)
ht = fh.Histories(
'../../utils/cgmf/tests/u235nf-time-events/histories.cgmf.parallel.0.reference',
nevents=nread)
self.assertTrue(h.getTimeCoincidenceWindow() > 0)
self.assertEqual(ht.getTimeCoincidenceWindow(), 1.5e-7)
def test_nubarg_for_events(self):
"""
test the nubargtot function calculation (average gammas per event)
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories('../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',nevents=nread)
nug = h.getNugtot()
self.assertEqual(np.mean(nug),h.nubargtot())
def test_gamma_ages(self):
"""
ensure the gamma times exist when they should be calculated
"""
nread = 10
#h = fh.Histories('92235_18MeV_timing.cgmf',nevents=nread)
h = fh.Histories(
'../../utils/cgmf/tests/u235nf-time-events/histories.cgmf.parallel.0.reference',
nevents=nread)
self.assertTrue(len(h.getGammaAges()) > 0)
def test_number_of_fragments(self):
"""
test that the number of fragments is twice the number of events
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories(
'../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',
nevents=nread)
self.assertEqual(h.numberFragments, 2 * nread)
def test_nevents_read(self):
"""
test that the correct number of events has been read
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories(
'../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',
nevents=nread)
self.assertEqual(h.numberEvents, nread)
def test_mean_list_calculation(self):
"""
test the calculation of the mean of a list
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories('../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',nevents=nread)
nE = h.getNeutronEcm()
l = []
for x in nE:
l+=x
self.assertEqual(np.mean(l),h.meanList(nE))
def test_nubar_for_fragments(self):
"""
test the nubar function calculation (average neutrons per fragment)
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories('../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',nevents=nread)
nuLF = h.getNuLF()
nuHF = h.getNuHF()
nutot = np.sum(nuLF)+np.sum(nuHF)
nutot = nutot/(2.*len(nuLF))
self.assertEqual(h.nubar(),nutot)
def test_nubar_for_events(self):
"""
test the nubartot function calculation (average neutrons per event)
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories('../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',nevents=nread)
nuLF = h.getNuLF()
nuHF = h.getNuHF()
nuPF = h.getPreFissionNu()
nutot = nuLF + nuHF + nuPF
nutot = np.mean(nutot)
self.assertEqual(h.nubartot(),nutot)
def test_neutron_multiplicity_distribution(self):
"""
test the P(nu) calculation
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories('../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',nevents=nread)
nu,pnu = h.Pnu()
nuvals = h.getNutot()
numax = np.max(nuvals)
self.assertEqual(nu[-1],numax)
for i in range(numax+1):
prob = (len(nuvals[nuvals==i]))/float(nread)
self.assertEqual(prob,pnu[i])
def test_gamma_multiplicity_distribution(self):
"""
test the gamma multiplicity distribution, P(N_gamma)
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories('../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',nevents=nread)
nug,pnug = h.Pnug()
nugvals = h.getNugtot()
nugmax = np.max(nugvals)
self.assertEqual(nug[-1],nugmax)
for i in range(nugmax+1):
prob = (len(nugvals[nugvals==i]))/float(nread)
self.assertEqual(prob,pnug[i])
def test_neutron_spectrum_energy_threshold(self):
"""
test the implementation of the neutron energy threshold
specificially in the PFNS, but is constructed the same elsewhere
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories('../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',nevents=nread)
Eth = 0.5
eng,pfns = h.pfns(Eth=Eth)
eng0,pfns0 = h.pfns()
emin = eng[pfns>0]
emin0 = eng0[pfns0>0]
if (Eth>0):
self.assertTrue(emin0[0]<emin[0])
else:
self.assertEqual(emin0[0],emin[0])
def test_mean_neutron_energy_in_lab_frame(self):
"""
test the calculation of the average neutron energy in the lab frame
"""
nread = 10
#h = fh.Histories('92235_18MeV.cgmf',nevents=nread)
h = fh.Histories('../../utils/cgmf/tests/u235nf-18MeV-events/histories.cgmf.parallel.0.reference',nevents=nread)
nE = h.getNeutronElab()
nPFE = h.getPreFissionNeutronElab()
nuPF = h.getPreFissionNu()
l = []
for x in nE:
l+=x
for i in range(len(nuPF)):
if (nuPF[i]>0):
l+=nPFE[i]
self.assertEqual(np.mean(l),h.meanNeutronElab())