def test_zeta_vs_cross_section(self): s = self.signal factors = [3, 5] method = 'zeta' intensities = s.get_lines_intensity() zfactors = utils_eds.edx_cross_section_to_zeta([3, 5], ['Al', 'Zn']) factors2 = utils_eds.zeta_to_edx_cross_section(zfactors, ['Al', 'Zn']) np.testing.assert_allclose(factors, factors2, atol=1e-3) res = s.quantification(intensities, method, factors=utils_eds.edx_cross_section_to_zeta( [22.402, 21.7132], ['Al', 'Zn'])) res2 = s.quantification(intensities, method='cross_section', factors=[22.402, 21.7132]) np.testing.assert_allclose(res[0][0].data, res2[0][0].data, atol=1e-3) np.testing.assert_allclose(res[0][0].data, np.ones((2, 2)) * 36.2969, atol=1e-3)
def test_edx_cross_section_to_zeta(self): cs = [3, 6] elements = ['Pt', 'Ni'] res = utils_eds.edx_cross_section_to_zeta(cs, elements) np.testing.assert_allclose(res, [1079.815272, 162.4378035], atol=1e-3)