set terminal postscript eps color enhanced solid font 'Helvetica,22' size 15cm,17cm
# set colors classic
set nobars
set loadpath 'data'

# set of labels used throughout the script
lab1 = ' 30 {/Symbol m}W'
lab2 = '130 {/Symbol m}W'
lab3 = '290 {/Symbol m}W'
lab4 = '420 {/Symbol m}W'
lab5 = '580 {/Symbol m}W'
lab6 = '640 {/Symbol m}W'

# general label offset
lab_off = 15.8

# overall style commands
set border 3
set xtics nomirror
set ytics nomirror

# first plot, rates
set output 'power_rates.eps'
set multiplot layout 2,1
set lmargin 10

set yrange [0:34e4]
set xrange [0:16]
th_ext = 14.8
set xtics 5
set key off

# top plots, signal rates
set ylabel font 'Helvetica,24' 'signal rate (1/s)'
# set key on outside top horizontal left
unset xlabel
set format x ''
set format y '%.0s%c'
set ytics 100000

load 'mypal.pal'
set label 'P_{780}' font 'Helvetica,24' center at first lab_off, graph 1.01
plot \
	'fit_signal.dat' u 2:($2<th_ext?$3:0/0) index 0 w lines ls 3 notitle,\
	'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):6:7 w yerror title lab1 ls 3,\
	'fit_signal.dat' u 2:($2<th_ext?$3:0/0) index 1 w lines ls 4 notitle,\
	'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):6:7 w yerror title lab2 ls 4,\
	'fit_signal.dat' u 2:($2<th_ext?$3:0/0) index 2 w lines ls 5 notitle,\
	'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):6:7 w yerror title lab3 ls 5,\
	'fit_signal.dat' u 2:($2<th_ext?$3:0/0) index 3 w lines ls 6 notitle,\
	'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):6:7 w yerror title lab4 ls 6,\
	'fit_signal.dat' u 2:($2<th_ext?$3:0/0) index 4 w lines ls 7 notitle,\
	'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):6:7 w yerror title lab5 ls 7,\
	'fit_signal.dat' u 2:($2<th_ext?$3:0/0) index 5 w lines ls 8 notitle,\
	'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):6:7 w yerror title lab6 ls 8,\
	 'plot_data_0.03.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab1, $1)) with labels notitle tc ls 11,\
	 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab2, $1)) with labels notitle tc ls 12,\
	 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab3, $1)) with labels notitle tc ls 13,\
	 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab4, $1)) with labels notitle tc ls 14,\
	 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab5, $1)) with labels notitle tc ls 15,\
	 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):($6+1e4):(sprintf(lab6, $1)) with labels notitle tc ls 16
unset format x
unset label


# bottom plot, idler rates
set xlabel font 'Helvetica,24' 'P_{776} (mW)'
set ylabel font 'Helvetica,24' 'idler rate (1/s)'
set key off
set yrange [0:47e4]

# set label 'P_{780}' right at first lab_off, graph 1.05
plot \
	'fit_idler.dat' u 2:($2<th_ext?$3:0/0) index 0 w lines ls 11 notitle,\
	'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):8:9 w yerror notitle ls 11,\
	'fit_idler.dat' u 2:($2<th_ext?$3:0/0) index 1 w lines ls 12 notitle,\
	'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):8:9 w yerror notitle ls 12,\
	'fit_idler.dat' u 2:($2<th_ext?$3:0/0) index 2 w lines ls 13 notitle,\
	'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):8:9 w yerror notitle ls 13,\
	'fit_idler.dat' u 2:($2<th_ext?$3:0/0) index 3 w lines ls 14 notitle,\
	'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):8:9 w yerror notitle ls 14,\
	'fit_idler.dat' u 2:($2<th_ext?$3:0/0) index 4 w lines ls 15 notitle,\
	'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):8:9 w yerror notitle ls 15,\
	'fit_idler.dat' u 2:($2<th_ext?$3:0/0) index 5 w lines ls 16 notitle,\
	'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):8:9 w yerror notitle ls 16,\
	 'plot_data_0.03.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab1, $1)) with labels notitle tc ls 11,\
	 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab2, $1)) with labels notitle tc ls 12,\
	 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab3, $1)) with labels notitle tc ls 13,\
	 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab4, $1)) with labels notitle tc ls 14,\
	 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab5, $1)) with labels notitle tc ls 15,\
	 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):($8 + 1e4):(sprintf(lab6, $1)) with labels notitle tc ls 16
