set term postscript color enhanced "Times-Roman" 14
set output "capture_vs_kick_four.ps"
set xlabel "kick [G]"
set ylabel "number of muons"
set multiplot
#set xtics 0,100,1000
#set xrange [150:300]
set origin 0.,0.
set title "Dependence on Scattering Through Inflector End Coils"
set label at 120,9000 "Number lost in ring"
set label at 120,2500 "Number stored"
plot  "20140317-1/fort.32" u ($19*10000):16 w lp ls 1.5 pt 6 lc 3 t 'no scattering', "" u ($19*10000):($15-$16) w l lc 3 not \
, "20140312/fort.32" u ($19*10000):16 w lp ps 1.5 pt 5 lc 1 t 'upstream scattering' , "" u ($19*10000):($15-$16) w l lc 1 not \
,   "20140313/fort.32" u ($19*10000):16 w lp ps 1.5 pt 4 lc 4 t 'downstream scattering', "" u ($19*10000):($15-$16) w l lc 4 not \
,  "20140316/fort.32" u ($19*10000):16 w lp ps 1.5 lc 7 t 'upstream and downstream scattering', "" u ($19*10000):($15-$16) w l lc 7 not