lm = self._logquasipdf(m, p, b)