<div dir="ltr"><div><div>Hi Johannes, <br><br></div>I guess I hadn't tested the Chantler results below about 1 keV, and certainly not at 100 eV. To be clear, Chantler claims the error bars on these value of "50 to 200%" below 200 eV.   But the spline interpolation is much, much worse than that!<br><br></div>Spline interpolation is definitely preferred at higher energies, and where the grid of tabulated data is fine enough.   But we should better explore when that is valid, and/or automatically use linear interpolation below some energy, and/or allow the user to select the interpolation order. <br><div><div><br></div><div>Anyway, thanks very much.   We should make this better at all energy ranges.<br><br></div><div>--Matt<br></div><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Oct 14, 2015 at 1:11 AM, Dr. Johannes Zellner <span dir="ltr"><<a href="mailto:johannes@zellner.org" target="_blank">johannes@zellner.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="font-size:12pt;font-family:Verdana,Geneva,sans-serif">
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace">Hi,</div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace">the larch xraydb plugin doesn't provide correct numbers for chantler f1 data around edges.</div>
<div style="margin:0;padding:0;font-family:monospace">Apparently this is due to incorrect interpolation in xray/xraydb.py in line 274.</div>
<div style="margin:0;padding:0;font-family:monospace">Below is a plot of "delta" from xraydb_plugin.xray_delta_beta() for Silicon which is terrribly wrong.</div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace"><img src="cid:1444803065561df1f9d6aa6741075083@zellner.org" alt=""></div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace">The more correct "delta" as from chantler is shown in the following plot:</div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace"><img src="cid:1444803065561df1f9d6e1b984800623@zellner.org" alt=""></div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace">--> I'd recommend to replace (line 274 in plugin xray/xraydb.py):</div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace">out = UnivariateSpline(te, ty, s=smoothing)(energy)</div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace">by a simple linear interpolation</div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace">out = np.interp(energy, te, ty)</div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace">The latter yields linear segments between chantler's data, but this is still better than the totally unphysical spline interpolation.</div>
<div style="margin:0;padding:0;font-family:monospace"> </div>
<div style="margin:0;padding:0;font-family:monospace">Best regards,</div><span class="HOEnZb"><font color="#888888">
<div style="margin:0;padding:0;font-family:monospace"><br>
<div>
<div style="margin:0;padding:0;font-family:monospace"><span>-- <br>Johannes</span></div>
</div>
</div>
</font></span></div>
<br>_______________________________________________<br>
Ifeffit mailing list<br>
<a href="mailto:Ifeffit@millenia.cars.aps.anl.gov">Ifeffit@millenia.cars.aps.anl.gov</a><br>
<a href="http://millenia.cars.aps.anl.gov/mailman/listinfo/ifeffit" rel="noreferrer" target="_blank">http://millenia.cars.aps.anl.gov/mailman/listinfo/ifeffit</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature">--Matt Newville <newville at <a href="http://cars.uchicago.edu" target="_blank">cars.uchicago.edu</a>> 630-252-0431<br></div>
</div>