set term postscript color enhanced "Times-Roman" 14
set output "capture_vs_hkick.ps"
set xlabel "kick [T]"
set ylabel "number of muons"
set multiplot
set xtics 0,0.002,0.275
set origin 0.,0.
set size 0.5, 0.5
plot "fort.32" u 19:16 w l lt -1 t 'muons captured', "" u 19:($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 19:4 w l lt -1 t 'CBO amplitude', "" u 19:10 w l lt 1 t "modulation depth" 


set origin 0.5,0.5
set ylabel "Arb"
plot "fort.32" u 19:($16/$4) w l lt -1 t 'muons/cbo', "" u 19:($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
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