main.dvi

/MOIRCS 2010 12 14

3 1 5 2 MOIRCS MOS 7 2.1 MOIRCS........................................ 8 2.2.................................... 11 2.3 MOIRCS MOS.................................. 12 2.3.1..................................... 13 2.3.2................................... 14 2.3.3.................................. 14 3 17 3.1 MOIRCS MOS........................... 18 3.2 MCSMDP........................................ 20 3.2.1 IRAF PyRAF................................ 20 3.2.2.............................. 21 3.2.3.................................. 22 3.2.4................................... 22 4 27 4.1...................................... 28 4.2......................................... 29 4.2.1.................... 29 4.2.2 A-B................................. 32 4.2.3................. 32 4.2.4................................... 34 4.3.................................... 35 4.4......................................... 36 4.4.1..................................... 36 4.4.2................................. 46 4.5....................................... 47 4.6..................................... 49

4 5 53 5.1.................................... 54 5.2....................... 55 A MCSMDP 59 A.1 MCSMDP............................ 60 A.2....................................... 61 A.2.1................................ 61 A.2.2........................... 61 A.2.3..................................... 62 A.2.4................................ 63 A.2.5..................................... 63 A.2.6......................... 64 A.2.7 1..................... 65

1 5

6 1 MOIRCS (MOS) IRAF (Image Reduction and Analysis Facility) python PyRAF UNIX PyRAF IRAF (R 500 2000) PyRAF MCSMDP (MOIRCS MOS Data Pipelines) MCSMDP 2010 S07A-083 (SMOKA ) PI ( ) (MDP Yoshikawa et al. 2010 [5] 2010 2010 12 14

2 MOIRCS MOS 7

8 2 MOIRCS MOS 2.1 MOIRCS MOIRCS (Multi-Object InfraRed Camera and Spectrograph) (0.9 2.5µm) YJHK s 4 7 VPH (R 700) (R 1500) (R 3000) 2.1 MOIRCS 40 FOCAS (Faint Object Camera and Spectrograph) MOIRCS 2µm 100K a 2.1: MOIRCS [µm] [R] b [Å/pixel] zj500 0.9 1.78 700 5.57 HK500 1.3 2.5 640 7.72 H 820 K c R1300 1.16 1.34 1500 1.91 J 4 c 1.45 1.80 1600 2.61 H 3 c 2.00 2.40 1500 3.88 K 2 d VPH J 1.23 3050 0.96 d VPH H 1.65 2940 1.31 a Web b 0. 5 c J, H, K 4 3 2 J K d Web MOIRCS HAWAII-2 HAWAII-2 2048pixel 2048pixel HgCdTe MOIRCS 0. 117/pixel MOIRCS

2.1. MOIRCS 9 2.1: MOIRCS x 1 x 2 VPH-J VPH-H 2.1 0. 5 ( 300pixel) 4 3. 5 2 2.1 2 (PA) 0 1 2 x y 1 x 2 VPH-J VPH-H MOIRCS (Suzukietal. 2008[4]) Web (http://naoj.org/observing/instruments/moircs/index.html)

10 2 MOIRCS MOS MOIRCS FITS (Flexible Image Transport System) 1 1 2 FITS ID ( ID).fits ID MCSA00057429 MCS MOIRCS A 8 MOIRCS 1 2 FITS FITS FITS 3 MOIRCS FITS ID FRAMEID ID MOIRCS FITS MOS 2.2 FITS IRAF/PyRAF FITS FITS FITS FITS FITS (http://www.fukuoka-edu.ac.jp/ kanamitu/fits/index.html)

2.2. 11 2.2: MOIRCS FITS FRAMEID EXP-ID DET-ID OBS-MOD DATA-TYP OBJECT SLIT DISPERSR EXPTIME K DITWID K DITCNT DATE-OBS HST-STR UT-STR MCSA00057430 ID MCSA00057429 ID 1 2 ID 1 2 SPEC MOS IMAG SPEC OBJECT DARK, DOMEFLAT, INSTFLAT CDFN MASK01 CDFN1 HK500 900.050 3.000 1 1 2 2010-12-14 UT 22:25:30.234 08:25:30.234 UT 2.2 1 FOCAS MOIRCS FOCAS Suprime-Cam CCD (Charge Coupled Device) 1µm MOIRCS HAWAII-2 HgCdTe MOIRCS HAWAII-2 CCD (H 2 O) (CO 2 ) 2.2 MOIRCS 1 MOIRCS 0.9 2.5µm

12 2 MOIRCS MOS HK500 1.3 2.5µm H K 1.9µm 2.2 MOIRCS OH O 3 FMOS OHS (OH Suppressor) OH saturate 10 15 A-B IRCS COMICS 3µm MOIRCS cold stop 2.5µm 8m Ian McLean Electronic Imaging in Astronomy; Detectors and Instrumentation [2] 2.3 MOIRCS MOS MOIRCS MOS FOCAS MOIRCS MOS

2.3. MOIRCS MOS 13 transmission 1.0 0.8 0.6 0.4 0.2 MOIRCS Ks MOIRCS H MOIRCS J atmosphere J H K s 9.0 8.0 7.0 6.0 5.0 4.0 3.0 2.0 A] [erg/s/cm 2 / F 17 10 1.0 0.0 0.0 10000 12000 14000 16000 18000 20000 22000 24000 [ A] 2.2: : MOIRCS J: H: K s : : MOIRCS 2.3.1 MOIRCS MOIRCS MOS (PA) PI PA PI 2 mdp FOCAS mdp wmdp moircs MDP wmdp moircs sbr 2 MOS

14 2 MOIRCS MOS 100K 2.3.2 6 2.3 2 3 1 2 2.3.3 6pixel 3 4

2.3. MOIRCS MOS 15 OH OH A0V MOS

16 2 MOIRCS MOS sky image sky + mask image sky + mask image Mask install Telescope shift/rotate Detect star positions Detect hole positions (Calcuration of shift/rotate) Rough pointing Telescope shift MOIRCS PA rotation detector 2 Spectroscopy start! detector 1 Precision pointing (Alignment stars come to each hole center) 2.3: MOS 2006,

3 17

18 3 3.1 MOIRCS MOS MOS (Hα) Hα 2 S/N x y specific flux (F λ [ergs 1 cm 2 Å 1 ]) IRAF MOS 3.1 A-B x y OH

3.1. MOIRCS MOS 19 3.1: MOIRCS MOS

20 3 A-B OH A-B 2 (F λ [ergsec 1 cm 2 Å 1 ]) 3.2 MCSMDP MCSMDP (MOIRCS MOS Data Pipelines) MOIRCS MOIRCS MOS PyRAF 1 Web 3.2.1 IRAF PyRAF IRAF (Image Reduction and Analysis Facility) NOAO (National Optical Astronomical Observatories) NOAO 2 IRAF (cl ) ( ) NOAO IRAF PyRAF IRAF Python NASA STScI (Space Telescope Science Institute) 3 PyRAF cl IRAF PyRAF Python Python IRAF MCSMDP PyRAF PyRAF MCSMDP IRAF version 2.14.1 (lnux ) PyRAF version 1.9 4 PyRAF 1 http://www.cc.kyoto-su.ac.jp/ tomohiro/mcsmdp/ 2 http://iraf.nao.ac.jp/iraf/web/ 3 http://www.stsci.edu/resources/software hardware/pyraf 4 2010 11 IRAF version 2.15 PyRAF

3.2. MCSMDP 21 Ubuntu Linux 10.4 IRAF/PyRAF 3.2.2 MCSMDP (Ubuntu Linux *-dev ) (apt, yum ) IRAF, Python stsci python stsci python PyRAF python STScI Web Web IRAF, Python PyRAF ds9, xpa Harvard-Smithonian Center for Astrophysics Web 5 ximtools RO Python Package RusselOwen 6 Python Python setuptools easy install matplotlib python 7 Linux scipy python 8 Linux ipython python PyRAF PyRAF 9 Linux R, rpy python 1011 Linux 5 http://hea-www.harvard.edu/rd/ds9/ 6 http://www.astro.washington.edu/users/rowen/ropackage/overview.html 7 http://matplotlib.sourceforge.net/ 8 http://www.scipy.org 9 http://ipython.scipy.org/moin 10 http://www.r-project.org/ 11 http://rpy.sourceforge.net/

22 3 3.2.3 MCSMDP MDPDB MCSMDP 12 home $ wget http://www. cc.kyoto su.ac. jp/ tomohiro/mcsmdp/mcsmdp v1 0 1. tgz $ wget http://www. cc.kyoto su.ac. jp/ tomohiro/mcsmdp/mdpdb v1 0 1. tgz $ tar xvfz MCSMDP v1 0 1. tgz $ tar xvfz MDPDB v1 0 1. tgz MCSMDP setup.sh $ cd MCSMDP v1 0 1 $./ setup.sh 13 bash, ksh, sh, zsh.bash profile (bash ).profile (ksh, sh, zsh ) source ${HOME}/MCSMDP v1 0 1/etc/mcsmdp.sh tcsh, csh.cshrc source ${HOME}/MCSMDP v1 0 1/etc/mcsmdp. csh 3.2.4 mcsmdp mcsmdp iraf mkiraf 14 pyraf (-->) $ mcsmdp mkiraf? (yes ):... > 12 $ bash 13 14 mkiraf iraf mkiraf no

3.2. MCSMDP 23 mcsmdp mcsmdp 15 > mcsmdp... > pyraf IRAF cl pyraf STScI The PyRAF Tutrial 16 mdpdisplay > lpar mdpdisplay file = HK500 MODS11 0390fl. fits Image ( list ) to be loaded (frame = 1) Display frame to be loaded (multiframe = no) Increment frame number for image list? ( scale = linear ) scale type (smode = zscale ) scale mode (z1 = INDEF) minimum greylevel to be displayed (z2 = INDEF) maximum greylevel to be displayed (cmap = BB ) color map ( invert = no) invert color map? (mode = al ) > epar mdpdisplay 3.2 () iraf pyraf epar Save & Quit pyraf Execute Unlearn Cancel Mdpdisplay Help MCSMDP iraf 15 iraf loginuser.cl pyraf 16 http://stsdas.stsci.edu/stsci python epydoc/docs/pyraf tutorial.pdf

24 3 3.2: mdpdisplay pyraf = 1 > mdpdisplay MCSA00057147. fits frame=2 mcsmdp.exit >. exit $ ipython mcsmdp --ipython pyraf ipython $ mcsmdp ipython... In [1]:

3.2. MCSMDP 25 ipython /.ipython/ipythonrc-pyraf ipython mcsmdp Ctrl+D import mod pylab execfile pylabinit.py ipython matplotlib python pyraf IRAF cl pyraf

4 27

28 4 4.1 1 $ mkdir p /mfs02a/tomohrys/mcsmdp work # $ cd /mfs02a/tomohrys/mcsmdp work $ cp /mfs01b/tomohrys/mcsmdp sample/.. mdp, fits 2 mdp MDP fits FITS mcsmdp hselect $ mcsmdp > mcsmdp > hselect. fits $I,OBS MOD,DATA TYP,OBJECT,DISPERSR yes MCSA00057015. fits SPEC MOS OBJECT DOMEFLAT HK500 MCSA00057016. fits SPEC MOS OBJECT DOMEFLAT HK500 MCSA00057017. fits SPEC MOS OBJECT DOMEFLAT HK500... MCSA00057114. fits SPEC OBJECT M53735(A0V:J8.9:H8.9:K8.9) HK500 MCSA00057116. fits SPEC OBJECT M53735(A0V:J8.9:H8.9:K8.9) HK500 MCSA00057147. fits SPEC MOS OBJECT CDFN MASK02 HK500 MCSA00057148. fits SPEC MOS OBJECT CDFN MASK02 HK500 MCSA00057149. fits SPEC MOS OBJECT CDFN MASK02 HK500... OBS-MOD (SPEC MOS) (SPEC) DATA-TYP DOMEFLAT OBJECT DISPERSR HK500 hselect 2 > hselect MCSA. fits $I OBJECT = DOMEFLAT & @ DET ID = 1 > flat1. lst > hselect MCSA. fits $I OBJECT = DOMEFLAT & @ DET ID = 2 > flat2. lst 1 http://www.cc.kyoto-su.ac.jp/ tomohiro/mcsmdp/mcsmdp sample.tbz 2

4.2. 29 > hselect MCSA. fits $I OBJECT = CDFN MASK02 & @ DET ID = 1 & K DITCNT = 1 > obj1a. lst > hselect MCSA. fits $I OBJECT = CDFN MASK02 & @ DET ID = 1 & K DITCNT = 2 > obj1b. lst > hselect MCSA. fits $I OBJECT = CDFN MASK02 & @ DET ID = 2 & K DITCNT = 1 > obj2a. lst > hselect MCSA. fits $I OBJECT = CDFN MASK02 & @ DET ID = 2 & K DITCNT = 2 > obj2b. lst flat1, flat2 1, 2 obj1a, obj1b 1 A B obj2a, obj2b 2 2 1 1 2 4.2 4.2.1 craverage 0 1 1 BPM.pl sed UNIX BPM >! sed s /\(. \)\. fits /BPM\/\1\.pl/ obj1a. lst > bpm1a. lst >! sed s /\(. \)\. fits /BPM\/\1\.pl/ obj1b. lst > bpm1b. lst > mkdir BPM MCSMDP 1 mdpdb$bpm/nlbpm1 FF64r.fits 2 mdpdb$bpm/nlbpm2 FF64r.fits craverage > imcopy mdpdb$bpm/nlbpm1 FF64r. fits,mdpdb$bpm/nlbpm1 FF64r. fits, mdpdb$bpm/nlbpm1 FF64r. fits,mdpdb$bpm/nlbpm1 FF64r. fits @bpm1a. lst > imcopy mdpdb$bpm/nlbpm1 FF64r. fits,mdpdb$bpm/nlbpm1 FF64r. fits, mdpdb$bpm/nlbpm1 FF64r. fits,mdpdb$bpm/nlbpm1 FF64r. fits @bpm1b. lst

30 4 epar craverage input = output = (crmask = ) (average = ) (sigma = ) List of input images List of output images List of output cosmic ray and object masks List of output block average filtered images List of output sigma images (navg = 15) (nrej = 100) (nbkg = 5) (nsig = 10) (var0 = 0.0) (var1 = 0.0) (var2 = 0.0) Block average box size Number of high pixels to reject from the average Background annulus width Box size for sigma calculation Variance coefficient for DNˆ0 term Variance coefficient for DNˆ1 term Variance coefficient for DNˆ2 term (crval = 1) Mask value for cosmic rays ( lcrsig = 100.0) Low cosmic ray sigma outside object ( hcrsig = 10.0) High cosmic ray sigma outside object (crgrow = 0.0) Cosmic ray grow radius (objval = 0) Mask value for objects ( lobjsig = 10.0) Low object detection sigma (hobjsig = 5.0) (objgrow = 0.0) High object detection sigma Object grow radius craverage > craverage @obj1a. lst crmask=@bpm1a. lst > craverage @obj1b. lst crmask=@bpm1b. lst fixpix Y cr >! sed s/\(. \)/ cr\1/ obj1a. lst > crobj1a. lst >! sed s/\(. \)/ cr\1/ obj1b. lst > crobj1b. lst fixpix > imcopy @obj1a. lst @crobj1a. lst > imcopy @obj1b. lst @crobj1b. lst

4.2. 31 epar fixpix images = masks = List of images to be fixed List of bad pixel masks ( linterp = INDEF ) Mask values for line interpolation (cinterp = 1 ) (verbose = no) Mask values for column interpolation Verbose output? ( pixels = no) List pixels? fixpix > fixpix @crobj1a. lst @bpm1a. lst > fixpix @crobj1b. lst @bpm1b. lst MCSMDP MCSMDP crrejection 1. ( 1) 2. A-B A A/B 2 3. 1, 2, OR 4. 5. 2 1 crrejection inimage1 = inimage2 = outimage1 = outimage2 = bpm = (navg = 15) input frame at A position input frame at B position output frame at A position output frame at B position badpixel mask Block average box size

32 4 (nrej = 100) (nbkg = 5) (nsig = 10) (var0 = 0.0) (var1 = 0.0) (var2 = 0.0) Number of high pixels to reject from the average Background annulus width Box size for sigma calculation Variance coefficient for DNˆ0 term Variance coefficient for DNˆ1 term Variance coefficient for DNˆ2 term ( lcrsig = 100.0) Low cosmic ray sigma outside object ( hcrsig = 10.0) High cosmic ray sigma outside object (crgrow = 0.0) Cosmic ray grow radius (bpmdir = BPM ) directory name to store bad pixel masks crrejection ( cr2* BPM2 ) >! sed s/\(. \)/ cr2\1/ obj1a. lst > cr2obj1a. lst >! sed s/\(. \)/ cr2\1/ obj1b. lst > cr2obj1b. lst > crrejection @obj1a. lst @obj1b. lst @cr2obj1a. lst @cr2obj1b. lst mdpdb$bpm/ nlbpm1 FF64r. f i t s bpmdir= BPM2 4.2.2 A-B A B imarith ab >! sed s/\(. \)/ab\1/ crobj1a. lst > abobj1a. lst > imarith @crobj1a. lst @crobj1b. lst @abobj1a. lst 4.2.3 imcombine input = output = (headers = ) (bpmasks = ) (rejmasks = ) (nrejmasks = ) (expmasks = ) (sigmas = ) (imcmb = $I ) List of images to combine List of output images List of header files (optional) List of bad pixel masks (optional) List of rejection masks (optional) List of number rejected masks (optional) List of exposure masks (optional) List of sigma images (optional) Keyword for IMCMB keywords

4.2. 33 ( logfile = STDOUT ) Log file (combine = median ) Type of combine operation ( reject = sigclip ) Type of rejection ( project = no) Project highest dimension of input images? (outtype = real ) Output image pixel datatype ( outlimits = ) Output limits (x1 x2 y1 y2...) ( offsets = none ) Input image offsets (masktype = none ) Mask type (maskvalue = 0 ) Mask value (blank = 0.0) Value if there are no pixels ( scale = exposure ) Image scaling (zero = none ) Image zero point offset (weight = none ) Image weights ( statsec = ) Image section for computing statistics (expname = EXPTIME ) Image header exposure time keyword ( lthreshold = INDEF) Lower threshold (hthreshold = INDEF) (nlow = 1) (nhigh = 1) (nkeep = 1) (mclip = yes) (lsigma = 3.0) (hsigma = 3.0) Upper threshold 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 (rdnoise = 0. ) ccdclip : CCD readout noise ( electrons ) (gain = 1. ) ccdclip : CCD gain ( electrons/dn) ( snoise = 0. ) ccdclip : Sensitivity noise ( fraction ) ( sigscale = 0.1) Tolerance for sigma clipping scaling corrections ( pclip = 0.5) pclip : Percentile clipping parameter (grow = 0.0) imcombine Radius ( pixels ) for neighbor rejection > imcombine @flat1. lst HK500 CDFN2 Domeflat1. fits 1 > imarith HK500 CDFN2 Domeflat1. fits / 10000. HK500 CDFN2 Domeflat1. fits

34 4 fixpix > fixpix HK500 CDFN2 Domeflat1. fits mdpdb$bpm/nlbpm1 FF64r. fits A-B fl >! sed s/\(. \)/ fl \1/ abobj1a. lst > flobj1a. lst > imarith @abobj1a. lst / HK500 CDFN2 Domeflat1. fits @flobj1a. lst 4.2.4 geotran input = output = database = transforms = Input data Output data Name of GEOMAP database file Names of coordinate transforms in database file (geometry = geometric ) Transformation type ( linear, geometric) (xin = INDEF) (yin = INDEF) X origin of input frame in pixels Y origin of input frame in pixels ( xshift = INDEF) X origin shift in pixels ( yshift = INDEF) Y origin shift in pixels (xout = INDEF) (yout = INDEF) (xmag = INDEF) (ymag = INDEF) (xrotation = INDEF) (yrotation = INDEF) (xmin = INDEF) (xmax = INDEF) (ymin = INDEF) X origin of output frame in reference units Y origin of output frame in reference units X scale of input picture in pixels per reference unit Y scale of input picture in pixels per reference unit X axis rotation in degrees Y axis rotation in degrees Minimum reference x value of output picture Maximum reference x value of output picture Minimum reference y value of output picture

4.3. 35 (ymax = INDEF) Maximum reference y value of output picture ( xscale = 1.0) X scale of output picture in reference units per pixel ( yscale = 1.0) Y scale of output picture in reference (ncols = INDEF) units per pixel Number of columns in the output picture ( nlines = INDEF) Number of lines in the output picture (xsample = 1.0) (ysample = 1.0) Coordinate surface sampling interval in x Coordinate surface sampling interval in y ( interpolant = linear ) Interpolant (boundary = constant ) Boundary extension (nearest, constant, (constant = 0.0) (fluxconserve = yes) (nxblock = 512) (nyblock = 512) (verbose = yes) reflect,wrap) Constant boundary extension Preserve image flux? X dimension of working block size in pixels Y dimension of working block size in pixels Print messages about the progress of the task geotran gc >! sed s/\(. \)/gc\1/ flobj1a. lst > gcobj1a. lst > geotran @flobj1a. lst @gcobj1a. lst mdpdb$geomap/mcsdistcrr1 feb07new.dbs mcsdistcrr1 feb07new.gmp geotran MDPDB mdpdb$geomap/mcsdistcrr1 feb07new.dbs mcsdistcrr1 feb07new.gmp 1 mdpdb$geomap/mcsdistcrr2 feb07new.dbs mcsdistcrr2 feb07new.gmp 2 database transform 3 4.3 3 MCSRED

36 4 MCSMDP maskplot MDP > maskplot CDFN MASK02.mdp image=gcflabcrmcsa00057147. fits raw+ MODS11-0390 imcopy >! sed s /\(. \)\. fits /\1 MODS11 0390\. fits / gcobj1a. lst > gcmods11 0390. l s t >! sed s/\(. \)/\1[,1755:1840]/ gcobj1a. lst > cut. lst > imcopy @cut. lst @gcmods11 0390. lst 4.4 4.4.1 X OH A-B >! sed s/\(. \)/ flsky \1/ crobj1a. lst > flsky1a. lst > imarith @crobj1a. lst / HK500 CDFN2 Domeflat1. fits @flsky1a. lst >! sed s/\(. \)/gc\1/ flsky1a. lst > gcsky1a. lst > geotran @flsky1a. lst @gcsky1a. lst mdpdb$geomap/mcsdistcrr2 feb07new.dbs mcsdistcrr2 feb07new.gmp >! sed s /\(. \)\. fits /\1 MODS11 0390\. fits / gcsky1a. lst > gcskymods11 0390. lst >! sed s/\(. \)/\1[,1755:1840]/ gcsky1a. lst > cut. lst > imcopy @cut. lst @gcskymods11 0390. lst identify identify X images = Images containing features to be identified

4.4. 37 crval = cdelt = Approximate coordinate (at reference pixel ) Approximate dispersion ( section = middle line ) Section to apply to two dimensional (database = database ) images Database in which to record feature data ( coordlist = mdpdb$ohlist/list NS HK500 ) User coordinate list (units = ) (nsum = 20 ) (match = 3.0) (maxfeatures = 50) (zwidth = 100.0) (ftype = emission ) (fwidth = 8.0) (cradius = 5.0) (threshold = 0.0) (minsep = 2.0) (function = chebyshev ) (order = 4) (sample = ) Coordinate units Number of lines/columns/bands to sum in 2D images Coordinate list matching limit Maximum number of features for automatic identification Zoom graph width in user units Feature type Feature width in pixels Centering radius in pixels Feature threshold for centering Minimum pixel separation Coordinate function Order of coordinate function Coordinate sample regions ( niterate = 10) Rejection iterations ( low reject = 3.0) Lower rejection sigma ( high reject = 3.0) Upper rejection sigma (grow = 0.0) (autowrite = no) (graphics = stdgraph ) (cursor = ) (aidpars = ) Rejection growing radius Automatically write to database Graphics output device Graphics cursor input Automatic identification algorithm parameters section nsum coordlist OH Rousselot et al. (2000) [3] HK500 MCSMDP 4.1 identify > identify gcflskycrmcsa00057147 MODS11 0390. fits

38 4 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 x1e+4 14000 16000 18000 20000 22000 Angstrom 14356.719 14518.961 14698.437 14887.699 15055.55 15187.14 15240.954 15287.789 15332.403 15432.156 15540.328 15656.963 15833.272 15972.596 16030.832 16128.609 16235.377 16502.365 16692.38 16903.679 16955.078 17008.757 17123.659 17247.926 17386.696 17449.967 17653.222 17880.299 17993.962 18118.494 19250.306 19350.119 19771.862 20008.163 20275.839 20339.497 20412.68 20499.364 20563.548 20729.015 20909.569 21115.857 21176.557 21249.592 21507.308 21711.17 21802.312 21873.518 21955.638 22052.366 22124.875 OH line list for HK500 grism (mdpdb$ohlist/list_ns_hk500) 22311.799 22517.969 4.1: HK500 OH

4.4. 39 4.2 X 20 4.1 4.1: identify m d c w f l q identify? b t x y X Y? m 4.1 5 m 4 5 f 4.3 residual 3 Chebychev 5 identify iraf 4.2 4 H 4 1 X 5 order=4 iraf order order-1

40 4 4.2: identify 4.2: identify h,i,j,k,l f d u c s z :order n n-1 :function :niterate :low reject :high reject q sigma rejection sigma rejection sigma rejection?

4.4. 41 4.3: identify q l f HK500 7.72[Å/pixel] RMS Å q q identify Write feature data to the database (yes)? return database (database/idgcflskycrmcsa00057147 MODS11-0390) reidentify identify reidentify

42 4 reference = Reference image images = Images to be reidentified answer = yes Fit dispersion function interactively? crval = Approximate coordinate (at reference pixel ) cdelt = Approximate dispersion ( interactive = no ) Interactive fitting? ( section = middle line ) Section to apply to two dimensional images (newaps = yes) Reidentify apertures in images not in reference? (override = yes) Override previous solutions? ( refit = yes) Refit coordinate function? (trace = yes) Trace reference image? (step = 20 ) Step in lines/columns/bands for tracing an image (nsum = 20 ) Number of lines/columns/bands to sum ( shift = 0. ) Shift to add to reference features (INDEF to search) (search = 0.0) Search radius ( nlost = 100) Maximum number of features which may be lost (cradius = 5.0) Centering radius (threshold = 0.0) Feature threshold for centering (addfeatures = no) Add features from a line list? ( coordlist = mdpdb$ohlist/list NS HK500 ) User coordinate list (match = 3.0) Coordinate list matching limit (maxfeatures = 50) Maximum number of features for automatic identification (minsep = 2.0) Minimum pixel separation (database = database ) Database ( logfiles = logfile ) List of log files ( plotfile = ) Plot file for residuals (verbose = yes) Verbose output? (graphics = stdgraph ) Graphics output device (cursor = ) Graphics cursor input

4.4. 43 (aidpars = ) Automatic identification algorithm parameters > reidentify gcflskycrmcsa00057147 MODS11 0390. fits gcflskycrmcsa00057147 MODS11 0390. f i t s 20 6 X,Y 2 X Y fitcoords fitcoords images = (fitname = ) Images whose coordinates are to be fit Name for coordinate fit in the database ( interactive = yes) Fit coordinates interactively? (combine = no) (database = database ) Combine input coordinates for a single fit? Database ( deletions = deletions.db ) Deletion list file (not used if null ) (function = chebyshev ) (xorder = 4) (yorder = 3) ( logfiles = STDOUT, logfile ) Log files Type of fitting function ( plotfile = plotfile ) Plot log file (graphics = stdgraph ) (cursor = ) X order of fitting function Y order of fitting function Graphics output device Graphics cursor input function xorder identify, reidentify > fitcoords gcflskycrmcsa00057147 MODS11 0390 identify fitcoords 4.3 fitcoords X Y 6 interactive=yes

44 4 4.3: fitcoords x[key],y[key] r f :order n X Y key n-1 :function q fitcoords? x y z s r [key] X Y (s-z) xxyy r X Y identify reidentify X x Y ( xxyr ) xorder RMS 2 X y Y xyyr Y q fitcoords 1 (MCSA00057147) MOIRCS 3 7 1 fitcoords transform transform input = output = fitnames = (minput = ) Input images Output images Names of coordinate fits in the database Input masks 7 MCSMDP

4.4. 45 4.4: 2 : (moutput = ) (database = database ) (interptype = linear ) Output masks Identify database Interpolation type (x1 = INDEF) Output starting x coordinate (x2 = INDEF) Output ending x coordinate (dx = INDEF) (nx = INDEF) (xlog = no) Output X pixel interval Number of output x pixels Logarithmic x coordinate? (y1 = INDEF) Output starting y coordinate (y2 = INDEF) Output ending y coordinate (dy = INDEF) (ny = INDEF) (ylog = no) Output Y pixel interval Number of output y pixels Logarithmic y coordinate? ( flux = yes) Conserve flux per pixel? (blank = INDEF) Value for out of range pixels ( logfiles = STDOUT, logfile ) List of log files transform tr >! sed s/\(. \)/ tr\1/ gcmods11 0390. lst > trmods11 0390. lst > transform @gcmods11 0390. lst @trmods11 0390. lst gcflskycrmcsa00057147 MODS11 0390 transform ds9 4.4 Y transform FITS ds9 information panel > mdpdisplay gcflabcrmcsa00057147 MODS11 0390. f i t s frame=1 > mdpdisplay trgcflabcrmcsa00057147 MODS11 0390. f i t s frame=2

46 4 4.4.2 A-B A-B transform Y Y background bg input = output = (axis = 2) Input images to be background subtracted Output background subtracted images Axis along which background is fit and subtracted ( interactive = yes) Set fitting parameters interactively? (sample = ) (naverage = 1) (function = chebyshev ) (order = 3) Sample of points to use in fit Number of points in sample averaging Fitting function Order of fitting function ( low reject = 3.0) Low rejection in sigma of fit ( high reject = 3.0) High rejection in sigma of fit ( niterate = 3) Number of rejection iterations (grow = 0.0) (graphics = stdgraph ) (cursor = ) Rejection growing radius Graphics output device Graphics cursor input >! sed s/\(. \)/bg\1/ trmods11 0390. lst > bgmods11 0390. lst > background @trmods11 0390. l s t @bgmods11 0390. l s t Fit column X ds9 1 background identify 2 (order=3) s

4.5. 47 1 q Fit column ds9 4.5 4 A-B A B A B 8 B A > hselect @bgmods11 0390. lst $I,K DITWID yes bgtrgcflabcrmcsa00057147 MODS11 0390. fits 3.000 bgtrgcflabcrmcsa00057151 MODS11 0390. fits 3.000 bgtrgcflabcrmcsa00057159 MODS11 0390. fits 3.000 bgtrgcflabcrmcsa00057163 MODS11 0390. fits 3.000 3 MOIRCS 0.117 arcsec/pixel 3.0/0.117 26pixel imshift >! sed s/\(. \)/sh\1/ bgmods11 0390. lst > shmods11 0390. lst > imshift @bgmods11 0390. lst @shmods11 0390. lst 0 26 imarith >! sed s/\(. \)/ng\1/ shmods11 0390. lst > ngmods11 0390. lst > imarith @shmods11 0390. lst 1 @ngmods11 0390. lst imcombine median imcombine input = output = (headers = ) (bpmasks = ) (rejmasks = ) (nrejmasks = ) List of images to combine List of output images List of header files (optional) List of bad pixel masks (optional) List of rejection masks (optional) List of number rejected masks

48 4 (optional) (expmasks = ) List of exposure masks (optional) (sigmas = ) List of sigma images (optional) (imcmb = $I ) Keyword for IMCMB keywords ( logfile = STDOUT ) Log file (combine = median ) Type of combine operation ( reject = sigclip ) Type of rejection ( project = no) Project highest dimension of input images? (outtype = real ) Output image pixel datatype ( outlimits = ) Output limits (x1 x2 y1 y2...) ( offsets = none ) Input image offsets (masktype = none ) Mask type (maskvalue = 0 ) Mask value (blank = 0.0) Value if there are no pixels ( scale = exposure ) Image scaling (zero = none ) Image zero point offset (weight = exposure ) Image weights ( statsec = ) Image section for computing statistics (expname = EXPTIME ) Image header exposure time keyword ( lthreshold = INDEF) Lower threshold (hthreshold = INDEF) Upper threshold (nlow = 1) minmax: Number of low pixels to reject (nhigh = 1) minmax: Number of high pixels to reject (nkeep = 1) Minimum to keep (pos) or maximum to reject (neg) (mclip = yes) Use median in sigma clipping algorithms? (lsigma = 3.0) Lower sigma clipping factor (hsigma = 3.0) Upper sigma clipping factor (rdnoise = 0. ) ccdclip : CCD readout noise ( electrons ) (gain = 1. ) ccdclip : CCD gain ( electrons/dn) ( snoise = 0. ) ccdclip : Sensitivity noise ( fraction ) ( sigscale = 0.1) Tolerance for sigma clipping scaling corrections ( pclip = 0.5) pclip : Percentile clipping parameter (grow = 0.0) Radius ( pixels ) for neighbor rejection (mode = al )

4.6. 49 > imcombine @bgmods11 0390. l s t,@ngmods11 0390. l s t HK500 MODS11 0390 4.6 1 1 ADU (F λ [ergsec 1 cm 2 Å 1 ]) N obs (λ) F λ,int R(λ) N obs (λ) = R(λ) F λ,int (4.1) A Castelli (2004) [1] 8 2 2 HK500 M53735.fits 2 apall 2 1 2 yes > apall HK500 M53735. fits Recenter apertures for HK500 M53735? ( yes ): Resize apertures for HK500 M53735? ( yes ): Edit apertures for HK500 M53735? ( yes ): 8

50 4 4.5: apall ( : : : : ) 4.5 b 4.5 s q Trace apertures for HK500 M53735? Fit traced positions for HK500 M53735 interactively? Fit curve to aperture 1 of HK500 M53735 interactively 4.5 s :order Write apertures for HK500 M53735 to database? Extract aperture spectra for HK500 M53735? Review extracted spectra from HK500 M53735? 4.5 HK500 M53735.ms.fits N(λ)

4.6. 51 R(λ) MCSMDP rescurve stddata = one D standard star data stdmag = Vega magnitude of the standard at magfilter output = output response curve ( magfilter = mdpdb$filter/2mass K. fits ) transmission curve of filter for stdmag (moddata = mdpdb$standard/a0v g50. spc ) standard star model spectrum (vegamodel = mdpdb$filter/alpha lyr stis 003. fits ) Vega model spectrum (A0V g50.spc) K Vega Vega Vega MDPDB HIPPARCOS A0 M53735 K 2MASS 8.856 > rescurve HK500 M53735.ms. fits 8.856 resc CDFN2 HK500. fits resc CDFN2 HK500.fits R(λ) 9 2 R(λ) 2 MCSMDP mdpfcalib 1 (= > mdpfcalib HK500 MODS11 0390. fits resc CDFN2 HK500. fits HK500 MODS11 0390fl. fits HK500 MODS11-0390fl.fits [ergsec 1 cm 2 Å 1 ] 9 FITS ds9

5 53

54 5 5.1 z 2 (Yoshikawa et al. 2010 [5]) 2 2 ds9 A-B 2 Y > mdpdisplay HK500 MODS11 0390fl. fits frame=1 splot 2 1 splot > splot HK500 MODS11 0390fl. fits splot 5.1 :nsum ds9 2 5 7 ( w ) a k k Gaussian FWHM HK500 1.3µm 2.5µm z 2 6000Å 5.2 AGN 1 2 1

5.2. 55 5.1: splot a w k s q a Gaussian k fitcoords? :nsum n n 5.2: [Å] [SII] Hα 6716, 6731 1:1 6563 [NII] 6548, 6583 Hα 1:2 [OIII] 4959, 5007 1:3 Hβ 4861 5.2 splot Hα (F) (L) L = F 4πd 2 L(z) (5.1) d L (z) MCSMDP 2 1 > cosmology 2.0 Luminosity Distance(m): 4.800e+26 Angular Diameter Distance(m): 5.333e+25 Look Back Time(yr ): 1.024e+10

56 5 Hα 1 NED Extinction Calculator 2 E(B V) = 0.011 2 http://nedwww.ipac.caltech.edu/forms/calculator.html

57 [1] Castelli, F., Kurucz, R. L. 2004, arxiv:astro-ph/0405087 [2] McLean, 1997, Electronic Imaging in Astronomy; Detectors and Instrumentation (John Wiley & Sons Ltd., 1997) [3] Rousselot et al. 2000, A&A, 354, 1134 [4] Suzuki, R., et al. 2008, PASJ, 60, 1347 [5] Yoshikawa, T., et al. 2010, ApJ, 718, 112

A MCSMDP 59

60 A MCSMDP A.1 MCSMDP MCSMDP 3.1 mdpproc, mdpcombine MCSMDP S/N A-B MDP : OH 1 OH - : OH 1 OH multiextention FITS (MEF) 2 extention Poisson 3 extention SPECFIT STSDAS specfit [NII]-Hα HK500 Gaussian Poisson MOIRCS Poisson g xi σ i = 1 g r = 1 xi r g, (A.1) x i g e ADU 1 r response curve ADUerg 1 seccm 2 Å 1 2Σ obj σi 2 σ obj =, (A.2) t exp t exp 2 A-B

A.2. 61 A.2 A.2.1 mkdflat infile outfile bpm (MDPDB ) offfile normilize imstat > mkdflat @flat1. lst HK500 CDFN2 Domeflat1. fits mdpdb$bpm/nlbpm1 FF64r. fits A.2.2 mktemplate pyraf mktemplate MCSMDP mktemplate.conf mktemplate > copy MCSMDP$doc/mktemplate. conf. >! vi mktemplate. conf > mktemplate mktemplate. conf > log.py mktemplate.conf bash = $ ## parameters coordlist= mdpdb\$ohlist/list NS HK500 mdpfile= CDFN MASK02.mdp prefix= HK500 # calibration data for chip1 bpm1= mdpdb\$bpm/nlbpm1 FF64r. fits

62 A MCSMDP flat1= HK500 CDFN2 Domeflat1. fits gdb1= mdpdb\$geomap/mcsdistcrr1 feb07new.dbs gmp1= mcsdistcrr1 feb07new.gmp # calibration data for chip2 bpm2= mdpdb\$bpm/nlbpm2 FF64r. fits flat2= HK500 CDFN2 Domeflat2. fits gdb2= mdpdb\$geomap/mcsdistcrr2 feb07new.dbs gmp2= mcsdistcrr2 feb07new.gmp coordlist mdpfile prefix 2 prefix + + MDP + ḟits bpm[12 ] flat[12 ] A.2.3 mktemplate.sh AB A (B ) B(A) hselect @files. list $I @ DET ID = 1 & K DITCNT = 1 > files1a. list hselect @files. list $I @ DET ID = 1 & K DITCNT = 2 > files1b. list hselect @files. list $I @ DET ID = 2 & K DITCNT = 1 > files2a. list hselect @files. list $I @ DET ID = 2 & K DITCNT = 2 > files2b. list mdpproc unlearn mdpproc iraf.mdpproc.bpmask= mdpdb$bpm/nlbpm1 FF64r. fits iraf.mdpproc.gdb= mdpdb$geomap/mcsdistcrr1 feb07new.dbs iraf.mdpproc.gmp= mcsdistcrr1 feb07new.gmp iraf.mdpproc. flat= HK500 CDFN2 Domeflat1. fits mdpproc @obj1a. lst @obj1b. lst Multi Extention FITS (MEF) 0 ab 1 ( )

A.2. 63 A.2.4 mdptran MDP FITS > mdptran CDFN MASK02.mdp CDFN MASK02tr1.mdp gcflabcrmcsa00057147. fits 1 MDP 2 MDP 3 FITS maskplot ds9 FITS select shape then return or q to quit : return click the object position : ds9 HK500 16692Å 3 4 5 q +return MDP ds9 type q to write mdp database to file : q +return MDP q MDP cutspec mdpproc > cutspec @gcobj1a. lst CDFN MASK02tr1.mdp name=ale03 205 A.2.5 mdpcombine mdptran ds9 4

64 A MCSMDP > unlearn mdpcombine > iraf.mdpcombine. coordlist= mdpdb$ohlist/list NS HK500 > iraf.mdpcombine.ymin=42 > iraf.mdpcombine.ymax=2003 > mdpcombine @gcflab1. list CDFN MASK02tr1.mdp HK500 MODS11 0390 MODS11 0390 id mode=manual 1 id mode=manual OH mktemplate sed MDP mdpcombine!awk {print $8} CDFN MASK02tr1.mdp sed s/\(. \)/mdpcombine\ \@gcflab1\. list \ CDFN MASK02 tr1 \.mdp\ HK500 \1 \1/ mdpcombine > iraf.mdpcombine.id mode= auto > iraf.mdpcombine. refname= MODS11 0390 > iraf.mdpcombine. bg inter=no > mdpcombine @gcflab1. list CDFN MASK02tr1.mdp HK500 MODS11 0094 MODS11 0094 A.2.6 MOS mdpproc mdpcombine MDP ymin, ymax > unlearn mdpproc > iraf.mdpproc.bpmask= mdpdb$bpm/nlbpm2 FF64r. fits > iraf.mdpproc.gdb= mdpdb$geomap/mcsdistcrr2 feb07new.dbs > iraf.mdpproc.gmp= mcsdistcrr2 feb07new.gmp > iraf.mdpproc. flat= HK500 CDFN2 Domeflat2. fits > mdpproc MCSA00057114. fits MCSA00057116. fits > unlearn mdpcombine > iraf.mdpcombine. coordlist= mdpdb$ohlist/list NS HK500 > mdpcombine gcflabmcsa00057114. f i t s INDEF HK500 M53735 M53735 id mode=manual coordlist= mdpdb$ohlist/list NS HK500 calframe= INDEF ymin=900 ymax=1030

A.2. 65 idenfity id section line # background s apall 1 rescurve mdpfcalib 2 apall MEF imcopy MEF mdpfcalib > imcopy HK500 M53735. fits [0] HK500 M53735 0. fits > apall HK500 M53735 0. fits > rescurve HK500 M53735 0.ms. fits 8.856 resc CDFN2 HK500. fits > mdpfcalib HK500 MODS11 0390. fits resc CDFN2 HK500. fits HK500 MODS11 0390fl. fits A.2.7 1 matplotlib mcsmdp ipython $ mcsmdp ipython mdp2dplot mdpextract 2 mdp2dplot 2 (center) (width) nsum mdpextract SPECFIT fittype Halpha ([NII]- Hα) OIII ([OIII]-Hβ) emission SPECFIT SPECFIT mdpextract mdp2dplot In [1]: mcsmdp In [2]: unlearn mdp2dplot In [3]: iraf.mdp2dplot. xlabel=r $\lambda {\rm obs} \left [\AA\right ]$ In [4]: iraf.mdp2dplot. ylabel= r $F \lambda\,[\times10ˆ{ 18}\,erg\,cmˆ{ 2}sˆ{ 1}\AAˆ{ 1}]$ In [5]: iraf.mdp2dplot. bottomplot=false

66 A MCSMDP In [6]: iraf.mdp2dplot. yscale=1e18 In [7]: iraf.mdp2dplot. linelist= mdpdb$linelist/nebulae tex. gaia In [8]: iraf.mdp2dplot. linewidth=1 In [9]: iraf.mdp2dplot. linealpha=1 In [10]: iraf.mdp2dplot. lstyle= k In [11]: iraf.mdp2dplot.ymin=0 In [12]: mdp2dplot HK500 MODS11 0390fl 59 21990 width=800 nsum=5 In [13]: mdpextract HK500 MODS11 0390fl 59 21990 HK500 MODS11 0390fl Halpha width=800 nsum=5 sub frac=0.5 sub fwhm=42 In [14]: mdp2dplot HK500 MODS11 0390fl Halpha 59 21990 plotfit+ ymin=0 width =800 In [15]: pylab. savefig ( HK500 MODS11 0390 fl Halpha.png ) A.1

A.2. 67 MODS11-0390 [21590.0A:22390.0A][57:61] 2.0 rest frame wavelength [A] 6444 6474 6503 6533 6563 6593 6623 6653 1 ] A 1.5 1.0 10 18 1 [N II ] 6548 H [N II ] 6583 2 s erg cm [ F 0.5 0.0 21600 21700 21800 21900 22000 22100 22200 22300 [ obs A ] A.1: mdp2dplot 1