# SM file to plot the average delpk spectra of the 1D DM "spectra".
#
erase
window 1 1 1 1
ctype default
expand 1.1
lweight 2
ltype 0
ptype 4 3
location 4500 31000 3500 31000
notation -4 4 -4 4
define TeX_strings 1
#
window 1 -2 1 1
limits -0.99 2.0 -2.0 1.99
ticksize -1 10 0 0
box
xlabel k (h/Mpc)
ylabel log_{10}[k P(k)/\pi]
expand 0.9
#
data OS_z3_Kf200.dat
read {k 1 pk 5} SET lk=LG(k) SET pk=LG(pk)
connect lk pk
data OS_z3_Kf100.dat
read {k 1 pk 5} SET lk=LG(k) SET pk=LG(pk)
connect lk pk
data OS_z3_Kf050.dat
read {k 1 pk 5} SET lk=LG(k) SET pk=LG(pk)
connect lk pk
data OS_z3_Kf025.dat
read {k 1 pk 5} SET lk=LG(k) SET pk=LG(pk)
connect lk pk
#
ctype default
relocate -0.8 -0.50 putlabel 6 128^3/ 256^3, 30Mpc/h
data "OS30A/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
#
ctype blue
relocate -0.8 -0.75 putlabel 6 256^3/ 512^3, 30Mpc/h
data "OS30B/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
#
ctype cyan
relocate -0.8 -1.00 putlabel 6 512^3/1024^3, 30Mpc/h
data "OS30C/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
#
ctype magenta
relocate -0.8 -1.25 putlabel 6 512^3/1024^3, 20Mpc/h
data "OS20A/delpk_1.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
#
ctype red
relocate -0.8 -1.50 putlabel 6 256^3/ 512^3, 15Mpc/h
data "OS15B/delpk_1.dat"
read {k 1 d2 2} SET lk=LG(k) SET d2=LG(d2)
ptype 4 0 points lk d2
#
ctype default
expand 1.1
window 1 -2 1 2
limits -0.99 2.0 0.00001 50
ticksize -1 10 0 0
box 0 2 0 0
ylabel k P(k)/\pi
expand 0.9
#
relocate -0.8 45 putlabel 6  O&S, z=3
relocate -0.8 42 putlabel 6  Subrandom ICs
#
ctype default
data "OS30A/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k)
ptype 4 0 points lk d2
#
ctype blue
data "OS30B/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k)
ptype 4 0 points lk d2
#
ctype cyan
data "OS30C/delpk_7.dat"
read {k 1 d2 2} SET lk=LG(k)
ptype 4 0 points lk d2
#
ctype magenta
data "OS20A/delpk_1.dat"
read {k 1 d2 2} SET lk=LG(k)
ptype 4 0 points lk d2
#
ctype red
data "OS15B/delpk_1.dat"
read {k 1 d2 2} SET lk=LG(k)
ptype 4 0 points lk d2
#