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 'Led_on_200us_noblind_self.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'
set label "(d)" right at graph 0.125, graph 0.96
# fname1 ='Led_on_20us_noblind.txt'
stats fname using 1:2 name "on"
# stats fname1 using 1:2 name "on"
set xrange [-0.6:30]
set yrange [0:on_max_y/on_sum_y*1.1]
plot fname  u 1:($2)/on_sum_y:(sqrt($2)/on_sum_y) ls 6 pt 6 ps 1 lw 1 with boxerrorbars title "",\
     # fname1  u 1:($2)/on_sum_y:(sqrt($2)/on_sum_y) ls 2  lw 1 with boxerrorbars title "No attack present ",\


set size 0.48,0.5
set origin 0.5,0.4
set obj 1 rectangle behind from graph 0.,0. to graph 1,1
set obj 1 fillstyle solid 1.0 fillcolor rgbcolor "white"
unset ylabel
unset xlabel
unset key
unset ytics
unset xtics
unset yrange
unset xrange
set xtics 30
set ytics 0.1
set ylabel ' Detection Probability' offset 2,0
set xlabel 'Time difference {/Symbol D}t (ns)'
overall_pulse = 7608
set yrange [0:0.15]
set xrange [0:71]
set ytics 0.1  offset 0.9,0
set format y '10^{%T}'
set ytics add('0' 0)
fname1 = "no_blinding_200us_g2.txt"
set bars small
plot fname1  u 1:($2 == 0 ? NaN : $2)/overall_pulse:(sqrt(($2)))/overall_pulse  ls 1 pt 6 ps 1 with errorbars ,\