code: def marginalization_bins(self):
          
          log_mean = np.log10(self.mean())
          # Default is to marginalize over two decades,
          # centered on mean, using 1000 bins
          return np.logspace(-1. + log_mean, 1. + log_mean, 1001)/self._j_ref
docs: Binning to use to do the marginalization integrals