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 = 0.05
b = 68
c = 5
d = 0.1
e = 60
f = 5

set noytics
set noxtics


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

set xrange[-80:80]
#plot f(x), 'distribution_ele.txt' using 1:2
fit f(x) 'distribution_vac.txt' using 1:2 via a,b,c
fit g(x) 'distribution_ele.txt' using 1:2 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,0.01 to c,0.01 as 8
set yrange[0:0.06]
set ylabel "{p_t}(x)" # font "Helvetica, 26" offset -1,0
set label "{{/Symbol s}_t}" at 18,0.015 # font "Helvetica, 26"
set label 2 "(a)"  at 65,0.042 font "Helvetica,18"

plot 'distribution_vac.txt' using ($1)-b:2 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,0.025 to f,0.025 as 8
unset label
set label "{{/Symbol s}_e}" at 10,0.04 # font "Helvetica, 26"
set label 2 "(b)"  at 65,0.12 font "Helvetica,18"
set yrange[0:0.16]
unset ylabel
set ylabel "{p_e}(x)" # font "Helvetica, 26" offset -1,0
set samples 1000
plot 'distribution_ele.txt' using ($1)-e:2 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 16,0.015  font "Helvetica, 16"
set label 2 "(c)"  at 65,0.042 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:0.055]
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^{15}" -63, "0" 0, "2^{15}-1" 64) 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'