set terminal postscript eps color enhanced solid font "Helvetica,20" size 13cm,17cm
set colorsequence classic
set output 'g2.eps'
set multiplot layout 4, 1

infile0='data_g2/75amp.dat'
infile0_s='sim_g2/75amp.dat'
infile1='data_g2/195amp.dat'
infile1_s='sim_g2/195amp.dat'
infile2='data_g2/265amp.dat'
infile2_s='sim_g2/265amp.dat'
infile3='data_g2/350amp.dat'
infile3_s='sim_g2/350amp.dat'
infile4='data_g2/455amp.dat'
infile4_s='sim_g2/455amp.dat'

unset ylabel
unset xlabel
unset key

#
set yrange [0:1.85]
set xrange [-210:210]
set ytics 0.5 nomirror
set xtics 100 nomirror
set format x ""
set lmargin 8
set rmargin 1.2
set tmargin at screen 0.98
set bmargin at screen 0.82
set label 1 at -203, 1.65
set label 1 "(a)" tc lt 7
plot infile0 u ($1-476):($4/147.08):(sqrt($4)/147.08) w yerr lc 1 ps 0.9 pt 7 t '195amp_data',\
infile0_s u ($1-476):($2/147.08) w lines lw 3 lc 3
#
unset label 1
set label 1 at -203, 1.65
set label 1 "(b)" tc lt 7
set label 2 at -200, 0.3
set label 2 "{/Symbol W} = 2{/Symbol p} x 15.8(2) MHz" tc lt 7
set tmargin at screen 0.80
set bmargin at screen 0.63
plot infile1 u ($1-479):($4/298.6):(sqrt($4)/298.6) w yerr lc 1 ps 0.9 pt 7 t '195amp_data',\
infile1_s u ($1-479):($2/298.6) w lines lw 3 lc 3
#
unset label 1
unset label 2
set label 1 at -203, 1.65
set label 1 "(c)" tc lt 7
set label 2 at -200, 0.3
set label 2 "{/Symbol W} = 2{/Symbol p} x 27.8(2) MHz" tc lt 7
set tmargin at screen 0.61
set bmargin at screen 0.44
set ylabel "g^{2} ({/Symbol t})" offset 1,0
plot infile2 u ($1-479):($4/251.58):(sqrt($4)/251.58) w yerr lc 1 ps 0.9 pt 7 t '195amp_data',\
infile2_s u ($1-479):($2/251.58) w lines lw 3 lc 3
#
unset label 1
unset label 2
set label 1 at -203, 1.65
set label 1 "(d)" tc lt 7
set label 2 at -200, 0.3
set label 2 "{/Symbol W} = 2{/Symbol p} x 43.9(2) MHz" tc lt 7
set tmargin at screen 0.42
set bmargin at screen 0.25
unset ylabel
plot infile3 u ($1-481):($4/222):(sqrt($4)/222) w yerr lc 1 ps 0.9 pt 7 t '195amp_data',\
infile3_s u ($1-481):($2/222) w lines lw 3 lc 3
#
unset label 1
unset label 2
set label 1 at -203, 1.65
set label 1 "(e)" tc lt 7
set label 2 at -200, 0.3
set label 2 "{/Symbol W} = 2{/Symbol p} x 64.3(2) MHz" tc lt 7
set format x "%.0f"
set tmargin at screen 0.23
set bmargin at screen 0.06
set xlabel "Time {/Symbol t} (ns)" offset 0,0.4
plot infile4 u ($1-481):($4/257):(sqrt($4)/257) w yerr lc 1 ps 0.9 pt 7 t '195amp_data',\
infile4_s u ($1-481):($2/257) w lines lw 3 lc 3
unset multiplot
pause -1