se

set terminal pdfcairo color enhanced font "Arial,12"
set samples 1000
set ylabel 'Intensity'
set xlabel 'energy transfer / meV'
#set style fill transparent solid 0.50 noborder
set style fill solid 0.50 noborder
set key vertical maxrows 2 right top
set noautoscale xy
set noytics

set output "se3.pdf"

f10(x)=(i10/(sqrt(2*pi)*b10))*exp((-(x-c10)**2)/(2*b10**2))+d10
i10 = 0.0291819
b10 = 0.614285
c10 = 6.37132
d10 = 0.0459446
fit [4.8:7.6][] f10(x) "10.txt" using 1:2 via i10,b10,c10,d10

f08(x)=(i08/(sqrt(2*pi)*b08))*exp((-(x-c08)**2)/(2*b08**2))+d08
i08 = 0.0366142
b08 = 0.567193
c08 = 6.41539
d08 = 0.0502077
fit [4.8:7.6][] f08(x) "08.txt" using 1:2 via i08,b08,c08,d08

f04(x)=(i04/(sqrt(2*pi)*b04))*exp((-(x-c04)**2)/(2*b04**2))+d04
i04 = 0.0109621
b04 = 0.39314
c04 = 6.46416
d04 = 0.0554256
fit [4.8:7.6][] f04(x) "04.txt" using 1:2 via i04,b04,c04,d04

f00(x)=(i00/(sqrt(2*pi)*b00))*exp((-(x-c00)**2)/(2*b00**2))+d00
i00 = 0.0180383
b00 = 0.446996
c00 = 6.52505
d00 = 0.0403075
fit [4.8:7.6][] f00(x) "00.txt" using 1:2 via i00,b00,c00,d00

set label sprintf("{/:Italic T} = 200K \n fit from 4.8 to 7.6") at graph 0.57,0.95 right 

plot [2:10][0:0.16] \
    "00.txt" u 1:($2+0.070) smooth csplines with filledc y=0 lc rgb "#FF4B00" notitle,\
    "04.txt" u 1:($2+0.030) smooth csplines with filledc y=0 lc rgb "#F6AA00" notitle,\
    "08.txt" u 1:($2-0.010) smooth csplines with filledc y=0 lc rgb "#03AF7A" notitle,\
    "10.txt" u 1:($2-0.030) smooth csplines with filledc y=0 lc rgb "#005AFF" notitle,\
    "00.txt" u 1:($2+0.070):3 w e pt 6 ps 0.7 lc rgb "#CC3900" title 'Se : 0.0',\
    "04.txt" u 1:($2+0.030):3 w e pt 6 ps 0.7 lc rgb "#C18400" title 'Se : 0.4',\
    "08.txt" u 1:($2-0.010):3 w e pt 6 ps 0.7 lc rgb "#027752" title 'Se : 0.8',\
    "10.txt" u 1:($2-0.030):3 w e pt 6 ps 0.7 lc rgb "#004ACC" title 'Se : 1.0',\
    f00(x)+0.070 w l lc rgb "#992B00" notitle,\
    f04(x)+0.030 w l lc rgb "#8E6100" notitle,\
    f08(x)-0.010 w l lc rgb "#014530" notitle,\
    f10(x)-0.030 w l lc rgb "#003899" notitle

save 'se.log'

temperature.g

set terminal pdfcairo color enhanced font "Arial,12"
set encoding iso
set size square
set lmargin 0
set rmargin 0
set terminal pdf size 5.4,5
set key right at graph 0.97,0.94
set xlabel "{/:Italic E} (meV)"
set ylabel "Intensity"

circ = 0.6
set linetype 1 pt 6 ps circ lc "#FF4B00"
set linetype 2 lc "#FF4B00"
set linetype 3 pt 6 ps circ lc "#005AFF"
set linetype 4 lc "#005AFF"
set linetype 5 pt 6 ps circ lc "#03AF7A"
set linetype 6 lc "#03AF7A"
set linetype 7 pt 6 ps circ lc "#4DC4FF"
set linetype 8 lc "#4DC4FF"

x = "10"

set label 1 right at graph 0.96,0.96 "{/:Italic x} = ".x[1:1].".".x[2:2].", {/:Italic Q} = 2.8 (\305)"
set output "x".x.".pdf"
plot\
"x".x."t10.txt" u 1:2:(sqrt($2)) title "{/:Italic T} = 10 K" w e,\
"" u 1:2 notitle w l,\
"x".x."t100.txt" u 1:2:(sqrt($2)) title "{/:Italic T} = 100 K" w e,\
"" u 1:2 notitle w l,\
"x".x."t200.txt" u 1:2:(sqrt($2)) title "{/:Italic T} = 200 K" w e,\
"" u 1:2 notitle w l,\
"x".x."t300.txt" u 1:2:(sqrt($2)) title "{/:Italic T} = 300 K" w e,\
"" u 1:2 notitle w l,\

set terminal pop

x.g

set terminal pdfcairo color enhanced font "Arial,12"
set encoding iso
set size square
set lmargin 0
set rmargin 0
set terminal pdf size 5.4,5
set key right at graph 0.97,0.94
set xlabel "{/:Italic E} (meV)"
set ylabel "Intensity"

circ = 0.6
set linetype 1 pt 6 ps circ lc "#FF4B00"
set linetype 2 lc "#FF4B00"
set linetype 3 pt 6 ps circ lc "#005AFF"
set linetype 4 lc "#005AFF"
set linetype 5 pt 6 ps circ lc "#03AF7A"
set linetype 6 lc "#03AF7A"
set linetype 7 pt 6 ps circ lc "#4DC4FF"
set linetype 8 lc "#4DC4FF"

t = "300"

set label 1 right at graph 0.96,0.96 "{/:Italic T} = ".t.", {/:Italic Q} = 2.8 (\305)"
set output "t".t.".pdf"
plot\
"x00t".t.".txt" u 1:2:(sqrt($2)) title "{/:Italic x} = 0.0" w e,\
"" u 1:2 notitle w l,\
"x08t".t.".txt" u 1:2:(sqrt($2)) title "{/:Italic x} = 0.8" w e,\
"" u 1:2 notitle w l,\
"x10t".t.".txt" u 1:2:(sqrt($2)) title "{/:Italic x} = 1.0" w e,\
"" u 1:2 notitle w l,\

set terminal pop