Commit 3baa8f2f authored by Plaszczynski Stephane's avatar Plaszczynski Stephane

check minimization and thetas pb

parent e6aa494a
......@@ -53,8 +53,6 @@ macro_append CAMEL_linkopts " -L$(CFITSIODIR)/lib -lcfitsio "
#application FitClassFromPico -group=exec -s=../src/camel exec/FitClassFromPico.cc
#CLASS FastPk access
library class_extra -no_share $(CAMELROOT)/src/class_extra/*.c
macro_append CAMEL_linkopts " -lclass_extra "
#CAMEL LIBRARIES
......@@ -112,9 +110,10 @@ macro_append testHiLLiPOP_dependencies " MinuitFit "
macro_append test_jla_dependencies " MinuitFit "
macro_append testPico_dependencies " MinuitFit "
#SP 7/6/19:
#Pk signature in CLASS changed in v2.7 : need to propagate a macro
# defined in camel_setup
macro_append cppflags "" CLASSV_ABOVE_2_7 " -DCLASSV_ABOVE_2_7"
macro_append cflags "" CLASSV_ABOVE_2_7 " -DCLASSV_ABOVE_2_7"
library class_extra -no_share $(CAMELROOT)/src/class_extra/*.c
......@@ -4,7 +4,6 @@
/* stuf for class below v2.7*/
#ifndef CLASSV_ABOVE_2_7
int camel_spectra_fast_pk_at_kvec_and_zvec(
struct background * pba,
struct spectra * psp,
......
###############################################################
#Boltzman Engine
engine=class
###############################################################
#Cosmological model
par H0 cosm 68 5 60 80
par omega_b cosm 0.02224 0.00027 0.017 0.027
par omega_cdm cosm 0.1192 0.0026 0.09 0.15
par tau_reio cosm 0.07 0.13E-01 0.01 0.20
par log(10^10A_s) cosm 3.07 0.025 2.7 3.5
par n_s cosm 0.96 0.0070 0.9 1.1
###############################################################
gauss1 tau_reio 0.058 0.012
###############################################################
#LIKELIHOODS
#-----------------------------
#CMB high-l (HILLIPOP)
Hillipop=HiLLiPOP/DX11dHM_superExt_CO_TT.lik
par A_planck nui 1. 0.001 0.9 1.1
par c0 nui 0. 0.001 -0.05 0.05
par c1 nui 0. 0.001 -0.05 0.05
fix c2 nui 0. 0.001 -0.05 0.05
par c3 nui 0. 0.001 -0.05 0.05
par c4 nui 0.004 0.001 -0.05 0.05
par c5 nui 0.004 0.001 -0.05 0.05
par Aradio nui 1. 1E-01 0.0 2
par Adusty nui 1. 1E-01 0.0 2
par Asz nui 1. 0.1 0.0 10
par Acib nui 1. 0.1 0.0 10
par AdustTT nui 1. 0.1 0.0 2
par Aksz nui 0. 1.0 0.0 10
par Aszxcib nui 0. 1.0 0.0 10
#PRIORS
gauss1 AdustTT 1 0.2
gauss1 Asz 1 0.2
gauss1 Acib 1 0.2
gauss1 Aradio 1 0.2
gauss1 Adusty 1 0.2
gauss1 A_planck 1 0.0025
gauss1 c0 0 2e-3
gauss1 c1 0 2e-3
gauss1 c3 0 2e-3
gauss1 c4 0.002 2e-3
gauss1 c5 0.002 2e-3
#-----------------------------
###############################################################
###############################################################
#CLASS setup
class non\ linear halofit
class N_ncdm 1
class m_ncdm 0.06
class N_eff 2.046
class k_pivot 0.05
class lensing yes
class sBBN\ file bbn/sBBN.dat
precisionFile=class_pre/hpjul3.pre
###############################################################
###############################################################
#fitter
nitermax=50000
set_stra=2
set_tol=0.00001
#output options
remove_cosmo_limits=false
doHesse=true
###############################################################
###############################################################
#Boltzman Engine
engine=class
###############################################################
#Cosmological model
par omega_b cosm 0.02224 0.00027 0.017 0.027
par omega_cdm cosm 0.1192 0.0026 0.09 0.15
par 100*theta_s cosm 1.0418 0.6E-04 1.03 1.05
par tau_reio cosm 0.07 0.13E-01 0.01 0.20
par log(10^10A_s) cosm 3.07 0.025 2.7 3.5
par n_s cosm 0.96 0.0070 0.9 1.1
###############################################################
gauss1 tau_reio 0.058 0.012
###############################################################
#LIKELIHOODS
#-----------------------------
#CMB high-l (HILLIPOP)
Hillipop=HiLLiPOP/DX11dHM_superExt_CO_TT.lik
par A_planck nui 1. 0.001 0.9 1.1
par c0 nui 0. 0.001 -0.05 0.05
par c1 nui 0. 0.001 -0.05 0.05
fix c2 nui 0. 0.001 -0.05 0.05
par c3 nui 0. 0.001 -0.05 0.05
par c4 nui 0.004 0.001 -0.05 0.05
par c5 nui 0.004 0.001 -0.05 0.05
par Aradio nui 1. 1E-01 0.0 2
par Adusty nui 1. 1E-01 0.0 2
par Asz nui 1. 0.1 0.0 10
par Acib nui 1. 0.1 0.0 10
par AdustTT nui 1. 0.1 0.0 2
par Aksz nui 0. 1.0 0.0 10
par Aszxcib nui 0. 1.0 0.0 10
#PRIORS
gauss1 AdustTT 1 0.2
gauss1 Asz 1 0.2
gauss1 Acib 1 0.2
gauss1 Aradio 1 0.2
gauss1 Adusty 1 0.2
gauss1 A_planck 1 0.0025
gauss1 c0 0 2e-3
gauss1 c1 0 2e-3
gauss1 c3 0 2e-3
gauss1 c4 0.002 2e-3
gauss1 c5 0.002 2e-3
#-----------------------------
###############################################################
###############################################################
#CLASS setup
class non\ linear halofit
class N_ncdm 1
class m_ncdm 0.06
class N_eff 2.046
class k_pivot 0.05
class lensing yes
class sBBN\ file bbn/sBBN.dat
precisionFile=class_pre/hpjul3.pre
###############################################################
###############################################################
#fitter
nitermax=50000
set_stra=2
set_tol=0.00001
#output options
remove_cosmo_limits=false
doHesse=true
###############################################################
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