IRAF

Transcription

1 IRAF

2 IRAF exercise

3 I tutorial iraf/misc exer211.tar.z login.cl % tar -zxvf exer211.tar.z % cl> # 3

4 II Exercise 1. intro IRAF xgterm ds9 xgterm IRAF intro % ds9 # ds9 % cl # IRAF IRAF login cl> cd exercises/intro # exercises IRAF login IRAF * 1 cl> stty # cl> stty xgtermc nlines=40 ncols=80 # xgterm 80, 40 [ ] RFITS cl> unlearn rfits cl> rfits fintro* "" junk old+ IRAF cl> diskspace cl> path cl> dir cl> dir *.fits cl> dir uparm cl> dir l+ cl> package cl> help images cl> phelp imheader cl> lpar imh # # # # # uparm # # # "images" # "imheader" # #? * 1 xgterm,xterm,gterm,vt640,vt100 4

5 cl> lpar imhead # cl> unlearn imhead # cl> imhead im010 # - cl> imhead # cl> imhead im010 l+ cl> ^ page #? cl> lpar imhead cl> imhead im010 l+ u- cl> epar imhead # : # eparam [crtl-c ] cl> lpar imhead # cl> imhead im011 page cl> dir uparm # uparm "imheader" cl> unlearn imhead # cl> dir uparm # cl> imhead im010 l+ page # "exptime" cl> unlearn hselect cl> lpar hselect cl> hselect im*.fits $I,exptime yes # hselect - cl> unlearn hedit cl> lpar hedit cl> hedit im*.fits notice "test data" add+ # (NOTICE) (test data) cl> imhead im010 l+ # cl> hselect im*.fits $I,notice yes (exercises/intro) (im010,im011) cl> dir im*.fits cl> imhead im*.fits # 2 # 5

6 cl> unlearn display cl> epar display cl> lpar display cl> display im010 1 # fill=yes # display # ds9 frame zoom ( View information panel ) panner box im cl> display im011 2 frame frame blink 2 Blink interval single ds9 cl> display im010 1 # im010 cl> unlearn imexamine cl> lpar imexamine # cl> imexamine im010 1 # ds9 - imexamine "interactive image cursor mode" a. "a" xgterm x y 6

7 FWHM - b. "?" - xgterm "" c. 3 im010 3 "a" - d. ds9 "d" xgterm "im011" e. ds9 "" "imexamine" 2 1 x=-0.53 y=-1.68 cl> lpar help cl> help imshift sec=description # description cl> help imshift sec=example cl> unlearn imshift cl> lpar imshift cl> imshift im011 s # im011 im010 s011 im010 s011 2 cl> unlearn imsum imarith a. cl> lpar imsum cl> imsum im010,s011 aver1 pixt=r calct=r option=average v+ <"epar imsum" ":go" > [ REAL( ) ] b. cl> lpar imarith 7

8 cl> imarith im010 + s011 aver2 pixt=r calct=r v+ <"epar imarith" ":go" > cl> imarith aver2 / 2 aver2 # c. [ ] cl> unlearn imstatistics cl> lpar imstat cl> imstat aver*.fits IRAF cl> history cl> ^ cl> e lpar cl> e cl> history 100 cl> history 100 > hfile cl> page hfile cl> history 100 page # # # lpar DEL # / # 100 # # # > IRAF cl> help plot cl> contour s011 # "plot" # ds9 8

9 cl> =gcur :.help # "interactive graphics cursor mode" # IRAF Z - "Z" A - C caps cl> surface s011 cl> surface s011 >G meta # # cl> implot s011 # "?" : command c - column plot, :l :c :.write meta # < "implot" > "Z" # cl> implot dev$wpix # IRAF (M51) :w world # RA DEC :f %H # hh:mm:ss :.write meta c # "c" :f %h # dd:mm:ss :.write meta 9

10 # meta cl> unlearn gkimosaic gkidir cl> gkidir meta # cl> lpar gkimosaic cl> gkimos meta # # cl> gkimos dev$vdm.gki # IRAF <space bar> # # IMEXAMINE cl> display dev$wpix 1 # IRAF M51 wpix pix wpix world coordinate cl> imexamine [ 224, 131 ]? # "" z # m # [ ] s # :epar # ds9 xgterm : l # g # "interactive graphics cursor mode" :naverage 10 # 10 i # "interactive image cursor mode" 10

