I II 3 III

Size: px

Start display at page:

Download "I 2 1 2 2 2 3 2 II 3 III 6 1 8 1.1.......................... 8 1.1.1...................... 8 1.1.2........................ 10 1.1.3..................."

きみつぐすみい
7 years ago
Views:

1 iraf

2 I II 3 III IV 20 1 login.cl 20 2 sky (?) 20 1

. 16 1.5......................... 18 2 19 2.1........... 19 2.2................... 20 IV 20 1 login.

3 I 1 1.3m GRB [pix] [count] [Å] [count] iraf iraf SAOimageds9 2 iraf :noao.twodspec.apextract.apall input = obj.fits List of input images 3 iraf tex iraf 2

4 (?) II :noao.twodspec.apextract.apall apall ap> epar apall input = obj.fits List of input images (output = out.fits) List of output spectra (apertur= ) Apertures (format = onedspec) Extracted spectra format (referen= ) List of aperture reference images (profile= ) List of aperture profile images input fits output format multispec ondspec 3

fits) List of output spectra (apertur= ) Apertures (format = onedspec) Extracted

5 # DEFAULT APERTURE PARAMETERS (lower = -5.) Lower aperture limit relative to center (upper = 5.) Upper aperture limit relative to center (apidtab= ) Aperture ID table (optional) lower upper lower upper # TRACING PARAMETERS (t_nsum = 10) Number of dispersion lines to sum (t_step = 10) Tracing step (t_nlost= 3) Number of consecutive times profile is lost before quitti (t_funct= legendre) Trace fitting function (t_order= 3) Trace fitting function order (t_sampl= *) Trace sample regions (t_naver= 1) Trace average or median (t_niter= 0) Trace rejection iterations (t_low_r= 3.) Trace lower rejection sigma (t_high_= 3.) Trace upper rejection sigma (t_grow = 0.) Trace rejection growing radius t order 3 (?)legendre 4

6 :go Find apertures for obj? : yes Number of apertures to be found automatically: 1 Resize apertures for obj? : yes Edit apertures for obj? : yes 1: 1 q d 5

7 d m q yes 2: 2 d d 3 ds9 f q yes 4 q out.0001.imh iraf out.0001.pix 6

8 3: 4: projection 7

9 III [Å] fil532.fits cl> imstat fil532.fits STDDEV σ cl> display fil532.fits 1 cl> imexamine SAOimage a FWHM ENCLOSED FWHM q 8

10 imarith :noao.digiphot.daophot.daofind Input image(s) (fil532.fits):fil532.fits Output coordinate file(s) (default: image.coo.?) (default):532.fits FWHM of features in scale units (2.5) (CR or value):fwhm New FWHM of features: 2.5 scale units 2.5 pixels Standard deviation of background in counts (0.) (CR or value):σ New standard deviation of background: 0. counts Detection threshold in sigma (4.) (CR or value): New detection threshold: 4. sigma 0. counts Minimum good data value (INDEF) (CR or value): New minimum good data value: INDEF counts Maximum good data value (INDEF) (CR or value): Input Output FWHM σ Output ** ** fil*.fits 0 [x i : x e, y i : y e ] daofind fil*.fits[x i : x e, y i : y e ] *.fits 9

) (CR or value):σ New standard deviation of background: 0. counts Detection threshold in sigma (4.) (CR or value): New detection threshold: 4. sigma 0.

11 **daofind FWHM ** :noao.digiphot.daophot.datapars FWHM datapars FWHM fwhm sigma :images.imgeom.imshift (p8) file:fil1.fits x +1.5 y -2.3 file:shiftfil1.fits im> imshift fil1.fits Output images: shiftfil1.fits Fractional pixel shift in x: 1.5 Fractional pixel shift in y: -2.3 shiftfil1.fits (p10) shiftfil1.fits shiftfil2.fits shiftfiln.fits 1 10

fits im> imshift fil1.fits Output images: shiftfil1.fits Fractional pixel shift in x: 1.

12 imarith calib.fits ** (?)** calib.dat dat shiftfil1.fits shiftfil2.fits shiftfiln.fits 1.2 II(p3) :noao.twodspec.apextract.apall 1.1.3(p1.1.3) calib.fits II(p3) out.0001.imh id.fits textitapall 11

