Catalogue selection methods

In this notebook, we will utilize the Galaxy class selection methods of galaxies stored in the Catalogue object to create boolean selection columns in the outputted fits table. These boolean columns can be loaded back in to new Catalogue objects directly from the fits catalogue upon load-in, handy for any subsequent plotting of cutouts, scaling relations, and UV LF / mass function generation. While the details on the available galaxy selection techniques are outlined in the Galaxy selection notebook, we will show the catalogue level implementation of the EPOCHS selection here. It is worth noting that all other galaxy selection options are also available on a catalogue scale.

[1]:

# imports
import astropy.units as u
from galfind import EAZY, Catalogue, EPOCHS_Selector
from galfind.Data import morgan_version_to_dir

WARNING:galfind:Could not change permissions of /raid/scratch/work/austind/GALFIND_WORK/Log_files/2024-12-11.log to 777.

Reading GALFIND config file from: /nvme/scratch/work/austind/GALFIND/galfind/../configs/galfind_config.ini

WARNING:galfind:Aperture corrections for VISTA not found in /nvme/scratch/work/austind/GALFIND/galfind/Aperture_corrections/VISTA_aper_corr.txt
WARNING:galfind:Aperture corrections for VISTA not found in /nvme/scratch/work/austind/GALFIND/galfind/Aperture_corrections/VISTA_aper_corr.txt

Failed to `import dust_attenuation`
Install from the repo with $ pip install git+https://github.com/karllark/dust_attenuation.git

[2]:

survey = "JOF"
version = "v11"
instrument_names = ["NIRCam"]
aper_diams = [0.32] * u.arcsec
forced_phot_band = ["F277W", "F356W", "F444W"]
min_flux_pc_err = 10.
SED_fit_params_arr = [
    {"templates": "fsps_larson", "lowz_zmax": 4.0},
    {"templates": "fsps_larson", "lowz_zmax": 6.0},
    {"templates": "fsps_larson", "lowz_zmax": None}
]

JOF_cat = Catalogue.pipeline(
    survey,
    version,
    instrument_names = instrument_names,
    version_to_dir_dict = morgan_version_to_dir,
    aper_diams = aper_diams,
    forced_phot_band = forced_phot_band,
    min_flux_pc_err = min_flux_pc_err
)
# load sextractor half-light radii
JOF_cat.load_sextractor_Re()

# load EAZY SED fitting results
for SED_fit_params in SED_fit_params_arr:
    EAZY_fitter = EAZY(SED_fit_params)
    EAZY_fitter(JOF_cat, aper_diams[0], load_PDFs = True, load_SEDs = True, update = True)

print(JOF_cat)

INFO:galfind:Loaded aper_diams=<Quantity [0.32] arcsec> for F277W+F356W+F444W
INFO:galfind:Combined mask for NIRCam/F277W+F356W+F444W already exists at /raid/scratch/work/austind/GALFIND_WORK/Masks/JOF/combined/JOF_F277W+F356W+F444W_auto.fits
WARNING: hdu= was not specified but multiple tables are present, reading in first available table (hdu=1) [astropy.io.fits.connect]
WARNING:astroquery:hdu= was not specified but multiple tables are present, reading in first available table (hdu=1)
WARNING:galfind:Aperture correction columns already in /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits
Calculating depths:   0%|          | 0/15 [00:00<?, ?it/s]
INFO:galfind:Calculated/loaded depths for JOF v11 NIRCam
INFO:galfind:Local depth columns already exist in /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits
INFO:galfind:Loaded 'has_data_mask' from /raid/scratch/work/austind/GALFIND_WORK/Masks/JOF/has_data_mask/JOF_MASTER_Sel-F277W+F356W+F444W_v11.h5
INFO:galfind:Making JOF v11 JOF_MASTER_Sel-F277W+F356W+F444W_v11 catalogue!
WARNING:galfind:cat_aper_diams not in kwargs.keys()=dict_keys(['ZP', 'min_flux_pc_err'])! Setting to aper_diams=<Quantity [0.32] arcsec>
WARNING:galfind:cat_aper_diams not in kwargs.keys()=dict_keys(['ZP', 'min_flux_pc_err'])! Setting to aper_diams=<Quantity [0.32] arcsec>
WARNING:galfind:cat_aper_diams not in kwargs.keys()=dict_keys([])! Setting to aper_diams=<Quantity [0.32] arcsec>
INFO:galfind:Made /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits catalogue!
INFO:galfind:Loaded FLUX_RADIUS from /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits saved as sex_Re for cat_band_properties[0].keys()=dict_keys(['F090W', 'F115W', 'F150W', 'F162M', 'F182M', 'F200W', 'F210M', 'F250M', 'F277W', 'F300M', 'F335M', 'F356W', 'F410M', 'F444W'])
INFO:galfind:Making .in file for EAZY_fsps_larson_zmax=4.0 SED fitting for JOF v11 NIRCam
INFO:galfind:Made .in file for EAZY_fsps_larson_zmax=4.0 SED fitting for JOF v11 NIRCam.

