reset
set term eps font "arial,12" enhanced size 8.6cm,10cm
set output "bristol_rho_vs_current.eps"

set multiplot 

# CODE FOR WATERFALL PLOT AT THE END OF TEXT IF NEEDED 

# alternatively only tracking the peak value
set size 1,0.5
set origin 0,0.5
set ylabel "power fraction"
set xrange [48:52]
set yrange [-0.05:1.05]
set ytics 0,0.2,1
set ytics auto format "%.2f"
set format x ""
set label 1 "r_{/Symbol a}" at 48.1,0.95
set label 2 "r_{/Symbol b}" at 48.1,0.1


# choice of colors
plot '780p07_780p34_ratio_points_removed.dat' using 1:2 with points pointtype 5 ps 0.5 linetype rgb "#000000" title "",'' using 1:3 with points pointtype 6 ps 0.5 linetype rgb "#000000" title ""



# fraction of coherent against current
set size 1,0.55
set origin 0,0
set xlabel "current (mA)"
set ylabel "{/Symbol r}" enhanced
set xrange [48:52]
set yrange [-0.05:1.05]
set ytics auto format "%.2f"
set bars large
unset format x
unset label

#plot  '780p07_ratio_max_points_removed.dat' using 1:2 with points pointtype 5 ps 0.5 linetype rgb "#000000" title "",'780p34_ratio_max_points_removed.dat' using 1:2 with points pointtype 6 ps 0.5 linetype rgb "#000000" title "",'rho_vs_current.dat' using 1:2:3:4 with yerr linestyle 7 ps 0.2 linetype rgb "#0000FF" title "", '' using 1:5:6:7 with yerr linestyle 7 ps 0.2 linetype rgb "#FF0000" title ""

plot 'rho_vs_current.dat' using 1:2:3:4 with yerr linestyle 7 ps 0.3 lw 2 linetype rgb "#0000FF" title "", '' using 1:5:6:7 with yerr linestyle 7 ps 0.3 lw 2 linetype rgb "#FF0000" title ""


unset multiplot


# g2X waterfall plot
#set size 1,0.5
#set origin 0,0.5
#set ylabel "wavelength [nm]"
#set xlabel "current [mA]"
#set cblabel "normalised intensity"
#set xrange [48:52]
#set yrange [779.8:782.2]
#set ytics 779.8,0.1,780.2
#set cbrange [0:0.4]
#set cbtics 0,0.1,0.4

# choice of colors
#set palette negative rgb 34,35,36
#plot "current_wavelength_intensity.dat" using 1:2:3 with image failsafe title ""