unset label
unset multiplot



# plot of the efficiencies
set output 'power_eff.eps'
set multiplot layout 2,1
# set xrange [0:16]
unset yrange
unset format y
set ytics autofreq

# top plots, signal efficiencies
set ylabel font 'Helvetica,24' 'signal heralded efficiency (%)'
unset xlabel
set format x ''
# set key on outside top horizontal left
unset key
set label 'P_{780}' font 'Helvetica,24' center at first lab_off, graph 1.01
plot \
	'fit_eff_s.dat' u 2:($2<th_ext?($3 * 100):0/0) index 0 w lines ls 3 notitle,\
	'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):($4/$6*100):((sqrt(($5/$6)**2 + ($7*$4/($6**2))**2))*100) w yerror title lab1 ls 3,\
	'fit_eff_s.dat' u 2:($2<th_ext?($3 * 100):0/0) index 1 w lines ls 4 notitle,\
	'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):($4/$6*100):((sqrt(($5/$6)**2 + ($7*$4/($6**2))**2))*100) w yerror title lab2 ls 4,\
	'fit_eff_s.dat' u 2:($2<th_ext?($3 * 100):0/0) index 2 w lines ls 5 notitle,\
	'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):($4/$6*100):((sqrt(($5/$6)**2 + ($7*$4/($6**2))**2))*100) w yerror title lab3 ls 5,\
	'fit_eff_s.dat' u 2:($2<th_ext?($3 * 100):0/0) index 3 w lines ls 6 notitle,\
	'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):($4/$6*100):((sqrt(($5/$6)**2 + ($7*$4/($6**2))**2))*100) w yerror title lab4 ls 6,\
	'fit_eff_s.dat' u 2:($2<th_ext?($3 * 100):0/0) index 4 w lines ls 7 notitle,\
	'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):($4/$6*100):((sqrt(($5/$6)**2 + ($7*$4/($6**2))**2))*100) w yerror title lab5 ls 7,\
	'fit_eff_s.dat' u 2:($2<th_ext?($3 * 100):0/0) index 5 w lines ls 8 notitle,\
	'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):($4/$6*100):((sqrt(($5/$6)**2 + ($7*$4/($6**2))**2))*100) w yerror title lab6 ls 8,\
	 'plot_data_0.03.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab1, $1)) with labels notitle tc ls 3,\
	 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab2, $1)) with labels notitle tc ls 4,\
	 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab3, $1)) with labels notitle tc ls 5,\
	 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab4, $1)) with labels notitle tc ls 6,\
	 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab5, $1)) with labels notitle tc ls 7,\
	 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab6, $1)) with labels notitle tc ls 8
unset format x
unset label

