mrpy.extra.likelihoods.SampleLike

class mrpy.extra.likelihoods.SampleLike(logm, logHs, alpha, beta, lnA, log_mmin=None, rhom=83265000000.0)

A subclass of mrpy.core.MRP which adds the likelihood (and derivatives) of a model given a sample of masses.

The likelihoods in this class are true under a number of simplifications. Firstly, the effective volume is constant down to a threshold minimum mass. Secondly, the masses have no measurement uncertainty. While these simplifications are rather restrictive, nevertheless the quantities here are suprisingly useful, for instance with estimating parameters for simulations, or as an ingredient in a more sophisticated analysis.

Parameters:

logm : array_like

Vector of log10 masses.

logHs, alpha, beta, lnA : array_like

The parameters of the MRP.

log_mmin : array_like, optional

Log-10 truncation mass of the MRP. By default is set to the minimum mass in logm.

rhom : float

Mass density of the universe. Only used if the normalisation is set to Arhom.

Methods

__init__(logm, logHs, alpha, beta, lnA[, …])
dndlog10m([log]) Return the MRP in log10 space at `m’.
dndm([log]) Return the MRP at m.
ngtm([log]) The number density greater than mmin.
rho_gtm([log]) The mass-weighted integral of the MRP, in reverse (ie.