# gnuplot -e 'filename=
# output from g2_tracking
#set term x11
set terminal pdf enhanced fontscale 0.75 size 6.0in, 3.5in
name = 'Beam_moments_scraping_vertical.pdf'
set output name
#d0=0.0025
#xlow=sprintf("%i",xmin)
#xhigh=sprintf("%i",xmax)
#print xlow, xhigh
#file ='EnergyDecoherence_x_'.xlow.'-'.xhigh.'-'.TimeBin.'.pdf'
#print file
#set output file
#filename ='xaverage_energy_'.TimeBin.'.dat'
#print filename
#set xrange [xmin:xmax]
set ylabel 'Vertical beam width ({/Symbol s}) [mm]'
set xlabel 'Time [{/Symbol m}s]'
#set label at graph 0.68,0.9 'Time Bin = '.TimeBin
set label 'plotting scrips/Beam\_moments\_scraping.gnu' at graph 1.02,0.02 rotate left font 'Verdana,8'
#delta0 = sprintf("%.4f",d0)
#title='{/Symbol D}_0 = '.delta0
set grid
 plot '20160913_155224/Beam_moments_tbt.dat' u ($5*1.e6):($9**.5*1000) w l t 'scraping off - 20160913\_155224', '20160913_084008/Beam_moments_tbt.dat' u ($5*1.e6):($9**.5*1000) w l t 'scraping on - 20160913\_084008'
name = 'Beam_moments_scraping_horizontal.pdf'
set output name
set ylabel 'Horizontal_beam width ({/Symbol s}) [mm]'
 plot '20160913_155224/Beam_moments_tbt.dat' u ($5*1.e6):($6**.5*1000) w l t 'scraping off - 20160913\_155224', '20160913_084008/Beam_moments_tbt.dat' u ($5*1.e6):($6**.5*1000) w l t 'scraping on - 20160913\_084008'