reset

set term postscript enhanced eps color size 8.6cm,6cm font "Helvetica,16pt"
set output 'distribution.eps'

f(x) = a*exp((-(x-b)**2)/(c**2))
g(x) = d*exp((-(x-e)**2)/(f**2))

a = 170000
b = 10
c = 700
d = 400000
e = 10
f = 300

set noytics
set noxtics


#set style arrow 3 nohead
set style fill transparent solid 0.8 border lc rgb 'grey'

set xrange[-5000:5000]
#plot f(x), 'distribution_ele.txt' using 1:2
fit f(x) 'histogram_30.txt' using 2:1 via a,b,c
fit g(x) 'histogram_0.txt' using 2:1 via d,e,f

set multiplot 
# layout 3, 1
# set lmargin at screen 0.05
# set tmargin at screen 0.95
# set bmargin at screen 0.7

#plot1
set size 1,0.35
set origin 0,0.66
unset xlabel
unset key 
set style arrow 8 heads size screen 0.008,90 lc rgb 'red' lw 1 front
set arrow from 0,40000 to c,40000 as 8
set yrange[0:400000]
set ylabel "{p_t}(x)" # font "Helvetica, 26" offset -1,0
set label "{{/Symbol s}_t}" at 1200,100000 # font "Helvetica, 26"
set label 2 "(a)"  at 3500,350000 font "Helvetica,18"

plot 'histogram_30.txt' using ($2)-b:1 with boxes lc rgb 'grey', f(x+b) w l lw 1 lc rgb 'red'

#plot f(x+b) w l lw 1 lc rgb 'red'


#plot2
unset key
unset arrow
#set tmargin at screen 0.665
#set bmargin at screen 0.414
set size 1,0.35
set origin 0,0.38
set style arrow 8 heads size screen 0.008,90 lc rgb 'red' lw 1 front
set arrow front from 0,50000 to f,50000 as 8
unset label
set label "{{/Symbol s}_e}" at 700,100000 # font "Helvetica, 26"
set label 2 "(b)"  at 3500,420000 font "Helvetica,18"
set yrange[0:600000]
unset ylabel
set ylabel "{p_e}(x)" # font "Helvetica, 26" offset -1,0
set samples 1000
plot 'histogram_0.txt' using ($2)-e:1 with boxes lc rgb 'grey', g(x+e) w l lw 1 lc rgb 'red'
#plot g(x+e) w l lw 1 lc rgb 'black'

#plot3
unset key
unset arrow
#set tmargin at screen 0.35
#set bmargin at screen 0.1
set size 1,0.45
set origin 0,0
#set xlabel 'x' font "Helvetica, 26" offset 40,0
set style arrow 8 heads size screen 0.008,90 lc rgb 'blue' lw 1 front
set arrow from 0,0.01 to sqrt(c**2-f**2),0.01 as 8
unset label
unset arrow 8
set label "{{/Symbol s}_q} = ({{{/Symbol s}_t}^2} - {{{/Symbol s}_e}^2})^{1/2}" at 1200,110000  font "Helvetica, 16"
set label 2 "(c)"  at 3500,350000 font "Helvetica,18"
#set label 3 "max{/Symbol \\133}p_q(x){/Symbol \\135} {/Symbol \\273} 9.38{/Symbol \\264}10^{-5}" at -75, 0.044 font "Helvetica, 16"
set style arrow 3 backhead filled  size screen 0.015,15,45 ls 1 lc rgb 'blue' lw 1 front
set arrow 3 from -3.5, 0.046 to -12,0.046 as 3
set yrange[0:400000]
unset ylabel
set ylabel "{p_q}(x)" # font "Helvetica, 26" offset -1,0
set xlabel "ADC value"

#set xtics ("0" -63, "32767" 0, "65535" 64) nomirror # font "Helvetica, 26" offset 0,-0.5

set xtics ("-2^{11}" -2048, "0" 0, "2^{11}-1" 2047) nomirror # font "Helvetica, 26" offset 0,-0.5

#plot (1/(2.507*sqrt(c**2-f**2)))*exp((-(x)**2)/(c**2-f**2)) lw 1 lc rgb 'blue'
plot a*exp((-(x)**2)/(c**2-f**2)) lw 1 lc rgb 'blue'