set terminal postscript eps enhanced color size 8.6cm, 5cm dl 3\
font "Helvetica, 16pt"
set output 's_vs_bin.eps'

set style line 1 lt 1 lc rgb "black" lw 2.5

set style line 11 lt 1 pt 6 ps .2 lc rgb "#034E7B" lw 3
set style line 12 lt 1 pt 4 ps 1.1 lc rgb "orange" lw 3
set style line 13 lt 1 pt 6 ps .7 lc rgb "#b40426" lw 3

set style line 3 lt 1 pt 6 ps .2 lc rgb "#034E7B" lw 3
# set style line 3 lt 1 pt 6 ps .2 lc rgb "#99baff" lw 3
set style line 4 lt 1 pt 6 ps .2 lc rgb "#f7a789" lw 3
set style line 5 lt 1 pt 6 ps .2 lc rgb "#e36a53" lw 3
set style line 6 lt 1 pt 6 ps .2 lc rgb '#b40426' lw 3
set style line 7 lt 1 pt 5 ps 1.1 lc rgb "orange" lw 3


unset bars
set xlabel 'bin width {/Symbol t} ({/Symbol m}s)'
set ylabel 'residual violation (S - 2)*10^{-3} ' #offset 2,0
set tics front
set ytics 5
set xrange [0:30]
# set yrange [1.995:2.013]
set yrange [-5:17]
set xtics 10
# set key  bottom left samplen .5 spacing 1
# set key on outside top horizontal left
unset key
set label 1 'experiment' at screen .6, .55
# set label 2 "numerical\nsimulation" at screen .7, .85
set label 2 "simulation" at screen .7, .85
set arrow 1 from .170,-5 to .170,0 nohead lc rgb 'red'

# fname = 'plot_data_20171127_203657_loophole_measuremen.dat'
fname = '20171214_171844_loophole_measurement_20k.dat'
fname_th = '20171214_171844_loophole_measurement_20k_th.dat'
plot 0 w filledcurves x1 ls 1 lc rgb "#e0e0e0" notitle ,\
     0 w l ls 1 notitle,\
     fname u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \
     fname_th u ($1 / 1e3):((abs($2)-2) * 1000) with line ls 12 title "simulation",\
     '20171214_171844_loophole_measurement_20k_0_to_0_1.dat' u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \
     '20171214_171844_loophole_measurement_20k_0_2_to_0_3.dat' u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \
     '20171214_171844_loophole_measurement_20k_0_3_to_0_4.dat' u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \
     '20171214_171844_loophole_measurement_20k_0_4_to_1.dat' u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment"#, \

# unset arrow 1

set output 'test.eps'
unset xrange
set xtics auto
unset yrange
unset label 1
unset label 2
set arrow from .170,-3 to .170,2 nohead lc rgb 'red'
plot \
     '20171214_171844_loophole_measurement_20k_0_to_0_1.dat' u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \
     '20171214_171844_loophole_measurement_20k_0_1_to_0_2.dat' u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \
     '20171214_171844_loophole_measurement_20k_0_2_to_0_3.dat' u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \
     '20171214_171844_loophole_measurement_20k_0_3_to_0_4.dat' u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \
     '20171214_171844_loophole_measurement_20k_0_4_to_1.dat' u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \


set output 'S_vs_bin_extended.eps'
set xlabel 'bin width {/Symbol t} ({/Symbol m}s)'
set ylabel 'residual violation (S - 2)*10^{-3} ' #offset 2,0
set tics front
# set ytics 5
unset yrange
set xrange [0:3000]
set yrange [-900:50]
set xtics auto
set ytics auto
set ytics 200
set key  bottom right samplen .5 spacing 1
# set label 1 'experiment' at screen .6, .55
# set label 2 "numerical\nsimulation" at screen .7, .85

fname = '20171214_171844_loophole_measurement_20k_long range.dat'
fname_th = '20171214_171844_loophole_measurement_20k_th_extended.dat'
plot 0 w filledcurves x1 ls 1 lc rgb "#e0e0e0" notitle ,\
     0 w l ls 1 notitle,\
     fname u ($1 / 1e3):((abs($2)-2) * 1000):3 with yerrorbars ls 11 title "experiment", \
     fname_th u ($1 / 1e3):((abs($2)-2) * 1000) with line ls 12 title "simulation"#,\



set output 'rnd_vs_bin.eps'
# set label 1 "numerical\nsimulation" at screen .3, .88
set label 1 "simulation" at screen .24, .88
set label 2 'n{/ZapfDingbats \325}{/Symbol \245}' at screen .24, .71
set label 10 'n=10^{10}' at screen .265, .55
set label 9 'n=10^{9}' at screen .30, .42
set label 8 'n=10^{8}' at screen .34, .25

set xlabel 'bin width {/Symbol t} ({/Symbol m}s)'
set ylabel 'rate of random bits (bits/s)' offset 2,0
set ytics 500
set xrange [0:30]
set yrange [0:2000]
set xtics 10
# set key  top right samplen .5 spacing 1
unset key
# set key on outside top horizontal left
# fname = 'plot_data_20171127_203657_loophole_measuremen.dat'
fname = '20171214_171844_loophole_measurement_20k_rnd.dat'
plot for [i=3:6] fname u 1:i ls i notitle,\
    fname u 1:7 with line ls 7