11 ds9 [ ] # ds9 cl> display dev$wpix 1 xmag=0.8 ymag=0.8 cl> wcslab dev$wpix 1 # "gflush" xmag, ymag 1 cl> dir cl> delete hfile,meta cl> imdelete *.fits ver+ # # im*.fits cl> logout IRAF A Beginner s Guide to Using IRAF, by Jeannette Barnes, August end of exercise

12 III Exercise 2. ccd1 IRAF CCD Kitt Peak George Jacoby 4 2 IRAF 2 2 xgterm( ) IRAF IRAF exercises phot cl> cd cl> cd exercises/phot rfits cl> unlearn rfits cl> rfits fm92* junk old+ m92*.fits V B 2 4 cl> dir cl> imhead m92*.fits 25 S/N 5 short 16 1 IMAGES imcombine m92*.fits 2 CCD 32 ( ) 12

13 implot cl> display m cl> implot m92006 ( ) :l m n, :c m n, C (m,n )? ?? ? colbias cl> noao cl> imred cl> bias cl> help colbias cl> unlearn colbias cl> epar colbias # unlearn # lpar colbias cl> lpar colbias input = "m92*.fits" Input images output = "%m%tr%92*.fits" Output images (bias = "[335:350,2:510]") Bias section (trim = "[1:318,2:510]") Trim section (median = no) Use median instead of average in column bias? (interactive = yes) Interactive? (function = "chebyshev") Fitting function (order = 1) Order of fitting function (low_reject = 3.) Low sigma rejection factor (high_reject = 3.) High sigma rejection factor (niterate = 1) Number of rejection iterations 13

14 (logfiles = "") (graphics = "stdgraph") (cursor = "") (mode = "l") Log files Graphics output device Graphics cursor input? section cl> sections %m%tr%92*.fits colbias cl> colbias cl> dir cl> imhead tr*.fits cl> implot tr92007 # # imarith imarith cl> files tr*.fits > zlist cl> edit zlist cl> cl> unlearn imarith cl> epar imarith # # lpar cl> lpar imarith operand1 = "@zlist" op = "-" operand2 = "tr92001" result = "@zlist" (title = "") (divzero = 0.) (hparams = "") Operand image or numerical constant Operator Operand image or numerical constant Resultant image Title for resultant image Replacement value for division by zero List of header parameters 14

15 (pixtype = "") (calctype = "") (verbose = yes) (noact = no) (mode = "l") Pixel type for resultant image Calculation data type Print operations? Print operations without performing them? imarith cl> imarith noao.imred.genric darksub 2 ( 1 ) imstatistics 2 imarith cl> phelp imstatistics cl> imstat tr92006,tr92007 fields= image,mode cl> imarith tr92006 / 1313 Bflat cl> imarith tr92007 / 1468 Vflat cl> implot Bflat # Vflat cl> display Bflat 1 # Vflat??? imarith 2 cl> imarith tr92010,tr92011 /??? n92010,n92011 cl> imarith tr92014,tr92015 /???????????????? cl> imhead n92*.fits display implot? ccdred mkskyflat cl> imdelete tr*.fits,n92*.fits,bflat,vflat ver+ cl> del zlist cl> imhead m92*.fits 15

16 2 cl> imhead m92*.fits ccdred noao.imred ccdred cl> package cl>???? # ccdred ccdred translation 3 IMAGETYP(object,flat,zero )EXP- TIME( )SUBSET( ) ccdlist cl> imheader m92015 l+ cl> unlearn ccdred cl> lpar ccdlist cl> ccdlist m92*.fits cl> ccdlist m92*.fits l+ # imagettyp,extime,subset KPNO translation cl> setinstrument? direct # translation [ ccdred epar verbose yes : ] [ ccdproc epar 1 : ] cl> ccdlist m92*.fits cl> type subsets 16

17 cl> dir ccddb$kpno # kpno translation files cl> type ccddb$kpno/direct.dat cl> lpar ccdred ccdred short pixeltype backup zerocombine flatcombine ccdproc biassec trimsec image epar cl> imhead m92015 l+ cl> epar ccdproc cl> lpar ccdproc cl> lpar ccdproc images = "m92*.fits" List of CCD images to correct (ccdtype = "object") CCD image type to correct (max_cache = 0) Maximum image caching memory (in Mbytes) (noproc = no) List processing steps only?\n (fixpix = no) Fix bad CCD lines and columns? (overscan = yes) Apply overscan strip correction? (trim = yes) Trim the image? (zerocor = yes) Apply zero level correction? (darkcor = no) Apply dark count correction? (flatcor = yes) Apply flat field correction? (illumcor = no) Apply illumination correction? (fringecor = no) Apply fringe correction? (readcor = no) Convert zero level image to readout correction? (scancor = no) Convert flat field image to scan correction?\n (readaxis = "line") Read out axis (column line) (fixfile = "") File describing the bad lines and columns (biassec = "[335:350,2:510]") Overscan strip image section (trimsec = "[1:318,2:510]") Trim data section (zero = "") Zero level calibration image 17

