docs: def get_observed_mmax_sigma(self, default=None):
          """
          :returns: the sigma for the maximum observed magnitude
          """
          if not isinstance(self.data['sigmaMagnitude'], np.ndarray):
              obsmaxsig = default
          else:
              obsmaxsig = self.data['sigmaMagnitude'][
                  np.argmax(self.data['magnitude'])]
          return obsmaxsig