FAST Spectrograph

Processing FAST Spectrograph Data
Jessica Mink, 2023-March-7

Telescope Data Center
FAST Spectrum
[TDC Home]  [TDC Search] [OIR Home]

FAST Data Directories

All FAST data since January 5, 2005 has been reduced more or less this way

Task

Commands to Reduce Last Night's Data

3. Align spectra for one configuration

In each reduction/working configuration directory for the date, extract and align 2-D spectra using fastpack.roadrunner in IRAF. When IRAF starts, it automatically moves to the FAST IRAF home directory, /home/fast/IRAF2.17. The FAST working directory, ared is linked from that directory, too, so you can use the same paths inside and outside of IRAF.

Sometimes it's worthwhile to run IRAF imstat on COMP, FLAT, or BIAS files to make sure that no files are way off statistically.

Operations are logged in the file road.log. To follow the log, run the shell script wroadlog in a different window but the same directory to follow the log.

1. Catalog FAST raw data archive

2. Transfer raw data to reduction/working directory

3. Align spectra for one configuration

4. Extract, Cross-correlate and Check spectra

5. Make sure that files are archived appropriately

6. Distribute FAST data from one night

 
> ecl
  IRAF2.17 with SAO packages 
  Maintained by Jessica Mink, jmink@cfa.harvard.edu
  January 17, 2023

  This product includes results achieved by the IRAF64 project in 2006-
  2009 directed by Chisato Yamauchi (C-SODA/ISAS/JAXA).

  Welcome to IRAF.  To list the available commands, type ? or ??.  To get
  detailed information about a command, type `help '.  To run  a
  command  or  load  a  package,  type  its name.   Type  `bye' to exit a
  package, or `logout' to get out  of the CL.

  Visit http://github.com/iraf-community/iraf/issues to report problems.
  The following commands or packages are currently defined:

      dataio.     fmisc.      images.     mike.       obslog.     proto.      softools.   svdfit.     tres.       
      fastpack.   hectospec.  language.   mscred.     obsolete.   rvsao.      starbase.   system.     utilities.  
      fitsutil.   hmisc.      lists.      noao.       plot.       saotdc.     stsdas.     tables.     wcstools.   

ecl> cd ared/2023/2023.0122/B2D600T752A2.0
ecl> pwd
/data/mc4/fast/reductions/2023/2023.0122/B2D600T752A2.0

ecl> roadrunner >& road.log
 
In another window, capture the progress in real time by running
 
> wroadlog
42  data images
Gain Report: BinBy = 2 Disp = 600 Tilt = 752. Aper = 2.


Blist = 20 
Flist = 20 
Picked for gain computation

bias1           bias2           flat1           flat2      
-----           -----           -----           -----      
0030.BIAS.fits  0036.BIAS.fits  0148.FLAT.fits  0154.FLAT.fits

       48880.15 + 47416.1 - 767.4073 - 767.6665
gain = ---------------------------------------
       351.8389 * 351.8389 - 8.989439 * 8.989439

noise= 8.989439 * 0.76599570139646 / 1.4142135623731

ChipGain

Gain            Noise
------          -----
0.766           4.869
Deleting /tmp/b123233uf.fits and /tmp/f123233tf.fits

...

Shift Computations for reidentify

COMP            SHIFT   FWHM
----            -----   ----
0141.COMP.fits  -8.188 4.15
Matchlist
COMP    HDRFILE                 DTmin   Sep             Adj
----    ---------------------   ------  -----------     ---
0141    0142.TIC87496804.fits    4.234  0:00:01.505       1
transform input="0142.TIC87496804" output="0142.TIC87496804TF" fitnames="0141.COMP"

Community IRAF V2.17 fast@cfa0 Mon 17:37:49 23-Jan-2023
  Transform 0142.TIC87496804 to 0142.TIC87496804TF.
  Conserve flux per pixel.
  User coordinate transformations:
    0141.COMP
  Interpolation is spline3.
  Using edge extension for out of bounds pixel values.
  Output coordinate parameters are:
    x1 =      5577., x2 =      7584., dx =     0.7489, nx = 2681, xlog = no
    y1 =         1., y2 =       161., dy =         1., ny =  161, ylog = no
Directory       /data/mc4/fast/reductions/2023/2023.0122/B2D600T752A2.0
Reidentify      Community IRAF V2.17 fast@cfa0 Mon 17:37:45 23-Jan-2023
Reference       HENEAR_TEMPLATE.600.2.20.752
         
Comp            MinL    MaxL    AvShift Smin    Smax    AvRMS   RMSmin  RMSmax
---------       ----    ----    ------- -----   -----   ------  ------  ------
0141.COMP         18      19      -6.22 -6.37   -6.06   0.0482   0.031   0.056