# 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.)