distances_bm Namespace Reference

Functions
def	do_all (z_arr, cpar, prefix)

Variables
float	TCMB = 2.725
int	PLOT_STUFF = 0
int	WRITE_STUFF = 1
int	FS = 16

Function Documentation

def distances_bm.do_all	(		z_arr,
			cpar,
			prefix
	)

Definition at line 13 of file distances_bm.py.

def do_all(z_arr,cpar,prefix) :
    pcs=csm.PcsPar()
    pcs.background_set(cpar['om'],cpar['ol'],cpar['ob'],cpar['w0'],cpar['wa'],cpar['hh'],TCMB)

    a_arr=1./(z_arr+1)
    chi_arr=np.array([pcs.radial_comoving_distance(a) for a in a_arr])

    if PLOT_STUFF==1 :
        plt.plot(z_arr,chi_arr); plt.xlabel('$z$',fontsize=FS); plt.ylabel('$\\chi(z)\\,[{\\rm Mpc}\\,h^{-1}]$',fontsize=FS); plt.show()

    if WRITE_STUFF==1 :
        np.savetxt(prefix+"_chi.txt",np.transpose([z_arr,chi_arr]),header="[1] z, [2] chi(z) (Mpc/h)")

z_arr=np.array([0.,1.,2.,3.,4.,5.])

cpar_model1={'om': 0.3,'ol': 0.7 ,'ob':0.05,'hh': 0.7,'s8': 0.8,'ns': 0.96,'w0': -1.0, 'wa': 0.0}
cpar_model2={'om': 0.3,'ol': 0.7 ,'ob':0.05,'hh': 0.7,'s8': 0.8,'ns': 0.96,'w0': -0.9, 'wa': 0.0}
cpar_model3={'om': 0.3,'ol': 0.7 ,'ob':0.05,'hh': 0.7,'s8': 0.8,'ns': 0.96,'w0': -0.9, 'wa': 0.1}
cpar_model4={'om': 0.3,'ol': 0.75,'ob':0.05,'hh': 0.7,'s8': 0.8,'ns': 0.96,'w0': -0.9, 'wa': 0.1}
cpar_model5={'om': 0.3,'ol': 0.65,'ob':0.05,'hh': 0.7,'s8': 0.8,'ns': 0.96,'w0': -0.9, 'wa': 0.1}

do_all(z_arr,cpar_model1,"model1")
do_all(z_arr,cpar_model2,"model2")
do_all(z_arr,cpar_model3,"model3")
do_all(z_arr,cpar_model3,"model4")
do_all(z_arr,cpar_model3,"model5")

Variable Documentation

int distances_bm.FS = 16

Definition at line 11 of file distances_bm.py.

int distances_bm.PLOT_STUFF = 0

Definition at line 9 of file distances_bm.py.

float distances_bm.TCMB = 2.725

Definition at line 8 of file distances_bm.py.

int distances_bm.WRITE_STUFF = 1

Definition at line 10 of file distances_bm.py.