12/14/2005   [[[ 1 ]]]

Slit 12

z=2.53
 [OII]= 13156A
[OIII]= 17675A
    Ha= 23167A

SXDS zJ
imarith ../data/MCSA00014965.fits - ../data/MCSA00014967.fits 14965-14967.fits
imarith ../data/MCSA00014966.fits - ../data/MCSA00014968.fits 14966-14968.fits
imarith ../data/MCSA00014576.fits - ../data/MCSA00014578.fits 14576-14578.fits

-----------------------------------
Slit 12  zJ  totak= 6600s = 110min
-----------------------------------
2日にまたがる
スリットの位置が異なる ~10pix
同じ領域を撮ってフラットを同一の使おうと思ったけど、
波長較正でこけるかもしれないので別々にやることにする。

day1:
imcopy ../data/MCSA00014576.fits[*,340:480] 12o_zJ_14576.fits zJ500 300s
imcopy ../data/MCSA00014578.fits[*,340:480] 12o_zJ_14578.fits
imcopy ../data/MCSA00014580.fits[*,340:480] 12o_zJ_14580.fits       600s
imcopy ../data/MCSA00014582.fits[*,340:480] 12o_zJ_14582.fits
imcopy ../data/MCSA00014584.fits[*,340:480] 12o_zJ_14584.fits
imcopy ../data/MCSA00014586.fits[*,340:480] 12o_zJ_14586.fits
imcopy ../data/MCSA00014588.fits[*,340:480] 12o_zJ_14588.fits 存在しない
imcopy ../data/MCSA00014590.fits[*,340:480] 12o_zJ_14590.fits
imcopy ../data/MCSA00014592.fits[*,340:480] 12o_zJ_14592.fits

day2:
imcopy ../data/MCSA00014966.fits[*,330:460] 12o_zJ_14966.fits  zJ500 600sec
imcopy ../data/MCSA00014968.fits[*,330:460] 12o_zJ_14968.fits
imcopy ../data/MCSA00014970.fits[*,330:460] 12o_zJ_14970.fits
imcopy ../data/MCSA00014972.fits[*,330:460] 12o_zJ_14972.fits

---------
差し引き
---------
day1:
imarith 12o_zJ_14576.fits - 12o_zJ_14578.fits 12o_zJ_76-78.fits 
imarith 12o_zJ_14580.fits - 12o_zJ_14582.fits 12o_zJ_80-82.fits 
imarith 12o_zJ_14584.fits - 12o_zJ_14586.fits 12o_zJ_84-86.fits 
imarith 12o_zJ_14590.fits - 12o_zJ_14592.fits 12o_zJ_90-92.fits 

day2:
imarith 12o_zJ_14966.fits - 12o_zJ_14968.fits 12o_zJ_66-68.fits 
imarith 12o_zJ_14970.fits - 12o_zJ_14972.fits 12o_zJ_70-72.fits

---------
波長較正
---------
for day1:
identify 12o_zJ_14580.fits section="middle line" coordli="OHlines_zJ500.dat" nsum=5 function="cheby"
14+l order=3 rms=0.8117
reidentify 12o_zJ_14580.fits 12o_zJ_14580.fits section="middle line" step=10 nsum=5 overrid+ interac+ nlost=3
fitcoord 12o_zJ_14580 function="cheby" xorder=1 yorder=1 interac+
rms=2.73
    X order = 3
    Y order = 1
    (1, 1) = 6425.517  (2048, 1) = 18224.98
    (1, 141) = 6438.404  (2048, 141) = 18237.86

transform 12o_zJ_76-78.fits 12o_zJ_76-78_tr.fits 12o_zJ_14580 interpt="linear"
transform 12o_zJ_80-82.fits 12o_zJ_80-82_tr.fits 12o_zJ_14580 interpt="linear"
transform 12o_zJ_84-86.fits 12o_zJ_84-86_tr.fits 12o_zJ_14580 interpt="linear"
transform 12o_zJ_90-92.fits 12o_zJ_90-92_tr.fits 12o_zJ_14580 interpt="linear"