18 (dark = "") (flat = "") (illum = "") (fringe = "") (minreplace = 1.) (scantype = "shortscan") (nscan = 1) (interactive = yes) (function = "chebyshev") (order = 1) (sample = "*") (naverage = 1) (niterate = 1) (low_reject = 3.) (high_reject = 3.) Dark count calibration image Flat field images Illumination correction images Fringe correction images Minimum flat field value Scan type (shortscan longscan) Number of short scan lines\n Fit overscan interactively? Fitting function Number of polynomial terms or spline pieces Sample points to fit Number of sample points to combine Number of rejection iterations Low sigma rejection factor High sigma rejection factor (grow = 0.) (mode = "l") Rejection growing radius zerocor flatcor cl> ccdproc cl> page logfile cl> imhead m92*.fits cl> imhead m92014 l+ # exercise Type "help ccdred" to see a list of the tasks in this package. Each task has an online help page. Also see the list of "Additional Help Topics". A User s Guide to CCD Reductions with IRAF, by Phil Massey, June There is a "demo" task in the IMRED.CCDRED.CCDTEST package that may be instructive to run as well. 18

19 end of exercise

20 IV Exercise 3. ccd2 ccdred twodspec.longslit Kitt Peak Coude Feed telescope rfits IRAF cl> cd cl> cd exercises/spec cl> unlearn rfits cl> help unlearn cl> rfits fsp* "" junk old+ cl> dir cl> imhead sp*.fits # IRAF # spectral data # rfits # rfits # # ccdred ( phelp ccdred ccdred???? ) cl> package cl>???? cl> unlearn ccdred cl> imhead sp0020 l+ cl> ccdlist sp*.fits # # sp0020 long header # CCD setinstrument setinstrument ccdproc cl> dir ccddb$kpno # CL cl> type ccddb$kpno/coude.cl # ccdproc cl> setinstrument? # Instrument ID? coude # 20

21 ccdred : ccdproc flatcor no : 10 1 zerocombine cl> lpar zerocombine # zerocombine cl> lpar zerocombine input = (output = "Zero") (combine = "average") (reject = "minmax") (ccdtype = "zero") (process = no) (delete = no) (clobber = no) (scale = "none") (statsec = "") (nlow = 0) (nhigh = 1) (nkeep = 1) (mclip = yes) (lsigma = 3.) (hsigma = 3.) (rdnoise = "0.") (gain = "1.") (snoise = "0.") (pclip = -0.5) (blank = 0.) (mode = "l") List of zero level images to combine Output zero level name Type of combine operation Type of rejection CCD image type to combine Process images before combining? Delete input images after combining? Clobber existing output image? Image scaling Image section for computing statistics minmax: Number of low pixels to reject minmax: Number of high pixels to reject Minimum to keep (pos) or maximum to reject (neg Use median in sigma clipping algorithms? Lower sigma clipping factor Upper sigma clipping factor ccdclip: CCD readout noise (electrons) ccdclip: CCD gain (electrons/dn) ccdclip: Sensitivity noise (fraction) pclip: Percentile clipping parameter Value if there are no pixels zerocombine ccdtype zero sp*.fits cl> zerocombine sp*.fits imstatistics 21

22 cl> imstat sp000?.fits,sp0010,zero # cl> display Zero 1 fill+ # cl> implot Zero # overscan region trim parameters cl> implot sp0014 overscan region biassec(overscan region) trimsec(trim parameters) imhead cl> imhead sp0020 l+ ccdproc epar biassec trimsec zerocor yes ccdtype cl> epar ccdproc ccdproc cc> lpar ccdproc images = "" List of CCD images to correct (ccdtype = "") CCD image type to correct (max_cache = 0) Maximum image caching memory (in Mbytes) (noproc = no) List processing steps only?\n (fixpix = no) Fix bad CCD lines and columns? (overscan = yes) Apply overscan strip correction? (trim = yes) Trim the image? (zerocor = yes) Apply zero level correction? (darkcor = no) Apply dark count correction? (flatcor = no) Apply flat field correction? (illumcor = no) Apply illumination correction? (fringecor = no) Apply fringe correction? (readcor = no) Convert zero level image to readout correction? (scancor = no) Convert flat field image to scan correction?\n (readaxis = "line") Read out axis (column line) (fixfile = "") File describing the bad lines and columns (biassec = "[105:130,1:1022]") Overscan strip image section (trimsec = "[34:74,1:1022]") Trim data section 22

