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 _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