print 'input: filename=dir/all_EndOfTracking_phase_space.dat; events=100; time_off=35e-6'
set terminal pdf enhanced fontscale 0.5 size 6.0in, 4.0in
outputfile = 'all_energy_vs_time.pdf'
print outputfile
set output outputfile

set timestamp offset graph 1.04,0.65 rotate font ',6'
set label 3 'plotting_scripts/energy_vs_time.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set xlabel  'Time [ns]'
set ylabel '{/Symbol D}p/p X 10^{-4}'
f(x) = a+b*x+c*x**2
 fit f(x) 'all_Energy_vs_time_0.dat' u ($3*1.e9):2 via a,b,c
set label at graph 0.05,0.95 'f(x) = a+bx+cx^2'
aword=sprintf("%4.3e",a)
set label at graph 0.05,0.9 'a = '.aword
bword=sprintf("%4.3e",b)
set label at graph 0.05,0.85 'b = '.bword
cword=sprintf("%4.3e",c)
set label at graph 0.05,0.8 'c = '.cword

plot 'all_Energy_vs_time_0.dat' u ($3*1.e9):2 pt 5 ps .25 not,f(x)
print outputfile