docs: def s_to_ev(offset_us, source_to_detector_m, array):
      """convert time (s) to energy (eV)
      Parameters:
      ===========
      numpy array of time in s
      offset_us: float. Delay of detector in us
      source_to_detector_m: float. Distance source to detector in m
  
      Returns:
      ========
      numpy array of energy in eV
      """
      lambda_a = 3956. * (array + offset_us * 1e-6) / source_to_detector_m
      return (81.787 / pow(lambda_a, 2)) / 1000.