logp

Import in python using from mcutils.scripts.non_central_chi.correct import logp

mcutils.scripts.non_central_chi.correct.logp(observed_signal, log_noise, ncoils, f_no_signal, dist_param1, dist_param2, use_dist='gamma')

Computes the logp of the total histogram

Parameters

• observed_signal – observed signal level

• log_noise – log of the noise variance

• ncoils – number of coils

• f_no_signal – fraction of observed signal with no true underlying signal

• dist_param1 – shape of the gamma distribution or log(mean) of lognorm distribution

• dist_param2 – rate of the gamma distribution or std of lognorm distribution

• use_dist – which distribution to use to model the true signal

Returns

log(P) for the observed signal