Unverified Commit 0034204e authored by Akke Viitanen's avatar Akke Viitanen
Browse files

add Suh hiz mbh

parent fbc4cc96

Makefile

+6 −0
Original line number Diff line number Diff line
@@ -226,6 +226,12 @@ fig/plambda_qu_sf_ctn_ctk_20251025.pdf: src/scripts/plots/plot_plambda_qu_sf_ctn 
fig/mbh_mstar.pdf: src/scripts/plots/plot_mbh_mstar.py

	python3 $< $@



fig/mbh_mstar_test_z5-5.5.pdf: src/scripts/plots/plot_mbh_mstar.py

	python3 $< $@ --redshifts 0.0 4.0 5.0 5.1 5.2 5.3 5.4 5.45 --zfmt "%.2f" --vmax=6.0



overleaf/fig/mbh_mstar_20260225.pdf: src/scripts/plots/plot_mbh_mstar.py

	python3 $< $@



fig/mbh_mstar_with_zou2024.pdf: src/scripts/plots/plot_mbh_mstar.py

	python3 $< $@ --plot_zou2024

src/lsst_inaf_agile/mbh.py

+16 −3
Original line number Diff line number Diff line
@@ -32,8 +32,21 @@ def get_log_mbh_pacucci2023(log_mstar): 
    return -2.43 + 1.06 * log_mstar





def get_log_mbh_suh2020(log_mstar, alpha=1.64, beta=10.29, dalpha=0.07, dbeta=0.04, add_error=False):

def get_log_mbh_suh2020(

    log_mstar,

    alpha=1.64,

    beta=10.29,

    dalpha=0.07,

    dbeta=0.04,

    add_error=False,

    is_hiz=False,

):

    """Return the first equation of Suh+2020 Sec. 5."""



    # NOTE: derived in ecogia 20260302

    if is_hiz:

        return 0.7804437564499483 * log_mstar + -0.3367970928572994



    if add_error:

        alpha += np.random.normal(scale=dalpha)

        beta += np.random.normal(scale=dbeta)
@@ -201,7 +214,7 @@ FUN_LOG_MBH_CONTINUITY_NEW2 = {} 




def get_log_mbh_continuity_new2(

    log_mstar, z, filename="data/AGILE_Mstar_Mbh/DECODEmeans_Mstar_Mbh_Zou2024cc.dat"

    log_mstar, z, filename="data/AGILE_Mstar_Mbh/DECODEmeans_Mstar_Mbh_Zou2024cc.dat", *args, **kwargs

):

    """

    Return logMBH according to the continuity equation with Zou+24 p(lambda).
@@ -228,7 +241,7 @@ def get_log_mbh_continuity_new2( 


        from scipy.interpolate import RegularGridInterpolator



        FUN_LOG_MBH_CONTINUITY_NEW2[filename] = RegularGridInterpolator((mvec, zvec), MBH)

        FUN_LOG_MBH_CONTINUITY_NEW2[filename] = RegularGridInterpolator((mvec, zvec), MBH, *args, **kwargs)



    # clip

    log_mstar = np.atleast_1d(log_mstar)