Commit 91e011aa authored by Plaszczynski Stephane's avatar Plaszczynski Stephane
Browse files

tricks to make profile look better

parent 9a6b9ca2
......@@ -15,13 +15,19 @@ l=label.split()
d=np.loadtxt(file,skiprows=1)
ichi2=l.index('chi2')
ivar=ichi2-1 # the scanned variable index is always below the chi2 value
ivar=ichi2-1 # the scanned variable index is always one index below chi2
chi2=d[:,ichi2]
var=d[:,ivar]
f=interp1d(var,chi2,kind='quadratic')
#no need to keep points above dchi2>5
dchi2=chi2-chi2.min()
w=dchi2<5
chi2=chi2[w]
var=var[w]
x=np.linspace(var[0],var[-1],100)
f=interp1d(var,chi2,kind='cubic')
x=np.linspace(var[0],var[-1],500)
y=f(x)
ymin=y.min()
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment