set terminal postscript eps color enhanced solid font 'Helvetica, 16pt' size 8.6cm, 5.5cm  dl .5
	set style line 1 lt 1 lc rgb "black" lw 2.5
	set style line 3 lt 1 pt 6 ps 1 lc rgb "#034E7B" lw 3
	set style line 4 lt 1 pt 6 ps 1 lc rgb "#f7a789" lw 3
	set style line 5 lt 1 pt 6 ps 1 lc rgb "#e36a53" lw 3
	set style line 6 lt 1 pt 6 ps 1 lc rgb '#b40426' lw 3
	set style line 7 lt 1 pt 5 ps 1.1 lc rgb "orange" lw 3
	set style line 10 lt 1 lc 0 lw 0.3
	set style line 11 lt 1 pt 6 ps 1 lc rgb "#034E7B" lw 3
	set style line 12 lt 1 pt 4 ps 1.1 lc rgb "orange" lw 6
	set style line 13 lt 1 pt 6 ps .7 lc rgb "#b40426" lw 3
	set style line 14 lt 1 pt 6 ps 0 lc rgb "#034E7B" lw 3
	set output 'combined_snspd.eps'
	unset bars
	set multiplot layout 3,1 margins 0.12,0.995,.14,.98 spacing 0,0.020
	set autoscale x2fix
	set xtics nomirror
	set ytics nomirror
	set xtics 100
	set ytics 2
	unset key
	fname810 = '810nm.txt'
	fname1550 = '1550nm.txt'
	fname550='550nm.txt'

	stats fname810 u 1:2 name "state"
	k1 = state_pos_max_y
	stats "1550nm.txt" u 1:2 nooutput
	k2 =  STATS_pos_max_y
	stats "550nm.txt" u 1:2 nooutput
	k3 = STATS_pos_max_y
	N1_min =  0.111310831749
	N2_min = 0.754643865846
	N3_min = 0.3073

	set yrange [0 : 7.9]
	set xrange [-150: 150]
	set label 1 "(b)" at -130,7.1
	set ytics 2.5
	set xtics  ('' -100, '' -50,' ' 0 ,'' 50,' ' 100 )
	plot fname810   every ::6::100 u  (($3)-k1):($4/N1_min)/1000:(sqrt($4)/N1_min)/1000  ls 6 pt 6 ps 1 lw 3 with yerrorbars,\
	     fname810  every ::2::1000  u (($1)-k1):(($2)/N1_min)/1000   ls 3 with line,\

	set yrange [0 : 4.7]
	set xrange [-150: 150]
	set label 1 "(c)" at -130,4.2
	set ytics (0, 2,4)
	unset arrow
	plot fname1550  every ::6::100 u (($3)-k2):($4/N2_min)/1000:(sqrt($4)/N2_min)/1000  ls 6 pt 6 ps 1 lw 3 with yerrorbars ,\
		 fname1550  every ::6::990 u (($1)-k2):(($2)/N2_min)/1000   ls 3 with line,\
	 

	set ylabel ' Normalized Cross-correlation (10^3) ' offset 0.4,5.5
	set xlabel ' {/Symbol D} t (ps) ' offset 0,.4
	set yrange [0:5.7]
	set xtics (-100,-50,0,50,100)
	set ytics 2
	set label 1 "(d)" at -130,5.1
	plot fname550   every ::6::95 u (($3)-k3):($4/N3_min)/1000:(sqrt($4)/N3_min)/1000  ls 6 pt 6 ps 1 lw 3 with yerrorbars ,\
		 fname550   every ::6::990 u (($1)-k3):(($2)/N3_min)/1000   ls 3 with line,\