for day2:
identify 12o_zJ_14966.fits section="middle line" coordli="OHlines_zJ500.dat" nsum=5 function="cheby"
17 order=3 0.9686
reidentify 12o_zJ_14966.fits 12o_zJ_14966.fits section="middle line" step=10 nsum=5 overrid+ interac+ nlost=3
fitcoord 12o_zJ_14966 function="cheby" xorder=1 yorder=1 interac+
2.538
    X order = 3
    Y order = 1
    (1, 1) = 6502.662  (2048, 1) = 18233.07
    (1, 131) = 6513.268  (2048, 131) = 18243.68

transform 12o_zJ_66-68.fits 12o_zJ_66-68_tr.fits 12o_zJ_14966 interpt="linear"
transform 12o_zJ_70-72.fits 12o_zJ_70-72_tr.fits 12o_zJ_14966 interpt="linear"

------
FLAT
------
DISPAXISが入ってなかった。objのときは入ってたのに。ヘッダはダークになっている。

for day1:
imcopy ../data/MCSA00014824.fits[*,340:480] 12f_zJ_14824_1.fits
imcopy ../data/MCSA00014826.fits[*,340:480] 12f_zJ_14826_1.fits
imcopy ../data/MCSA00014828.fits[*,340:480] 12f_zJ_14828_1.fits
imcombine 12f_zJ_14824_1.fits,12f_zJ_14826_1.fits,12f_zJ_14828_1.fits 12f_zj_1_p25.fits combine="median" reject="none"

hedit 12f_zj_1_p25.fits DISPAXIS 1 add+ ver-
transform 12f_zj_1_p25.fits 12f_zj_1_p25_tr.fits 12o_zJ_14580 interpt="linear"

blkavg 12f_zj_1_p25_tr.fits 12f_zj_1_p25_trblk.fits 1 141 option="average"
imcopy 12f_zj_1_p25_trblk.fits[*,1] 12f_zj_1_p25_trblk1.fits
imarith 12f_zj_1_p25_tr.fits / 12f_zj_1_p25_trblk1.fits 12f_zj_1_p25_trn.fits
ccdmask 12f_zj_1_p25_trn.fits 12f_zj_1_p25_trn
imcopy 12f_zj_1_p25_trn.fits 12f_zj_1_p25_trnm.fits 
fixpix 12f_zj_1_p25_trnm.fits 12f_zj_1_p25_trn.pl   << flat for day1


for day2:
imcopy ../data/MCSA00014824.fits[*,330:460] 12f_zJ_14824_2.fits
imcopy ../data/MCSA00014826.fits[*,330:460] 12f_zJ_14826_2.fits
imcopy ../data/MCSA00014828.fits[*,330:460] 12f_zJ_14828_2.fits
imcombine 12f_zJ_14824_2.fits,12f_zJ_14826_2.fits,12f_zJ_14828_2.fits 12f_zj_2_p25.fits combine="median" reject="none"

hedit 12f_zj_2_p25.fits DISPAXIS 1 add+ ver-
transform 12f_zj_2_p25.fits 12f_zj_2_p25_tr.fits 12o_zJ_14966 interpt="linear"

blkavg 12f_zj_2_p25_tr.fits 12f_zj_2_p25_trblk.fits 1 131 option="average"
imcopy 12f_zj_2_p25_trblk.fits[*,1] 12f_zj_2_p25_trblk1.fits
imarith 12f_zj_2_p25_tr.fits / 12f_zj_2_p25_trblk1.fits 12f_zj_2_p25_trn.fits
ccdmask 12f_zj_2_p25_trn.fits 12f_zj_2_p25_trn
imcopy 12f_zj_2_p25_trn.fits 12f_zj_2_p25_trnm.fits 
fixpix 12f_zj_2_p25_trnm.fits 12f_zj_2_p25_trn.pl   << flat for day2

