Kromek GR-1付属のSpectrometer Test Data Sheetを元に
E(keV) = I + k \times ch
の I と k を求める。
Cs137Actual = 879.76
Am241Actual = 76.06
|
|
import numpy, scipy.optimize
xdata = numpy.array([
Am241Actual,
Cs137Actual,
])
ydata = numpy.array([
59.54,
661.66,
])
def energy(x, i, k):
return i + k * x
(popt, pcov) = scipy.optimize.curve_fit(energy, xdata, ydata, p0 = [0.0, 0.75])
|
|
formula_energy = "$ E = (%.3e) + (%.3e)ch $" % tuple(popt)
def fit_energy(x):
return energy(x, *popt)
plot_data = list_plot(zip(xdata, ydata), rgbcolor='gray')
plot_fit = plot(fit_energy, rgbcolor='red', xmax=4096, legend_label=formula_energy)
(plot_data + plot_fit).show(axes_labels=["ch", "E(kev)"])
|
|
