#gnuplot -e 'title ="directory/subdir"' $d9/plotting_scripts/Vparam2.gnu
#set term postscript color enhanced "Times-Roman" 26
set terminal pdf enhanced fontscale 0.75 size 6.0in, 4.0in
set output "Vparam2.pdf"
#set title title noenhance
load "VparamDependence.dat"
reset
set border linewidth 1.5
set key bottom left
set ylabel "% captured" offset 2,0,0
#set multiplot		
#set xtics nomirror
set yrange [0:4]
set ytics 0,1
set mytics 2
set label 'plotting scripts/Vparam2.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6'
print column1
plot "VparamDependence.dat" u column1:($16/nmuons)*100:(($16)**0.5/nmuons*100) i 1:1 w errorlines lw 1.5 lt -1 t '% muons captured'
reset
set ylabel "CBO amplitude [mm]" offset -1,0,0
set ytics nomirror
set yrange [0:] 
#set ytics 0,5
plot "VparamDependence.dat" u column1:4:($4*($16/nmuons)**0.5)  i  1:1  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
#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