def set_temporary_blank(self, k, v, e, f, verbose=False):
if verbose:
self.debug('temp blank {}({:0.4f}+/-{:0.4f}) fit={}'.format(
k, v, e, f))
if k in self.isotopes:
iso = self.isotopes[k]
tb = iso.temporary_blank
if tb is None:
iso.temporary_blank = tb = Blank(iso.name, iso.detector)
tb.value, tb.error, tb.fit = v, e, f
def set_temporary_blank(self, k, v, e, f):
tol = 0.00001
if k in self.isotopes:
iso = self.isotopes[k]
if iso.temporary_blank is not None:
tb = iso.temporary_blank
if abs(tb.value - v) < tol and abs(tb.error - e) < tol:
return
else:
self.debug('temp blank {}({:0.4f}+/-{:0.4f}) fit={}'.format(k, v, e, f))
tb.value, tb.error, tb.fit = v, e, f
else:
self.debug('temp blank {}({:0.4f}+/-{:0.4f}) fit={}'.format(k, v, e, f))
iso.temporary_blank = b = Blank(k, iso.detector)
b.value = v
b.error = e
b.fit = f
def _sync(self, obj):
for dbiso in obj.isotopes:
r = dbiso.results[-1]
uv = r.Iso
ee = r.IsoEr
bv = r.Bkgd
be = r.BkgdEr
key = dbiso.Label
iso = Isotope(name=key, value=uv, error=ee)
iso.baseline = Baseline(name=key,
reverse_unpack=True,
dbrecord=dbiso.baseline,
unpacker=lambda x: x.PeakTimeBlob,
fit='average_SEM')
iso.blank = Blank(name=key, value=bv, error=be)
self.isotopes[key] = iso
def _sync(self, obj):
arar = obj.araranalyses[-1]
if arar:
self.j = ufloat(arar.JVal, arar.JEr)
self.age = arar.Age
self.age_err = arar.ErrAge
self.age_err_wo_j = arar.ErrAgeWOErInJ
self.rad40_percent = ufloat(arar.PctRad, arar.PctRadEr)
for dbiso in obj.isotopes:
r = dbiso.results[-1]
uv = r.Iso
ee = r.IsoEr
bv = r.Bkgd
be = r.BkgdEr
key = dbiso.Label
n = dbiso.NumCnts
iso = Isotope(name=key, value=uv, error=ee, n=n)
det =dbiso.detector
iso.ic_factor=ufloat(det.ICFactor, det.ICFactorEr)
iso.fit = r.fit.Label.lower()
iso.baseline = Baseline(name=key,
reverse_unpack=True,
dbrecord=dbiso.baseline,
unpack=True,
unpacker=lambda x: x.PeakTimeBlob,
error_type='SEM',
fit='average')
iso.baseline.set_filter_outliers_dict()
iso.blank = Blank(name=key, value=bv, error=be)
self.isotopes[key] = iso
def set_temporary_blank(self, k, v, e):
self.debug('setting temporary blank iso={}, v={}, e={}'.format(
k, v, e))
if self.isotopes.has_key(k):
iso = self.isotopes[k]
iso.temporary_blank = Blank(value=v, error=e)
def set_temporary_blank(self, k, v, e):
if self.isotopes.has_key(k):
iso = self.isotopes[k]
iso.temporary_blank = Blank(value=v, error=e)