set term postscript color enhanced "Times-Roman" 14
set output "capture_vs_betay.ps"
set xlabel "beta y [m]"
set ylabel "number of muons"
set multiplot
set xtics 0,1,20
set origin 0.,0.
set size 0.5, 0.5
plot "fort.32" u 20:16 w l lt -1 t 'muons captured', "" u 20:($15-$16) w l lt 1 t "muons lost in ring" 

set origin 0.,0.5
set size 0.5,0.5
set ylabel "Amplitude [mm]"
plot "fort.32" u 20:4 w l lt -1 t 'CBO amplitude', "" u 20:10 w l lt 1 t "modulation depth" 


set origin 0.5,0.5
set ylabel "Arb"
plot "fort.32" u 20:($16/$4) w l lt -1 t 'muons/cbo', "" u 20:($16/($4*$10)) w l lt 1 t "muons/(cbo X mod)" 
set origin 0.5,0.
unset border
unset xtics
unset ytics
unset ylabel
unset xlabel
set yrange [0:10]
set xrange [1.:3]
set label font ",10"
plot "gnu_input.dat" u 1:2 :3 w labels left font ",10" not