23 (zero = "Zero") (dark = "") (flat = "") (illum = "") (fringe = "") (minreplace = 1.) (scantype = "shortscan") (nscan = 1) (interactive = yes) (function = "chebyshev") (order = 1) (sample = "*") (naverage = 1) (niterate = 1) (low_reject = 3.) (high_reject = 3.) (grow = 0.) (mode = "l") Zero level calibration image Dark count calibration image Flat field images Illumination correction images Fringe correction images Minimum flat field value Scan type (shortscan longscan) Number of short scan lines\n Fit overscan interactively? Fitting function Number of polynomial terms or spline pieces Sample points to fit Number of sample points to combine Number of rejection iterations Low sigma rejection factor High sigma rejection factor Rejection growing radius 4 flatcombine cl> implot sp0014 :i sp0015 o :l 512 :i sp0016 o :l 512 :i sp0017 o :l 512 flatcombine cl> epar flatcombine > lpar flatcombine input = (output = "Flat") (combine = "average") List of flat field images to combine Output flat field root name Type of combine operation 23

24 (reject = "crreject") (ccdtype = "flat") (process = yes) (subsets = no) (delete = no) (clobber = no) (scale = "mode") (statsec = "") (nlow = 1) (nhigh = 1) (nkeep = 1) (mclip = yes) (lsigma = 3.) (hsigma = 3.) (rdnoise = "rdnoise") (gain = "gain") (snoise = "0.") (pclip = -0.5) (blank = 1.) (mode = "l") Type of rejection CCD image type to combine Process images before combining? Combine images by subset parameter? Delete input images after combining? Clobber existing output image? Image scaling Image section for computing statistics minmax: Number of low pixels to reject minmax: Number of high pixels to reject Minimum to keep (pos) or maximum to reject (neg Use median in sigma clipping algorithms? Lower sigma clipping factor Upper sigma clipping factor ccdclip: CCD readout noise (electrons) ccdclip: CCD gain (electrons/dn) ccdclip: Sensitivity noise (fraction) pclip: Percentile clipping parameter Value if there are no pixels flatcombine cl> flatcombine sp*.fits cl> imhead Flat cl> implot Flat overscan trim zero correct cl> files sp*.fits > plist cl> edit plist cl> cl> cl> cl> imhead sp0025 l+ # # bias flat # CCD # CCD # notice reduction flags # notice reduction history twodspec.longslit response ( ) (response) response Interactive Curve Fitting (ICFIT) 24

25 cl> twodspec cl> longslit lpar longslit dispaxis 2 cl> unlearn response cl> phelp response cl> lpar response cl> response Flat Flat Resp? # look at cursor options :function cheby # chebyshef :order 5 # f k h cl> implot Resp cl> lpar combine cl> combine sp0011,sp0012,sp0013 Sky reject=crreject \ scale=mode weight=!exptime rdnoise=rdnoise gain=gain cl> implot Sky cl> unlearn imarith cl> imarith Sky / Resp SkyFlat cl> display SkyFlat 1 fill+ cl> implot SkyFlat illumination illumination function ( ) (default 5) illumination function cl> phelp illumination cl> unlearn illumination cl> lpar illum nbin=5 cl> illum SkyFlat Illum (5 ) 25

26 ? # display cursor options :function cheby :order 2 # straight line :high 1 # set sigma rejection to 2 :grow 2 # radius of pixel rejection :niterate 2 # iterate the solution f # fit new parameters # accept and go to next bin f cl> implot Illum cl> imarith SkyFlat / Illum test cl> implot test # Illum # illumination function cl> imarith Resp * Illum SuperFlat cl> implot SuperFlat :l :i Resp o :l # Resp SuperFlat cl> epar ccdproc # flatcor yes flat SuperFlat cl> cl> imhead sp0020 l+ # Exercise Type "help ccdred" to see a list of the tasks in this package. Each task has an online help page. Also see the list of "Additional Help Topics". There is online help for both RESPONSE and ILLUMINATION. A User s Guide to CCD Reductions with IRAF, by Phil Massey, June

27 A User s Guide to Reducing Slit Spectra with IRAF, by Phil Massey, Frank Valdes, Jeannette Barnes, April end of exercise

28 V Exercise 4. photom APPHOT PHOTCAL exercise 2 M92 CCDPROC xgterm IRAF login xgterm (xterm, kterm aperture photometry IRAF exercises/phot cl> cd cl> cd exercises/phot m92*.fits 4 B,V 2 cl> dir cl> imhead m92*.fits cl> imhead m92010 l+ APERTURE PHOTOMETRY EXPTIME FILTERS AIRMASS ASTUTIL SETAIRMASS * 2 cl> astutil cl> phelp setairmass cl> unlearn setairmass cl> lpar setairmass cl> setairmass m92*.fits updatecl> setairmass m92*.fits cl> imhead m92014 l+ # KPNO # * 2 effective airmass help 28

29 cl> bye # astutil aperture photometry DIGIPHOT APPHOT cl> digiphot cl> apphot (FWHM) FWHM IMEXAMINE cl> display m cl> imexamine r # # FWHM Gaussian FWHM FWHM FWHM IMEXAMINE d r d r d r # m92011 # # m92014 # # m92015 # # (uit) QPHOT* 3 EPAR cl> unlearn apphot cl> epar phot ap> lpar phot image = "" cbox = 5. annulus = Input image Centering box width in pixels Inner radius of sky annulus in pixels * 3 Quick aperture PHOTometry 29

30 dannulus = apertures = "10" (coords = "") (output = "default") (plotfile = "") (zmag = 26.) (exposure = "exptime") (airmass = "airmass") (filter = "filters") (obstime = "ut") (epadu = 14.) (interactive = yes) (radplots = no) (verbose = no) (graphics = "stdgraph") (display = "stdimage") (icommands = "") (gcommands = "") (mode = "l") Width of the sky annulus in pixels List of photometry apertures Coordinate list Results file Plot metacode file Zero point of magnitude scale Exposure time image header keyword Airmass image header keyword Filter image header keyword Time of observation image header keyword Instrument gain Interactive mode Plot the radial profiles in interactive mode Print messages Graphics device Display device Image cursor: [x y wcs] key [cmd] Graphics cursor: [x y wcs] key [cmd] QPHOT cbox=5, annulus=15, dannulus=10, apertures=10 annulus sky annulus dannulus ( ) docs/m92.ps PostScript cl> display m cl> phot m92010? # i # 6 [extraction box ] ["v" ] [centering box 2.5 "return" ] 30

31 [inner sky radius 15 "return" 15 ] [outer sky radius 25 "return" 10 ] [ 10 "return" ] # [ - center, sky, magnitude, error ] [ radial plot ] w # :radplot yes # radial profile option [ ] 6, 9, 5, 8, 7, 23, 13, 20, 18 -.mag.1 TXDUMP cl> dir cl> page m92010.mag.1 cl> lpar txdump cl> txdump m92010.mag.1 image,xcenter,ycenter,mag,msky,stdev yes cl> txdump m92010.mag.1 mag,merr yes graph point+ QPHOT V B cl> display m [ 6 31

32 - x 4.5 ] TXDUMP LISTS LINTRAN cl> txdump m92010.mag.1 xcenter,ycenter yes > coords cl> type coords cl> display m cl> tvmark 1 coords mark=circle radii=10 color=205 QPHOT B,V cbox cl> epar phot ap> lpar phot image = "m92011,m92014,m92015" Input image cbox = 7. Centering box width in pixels annulus = 15. Inner radius of sky annulus in pixels dannulus = 10. Width of the sky annulus in pixels apertures = "10" List of photometry apertures (coords = "coords") Coordinate list (output = "default") Results file (plotfile = "") Plot metacode file (zmag = 26.) Zero point of magnitude scale (exposure = "exptime") Exposure time image header keyword (airmass = "airmass") Airmass image header keyword (filter = "filters") Filter image header keyword (obstime = "ut") Time of observation image header keyword (epadu = 14.) Instrument gain (interactive = no) Interactive mode (radplots = no) Plot the radial profiles in interactive mode (verbose = no) Print messages (graphics = "stdgraph") Graphics device (display = "stdimage") Display device (icommands = "") Image cursor: [x y wcs] key [cmd] (gcommands = "") Graphics cursor: [x y wcs] key [cmd] (mode = "l") 32

33 cl> phot cl> dir *.mag* # mag? cl> txdump m92*.mag.1 mag,merr yes graph point+ cl> txdump *.mag* xcenter,ycenter,mag,merr,ifilter yes Merr - V=60 B=50 PHOTOMETRIC CALIBRATIONS PHOTOCAL cl> photcal cl> unlearn photcal MKCATALOG stds fstds.dat stds MKCATALOG cl> page stds cl> page fstds.dat QPHOT V B-V cl> edit imsets cl> type imsets M92 : m92010 m92014 M92 : m92011 m92015 *.mag* imsets V=60 B=50 cl> phelp mknobsfile cl> epar mknobsfile ph> lpar mknobsfile photfiles = "*.mag*" idfilters = "60,50" imsets = "imsets" observations = "sobs" The input list of APPHOT/DAOPHOT databases The list of filter ids The input image set file The output observations file 33

34 (obsparams = "") The input observing parameters file (obscolumns = "2,3,4") The format of obsparams (minmagerr = 0.001) The minimum error magnitude (shifts = "") The input x and y coordinate shifts file (apercors = "") The input aperture corrections file (aperture = 1) The aperture number of the extracted magnitude (tolerance = 8.) The tolerance in pixels for position matching (allfilters = yes) Output only objects matched in all filters (verify = no) Verify interactive user input? (verbose = yes) Print status, warning and error messages? (mode = "l") cl> mknobsfile sobs stds 6, 9, 5, 8, 7, 23, 13, 20, 18 9 cl> edit sobs MKCONFIG comfiguration file stds, sobs m92fig MKCONFIG cl> phelp mkconfig cl> type m92fig V, B-V FITPARM cl> phelp fitparams cl> epar fitparams ph> lpar fitparams observations = "sobs" List of observations files catalogs = "stds" List of standard catalog files config = "m92fig" Configuration file parameters = "calib" Output parameters file (weighting = "uniform") Weighting type (uniform,photometric,euations) (addscatter = yes) Add a scatter term to the weights? (tolerance = E-5) Fit convergence tolerance (maxiter = 15) Maximum number of fit iterations (nreject = 0) Number of rejection iterations 34

35 (low_reject = 3.) Low sigma rejection factor (high_reject = 3.) High sigma rejection factor (grow = 0.) Rejection growing radius (interactive = yes) Solve fit interactively? (logfile = "STDOUT") Output log file (log_unmatche = yes) Log any unmatched stars? (log_fit = no) Log the fit parameters and statistics? (log_results = no) Log the results? (catdir = "") The standard star catalog directory (graphics = "stdgraph") Output graphics device (cursor = "") Graphics cursor input (mode = "l") cl> fitparams? [ V filter ] :vshow :results # # # (B ) [V ] cl> page calib # FITPARAMS sobs INVERTFIT cl> epar invertfit 35

36 cl> lpar invertfit observations = "sobs" config = "m92fig" parameters = "calib" calib = "results" (catalogs = "stds") (errors = "obserrors") (objects = "all") (print = "") (format = "") List of observations files Configuration file Fitted parameters file Output calibrated standard indices file List of standard catalog files Error computation type (undefined,obserrors,eu Objects to be fit (all,program,standards) Optional list of variables to print Optional output format string (append = no) Append output to an existing file? (catdir = "") The standard star catalog directory (mode = "l") cl> invertfit cl> page results aperture QPHOT aperture MKOBSFILE aperture A User s Guide to Stellar CCD Photometry with IRAF, by Philip Massey and Lindsey E. Davis, April A User s Guide to the IRAF Apphot Package, by Lindsey Elspeth Davis, May Specifications for the Aperture Photometry Package, Lindsey Davis, October Online help is available for all tasks. and "phelp pcintro" in the PHOTCAL package. Especially see "phelp config" end of exercise

37 VI Exercise 5. spectra Exersice 3 cl> cd cl> cd exercises/spec sp0018 sp0027 cl> dir cl> imhead sp*.fits cl> imhead sp0018 l+ #? cl> imhead sp*.fits sp0018.fits[41,1022][real]: comp 6707 start of night sp0020.fits[41,1022][real]: DHCep 6707 sp0021.fits[41,1022][real]: DHCep 6707 sp0022.fits[41,1022][real]: DHCep 6707 sp0023.fits[41,1022][real]: AHCep 6707 sp0024.fits[41,1022][real]: AHCep 6707 sp0025.fits[41,1022][real]: AHCep 6707 sp0027.fits[41,1022][real]: comp 6707 noao.twodspec.longslit cl> imred cl> kpnoslit apall apall FWFM(Full Width Full Maximum) implot cl> implot sp0020 # C FWHM 50 "C" 37

38 "C" 2 [ FWHM ] o # :l 50 # 50 o :l 1000 c # "c" :i sp0025 # :l 511 [ implot ] apall cl> unlearn apall cl> epar apall # : # cl> lpar apall page input = "@extract" nfind = 1 (output = "") (format = "multispec") (references = "") (profiles = "") (interactive = yes) (find = yes) (recenter = yes) (resize = yes) (edit = yes) List of input images Number of apertures to be found automatically List of output spectra Extracted spectra format List of aperture reference images List of aperture profile images\n Run task interactively? Find apertures? Recenter apertures? Resize apertures? Edit apertures? 38

39 (trace = yes) (fittrace = yes) (extract = yes) (extras = no) (review = yes) (line = INDEF) (nsum = 10) (dispaxis = 2) (lower = -5.) (upper = 5.) (apidtable = "") (b_function = "chebyshev") (b_order = 2) (b_sample = "-10:-6,6:10") (b_naverage = -10) (b_niterate = 0) (b_low_reject = 3.) (b_high_rejec = 3.) (b_grow = 0.) (width = 6.) (radius = 10.) (threshold = 0.) (minsep = 5.) (maxsep = 1000.) (order = "increasing") (apertures = "") (npeaks = INDEF) (shift = yes) (llimit = INDEF) (ulimit = INDEF) (ylevel = 0.1) (peak = yes) (bkg = yes) (r_grow = 0.) (avglimits = no) (t_nsum = 10) (t_step = 10) (t_nlost = 3) (t_function = "legendre") Trace apertures? Fit the traced points interactively? Extract spectra? Extract sky, sigma, etc.? Review extractions?\n Dispersion line Number of dispersion lines to sum\n\n# DEFAULT Dispersion axis (1=along lines, 2=along columns Lower aperture limit relative to center Upper aperture limit relative to center Aperture ID table (optional)\n\n# DEFAULT BACKG Background function Background function order Background sample regions Background average or median Background rejection iterations Background lower rejection sigma Background upper rejection sigma Background rejection growing radius\n\n# APERTU Profile centering width Profile centering radius Detection threshold for profile centering\n\n# Minimum separation between spectra Maximum separation between spectra Order of apertures\n\n# RECENTERING PARAMETERS\n Select apertures Select brightest peaks Use average shift instead of recentering?\n\n# Lower aperture limit relative to center Upper aperture limit relative to center Fraction of peak or intensity for automatic wid Is ylevel a fraction of the peak? Subtract background in automatic width? Grow limits by this factor Average limits over all apertures?\n\n# TRACING Number of dispersion lines to sum Tracing step Number of consecutive times profile is lost bef Trace fitting function 39

40 (t_order = 2) (t_sample = "*") (t_naverage = 1) (t_niterate = 0) (t_low_reject = 3.) (t_high_rejec = 3.) (t_grow = 0.) (background = "fit") (skybox = 1) (weights = "none") (pfit = "fit1d") (clean = no) (saturation = INDEF) (readnoise = "0.") (gain = "1.") (lsigma = 4.) (usigma = 4.) (nsubaps = 1) (mode = "l") Trace fitting function order Trace sample regions Trace average or median Trace rejection iterations Trace lower rejection sigma Trace upper rejection sigma Trace rejection growing radius\n\n# EXTRACTION Background to subtract Box car smoothing length for sky Extraction weights (none variance) Profile fitting type (fit1d fit2d) Detect and replace bad pixels? Saturation level Read out noise sigma (photons) Photon gain (photons/data number) Lower rejection threshold Upper rejection threshold Number of subapertures per aperture extras off b order=2( ) b naver=-10 width=6( 6 ) nfind=1 background=fit apall cl> files sp002?.fits > extract cl> edit extract # sp cl> # cl> apall yes (sp0020.fits)? # b # 40

41 ? # t # s # s # s # s # f # :sample # # yes :func spline3 :order 3 :niter 1 f # yes 41

42 # database cl> imhead *.ms* cl> imhead sp0025.ms l+ page cl> dir database cl> page database/apsp0021 # # apsum cl> unlearn apsum cl> epar apsum cl> lpar apsum input = "sp0018,sp0027" List of input images (output = "") List of output spectra (format = "multispec") Extracted spectra format (refinput = "sp0018,sp0027" List of input images (output = "") List of output spectra (format = "multispec") Extracted spectra format (references = "sp0025") List of aperture reference images (profiles = "") List of aperture profile images\n (interactive = no) Run task interactively? (find = no) Find apertures? (recenter = no) Recenter apertures? (resize = no) Resize apertures? (edit = no) Edit apertures? (trace = no) Trace apertures? (fittrace = no) Fit the traced points interactively? (extract = yes) Extract apertures? (extras = no) Extract sky, sigma, etc.? 42

43 (review = yes) (line = INDEF) (nsum = 10) (background = "none") (weights = "none") (pfit = "fit1d") (clean = no) (skybox = 1) (saturation = INDEF) (readnoise = "0.") (gain = "1.") (lsigma = 4.) (usigma = 4.) (nsubaps = 1) (mode = "l") Review extractions?\n Dispersion line Number of dispersion lines to sum\n Background to subtract (none average fit) Extraction weights (none variance) Profile fitting type (fit1d fit2d) Detect and replace bad pixels? Box car smoothing length for sky Saturation level Read out noise sigma (photons) Photon gain (photons/data number) Lower rejection threshold Upper rejection threshold Number of subapertures per aperture cl> apsum cl> imhead *.ms.fits cl> implot sp0018.ms :i sp0027.ms l # docs spec.eps cl>!lpr spec.eps identify reidentify cl> unlearn identify reidentify cl> implot sp0027.ms identify fwidth=4 coordlist linelists$thorium.dat identify database cl> lpar identify cl> epar identify cl> identify sp0027.ms 43

44 ? # w #? # w f # "m" "m" "w" "e" "e" "m" "w" "a" f # ICFIT? # :func cheby # :order 3 # f # # IDENTIFY w # a l # thorium.dat f # ICFIT :nite 2 # f # d # "d" f # # IDENTIFY d # ICFIT RMS 0.01 A f # l # # IDENTIFY # IDENTIFY 44

45 cl> dir database # idsp0027.ms cl> page database/idsp0027.ms sp0027.ms.fits sp0027.ms reidentify cl> reidentify sp0027.ms sp0018.ms \ coord=linelists$thorium.dat v+ inter+ cl> imhead sp0027.ms l+ refspec1 refspectra UT(Universal Time) UT setairmass UT cl> hselect sp*.ms.fits $I,ut,exptime yes cl> lpar setairmass cl> setairmass sp*.ms.fits updatecl> setairmass sp*.ms.fits cl> imhead sp0027.ms l+ refspectra cl> epar refspectra cl> lpar refspectra input = "sp*.ms.fits" answer = "no" (references = "sp*.ms.fits") (apertures = "") (refaps = "") (ignoreaps = yes) (select = "interp") (sort = "utmiddle") (group = "") (time = yes) (timewrap = 17.) (override = no) (confirm = yes) (assign = yes) List of input spectra Accept assignment? List of reference spectra Input aperture selection list Reference aperture selection list Ignore input and reference apertures? Selection method for reference spectra Sort key Group key Is sort key a time? Time wrap point for time sorting Override previous assignments? Confirm reference spectrum assignments? Assign the reference spectra to the input spect 45

46 (logfiles = "STDOUT,logfile") List of logfiles (verbose = no) Verbose log output? (mode = "l") cl> refspectra # "yes" "YES" cl> imhead sp0025.ms l+ # dispcor refspec dispcor cl> epar dispcor cl> lpar dispcor input = "sp*.ms.fits" List of input spectra output = "%sp%l%*.ms.fits" List of output spectra (linearize = yes) Linearize (interpolate) spectra? (database = "database") Dispersion solution database (table = "") Wavelength table for apertures (w1 = INDEF) Starting wavelength (w2 = INDEF) Ending wavelength (dw = INDEF) Wavelength interval per pixel (nw = INDEF) Number of output pixels (log = no) Logarithmic wavelength scale? (flux = yes) Conserve flux? (samedisp = no) Same dispersion in all apertures? (global = yes) Apply global defaults? (ignoreaps = no) Ignore apertures? (confirm = yes) Confirm dispersion coordinates? (listonly = yes) List the dispersion coordinates only? (verbose = yes) Print linear dispersion assignments? (logfile = "") Log file (mode = "l") cl> dispcor (changing assignments) NO cl> dispcor listonly- 46

47 cl> imhead l0025.ms l+ cl> splot l0025.ms # SPLOT SPLOT cl> specplot l0020.ms,l0021.ms,l0022.ms 3 cl> scombine l0020.ms,l0021.ms,l0022.ms spec1 \ combine=median scale=mode cl> splot spec1 o # g l0020.ms o g l0021.ms o g l0022.ms g spec1 splot bad pixels phelp A User s Guide to Reducing Slit Spectra with IRAF, by Phil Massey, Frank Valdes, and Jeannette Barnes, April Online help is available for all tasks. Especially see "phelp onedspec.package" and "phelp onedspec.specwcs". See "phelp doslit" for doing it all in one big gulp, and better too, since it extracts the arcs for each stellar spectrum. 47

48 end of exercise