set term pdf enhanced  fontscale 0.75 size 6.0in, 4.0in
set output "Vparam.pdf"
load "VparamDependence.dat"
set key left Left
set label 'mytest/quad_scan_bradial/20170616_223024/' at graph 0.55,1.02 font 'Verdana,6' noenhance
set label '/plotting_scripts/Vparampdf.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set timestamp offset graph 0.96,0.7 rotate font 'Verdana,6'
set ylabel "captured [%]" offset 2,0
set ytics nomirror
set y2label "CBO amplitude [mm]" offset -1,0
set y2tics nomirror 
plot "VparamDependence.dat" u column1:(100*$16/nmuons):(($16)**0.5/nmuons*100) axis x1y1 w errorlines lw 1.5 lt -1 t '% muons captured' , \
"" u column1:4:($4*($16/nmuons)**0.5)  axis x1y2 w errorl lw 1.5 lt 1 t "CBO amplitude"
#plot "VparamDependence.dat" u column1:16 w l lt -1 t 'muons captured', "" u column1:15 w l lt 1 t "muons exiting inflector" 
unset y2label
unset y2tics
unset terminal

#plot "gnu_input.dat" u 1:2 :3 w labels left font ",10" not