13 input = calib.fits List of input images (output = id.fits) List of output spectra (apertur= ) Apertures (format = onedspec) Extracted spectra format (referen= out.fits) List of aperture reference images (profile= ) List of aperture profile images (interac= yes) Run task interactively? (find = no) Find apertures? (recente= no) Recenter apertures? (resize = no) Resize apertures? (edit = no) Edit apertures? (trace = no) Trace apertures? (fittrac= yes) Fit the traced points interactively? (extract= yes) Extract spectra? (extras = yes) Extract sky, sigma, etc.? (review = yes) Review extractions? referen out.fits find recente resize edit trace no :go Write apertures for calib to database (yes): Extract aperture spectra for calib? (yes): Review extracted spectra from calib? (yes): Review extracted spectrum for aperture 1 from calib? (yes): yes id.fits out.fits calib.fits q id.fits? id.0001.imh iraf 12

(resize = no) Resize apertures? (edit = no) Edit apertures? (trace = no) Trace apertures? (fittrac= yes) Fit the traced points interactively? (extract= yes) Extract spectra?

14 1.3 :noao.twodspec.longslit.identify lo> identify id.0001.imh 5: 5 m Å l 6 f ( spline3 ) 7 13

15 6: 7: 14

16 8: 9: 15

17 h 8 q 9 q Write feature data to the database (yes)? yes identify emission absorption emission absorption **1.1.3 (?)** lo> epar identify (coordli= calib.dat) User coordinate list coordli dat 1.4 :out.0001.imh :id.0001.imh :noao.onedspec.refspectra on> epar refspectra 16

18 input = out.0001.imh List of input spectra (referen= id.0001.imh) List of reference spectra (apertur= ) Input aperture selection list (refaps = ) Reference aperture selection list (ignorea= yes) Ignore input and reference apertures? (select = interp) Selection method for reference spectra (sort = none) Sort key (group = none) Group key (time = no) Is sort key a time? (timewra= 17.) Time wrap point for time sorting (overrid= no) Override previous assignments? (confirm= yes) Confirm reference spectrum assignments? (assign = yes) Assign the reference spectra to the input spectr (logfile= STDOUT,logfile) List of logfiles (verbose= no) Verbose log output? answer = Accept assignment? (mode = ql) input referen sort group none :go [out.0001] refspec1= id.0001 Accept assignment? (no yes YES): yes 17

) Time wrap point for time sorting (overrid= no) Override previous assignments? (confirm= yes) Confirm reference spectrum assignments?

19 1.5 :noao.onedspec.dispcor on> epar dispcor input = out.0001.imh List of input spectra output = spect.fits List of output spectra (lineari= yes) Linearize (interpolate) spectra? (databas= 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? (samedis= no) Same dispersion in all apertures? (global = no) Apply global defaults? (ignorea= no) Ignore apertures? (confirm= no) Confirm dispersion coordinates? (listonl= no) List the dispersion coordinates only? (verbose= yes) Print linear dispersion assignments? (logfile= ) Log file (mode = ql) input output :go spect.fits 18

(nw = INDEF) Number of output pixels (log = no) Logarithmic wavelength scale? (flux = yes) Conserve flux? (samedis= no) Same dispersion in all apertures? (global = no) Apply global defaults?

20 2 2.1 :noao.onedspec.longslit.splot 2 on> splot spect.fits /reffig:spect1 10: z c k k ( ) g ( x FWHM eqw 2 twodspec 19

21 ) e e q iraf 2.2 sarith imarith IV 1 login.cl login.cl! iraf UNIX rm /textitlogin.cl package user task $rm = $foreign OK 2 sky (?) 1.3mφ GRB sky sky :noao.imred.generic.background [x i : x e, y i : y e ] 20

22 (sample = x i : xe, y i : y e ) Sample of points to use in fit sky 21

longslit.dvi

longslit.dvi IRAF [ ] HDS A 1 HDS A.1 HDS p4 1 1 A.2 p9 p13 p19 A.3 HH46-47 2 1 1: 2: 2 A.4 p7-8 IRAF ecl> imhead *.fits HDSA ID FLAT, BIAS, COMPARISON HH46-47 HDS 2 CCD ID fits p9 CCD ID B CCD p9-p11 IRAF overscan.cl