# bottom plot, idler efficiency
set xlabel font 'Helvetica,24' 'P_{776} (mW)'
set ylabel font 'Helvetica,24' 'idler heralded efficiency (%)'
set key off
plot \
	'fit_eff_i.dat' u 2:($2<th_ext?($3 * 100):0/0) index 0 w lines ls 3 notitle,\
	'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):($4/$8*100):((sqrt(($5/$8)**2 + ($9*$4/($8**2))**2))*100) w yerror title lab1 ls 11,\
	'fit_eff_i.dat' u 2:($2<th_ext?($3 * 100):0/0) index 1 w lines ls 4 notitle,\
	'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):($4/$8*100):((sqrt(($5/$8)**2 + ($9*$4/($8**2))**2))*100) w yerror title lab2 ls 12,\
	'fit_eff_i.dat' u 2:($2<th_ext?($3 * 100):0/0) index 2 w lines ls 5 notitle,\
	'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):($4/$8*100):((sqrt(($5/$8)**2 + ($9*$4/($8**2))**2))*100) w yerror title lab3 ls 13,\
	'fit_eff_i.dat' u 2:($2<th_ext?($3 * 100):0/0) index 3 w lines ls 6 notitle,\
	'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):($4/$8*100):((sqrt(($5/$8)**2 + ($9*$4/($8**2))**2))*100) w yerror title lab4 ls 14,\
	'fit_eff_i.dat' u 2:($2<th_ext?($3 * 100):0/0) index 4 w lines ls 7 notitle,\
	'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):($4/$8*100):((sqrt(($5/$8)**2 + ($9*$4/($8**2))**2))*100) w yerror title lab5 ls 15,\
	'fit_eff_i.dat' u 2:($2<th_ext?($3 * 100):0/0) index 5 w lines ls 8 notitle,\
	'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):($4/$8*100):((sqrt(($5/$8)**2 + ($9*$4/($8**2))**2))*100) w yerror title lab6 ls 16,\
	 'plot_data_0.03.dat' every ::9 u ($0*0+lab_off):($4/$8*100):(sprintf(lab1, $1)) with labels notitle tc ls 3,\
	 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):($4/$8*100):(sprintf(lab2, $1)) with labels notitle tc ls 4,\
	 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):($4/$8*100):(sprintf(lab3, $1)) with labels notitle tc ls 5,\
	 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):($4/$8*100):(sprintf(lab4, $1)) with labels notitle tc ls 6,\
	 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):($4/$8*100):(sprintf(lab5, $1)) with labels notitle tc ls 7,\
	 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):($4/$8*100-.4):(sprintf(lab6, $1)) with labels notitle tc ls 8
unset multiplot


# plot of pair rate and coherence times
set terminal postscript eps color enhanced solid font 'Helvetica,22' size 15cm,8.5cm
set output 'power_coh.eps'
set lmargin 10

unset format x
# set yr [0:22]
set yr [6:22]
set ytics autofreq
set xlabel font 'Helvetica,24' 'P_{776} (mW)'
set ylabel font 'Helvetica,24' 'coherence time (ns)'
set key off
set label 'P_{780}' font 'Helvetica,24' center at first lab_off, graph 1.01
plot \
	'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerrorlines title lab1 ls 3,\
	'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerrorlines title lab2 ls 4,\
	'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerrorlines title lab3 ls 5,\
	'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerrorlines title lab4 ls 6,\
	'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerrorlines title lab5 ls 7,\
	'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerrorlines title lab6 ls 8,\
	 'plot_data_0.03.dat' every ::9 u ($0*0+lab_off):($10-.7):(sprintf(lab1, $1)) with labels notitle tc ls 3,\
	 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab2, $1)) with labels notitle tc ls 4,\
	 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab3, $1)) with labels notitle tc ls 5,\
	 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab4, $1)) with labels notitle tc ls 6,\
	 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab5, $1)) with labels notitle tc ls 7,\
	 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab6, $1)) with labels notitle tc ls 8
unset multiplot

# plot of pair rate
set terminal postscript eps color enhanced solid font 'Helvetica,22' size 15cm,8.5cm
set nobars
set output 'power_pairs.eps'
unset yr
unset label
set lmargin 10
set ylabel font 'Helvetica,24' 'pair rate (1/s)'
set format y '%.0s%c'
set key on outside top horizontal right
load 'mypal.pal'
plot \
	'fit_pairs.dat' u 2:3 index 0 w lines ls 3 notitle,\
	'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab1 ls 3,\
	'fit_pairs.dat' u 2:3 index 1 w lines ls 4 notitle,\
	'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab2 ls 4,\
	'fit_pairs.dat' u 2:3 index 2 w lines ls 5 notitle,\
	'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab3 ls 5,\
	'fit_pairs.dat' u 2:3 index 3 w lines ls 6 notitle,\
	'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab4 ls 6,\
	'fit_pairs.dat' u 2:3 index 4 w lines ls 7 notitle,\
	'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab5 ls 7,\
	'fit_pairs.dat' u 2:3 index 5 w lines ls 8 notitle,\
	'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab6 ls 8

# set xlabel font 'Helvetica,24' 'P_{776} (mW)'