Minimal full run
We have developed some high level functions which allow the various stages of running LePHARE to be abstracted via a minimal interface.
In this notebook we provide a simple example of running this minimal set up with a table that is already in an acceptable format.
We use a very sparese redhsift grid to allow quick computation. In order to use the results scientically a higher redshift density should be used.
The notebook can be downloaded here.
[1]:
import lephare as lp
from astropy.table import Table
import numpy as np
from matplotlib import pylab as plt
%matplotlib inline
User defined LEPHAREDIR is set. Code runs depend on all required
auxiliary data being present at /home/docs/.cache/lephare/data.
User defined LEPHAREWORK is set. All intermediate files will be written to:
/home/docs/.cache/lephare/work
OMP_NUM_THREADS set to 1 (max available CPUs: 2)
[2]:
# Open the config file
config = lp.default_cosmos_config.copy()
[3]:
# this line reduced the zgrid density from the default to make the notebook run faster.
# Comment this out for better science results
config["Z_STEP"] = "0.1,0.,7."
[4]:
# Download the required data
lp.data_retrieval.get_auxiliary_data(keymap=config, additional_files=["examples/COSMOS.in"])
Registry file downloaded and saved as data_registry.txt.
/home/docs/checkouts/readthedocs.org/user_builds/lephare/envs/latest/lib/python3.10/site-packages/lephare/data_retrieval.py:431: UserWarning: /home/docs/.cache/lephare/data/sed/STAR/STAR_MOD_ALL.list is present locally and will not be overwritten.
warnings.warn(
/home/docs/checkouts/readthedocs.org/user_builds/lephare/envs/latest/lib/python3.10/site-packages/lephare/data_retrieval.py:431: UserWarning: /home/docs/.cache/lephare/data/sed/GAL/COSMOS_SED/COSMOS_MOD.list is present locally and will not be overwritten.
warnings.warn(
/home/docs/checkouts/readthedocs.org/user_builds/lephare/envs/latest/lib/python3.10/site-packages/lephare/data_retrieval.py:431: UserWarning: /home/docs/.cache/lephare/data/sed/QSO/SALVATO09/AGN_MOD.list is present locally and will not be overwritten.
warnings.warn(
Checking/downloading 442 files...
442 completed.
All files downloaded successfully and are non-empty.
Checking/downloading 1 files...
1 completed.
All files downloaded successfully and are non-empty.
[5]:
# Run the prepare stage to process filters, build the SED libraries, and theoretical magnitudes
lp.prepare(config)
# NAME IDENT Lbda_mean Lbeff(Vega) FWHM AB-cor TG-cor VEGA M_sun(AB) CALIB Lb_eff Fac_corr
u_cfht.lowres 1 3844.1 3908.3 537.9 0.3150 -0.3891 -20.6345 6.0327 0 0.3815 1.0000
u_new.pb 2 3689.8 3749.5 456.0 0.6195 -0.2745 -20.8527 6.3135 0 0.3668 1.0000
gHSC.pb 3 4850.6 4760.2 1194.2 -0.0860 -0.2458 -20.7272 5.0764 0 0.4780 1.0000
rHSC.pb 4 6241.4 6142.0 1538.5 0.1466 0.2580 -21.5143 4.6523 0 0.6178 1.0000
iHSC.pb 5 7716.5 7636.7 1476.1 0.3942 0.6138 -22.2286 4.5323 0 0.7666 1.0000
zHSC.pb 6 8915.2 8906.7 768.0 0.5169 0.7625 -22.6733 4.5147 0 0.8903 1.0000
yHSC.pb 7 9801.2 9770.6 796.8 0.5534 0.7763 -22.9145 4.5081 0 0.9782 1.0000
Y.lowres 8 10222.2 10195.8 919.2 0.6043 0.8180 -23.0574 4.5130 0 1.0206 1.0000
J.lowres 9 12554.9 12480.6 1712.3 0.9228 -99.0000 -23.8194 4.5638 0 1.2514 1.0000
H.lowres 10 16497.4 16351.6 2892.7 1.3701 -99.0000 -24.8565 4.7045 0 1.6409 1.0000
K.lowres 11 21576.6 21434.8 2926.3 1.8335 -99.0000 -25.9057 5.1316 0 2.1502 1.0000
IB427.lowres 12 4263.9 4255.9 206.5 -0.1446 -0.4942 -20.4117 5.5152 0 0.4262 1.0000
IB464.lowres 13 4636.1 4633.3 218.1 -0.1520 -0.3463 -20.5860 5.0658 0 0.4634 1.0000
IB484.lowres 14 4850.7 4846.4 228.2 -0.0241 -0.1770 -20.8122 4.9880 0 0.4849 1.0000
IB505.lowres 15 5063.6 5060.6 230.9 -0.0656 -0.1366 -20.8639 4.9423 0 0.5061 1.0000
IB527.lowres 16 5262.3 5258.9 242.1 -0.0260 -0.0464 -20.9871 4.8937 0 0.5260 1.0000
IB574.lowres 17 5766.0 5762.0 272.1 0.0657 0.1377 -21.2773 4.7042 0 0.5763 1.0000
IB624.lowres 18 6234.5 6230.0 300.6 0.1527 0.2768 -21.5339 4.6386 0 0.6232 1.0000
IB679.lowres 19 6782.6 6778.6 336.5 0.2542 0.4288 -21.8183 4.5709 0 0.6779 1.0000
IB709.lowres 20 7075.0 7070.6 315.6 0.2982 0.4968 -21.9541 4.5558 0 0.7072 1.0000
IB738.lowres 21 7362.7 7358.5 323.5 0.3460 0.5577 -22.0886 4.5449 0 0.7360 1.0000
IB767.lowres 22 7686.6 7681.2 363.8 0.3992 0.6164 -22.2351 4.5243 0 0.7683 1.0000
IB827.lowres 23 8245.9 8241.0 343.7 0.4891 0.7300 -22.4777 4.5161 0 0.8243 1.0000
NB711.lowres 24 7119.9 7119.4 72.6 0.3072 0.5084 -21.9774 4.5542 0 0.7120 1.0000
NB816.lowres 25 8149.9 8149.2 119.8 0.4713 0.7098 -22.4349 4.5154 0 0.8149 1.0000
NB118.lowres 26 11909.2 11908.6 112.3 0.8376 -99.0000 -23.6250 4.5554 0 1.1909 1.0000
irac_ch1.lowres 27 35763.0 35263.7 7410.9 2.7951 -99.0000 -27.9585 6.0679 1 3.5634 1.0036
irac_ch2.lowres 28 45289.5 44608.6 10104.6 3.2634 -99.0000 -28.9384 6.5680 1 4.5111 1.0040
irac_ch3.lowres 29 57873.0 56764.7 13508.5 3.7537 -99.0000 -29.9581 7.0472 1 5.7592 1.0050
irac_ch4.lowres 30 80442.3 77033.1 28394.0 4.3959 -99.0000 -31.2962 7.6701 1 7.9590 1.0110
#######################################
# It s translating SEDs to binary library #
# with the following options :
# Config file :
# Library type : STAR
# STAR_SED :/home/docs/.cache/lephare/data/sed/STAR/STAR_MOD_ALL.list
# STAR_LIB :LIB_STAR
# STAR_LIB doc:/home/docs/.cache/lephare/work/lib_bin/LIB_STAR.doc
# STAR_FSCALE :0.0000
#######################################
Number of templates in the list 254
Number of SED in the list (including different ages) 254
WRONG NUMBER OF ARGUMENTS FOR OPTION MOD_EXTINC
We have 2 instead of 8
Use default value 0,0 for all filters
#######################################
# It s computing the SYNTHETIC MAGNITUDES #
# For Gal/QSO libraries with these OPTIONS #
# with the following options :
# Config file :
# Filter file : filter_cosmos
# Magnitude type : AB
# COSMOLOGY :70.0000,0.3000,0.7000
# STAR_LIB_IN :/home/docs/.cache/lephare/work/lib_bin/LIB_STAR(.doc & .bin)
# STAR_LIB_OUT :/home/docs/.cache/lephare/work/lib_mag/STAR_COSMOS(.doc & .bin)
# LIB_ASCII NO
# MW_GALAMETZ NO
# EXT_MW_CURVE CARDELLI
# CREATION_DATE Wed Oct 28 09:41:40 2026
#############################################
#######################################
# It s translating SEDs to binary library #
# with the following options :
# Config file :
# Library type : GAL
# GAL_SED :/home/docs/.cache/lephare/data/sed/GAL/COSMOS_SED/COSMOS_MOD.list
# GAL_LIB :LIB_CE
# GAL_LIB doc:/home/docs/.cache/lephare/work/lib_bin/LIB_CE.doc
# GAL_LIB phys:/home/docs/.cache/lephare/work/lib_bin/LIB_CE.phys
# SEL_AGE :none
# GAL_FSCALE :1.0000
# AGE_RANGE 0.0000 15000000000.0000
#######################################
Number of templates in the list 31
Number of SED in the list (including different ages) 31
WRONG NUMBER OF ARGUMENTS FOR OPTION MOD_EXTINC
We have 2 instead of 8
Use default value 0,0 for all filters
#######################################
# It s computing the SYNTHETIC MAGNITUDES #
# For Gal/QSO libraries with these OPTIONS #
# with the following options :
# Config file :
# Filter file : filter_cosmos
# Magnitude type : AB
# GAL_LIB_IN :/home/docs/.cache/lephare/work/lib_bin/LIB_CE(.doc & .bin)
# GAL_LIB_OUT :/home/docs/.cache/lephare/work/lib_mag/CE_COSMOS(.doc & .bin)
# Z_STEP :0.1000 0.0000 7.0000
# COSMOLOGY :70.0000,0.3000,0.7000
# EXTINC_LAW :SMC_prevot.dat SB_calzetti.dat SB_calzetti_bump1.dat SB_calzetti_bump2.dat
# MOD_EXTINC :0 0 0 0 0 0 0 0
# EB_V :0.0000 0.0500 0.1000 0.1500 0.2000 0.2500 0.3000 0.3500 0.4000 0.5000
# EM_LINES EMP_UV
# EM_DISPERSION 0.5000,0.7500,1.0000,1.5000,2.0000,
# MW_GALAMETZ NO
# EXT_MW_CURVE CARDELLI
# LIB_ASCII NO
# CREATION_DATE Wed Oct 28 09:41:40 2026
#############################################
#######################################
# It s translating SEDs to binary library #
# with the following options :
# Config file :
# Library type : QSO
# QSO_SED :/home/docs/.cache/lephare/data/sed/QSO/SALVATO09/AGN_MOD.list
# QSO_LIB :LIB_QSO
# QSO_LIB doc:/home/docs/.cache/lephare/work/lib_bin/LIB_QSO.doc
# QSO_FSCALE :1.0000
#######################################
Number of templates in the list 30
Number of SED in the list (including different ages) 30
WRONG NUMBER OF ARGUMENTS FOR OPTION MOD_EXTINC
We have 2 instead of 8
Use default value 0,0 for all filters
#######################################
# It s computing the SYNTHETIC MAGNITUDES #
# For Gal/QSO libraries with these OPTIONS #
# with the following options :
# Config file :
# Filter file : filter_cosmos
# Magnitude type : AB
# QSO_LIB_IN :/home/docs/.cache/lephare/work/lib_bin/LIB_QSO(.doc & .bin)
# QSO_LIB_OUT :/home/docs/.cache/lephare/work/lib_mag/QSO_COSMOS(.doc & .bin)
# Z_STEP :0.1000 0.0000 7.0000
# COSMOLOGY :70.0000,0.3000,0.7000
# EXTINC_LAW :SMC_prevot.dat SB_calzetti.dat SB_calzetti_bump1.dat SB_calzetti_bump2.dat
# MOD_EXTINC :0 0 0 0 0 0 0 0
# EB_V :0.0000 0.0500 0.1000 0.1500 0.2000 0.2500 0.3000 0.3500 0.4000 0.5000
# MW_GALAMETZ NO
# EXT_MW_CURVE CARDELLI
# LIB_ASCII NO
# CREATION_DATE Wed Oct 28 09:41:43 2026
#############################################
[6]:
# Load the COSMOS example table. Lets just look at the first 100 objects to keep the notebook fast
input_table = Table.read(f"{lp.LEPHAREDIR}/examples/COSMOS.in", format="ascii")[:100]
input_table[input_table.colnames[-1]] = input_table[input_table.colnames[-1]].astype(
str
) # Make the final comments a str
[7]:
# Calculate the photometric redshifts
output, photozlist = lp.process(config, input_table)
Using user columns from input table assuming they are in the standard MEME order.
Processing 100 objects with 30 bands
#######################################
# PHOTOMETRIC REDSHIFT with OPTIONS #
# Config file :
# CAT_IN : change_me_to_output_filename_required.ascii
# CAT_OUT : zphot.out
# CAT_LINES : 0 1000000000
# PARA_OUT : /home/docs/.cache/lephare/data/examples/output.para
# INP_TYPE : F
# CAT_FMT[0:MEME 1:MMEE] : 0
# CAT_MAG : AB
# ZPHOTLIB : CE_COSMOS STAR_COSMOS QSO_COSMOS
# APPLY_MW_EXTINCTION : NO
# FIR_LIB :
# FIR_LMIN : 7.000000
# FIR_CONT : -1.000000
# FIR_SCALE : -1.000000
# FIR_FREESCALE : YES
# FIR_SUBSTELLAR : NO
# ERR_SCALE : 0.020000 0.020000 0.020000 0.020000 0.020000 0.020000 0.020000 0.050000 0.050000 0.050000 0.050000 0.020000 0.020000 0.020000 0.020000 0.020000 0.020000 0.020000 0.020000 0.020000 0.020000 0.020000 0.020000 0.050000 0.050000 0.050000 0.050000 0.100000 0.2000# Offsets added to the modeled magnitudes (or substracted to the observed): 0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0
00 0.300000
# ERR_FACTOR : 1.500000
# GLB_CONTEXT : 0
# FORB_CONTEXT : -1
# DZ_WIN : 1.000000
# MIN_THRES : 0.020000
# MAG_ABS : -24.000000 -5.000000
# MAG_ABS_AGN : -30.000000 -10.000000
# MAG_REF : 3
# NZ_PRIOR : -1 -1
# Z_INTERP : YES
# Z_METHOD : BEST
# MABS_METHOD : 1
# MABS_CONTEXT : 33556478
# MABS_REF : 11
# AUTO_ADAPT : NO
# ADAPT_BAND : 5
# ADAPT_LIM : 1.500000 23.000000
# ADAPT_ZBIN : 0.010000 6.000000
# ZFIX : NO
# SPEC_OUT : NO
# CHI_OUT : NO
# PDZ_OUT : test
#######################################
Reading input librairies ...
Read lib
Number of keywords to be read in the doc: 15
Number of keywords read at the command line (excluding -c config): 0
Reading keywords from /home/docs/.cache/lephare/work/lib_mag/QSO_COSMOS.doc
Number of keywords read in the config file: 18
Keyword NUMBER_ROWS not provided
Keyword NUMBER_SED not provided
Keyword Z_FORM not provided
Reading library: /home/docs/.cache/lephare/work/lib_mag/QSO_COSMOS.bin
Done with the library reading with 2130 SED read.
Number of keywords to be read in the doc: 15
Number of keywords read at the command line (excluding -c config): 0
Reading keywords from /home/docs/.cache/lephare/work/lib_mag/STAR_COSMOS.doc
Number of keywords read in the config file: 18
Keyword NUMBER_ROWS not provided
Keyword NUMBER_SED not provided
Keyword Z_FORM not provided
Reading library: /home/docs/.cache/lephare/work/lib_mag/STAR_COSMOS.bin
Done with the library reading with 2384 SED read.
Number of keywords to be read in the doc: 15
Number of keywords read at the command line (excluding -c config): 0
Reading keywords from /home/docs/.cache/lephare/work/lib_mag/CE_COSMOS.doc
Number of keywords read in the config file: 18
Keyword NUMBER_ROWS not provided
Keyword NUMBER_SED not provided
Keyword Z_FORM not provided
Reading library: /home/docs/.cache/lephare/work/lib_mag/CE_COSMOS.bin
Done with the library reading with 13389 SED read.
Read lib out
Done reading the FIR lib
# NAME IDENT Lbda_mean Lbeff(Vega) FWHM AB-cor VEGA CALIB Fac_corr
u_cfht.lowres 1 3844.1 3908.3 537.9 0.3150 -20.6345 0 1.0000
u_new.pb 2 3689.8 3749.5 456.0 0.6195 -20.8527 0 1.0000
gHSC.pb 3 4850.6 4760.2 1194.2 -0.0860 -20.7272 0 1.0000
rHSC.pb 4 6241.4 6142.0 1538.5 0.1466 -21.5143 0 1.0000
iHSC.pb 5 7716.5 7636.7 1476.1 0.3942 -22.2286 0 1.0000
zHSC.pb 6 8915.2 8906.7 768.0 0.5169 -22.6733 0 1.0000
yHSC.pb 7 9801.2 9770.6 796.8 0.5534 -22.9145 0 1.0000
Y.lowres 8 10222.2 10195.8 919.2 0.6043 -23.0574 0 1.0000
J.lowres 9 12554.9 12480.6 1712.3 0.9228 -23.8194 0 1.0000
H.lowres 10 16497.4 16351.6 2892.7 1.3701 -24.8565 0 1.0000
K.lowres 11 21576.6 21434.8 2926.3 1.8335 -25.9057 0 1.0000
IB427.lowres 12 4263.9 4255.9 206.5 -0.1446 -20.4117 0 1.0000
IB464.lowres 13 4636.1 4633.3 218.1 -0.1520 -20.5860 0 1.0000
IB484.lowres 14 4850.7 4846.4 228.2 -0.0241 -20.8122 0 1.0000
IB505.lowres 15 5063.6 5060.6 230.9 -0.0656 -20.8639 0 1.0000
IB527.lowres 16 5262.3 5258.9 242.1 -0.0260 -20.9871 0 1.0000
IB574.lowres 17 5766.0 5762.0 272.1 0.0657 -21.2773 0 1.0000
IB624.lowres 18 6234.5 6230.0 300.6 0.1527 -21.5339 0 1.0000
IB679.lowres 19 6782.6 6778.6 336.5 0.2542 -21.8183 0 1.0000
IB709.lowres 20 7075.0 7070.6 315.6 0.2982 -21.9541 0 1.0000
IB738.lowres 21 7362.7 7358.5 323.5 0.3460 -22.0886 0 1.0000
IB767.lowres 22 7686.6 7681.2 363.8 0.3992 -22.2351 0 1.0000
IB827.lowres 23 8245.9 8241.0 343.7 0.4891 -22.4777 0 1.0000
NB711.lowres 24 7119.9 7119.4 72.6 0.3072 -21.9774 0 1.0000
NB816.lowres 25 8149.9 8149.2 119.8 0.4713 -22.4349 0 1.0000
NB118.lowres 26 11909.2 11908.6 112.3 0.8376 -23.6250 0 1.0000
irac_ch1.lowres 27 35763.0 35263.7 7410.9 2.7951 -27.9585 1 1.0036
irac_ch2.lowres 28 45289.5 44608.6 10104.6 3.2634 -28.9384 1 1.0040
irac_ch3.lowres 29 57873.0 56764.7 13508.5 3.7537 -29.9581 1 1.0050
irac_ch4.lowres 30 80442.3 77033.1 28394.0 4.3959 -31.2962 1 1.0110
Source 95 // Band 19 removed to improve the chi2, with old and new chi2 1349.9118 941.3387
Source 95 // Band 16 removed to improve the chi2, with old and new chi2 1349.9118 636.6527
[8]:
# Have a look at the output table
output[:5]
[8]:
Table length=5
| IDENT | Z_BEST | Z_BEST68_LOW | Z_BEST68_HIGH | Z_MED | Z_MED68_LOW | Z_MED68_HIGH | CHI_BEST | MOD_BEST | EXTLAW_BEST | EBV_BEST | Z_SEC | CHI_SEC | MOD_SEC | EBV_SEC | ZQ_BEST | CHI_QSO | MOD_QSO | MOD_STAR | CHI_STAR | MAG_OBS() | ERR_MAG_OBS() | K_COR() | MAG_ABS() | EMAG_ABS() | MABS_FILT() | SCALE_BEST | NBAND_USED | CONTEXT | ZSPEC | AGE_BEST | AGE_INF | AGE_MED | AGE_SUP | LDUST_BEST | LDUST_INF | LDUST_MED | LDUST_SUP | LUM_TIR_BEST | LUM_TIR_INF | LUM_TIR_MED | LUM_TIR_SUP | MASS_BEST | MASS_INF | MASS_MED | MASS_SUP | SFR_BEST | SFR_INF | SFR_MED | SFR_SUP | SSFR_BEST | SSFR_INF | SSFR_MED | SSFR_SUP | COL1_INF | COL1_MED | COL1_SUP | COL2_INF | COL2_MED | COL2_SUP | LUM_NUV_BEST | LUM_R_BEST | LUM_K_BEST | MAG_PRED() | STRING_INPUT | PDF_BAY_ZG() |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| str3 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | int64 | float64 | float64 | float64 | float64 | int64 | float64 | float64 | float64 | int64 | int64 | float64 | float64[30] | float64[30] | float64[30] | float64[30] | float64[30] | float64[30] | float64 | int64 | int64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64 | float64[0] | str7 | float64[71] |
| 1 | 4.827588053176267 | 4.792864162484217 | 4.8247035987094815 | 4.809353639529532 | 4.735213468920509 | 4.877432705568992 | 78.4360919474212 | 29 | 1.0 | 0.0 | -99.9 | 1000000000.0 | -99 | -99.0 | 4.761192709900175 | 109.61132233693975 | 29 | 204 | 149.7851058075376 | 1000.0 .. 1000.0 | 1000.0 .. 1000.0 | 15.395790170924517 .. -2.2239838725661887 | -20.582116047452615 .. -20.07959855410581 | -1.0 .. -1.0 | -1.0 .. -1.0 | 233.50331104337775 | 27 | 0 | 4.839 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -999.0 | -999.0 | -999.0 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | 9.565328333941267 | 9.178277583553516 | 8.144090864033764 | 3004.0 | 0.0 .. 5.288727950493209e-66 | |
| 2 | 1.088886232875882 | 1.0946185834175848 | 1.103424155117143 | 1.0999947373384336 | 1.031994716536923 | 1.1679983366622735 | 93.02041227499744 | 15 | 1.0 | 0.0 | -99.9 | 1000000000.0 | -99 | -99.0 | 1.078838993975888 | 84.17796598747135 | 5 | 136 | 609.0491574281658 | 24.823528777031047 .. 1000.0 | 0.0964327877384378 .. 1000.0 | 0.8272439792088108 .. 0.06739097021029379 | -20.677429250740463 .. -22.863215297164704 | 0.3347937957432432 .. 2.5620645127976367 | 4.0 .. 9.0 | 308.9564059450365 | 27 | 0 | 1.166 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -999.0 | -999.0 | -999.0 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | 9.044273125660595 | 9.758733581455974 | 9.274099866183231 | 11004.0 | 2.6647756887591386e-163 .. 0.0 | |
| 3 | 0.17846490772410858 | 0.19827523183695872 | 0.20068631164152417 | 0.1999999999999743 | 0.13199999928471856 | 0.26800000071524754 | 166.77846176528962 | 6 | 1.0 | 0.0 | -99.9 | 1000000000.0 | -99 | -99.0 | 0.14726583473678412 | 225.10533019688802 | 2 | 139 | 1153.7975332198687 | 20.455816150467925 .. 1000.0 | 0.030183067850864845 .. 1000.0 | 0.7314680606005094 .. -0.4008057083232538 | -20.19477084038172 .. -20.773707996980512 | 0.409020137439148 .. 2.542738489332491 | 2.0 .. 5.0 | 214.28134393592472 | 28 | 0 | 0.1649 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -999.0 | -999.0 | -999.0 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | 8.022664495154324 | 9.751027607000434 | 9.24549800983544 | 11004.0 | 2.439573184376587e-301 .. 0.0 | |
| 4 | 5.65234530640983 | 5.495987902917178 | 5.709314255608968 | 5.551861054559076 | 5.415850160317101 | 5.688834397248671 | 27.214339736500367 | 29 | 1.0 | 0.0 | -99.9 | 1000000000.0 | -99 | -99.0 | 5.470470735043736 | 39.96775754447561 | 29 | 196 | 62.86395512171589 | 28.303047821698073 .. 1000.0 | 1.9744363436325476 .. 1000.0 | 17.004864816041504 .. -2.3987749671381082 | -21.05603593341305 .. -20.593117782458563 | -1.0 .. -1.0 | -1.0 .. -1.0 | 374.8115308873582 | 28 | 0 | 5.742 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -999.0 | -999.0 | -999.0 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | 9.770848234033718 | 9.383797483645967 | 8.349610764126215 | 3004.0 | 0.0 .. 3.833635420080481e-32 | |
| 5 | 0.9762834846739157 | 0.8984330128306971 | 1.0382097047745622 | 0.963203048464267 | 0.8329732360868706 | 1.094240114721352 | 43.68947867421007 | 24 | 1.0 | 0.0 | 2.7 | 50.82152956860226 | 25 | 0.0 | 0.6270940501343951 | 53.521608136454546 | 4 | 136 | 108.46789778535693 | 28.54446872020207 .. 1000.0 | 2.1272798742731407 .. 1000.0 | -0.14490217829212163 .. 0.3995215285798217 | -16.664482956205358 .. -17.655057721704065 | 0.03997003948488498 .. -1.0 | 4.0 .. -1.0 | 8.439872870532358 | 28 | 0 | 5.745 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -999.0 | -999.0 | -999.0 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | -99.9 | 7.888820165090046 | 7.94419345157095 | 7.121595851693947 | 3003.0 | 8.23266813272922e-22 .. 8.777207959254923e-35 |
[9]:
# Plot the predicted redhsifts against the spectroscopic redshifts
plt.scatter(np.log10(1 + output["ZSPEC"]), np.log10(1 + output["Z_BEST"]))
[9]:
<matplotlib.collections.PathCollection at 0x738e9f4d7070>
[10]:
# Look at one example posterior redshift distribution
[11]:
plt.plot(photozlist[0].pdfmap[11].xaxis, photozlist[0].pdfmap[11].vPDF)
[11]:
[<matplotlib.lines.Line2D at 0x738e841defe0>]
[ ]: