unset multiplot
!rm VparamE821_E989CBO.ps
set term postscript color enhanced "Times-Roman" 22
set output "VparamE821_E989CBO.ps"

column1=19
nmuons=10000
set border linewidth 1.5
set xlabel 'Kick [G]"
set ylabel "CBO Amplitude [mm]" offset 2,0,0
set xtics nomirror
set yrange [0:]
set origin 0.,0.
set style arrow 1 head filled size screen 0.025,30,45 lt 1 lw 3 lc 0
set arrow from 166,29 to 166,33 as 1
set label at 125.,30. "E821" # font "TimesRoman" 20
set label at 125.,27. "Voltage"
set label at 125.,24. "Maximum"
plot "821kick-20150614/VparamDependence.dat" u column1:4 i 1:1 w l lw 2 lt 1 lc 1 t 'E821 kicker',\
 "" u column1:4 i 2:2 w l lw 2 lt 2 lc 1  not '% muons captured' ,\
 "127cm_table_20150507/Kick1-3Scan/VparamDependence.dat" u column1:4 i 1:3 w l lw 2 lt 1 lc 3 t 'E989 kicker' ,\
 "127cm_table_20150507/Kick1-3Scan/VparamDependence.dat" u column1:4:($4*($16/nmuons)**0.5) i 2:2  w errorlines lw 2 lt 1 lc 3 not 'E989 kicker'

#unset y2label
#unset y2tics
#set origin 0.,0.5
#set size 0.5,0.5
#set ylabel "Amplitude [mm]"
#plot "VparamDependence.dat" u column1:4 w l lt -1 t 'CBO amplitude', "" u column1:10 w l lt 1 t "modulation depth" 


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