ESI data reduction - isobel Hook Data are binned 2x1 #Ran esibias on bias frames #esibias @bias.list # not sure this is working for binned data. skip this step. Made a combined bias (no overscan correction) using imcombine. imcombine @bias.list ebias_n1 combine=average reject=avsigclip scale=none # run esibias on the flats - decided not to because it doesn't seem to work for binned data. #esibias @flat.list #ccdproc the flats to take off the bias ccdproc @flat.list ccdtype="" fixpix- oversc- trim- zero+ dark- flat- zero="ebias_n1" #Make a combined flat imcombine @flat.list output=eflat0.7_n1 combine=average reject=avsigclip \ scale=median statsec=[460:500,2150:2200] brun = "1037" run no? cng = "1" night? (biasfile = "ebias_n1") bias file? (flatfile = "eflat0.7_n1") flat file? (reffile = ) trace ref file? (refarc = "e22_23.csx") Arc ref file? (bsens = "home$../testdata/e2049std_2sens") Sens function file? (dobias = yes) Bias correct data? (doflat = yes) Flatten data? (dotrim = yes) Trim data? (dodisp = yes) Display image? (dospec = yes) Reduce as a spectrum? (mode = "q") # arc: imcop n1raw$esi1017 e1017 #copied database file apidata1_e1058 ape1058 apall e1017 out=e1017arc refer=e1058 trace- back- recen- intera- ecreidentify e1017arc reference=e22_23.csx #ecreident doesn;y have an interactive mode. ecidentify e1017arc # reduce the standard brun = "1059" run no? cng = "1" night? (biasfile = "ebias_n1") bias file? (flatfile = "eflat0.7_n1") flat file? (reffile = ) trace ref file? (refarc = "e1017arc") Arc ref file? (bsens = "home$../testdata/e2049std_2sens") Sens function file? (dogreg = no) File from Greg? (dodisp = yes) Display image? (dospec = yes) Reduce as a spectrum? (mode = "ql") use e1058_ww standard n1_1059ww out=e1059std extinct=onedstds$mkoextinct.dat caldir=onedstds$oke1990/ star_nam=bd75325 observa=keck sensfunc e1059std e1059sens extinct=onedstds$mkoextinct.dat observa=keck standard n1_1063ww out=e1063std extinct=onedstds$mkoextinct.dat caldir='/home/esi8/Greg/iraf/spstds/' star_nam=hd93521_stis_001x observa=keck sensfunc e1063std e1063sens extinct=onedstds$mkoextinct.dat observa=keck calibrate n1_1063ww n1_1063c sens=e1063sens observat=keck ## WHEN RUNNING ESI_IS : CAN APPLY A SHIFT TO ALL APERTURES BY ## PRESSING 'a' FIRST ###################################################################' Setup to run Greg/ Lifan's code: .run esi_red .run la_cosspec in my n1 directory idata1='/h/scratch41/esi8/Isobel/n1 cd,idata1 the esi_red,..... Preparation of standard using their code esi_red,'../../night1/esi1063.fits','Eesi1063.fits',flat=flatfile,niter=0,sorder=2 ###################################################################' Object S01-027 Frames 1065, 1066 Reduced with esi_is esi_is 1065 1 biasfile=ebias_n1 flatfile=eflat0.7_n1 reffile=idata1$e1063 refarc=e1017arc bsens=e1063sens re-reduced the next night: esi_is 1065 1 biasfile=ebias_n1 flatfile=eflat0.7_n1 reffile=idata1$e1063 refarc=e1017arc bsens=e1063sens flatfile='/h/scratch41/esi8/lifan/night1/NFlatn1.fits' esi_red,'../../night1/esi1065.fits','Eesi1065.fits',flat=flatfile Looks like we missed the object. ###################################################################' Object S01-016 Frame 1070 Reduced with esi_is z=0.5612 (OII doublet, OIII, Hb, Hg) Possibly some SN features. Try greg's version on 1070 - cant get it to work? floating point overflow. Files: S02-016_snz.ps, n1_1070cc_comb.asc, n1_1070cc_comb.fits ###################################################################' Object S01-035 Frame 1074,1075,1076 Reduced with esi_is z=1.37?? (OII) 'SN' is red, could match maybe. Trace is weaker in 1075, almost invisible in 1076. Combination of all three: n1_S02_35c.asc, n1_S02_35c.ps Try greg's version on 1074 - cant get it to work? floating point overflow. flatfile='/h/scratch41/esi8/lifan/night1/NFlatn1.fits' (normalised, overscan corrected flat with 0.75" slit) esi_red,'../../night1/esi1074.fits','Eesi1074.fits',flat=flatfile ccdproc Eesi1074.fits fixpix+ fixfile='/home/esi8/lifan/iraf/echelle/esimask1.dat' - all other flags to no. Files (from combined spectrum) S02-035c.asc, S02-035c.fits, S02-035c.ps ###################################################################' # NIGHT 2 # ###################################################################' Made a flat using imcombine - no other steps done (no overscan correction etc) = eflat1.0_n2 Object SNF : 2002-031 Frame 2034 Blue objects. Broad p-cygni type lines at z=0: strange star? ###################################################################' Object SNF : 2002-028 Frame 2035 Ran my reduction - got very odd wiggles in the blue. Otherwise featureless. esi_red,'../../night2/esi2035.fits','Eesi2035.fits',flat='eflat1.0_n2.fits',/fproc Its subtracting the object. Need lower order for sky subtraction esi_red,'../../night2/esi2035.fits','Eesi2035.fits',flat='eflat1.0_n2.fits',/fproc,sorder=3 ############################################# Next object S02-049 Frames 2051,2052 esi_red,'../../night2/esi2051.fits','Eesi2051.fits',flat='eflat1.0_n2.fits',/fproc,sorder=3 esi_red,'../../night2/esi2052.fits','Eesi2051.fits',flat='eflat1.0_n2.fits',/fproc,sorder=3 Fixpixed both. Extracted both using apertures from 2051 Ran Lifan's code using flux calibration from tonight with 1,0arcsec slit. Files: S02-049c.asc, S02-049c.fits, S02-049c.ps ############################################# S02-032 Frame 2058 2059 (shortned) Dont see anything in sky-subtracted image esi_red,'../../night2/esi2058.fits','Eesi2058.fits',flat='eflat1.0_n2.fits',/fproc,sorder=3 Maybe something in order 9 - but may also be in the flat - ################################################################## C02-000 Frame 2060 Done with my reduction. Probably z=0.26 at day +11, from SN features only. Files: C02-000_gal.fits, C02-000_gal.asc (host) C02-000_sn.fits, C02-000_sn.asc, C02-000_sn.ps (SN) ################################################################## NIGHT 3 ################################################################## C02-008 Frame 3046 esi_red,'../../night3/esi3046.fits','Eesi3046.fits',flat='../n2/eflat1.0_n2.fits',/fproc,sorder=3 SN should be in the middle of the slit. z=0.28 Type I (from Sn features). No SiII dip or galaxy lines. Files: C02-008c.fits C02-008c.asc, ################################################################## Summary S01-027 No trace in the spectrum in either file. S01-016 z=0.5612 galaxy lines: OII doublet, OIII, Hb, Hg. Possible SN features? S01-035 z=1.37? OII in neighbour object 5" away. Possibly SN features at that z. S02-049 z~0.85?? Maybe a SN bump around 8000A. Needs redoing. Flux in blue though. S02-032 Nothing in the trace on centre. (flat-fielding error to side?) C02-000 z=0.26 Type Ia at day around +11, from SN features only. Si feature visible. C02-008 z=0.28 Type I (from SN features). No SiII dip or galaxy lines.