Peak Finding and Fitting¶

Peak search¶

PeakSearch locates peaks in a spectrum given the detector resolution (FWHM) and a minimum signal-to-noise ratio.

from wara import Spectrum, PeakSearch

sp = Spectrum.from_file("data.csv")
ps = PeakSearch(sp, fwhm=10, snr=3)
ps.find_peaks()
ps.plot()

Peak fitting¶

Once peaks are found, PeakFit (and AdvancedFit for overlapping peaks) decompose the spectrum into Gaussian components on top of a background model.

from wara import PeakFit

pf = PeakFit(sp, ps.peaks)
pf.fit()
pf.plot()
print(pf.results)

For full API details see the API reference for wara.peaksearch and wara.peakfit.