NB making postscript files: - doesn't seem to work? .snap psdump A file /tmp/irafdmpNNNNN.ps (where NNNNN is some number) will be created. (j$=/home/astro23/deepsearch/julia/reductions/arcs/) l071 is a standard. Isobel extracted standard & copied ap file from database to arc directory (had to rename it to fool IRAF). cd k$is/arcs/ apall j$arcs/l068 out=068arc refer=l071 trace- back- intera- (Hg+Kr+Ne+Ar) apall j$arcs/l078 out=078arc refer=l071 trace- back- intera- (Hg+Kr+Ne+Ar) - j. I used ap#2 from 072 as ref apall j$arcs/l078 out=078arc refer=l072 trace- back- intera- nfind=1 (Hg+Kr+Ne+Ar) - I. apall j$arcs/l084 out=084arc refer=l071 trace- back- intera- (Ne) apall j$arcs/l085 out=085arc refer=l071 trace- back- intera- (Ar) apall j$arcs/l086 out=086arc refer=l071 trace- back- intera- (Hg+Kr) apall j$arcs/l087 out=087arc refer=l071 trace- back- intera- (Hg+Kr+Ne+Ar) apall j$arcs/l088 out=088arc refer=l071 trace- back- intera- (Ne+Ar) All these need extracting. Use object aperture ? Is standard good enough ? 095 Ne+Ar apall l095 out=al095arc refer=l091sig8cr trace- back- intera- nfind=1 096 Hg+Kr apall l096 out=al096arc refer=l091sig8cr trace- back- intera- nfind=1 redid this using ap#2 (peak) of obj 090 as ref - I 095 Ne+Ar apall j$arcs/l095 out=095arc refer=l090 trace- back- intera- nfind=1 096 Hg+Kr apall j$arcs/l096 out=096arc refer=l090 trace- back- intera- nfind=1 106 Hg+Kr apall j$arcs/l106 out=106arc refer=l097 trace- back- intera- nfind=1 107 Ne+Ar apall j$arcs/l107 out=107arc refer=l097 trace- back- intera- nfind=1 125 Hg+Kr apall l125arc out=al125arc refer=l134 trace- back- intera- nfind=1 126 Ne+Ar apall l126arc out=al126arc refer=l134 trace- back- intera- nfind=1 Shouldn't use 135 as ref here. 136 Hg+Kr apall l136arc out=al136arc refer=l135 trace- back- intera- nfind=1 137 Ne+Ar apall l137arc out=al137arc refer=l135 trace- back- intera- nfind=1 I redid it using ap #2 of 141 as ref - I 136 Hg+Kr apall j$l136 out=136arc refer=l141 trace- back- intera- nfind=1 137 Ne+Ar apall j$l137 out=137arc refer=l141 trace- back- intera- nfind=1 151 Hg+Kr apall l151 out=al151arc refer=l145sig8cr trace- back- intera- nfind=1 152 Ne+Ar apall l152 out=al152arc refer=l145sig8cr trace- back- intera- nfind=1 Isobel's version: 151 Hg+Kr apall j$arcs/l151 out=151arc refer=l145 trace- back- intera- nfind=1 152 Ne+Ar apall j$arcs/l152 out=152arc refer=l145 trace- back- intera- nfind=1 Use ap #2 of l162 as ref apall j$arcs/l163 out=163arc refer=l162 trace- back- intera- nfind=1 apall j$arcs/l164 out=164arc refer=l162 trace- back- intera- nfind=1 apall j$arcs/l174 out=174arc refer=l172 trace- back- intera- nfind=1 apall j$arcs/l175 out=175arc refer=l172 trace- back- intera- nfind=1 apall j$arcs/l233 out=233arc refer=l239 trace- back- intera- nfind=1 -I apall j$arcs/l234 out=234arc refer=l239 trace- back- intera- nfind=1 apall j$arcs/l243 out=243arc refer=l241 upper=3.0 lower=-3.0 trace- back- intera- nfind=1 apall j$arcs/l244 out=244arc refer=l241 upper=3.0 lower=-3.0 trace- back- intera- nfind=1 apall j$arcs/l250 out=250arc refer=l248 trace- back- intera- nfind=1 apall j$arcs/l251 out=251arc refer=l248 trace- back- intera- nfind=1 use ap2 0f 248 as trace (edited database/apfile) 259 ] UV setup 260 ] apall j$arcs/l267 out=267arc refer=l273 trace- back- intera- nfind=1 apall j$arcs/l268 out=268arc refer=l273 trace- back- intera- nfind=1 Identification of arcs ---------------------- Used linelist ~imh/iraf/keck_nov95/cuarne Isobels arc spectra and identifications are in directory k$is/arcs/ Julias arc spectra and identifications are in directory s$ (=/home/astro23/deepsearch/keck/spectra) linelist=~imh/iraf/keck_nov95/cuarne combined arcs (233_234arc) combined together using with 'sarith': sarith l233 op=+ l234 out=233_234arc 068 done by hand: 4th order 078 reidentified using 068 as ref: 4th order -legendre - j. 5th order legendre: rms=0.35 - I 087 reidentified using 068 as ref: 5th order -legendre 095_096arc reidentified using 068 sd ref: 4th order -legendre rms=0.2879 -j. 095_096arc reidentified using 068 as ref: 5th order -legendre rms=0.35 -I. 267_268arc " 5th order. rms=0.34 174_175arc 250_251arc reidentified using 267_268arc : 5th order - legendre 243_244arc reidentified using 267_268arc : 106_107arc reidentified using 068arc( isobel's arc) - j. 106_107arc reidentified using 095_096arc rms=0.4, 5th order. - I. 125_126arc reidentified using l106_l107arc: 4th order spline3 - j. 136_137arc reidentified using l125_l126arc: " RMS=0.2403 - j. where is this ? I redid it using 267_268 as ref : 5th order legendre, rms=0.28-I 151_152arc reidentified using l125_l126arc: " legendre rms=0.1512 - j reident using 095_096arc - useless. did by hand: legendre order 5 rms=0.25 - I 163_164arc reidentified using 174_175arc: rms=0.244, 4th order legendre 233_234arc reidentified using 243_244arc as ref. rms=0.35, order=4 ------------------------------------------------------------------------ DISPERSION CORRECTION (Wavelength Calibration) *note before dispcor, put reference spectrum in header using hedit: hedit al070 REFSPEC1 l078arc (using hedit in place of refspec) - j. dispcor al070 dl070 - to wavelength calibrate - j. l078arc - for standards 069-071 (HD19445) - j. 095_096arc - for obj 090-091 (SN95103) - j. 106_107arc - for standards 108,109 (HD84937) -j. 125_126arc - for standards 134,135 (g158-100) - may not use; not in onedstds$ - j. 136_137arc - for 141,142 (SN9570) - j. 151_152arc - for 143-150 (SN 95116, 9571, 9579) - j. 163_164arc - for 157-170 (SN 9568, 95115, 95104) - j. 174_175arc - for standard 171,172 - j. 233_234arc - 243_244arc - for 241,242 (SN 9569) - j. 250_251arc - for 248-256 (SN 9528, 9591) not using 259_260arc for SN9591 cuz UV setup -j. 267_268arc - for 265-273 (SN R42, R47) - j. Standards --------- Have to use very high order (e.g. 25) in sensfunc because of funny shape of flat-field. 2D flat-fielded standard images are in j$star/ Isobel's standards are in /home/astro23/deepsearch/keck/is/stand/ sensl070.imh - hd19445 spline3 order=13 - j. calibrated standard has peak at 9200 - bad. -I I redid it: ap070, W calib with 068arc 24th order legendre. Added points. sensl071.imh - hd19445 spline3 order=13 RMS=0.0275 - j. sensl108.imh - hd84937 exp=1s spline3 order=17 RMS=0.0240 sensl109.imh - hd84937 exp=1s spline3 order=15 RMs=0.0183 sensl134.imh - g158-100 extracted by julia: sensl135.imh - g158-100 " s$dl135. added lots of points and fitted 23 legendre. >> 135sens sensl171.imh - hd84937 sensl172.imh - hd84937 spline3 order=11 (not good enough try higher) - j. 172sens - v. high order legendre (25 ?) sensl173.imh - hd84937 spline3 order=15 rms=0.0159 sens sens108_109 for FLUX CALIBRATION: ---------------- * all stars calibrated with respect to themselves. calib dl071 cl071 sens=sensl071 I decided to sum sensl118 and sensl109 together with an average. - j. imsum sensl108,sensl109 sens108_109 opt=average sensl071 used for calibration of objects l072 to l088 (1st half of night) sens108_109 used for calibration of objects l089 to end of night (2nd half of night) sens173 used for calibration of all objects for night 2 (171 & 172 - 1 sec exposures; 173 - 2s) Extraction of objects --------------------- Flat_fielded data is in Flat-fielded, CR-cleaned data is in /home/astro23/deepsearch/spectra/ moved to /home/sp/fall95k/flt/. moved reductions in /home/astro23/deepsearch/isobel/keck_nov95/reduced to /home/sp/fall95k/rdx/isobel/old/. Isobel's extractions are in /home/astro23/deepsearch/spectra/is/ moved these to /home/sp/fall95k/rdx/isobel/ ascii files in /home/sp/fall95k/rdx/isobel/asc -------------------------------------------------------------------- R47 273 one object only 272 missing - will move -------------------------------------------------------------------- R42 265 Faint object in ap#1, SN is bright object @ row 187, in ap#2 266 " Combined to make apR42 W calib with 267_268arc F calib with k$is/stand/172sens -------------------------------------------------------------------- SN9591 255 4 objects. SN is in ap#2 ? (ap 1 has CR - only use sky to left) 256 " ap 1 OK this time Combined >> ap9591 W calib 250_251arc F calib l172 ap 1 ap 2 SN - type II ? z=0.314 ap 3 bumpy spectrum ? ap 4 wrote as ascii using specconv : ap9591c_ap2.spec wspectext ap9591c[*,1] asc/ap9591c_ap1.asc wspectext ap9591c[*,3] asc/ap9591c_ap3.asc wspectext ap9591c[*,4] asc/ap9591c_ap4.asc wfits ap9591c[*,1] fits/ap9591c_ap1.fits wfits ap9591c[*,2] fits/ap9591c_ap2.fits wfits ap9591c[*,3] fits/ap9591c_ap3.fits wfits ap9591c[*,4] fits/ap9591c_ap4.fits -------------------------------------------------------------------- SN9528 248 4 objects (SN on shoulder of gal) 249 " ap 1 faint gal ap 2 z=0.160 ap 3 SN, clear wiggles ap 4 galaxy z=0.425 Combined >> ap9528 W calib 250_251arc F calib l172 wrote as ascii with specconv: ap9528c_ap3.spec wspectext ap9528c[*,1] asc/ap9528c_ap1.asc wspectext ap9528c[*,2] asc/ap9528c_ap2.asc wspectext ap9528c[*,4] asc/ap9528c_ap4.asc wfits ap9528c[*,1] fits/ap9528c_ap1.fits wfits ap9528c[*,2] fits/ap9528c_ap2.fits wfits ap9528c[*,3] fits/ap9528c_ap3.fits wfits ap9528c[*,4] fits/ap9528c_ap4.fits -------------------------------------------------------------------- SN9569 241 one object- possible shoulder ? difficultish to trace 242 Combined >> ap9569 W calib 243_244arc F calib l172 ap1 SN+g z=0.498 from OII, Ca H&K in galaxy wrote as asc using specconv: ap9569c_ap1.spec wrote as fits: wfits ap9569c fits/ap9569c.fits tried to extract shoulder >> ap9569_2 ap 1 = peak ap 2 = shoulder which is SN ? - ap 2 ?? -------------------------------------------------------------------- SN95126 238 cant see any object in 2D image ! 239 2 objects: bright thing way off, SN near centre (merged with gal) 240 "" W calib = 233_244 F calib 172 (24 pixels outside range) 072 Can't do CR removal because 073 screwed up. - extracted 2 objects as above W calib=078arc F calib: 070sens 073 Short exposure. Jog in frame (missing pixel ?) - dont even try to extract combined 072,239,240 >> ap95126c ap#1 bright gal ? ap#2 SN + host difference in slope between 072 and 239, esp for bright gal. ascii files: ap2 > ap95126c_ap2.spec (using specconv) wspectext ap95126c[*,1,1] asc/ap95126c_ap1.asc wfits ap95126c[*,1,1] fits/ap95126c_ap1.fits wfits ap95126c[*,2,1] fits/ap95126c_ap2.fits -------------------------------------------------------------------- SN95103 169 attempted to extract 2: #1= lh shoulder, #2= peak 170 (wasn't cr removed- near dawn) " W calib: 174_175arc -I F calib: 172sens - I 089 missing - short exp time. 090 extracted 2 objects: #1 = LH shoulder, #2 = peak - I 091 extracted 2 objects: #1 = LH shoulder, #2 = peak - I W calib: 095_096arc (on individual spectra. Add all later) -I F calib: 070sens - I Combined all 4 using scombine & median weighting: >>ap95103c Mar 5 19:04: SCOMBINE combine = average, scale = median, zero = none, weight = median reject = avsigclip, mclip = yes, nkeep = 1 lsigma = 3., hsigma = 3. blank = 0. sample = [6000:8000] Images Median Scale Weight ap169c[ 1] 4.7E-19 1.584 0.158 ap170c[ 1] 3.0E-19 2.448 0.102 ap090c[ 1] 1.2E-18 0.613 0.408 ap091c[ 1] 9.9E-19 0.753 0.332 should give 169, 170 more weight ?? ap#1 SN - sky @5577 ? ap#2 galaxy: OII, Hb, H\& K narrow lines removed: ap95103_nolines Rob wrote as ascii: ap95103c_nolines > 95103.sn.keck.951126.asc I (IMH) wrote as fits wfits ap95103c_nolines[*,1] fits/ap95103c_nolines_ap1.fits wfits ap95103c_nolines[*,2] fits/ap95103c_nolines_ap2.fits -------------------------------------------------------------------- SN95104 161 problems with trace. Something dragging it off = emission line in gal ? Use a reference - try obj 4. Not great. SN is bright object 162 crs not removed - needs moving (near dawn) also difficult to trace. SN ok though ? W calib: used 163_164 arc -I F calib: 172sens - GOT WARNING: W calib is different from 172 ? ap161c: bias WARNING: 32 pixels outside of flux calibration limits Extinction correction applied Flux calibration applied ap162c: bias WARNING: 32 pixels outside of flux calibration limits Extinction correction applied Flux calibration applied combined after flux calibration >> ap95104c ap#1 ap#2 SN - emission lines: z=0.48 OII, OIII Hb (Ha ?) H&K, ap#3 ap#4 ap#2 > ap95104c_ap2.spec. Looks like Rknop copied to 95104.sn.keck.951126.asc wspectext ap95104c[*,1] asc/ap95104c_ap1.asc wspectext ap95104c[*,3] asc/ap95104c_ap3.asc wspectext ap95104c[*,4] asc/ap95104c_ap4.asc wfits ap95104c[*,1] fits/ap95104c_ap1.fits wfits ap95104c[*,2] fits/ap95104c_ap2.fits wfits ap95104c[*,3] fits/ap95104c_ap3.fits wfits ap95104c[*,4] fits/ap95104c_ap4.fits -------------------------------------------------------------------- SN9568 157 - julia has extracted these - no notes - 4 objects extracted. All bright enough to trace. - I - flatfielding problem ?? - I 158 " - 4 objects extracted. All bright enough to trace. - I W calib: 163_164arc F calib: 172sens - WARNING ap157c: bias WARNING: 32 pixels outside of flux calibration limits Extinction correction applied Flux calibration applied ap158c: bias WARNING: 32 pixels outside of flux calibration limits Extinction correction applied Flux calibration applied combined >> ap9568c ap#1 - SN ap#2 - gal @ z=0.66 ? (OII, Ca H & K plus other H lines) ap#3 ap#4 ap1 written to ascii using specconv >> ap9568c_ap1.spec wspectext ap9568c[*,2] asc/ap9568c_ap2.asc wspectext ap9568c[*,3] asc/ap9568c_ap3.asc wspectext ap9568c[*,4] asc/ap9568c_ap4.asc wfits ap9568c[*,1] fits/ap9568c_ap1.fits wfits ap9568c[*,2] fits/ap9568c_ap2.fits wfits ap9568c[*,3] fits/ap9568c_ap3.fits wfits ap9568c[*,4] fits/ap9568c_ap4.fits -------------------------------------------------------------------- SN95115 - julia has extracted these - no notes ? >> sn95115sum. 159 3 objects 160 3 objects W calib: used 163_164arc. - I F calib: used 172sens >> ap95115c. scales look different from Julia's ap#1 APNUM1 = '1 1 133.27 139.27' z=0.204 gal OIII, Hb, Ha ap#2 APNUM2 = '2 2 183.79 189.79' - target ap#3 APNUM3 = '3 3 319.32 325.32' z=0.435 OII, H&K ascii file: ap95115c_ap2.spec (made using specconv) wspectext ap95115c[*,1,1] asc/ap95115c_ap1.asc wspectext ap95115c[*,3,1] asc/ap95115c_ap3.asc wfits ap95115c[*,1,1] fits/ap95115c_ap1.fits wfits ap95115c[*,2,1] fits/ap95115c_ap2.fits wfits ap95115c[*,3,1] fits/ap95115c_ap3.fits -------------------------------------------------------------------- SN9579 149 6 objects: SN is middle of 3 @row 188. - using Ap3 to trace Ap4 -j. - did same thing - I object order got messed up ? - I 150 tracing problems for both of these : objects merge at low pixel end -I W calib: used 151_152arc -I F calib: used 135sens combined >> ap9579c - I #1 2 narrow lines. can't get z #2 galaxy z=0.466 (OII, OIII), H&K #3 SN (sky line @5577 ?) #4 galaxy z=0.464 H&K #5 one narrow line - z=? #6 very red object 28-jul-1997 IMH noticed aperture numbers are out of synch with IRAF's numbering scheme. Re-wrote ascii files so they match the above description. wspectext ap9579c[*,1] asc/ap9579_ap1.asc wspectext ap9579c[*,2] asc/ap9579_ap4.asc wspectext ap9579c[*,3] asc/ap9579_ap2.asc wspectext ap9579c[*,4] asc/ap9579_ap3.asc wspectext ap9579c[*,5] asc/ap9579_ap5.asc wspectext ap9579c[*,6] asc/ap9579_ap6.asc wfits ap9579c[*,1] fits/ap9579_ap1.fits wfits ap9579c[*,2] fits/ap9579_ap4.fits wfits ap9579c[*,3] fits/ap9579_ap2.fits wfits ap9579c[*,4] fits/ap9579_ap3.fits wfits ap9579c[*,5] fits/ap9579_ap5.fits wfits ap9579c[*,6] fits/ap9579_ap6.fits -------------------------------------------------------------------- SN9571 146 4 objects 147 star @row 71, other objects @ row 173, 200(SN?), 275. tracing seems to be better for 147. - using spline3 order 5. -j. 147 has header info missing - add by hand hedit ap147w AIRMASS 1.09 add+ hedit ap147w EXPOSURE 1500 add+ 5 objects - all bright enough to trace - I W calib: used 136_137 F calib: used 135sens - WARNING: 20 pix outside range >>ap9571c ap#1 bright star ap#2 possibly split into 2 ? - this is SN - see Ivan's charts emission line @ 6954.8 - oII ? >> z=0.86 ap#3 emission line @5218- oII ? ap#4 wiggles: M star ? ap#5 wrote ap 2 using specconv ? > ap9571c_ap2.spec wspectext ap9571c[*,1,1] asc/ap9571c_ap1.asc wspectext ap9571c[*,3,1] asc/ap9571c_ap3.asc wspectext ap9571c[*,4,1] asc/ap9571c_ap4.asc wspectext ap9571c[*,5,1] asc/ap9571c_ap5.asc wfits ap9571c[*,1,1] fits/ap9571c_ap1.fits wfits ap9571c[*,2,1] fits/ap9571c_ap2.fits wfits ap9571c[*,3,1] fits/ap9571c_ap3.fits wfits ap9571c[*,4,1] fits/ap9571c_ap4.fits wfits ap9571c[*,5,1] fits/ap9571c_ap5.fits --------------------------------------- tried to extract SN alone - not really separated ?. >> ap9571_2c ?? >> ap9571_2_ap2 >> ap9571_nolines >> ap9571_n1 & ap9571_n2 ap#1 ? (SN/ gal- WHICH IS WHICH ? ap#2 ? -------------------------------------------------------------------- SN 95116 : julia has reduced it but no notes -I 145 1 object - straightforward extraction 144 143 W calib 136_137arc -I F calib: 135sens (g158). >> (combined x files then calibrated) ap95116c -I since then I calibrated individual extractions for s/n calc What did rob write out as ascii ? fits: wfits ap95116c fits/ap95116c.fits wspectext ap95116c asc/ap95116c.asc -------------------------------------------------------------------- SN9570 141 3 objects (1& 2 joined) APNUM1 = '1 1 184.54 188.33' APNUM2 = '2 2 188.32 194.10' APNUM3 = '3 3 321.98 327.98' 142 2 objects (merged- clearly separated for some of spectrum - but hard to trace both) APNUM1 = '1 1 182.77 187.47' APNUM2 = '2 2 187.36 193.19' W calib: 136_137 F calib: 135sens ap#1 - emission line seems stronger in Ap1- is this the SN or gal? flux higher in ap #1 by x3 ap#2 18/11/97 summed ap141c and ap142c together. Ap1 has stringer emission lines OII & OIII wspectext sn9570c_sum[*,1,1] asc/sn9570c_sum_ap1.asc wspectext sn9570c_sum[*,2,1] asc/sn9570c_sum_ap2.asc wspectext sn9570c_sum[*,3,1] asc/sn9570c_sum_ap3.asc wfits sn9570c_sum[*,1,1] fits/sn9570c_sum_ap1.fits wfits sn9570c_sum[*,2,1] fits/sn9570c_sum_ap2.fits wfits sn9570c_sum[*,3,1] fits/sn9570c_sum_ap3.fits -------------------------------------------------------------------- SN95110 97 - missing - screwup. need to do CR removal. 98 " 99 " 4 objects: only 3 extracted - faint galaxy at bottom but CRs make it hard to locate trace W calib:106_107arc F calib: 108 ?? - but no arcs for 108. Use 070. >> ap95110c ap#1 - emission line ap#2 ap#3 - SN. ascii: ap95110c_ap3.spec Sn extraction written to 95110_sn. Rob made ascii file of this:95110.sn.keck.951126.asc. IMH made fits file wfits 95110_sn fits/95110_sn.fits -------------------------------------------------------------------- NB IRAF's convention for naming spectra: e.g: sn95116.143.3001.imh | | / \ type of spectrum aperture (object) number i.e. optimal, unweighted, sky or error N.B starts at 0 specconv ap95110c 3 9570 not done yet specconv ap95116c 1 specconv ap9571c 2 specconv ap9579c 3 specconv ap95115c 2 specconv ap9568c 1 specconv ap95104c 2 specconv ap95103c 1 specconv ap95126c 2 specconv ap9569c 1 specconv ap9528c 3 specconv ap9591c 2