src_cross="out_cross.dat"
src_seedpower="030512c_g2_seedpower.dat"
src_auto_795="g2_auto_795.dat"

set terminal postscript enhanced eps color size 10cm,10cm font "Helvetica,16pt"
set output "g2_combined.eps"
set isosamples 40,40
set samples 1000,1000
set border 3
set nobar
set multiplot
set bar 0

set size 1,0.5
set origin 0.0,0.48


set ylabel font "Helvetica,16" "g^{(2)}"
set xlabel font "Helvetica,16" "{/Symbol D}t (ns)"

set ytics 0.5 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})"
#fit [-75:75] f(x) src_auto_795 u ($1-325):($2/55):(sqrt($2)/55) via a,b,c,x0
set xtics nomirror
set ytics nomirror
f(x)=a*(1+b*(exp(-abs(x-x0)/c)))
a=1;b=1;c=9.9;x0=0
fit [-75:75] f(x) src_auto_795 u ($1-328):($2/55.5):(sqrt($2)/55.5) via a,b,c
plot [-75:75] src_auto_795 u ($1-328):($2/55.5):(sqrt($2)/55.5) w yerrorbars notitle lt 1 pt 7 ps 1, f(x) w l lc rgb "blue" lw 2 lt 1 notitle


 
#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
unset label
set size 1,0.48
set origin 0.01,0.01

set xlabel font "Helvetica,16" "seed power (mW)"
set ylabel font "Helvetica,16" "g^{(2)}"
#set format x ""
#set label font "Helvetica,18" "G_{SS}^{(2)} ({/Symbol t}) at first -80, first 75 rotate #by 90"
#set label font "Helvetica,18" "g_{SS}^{(2)} ({/Symbol t}) at first 75, first 75 rotate by #90"
set ytics 0.5 nomirror font "Helvetica,16"

#set y2label font "Helvetica,18" "g_{SS}^{(2)} ({/Symbol t})
set xtics nomirror
set ytics nomirror
#set style line 1 lc rgb '#8b1a0e' pt 1 ps 1 lt 1 lw 2
plot src_seedpower u 1:2:3 w yerrorbars notitle pt 7 ps 1