#set term postscript color enhanced "Times-Roman" 26
set terminal pdf enhanced fontscale 0.75 size 6.0in, 4.0in
outfile = "Vparam2_compare".ncompare.".pdf"
set output outfile
#set title title noenhance
#dir1 = 'kick_scan_20150328/20171118_232900' 
#ncompare = 4
#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'

reset
set border linewidth 1.5
set key top left
set ylabel "% captured" offset 2,0,0
#set multiplot		
#set xtics nomirror
set yrange [0:4]
set ytics 0,1
set mytics 2
if (nokey==1){set key off}
if(ncompare >= 1){
    set label at graph 0.,1.02 dir1  font 'Verdana, 6' noenhance
    load dir1."/VparamDependence.dat"
    nmuons1=nmuons
}
if(ncompare >= 2){
    set label at graph 0.,1.04 dir2  font 'Verdana, 6' noenhance
    load dir2."/VparamDependence.dat"
    nmuons2=nmuons
}
if(ncompare >= 3){
    set label at graph 0.,1.06 dir3  font 'Verdana, 6' noenhance
    load dir3."/VparamDependence.dat"
    nmuons3=nmuons
}
if(ncompare >= 4){
    set label at graph 0.,1.08 dir4  font 'Verdana, 6' noenhance
    load dir1."/VparamDependence.dat"
    nmuons4=nmuons
}
if(ncompare > 4) {print 'There are too many files'}


set label 'plotting scripts/Vparam2_compare_beta.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
if(ncompare == 1){
plot dir1."/VparamDependence.dat" u column1:($16/nmuons1)*100:(($16)**0.5/nmuons1*100) i 1:1 w errorlines lw 1.5 lt 2 pt 5 ps .5 t t1
}
if(ncompare == 2){
plot dir2."/VparamDependence.dat" u column1:($16/nmuons2)*100:(($16)**0.5/nmuons2*100) i 1:1 w errorlines lw 1.5 lt -1 t t2,dir1."/VparamDependence.dat" u column1:($16/nmuons1)*100:(($16)**0.5/nmuons1*100) i 1:1 w errorlines lw 1.5 lt 2 pt 5 ps .5 t t1
}
if(ncompare == 3){
plot dir2."/VparamDependence.dat" u column1:($16/nmuons2)*100:(($16)**0.5/nmuons2*100) i 1:1 w errorlines lw 1.5 lt -1 t t2,dir1."/VparamDependence.dat" u column1:($16/nmuons1)*100:(($16)**0.5/nmuons1*100) i 1:1 w errorlines lw 1.5 lt 2 pt 5 ps .5 t t1,dir3."/VparamDependence.dat" u column1:($16/nmuons3)*100:(($16)**0.5/nmuons3*100) i 1:1 w errorlines lw 1.5 lt 3 pt 5 ps .5 t t3}
if(ncompare == 4){
plot dir2."/VparamDependence.dat" u column1:($16/nmuons2)*100:(($16)**0.5/nmuons2*100) i 1:1 w errorlines lw 1.5 lt -1 t t2,dir1."/VparamDependence.dat" u column1:($16/nmuons1)*100:(($16)**0.5/nmuons1*100) i 1:1 w errorlines lw 1.5 lt 2 pt 5 ps .5 t t1,dir3."/VparamDependence.dat" u column1:($16/nmuons3)*100:(($16)**0.5/nmuons3*100) i 1:1 w errorlines lw 1.5 lt 3 pt 5 ps .5 t t3,dir4."/VparamDependence.dat" u column1:($16/nmuons4)*100:(($16)**0.5/nmuons4*100) i 1:1 w errorlines lw 1.5 lt 4 pt 5 ps .5 t t4}
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
if(ncompare ==1){
plot dir1."/VparamDependence.dat" u column1:4:($4*($16/nmuons1)**0.5)  i  1:1  w errorl lw 1.5 lt 1 t t1
}
if(ncompare ==2){
plot dir2."/VparamDependence.dat" u column1:4:($4*($16/nmuons2)**0.5)  i  1:1  w errorl lw 1.5 lt 1 t t2 ,dir1."/VparamDependence.dat" u column1:4:($4*($16/nmuons1)**0.5)  i  1:1  w errorl lw 1.5 lt 2 pt 5 ps .5 t t1
}
if(ncompare ==3){
plot dir2."/VparamDependence.dat" u column1:4:($4*($16/nmuons2)**0.5)  i  1:1  w errorl lw 1.5 lt 1 t t2 ,dir1."/VparamDependence.dat" u column1:4:($4*($16/nmuons1)**0.5)  i  1:1  w errorl lw 1.5 lt 2 pt 5 ps .5 t t1 ,dir3."/VparamDependence.dat" u column1:4:($4*($16/nmuons3)**0.5)  i  1:1  w errorl lw 1.5 lt 3 pt 5 ps .5 t t3
}
if(ncompare ==4){
plot dir2."/VparamDependence.dat" u column1:4:($4*($16/nmuons2)**0.5)  i  1:1  w errorl lw 1.5 lt 1 t t2 ,dir1."/VparamDependence.dat" u column1:4:($4*($16/nmuons1)**0.5)  i  1:1  w errorl lw 1.5 lt 2 pt 5 ps .5 t t1 ,dir3."/VparamDependence.dat" u column1:4:($4*($16/nmuons3)**0.5)  i  1:1  w errorl lw 1.5 lt 3 pt 5 ps .5 t t3 ,dir4."/VparamDependence.dat" u column1:4:($4*($16/nmuons4)**0.5)  i  1:1  w errorl lw 1.5 lt 4 pt 5 ps .5 t t4
}
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