code: def _pdf(self, x, dist, cache):
          
          output = evaluation.evaluate_density(dist, numpy.arccosh(x), cache=cache)
          output /= numpy.where(x != 1, numpy.sqrt(x*x-1), numpy.inf)
          return output
docs: Probability density function.