You can read a hungarian version here: media:EchelleSpectra.pdf
In this type of convention the starting letters of the imagename refers to the corrections. So forexample a bdfobj.fits is a scientific image corrected with bias dark and flat. We are renaming the images from .fit to .fits, because the IRAF command doecslit only accepts .fits (maybe this depends on the version of IRAF).
We have to do the correction for dark, flat and scientific images too..
We have to do the correction for flat and scientific images.
After this with the implot task open flat.fits and select the useful part of the aperture.
Setting the RCC telescope at Piszkéstető Mountain
If our images do not contain headers like: RA,DEC,EPOCH,IMTYP we have to add them manually.
We can look for coordinates in the SIMBAD database. You have to set the RA coordinates like this: "'44:58:06'", and the DEC coordinates like this: '15:45:34' (no idea why, maybe it depends on the version of IRAF too, but this is the only way the setjd command will work). IMTYP is object for scientific images, and comp for ThAr images.
For the ThAr images you have to leave the hjd field empty!!!
We have to meassure the width of our apertures. One way of doing this is to open an image with DS9 and from the menu use Region => Shape => Projection option and mark a line perpendicular to the spectral lines. In the popup Projection window we can read the width of the apertures.
The doecslit task has basic parameters, which you can set by typing :e in the field of sparams. First run the task, then type a capital i. This will cancel the task and put you back to the xgterm window from the irafterm window. This is necessary, because, you cant change the size of the irafterm window while it is active. Enlarge the window so you can work with it comfortably, then run doecslit again.
WARNING!!! DO NOT CLOSE THE IRAFTERM WINDOW, OR IRAF WILL FREEZE!!! (at least my version of IRAF can't handle this)
After this we can look at the apertures. We can delete undesired orders with hovering over them with the cursor and pressing d, and mark orders with m. You van define the order number with o.
With pressing q we can move forward to the next step of the task. It will ask if we would like to fit the aperture interactively. If you accept, you can look at all the aperture fit moving forward with q and RETURN. You can delete points from the fits hovering over it your cursor and pushing d.
If we have multiple object images, and the fitting our acceptable, you can select no for the interactive question.
After tracking the apertures, the next step in the doecslit task is to identify the lines in the spectral lamp image. You can use the spectral maps of Balázs Csák (http://tesla.gothard.hu/~obs/.po/PO_ThAr.pdf.pdf) and Levente Kriskovics (media:tharmap.pdf), or using an online site like: http://iraf.noao.edu/specatlas/thar/thar.html
By pressing m on a line you can set its wavelength, with d you can delete if you mistyped something. With the + and - signs (on the numpad) you can move between the labeled lines. If you are ready, you can move to the next order by pushing k or to the previous order with j. If we are ready with all the orders, push y and IRAF will try to identify more lines automatically (you have to do this on every order).
After this push f to fit the dispersion function, where you can delete fitted points with d. Here are the identified lines by doecslit:
After this you will have the .ec files, you can view these with the splot task
You can move between the orders with the () buttons, you can zoom in an area with a, zoom out with w and a. The splot task: at the left side of the spectra you can see the H alpha line
Making list files for the wavelength calibrated images:
The fitted curve of the continuum task:
If you would like to exclude a bigger absoprtion region you can mark fitted regions by pressing s at the sides, and pushing f for a refit. This will be valid for the next order too, so you have to push t and f. After this we will have the continuum normalized orders.
A continuum normalized order with the H alpha line on its left side:
Making list files:
Making the 1 dimensional spectra:
You can view the spectrum with the splot task
You can see that the different orders wont overlap over 7100 angstroms an that the edges are not perfect.
We can trim the spectra over 1 with the imreplace task.
After this you can examine your spectra with the splot task. For example you can zoom in on the H alpha line by pressing a twice at the edges of the desired area, and press k two times at the edges of a line to fit a Gauss-curve, k and v for a Voigt-profile, k and l for a Lorentz-profile, you can read the equivalent width in the irafterm window's lower left corner.
The fitted Voigt-profile with the splot task