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'
set output outfile

set timestamp offset graph 1.16,0.7 rotate font 'Verdana,6'
set label 'plotting_scripts/phase_space_to_hist.gnu' at graph 1.02,0.02 rotate left font 'Verdana,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 'Verdana,6' noenhance

filename = 'bin_average_rms_7.dat'     #(4->7 for acos(s.p) instead of atan(sx/sz)
fit_err = '2D_hist_fit_coef_err_7.dat'
load fit_err
f(x) = a*x+b

#set xrange [-0.0:0.05]
#set xrange [-0.05:0.05]   #####
set xlabel 'x [m]'
#set yrange [-0.015:0.015]
set ylabel 'spin angle [rad]'
scale = 1
a=0.;b=0.
print filename
fit f(x)  filename u ($1/scale):($5>0 ? $2:1/0):(1./$5**.5) yerror via a,b

#set title filename.' - '.file noenhance
set label 'fit binned data' at graph 0.2, 0.95
set label at graph 0.2,0.9 'slope= '.sprintf("%5.4f",a).'{/Symbol \261}'.sprintf("%5.4f",a_err).' [rad/m]' font ',10'
set label at graph 0.2,0.85 'offset= '.sprintf("%5.4f",b).'{/Symbol \261}'.sprintf("%5.4f",b_err).'[rad]' font ',10'

plot  filename u ($1/scale):2:4 w err not 'data', f(x) not
unset label
set label 'plotting_scripts/phase_space_to_hist.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set label at graph 0.0,1.02 name font 'Verdana,6' noenhance
set label 'fit unbinned data' at graph 0.2,0.95
set label at graph 0.2,0.9 'slope= '.sprintf("%5.4f",coef1).'{/Symbol \261}'.sprintf("%5.4f",coef1_err).' [rad/m]' font ',10'
set label at graph 0.2,0.85 'offset= '.sprintf("%5.4f",coef0).'{/Symbol \261}'.sprintf("%5.4f",coef0_err).'[rad]' font ',10'
g(x) = coef0 + coef1*x
plot  filename u ($1/scale):2:4 w err not 'data', g(x) not


unset label
set label 'plotting_scripts/phase_space_to_hist.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set label at graph 0.0,1.02 name font 'Verdana,6' noenhance
filename='bin_average_rms_8.dat'
fit_err = '2D_hist_fit_coef_err_8.dat'
load fit_err
#set xrange [-0.015:0.015]
set xlabel "x' [rad]"
#set yrange [0.6:]
set ylabel 'spin angle [rad]'

a=1;b=1;
print filename
fit f(x)  filename u ($1):($5>0 ? $2:1/0):(1./$5**.5) yerror via a,b
set label 'fit binned data' at graph 0.2,0.95
set label at graph 0.2,0.9 'slope= '.sprintf("%5.4f",a).'{/Symbol \261}'.sprintf("%5.4f",a_err) font ',10'
set label at graph 0.2,0.85 'offset= '.sprintf("%5.4f",b).'{/Symbol \261}'.sprintf("%5.4f",b_err).' [rad]' font ',10'
plot  filename u ($1):2:4 w err not, f(x) not

unset label
set label 'fit unbinned data' at graph 0.2,0.95
set label at graph 0.2,0.9 'slope= '.sprintf("%5.4f",coef1).'{/Symbol \261}'.sprintf("%5.4f",coef1_err) font ',10'
set label at graph 0.2,0.85 'offset= '.sprintf("%5.4f",coef0).'{/Symbol \261}'.sprintf("%5.4f",coef0_err).'[rad]' font ',10'
g(x) = coef0 + coef1*x
plot  filename u ($1/scale):2:4 w err not 'data', g(x) not


#--------------------------------
if(number == 10){

unset label
set label 'plotting_scripts/phase_space_to_hist.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set label at graph 0.0,1.02 name font 'Verdana,6' noenhance
filename='bin_average_rms_10.dat'
fit_err = '2D_hist_fit_coef_err_10.dat'
load fit_err
#set xrange [-0.015:0.015]
set xlabel "amplitude [m-rad]"
#set yrange [0.6:]
set ylabel 'spin angle [rad]'

a=1;b=1;
print filename
fit f(x)  filename u ($1):($5>0 ? $2:1/0):(1./$5**.5) yerror via a,b
set label 'fit binned data' at graph 0.2,0.95
set label at graph 0.2,0.9 'slope= '.sprintf("%5.4f",a).'{/Symbol \261}'.sprintf("%5.4f",a_err) font ',10'
set label at graph 0.2,0.85 'offset= '.sprintf("%5.4f",b).'{/Symbol \261}'.sprintf("%5.4f",b_err).' [rad]' font ',10'
plot  filename u ($1):2:4 w err not, f(x) not

unset label
set label 'fit unbinned data' at graph 0.2,0.95
set label at graph 0.2,0.9 'slope= '.sprintf("%5.4f",coef1).'{/Symbol \261}'.sprintf("%5.4f",coef1_err) font ',10'
set label at graph 0.2,0.85 'offset= '.sprintf("%5.4f",coef0).'{/Symbol \261}'.sprintf("%5.4f",coef0_err).'[rad]' font ',10'
g(x) = coef0 + coef1*x
plot  filename u ($1/scale):2:4 w err not 'data', g(x) not
}
#----------------------------------
unset label
set label 'plotting_scripts/phase_space_to_hist.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set label at graph 0.0,1.02 name font 'Verdana,6' noenhance
filename='bin_average_rms_6.dat'
fit_err = '2D_hist_fit_coef_err_6.dat'
load fit_err
set xrange [-0.008:0.008]
set xlabel '{/Symbol D}p/p'
#set yrange [-0.04:0.04]
set ylabel 'spin angle [rad]'

a=1.;b=1.
print filename
#fit f(x)  filename u ($1):($5>0 ? $2:1/0):(1./$5**.5) yerror via a,b

#set title filename noenhance
set label 'fit binned data' at graph 0.2,0.25
set label at graph 0.2,0.2 'slope= '.sprintf("%5.4f",a).'{/Symbol \261}'.sprintf("%5.4f",a_err).' [rad]' font ',10'
set label at graph 0.2,0.15 'offset= '.sprintf("%5.4f",b).'{/Symbol \261}'.sprintf("%5.4f",b_err) font ',10'

plot  filename u ($1):2:4 w err not 'data', f(x) not

unset label
set label 'fit unbinned data' at graph 0.2,0.25
set label at graph 0.2,0.2 'slope= '.sprintf("%5.4f",coef1).'{/Symbol \261}'.sprintf("%5.4f",coef1_err).' [rad]' font ',10'
set label at graph 0.2,0.15 'offset= '.sprintf("%5.4f",coef0).'{/Symbol \261}'.sprintf("%5.4f",coef0_err) font ',10'
g(x) = coef0 + coef1*x
plot  filename u ($1/scale):2:4 w err not 'data' #, g(x) not

print outfile