set term postscript color enhanced "Times-Roman" 14
set output "capture_vs_width.ps"
set xlabel "width [ns]"
set ylabel "number of muons"
set multiplot
#set xtics 0,100,1000
#set xrange [10:300]
set origin 0.,0.
set size 0.5, 0.5
plot "fort.32" u ($21*1.e9):16 w l lt -1 t 'muons captured', "" u ($21*1.e9):($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 ($21*1.e9):4 w l lt -1 t 'CBO amplitude', "" u ($21*1.e9):10 w l lt 1 t "modulation depth" 


set origin 0.5,0.5
set ylabel "Arb"
plot "fort.32" u ($21*1.e9):($16/$4) w l lt -1 t 'muons/cbo', "" u ($21*1.e9):($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]
unset xlabel
set xrange [1.:3]
set label font ",10"
plot "gnu_input.dat" u 1:2 :3 w labels left font ",10" not