Running SED fitting took 0.1s

INFO:galfind:Loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam

Loading properties and associated errors took 1.4s

Constructing redshift PDFs: 100%|██████████| 16335/16335 [00:00<00:00, 57825.10it/s]
INFO:galfind:Finished loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam
INFO:galfind:Loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
Constructing SEDs: 100%|██████████| 16335/16335 [00:01<00:00, 12404.10it/s]
INFO:galfind:Finished loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
INFO:galfind:Updating SED results in galfind catalogue object
Updating galaxy SED results: 100%|██████████| 16335/16335 [00:00<00:00, 157758.30it/s]
INFO:galfind:Making .in file for EAZY_fsps_larson_zmax=6.0 SED fitting for JOF v11 NIRCam
INFO:galfind:Made .in file for EAZY_fsps_larson_zmax=6.0 SED fitting for JOF v11 NIRCam.

Running SED fitting took 0.1s

INFO:galfind:Loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam

Loading properties and associated errors took 1.3s

Constructing redshift PDFs: 100%|██████████| 16335/16335 [00:00<00:00, 58049.29it/s]
INFO:galfind:Finished loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam
INFO:galfind:Loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
Constructing SEDs: 100%|██████████| 16335/16335 [00:01<00:00, 8771.30it/s]
INFO:galfind:Finished loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
INFO:galfind:Updating SED results in galfind catalogue object
Updating galaxy SED results: 100%|██████████| 16335/16335 [00:00<00:00, 114869.76it/s]
INFO:galfind:Making .in file for EAZY_fsps_larson_zfree SED fitting for JOF v11 NIRCam
INFO:galfind:Made .in file for EAZY_fsps_larson_zfree SED fitting for JOF v11 NIRCam.

Running SED fitting took 0.1s

INFO:galfind:Loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam

Loading properties and associated errors took 1.3s

Constructing redshift PDFs: 100%|██████████| 16335/16335 [00:00<00:00, 45896.80it/s]
INFO:galfind:Finished loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam
INFO:galfind:Loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
Constructing SEDs: 100%|██████████| 16335/16335 [00:02<00:00, 7851.67it/s]
INFO:galfind:Finished loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
INFO:galfind:Updating SED results in galfind catalogue object
Updating galaxy SED results: 100%|██████████| 16335/16335 [00:00<00:00, 121496.96it/s]
WARNING: hdu= was not specified but multiple tables are present, reading in first available table (hdu=1) [astropy.io.fits.connect]
WARNING:astroquery:hdu= was not specified but multiple tables are present, reading in first available table (hdu=1)

****************************************
CATALOGUE(JOF,v11,NIRCam):
----------
CAT PATH = /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits
TOTAL GALAXIES = 16335
RA RANGE = [53.01070689 53.11059594] deg
DEC RANGE = [-27.91226173 -27.83206063] deg
----------
****************************************
MULTIPLE_FILTER
----------
FACILITY: JWST
INSTRUMENT: NIRCam
FILTERS: ['F090W', 'F115W', 'F150W', 'F162M', 'F182M', 'F200W', 'F210M', 'F250M', 'F277W', 'F300M', 'F335M', 'F356W', 'F410M', 'F444W']
****************************************

Now that we have loaded the blank catalogue, we will perform the selection. Since we are running on a Catalogue this time around instead of a Galaxy (as in Galaxy Selection, Example 4), we do not need to insert the catalogue filterset when instantiating EPOCHS_Selector as this information is already stored in the catalogue we are running.

[3]:

# perform EPOCHS selection
epochs_selector = EPOCHS_Selector(aper_diams[0], EAZY_fitter, allow_lowz = False, unmasked_instruments = "NIRCam")
EPOCHS_JOF_cat = epochs_selector(JOF_cat, return_copy = True)

