# Choose fiducial values for parameters (f(z) is defined later)
fiducial = {
'aperp': 1.,
'apar': 1.,
'bHI0': 0.702,
'A': 1.,
'sigma_nl': 7. #14. #0.5
}
# Load cosmology and experimental settings
cosmo = rf.experiments.cosmo
# Precompute cosmological functions and derivs.
camb_matterpower = "/home/phil/oslo/iswfunction/cosmomc/camb/testX_matterpower.dat"
cosmo_fns, cosmo = rf.precompute_for_fisher(cosmo, camb_matterpower)
H, r, D, f = cosmo_fns
# Get wiggles function and smoooth P(k)
k = np.logspace(np.log10(0.02), np.log10(0.4), 2000)
fbao = cosmo['fbao'](k)
pksmooth = cosmo['pk_nobao'](k)
# Plot results
P.subplot(111)
#P.plot(k, k**1. * pksmooth)
#P.plot(k, k**1. * pksmooth * (1. + fbao))
P.plot(k, fbao, 'k-', lw=1.5)
P.xlim((np.min(k), np.max(k)))
P.xscale('log')
import euclid
cosmo = rf.experiments.cosmo
#names = ["GBT", "BINGO", "WSRT", "APERTIF", "JVLA", "ASKAP", "KAT7", "MeerKAT", "SKA1", "SKAMID", "SKAMID_COMP", "iSKAMID", "iSKAMID_COMP", "SKA1_CV"]
#names = ["SKA1", "SKAMID", "SKAMID_COMP", "iSKAMID", "iSKAMID_COMP"]
names = [
"SKAMID",
] # "SKA1"] # , "iSKAMID_COMP"]
#ls = ['k-', 'r-', 'b--', 'm-', 'c--']
cols = ['r', 'g', 'c']
colours = ['#22AD1A', '#3399FF', '#ED7624']
cosmo_fns, cosmo = rf.precompute_for_fisher(rf.experiments.cosmo,
"camb/rf_matterpower.dat")
H, r, D, f = cosmo_fns
# Fiducial value and plotting
fig = P.figure()
ax1 = fig.add_subplot(111)
for k in range(len(names)):
root = "output/" + names[k]
# Load cosmo fns.
dat = np.atleast_2d(np.genfromtxt(root + "-cosmofns-zc.dat")).T
zc, Hc, dAc, Dc, fc = dat
zs, Hs, dAs, Ds, fs = np.genfromtxt(root + "-cosmofns-smooth.dat").T
kc = np.genfromtxt(root + "-fisher-kc.dat").T
#!/usr/bin/python
"""
OBSOLETE
"""
import numpy as np
import pylab as P
from rfwrapper import rf
from radiofisher import experiments as e
from units import *
expt = e.SKAMID
expt['Sarea'] /= 6.
cosmo_fns, cosmo = rf.precompute_for_fisher(e.cosmo, "camb/rf_matterpower.dat")
H, r, D, f = cosmo_fns
z = np.linspace(1e-2, 3., 300)
rnu = C * (1. + z)**2. / H(z) # Perp/par. dist. scales
Dmax = 100e3 # 100 km max. baseline
Dmin = 15. # 15m dish diameter (actually, this would give FOV, not Dmin)
# INTERFEROM.
kmin_int = 2. * np.pi * Dmin * (1420.0e6) / (3e8 * r(z) * (1. + z))
kmax_int = 2. * np.pi * Dmax * (1420.0e6) / (3e8 * r(z) * (1. + z))
# SINGLE-DISH
Vphys = expt['Sarea'] * (expt['survey_dnutot'] /
expt['nu_line']) * r(z)**2. * rnu
kmin = 2. * np.pi / Vphys**(1. / 3.)