class IsotopeTestCase(unittest.TestCase):
def setUp(self):
self.iso = Isotope()
xs = linspace(10, 410, 400)
c = 1000.0
b = -2.0
a = 0.0
ys = a * xs * xs + b * xs + c
self.iso.trait_set(xs=xs, ys=ys)
self.iso.fit = 'parabolic'
def test_value(self):
v = self.iso.value
self.assertEqual(v, 999.999999999995)
def test_filtered_value(self):
self.iso.ys[[1, 2, 3]] = [0, 1, 1]
d = {
'filter_outliers': True,
'filter_outliers_iteration': 1,
'filter_outliers_std_dev': 2
}
self.iso.set_filtering(d)
v = self.iso.value
self.assertEqual(v, 999.99999999997499)
class IsotopeTestCase(unittest.TestCase):
def setUp(self):
self.iso = Isotope()
xs = linspace(10, 410, 400)
c = 1000.0
b = -2.0
a = 0.0
ys = a * xs * xs + b * xs + c
self.iso.trait_set(xs=xs, ys=ys)
self.iso.fit = 'parabolic'
def test_value(self):
v = self.iso.value
self.assertEqual(v, 999.999999999995)
def test_filtered_value(self):
self.iso.ys[[1, 2, 3]] = [0, 1, 1]
d = {'filter_outliers': True,
'filter_outliers_iteration': 1,
'filter_outliers_std_dev': 2}
self.iso.set_filtering(d)
v = self.iso.value
self.assertEqual(v, 999.99999999997499)
def set_isotope(self, iso, v, **kw):
if not self.isotopes.has_key(iso):
niso = Isotope(name=iso)
self.isotopes[iso] = niso
else:
niso = self.isotopes[iso]
niso.set_uvalue(v)
niso.trait_set(**kw)
return niso
def set_isotope(self, iso, v, **kw):
if not self.isotopes.has_key(iso):
niso = Isotope(name=iso)
self.isotopes[iso] = niso
else:
niso = self.isotopes[iso]
niso.set_uvalue(v)
niso.trait_set(**kw)
return niso
def _sync_isotopes(self, yd, meas_analysis, unpack):
"""
load isotopes from file
"""
isos = {}
# print yd['isotopes']
for iso in yd['isotopes']:
ii = Isotope(name=iso['name'],
detector=iso['detector'],
ic_factor=self._to_ufloat(iso, 'ic_factor'),
discrimination=self._to_ufloat(iso, 'discrimination'))
ii.trait_set(_value=iso['value'], _error=iso['error'])
ii.set_blank(iso['blank'], iso['blank_err'])
ii.set_baseline(iso['baseline'], iso['baseline_err'])
isos[iso['name']] = ii
self.isotopes = isos
def _sync_isotopes(self, yd, meas_analysis, unpack):
"""
load isotopes from file
"""
isos = {}
# print yd['isotopes']
for iso in yd['isotopes']:
ii = Isotope(name=iso['name'],
detector=iso['detector'],
ic_factor=self._to_ufloat(iso, 'ic_factor'),
discrimination=self._to_ufloat(iso, 'discrimination'))
ii.trait_set(_value=iso['value'], _error=iso['error'])
ii.set_blank(iso['blank'], iso['blank_err'])
ii.set_baseline(iso['baseline'], iso['baseline_err'])
isos[iso['name']] = ii
self.isotopes = isos
