set term post eps enh color 12
set out "D_Ekin_cubic_s_DMFT_QMC_sN.eps"
set size 0.5,0.54
set samples 1000
set multiplot

set size 0.49,0.23
set origin 0.01,0.06
set rmargin 0.5
set lmargin 7
set tmargin 0
set bmargin 0

set ytics 0.5
set mytics 2
set xtics 0,0.2,1.4
set mxtics 2

fac=1
PS=0.8
PS1=0.8
PS2=0.6

s_DMFT_U2=0.0368944
s_DMFT_U4=0.233629
s_DMFT_U6=0.471987
s_DMFT_U8=0.631025
s_DMFT_U10=0.692283
s_DMFT_U12=0.6920
s_DMFT_U15=0.6921
s_QMC_U4=0.161456
s_QMC_U6=0.319649
s_QMC_U8=0.350211
s_QMC_U10=0.3072
s_QMC_U12=0.3409
s_QMC_U15=0.3466


set xlabel "s"  offset 0,0.5
set ylabel "E_{kin} / t" offset 1

set nokey

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

set label 1 "{/=16 (a)}" at screen 0,0.5
set label 2 "{/=16 (b)}" at screen 0,0.27


plot [0:1.95][-2.1:0] \
"Data/DMFT/DMFT_TN_SN_D_Ekin.dat" u 3:5 t"" w l lt 1 lc -1 lw 1,\
"Data/DQMC_3D/TN_SN_D_Ekin_QMC.dat" u 3:5 t"" w l lt 1 lc -1 lw 1,\
"Data/DMFT/U0_cubic_E_s.dat" u 3:2 t"{/=10 U/t = 0}" w l lt -1,\
"Data/DMFT/results_U020mu0_ent_all.dt0" u 18:13 t"" w l lt 1 lc 1 lw 3,\
"Data/DQMC_3D/threedU2.dat" u 3:($4-2*$8) t"" w p lc 1 lw 3 pt 4 ps PS2,\
"Data/DMFT/results_U040mu0_ent_all.dt0" u 18:($18>0.03?$13:1/0) t"" w l ls 2,\
"Data/DQMC_3D/threedU4.dat" u 3:($4-4*$8) t"" w p lc rgb "#00BB00" lw 3 pt 2 ps PS1,\
"Data/DMFT/results_U060mu0_ent_all.dt0" u 18:($18>0.05?$13:1/0) t"" w l lt 3 lw 3,\
"Data/DQMC_3D/threedU6.dat" u 3:($4-6*$8) t"" w p lt 3 lw 3 pt 6 ps PS1,\
"Data/DMFT/results_U080mu0_ent_all.dt0" u 18:($18>0.065?$13:1/0) t"" w l lt 4 lw 3,\
"Data/DQMC_3D/threedU8.dat" u 3:($3>0.03?($4-8*$8):1/0) t"" w p lt 4 lw 3 pt 10 ps PS1,\
"Data/DMFT/results_U100mu0_ent_all.dt0" u 18:($18>0.075?$13:1/0) t"" w l lt 5 lw 3,\
"Data/DQMC_3D/threedU10_Th.dat" u ($1):7 t"" w p lt 5 lw 3 pt 8 ps PS1,\
"Data/DMFT/results_U150mu0_ent_all.dt0" u 18:($18>0.1?$13:1/0) t"" w l lt 6 lc 8 lw 3,\
"Data/DQMC_3D/U15L6dat0.dat" u 3:6 t "" w p lt 6 lc 8 lw 3 pt 12 ps PS1


set size 0.49,0.23
set origin 0.01,0.29
set tmargin 0
set bmargin 0

set xlabel ""
set format x ""
set ylabel "D" offset 1.8

#set key samplen 2 spacing 1.2 at 1.79,0.205
set key samplen 2 spacing 1.2 

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


