set terminal postscript eps color enhanced font "Helvetica,22" size 12.6cm,8cm
set output 'hom_dip_150ns.eps'

tau_s=26.18 #single atom decay time
tau_f=13.61 #fwm decay time

Pn(x)=0.5 #non-interfering part
Pint(x)=( (x>=0)? (2*tau_s*tau_f/(tau_s+tau_f)**2)*exp(-x/tau_s) : (2*tau_s*tau_f/(tau_s+tau_f)**2)*exp(x/tau_f) ) #interfering part

set samples 1000
set border 3
set bars 0
set xtics 10 nomirror
set ytics 0.2 nomirror
#set xlabel "Excitation pulse delay (ns)"
set xlabel "Excitation time delay {/Symbol D}T (ns)"
set ylabel "Normalized coincidence probability"

plot [-55:55][0:1.1] 'hom_dip_150ns_window.dat' u ($2):(1-$4):(1.5):5 w xyerrorbars lc rgb "red" lt 1 pt 6 ps 1.5 lw 2 notitle,\
1-Pint(x)/Pn(x) w l lc rgb "blue" lt 1 lw 3 notitle
#'hom_dip_theory_150ns_window.dat' u 1:2 w l lc rgb "blue" lt 1 notitle
#'hom_dip_nocorrection.dat' u 1:2:(1):3 w xyerrorbars lc rgb "blue" lt 1 pt 4 notitle,\