1 / 39

MOS Data Reduction

MOS Data Reduction. Michael Balogh University of Durham. Outline. (Automatic) identification of slits and galaxies Distortion correction Background subtraction Wavelength calibration Flat fields and flux calibration. Data Reduction software.

caitir
Download Presentation

MOS Data Reduction

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. MOS Data Reduction Michael Balogh University of Durham

  2. Outline • (Automatic) identification of slits and galaxies • Distortion correction • Background subtraction • Wavelength calibration • Flat fields and flux calibration

  3. Data Reduction software • IRAF: Can deal with multiobject spectroscopy, but handles the following inelegantly: • wavelength calibration • distortion corrections • Dan Kelson’s recently public software: http://www.ociw.edu/~kelson/ • designed for use with MOS data • handles wavelength calibration and distortion corrections easily • Employs new technique for optimal background subtraction • But is somewhat obscure Note: neither software package deals easily with ultraplex data

  4. MOS data: the spectra

  5. MOS data: flats

  6. MOS data: arc lamps

  7. Ultraplex data

  8. Identification of Objects

  9. Identification of Objects: IRAF

  10. Identification of Objects: IRAF Interactively identify object(s) in each slit Specify extent to extract in 1D spectrum Can be tricky for faint spectra because optimal columns to extract will vary from slit to slit (in some cases will hit bright sky lines, in other cases miss bright part of spectrum)

  11. Identification of Objects: Kelson Identify slits in flat field image Laplacian filter helps define slit edges Pick object location on 2D image (using ds9, for example) Kelson 2003

  12. Distortion correction NIRSPEC: Kelson 2003 Spectra are usually curved, due to instrument distortions

  13. Distortion correction Two options: 1. Rectify image before extracting spectra. Makes reduction easier, but introduces residuals in sky subtraction. 2. Measure distortion, but extract spectra from original frame and map to rectified coordinate frame.

  14. Distortion correction: IRAF Curvature in spatial direction is tricky to correct; not easily implemented. Curvature in spectral direction can be traced when extracting spectrum. Must be done interactively and probably not used when extracting arc spectrum Need to be able to see the spectrum… d d

  15. Distortion Corrections: Kelson • Trace slits in flat field to map distortion in spectral direction • For each slit, trace sky lines (or arc lines) to map distortion in spatial direction

  16. Kelson 2003

  17. Background Correction

  18. Background Correction Usual procedure: Define background region on either side of object Fit polynomial across dispersion Assumes no distortion in spatial direction, so must correct first

  19. Background Correction Rebinning: introduces correlated noise, smears bad pixels, produces artifacts/residuals, and forces sky spectrum to have common pixelization Instead: perform least-squares fit to sky spectrum in original coordinates. This provides better sampling in rectified coordinates. Kelson, PASP, in press

  20. Background Correction Rebinning: introduces correlated noise, smears bad pixels, produces artifacts/residuals, and forces sky spectrum to have common pixelization Instead: perform least-squares fit to sky spectrum in original coordinates. This provides better sampling in rectified coordinates. Kelson, PASP, in press

  21. Background Subtraction 2D LRIS spectrum Spectrum profile in rectified coordinates Compare smoothed version of above with profile from single pixel width Kelson, PASP in press

  22. Background Correction • Define sky regions (either directly, or using s-clipping techniques) • Fit bivariate B-spline (Dierckx 1993) as a function of rectified coordinates • Essentially approximates an interpolating spline along the wavelength coordinate, but with much finer sampling than available in a single CCD row • Can generalize further and fit simultaneously to all spectra in a frame. Thus get improved resolution even if distortions are small. Kelson, PASP, in press

  23. LRIS Raw Sky model Background subtracted rms-smoothed, divided by noise: no residuals! Kelson, PASP, in press

  24. NIRSPEC Raw Sky model Background subtracted rms-smoothed, divided by noise: no residuals! Kelson, PASP, in press

  25. Wavelength calibration

  26. Wavelength calibration I Extract arc lamp spectrum for each slit IRAF: identify a few lines and fit low-order function. Then easy to use this fit to find more lines and improve quality of the fit. Task reidentify to find arc lines in other slits on same image does not work well. Usually have to do each slit separately. Not clear to me if this uses trace information from spectrum.

  27. Wavelength calibration II Kelson (2003) software Automatically identify lines in all slits, and computes pixel-wavelength transformation Don’t know how it works, but it does! Can do in minutes what used to take me days with IRAF.

  28. Flat fielding

  29. Flat fielding • Remove the “slit function”: variation in sensitivity along the slit Needed to correct for uneven slits

  30. Flat fielding • Remove the “blaze”: variation in sensitivity in dispersion direction Needed for flux calibration, unless star observed in every slit

  31. Flat fielding • Remove pixel-to-pixel sensitivity variations. Usually introduces a lot of noise

  32. Flux Calibration • Observe photometric standard through one (or more) slits • Reduce normally, and flat field (remove “blaze” function) • Divide by known spectral shape to get detector response as function of wavelength.

  33. Conclusions For LDSS2 spectra, I find both give similar quality results

More Related