set terminal postscript eps color solid enhanced font "Helvetica,16"
set output 'plot_g2.eps'

set isosamples 20,20
set samples 1000,1000

set multiplot


# Dark counts of the APDs used to detect the photons
d1=175;d2=520;

# Pairs to singles efficiency of the source.
eff1=0.184;
eff2=0.14;

# ratio of total pairs to peak coincidences in 1 ns bin 
t=8.4;

# Expected g(2) as a function of pair rate x
g(x)=((((x/eff1)+d1)*((x/eff2)+d2)*1e-9)+(x/t))/(((x/eff1)+d1)*((x/eff2)+d2)*1e-9);

set size 1,1
set origin 0,0

set xlabel "Pair Rate"
set ylabel "g^{(2)}(0)"

plot [0:4500] 'rate_g2_singles.dat' u 1:2:5 w yerrorbars notitle ps 1 pt 5, g(x) lc rgb "blue" notitle

set size 0.6,0.6
set origin 0.3,0.3

set xtics 20 font "Helvetica,12"
set ytics 2000 font "Helvetica,12"

unset xlabel
unset ylabel

plot [0:200] 'rate_g2_singles.dat' u 1:2:5 w yerrorbars notitle ps 1 pt 5, g(x) lc rgb "blue" notitle