def read(file,
bad_mca_idx=[],
total=True,
align=True,
correct=True,
tau=None,
det_idx=0,
emin=-1.0,
emax=-1.0,
xrf_params={},
lines=None,
fmt='CARS'):
"""
Read detector files
Examples:
---------
>>m = xrf.read(file="file_name",bad_mca_idx=[],
total=True,align=True,tau=None)
Notes:
------
* xrf_params should have the same format as returned
by Xrf.get_params
* if xrf_params == None: returns a list of med objects
otherwise: returns a list of xrf objects (default)
"""
med = read(file,
bad_mca_idx=bad_mca_idx,
total=total,
align=align,
correct=correct,
tau=tau,
fmt=fmt)
if med == None: return None
xrf = med2xrf(med,
xrf_params=xrf_params,
lines=lines,
det_idx=det_idx,
emin=emin,
emax=emax)
return xrf
def read(
file,
bad_mca_idx=[],
total=True,
align=True,
correct=True,
tau=None,
det_idx=0,
emin=-1.0,
emax=-1.0,
xrf_params={},
lines=None,
fmt="CARS",
):
"""
Read detector files
Examples:
---------
>>m = xrf.read(file="file_name",bad_mca_idx=[],
total=True,align=True,tau=None)
Notes:
------
* xrf_params should have the same format as returned
by Xrf.get_params
* if xrf_params == None: returns a list of med objects
otherwise: returns a list of xrf objects (default)
"""
med = read(file, bad_mca_idx=bad_mca_idx, total=total, align=align, correct=correct, tau=tau, fmt=fmt)
if med == None:
return None
xrf = med2xrf(med, xrf_params=xrf_params, lines=lines, det_idx=det_idx, emin=emin, emax=emax)
return xrf
"""
lines = []
for pk in self.xrf[0].peaks:
lines.append(pk.label)
self.lines = lines
self.peaks = {}
for l in lines:
p = peak_areas(self.xrf, l)
self.peaks[l] = p
##############################################################################
##############################################################################
if __name__ == "__main__":
from matplotlib import pyplot
xrf = read(file='_test.xrf', bad_mca_idx=[0, 2, 13], emin=4., emax=9.)
pyplot.plot(xrf.get_energy(), xrf.get_data(), '-k')
#pyplot.show()
#
xrf.init_lines(['Fe ka', 7.12])
#xrf.init_bgr()
#
xrf.calc()
pyplot.plot(xrf.get_energy(), xrf.predicted, '-r')
#
xrf.fit()
en = xrf.get_energy()
pyplot.plot(en, xrf.predicted, '-b')
#
#for peak in xrf.peaks:
# pyplot.plot(en,1+peak.calc(en),'.')
lines = []
for pk in self.xrf[0].peaks:
lines.append(pk.label)
self.lines = lines
self.peaks = {}
for l in lines:
p = peak_areas(self.xrf, l)
self.peaks[l] = p
##############################################################################
##############################################################################
if __name__ == "__main__":
from matplotlib import pyplot
xrf = read(file="_test.xrf", bad_mca_idx=[0, 2, 13], emin=4.0, emax=9.0)
pyplot.plot(xrf.get_energy(), xrf.get_data(), "-k")
# pyplot.show()
#
xrf.init_lines(["Fe ka", 7.12])
# xrf.init_bgr()
#
xrf.calc()
pyplot.plot(xrf.get_energy(), xrf.predicted, "-r")
#
xrf.fit()
en = xrf.get_energy()
pyplot.plot(en, xrf.predicted, "-b")
#
# for peak in xrf.peaks:
# pyplot.plot(en,1+peak.calc(en),'.')
