WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe...

26
WiFeS workshop, RSAA 28 April 2010 WiFeS IRAF DATA REDUCTION PACKAGE

Transcript of WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe...

Page 1: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

WiFeS workshop, RSAA 28 April 2010

WiFeS IRAF DATA

REDUCTION PACKAGE

Page 2: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Turn raw files into calibrated 3D images

λ

yx

Page 3: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

( http://rsaa.anu.edu.au/observing/ssowiki/index.php/WiFeS_data_reduction )

Link to scripts here:RSAA home page ! WiFeS in menu ! ‘Data Reduction’ at bottom of wiki page

WiFeS package

Gemini package (NIFS)

IRAF/NOAO spectroscopy tasks

Requires IRAF (version 2.14)

+ external packages: Gemini, TABLES & STSDAS

An external IRAF package based on Gemini for NIFS

Page 4: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Reformat files

Process calibration data

The package is set up for this general process:

wftable

Create image lists

wfcal

wftable

wfreduceApply calibrations &

construct data cubes

file lists

calibration solutions

Set I/O directorieswifes

Smooth-spectrum standard

Flux standard

Science data

Page 5: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Raw files - simple FITS format

Multi-extension FITS files

Convert to multi-extension FITS

T2m3wr_20090405_124930[0][0]

T2m3wr_20090405_124930[1][4138,4096]image data + header

primary header unit⤺⤺

T2m3wr.20090405.134344.123[0][4138,4096]

image data + header⤺

wftable: Convert images to MEF & create file lists

Create image lists for grating+dichroic

wftable

wftable

Images assigned to text list files by the IMAGETYP header keyword value

IMAGETYP FileARC list_arc

FLAT list_gflat

SKYFLAT list_sflat

WIRE list_wire

ZERO list_zero

OBJECT list_obj

wfcal

wfreduce

}list files

Page 6: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Prepare images

wavelength solution

spatialsolution

master flat field frame

median bias

File listsCalibration files

Combine bias frames

Normalise flat field Identify & fit arc lines Find slice centres

Fit & correct slit profile

Measure image shift

wfreduce

dorefflat

dowfprepare

dozero

doarc dowiredogflat

dosflat

wfcal: Steps for processing calibration files

reference frame

Page 7: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Expand binned data

Overscan bias subtraction

Bias frame surface fitting

`Prepare’

wfreduce: Apply wfcal results & construct cubes

Bias subtractionCombine frames

Sky subtraction

Rescale & form data cube

DAR correction

`Reduce’

dowfprepare

docombine

doreduce

dofitcoorddotransform

dodar

dosubskydoslitsky

overscan or fitbias

doskyframe

dozero

Combine shuffled slices

dosubshuffled

data cube

Subtract median bias frame (from wfcal)

Page 8: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

`Prepare’

wfreduce: Apply wfcal results & construct cubes

Bias subtractionCombine frames

Sky subtraction

Rescale & form data cube

DAR correction

`Reduce’

dowfprepare

docombine

doreduce

dofitcoorddotransform

dodar

dosubskydoslitsky

overscan or fitbias

doskyframe

dozero

Median-combine multiple frames

Combine shuffled slices

dosubshuffled

data cube

Page 9: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

`Prepare’

wfreduce: Apply wfcal results & construct cubes

Bias subtractionCombine frames

Sky subtraction

Rescale & form data cube

DAR correction

`Reduce’

dowfprepare

docombine

doreduce

dofitcoorddotransform

dodar

dosubskydoslitsky

overscan or fitbias

doskyframe

dozero

Combine shuffled slices

dosubshuffled

data cube

Subtract a separate blank sky frame

Subtract a sub-aperture average sky spectrum

Subtract `long-slit’ spectrum in each slice

Page 10: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

`Prepare’

wfreduce: Apply wfcal results & construct cubes

Bias subtractionCombine frames

Sky subtraction

Rescale & form data cube

DAR correction

`Reduce’

dowfprepare

docombine

doreduce

dofitcoorddotransform

dodar

dosubskydoslitsky

overscan or fitbias

doskyframe

dozero

Combine shuffled slices

dosubshuffled Shift and combine nod- or subnod-and-shuffled data

data cube

Page 11: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

`Prepare’

wfreduce: Apply wfcal results & construct cubes

Bias subtractionCombine frames

Sky subtraction

Rescale & form data cube

DAR correction

`Reduce’

dowfprepare

docombine

doreduce

dofitcoorddotransform

dodar

dosubskydoslitsky

overscan or fitbias

doskyframe

dozero

Combine shuffled slices

dosubshuffled

Extract 25 slices

Apply flat field

data cube

Page 12: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

`Prepare’

wfreduce: Apply wfcal results & construct cubes

Bias subtractionCombine frames

Sky subtraction

Rescale & form data cube

DAR correction

`Reduce’

dowfprepare

docombine

doreduce

dofitcoorddotransform

dodar

dosubskydoslitsky

overscan or fitbias

doskyframe

dozero

Combine shuffled slices

dosubshuffled

Transform image scales

Stack slices into 3D image

data cube

Page 13: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

`Prepare’

wfreduce: Apply wfcal results & construct cubes

Bias subtractionCombine frames

Sky subtraction

Rescale & form data cube

DAR correction

`Reduce’

dowfprepare

docombine

doreduce

dofitcoorddotransform

dodar

dosubskydoslitsky

overscan or fitbias

doskyframe

dozero

Combine shuffled slices

dosubshuffled

data cube

Page 14: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

Extract telluric star spectrum

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

λx

y

(dotelcorrect)

dotelextract

Page 15: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

Extract telluric star spectrum

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)

dotelextract

Page 16: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)

Na D

Ca II

Remove stellar absorption

dofixtel

Page 17: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)

Remove stellar absorption

dofixtel

Page 18: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)

Remove stellar absorption

dofixtel

Page 19: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)

Remove stellar absorption

Na D

Ca II

dofixtel

Page 20: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)

Normalise telluric spectrum

Extract and normalise example object spectrum

dotelcorrect

Page 21: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

Correct telluric features

wfreduce: Further steps for standard star data

(dotelcorrect)dotelcorrect

Page 22: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)Extract flux spectrum

dofluxextract

Page 23: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)

Page 24: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)Calculate flux

scale

dofluxcal

Page 25: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

Process smooth-spectrum standard

star (red)

Process flux standard star

Process science frames & apply standard star

calibrations

telluric spectrum

flux calibration

dotelextract

dofixtel

dofluxextract

(dotelcorrect)

dofluxcal

wfreduce: Further steps for standard star data

(dotelcorrect)Calculate flux

scale

dofluxcal

Page 26: WiFeS IRAF DATA R P - Intranet - ANU...dodar dosubsky doslitsky overscan or fitbias doskyframe dozero Combine shuffled slices dosubshuffled data cube Subtract a separate blank sky

wfreduce: a few last tasks...

doscience

dobincube

dosimcube

output datacube T2m3wr_20090405_124930_R3000_RT560_res_t.fits

..._res_bin.fits

..._res_tds9.fits

..._sci.fits

..._res_c.fits

Form projected image andintegrated spectrum

Re-bin cube spatial axes

Save simple FITS cube file