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 = 'differential_decay_dependencies.pdf'
set output outfile

set timestamp offset graph 1.16,0.7 rotate font 'Verdana,6'
set label 'plotting_scripts/differential_decay_dependencies.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

set print 'differential_decay_dependencies.dat'

do for [i=1:5]{



if(i ==1){
number=6
set xlabel  '{/Symbol D}p/p'
set ylabel  '{/Symbol f} [rad]'
word ='phi vs Dp/p'
units='[rad]'
}
if(i ==2){
number=7
set xlabel  'x [m]'
set ylabel  '{/Symbol f} [rad]'
word ='phi vs x'
units='[rad/m]'
}
if(i ==3){
number=8
set xlabel  "x' [rad]"
set ylabel  '{/Symbol f} [rad]'
word ="phi vs x'"
units=''
}
if(i ==4){
number='13'
set xlabel  "y [m]"
set ylabel  '{/Symbol f} [rad]'
word ="phi vs y"
units='[rad/m]'
}
if(i ==5){
number='12'
set xlabel  "y [m]"
set ylabel  '{/Symbol f} [rad]'
word ="phi vs y'"
units=' '
}

filename = 'bin_average_rms_'.number.'.dat'
fit_err = '2D_hist_fit_coef_err_'.number.'.dat'

load fit_err
print word
print 'p0=',coef0,'+/-',coef0_err,'[rad]'
print 'p1=',coef1,'+/-',coef1_err,units
print ' '

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

print outfile