def _get_isotope(self, f, name, ncnts):
iso = Isotope(name, 'Detector1')
iso.name = name
rs = (f.next() for i in xrange(ncnts))
xs, ys = zip(*((float(r[0]), float(r[1])) for r in rs))
iso.xs = xs
iso.ys = ys
return iso
def _get_isotope(self, f, name, ncnts, bk):
iso = Isotope(name, 'Detector1')
iso.set_ublank(bk)
iso.name = name
iso.set_fit('linear')
iso.set_fit_error_type('SEM')
rs = (next(f) for i in range(ncnts))
ys, xs = list(zip(*((float(r[0]), float(r[1])) for r in rs)))
iso.xs = array(xs)
iso.ys = array(ys)
return iso
def _get_isotope(self, f, name, ncnts, bk):
iso = Isotope(name, 'Detector1')
iso.set_ublank(bk)
iso.name = name
iso.set_fit('linear')
iso.set_fit_error_type('SEM')
rs = (next(f) for i in range(ncnts))
ys, xs = list(zip(*((float(r[0]), float(r[1])) for r in rs)))
iso.xs = array(xs)
iso.ys = array(ys)
return iso
def get_analysis_import_spec(self, delimiter=None):
f = self.file_gen(delimiter)
pspec = self.new_persistence_spec()
ident = os.path.splitext(os.path.basename(self.path))[0]
pspec.run_spec.uuid = ident
# labnumber = '100000'
# irradiation_level = 'A'
# irradiation_position = 1
aliquot = int(ident[-4:])
# pspec.run_spec.labnumber = labnumber
pspec.run_spec.aliquot = aliquot
# pspec.run_spec.irradiation_level = irradiation_level
# pspec.run_spec.irradiation_position = irradiation_position
version = next(f)
ncycles = next(f)
total_analysis_time = get_int(f, 1)
start = next(f)
end = next(f)
nzeros = next(f)
nanalysis_cycles = next(f)
toffset = next(f)
pspec.timestamp = datetime.strptime(get_next(f, 1), '#%Y-%m-%d %H:%M:%S#')
collector_gains = next(f)
print('casd', collector_gains)
# int_posts = [next(f) for i in range(41)]
# print(int_posts)
while 1:
i = next(f)
if i[0] == '"Number of peaks centred"':
npeakscentered = i[1]
break
print('fff', npeakscentered)
ndeflectors = next(f)[1]
source_ht = next(f)
half_plate_v = next(f)
trap = next(f)
trap_voltage = next(f)
repeller = next(f)
filament_v = next(f)
delta_hp = next(f)
z_lens = next(f)
delta_z = next(f)
max_current = next(f)
quad_1 = next(f)
cubic_1 = next(f)
lin_1 = next(f)
q18_cor = next(f)
q19_cor = next(f)
quad_2 = next(f)
cubic_2 = next(f)
lin_2 = next(f)
q28_cor = next(f)
q29_cor = next(f)
suppressor = next(f)
Deflect_IC1 = next(f)
Filter_IC0 = next(f)
Deflect_IC0 = next(f)
Deflect_IC2 = next(f)
Filter_IC3 = next(f)
Deflect_IC3 = next(f)
mdfpath = next(f)
analysis_type_info = next(f)
# trim off quotes
ati = analysis_type_info[1:-1]
ati = ati.split(' ')
if ati[0] == 'Blank':
pspec.run_spec.labnumber = 'ba-01'
pspec.run_spec.irradiation = 'NoIrradiation'
pspec.run_spec.irradiation_level = 'A'
pspec.run_spec.irradiation_position = 1
elif ati[0] == 'Air':
a = int(ati[1])
pspec.run_spec.labnumber = 'a-{:02n}'.format(a)
pspec.run_spec.irradiation = 'NoIrradiation'
pspec.run_spec.irradiation_level = 'A'
pspec.run_spec.irradiation_position = a + 1
ics = next(f)
discs = next(f)
[_ for _ in range(9)]
signals = []
for _ in range(total_analysis_time):
line = next(f)
_type = int(line[-1])
if _type:
signals.append([float(li) for li in line])
signals = array(signals)
with open(self.nice_path, 'r') as nice:
isotopes = {}
parser = NiceParser(signals)
for line in nice:
try:
lhs, rhs = list(map(str.strip, line.split('=')))
except ValueError:
continue
if lhs.startswith('Result'):
break
parser.set_tokens(rhs.split(' '))
ret, det_idx = parser.exp()
iso = Isotope(lhs, 'IC{}'.format(det_idx))
iso.name = lhs
iso.xs = ret.xs
iso.ys = ret.ys
isotopes[lhs] = iso
pspec.isotope_group = IsotopeGroup(isotopes=isotopes)
return pspec
