Commit ef0f5225 authored by Plaszczynski Stephane's avatar Plaszczynski Stephane
Browse files

add some data for testig stability of Engine versions

parent 6c464673
This source diff could not be displayed because it is too large. You can view the blob instead.
Tcmb= 2.7255
file1='cl_classv2.3.2.fits'
t1=mrdfits(file1,1)
file2='cl_classv2.4.3.fits'
t2=mrdfits(file2,1)
file3='plikTT_bflike_pico.fits'
;file3='cl_classv2.4.3_nl.fits'
t3=mrdfits(file3,1)
;camb
readcol,'base_plikHM_TT_lowTEB.minimum.theory_cl',ell,dl,skip=1
;class
;file3='output/plikbf_cl_lensed.dat'
;readcol,file3,ell,dltt,skip=7
;restart @lmin=2
l=dindgen(n_elements(t1.tt))
llp=l*(l+1)/(2*!dpi)
dl1=llp*t1.tt*1e12
dl1=dl1[2:*]
dl2=llp*t2.tt*1e12
dl2=dl2[2:*]
dl3=llp*t3.tt*1e12
dl3=dl3[2:*]
;dl3=dltt*1e12*Tcmb^2
cols=['black','blue','red']
cgplot,ell,dl1-dl,xtit='l',/xs,ytit='$\Delta$Cl',yr=[-8,8]
cgplot,ell,dl2-dl,/over,col=cols[1]
cgplot,ell,dl3-dl,/over,col=cols[2]
cgplot,!x.crange,[0,0],/over
al_legend,[file1,file2,file3],textcol=cgcolor(cols),/top,/right
engine=class
par omega_b cosm 0.2224156E-01 0.00027 0.017 0.027
par omega_cdm cosm 0.1197735E+00 0.0026 0.09 0.15
#par 100*theta_s cosm 0.1040862E+01 0.6E-04 1.0170 1.0617
par H0 cosm 0.6730713E+02 1 50 80
par tau_reio cosm 0.7814866E-01 0.13E-01 0.01016 0.18221
par log(10^10A_s) cosm 0.3090373E+01 0.025 2.7 3.5
par n_s cosm 0.9658314E+00 0.0070 0.9 1.1
par A_planck nui 0.1000296E+01 0.01 0.9 1.1 #ycal
par A_cib_217 nui 0.6662739E+02 1 0 200 #OK Acib217
par xi_sz_cib nui 0.4657732E-01 .01 0 1 #xi_sz_cib
par A_sz nui 0.7141324E+01 1 0 10 #OK Atsz143
par ps_A_100_100 nui 0.2518416E+03 10 0 400 #OK Aps100
par ps_A_143_143 nui 0.3916351E+02 1 0 400 #OK Aps143
par ps_A_143_217 nui 0.3360542E+02 1 0 400 #OK Aps143x217
par ps_A_217_217 nui 0.9783567E+02 1 0 400 #OK Aps217
par ksz_norm nui 0.8641823E-03 1 0 10 #OK Aksz
par gal545_A_100 nui 0.7412626E+01 .2 0 50 #OK AdustTT100
par gal545_A_143 nui 0.8981278E+01 9.41 0 50 #OK AdustTT143
par gal545_A_143_217 nui 0.1753156E+02 19.69 0 100 #OK AdustTT143x217
par gal545_A_217 nui 0.8202580E+02 82.4 0 100 #OK AdustTT217
par calib_100T nui 0.9978923E+00 .2 0 3 #OK c100
par calib_217T nui 0.9959253E+00 .2 0 3 #OK c217
fix cib_index nui -1.3 .2 -5 5 #OK ncib
gausssum2 A_sz ksz_norm 1.6 1. 9.5 3
gauss1 calib_100T 0.9990004 0.001
gauss1 calib_217T 0.99501 0.002
gauss1 gal545_A_100 7 2
gauss1 gal545_A_143 9 2
gauss1 gal545_A_143_217 21 8.5
gauss1 gal545_A_217 80 20
gauss1 A_planck 1 0.0025
#class setup
class k_pivot 0.05
class lensing yes
class sBBN\ file bbn/sBBN.dat
precisionFile=class_pre/hpjul2.pre
clikfile=planck_data/hi_l/plik/plik_dx11dr2_HM_v18_TT.clik
clikfile2=planck_data/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5c_Ap.clik
#fitter
nitermax=50000
set_stra=2
set_tol=0.00001
#output options
computecontour=false
profile=false
remove_cosmo_limits=false
remove_limits=false
class N_ncdm 1
fix m_ncdm cosmo 0.06
class N_eff 2.046
engine=class
par omega_b cosm 0.2224156E-01 0.00027 0.017 0.027
par omega_cdm cosm 0.1197735E+00 0.0026 0.09 0.15
#par 100*theta_s cosm 0.1040862E+01 0.6E-04 1.0170 1.0617
par H0 cosm 0.6730713E+02 1 50 80
par tau_reio cosm 0.7814866E-01 0.13E-01 0.01016 0.18221
par log(10^10A_s) cosm 0.3090373E+01 0.025 2.7 3.5
par n_s cosm 0.9658314E+00 0.0070 0.9 1.1
par A_planck nui 0.1000296E+01 0.01 0.9 1.1 #ycal
par A_cib_217 nui 0.6662739E+02 1 0 200 #OK Acib217
par xi_sz_cib nui 0.4657732E-01 .01 0 1 #xi_sz_cib
par A_sz nui 0.7141324E+01 1 0 10 #OK Atsz143
par ps_A_100_100 nui 0.2518416E+03 10 0 400 #OK Aps100
par ps_A_143_143 nui 0.3916351E+02 1 0 400 #OK Aps143
par ps_A_143_217 nui 0.3360542E+02 1 0 400 #OK Aps143x217
par ps_A_217_217 nui 0.9783567E+02 1 0 400 #OK Aps217
par ksz_norm nui 0.8641823E-03 1 0 10 #OK Aksz
par gal545_A_100 nui 0.7412626E+01 .2 0 50 #OK AdustTT100
par gal545_A_143 nui 0.8981278E+01 9.41 0 50 #OK AdustTT143
par gal545_A_143_217 nui 0.1753156E+02 19.69 0 100 #OK AdustTT143x217
par gal545_A_217 nui 0.8202580E+02 82.4 0 100 #OK AdustTT217
par calib_100T nui 0.9978923E+00 .2 0 3 #OK c100
par calib_217T nui 0.9959253E+00 .2 0 3 #OK c217
fix cib_index nui -1.3 .2 -5 5 #OK ncib
gausssum2 A_sz ksz_norm 1.6 1. 9.5 3
gauss1 calib_100T 0.9990004 0.001
gauss1 calib_217T 0.99501 0.002
gauss1 gal545_A_100 7 2
gauss1 gal545_A_143 9 2
gauss1 gal545_A_143_217 21 8.5
gauss1 gal545_A_217 80 20
gauss1 A_planck 1 0.0025
#class setup
class non\ linear halofit
class k_pivot 0.05
class lensing yes
class sBBN\ file bbn/sBBN.dat
precisionFile=class_pre/hpjul2.pre
clikfile=planck_data/hi_l/plik/plik_dx11dr2_HM_v18_TT.clik
clikfile2=planck_data/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5c_Ap.clik
#fitter
nitermax=50000
set_stra=2
set_tol=0.00001
#output options
computecontour=false
profile=false
remove_cosmo_limits=false
remove_limits=false
class N_ncdm 1
fix m_ncdm cosmo 0.06
class N_eff 2.046
engine=pico
par omega_b cosm 0.2224156E-01 0.00027 0.017 0.027
par omega_cdm cosm 0.1197735E+00 0.0026 0.09 0.15
par 100*theta_s cosm 0.1040862E+01 0.6E-04 1.0170 1.0617
par tau_reio cosm 0.7814866E-01 0.13E-01 0.01016 0.18221
par log(10^10A_s) cosm 0.3090373E+01 0.025 2.7 3.5
par n_s cosm 0.9658314E+00 0.0070 0.9 1.1
par A_planck nui 0.1000296E+01 0.01 0.9 1.1 #ycal
par A_cib_217 nui 0.6662739E+02 1 0 200 #OK Acib217
par xi_sz_cib nui 0.4657732E-01 .01 0 1 #xi_sz_cib
par A_sz nui 0.7141324E+01 1 0 10 #OK Atsz143
par ps_A_100_100 nui 0.2518416E+03 10 0 400 #OK Aps100
par ps_A_143_143 nui 0.3916351E+02 1 0 400 #OK Aps143
par ps_A_143_217 nui 0.3360542E+02 1 0 400 #OK Aps143x217
par ps_A_217_217 nui 0.9783567E+02 1 0 400 #OK Aps217
par ksz_norm nui 0.8641823E-03 1 0 10 #OK Aksz
par gal545_A_100 nui 0.7412626E+01 .2 0 50 #OK AdustTT100
par gal545_A_143 nui 0.8981278E+01 9.41 0 50 #OK AdustTT143
par gal545_A_143_217 nui 0.1753156E+02 19.69 0 100 #OK AdustTT143x217
par gal545_A_217 nui 0.8202580E+02 82.4 0 100 #OK AdustTT217
par calib_100T nui 0.9978923E+00 .2 0 3 #OK c100
par calib_217T nui 0.9959253E+00 .2 0 3 #OK c217
fix cib_index nui -1.3 .2 -5 5 #OK ncib
gausssum2 A_sz ksz_norm 1.6 1. 9.5 3
gauss1 calib_100T 0.9990004 0.001
gauss1 calib_217T 0.99501 0.002
gauss1 gal545_A_100 7 2
gauss1 gal545_A_143 9 2
gauss1 gal545_A_143_217 21 8.5
gauss1 gal545_A_217 80 20
gauss1 A_planck 1 0.0025
clikfile=planck_data/hi_l/plik/plik_dx11dr2_HM_v18_TT.clik
clikfile2=planck_data/low_l/bflike/lowl_SMW_70_dx11d_2014_10_03_v5c_Ap.clik
#fitter
nitermax=50000
set_stra=2
set_tol=0.00001
#output options
computecontour=false
profile=false
remove_cosmo_limits=false
remove_limits=false
fix m_ncdm cosmo 0.06
fix N_nu cosmo 3.046
fix m_ncdm cosm 0.06
fix Alens cosm 1
fix YHe cosm 0.248
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment