# STACK THEM ALL
set terminal postscript eps color enhanced solid font "Helvetica,20" size 12cm,13cm
set output 'pddeg2.eps'
set multiplot layout 3, 1
set colorsequence classic


infile3 = 'data/13pi_PDD_2019-01-11.dat'
infile0 = 'data/3pi_PDD_2019-01-10.dat'
infile1 = 'data/5pi_PDD_2019-01-10.dat'
infile2 = 'data/7pi_PDD_2019-01-11.dat'
infile4 = 'data/spinecho_2019-01-10.dat'
infile5 = 'data/15pi_PDD_2019-01-11.dat'
infiles0 = 'data/3pi_simulations.dat'
infiles1 = 'data/5pi_simulations.dat'
infiles2 = 'data/13pi_simulations.dat'

a0=0.5
b0 = -0.4
tau = 1

decay1(x) = a0 - exp(-x**2/tau**2)*b0
#fit decay1(x) infile0 u (4*$1/1000):5:6 via a0,b0,tau

unset ylabel
unset xlabel
unset key
#
set xrange [0:3]
set yrange [50:100]
set ytics 25 nomirror
set xtics 1 nomirror
set format x ""
set lmargin 8
set rmargin 1
set tmargin at screen 0.98
set bmargin at screen 0.71
plot infiles0 u ($1/1000):(43*$2+50) w lines lc 3 lw 3 not, infile0 u (4*$1/1000):(100*$5):(100*$6) w yerr lc 1 ps 1.2 pt 7 t '3 pi'
#
set tmargin at screen 0.66
set bmargin at screen 0.39
set ylabel "Probability in F=2 (%)" offset 0,0

a2=0.5
b2 = -0.2
tau2 = 1
decay2(x) = a2 - exp(-x**2/tau2**2)*b2
#fit decay2(x) infile1 u (6*$1/1000):5:6 via a2,b2,tau2

plot infiles1 u ($1/1000):(43*$2+50)  w lines lc 3 lw 3 not, infile1 u (6*$1/1000):(100*$5):(100*$6) w yerr lc 1 ps 1.2 pt 7 t '5 pi'

#
unset ylabel
set format x "%.0f"
set tmargin at screen 0.34
set bmargin at screen 0.07
set xlabel "Free evolution time {/Symbol t} (ms)" offset 0,0.5

a4= 0.3
b4 = -0.2
tau4 = 1.3
decay4(x) = a4 - exp(-x**2/tau4**2)*b4
#fit decay4(x) infile5 u (14*$1/1000):5:6 via a4,b4

plot infiles2 u ($1/1000):(43*$2+50) w lines lc 3 lw 3 not, infile3 u (14*$1/1000):(100*$5):(100*$6) w yerr lc 1 ps 1.2 pt 7 t '13 pi'
unset multiplot
pause -1