docs: def get_magnitude_time_distribution(self, magnitude_bins, time_bins,
                                          normalisation=False, bootstrap=None):
          '''
          Returns a 2-D histogram indicating the number of earthquakes in a
          set of time-magnitude bins. Time is in decimal years!
  
          :param numpy.ndarray magnitude_bins:
               Bin edges for the magnitudes
  
          :param numpy.ndarray time_bins:
              Bin edges for the times
  
          :param bool normalisation:
              Choose to normalise the results such that the total contributions
              sum to 1.0 (True) or not (False)
  
          :param int bootstrap:
              Number of bootstrap samples
  
          :returns:
              2D histogram of events in magnitude-year bins
          '''
          return bootstrap_histogram_2D(
              self.get_decimal_time(),
              self.data['magnitude'],
              time_bins,
              magnitude_bins,
              xsigma=np.zeros(self.get_number_events()),
              ysigma=self.data['sigmaMagnitude'],
              normalisation=normalisation,
              number_bootstraps=bootstrap)