set terminal postscript eps color enhanced font "Helvetica,16pt" size 8.6cm,12cm 
set output "g2_sun_noon2.eps"

# for fit function
N(x)=a+b*(1+(exp(-abs((x-d)/tc)))**2)

# common settings
set xrange [-11.87:-7.87]
set format y "%4g"
set format x "%2.1f"
# set tics scale 2
set tics nomirror
set format y2 "%2.1f"

# bottom figure
set origin 0,0
set size 1,0.44

a=430
b=260
tc=0.32
d=-9.87

fit N(x) 'g2_sun_noon2.dat' using 1:2:3 via a,b,tc,d

set y2r [600/(a+b):1000/(a+b)]
set y2tics 0.1

set xlabel "Timing separation {/Symbol t} [ns]"
set yrange [600:1000]
set ylabel "Two-Photon Coincidences"

plot 'g2_sun_noon2.dat' using 1:2:3 w yerrorbars notitle,\
N(x) w l ls -1 notitle

# d = -9.869 +/- 0.009 ns
# tc = 0.32 +/- 0.04 ns
# a = 430 +/- 20
# b = 260 +/- 20
# g2 peak = 1.37 +/- 0.01
# reduced chisq = 0.995616