colours = [ ['#CC0000', '#F09B9B'],
['#1619A1', '#B1C9FD'],
['#5B9C0A', '#BAE484'],
['#FFB928', '#FFEA28'] ]
# Define values for interferom./single-dish
Ddish = 15.
Dmin = 15.
Dmax = 1000.
dnutot = 600. / 1420. # Dimensionless
Sarea = 25e3 * (D2RAD)**2.
# Get r(z)
HH, rr, DD, ff = rf.background_evolution_splines(cosmo)
zz = np.linspace(0., 3., 1000)
r = rr(zz)
rnu = C*(1.+zz)**2. / HH(zz)
l = (1.+zz) * 3e8 / 1420e6 # metres
# Calculate horizon size
zeq = 3265. # From WMAP 9; horizon size is pretty insensitive to this
om = cosmo['omega_M_0']
ol = cosmo['omega_lambda_0']
h = cosmo['h']
orad = om / (1. + zeq)
aa = np.linspace(0., 1., 10000)
_z = 1./aa - 1.
integ = 1. / np.sqrt(om*aa + ol*aa**4. + orad)
import numpy as np
import pylab as P
from rfwrapper import rf
import matplotlib.patches
import matplotlib.cm
import os
from radiofisher import euclid
cosmo = rf.experiments.cosmo
names = ["EuclidRef_paper", "exptL_paper", "aexptM_paper", "exptS_paper"]
colours = ['#CC0000', '#1619A1', '#5B9C0A', '#990A9C'] # DETF/F/M/S
labels = ['DETF IV', 'Facility', 'Stage II', 'Stage I']
# Get f_bao(k) function
cosmo_fns = rf.background_evolution_splines(cosmo)
cosmo = rf.load_power_spectrum(cosmo, "cache_pk.dat", force_load=True)
fbao = cosmo['fbao']
# Fiducial value and plotting
fig = P.figure()
axes = [
fig.add_subplot(411),
fig.add_subplot(412),
fig.add_subplot(413),
fig.add_subplot(414)
]
for k in range(len(names)):
root = "output/" + names[k]
#!/usr/bin/python
"""
Calculate and plot the comoving volumes of some surveys.
"""
import numpy as np
import pylab as P
import scipy.integrate
import scipy.interpolate
from rfwrapper import rf
C = 3e5
cosmo = rf.experiments.cosmo
# Precalculate background evolution
H, r, D, f = rf.background_evolution_splines(cosmo, zmax=10., nsamples=500)
_z = np.linspace(0., 10., 1000)
_vol = C * scipy.integrate.cumtrapz(r(_z)**2. / H(_z), _z, initial=0.)
_vol *= (np.pi / 180.)**2. / 1e9 # per deg^2, in Gpc^3
vol = scipy.interpolate.interp1d(_z, _vol, kind='linear', bounds_error=False)
def Vsurvey(zmin, zmax, sarea):
return sarea * (vol(zmax) - vol(zmin))
# zmin, zmax, sarea, name
surveys = [
[0., 0.8, 10e3, "BOSS"],
[1.9, 3.5, 420., "HETDEX"],
[0.1, 1.9, 14e3, "DESI"],
[0.6, 2.1, 15e3, "Euclid"],
