set terminal postscript eps color enhanced solid font 'Helvetica, 16pt' size 8.6cm, 6cm  dl .5
set style line 1 lt 1 lc rgb "black" lw 2.5
set style line 2 lt 5 lc rgb "grey" lw 2
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 'selfb_200us_combined.eps'
# the histgram boxes setting
set style fill solid 0.6 
set boxwidth 1
set multiplot
set ylabel 'Probability Distribution ' offset 1.6,0
set xlabel 'Number of Photon Counts (Time interval 200 {/Symbol m}s) '
set ytics 0.5
set xtics 5
fname = 'normal_operation_Led_on_200us_self_blind.txt'
fname1 = 'Led_on_200us_just_b.txt'
set label "(e)" right at graph 0.125, graph 0.96
stats fname using 1:2 name "noblind"
stats fname1 using 1:2 name "blind"
set xrange [-0.6:30]
set yrange [0:1.05]
set key title " "
set key top left Left reverse samplen 2
plot fname  u 1:($2)/noblind_sum_y:(sqrt($2)/noblind_sum_y) ls 6 pt 6 ps 1 lw 1 lc "blue" with boxerrorbars title "Unblinded",\
     fname1  u 1:($2)/blind_sum_y:(sqrt($2)/blind_sum_y) ls 6  lw 1   with boxerrorbars title "Blinded ",\


set size 0.58,0.6
set origin 0.38,0.3
unset ylabel
unset xlabel
unset xtics
unset yrange
unset xrange
unset label
set ylabel ' Detection Probability' offset 3,0
set xlabel 'Time difference {/Symbol D}t (ns)'
set yrange [0:0.14]
set xrange [0:100]
set ytics 0.1  offset 0.9,0
set xtics 20
overall_pulse = 7659
fname = "no_blinding_200us_g2.txt"
fname1 = "just_blinding_200us_g2.txt"
set key title " "
set key at 36,0.14   samplen 2
plot fname  u 1:($2 == 0 ? NaN : $2)/overall_pulse:(sqrt(($2)))/overall_pulse  ls 6 pt 6 ps 1.5 lw 3 lc "blue" with errorbars title "Unblinded",\
     fname1  u 1:($2 == 0 ? NaN : $2)/overall_pulse:(sqrt(($2)))/overall_pulse  ls 6 pt 2 ps 1.5 lw 3 with errorbars title "Blinded",\