set term pdf fontscale 0.5 size 5in,3in
#filename='bin_average_rms_7.dat'       #(4->7 for acos(s.p) instead of atan(sx/sz)
#load filename
#outfile = file.'.pdf'
outfile = 'spin_debug.pdf'
set output outfile

set timestamp offset graph 1.175,0.7 rotate font ',6'
set label 'plotting_scripts/spin_debug.gnu' at graph 1.02,0.02 rotate left font ',6' noenhance
!pwd > directory
!sed -n -e 's/^/name="/' -e 's/$/"/p' directory > name
load 'name'
print 'directory='.name
set label at graph 0.0,1.02 name font ',6' noenhance

load 'domega_time_dependence.dat'
print 'dg_dt=', delta_2
set xlabel 'time [us]'
set ylabel 'Fractional momentum offset [%]'
h(x) = delta_1 + delta_2 * x
set xrange [0:1000]
plot h(x)*100 not

filename = '2D_hist_fit_coef_err_4.dat'
load filename
set label 3  filename at graph 0.02,0.95   noenhance
f(x) = coef0 + coef1*x+coef2*x**2+coef3*x**3
set xlabel 'fractional momentum offset'
set ylabel 'displacement [m]'
set xrange [-0.004:0.004]
plot f(x)
#g(x) = 
f(x) = coef0 + coef1*x+coef2*x**2+coef3*x**3
set xrange [-0.004:0.004]
plot f(x)
unset label 3

filename = '2D_hist_fit_coef_err_5.dat'
load filename

set ylabel "angle x' [rad]"
set label 3  filename at graph 0.02,0.95 noenhance
set xrange [-0.004:0.004]
plot f(x)
unset label 3
set timestamp offset graph 1.2,0.7 rotate font ',6'
g(x) = (delta_omega_0 + delta_omega_1* x + delta_omega_2 * x**2)*delta_2/omega_a
g_err(x) = sqrt((do0_err)**2 + (do1_err * x/2.)**2 + (do2_err * x**2/3)**2) * delta_2/omega_a
set xrange [0:1000]
set xlabel 'time [{/Symbol m}s]'
set ylabel '{/Symbol Dw_a/w_a} [ppb]'
plot g(x)*1.e9 not
plot (g(x)+g_err(x))*1.e9 lc -1 dt 2 not , g(x)*1.e9 not, (g(x)-g_err(x))*1.e9 lc -1 dt 2 not
set timestamp offset graph 1.15,0.7 rotate font ',6'
gavg(x) = (delta_omega_0 + delta_omega_1* x/2. + delta_omega_2/3 * x**2)*delta_2/omega_a
set ylabel '{/Symbol Dw_a/w_a} [ppb]'
set key left Left
plot g(x)*1.e9 t '{/Symbol Dw_a/w_a}', gavg(x)*1.e9 t '<{/Symbol Dw_a/w_a}>_T'



print outfile