docs: def _get_inter_event_tau(self, C, mag, num_sites):
          """
          Returns the inter-event standard deviation (tau), which is dependent
          on magnitude
          """
          base_vals = np.zeros(num_sites)
          if mag <= 4.5:
              return base_vals + C["t1"]
          elif mag >= 5.5:
              return base_vals + C["t2"]
          else:
              return base_vals + C["t1"] + (C["t2"] - C["t1"]) * (mag - 4.5)