Selecting bluewards_Lya_SNR<2.0_EAZY_fsps_larson_zfree_0.32as:   0%|          | 0/16335 [00:00<?, ?it/s]Selecting bluewards_Lya_SNR<2.0_EAZY_fsps_larson_zfree_0.32as: 100%|██████████| 16335/16335 [00:00<00:00, 82465.86it/s]
WARNING:galfind:property_name='bluewards_Lya_SNR<2.0_EAZY_fsps_larson_zfree_0.32as' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting redwards_Lya_SNR>5.0,5.0_widebands_EAZY_fsps_larson_zfree_0.32as: 100%|██████████| 16335/16335 [00:00<00:00, 27624.17it/s]
WARNING:galfind:property_name='redwards_Lya_SNR>5.0,5.0_widebands_EAZY_fsps_larson_zfree_0.32as' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting ALL_redwards_Lya_SNR>2.0_widebands_EAZY_fsps_larson_zfree_0.32as: 100%|██████████| 16335/16335 [00:00<00:00, 73698.42it/s]
WARNING:galfind:property_name='ALL_redwards_Lya_SNR>2.0_widebands_EAZY_fsps_larson_zfree_0.32as' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting red_chi_sq<3.0_EAZY_fsps_larson_zfree_0.32as: 100%|██████████| 16335/16335 [00:00<00:00, 90268.24it/s]
WARNING:galfind:property_name='red_chi_sq<3.0_EAZY_fsps_larson_zfree_0.32as' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting chi_sq_diff>4.0,dz>0.5_EAZY_fsps_larson_zfree_0.32as: 100%|██████████| 16335/16335 [00:00<00:00, 78391.61it/s]
WARNING:galfind:property_name='chi_sq_diff>4.0,dz>0.5_EAZY_fsps_larson_zfree_0.32as' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting zPDF>60%,|dz|/z<0.1_EAZY_fsps_larson_zfree_0.32as: 100%|██████████| 16335/16335 [00:00<00:00, 90060.12it/s]
WARNING:galfind:property_name='zPDF>60%,|dz|/z<0.1_EAZY_fsps_larson_zfree_0.32as' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F090W: 100%|██████████| 16335/16335 [00:00<00:00, 407806.60it/s]
WARNING:galfind:property_name='unmasked_F090W' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F115W: 100%|██████████| 16335/16335 [00:00<00:00, 428560.43it/s]
WARNING:galfind:property_name='unmasked_F115W' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F150W: 100%|██████████| 16335/16335 [00:00<00:00, 457834.09it/s]
WARNING:galfind:property_name='unmasked_F150W' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F162M: 100%|██████████| 16335/16335 [00:00<00:00, 465811.99it/s]
WARNING:galfind:property_name='unmasked_F162M' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F182M: 100%|██████████| 16335/16335 [00:00<00:00, 436605.97it/s]
WARNING:galfind:property_name='unmasked_F182M' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F200W: 100%|██████████| 16335/16335 [00:00<00:00, 416017.71it/s]
WARNING:galfind:property_name='unmasked_F200W' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F210M: 100%|██████████| 16335/16335 [00:00<00:00, 406425.29it/s]
WARNING:galfind:property_name='unmasked_F210M' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F250M: 100%|██████████| 16335/16335 [00:00<00:00, 350003.86it/s]
WARNING:galfind:property_name='unmasked_F250M' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F277W: 100%|██████████| 16335/16335 [00:00<00:00, 337517.15it/s]
WARNING:galfind:property_name='unmasked_F277W' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F300M: 100%|██████████| 16335/16335 [00:00<00:00, 339950.46it/s]
WARNING:galfind:property_name='unmasked_F300M' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F335M: 100%|██████████| 16335/16335 [00:00<00:00, 374540.56it/s]
WARNING:galfind:property_name='unmasked_F335M' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F356W: 100%|██████████| 16335/16335 [00:00<00:00, 345510.07it/s]
WARNING:galfind:property_name='unmasked_F356W' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F410M: 100%|██████████| 16335/16335 [00:00<00:00, 400964.20it/s]
WARNING:galfind:property_name='unmasked_F410M' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_F444W: 100%|██████████| 16335/16335 [00:00<00:00, 365014.52it/s]
WARNING:galfind:property_name='unmasked_F444W' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting unmasked_NIRCam: 100%|██████████| 16335/16335 [00:00<00:00, 474989.81it/s]
WARNING:galfind:property_name='unmasked_NIRCam' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting bluest_band_SNR<2.0_0.32as: 100%|██████████| 16335/16335 [00:00<00:00, 69916.41it/s]
WARNING:galfind:property_name='bluest_band_SNR<2.0_0.32as' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting sex_Re_F277W>45.0mas: 100%|██████████| 16335/16335 [00:00<00:00, 97083.08it/s]
WARNING:galfind:property_name='sex_Re_F277W>45.0mas' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting sex_Re_F356W>45.0mas: 100%|██████████| 16335/16335 [00:00<00:00, 94243.72it/s]
WARNING:galfind:property_name='sex_Re_F356W>45.0mas' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting sex_Re_F444W>45.0mas: 100%|██████████| 16335/16335 [00:00<00:00, 83504.21it/s]
WARNING:galfind:property_name='sex_Re_F444W>45.0mas' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting sex_Re_F277W+F356W+F444W>45.0mas: 100%|██████████| 16335/16335 [00:00<00:00, 349934.14it/s]
WARNING:galfind:property_name='sex_Re_F277W+F356W+F444W>45.0mas' already appended to hdu='SELECTION' .fits table, not overwriting!
Selecting EPOCHS_NIRCam_EAZY_fsps_larson_zfree_0.32as: 100%|██████████| 16335/16335 [00:00<00:00, 88523.67it/s]
WARNING:galfind:property_name='EPOCHS_NIRCam_EAZY_fsps_larson_zfree_0.32as' already appended to hdu='SELECTION' .fits table, not overwriting!
INFO:galfind:deepcopy(Catalogue)

