For almost all purposes we recommend the use of the
Spectral line indices (Lick indices)
Separate continuum and flux windows are used to determine spectral
line indices analogous to the Lick indices. These are stored in Line index HDU (HDU 5) of
Details of continuum fitting and line measurements
All of the line parameters are measured in the observed frame, and no correction has been made for the instrumental resolution. Line indices are invariant under transforming from observed to rest frame.
The continuum is fit using a median/mean filter. A sliding window is created of length 300 pixels for galaxies and stars or 1000 pixels for quasars. Pixels closer than 8 pixels(560km/s) for galaxies and stars or 30 pixels (2100 km/s) for QSOs to any reference line are masked and not used in the continuum measurement. The remaining pixels in the filter are ordered and the values between the 40th and 60th percentile are averaged to give the continuum. The category=1 lines are fit with a cruder continuum which is given by a fifth order polynomial fit which iteratively rejects outlying points.
Reference Line List
The list of lines which are fit are given as ASCII line list and as HTML line table below. Note that many times that a single line in the table actually represents multiple features. Since the line fits are allowed to drift in wavelength somewhat, the exact precision of the lines are not important. The wavelength precision does become important for the emission line determination.
Every line in the reference list is fit as a single Gaussian on top of the continuum subtracted spectrum. Lines that are deemed close enough are fitted simultaneously as a blend. The basic line fitting is performed by the SLATEC common mathematical library routine SNLS1E which is based on the Levenberg-Marquardt method. Parameters are constrained to fall within certain values by multiplying the returned chi-squared values by a steep function. Any lines with parameters falling close to these constraints should be treated with caution. The constraints are: sigma > 0.5 Angstrom, sigma < 100 Angstrom, and the center wavelength is allowed to drift by no more than 450 km/sec for stars and galaxies or 1500 km/sec for QSOs, except for the CIV line which is allowed to be shifted by as much as 3000 km/sec.
Testing the results
There are a number of ways that the line fitting can fail. If the continuum is bad the line fits will be compromised. The median/mean filtering routine will always fail for white dwarfs, some A stars as well as late-type stars. In addition is has trouble for galaxies with a strong 4000 Angstrom break. Likewise the line fitting will have trouble when the lines are not really Gaussian. The Levenberg-Marquardt routine can fall into local minima, which can happen when there is self-absorption in a QSO line or both a narrow and broad component for example. One should always check the chi-squared values to evaluate the quality of the fit. Examples of the quality of the line fits for over 100 spectra can be found in ps file with sample line fits (.ps.gz, 16 MB). In these plots, the red line is the continuum and the blue line is the continuum plus the sum of the Gaussian fits.
You can download the line list as ASCII table.
Last modified: Mon Sep 27 13:29:36 CDT 2004
|Send your questions to firstname.lastname@example.org|