--------------
MASK for OBJ
--------------
ccdmask 12o_zJ_76-78_tr.fits 12o_zJ_76-78_tr
ccdmask 12o_zJ_80-82_tr.fits 12o_zJ_80-82_tr
ccdmask 12o_zJ_84-86_tr.fits 12o_zJ_84-86_tr
ccdmask 12o_zJ_90-92_tr.fits 12o_zJ_90-92_tr

imcopy 12o_zJ_76-78_tr.fits 12o_zJ_76-78_trm.fits
imcopy 12o_zJ_80-82_tr.fits 12o_zJ_80-82_trm.fits
imcopy 12o_zJ_84-86_tr.fits 12o_zJ_84-86_trm.fits
imcopy 12o_zJ_90-92_tr.fits 12o_zJ_90-92_trm.fits

fixpix 12o_zJ_76-78_trm.fits 12o_zJ_76-78_tr.pl
fixpix 12o_zJ_80-82_trm.fits 12o_zJ_80-82_tr.pl
fixpix 12o_zJ_84-86_trm.fits 12o_zJ_84-86_tr.pl
fixpix 12o_zJ_90-92_trm.fits 12o_zJ_90-92_tr.pl

ccdmask 12o_zJ_66-68_tr.fits 12o_zJ_66-68_tr
imcopy 12o_zJ_66-68_tr.fits 12o_zJ_66-68_trm.fits
fixpix 12o_zJ_66-68_trm.fits 12o_zJ_66-68_tr.pl

ccdmask 12o_zJ_70-72_tr.fits 12o_zJ_70-72_tr
imcopy 12o_zJ_70-72_tr.fits 12o_zJ_70-72_trm.fits
fixpix 12o_zJ_70-72_trm.fits 12o_zJ_70-72_tr.pl

------------
OBJ / FLAT
------------
for day1:
imarith 12o_zJ_76-78_trm.fits / 12f_zj_1_p25_trnm.fits 12o_zJ_76-78_trmf.fits
imarith 12o_zJ_80-82_trm.fits / 12f_zj_1_p25_trnm.fits 12o_zJ_80-82_trmf.fits 
imarith 12o_zJ_84-86_trm.fits / 12f_zj_1_p25_trnm.fits 12o_zJ_84-86_trmf.fits
imarith 12o_zJ_90-92_trm.fits / 12f_zj_1_p25_trnm.fits 12o_zJ_90-92_trmf.fits 
for day2:
imarith 12o_zJ_66-68_trm.fits / 12f_zj_2_p25_trnm.fits 12o_zJ_66-68_trmf.fits
imarith 12o_zJ_70-72_trm.fits / 12f_zj_2_p25_trnm.fits 12o_zJ_70-72_trmf.fits


-----------
重ね合わせ
-----------
day1:
imarith 12o_zJ_76-78_trmf.fits * -1.0 12o_zJ_76-78_trmf-1.fits
imshift 12o_zJ_76-78_trmf-1.fits 12o_zJ_76-78_trmf-1sh.fits 0 44
imarith 12o_zJ_76-78_trmf.fits + 12o_zJ_76-78_trmf-1sh.fits 12o_zJ_76-78_trmf_comb.fits

imarith 12o_zJ_80-82_trmf.fits * -1.0 12o_zJ_80-82_trmf-1.fits
imshift 12o_zJ_80-82_trmf-1.fits 12o_zJ_80-82_trmf-1sh.fits 0 44
imarith 12o_zJ_80-82_trmf.fits + 12o_zJ_80-82_trmf-1sh.fits 12o_zJ_80-82_trmf_comb.fits

imarith 12o_zJ_84-86_trmf.fits * -1.0 12o_zJ_84-86_trmf-1.fits
imshift 12o_zJ_84-86_trmf-1.fits 12o_zJ_84-86_trmf-1sh.fits 0 44
imarith 12o_zJ_84-86_trmf.fits + 12o_zJ_84-86_trmf-1sh.fits 12o_zJ_84-86_trmf_comb.fits

