src_cross="out_cross.dat"
src_auto_762="g2_auto_762.dat"
src_auto_795="g2_auto_795.dat"

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

set multiplot

f(x)=a*(1+b*(exp(-abs(x-x0)/c)))
a=54.3;b=0.96;c=9.9;x0=0;
#fit [-125:125] f(x) src_auto_795 u ($1-325):2:(sqrt($2)) via a,b,c,x0

set size 0.985,0.39
set origin 0.02,0

set ylabel font "Helvetica,18" "G_{II}^{(2)} ({/Symbol t})"
set xlabel font "Helvetica,18" "{/Symbol t} (ns)"

set ytics 40 nomirror font "Helvetica,16"
set y2range [0.77:2.59]
set y2tics 0.5 font "Helvetica,16"
set y2label font "Helvetica,18" "g_{II}^{(2)} ({/Symbol t})"

plot [-75:75][40:135] src_auto_795 u ($1-328):2:(sqrt($2)) w yerrorbars notitle pt 7 ps 0.5, f(x) w l lc rgb "blue" lw 2 lt 1 notitle

set xlabel ""
set format x ""

 
a=105.5;b=0.94;c=17.8;x0=0;

#fit [-125:125] f(x) src_auto_762 u ($1-325):2 via a,b,c,x0

set size 0.985,0.33
set origin 0.02,0.355

set ylabel font "Helvetica,18" "G_{SS}^{(2)} ({/Symbol t})"
set ytics 50 nomirror font "Helvetica,16"
set y2range [0.82:2.47]
set y2tics 0.5 font "Helvetica,16"
set y2label font "Helvetica,18" "g_{SS}^{(2)} ({/Symbol t})

plot [-75:75][80:240] src_auto_762 u ($1-328):2:(sqrt($2)) w yerrorbars notitle pt 7 ps 0.5,f(x) w l lc rgb "blue" lw 2 lt 1 notitle

set size 1,0.33
set origin 0.005,0.65
set format x ""
set ylabel font "Helvetica,18" "G_{SI}^{(2)} ({/Symbol t})"
set ytics 600 nomirror font "Helvetica,16"
set y2range [0:6.3]
set y2tics 1.5 font "Helvetica,16"
set y2label font "Helvetica,18" "g_{SI}^{(2)} ({/Symbol t})  ({/Symbol \664} 10^3)"

set label font "Helvetica,14" "g_{II}^{(2)}(0) = \n\n 2.03 {/Symbol \661} 0.08" at first -60, first -3400
set label font "Helvetica,14" "g_{SS}^{(2)}(0) = \n\n 2.06 {/Symbol \661} 0.06" at first -60, first -750
set label font "Helvetica,14" "g_{SI}^{(2)}(0) = \n\n 5800 {/Symbol \661} 76 " at first -60, first 1800
#set label font "Helvetica,8" "{/Symbol t}_0 = 6.7(2) ns" at first 45, first 700
#set label font "Helvetica,8" "{/Symbol t}_0 = 21(1) ns" at first 45, first -950
#set label font "Helvetica,8" "{/Symbol t}_0 = 10.8(8) ns" at first 45, first -2600
set label font "Helvetica,16" "(c)" at first 50, first -3400
set label font "Helvetica,16" "(b)" at first 50, first -750
set label font "Helvetica,16" "(a)" at first 50, first 1800
#set label font "Helvetica,8" "O.D = 25" at first -65, first 1100
#set label font "Helvetica,8" '{/Symbol D} \~ 50 MHz' at first -65, first 750
#set label font "Helvetica,8" "O.D = 23.4" at first -65, first -650
#set label font "Helvetica,8" '{/Symbol D} \~ 0' at first -65, first -900
#set label font "Helvetica,8" "O.D = 23.4" at first -65, first -2300
#set label font "Helvetica,8" '{/Symbol D} \~ 0' at first -65, first -2550
f(x)=B+A*exp(-(x-x0)/C)
B=0;A=2390;C=6.7;x0=0;
#fit [0:100] f(x) src_cross u ($1-175):2 via A,B,C,x0
plot [-75:75][0:2400] src_cross u ($1-179):2:3 w yerrorbars notitle pt 7 ps 0.5, f(x) w l lc rgb "blue" lw 2 lt 1 notitle, 'coincidence_left.dat' u ($1-3):($2*2400) w l lc rgb "#008000" lw 2 lt 3 notitle,'coincidence_right.dat' u ($1-3):($2*2400) w l lc rgb "#008000" lw 2 lt 3 notitle