set term postscript color enhanced "Times-Roman" 14
set output "capture_vs_offset.ps"
set xlabel "initial position offset"
set ylabel "% muons captured"
set y2label "% muons lost in ring"
set multiplot
#set xtics 0,100,1000
#set xrange [150:350]
set origin 0.,0.
set size 0.5, 0.5
set autoscale y
set autoscale y2
set y2tics 
plot "fort.32" every ::1   u ($24*1000):($16/100) axis x1y1 w l lt -1 t 'muons captured' , ""  every ::1  u ($24*1000):(($15-$16)/100) axis x1y2 w l lt 1 t "muons lost in ring" 

unset y2label
unset y2tics
set origin 0.,0.5
set size 0.5,0.5
set ylabel "Amplitude [mm]"
plot "fort.32" u ($24*1000):4 w l lt -1 t 'CBO amplitude', "" u ($24*1000):10 w l lt 1 t "modulation depth" 


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