imarith 12o_zJ_90-92_trmf.fits * -1.0 12o_zJ_90-92_trmf-1.fits
imshift 12o_zJ_90-92_trmf-1.fits 12o_zJ_90-92_trmf-1sh.fits 0 44
imarith 12o_zJ_90-92_trmf.fits + 12o_zJ_90-92_trmf-1sh.fits 12o_zJ_90-92_trmf_comb.fits

imcombine 12o_zJ_76-78_trmf_comb.fits,12o_zJ_80-82_trmf_comb.fits,12o_zJ_84-86_trmf_comb.fits,12o_zJ_90-92_trmf_comb.fits 12o_zJ_day1.fits combine="sum" reject="none"
imcopy 12o_zJ_day1.fits[*,70:140] 12o_zJ_day1s.fits

day2:
imarith 12o_zJ_66-68_trmf.fits * -1.0 12o_zJ_66-68_trmf-1.fits
imshift 12o_zJ_66-68_trmf-1.fits 12o_zJ_66-68_trmf-1sh.fits 0 44
imarith 12o_zJ_66-68_trmf.fits + 12o_zJ_66-68_trmf-1sh.fits 12o_zJ_66-68_trmf_comb.fits

imarith 12o_zJ_70-72_trmf.fits * -1.0 12o_zJ_70-72_trmf-1.fits
imshift 12o_zJ_70-72_trmf-1.fits 12o_zJ_70-72_trmf-1sh.fits 0 44
imarith 12o_zJ_70-72_trmf.fits + 12o_zJ_70-72_trmf-1sh.fits 12o_zJ_70-72_trmf_comb.fits

imcombine 12o_zJ_66-68_trmf_comb.fits,12o_zJ_70-72_trmf_comb.fits 12o_zJ_day2.fits combine="sum" reject="none"
imcopy 12o_zJ_day2.fits[*,70:130] 12o_zJ_day2s.fits

見えてる。14900Aあたりに輝線?でもz=2.53なので、rest-frame=4220A?



----------------------------------------------------------------------------------------
SXDS HK
imarith ../data/MCSA00014988.fits - ../data/MCSA00014990.fits 14988-14990.fits


----------- 
Slit 12  HK        total = 6000s = 100min
-----------
imcopy ../data/MCSA00014988.fits[*,320:450] 12o_HK_14988.fits  HK500 600sec
imcopy ../data/MCSA00014990.fits[*,320:450] 12o_HK_14990.fits
imcopy ../data/MCSA00014992.fits[*,320:450] 12o_HK_14992.fits
#imcopy ../data/MCSA00014994.fits[*,320:450] 12o_HK_14994.fits 94データなし
imcopy ../data/MCSA00014996.fits[*,320:450] 12o_HK_14996.fits
imcopy ../data/MCSA00014998.fits[*,320:450] 12o_HK_14998.fits
imcopy ../data/MCSA00015000.fits[*,320:450] 12o_HK_15000.fits
imcopy ../data/MCSA00015002.fits[*,320:450] 12o_HK_15002.fits
imcopy ../data/MCSA00015004.fits[*,320:450] 12o_HK_15004.fits
imcopy ../data/MCSA00015006.fits[*,320:450] 12o_HK_15006.fits
imcopy ../data/MCSA00014984.fits[*,320:450] 12o_HK_14984.fits  HK500 300sec
imcopy ../data/MCSA00014986.fits[*,320:450] 12o_HK_14986.fits  

imarith 12o_HK_14988.fits - 12o_HK_14990.fits 12o_HK_88-90.fits
imarith 12o_HK_14996.fits - 12o_HK_14998.fits 12o_HK_96-98.fits
imarith 12o_HK_15000.fits - 12o_HK_15002.fits 12o_HK_00-02.fits
imarith 12o_HK_15004.fits - 12o_HK_15006.fits 12o_HK_04-06.fits
imarith 12o_HK_14984.fits - 12o_HK_14986.fits 12o_HK_84-86.fits 300s
imarith 12o_HK_14992.fits - 12o_HK_14990.fits 12o_HK_92-90.fits 600s 上のスペクトルだけ使う



