Add test MBH Graham23

#!/usr/bin/env python3

# Author: Akke Viitanen

# Email: akke.viitanen@helsinki.fi

# Date: 2025-12-14 20:59:01

"""

Test MBH

"""

from unittest import TestCase

from lsst_inaf_agile import mbh

class TestMBH(TestCase):

    def test_graham2023(self):

        import matplotlib.pyplot as plt

        import numpy as np

        for row in mbh.GRAHAM2023_ROWS:

            log_mstar = np.linspace(9, 12)

            log_mbh = mbh.get_log_mbh_graham2023(log_mstar, which=row)

            plt.plot(log_mstar, log_mbh, label=row)

        plt.xlabel(r"$\log M_\mathrm{star}$")

        plt.ylabel(r"$\log M_\mathrm{BH}$")

        log_mbh_kh13 = mbh.get_log_mbh_kormendy_ho2013(log_mstar)

        plt.plot(log_mstar, log_mbh_kh13, color="k", lw=2, label="KH13")

        log_mbh_kh13 = mbh.get_log_mbh_kormendy_ho2013(log_mstar, from_rv15=True)

        plt.plot(log_mstar, log_mbh_kh13, color="k", lw=2, label="KH13 (RV15)")

        plt.legend(loc="upper left")

        plt.show()