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"],