identify 12o_HK_14988.fits section="middle line" coordli="OHlines_HK500.dat" nsum=5 function="cheby"
13+r order=4 2.864
reidentify 12o_HK_14988.fits 12o_HK_14988.fits section="middle line" step=10 nsum=5 overrid+ interac+ nlost=3
fitcoord 12o_HK_14988 function="cheby" xorder=1 yorder=1 interac+
4.988
    X order = 4
    Y order = 2
    (1, 1) = 13188.79  (2048, 1) = 28669.61
    (1, 131) = 13164.85  (2048, 131) = 28721.69

transform 12o_HK_88-90.fits 12o_HK_88-90_tr.fits 12o_HK_14988 interpt="linear"
transform 12o_HK_96-98.fits 12o_HK_96-98_tr.fits 12o_HK_14988 interpt="linear"
transform 12o_HK_00-02.fits 12o_HK_00-02_tr.fits 12o_HK_14988 interpt="linear"
transform 12o_HK_04-06.fits 12o_HK_04-06_tr.fits 12o_HK_14988 interpt="linear"
transform 12o_HK_84-86.fits 12o_HK_84-86_tr.fits 12o_HK_14988 interpt="linear"
transform 12o_HK_92-90.fits 12o_HK_92-90_tr.fits 12o_HK_14988 interpt="linear"

------
FLAT
------
imcopy ../data/MCSA00014870.fits[*,320:450] 12f_HK_14870.fits
imcopy ../data/MCSA00014872.fits[*,320:450] 12f_HK_14872.fits
imcopy ../data/MCSA00014874.fits[*,320:450] 12f_HK_14874.fits
imcombine 12f_HK_14870.fits,12f_HK_14872.fits,12f_HK_14874.fits 12f_HK_p-25.fits combine="median" reject="none"

hedit 12f_HK_p-25.fits DISPAXIS 1 add+ ver-
transform 12f_HK_p-25.fits 12f_HK_p-25_tr.fits 12o_HK_14988 interpt="linear"

blkavg 12f_HK_p-25_tr.fits 12f_HK_p-25_trblk.fits 1 131 option="average"
imcopy 12f_HK_p-25_trblk.fits[*,1] 12f_HK_p-25_trblk1.fits
imarith 12f_HK_p-25_tr.fits / 12f_HK_p-25_trblk1.fits 12f_HK_p-25_trn.fits
ccdmask 12f_HK_p-25_trn.fits 12f_HK_p-25_trn
imcopy 12f_HK_p-25_trn.fits 12f_HK_p-25_trnm.fits 
fixpix 12f_HK_p-25_trnm.fits 12f_HK_p-25_trn.pl   << flat

ccdmask 12o_HK_88-90_tr.fits 12o_HK_88-90_tr
imcopy 12o_HK_88-90_tr.fits 12o_HK_88-90_trm.fits
fixpix 12o_HK_88-90_trm.fits 12o_HK_88-90_tr.pl

ccdmask 12o_HK_96-98_tr.fits 12o_HK_96-98_tr
imcopy 12o_HK_96-98_tr.fits 12o_HK_96-98_trm.fits
fixpix 12o_HK_96-98_trm.fits 12o_HK_96-98_tr.pl

ccdmask 12o_HK_00-02_tr.fits 12o_HK_00-02_tr
imcopy 12o_HK_00-02_tr.fits 12o_HK_00-02_trm.fits
fixpix 12o_HK_00-02_trm.fits 12o_HK_00-02_tr.pl

ccdmask 12o_HK_04-06_tr.fits 12o_HK_04-06_tr
imcopy 12o_HK_04-06_tr.fits 12o_HK_04-06_trm.fits
fixpix 12o_HK_04-06_trm.fits 12o_HK_04-06_tr.pl

