Feels like I'm cheating - a bit.
So I have;
BKGFIT - which fits the background alone
FSPFIT - which fits he complete spectrum
RAWDD - displays the raw`data
XCALIBER - x-axis calibration to energy
FSPDIS - display the full spectrum
STGEN - generate a standard peak
The order of events for the way the code works is as follows;
- Read and display the raw data (RAWDD)
- Generate a 'standard' peak from the peak data (STGEN)
- Fit the background (BKGFIT)
- Search for regions where there may be peaks (FSPFIT)
- Add a peak to the region (FSPFIT)
- Repeat 3,4,5 until there are no further peak regions identified
- Convert the x-axis channel numbers into energies (XCALIBER)
- Display the fitting results (FSPDIS)
During the fitting the user has full control over the level to which the code will identify and attempt to fit smaller and smaller peaks.
Can't believe we did this on machines available at the time!