lsst_inaf_agile.kelly2013

Implement Kelly+2013 BHMF.

Attributes

REDSHIFTS

Functions

`get_bhmf`(→ tuple[float, float, float, float])	Return Kelly & Shen 2013 BHMF at the given redshift z.
`_get_figax`()	Return Fig. 4 figure and axes.
`_plot_bhmf`(fig, axes, redshifts)	Plot Fig. 4 style BHMF given axes and redshifts.

Module Contents

REDSHIFTS = (0.4, 0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 1.8, 2.15, 2.65, 3.2, 3.75, 4.25, 4.75)[source]

get_bhmf(z: float, filename: str = 'data/kelly2013/table1.dat') → tuple[float, float, float, float][source]

Return Kelly & Shen 2013 BHMF at the given redshift z.

The BHMF is based on Table 1 of Kelly & Shen 2013. See ‘See Also’ below. The closest redshift bin to the input redshift is returned without any interpolation.

Parameters:

z (float) – Redshift of the black hole mass function.
filename (str) – Path to the machine-readable Kelly & Shen 2013 table 1.

Returns:

log_mbh, log_phi, log_phi_lo log_phi_hi – A 4-tuple containing the black hole mass function measurement.

Return type:

tuple[float, float, float, float]

Raises:

ValueError – If redshift is not within the Kelly & Shen 2013 range (z!=0-5).