path_git, path_datos_global = definir_path()
os.chdir(path_git)
sys.path.append('./Software/Funcionales/')
from funciones_sampleo import MCMC_sampler
from funciones_data import leer_data_pantheon, leer_data_cronometros, leer_data_BAO, leer_data_AGN
from funciones_alternativos import params_to_chi2
#ORDEN DE PRESENTACION DE LOS PARAMETROS: Mabs,omega_m,b,H_0,n
#%%
## Supernovas
os.chdir(path_git + '/Software/Estadística/Datos/Datos_pantheon/')
ds_SN = leer_data_pantheon('lcparam_full_long_zhel.txt')
# Cronómetros
os.chdir(path_git + '/Software/Estadística/Datos/')
ds_CC = leer_data_cronometros('datos_cronometros.txt')
# BAO
os.chdir(path_git + '/Software/Estadística/Datos/BAO/')
ds_BAO = []
archivos_BAO = [
'datos_BAO_da.txt', 'datos_BAO_dh.txt', 'datos_BAO_dm.txt',
'datos_BAO_dv.txt', 'datos_BAO_H.txt'
]
for i in range(5):
aux = leer_data_BAO(archivos_BAO[i])
ds_BAO.append(aux)
# AGN
os.chdir(path_git + '/Software/Estadística/Datos/Datos_AGN')
ds_AGN = leer_data_AGN('table3.dat')
M_true = -19.2
omega_m_true = 0.4
b_true = -0.5
H0_true = 73.48 #Unidades de (km/seg)/Mpc
n = 1
params_fijos = [H0_true,n]
#%%
#Datos de SN
os.chdir(path_git+'/Software/Estadística/Datos/Datos_pantheon/')
zcmb,zhel, Cinv, mb = leer_data_pantheon('lcparam_full_long_zhel.txt')
#Datos de crnómetros
os.chdir(path_git+'/Software/Estadística/Datos/')
z_data, H_data, dH = leer_data_cronometros('datos_cronometros.txt')
#%%
#Parametros a ajustar
nll = lambda theta: params_to_chi2_taylor(theta,params_fijos,zcmb, zhel, Cinv,
mb,z_data, H_data, dH,chi_riess=False)
initial = np.array([M_true,omega_m_true,b_true])
soln = minimize(nll, initial, options = {'eps': 0.01}, bounds =((-20,-18),(0.1,0.5),(-1, 1)))
M_ml, omega_m_ml, b_ml = soln.x
print(M_ml,omega_m_ml,b_ml)
os.chdir(path_git + '/Software/Estadística/Resultados_simulaciones')
np.savez('valores_medios_HS_CC+SN_3params_taylor', sol=soln.x)
