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
LEPHAREDIR is being set to the default cache directory:
/home/docs/.cache/lephare/data
More than 1Gb may be written there.
LEPHAREWORK is being set to the default cache directory:
/home/docs/.cache/lephare/work
Default work cache is already linked.
This is linked to the run directory:
/home/docs/.cache/lephare/runs/20260506T145844
[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.
Checking/downloading 445 files...
445 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 0.3844 0.3908 0.0538 0.3150 -0.3891 -20.6345 6.0327 0 0.3815 1.0000
u_new.pb 2 0.3690 0.3750 0.0456 0.6195 -0.2745 -20.8527 6.3135 0 0.3668 1.0000
gHSC.pb 3 0.4851 0.4760 0.1194 -0.0860 -0.2458 -20.7272 5.0764 0 0.4780 1.0000
rHSC.pb 4 0.6241 0.6142 0.1539 0.1466 0.2580 -21.5143 4.6523 0 0.6178 1.0000
iHSC.pb 5 0.7716 0.7637 0.1476 0.3942 0.6138 -22.2286 4.5323 0 0.7666 1.0000
zHSC.pb 6 0.8915 0.8907 0.0768 0.5169 0.7625 -22.6733 4.5147 0 0.8903 1.0000
yHSC.pb 7 0.9801 0.9771 0.0797 0.5534 0.7763 -22.9145 4.5081 0 0.9782 1.0000
Y.lowres 8 1.0222 1.0196 0.0919 0.6043 0.8180 -23.0574 4.5130 0 1.0206 1.0000
J.lowres 9 1.2555 1.2481 0.1712 0.9228 -99.0000 -23.8194 4.5638 0 1.2514 1.0000
H.lowres 10 1.6497 1.6352 0.2893 1.3701 -99.0000 -24.8565 4.7045 0 1.6409 1.0000
K.lowres 11 2.1577 2.1435 0.2926 1.8335 -99.0000 -25.9057 5.1316 0 2.1502 1.0000
IB427.lowres 12 0.4264 0.4256 0.0207 -0.1446 -0.4942 -20.4117 5.5152 0 0.4262 1.0000
IB464.lowres 13 0.4636 0.4633 0.0218 -0.1520 -0.3463 -20.5860 5.0658 0 0.4634 1.0000
IB484.lowres 14 0.4851 0.4846 0.0228 -0.0241 -0.1770 -20.8122 4.9880 0 0.4849 1.0000
IB505.lowres 15 0.5064 0.5061 0.0231 -0.0656 -0.1366 -20.8639 4.9423 0 0.5061 1.0000
IB527.lowres 16 0.5262 0.5259 0.0242 -0.0260 -0.0464 -20.9871 4.8937 0 0.5260 1.0000
IB574.lowres 17 0.5766 0.5762 0.0272 0.0657 0.1377 -21.2773 4.7042 0 0.5763 1.0000
IB624.lowres 18 0.6234 0.6230 0.0301 0.1527 0.2768 -21.5339 4.6386 0 0.6232 1.0000
IB679.lowres 19 0.6783 0.6779 0.0336 0.2542 0.4288 -21.8183 4.5709 0 0.6779 1.0000
IB709.lowres 20 0.7075 0.7071 0.0316 0.2982 0.4968 -21.9541 4.5558 0 0.7072 1.0000
IB738.lowres 21 0.7363 0.7358 0.0323 0.3460 0.5577 -22.0886 4.5449 0 0.7360 1.0000
IB767.lowres 22 0.7687 0.7681 0.0364 0.3992 0.6164 -22.2351 4.5243 0 0.7683 1.0000
IB827.lowres 23 0.8246 0.8241 0.0344 0.4891 0.7300 -22.4777 4.5161 0 0.8243 1.0000
NB711.lowres 24 0.7120 0.7119 0.0073 0.3072 0.5084 -21.9774 4.5542 0 0.7120 1.0000
NB816.lowres 25 0.8150 0.8149 0.0120 0.4713 0.7098 -22.4349 4.5154 0 0.8149 1.0000
NB118.lowres 26 1.1909 1.1909 0.0112 0.8376 -99.0000 -23.6250 4.5554 0 1.1909 1.0000
irac_ch1.lowres 27 3.5763 3.5264 0.7411 2.7951 -99.0000 -27.9585 6.0679 1 3.5634 1.0036
irac_ch2.lowres 28 4.5290 4.4609 1.0105 3.2634 -99.0000 -28.9384 6.5680 1 4.5111 1.0040
irac_ch3.lowres 29 5.7873 5.6765 1.3509 3.7537 -99.0000 -29.9581 7.0472 1 5.7592 1.0050
irac_ch4.lowres 30 8.0442 7.7033 2.8394 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 SED in the list 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
# CREATION_DATE Wed May 6 15:09:22 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 SED in the list 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,
# LIB_ASCII NO
# CREATION_DATE Wed May 6 15:09:22 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 SED in the list 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 # LIB_ASCII NO
# CREATION_DATE Wed May 6 15:09:28 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 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
# 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.200000 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: 13
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: 16
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: 13
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: 16
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: 13
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: 16
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
Read filt
# NAME IDENT Lbda_mean Lbeff(Vega) FWHM AB-cor VEGA CALIB Fac_corr
u_cfht.lowres 1 0.3844 0.3908 0.0538 0.3150 -20.6300 0 1.0000
u_new.pb 2 0.3690 0.3750 0.0456 0.6195 -20.8500 0 1.0000
gHSC.pb 3 0.4851 0.4760 0.1194 -0.0860 -20.7300 0 1.0000
rHSC.pb 4 0.6241 0.6142 0.1539 0.1466 -21.5100 0 1.0000
iHSC.pb 5 0.7716 0.7637 0.1476 0.3942 -22.2300 0 1.0000
zHSC.pb 6 0.8915 0.8907 0.0768 0.5169 -22.6700 0 1.0000
yHSC.pb 7 0.9801 0.9771 0.0797 0.5534 -22.9100 0 1.0000
Y.lowres 8 1.0220 1.0200 0.0919 0.6043 -23.0600 0 1.0000
J.lowres 9 1.2550 1.2480 0.1712 0.9228 -23.8200 0 1.0000
H.lowres 10 1.6500 1.6350 0.2893 1.3700 -24.8600 0 1.0000
K.lowres 11 2.1580 2.1430 0.2926 1.8330 -25.9100 0 1.0000
IB427.lowres 12 0.4264 0.4256 0.0207 -0.1446 -20.4100 0 1.0000
IB464.lowres 13 0.4636 0.4633 0.0218 -0.1520 -20.5900 0 1.0000
IB484.lowres 14 0.4851 0.4846 0.0228 -0.0241 -20.8100 0 1.0000
IB505.lowres 15 0.5064 0.5061 0.0231 -0.0656 -20.8600 0 1.0000
IB527.lowres 16 0.5262 0.5259 0.0242 -0.0260 -20.9900 0 1.0000
IB574.lowres 17 0.5766 0.5762 0.0272 0.0657 -21.2800 0 1.0000
IB624.lowres 18 0.6234 0.6230 0.0301 0.1527 -21.5300 0 1.0000
IB679.lowres 19 0.6783 0.6779 0.0336 0.2542 -21.8200 0 1.0000
IB709.lowres 20 0.7075 0.7071 0.0316 0.2982 -21.9500 0 1.0000
IB738.lowres 21 0.7363 0.7358 0.0323 0.3460 -22.0900 0 1.0000
IB767.lowres 22 0.7687 0.7681 0.0364 0.3992 -22.2400 0 1.0000
IB827.lowres 23 0.8246 0.8241 0.0344 0.4891 -22.4800 0 1.0000
NB711.lowres 24 0.7120 0.7119 0.0073 0.3072 -21.9800 0 1.0000
NB816.lowres 25 0.8150 0.8149 0.0120 0.4713 -22.4300 0 1.0000
NB118.lowres 26 1.1910 1.1910 0.0112 0.8376 -23.6200 0 1.0000
irac_ch1.lowres 27 3.5760 3.5260 0.7411 2.7950 -27.9600 1 1.0040
irac_ch2.lowres 28 4.5290 4.4610 1.0100 3.2630 -28.9400 1 1.0040
irac_ch3.lowres 29 5.7870 5.6760 1.3510 3.7540 -29.9600 1 1.0050
irac_ch4.lowres 30 8.0440 7.7030 2.8390 4.3960 -31.3000 1 1.0110
AUTO_ADAPT is set to YES. Computing offsets.
# 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
Source 95 // Band 18 removed to improve the chi2, with old and new chi2 1411.2432 968.7137
Source 95 // Band 14 removed to improve the chi2, with old and new chi2 1411.2432 637.8812
[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.829661886141765 | 4.793146933340135 | 4.826842617746427 | 4.810818628430521 | 4.735791721885345 | 4.87889777294346 | 78.14754768561784 | 29 | 1.0 | 0.0 | -99.9 | -99.0 | -99 | -99.0 | 4.76119270990018 | 109.61132233693942 | 29 | 204 | 149.7851058075375 | 1000.0 .. 1000.0 | 1000.0 .. 1000.0 | 966.251752740498 .. -2.2157305413111814 | -20.527785822720272 .. -20.097549849703952 | -1.0 .. -1.0 | -1.0 .. -1.0 | 236.81092817053948 | 27 | 0 | 4.839 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -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.572534241351857 | 9.184826330566931 | 8.1504144856926 | 3004.0 | 0.0 .. 5.741061854427692e-66 | |
| 2 | 1.083564531470666 | 1.094135192245143 | 1.1029630185323431 | 1.0999882512657646 | 1.0319882496041182 | 1.1679962416251295 | 93.46554905827557 | 15 | 1.0 | 0.0 | -99.9 | -99.0 | -99 | -99.0 | 1.0788389939758893 | 84.17796598747087 | 5 | 136 | 609.0491574281652 | 24.823528777031047 .. 1000.0 | 0.0964327877384378 .. 1000.0 | 0.823764788572272 .. 0.07419490624515242 | -20.654443788703727 .. -22.848686032630905 | 0.3119043586528303 .. 2.5621733087906122 | 4.0 .. 9.0 | 305.7656261849724 | 27 | 0 | 1.166 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -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.0409168284687 | 9.754650757625313 | 9.269802386023784 | 11004.0 | 2.3922346129588668e-163 .. 0.0 | |
| 3 | 0.17846490772411044 | 0.19827523183695886 | 0.2006863116415242 | 0.1999999999999743 | 0.13199999928471856 | 0.26800000071524754 | 166.77846176528413 | 6 | 1.0 | 0.0 | -99.9 | -99.0 | -99 | -99.0 | 0.14726583473678323 | 225.10533019688646 | 2 | 139 | 1153.7975332198678 | 20.455816150467925 .. 1000.0 | 0.030183067850864845 .. 1000.0 | 0.7314680606005165 .. -0.4008057083232578 | -20.194770840381747 .. -20.773707996980534 | 0.3861506934648933 .. 2.5428473843025685 | 2.0 .. 5.0 | 214.2813439359289 | 28 | 0 | 0.1649 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -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.023789314488688 | 9.751451019864028 | 9.245709899042403 | 11004.0 | 2.4395731843499614e-301 .. 0.0 | |
| 4 | 5.522630484838172 | 5.451015540388572 | 5.704357759502074 | 5.537830801543679 | 5.407604870133837 | 5.681757040439759 | 27.819407800209085 | 28 | 1.0 | 0.0 | -99.9 | -99.0 | -99 | -99.0 | 5.470470735043707 | 39.96775754447567 | 29 | 196 | 62.86395512171592 | 28.303047821698073 .. 1000.0 | 1.9744363436325476 .. 1000.0 | 965.8996691960407 .. -2.3434160669782074 | -21.053604143741875 .. -20.847790592306136 | -1.0 .. -1.0 | -1.0 .. -1.0 | 386.36301096697537 | 28 | 0 | 5.742 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -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.716167844028593 | 9.455481676185013 | 8.45774081482083 | 3004.0 | 0.0 .. 4.619342515518757e-32 | |
| 5 | 0.9745359630199829 | 0.899251207791433 | 1.0336485645860465 | 0.9603147412975389 | 0.8409376217859641 | 1.0850663557708564 | 43.396413985529065 | 24 | 1.0 | 0.0 | 2.7 | 25.0 | 25 | 0.0 | 0.6270940501343955 | 53.52160813645468 | 4 | 136 | 108.46789778535701 | 28.54446872020207 .. 1000.0 | 2.1272798742731407 .. 1000.0 | -0.15526096089426908 .. 0.4022187058951644 | -16.66001548110589 .. -17.653472268204332 | 0.031133529486414346 .. -1.0 | 4.0 .. -1.0 | 8.42653470979814 | 28 | 0 | 5.745 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -999.0 | -99.9 | -99.9 | -99.9 | -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.889236481214347 | 7.943940168363852 | 7.121132223822482 | 3003.0 | 8.33409720172609e-22 .. 9.696614608841358e-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 0x7ed33efaadd0>
[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 0x7ed33ee7bb80>]
[ ]: