set term post eps enh color 12
set out "D_vs_S_U15etc.eps"
set size 0.5,0.45
set multiplot
set rmargin 1
set lmargin 8
set tmargin 1
set samples 1000



PS=0.8
PS1=0.7
PS2=0.6


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 style line 7 lt 1 lc rgb "#FFFF7E" lw 15
set style line 8 lt 1 lc rgb "#00FFFF" lw 15

set xlabel "s"
set ylabel "D" offset 1
set key samplen 2 spacing 1.3 left Left reverse
set xtics 0.1
set ytics 0.01
set mytics 2
set arrow 1 from log(2),0.025 to log(2),0.01 nohead lt 4 lc -1 front

set arrow 2 from 0.1,0.01 to 0.1,0.0112 nohead lt -1 front
set arrow 3 from 0.2,0.01 to 0.2,0.0112 nohead lt -1 front
set arrow 4 from 0.3,0.01 to 0.3,0.0112 nohead lt -1 front
set arrow 5 from 0.4,0.01 to 0.4,0.0112 nohead lt -1 front
set arrow 6 from 0.5,0.01 to 0.5,0.0112 nohead lt -1 front
set arrow 7 from 0.6,0.01 to 0.6,0.0112 nohead lt -1 front
set arrow 8 from 0.7,0.01 to 0.7,0.0112 nohead lt -1 front
set arrow 9 from 0.8,0.01 to 0.8,0.0112 nohead lt -1 front
set arrow 10 from 0.9,0.01 to 0.9,0.0112 nohead lt -1 front

set arrow 11 from log(2),0.013 to 0.885,0.013 lt -1 lw 2 size 0.03,30 front
set label 1 "charge" at 0.89,0.022 right front
set label 2 "physics" at 0.89,0.019 right front
set label 3 "dominant" at 0.89,0.016 right front

#set arrow from 0,0.0051 to 0.35,0.0051 nohead ls 7
#set arrow from 0.35,0.0051 to 0.693,0.0051 nohead ls 8

#set style fill pattern

plot [0:0.925][0.01:0.05]\
x<log(2)?0.0143:1/0 t"" w filledcurve y1=0.01 lc rgb "#cccccc",\
"Data/DMFT/results_U150mu0_ent_all.dt0" u 18:($18>log(2)?$10:1/0) t"" w filledcurve y1=0.01 lc rgb "#cccccc",\
"Data/BA/BA_U866_allT_D.dat" u 4:8 t"d = 1" w l lt 4 lw 3,\
"Data/DQMC_2D/U12.25L6_beta_s_Ek_D.dat" u 2:4 t"d = 2" w lp ls 2,\
"Data/DQMC_3D/U15L6dat0.dat" u 3:8  t "d = 3" w lp lt 3 lw 3 pt 12 ps PS1,\
"Data/DMFT/results_U150mu0_ent_all.dt0" u 18:($18>0.1?$10:1/0) t"DMFT" w lp lt 1 lw 3 pt 7 ps PS


#($1<=2.1?$3-sfit($1):sDMFTh_T($1)

set size 0.19,0.19
set origin 0.24,0.24
set xlabel "{/=10 s}" offset 0,0.5
set ylabel "{/=10 D}" offset 1.5
set format "{/=10 %g}"
set nokey
set mytics 1

set label 1 "d&{/=4 x}=&{/=4 x}1" at 0.91,0.009 right
unset label 2
unset label 3

set arrow 1 from log(2),0 to log(2),0.0235 nohead lt 4 lc -1 front
unset arrow 2
unset arrow 3
unset arrow 4
unset arrow 5
unset arrow 6
unset arrow 7
unset arrow 8
unset arrow 9
unset arrow 10
unset arrow 11

set lmargin 2
set bmargin 2

set label "{/=8 U/&{/=1 x}t = 8.66}" at 0.54,0.0265 tc lt 4 #rotate by -12
set label "{/=8 U/&{/=1 x}t = 12}" at 0.54,0.016 tc lt 3 #rotate by -10
set label "{/=8 U/&{/=1 x}t = 25}" at 0.54,0.0065 tc lt 1 # rotate by -2

plot [0.45:0.95][0:0.0325]\
"Data/BA/BA_U866_allT_D.dat" u 4:8 t"U/&{/=4 x}t = 8.66" w l lt 4 lw 3,\
"Data/BA/BA_U1200_allT_D.dat" u 4:8 w l lt 3 lw 3,\
"Data/BA/BA_U2500_allT_D.dat" u 4:8 w l lt 1 lw 3