Function to calculate from a set of three (or optionally five or two) output files from the spectrometer, as produced by SpectraSuite, for each of two calibration lamps (D2 and FEL), and calibration data for the lamps, a set of calibration multipliers for the instrument. This function applies a non-linearity correction, an slit-function correction, and stray-light correction, and optionally smoothing by means of process_maya_arrays(). Optionaly a HDR ("high dynamic range") merging of spectra (based on integration time bracketing) using two different integration times is done.

## Arguments

- FEL.raw.counts
source_mscpt The raw counts and integration time data for the FEL light source.

- D2.raw.counts
source_mscpt The raw counts and integration time data for the D2 light source.

- pix.wavelengths
numeric vector The "true" wavelengths in nanometres at each pixel in the detector array.

- wl.range
numeric vector of length two Range of wavelengths for which to compute the calibration.

- FEL.k
numeric vector a numeric vector with n constants for the polynomial for the FEL lamp.

- D2.k
numeric vector a numeric vector with n constants for the polynomial for the D2 lamp.

- method
list Method variant to be used.

- verbose
Logical indicating the level of warnings wanted. Defaults to

`FALSE`

.