# uses entropy based on Paiva's QMC data
# uses Rost's 2d QMC data
set term post eps enh color 12
set out "D_S_vs_T_U15etc.eps"
set size 0.5,0.45
set rmargin 1
set lmargin 7
set samples 1000

set multiplot

set style line 2 lt 2 lc rgb "#00BB00" lw 3


set size 0.5,0.20
set origin 0,0.05
set bmargin 1
set tmargin 0
set ytics 0.2
set mytics 2
set xtics 0.5
set mxtics 2

set xlabel "T/t{/Symbol \327}(Z/6)^{1/2}" offset 0,0.2
set ylabel "s" offset 1
set key samplen 2 spacing 1.2 Left reverse at 2.65,0.62

#set label 1 "{/=16 b)}" at 1.8,0.65
set label 1 "{/=16 b)}" at 0.08,0.79
#set label 2 "{/=11 Heisenberg}" at 1.53,0.62
set label 2 "{/=11 Heisenberg}" at 1.79,0.695

PS=0.8
PS1=0.6
PS2=0.5



set style line 2 lt 2 lc rgb "#00BB00" lw 3 pt 4 ps PS2
set style line 5 lt 5 lc rgb "#00BBBB" lw 3
set arrow from 0,log(2) to 0.4,log(2) nohead lt 1 lc -1 lw 1
#lc rgb "#404040" lw 2


#set arrow from 1.77,0.1 to 1.91,0.1 nohead lt -1

plot [0:2.25][0:0.875]\
"Data/BA/BA_U866_allT_D.dat" u ($2*sqrt(3)):4 t"d = 1 {/=11(chain, U/t = 8.66)}" w l lt 4 lw 3,\
"Data/DQMC_2D/U12.25L6_beta_s_Ek_D.dat" u (sqrt(1.5)/$1):2 t"d = 2 {/=11(square, U/t = 12.25)}" w lp ls 2,\
"Data/DQMC_3D/U15L6dat0.dat" u 2:(2<0?3:1/0) t "d = 3 {/=11(cubic, U/t = 15)}" w lp lt 3 lw 3 pt 12 ps PS1,\
"Data/DMFT/results_U150mu0_ent_all.dt0" u (1/$2):($18>0.1?$18:1/0) t"DMFT {/=11(cubic, U/t = 15)}" w lp lt 1 lw 3 pt 7 ps PS,\
"Data/DQMC_3D/U15L6dat0.dat" u 2:3 t "" w lp lt 3 lw 3 pt 12 ps PS1,\
"Data/Heisenberg/fig14_heis_corr.dat" u (6/$1):(6/$1<1.8?$2:1/0) t"" w l lt 7 lw 3



set size 0.5,0.20
set origin 0,0.265
set bmargin 0
set tmargin 1
set xlabel ""
set format x ""
set ylabel "D" offset 2
set nokey
#set label "{/=14 (hyper)cubic lattice, n=1}" at 0.4,0.031

#set label 1 "{/=16 a)}" at 1.8,0.015
set label 1 "{/=16 a)}" at  0.08,0.035
set label 2 ""

set ytics 0.01
set mytics 2

fac=1
PS=0.8
PS1=0.6
PS2=0.5
f(beta)=1/(beta+1)

set arrow from 0,0.0145 to 0.4,0.0145 nohead lt 4 lc rgb "#404040" lw 2

plot [0:2.25][0.01:0.0375]\
"Data/BA/BA_U866_allT_D.dat" u ($2*sqrt(3)):8 t"BA: d=1, U/t=8.66" w l lt 4 lw 3,\
"Data/DQMC_2D/U12.25L6_beta_s_Ek_D.dat" u (sqrt(1.5)/$1):4 t"d=2, U/t=12.25" w lp ls 2,\
"Data/DMFT/results_U150mu0_ent_all.dt0" u (1/$2):($18>0.1?$10:1/0) t"DMFT: d=3, U/t=15" w lp lt 1 lw 3 pt 7 ps PS,\
"Data/DQMC_3D/U15L6dat0.dat" u 2:8 t "QMC: d=3, U/t=15" w lp lt 3 lw 3 pt 12 ps PS1