D_U15(x) = 0.0145 + exp ( -1.4233489     -4.6252736  *x+x**2*   -10.261082 )
D_U10(x) = 0.044 + exp ( -1.5714512     -3.1333736   *x+x**2*  -8.5874696 )
D_U8(x) = 0.072 + exp (  -1.7367676     -2.4327978  *x+x**2*   -9.8342685)
D_U6(x) = 0.108 + exp ( -1.963635      -2.0944032  *x+x**2*    -10.308031)
D_U4(x) = 0.150 + exp ( -2.3089814     -2.0638882   *x+x**2*   -9.5611966)
D_U2(x) = 0.198 + exp (-2.9578806     -2.0634092 *x+x**2*    -8.8929853 )

# plot [0:1] "Data/DMFT/results_U040mu0_ent_all.dt0" u 18:10 w lp, x<0.84?0.1565:1/0, "Data/DMFT/results_U060mu0_ent_all.dt0" u 18:10 w lp, x<0.81?0.1138:1/0, "Data/DMFT/results_U020mu0_ent_all.dt0" u 18:10 w lp, x<0.85?0.2027:1/0, "Data/DMFT/results_U080mu0_ent_all.dt0" u 18:10 w lp, x<0.76?0.076:1/0, "Data/DMFT/results_U100mu0_ent_all.dt0" u 18:10 w lp, x<0.715?0.0454:1/0, "Data/DMFT/results_U150mu0_ent_all.dt0" u 18:10 w lp, x<0.694?0.0145:1/0


set ytics 0.05
set mytics 2
set mxtics 2

DY=0.02
#set style arrow 1 head  size 0.03,30 lt 1 lw 2 lc -1 front
#set arrow from s_DMFT_U2,0.215+DY to s_DMFT_U2,0.215 as 1
#set arrow from s_DMFT_U2,0.21-DY to s_DMFT_U2,0.21 as 1
#set arrow from s_DMFT_U4,0.179+DY to s_DMFT_U4,0.179 as 1
#set arrow from s_QMC_U4,0.165-DY to s_QMC_U4,0.165 as 1
#set arrow from s_DMFT_U6,0.132+DY to s_DMFT_U6,0.132 as 1
#set arrow from s_QMC_U6,0.115-DY to s_QMC_U6,0.115 as 1
#set arrow from s_DMFT_U8,0.084+DY to s_DMFT_U8,0.084 as 1
#set arrow from s_QMC_U8,0.07-DY to s_QMC_U8,0.07 as 1
#set arrow from s_DMFT_U10,0.05+DY to s_DMFT_U10,0.05 as 1
#set arrow from s_QMC_U15,0.044-DY to s_QMC_U15,0.044 as 1
#set arrow from s_DMFT_U15,0.02+DY to s_DMFT_U15,0.02 as 1
#set arrow from s_QMC_U15,0.018-DY to s_QMC_U15,0.018 as 1