Of course we see that returning a deep copy of the catalogue object takes longer than not. Let’s have a look at how this changes the Catalogue print statement.

[4]:

print(EPOCHS_JOF_cat)

****************************************
CATALOGUE(JOF,v11,NIRCam):
----------
CAT PATH = /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits
TOTAL GALAXIES = 16335
RA RANGE = [53.01070689 53.11059594] deg
DEC RANGE = [-27.91226173 -27.83206063] deg
----------
****************************************
MULTIPLE_FILTER
----------
FACILITY: JWST
INSTRUMENT: NIRCam
FILTERS: ['F090W', 'F115W', 'F150W', 'F162M', 'F182M', 'F200W', 'F210M', 'F250M', 'F277W', 'F300M', 'F335M', 'F356W', 'F410M', 'F444W']
****************************************

Fantastic! We now have a catalogue that has been cropped to only the EPOCHS sample.

Once the selection has been run on the Catalogue object, it is saved in the fits catalogue and is automatically loaded back in immediately when re-initializing the catalogue.

[5]:

JOF_cat_new = Catalogue.pipeline(
    survey,
    version,
    instrument_names = instrument_names,
    version_to_dir_dict = morgan_version_to_dir,
    aper_diams = aper_diams,
    forced_phot_band = forced_phot_band,
    min_flux_pc_err = min_flux_pc_err
)
# load sextractor half-light radii
JOF_cat_new.load_sextractor_Re()

# load EAZY SED fitting results
for SED_fit_params in SED_fit_params_arr:
    EAZY_fitter = EAZY(SED_fit_params)
    EAZY_fitter(JOF_cat_new, aper_diams[0], load_PDFs = True, load_SEDs = True, update = True)

INFO:galfind:Loaded aper_diams=<Quantity [0.32] arcsec> for F277W+F356W+F444W
INFO:galfind:Combined mask for NIRCam/F277W+F356W+F444W already exists at /raid/scratch/work/austind/GALFIND_WORK/Masks/JOF/combined/JOF_F277W+F356W+F444W_auto.fits
WARNING:galfind:Aperture correction columns already in /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits
Calculating depths:   0%|          | 0/15 [00:00<?, ?it/s]
INFO:galfind:Calculated/loaded depths for JOF v11 NIRCam
INFO:galfind:Local depth columns already exist in /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits
INFO:galfind:Loaded 'has_data_mask' from /raid/scratch/work/austind/GALFIND_WORK/Masks/JOF/has_data_mask/JOF_MASTER_Sel-F277W+F356W+F444W_v11.h5
INFO:galfind:Making JOF v11 JOF_MASTER_Sel-F277W+F356W+F444W_v11 catalogue!
WARNING:galfind:cat_aper_diams not in kwargs.keys()=dict_keys(['ZP', 'min_flux_pc_err'])! Setting to aper_diams=<Quantity [0.32] arcsec>
WARNING:galfind:cat_aper_diams not in kwargs.keys()=dict_keys(['ZP', 'min_flux_pc_err'])! Setting to aper_diams=<Quantity [0.32] arcsec>
WARNING:galfind:cat_aper_diams not in kwargs.keys()=dict_keys([])! Setting to aper_diams=<Quantity [0.32] arcsec>
INFO:galfind:Made /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits catalogue!
INFO:galfind:Loaded FLUX_RADIUS from /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits saved as sex_Re for cat_band_properties[0].keys()=dict_keys(['F090W', 'F115W', 'F150W', 'F162M', 'F182M', 'F200W', 'F210M', 'F250M', 'F277W', 'F300M', 'F335M', 'F356W', 'F410M', 'F444W'])
INFO:galfind:Making .in file for EAZY_fsps_larson_zmax=4.0 SED fitting for JOF v11 NIRCam
INFO:galfind:Made .in file for EAZY_fsps_larson_zmax=4.0 SED fitting for JOF v11 NIRCam.