ccdmask 12o_HK_84-86_tr.fits 12o_HK_84-86_tr
imcopy 12o_HK_84-86_tr.fits 12o_HK_84-86_trm.fits
fixpix 12o_HK_84-86_trm.fits 12o_HK_84-86_tr.pl

ccdmask 12o_HK_92-90_tr.fits 12o_HK_92-90_tr
imcopy 12o_HK_92-90_tr.fits 12o_HK_92-90_trm.fits
fixpix 12o_HK_92-90_trm.fits 12o_HK_92-90_tr.pl



------------
OBJ / FLAT
------------
imarith 12o_HK_88-90_trm.fits / 12f_HK_p-25_trnm.fits 12o_HK_88-90_trmf.fits
imarith 12o_HK_96-98_trm.fits / 12f_HK_p-25_trnm.fits 12o_HK_96-98_trmf.fits
imarith 12o_HK_00-02_trm.fits / 12f_HK_p-25_trnm.fits 12o_HK_00-02_trmf.fits
imarith 12o_HK_04-06_trm.fits / 12f_HK_p-25_trnm.fits 12o_HK_04-06_trmf.fits
imarith 12o_HK_84-86_trm.fits / 12f_HK_p-25_trnm.fits 12o_HK_84-86_trmf.fits
imarith 12o_HK_92-90_trm.fits / 12f_HK_p-25_trnm.fits 12o_HK_92-90_trmf.fits

imarith 12o_HK_88-90_trmf.fits * -1.0 12o_HK_88-90_trmf-1.fits
imshift 12o_HK_88-90_trmf-1.fits 12o_HK_88-90_trmf-1sh.fits 0 44
imarith 12o_HK_88-90_trmf.fits + 12o_HK_88-90_trmf-1sh.fits 12o_HK_88-90_trmf_comb.fits

imarith 12o_HK_96-98_trmf.fits * -1.0 12o_HK_96-98_trmf-1.fits
imshift 12o_HK_96-98_trmf-1.fits 12o_HK_96-98_trmf-1sh.fits 0 44
imarith 12o_HK_96-98_trmf.fits + 12o_HK_96-98_trmf-1sh.fits 12o_HK_96-98_trmf_comb.fits

imarith 12o_HK_00-02_trmf.fits * -1.0 12o_HK_00-02_trmf-1.fits
imshift 12o_HK_00-02_trmf-1.fits 12o_HK_00-02_trmf-1sh.fits 0 44
imarith 12o_HK_00-02_trmf.fits + 12o_HK_00-02_trmf-1sh.fits 12o_HK_00-02_trmf_comb.fits

imarith 12o_HK_04-06_trmf.fits * -1.0 12o_HK_04-06_trmf-1.fits
imshift 12o_HK_04-06_trmf-1.fits 12o_HK_04-06_trmf-1sh.fits 0 44
imarith 12o_HK_04-06_trmf.fits + 12o_HK_04-06_trmf-1sh.fits 12o_HK_04-06_trmf_comb.fits

imarith 12o_HK_84-86_trmf.fits * -1.0 12o_HK_84-86_trmf-1.fits
imshift 12o_HK_84-86_trmf-1.fits 12o_HK_84-86_trmf-1sh.fits 0 44
imarith 12o_HK_84-86_trmf.fits + 12o_HK_84-86_trmf-1sh.fits 12o_HK_84-86_trmf_comb.fits

12o_HK_92-90_trmf.fits はこのまま。

imcombine 12o_HK_88-90_trmf_comb.fits,12o_HK_96-98_trmf_comb.fits,12o_HK_00-02_trmf_comb.fits,12o_HK_04-06_trmf_comb.fits,12o_HK_84-86_trmf_comb.fits,12o_HK_92-90_trmf.fits 12o_HK_trmf_all.fits combine="sum" reject="none"
imcopy 12o_HK_trmf_all.fits[1:1400,70:130] 12o_HK_trmf_alls.fits

14900のところにはやっぱり輝線がみえる。
z=2.5がまちがい?
14900がHaだとすると、z=1.27になる。