plot [0:1.95][0:0.25] \
"Data/DMFT/DMFT_TN_SN_D_Ekin.dat" u 3:4 t"" w l lt 1 lc -1 lw 1,\
"Data/DQMC_3D/TN_SN_D_Ekin_QMC.dat" u 3:4 t"" w l lt 1 lc -1 lw 1,\
"Data/DMFT/U0_cubic_E_s.dat" u 3:(1/0) t"U/t = 0" w l lt -1,\
"Data/DQMC_3D/threedU2.dat" u 3:(1/0) t"U/t = 2" w lp lt 1 lc 1 lw 3 pt 4 ps PS2,\
"Data/DQMC_3D/threedU4.dat" u 3:(1/0) t"U/t = 4" w lp lt 2 lc rgb "#00BB00" lw 3 pt 2 ps PS1,\
"Data/DQMC_3D/threedU6.dat" u 3:(1/0) t"U/t = 6" w lp lt 3 lc 3 lw 3 pt 6 ps PS1,\
"Data/DQMC_3D/threedU8.dat" u 3:(1/0) t"U/t = 8" w lp lt 4 lc 4 lw 3 pt 10 ps PS1,\
"Data/DQMC_3D/threedU10_Th.dat" u ($1):(1/0) t"U/t = 10" w lp lt 5 lc 5 lw 3 pt 8 ps PS1,\
"Data/DQMC_3D/U15L6dat0.dat" u 3:(1/0) t"U/t = 15" w lp lt 6 lc 8 lw 3 pt 12 ps PS1,\
"Data/DMFT/results_U020mu0_ent_all.dt0" u 18:10 t"" w l lt 1 lc 1 lw 3,\
"" u 18:(D_U2($2/(1+$2))) t"" smooth csplines lt 4 lc rgb "#404040",\
"Data/DQMC_3D/threedU2.dat" u 3:8 t"" w p lc 1 lw 3 pt 4 ps PS2,\
"Data/DMFT/results_U040mu0_ent_all.dt0" u 18:($18>0.03?$10:1/0) t"" w l ls 2,\
"" u 18:(D_U4($2/(1+$2))) t"" smooth csplines lt 4 lc rgb "#404040",\
"Data/DQMC_3D/threedU4.dat" u 3:8 t"" w p lc rgb "#00BB00" lw 3 pt 2 ps PS1,\
"Data/DMFT/results_U060mu0_ent_all.dt0" u 18:($18>0.05?$10:1/0) t"" w l lt 3 lw 3,\
"" u 18:(D_U6($2/(1+$2))) t"" smooth csplines lt 4 lc rgb "#404040",\
"Data/DQMC_3D/threedU6.dat" u 3:8 t"" w p lt 3 lw 3 pt 6 ps PS1,\
"Data/DMFT/results_U080mu0_ent_all.dt0" u 18:($18>0.065?$10:1/0) t"" w l lt 4 lw 3,\
"" u 18:(D_U8($2/(1+$2))) t"" smooth csplines lt 4 lc rgb "#404040",\
"Data/DQMC_3D/threedU8.dat" u 3:($3>0.03?$8:1/0) t"" w p lt 4 lw 3 pt 10 ps PS1,\
"Data/DMFT/results_U100mu0_ent_all.dt0" u 18:($18>0.075?$10:1/0) t"" w l lt 5 lw 3,\
"" u 18:(D_U10($2/(1+$2))) t"" smooth csplines lt 4 lc rgb "#404040",\
"Data/DQMC_3D/threedU10_Th.dat" u ($1):5 t"" w p lt 5 lw 3 pt 8 ps PS1,\
"Data/DMFT/results_U150mu0_ent_all.dt0" u 18:($18>0.1?$10:1/0) t"" w l lt 6 lc 8 lw 3,\
"" u 18:(D_U15($2/(1+$2))) t"" smooth csplines lt 4 lc rgb "#404040",\
"Data/DQMC_3D/U15L6dat0.dat" u 3:8 t "" w p lt 6 lc 8 lw 3 pt 12 ps PS1

set size 0.18,0.12
set origin 0.308,0.10

YOFF=-0.06
set label "{/=9 0}" at 0,YOFF center
set label "{/=9 5}" at 5,YOFF center
set label "{/=9 10}" at 10,YOFF center
set label "{/=9 15}" at 15,YOFF center

XOFF=-0.8
set label "{/=9 0.6}" at XOFF,0.6 right
set label "{/=9 0.4}" at XOFF,0.4 right
set label "{/=9 0.2}" at XOFF,0.2 right
set label "{/=9 0}" at XOFF,0 right

set xlabel "{/=9 U / t}" offset 0,0.4
#set ylabel "{/=9 s_N}" offset 0.8,-0.2
set ylabel "{/=9 s_N}" offset -0.8,-0.2

set ytics 0.2
set mytics 2
set nokey
set xtics 5
set mxtics 1
set format ""

PS2=0.7

set label "{/=8 DMFT}" at 10,0.62 tc lt 1
set label "{/=8 DQMC}" at 10,0.24 tc lt 3

plot [0:17.5][0:0.75]\
log(2)/2 t"" lt 5 lc rgb "#404040",\
log(2) t"" lt 9,\
"Data/DMFT/AF_mu0_TN_SN.dat" u 1:3 t"DMFT" w lp lt 1 lw 2 pt 6 ps PS2,\
"Data/DQMC_3D/TN_SN_QMC.dat" u 1:3 t"3D QMC" w lp lt 3 lw 2 pt 12 ps PS2
