Finding the peak positions

[Wile_E_Coyote]

Try to play with this:

Trying to find peak areas using example FitAwmi.C

@axel The standard version of the ${ROOTSYS}/tutorials/spectrum/peaks.C script fits “peaks’ heights”, while the version attached below can fit “peaks’ areas” instead (and you can run this version with ROOT 5 or ROOT 6). peaks.C (4.3 KB)

and that:

Trying to find peak areas using example FitAwmi.C

The attached “SpectrumFit.cxx” macro fits 1D spectra using either the AWMI method (an algorithm without matrix inversion) or the Stiefel-Hestens method (a conjugate gradient algorithm) from the TSpectrumFit class. The TSpectrum class is used to find peaks. To try this macro, in a ROOT (5 or 6) prompt, do … … for the “AWMI” method … .x SpectrumFit.cxx ... or ... .x SpectrumFit.cxx++ SpectrumFit(); another random set of peaks (integer areas) SpectrumFit(0, kFALSE); another random set of peaks (…