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' lab1c = ' 30 {/Symbol m}W' lab2c = '130 {/Symbol m}W' lab3c = '290 {/Symbol m}W' lab4c = '420 {/Symbol m}W' lab5c = '580 {/Symbol m}W' lab6c = '640 {/Symbol m}W' tc 'red' # 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:] # set xrange [0:16] set xrange [0:25] 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:3 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:3 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:3 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:3 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:3 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:3 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,\ 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab2, $1)) with labels notitle,\ 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab3, $1)) with labels notitle,\ 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab4, $1)) with labels notitle,\ 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab5, $1)) with labels notitle,\ 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):6:(sprintf(lab6, $1)) with labels notitle 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 label 'P_{780}' right at first lab_off, graph 1.05 plot \ 'fit_idler.dat' u 2:3 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:3 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:3 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:3 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:3 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:3 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,\ 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab2, $1)) with labels notitle,\ 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab3, $1)) with labels notitle,\ 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab4, $1)) with labels notitle,\ 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab5, $1)) with labels notitle,\ 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):8:(sprintf(lab6, $1)) with labels notitle unset label unset multiplot # plot of the efficiencies set output 'power_eff.eps' set multiplot layout 2,1 set xrange [0:16] 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:($3 * 100) 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:($3 * 100) 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:($3 * 100) 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:($3 * 100) 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:($3 * 100) 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:($3 * 100) 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,\ 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab2, $1)) with labels notitle,\ 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab3, $1)) with labels notitle,\ 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab4, $1)) with labels notitle,\ 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab5, $1)) with labels notitle,\ 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):($4/$6*100):(sprintf(lab6, $1)) with labels notitle 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:($3 * 100) 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:($3 * 100) 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:($3 * 100) 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:($3 * 100) 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:($3 * 100) 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:($3 * 100) 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,\ 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):($4/$8*100):(sprintf(lab2, $1)) with labels notitle,\ 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):($4/$8*100):(sprintf(lab3, $1)) with labels notitle,\ 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):($4/$8*100):(sprintf(lab4, $1)) with labels notitle,\ 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):($4/$8*100):(sprintf(lab5, $1)) with labels notitle,\ 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):($4/$8*100-.4):(sprintf(lab6, $1)) with labels notitle 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 multiplot layout 2,1 set xrange [0:16] set lmargin 10 # top plots, pair rate # set ylabel font 'Helvetica,24' 'pair rate (1/s)' # unset xlabel # set format x '' # set format y '%.0s%c' # set key on outside top horizontal right # # load 'greens.pal' # plot \ # 'plot_data_0.03.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab1 ls 3,\ # 'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab2 ls 4,\ # 'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab3 ls 5,\ # 'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab4 ls 6,\ # 'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab5 ls 7,\ # 'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):4:5 w yerror title lab6 ls 8 unset format x set yr [0: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 yerror title lab1 ls 3,\ 'plot_data_0.06.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerror title lab2 ls 4,\ 'plot_data_0.09.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerror title lab3 ls 5,\ 'plot_data_0.12.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerror title lab4 ls 6,\ 'plot_data_0.15.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerror title lab5 ls 7,\ 'plot_data_0.18.dat' u (250*($1-19)/(4750*0.5)):10:11 w yerror title lab6 ls 8,\ 'plot_data_0.03.dat' every ::9 u ($0*0+lab_off):($10-.7):(sprintf(lab1, $1)) with labels notitle,\ 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab2, $1)) with labels notitle,\ 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab3, $1)) with labels notitle,\ 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab4, $1)) with labels notitle,\ 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab5, $1)) with labels notitle,\ 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab6, $1)) with labels notitle unset multiplot # plot of pair rate set terminal postscript eps color enhanced solid font 'Helvetica,22' size 15cm,8.5cm # set colors classic set nobars set output 'power_pairs.eps' # set multiplot layout 2,1 set xrange [0:16] unset yr unset label set lmargin 10 # top plots, pair rate set ylabel font 'Helvetica,24' 'pair rate (1/s)' # unset xlabel # set format x '' 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 # unset format x # set yr [0:22] # set ytics autofreq set xlabel font 'Helvetica,24' 'P_{776} (mW)' set output 'power_times.eps' set ylabel font 'Helvetica,24' 'coherence time (ns)' set key off # load 'greys.pal' 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,\ 'plot_data_0.06.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab2, $1)) with labels notitle,\ 'plot_data_0.09.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab3, $1)) with labels notitle,\ 'plot_data_0.12.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab4, $1)) with labels notitle,\ 'plot_data_0.15.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab5, $1)) with labels notitle,\ 'plot_data_0.18.dat' every ::9 u ($0*0+lab_off):10:(sprintf(lab6, $1)) with labels notitle