#set term postscript color enhanced "Times-Roman" 26
set terminal pdf enhanced fontscale 0.75 size 6.0in, 4.0in
outfile = "Vparam2_compare3_measured.pdf"
set output outfile
#set title title noenhance
#dir1 = 'kick_scan_20150328/20171118_232900' 
dir1 ='kick_scan_20171118_halftail/20171119_124116'
dir2 = 'kick_scan_20171113/20171114_101423'
dir3 = 'kick_scan_20171119_05_tail/20171119_235059'
dir4 = 'kick_scan_20150328/20171115_084156'
dir5='kick_scan_20171113/20180318_192400'
dir6='kick_scan_20171113/20180320_131119' #/20180319_092127' #kicker1data_flattop.dat 
load dir1."/VparamDependence.dat"
nmuons5=nmuons
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 at graph 0.,1.02 dir1  font 'Verdana, 6' noenhance
set label at graph 0.,1.04 dir2  font 'Verdana, 6' noenhance
set label at graph 0.,1.06 dir3  font 'Verdana, 6' noenhance
set label at graph 0.,1.08 dir4  font 'Verdana, 6' noenhance
set label at graph 0.,1.1 dir5  font 'Verdana, 6' noenhance
set label 'plotting scripts/Vparam2_compare.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set timestamp offset graph 1.07,0.7 rotate font 'Verdana,6'
print column1
plot dir2."/VparamDependence.dat" u column1:($16/nmuons)*100:(($16)**0.5/nmuons*100) i 1:1 w errorlines lw 1.5 lt -1 t 'measured 11-17',dir1."/VparamDependence.dat" u column1:($16/nmuons)*100:(($16)**0.5/nmuons*100) i 1:1 w errorlines lw 1.5 lt 2 pt 5 ps .5 t 'halftail',dir3."/VparamDependence.dat" u column1:($16/nmuons)*100:(($16)**0.5/nmuons*100) i 1:1 w errorlines lw 1.5 lt 3 pt 5 ps .5 t '5% tail',dir4."/VparamDependence.dat" u column1:($16/nmuons)*100:(($16)**0.5/nmuons*100) i 1:1 w errorlines lw 1.5 lt 4 pt 5 ps .5 t 'ideal',dir5."/VparamDependence.dat" u column1:($16/nmuons)*100:(($16)**0.5/nmuons*100) i 1:1 w errorlines lw 1.5 lt 7 pt 6 ps .5 t 'check',dir6."/VparamDependence.dat" u column1:($16/nmuons)*100:(($16)**0.5/nmuons*100) i 1:1 w errorlines lw 1.5 lt 6 pt 6 ps .5 t 'flattop'
reset
set ylabel "CBO amplitude [mm]" offset -1,0,0
set ytics nomirror
set label at graph 0.,1.02 dir1  font 'Verdana, 6' noenhance
set label at graph 0.,1.04 dir2  font 'Verdana, 6' noenhance
set label at graph 0.,1.06 dir3  font 'Verdana, 6' noenhance
set label at graph 0.,1.08 dir4  font 'Verdana, 6' noenhance
set label 'plotting scripts/Vparam2_compare.gnu' at graph 1.02,0.02 rotate left font 'Verdana,6' noenhance
set timestamp offset graph 1.07,0.7 rotate font 'Verdana,6'
set yrange [0:] 
#set ytics 0,5
plot dir2."/VparamDependence.dat" u column1:4:($4*($16/nmuons)**0.5)  i  1:1  w errorl lw 1.5 lt 1 t 'measured 11-17' ,dir1."/VparamDependence.dat" u column1:4:($4*($16/nmuons)**0.5)  i  1:1  w errorl lw 1.5 lt 2 pt 5 ps .5 t "ideal" ,dir3."/VparamDependence.dat" u column1:4:($4*($16/nmuons)**0.5)  i  1:1  w errorl lw 1.5 lt 3 pt 5 ps .5 t "5% tail" ,dir4."/VparamDependence.dat" u column1:4:($4*($16/nmuons)**0.5)  i  1:1  w errorl lw 1.5 lt 4 pt 5 ps .5 t "ideal" 

print outfile
#plot "VparamDependence.dat" u column1:16 w l lt -1 t 'muons captured', "" u column1:15 w l lt 1 t "muons exiting inflector" 

# 'measured 11-17',\
#t 'measured 11-17',\ 
#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