imap_processing.idex.idex_l2a.analyze_peaks

imap_processing.idex.idex_l2a.analyze_peaks(tof_high: DataArray, high_sampling_time: DataArray, mass_scale: DataArray, event_num: int, peaks_2d: ndarray) → tuple[NDArray, NDArray, NDArray, NDArray]

Fit an EMG curve to the Time of Flight data around each peak.

Parameters:

• tof_high (xarray.DataArray) – The time of flight array.

• high_sampling_time (xarray.DataArray) – The high sampling time array.

• mass_scale (xarray.DataArray) – Time to mass scale.

• event_num (int) – Dust event number (for debugging purposes).

• peaks_2d (numpy.ndarray) – Nested list of peak indices.

Returns:

• params (numpy.ndarray) – Array of the EMG fit parameters (mu, sigma, lambda) at the corresponding mass. Empty mass slots contain zeros.

• area_under_emg (numpy.ndarray) – Array of the area under the EMG curve at that mass. Empty mass slots contain zeros.