LePhare

In addition to SED fitting with EAZY, we also include functionality to run SED fitting with LePhare, which is an older code based in Fortran.

import os

import astropy.units as u

os.environ["GALFIND_CONFIG_DIR"] = "/nvme/scratch/work/austind/GALFIND/testing"
os.environ["GALFIND_CONFIG_NAME"] = "test_galfind_config.ini"

%load_ext autoreload
%autoreload 2
from galfind.catalogues.Catalogue import Catalogue
from galfind.sed_fitting.LePhare import LePhare

# Load in a JOF data object
# survey = "JOF"
# version = "v11"
# instrument_names = ["NIRCam"]
# forced_phot_band = ["F277W", "F356W", "F444W"]
survey = "test"
version = "v0"
instrument_names = ["ACS_WFC", "NIRCam"]
forced_phot_band = ["F200W", "F444W"]
aper_diams = [0.32] * u.arcsec
min_flux_pc_err = 10.0

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,
    im_str = ["test"],
)

SED_fit_params = {"GAL_TEMPLATES": "BC03_Chabrier2003_Z(m42_m62)"}
LePhare_fitter = LePhare(SED_fit_params)
print(LePhare_fitter)

LePhare_fitter.compile_filters(cat.filterset)

for type in ["STAR", "QSO", "GAL"]:
    LePhare_fitter.compile_binary(type)
    LePhare_fitter.compile_templates(cat.filterset, type, save_suffix = "")