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

unset key
unset bars
set xlabel 'Area (nVs)' enhanced
set ylabel 'Number of traces' offset 2,0
set xrange [0:40]
set yrange [4:6500]
# set yrange [5:6000]
set logscale y 10
# set ytics 1000

set label 10 'n=1' font 'Helvetica,24' center at 11.5, 10
set label 11 'n=2' font 'Helvetica,24' center at 22.3, 10
set label 12 'n=3' font 'Helvetica,24' center at 32.7, 10
# set label 13 'n=3' font 'Helvetica,24' center at 2.35, 300
# set label 14 'n=4' font 'Helvetica,24' center at 3.05, 200
# set label 15 'n{/Symbol \263}5' font 'Helvetica,24' center at 4, 100

# plt.text(0.6, 600, 'n=0', ha='center', weight='bold')
# plt.text(1.1, 450, 'n=1', ha='center', weight='bold')
# plt.text(1.7, 250, 'n=2', ha='center', weight='bold')
# plt.text(2.4, 100, 'n=3', ha='center', weight='bold')
# plt.text(3, 22, 'n=4', ha='center', weight='bold')

# plot with errorbars
# plot 'area_histo_with_fit.dat' u 1:4 w lines lc 'blue' lw 4 notitle, \
    '' u 1:5 w lines lc 'red' lw 4 notitle, \
    '' u 1:6 w lines lc 'green' lw 4 notitle, \
    '' u 1:($4 +$5 + $6) w lines lc 'orange' lw 4 notitle, \
    '' u 1:2:3 w yerrorbar pt 6 lc 'black' lw 2 ps .7 notitle#, \
    # using 2:(700) with impulses lc rgb 'grey' lw 1.5 #,\
    # 'threshold.dat' using ($2 - .7 ):(650) w label
    # for [cf=2:6] 'fit.dat' u 1:(column(cf)) w lines lc 'black' lw 1 notitle

# plot without errorbars
plot 'area_histo_with_fit.dat' u 1:4 w lines lc 'blue' lw 4 notitle, \
    '' u 1:5 w lines lc 'red' lw 4 notitle, \
    '' u 1:6 w lines lc 'green' lw 4 notitle, \
    '' u 1:($4 +$5 + $6) w lines lc 'orange' lw 4 notitle, \
    '' u 1:2 pt 6 lc 'black' lw 2 ps .7 notitle#, \
    # using 2:(700) with impulses lc rgb 'grey' lw 1.5 #,\
    # 'threshold.dat' using ($2 - .7 ):(650) w label
    # for [cf=2:6] 'fit.dat' u 1:(column(cf)) w lines lc 'black' lw 1 notitle