set terminal postscript eps color enhanced solid font "Helvetica,22" size 15cm,17cm
set colors classic
set nobars
set style line 1 pt 4 ps 2 lc 'blue' lw 5
set style line 2 pt 6 ps 2 lc 'red' lw 5
set style line 3 pt 10 ps 2 lc rgb '#006400' lw 5
set style line 4 pt 12 ps 2 lc rgb '#2F4F4F' lw 5

set border 3
set xtics nomirror
set ytics nomirror

set output 'detuning_rates.eps'
set multiplot layout 2,1
set lmargin 11
set xrange [-17:15]

f(x) = c/((x-a)**2+b)
# f(x) = c * exp(-(x-a)**2/sqrt(2)/b**2)
a=0.1
b = 1
c = 100000
fit f(x) 'plot_data.dat' u (($1 - 68) * 2):($3/1.875):(sqrt($3)/1.875) via a,b,c
# top plots, single rates
set key t l
set ylabel font "Helvetica,24" "count rate (1/s)"
unset xlabel
set format x ""
set format y '%.0s%c'
plot \
	'plot_data.dat' u (($1 - 68) * 2):($3/1.875):(sqrt($3)/1.875)  title "signal" ls 1, \
	'plot_data.dat' u (($1 - 68) * 2):($4/1.875):(sqrt($4)/1.875) title "idler" ls 2
	# f(x)

# bottom plot, pair rates
unset format x
set format y '%.0s%c'
set xlabel font "Helvetica,24" "two-photon detuning (MHz)"
set ylabel font "Helvetica,24" "pair rate (1/s)"
plot \
	'plot_data.dat' u (($1 - 68) * 2):($2/1.875):(sqrt($2)/1.875) w yerrorbar notitle ls 4

########## plot of the efficiencies and coherence times
unset multiplot
set terminal postscript eps color enhanced solid font "Helvetica,22" size 15cm,8.5cm
set colors classic
set output 'detuning_eff.eps'
# set multiplot layout 2,1
set xrange [-17:15]
set ytics 2

# top plots, efficiencies
set key t l
set ylabel font "Helvetica,24" "efficiency (%)"
# unset xlabel
# set format x ""
plot \
	'plot_data.dat' u (($1 - 68) * 2):(($2/$3)*100) : ($2/$3**2 * sqrt($2*(1+$2/$3)) *100) w yerrorbar title "signal"  ls 1, \
	'plot_data.dat' u (($1 - 68) * 2):(($2/$4)*100) : ($2/$4**2 * sqrt($2*(1+$2/$4)) *100) w yerrorbar title "idler" ls 2, \
    'eff_s.dat' u 1:($2 * 100) w lines notitle ls 1, \
    'eff_i.dat' u 1:($2 * 100) w lines notitle ls 2

# set terminal postscript eps color enhanced solid font "Helvetica,22" size 15cm,8.5cm
# set colors classic
# set output 'detuning_eff.eps'

# unset format x
# set ytics 1
# set xlabel font "Helvetica,24" "two-photon detuning (MHz)"
# set ylabel font "Helvetica,24" "coherence time {/Symbol t} (ns)"
# plot \
# 	'plot_data_tau.dat' u (($1 - 68) * 2):6:7 notitle ls 4