Running SED fitting took 0.1s

INFO:galfind:Loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam

Loading properties and associated errors took 1.3s

Constructing redshift PDFs: 100%|██████████| 16335/16335 [00:00<00:00, 52702.52it/s]
INFO:galfind:Finished loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam
INFO:galfind:Loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
Constructing SEDs: 100%|██████████| 16335/16335 [00:01<00:00, 12606.34it/s]
INFO:galfind:Finished loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
INFO:galfind:Updating SED results in galfind catalogue object
Updating galaxy SED results: 100%|██████████| 16335/16335 [00:00<00:00, 109341.89it/s]
INFO:galfind:Making .in file for EAZY_fsps_larson_zmax=6.0 SED fitting for JOF v11 NIRCam
INFO:galfind:Made .in file for EAZY_fsps_larson_zmax=6.0 SED fitting for JOF v11 NIRCam.

Running SED fitting took 0.1s

INFO:galfind:Loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam

Loading properties and associated errors took 2.7s

Constructing redshift PDFs: 100%|██████████| 16335/16335 [00:00<00:00, 50266.21it/s]
INFO:galfind:Finished loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam
INFO:galfind:Loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
Constructing SEDs: 100%|██████████| 16335/16335 [00:01<00:00, 10613.19it/s]
INFO:galfind:Finished loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
INFO:galfind:Updating SED results in galfind catalogue object
Updating galaxy SED results: 100%|██████████| 16335/16335 [00:00<00:00, 109709.92it/s]
INFO:galfind:Making .in file for EAZY_fsps_larson_zfree SED fitting for JOF v11 NIRCam
INFO:galfind:Made .in file for EAZY_fsps_larson_zfree SED fitting for JOF v11 NIRCam.

Running SED fitting took 0.1s

INFO:galfind:Loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam

Loading properties and associated errors took 1.3s

Constructing redshift PDFs: 100%|██████████| 16335/16335 [00:00<00:00, 45106.84it/s]
INFO:galfind:Finished loading EAZY_fsps_larson property PDFs into JOF v11 NIRCam
INFO:galfind:Loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
Constructing SEDs: 100%|██████████| 16335/16335 [00:01<00:00, 10576.12it/s]
INFO:galfind:Finished loading EAZY_fsps_larson SEDs into JOF v11 NIRCam
INFO:galfind:Updating SED results in galfind catalogue object
Updating galaxy SED results: 100%|██████████| 16335/16335 [00:00<00:00, 107419.12it/s]

[6]:

if JOF_cat == JOF_cat_new:
    print("Catalogues are the same")
else:
    print("Catalogues are different")

Catalogues are different

Should you notice an error in, for instance, a custom selector, you can always delete the selection fits extension so that you don’t keep re-loading the dodgy results.

[7]:

JOF_cat_new.del_hdu(hdu = "SELECTION")

print(JOF_cat_new)

INFO:galfind:Writing table to /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits
INFO:galfind:Deleted hdu.upper()='SELECTION' from self.cat_path='/raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits'!
WARNING: hdu= was not specified but multiple tables are present, reading in first available table (hdu=1) [astropy.io.fits.connect]
WARNING:astroquery:hdu= was not specified but multiple tables are present, reading in first available table (hdu=1)

****************************************
CATALOGUE(JOF,v11,NIRCam):
----------
CAT PATH = /raid/scratch/work/austind/GALFIND_WORK/Catalogues/v11/NIRCam/JOF/(0.32)as/JOF_MASTER_Sel-F277W+F356W+F444W_v11.fits
TOTAL GALAXIES = 16335
RA RANGE = [53.01070689 53.11059594] deg
DEC RANGE = [-27.91226173 -27.83206063] deg
----------
****************************************
MULTIPLE_FILTER
----------
FACILITY: JWST
INSTRUMENT: NIRCam
FILTERS: ['F090W', 'F115W', 'F150W', 'F162M', 'F182M', 'F200W', 'F210M', 'F250M', 'F277W', 'F300M', 'F335M', 'F356W', 'F410M', 'F444W']
****************************************