# 
#  Compares result from Fourier Filter with Path from exp2feff

#  Read chi(k) data
read_data (file = ../data/cu50k.chi, group=cu, type=chi)

# Fourier filter
fftf(cu.chi, kmin = 1.5, kmax = 18.5, dk = 1.0, kweight = 2)
fftr(real=cu.chir_re, imag=cu.chir_im, rmin = 1.3, rmax = 3.0)

newplot cu.q, cu.chiq_im, color=blue
path(1, file=feff_exp.dat, degen=12, s02=1.0)

ff2chi(1, group=fexp)

plot fexp.k, fexp.chi*fexp.k^2 , color=red
