#!/usr/bin/gnuplot
infile = 'sum/g2sum_5ns.dat'
#infile = 'sum/g2sum_10ns.dat'
infile2 = 'sum/g2sum_5ns_long.dat'
infile3 = 'sum/g2sum_5ns_long_bg.dat'
set terminal postscript eps color enhanced solid font "Helvetica,23" size 20cm,9cm
#set term pdf
#set term pngcairo
set output "g2_sup.eps"

set border linewidth 1.5
set samples 5000
set multiplot 
unset bar
unset key
set title ' '

#####  ####
set lmargin at screen 0.10
set rmargin at screen 0.37

delay = 482.5
gamma0 = 2*pi*7 *1e-3
lam = 0.08077
yoffset = 0

set tics front
set xlabel "Time {/Symbol t} (ns)"
set ylabel 'g^{(2)}({/Symbol t})' offset 1,0
set pointsize 1
set xrange [-231:231]
set yrange [0.917:1.038]
set ytics 0.92,0.04,1 nomirror
set xtics 100 nomirror

set label 1 sprintf("{/:Bold a}") at graph -0.1,1.1 font "Helvetica, 24"

plot infile u ($1-delay):($6):5 w ye pt 7 lw 2 lc rgbcolor "#00c000"

#################################
#######

set lmargin at screen 0.40
set rmargin at screen 0.67

y0=1
amp=0.02
tau=800
x0=485

exp_decay(x) = (  x > 486 ) ? amp*exp(-(x-x0)/tau ) + y0 : 1/0  

FIT_LIMIT=1.e-10
fit [685:4500] exp_decay(x) infile2 u 1:6:5 via amp,tau


set xlabel "Time {/Symbol t} (ns)"
#set ylabel 'g^{(2)}({/Symbol t})' offset 2,0
unset ylabel
set pointsize 1
set xrange [-231:4500]
set yrange [0.917:1.038]
#set ytics 0.92,0.04,1 nomirror
set ytics format " "
set xtics 1000 nomirror

set label 1 sprintf("{/:Bold b}") at graph -0.1,1.1 font "Helvetica, 24"

plot infile2 u ($1-delay):($6):5 w l lc rgbcolor "#00c000" , exp_decay(x+x0) t 'exp decay' lc rgb "black" lw 4

#####
set lmargin at screen 0.70
set rmargin at screen 0.97
set label 1 sprintf("{/:Bold c}") at graph -0.1,1.1 font "Helvetica, 24"
plot infile3 u ($1-delay):($6):5 w l lc rgbcolor "#00c000"

unset multiplot
unset output
pause -1