12/15/2005 [[[ 3 ]]]
Slit 17 (ch2)
SCUBA 850.1
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 17 zJ totak= 6600s = 110min
-----------------------------------
day1:
imcopy ../data/MCSA00014576.fits[*,1265:1395] 17o_zJ_14576.fits zJ500 300s
imcopy ../data/MCSA00014578.fits[*,1265:1395] 17o_zJ_14578.fits
imcopy ../data/MCSA00014580.fits[*,1265:1395] 17o_zJ_14580.fits 600s
imcopy ../data/MCSA00014582.fits[*,1265:1395] 17o_zJ_14582.fits
imcopy ../data/MCSA00014584.fits[*,1265:1395] 17o_zJ_14584.fits
imcopy ../data/MCSA00014586.fits[*,1265:1395] 17o_zJ_14586.fits
imcopy ../data/MCSA00014588.fits[*,1265:1395] 17o_zJ_14588.fits 存在しない
imcopy ../data/MCSA00014590.fits[*,1265:1395] 17o_zJ_14590.fits
imcopy ../data/MCSA00014592.fits[*,1265:1395] 17o_zJ_14592.fits
day2:
imcopy ../data/MCSA00014966.fits[*,1250:1380] 17o_zJ_14966.fits zJ500 600sec
imcopy ../data/MCSA00014968.fits[*,1250:1380] 17o_zJ_14968.fits
imcopy ../data/MCSA00014970.fits[*,1250:1380] 17o_zJ_14970.fits
imcopy ../data/MCSA00014972.fits[*,1250:1380] 17o_zJ_14972.fits
---------
差し引き
---------
day1:
imarith 17o_zJ_14576.fits - 17o_zJ_14578.fits 17o_zJ_76-78.fits
imarith 17o_zJ_14580.fits - 17o_zJ_14582.fits 17o_zJ_80-82.fits
imarith 17o_zJ_14584.fits - 17o_zJ_14586.fits 17o_zJ_84-86.fits
imarith 17o_zJ_14590.fits - 17o_zJ_14592.fits 17o_zJ_90-92.fits
day2:
imarith 17o_zJ_14966.fits - 17o_zJ_14968.fits 17o_zJ_66-68.fits
imarith 17o_zJ_14970.fits - 17o_zJ_14972.fits 17o_zJ_70-72.fits
---------
波長較正
---------
for day1:
identify 17o_zJ_14580.fits section="middle line" coordli="OHlines_zJ500.dat" nsum=5 function="cheby"
17+l order=3 rms=0.9801
reidentify 17o_zJ_14580.fits 17o_zJ_14580.fits section="middle line" step=10 nsum=5 overrid+ interac+ nlost=3
fitcoord 17o_zJ_14580 function="cheby" xorder=1 yorder=1 interac+
rms=1.908
X order = 3
Y order = 1
(1, 1) = 5069.082 (2048, 1) = 16773.22
(1, 131) = 5079.231 (2048, 131) = 16783.37
transform 17o_zJ_76-78.fits 17o_zJ_76-78_tr.fits 17o_zJ_14580 interpt="linear"
transform 17o_zJ_80-82.fits 17o_zJ_80-82_tr.fits 17o_zJ_14580 interpt="linear"
transform 17o_zJ_84-86.fits 17o_zJ_84-86_tr.fits 17o_zJ_14580 interpt="linear"
transform 17o_zJ_90-92.fits 17o_zJ_90-92_tr.fits 17o_zJ_14580 interpt="linear"
for day2:
identify 17o_zJ_14966.fits section="middle line" coordli="OHlines_zJ500.dat" nsum=5 function="cheby"
14+l order=3 0.7537
reidentify 17o_zJ_14966.fits 17o_zJ_14966.fits section="middle line" step=10 nsum=5 overrid+ interac+ nlost=3
fitcoord 17o_zJ_14966 function="cheby" xorder=1 yorder=1 interac+
2.207
X order = 3
Y order = 1
(1, 1) = 4879.997 (2048, 1) = 16765.23
(1, 131) = 4889.625 (2048, 131) = 16774.86
transform 17o_zJ_66-68.fits 17o_zJ_66-68_tr.fits 17o_zJ_14966 interpt="linear"
transform 17o_zJ_70-72.fits 17o_zJ_70-72_tr.fits 17o_zJ_14966 interpt="linear"
------
FLAT p_25mm
------
DISPAXISが入ってなかった。objのときは入ってたのに。ヘッダはダークになっている。
スリットごとにCALフラットがその位置で撮れているか確認する。
for day1:
imcopy ../data/MCSA00014836.fits[*,1265:1395] 17f_zJ_14836_1.fits
imcopy ../data/MCSA00014838.fits[*,1265:1395] 17f_zJ_14838_1.fits
imcopy ../data/MCSA00014840.fits[*,1265:1395] 17f_zJ_14840_1.fits
imcombine 17f_zJ_14836_1.fits,17f_zJ_14838_1.fits,17f_zJ_14840_1.fits 17f_zj_1_p25.fits combine="median" reject="none"
hedit 17f_zj_1_p25.fits DISPAXIS 1 add+ ver-
transform 17f_zj_1_p25.fits 17f_zj_1_p25_tr.fits 17o_zJ_14580 interpt="linear"
blkavg 17f_zj_1_p25_tr.fits 17f_zj_1_p25_trblk.fits 1 131 option="average"
imcopy 17f_zj_1_p25_trblk.fits[*,1] 17f_zj_1_p25_trblk1.fits
imarith 17f_zj_1_p25_tr.fits / 17f_zj_1_p25_trblk1.fits 17f_zj_1_p25_trn.fits
ccdmask 17f_zj_1_p25_trn.fits 17f_zj_1_p25_trn
imcopy 17f_zj_1_p25_trn.fits 17f_zj_1_p25_trnm.fits
fixpix 17f_zj_1_p25_trnm.fits 17f_zj_1_p25_trn.pl << flat for day1
for day2:
imcopy ../data/MCSA00014836.fits[*,1250:1380] 17f_zJ_14836_2.fits
imcopy ../data/MCSA00014838.fits[*,1250:1380] 17f_zJ_14838_2.fits
imcopy ../data/MCSA00014840.fits[*,1250:1380] 17f_zJ_14840_2.fits
imcombine 17f_zJ_14836_2.fits,17f_zJ_14838_2.fits,17f_zJ_14840_2.fits 17f_zj_2_p25.fits combine="median" reject="none"
hedit 17f_zj_2_p25.fits DISPAXIS 1 add+ ver-
transform 17f_zj_2_p25.fits 17f_zj_2_p25_tr.fits 17o_zJ_14966 interpt="linear"
blkavg 17f_zj_2_p25_tr.fits 17f_zj_2_p25_trblk.fits 1 131 option="average"
imcopy 17f_zj_2_p25_trblk.fits[*,1] 17f_zj_2_p25_trblk1.fits
imarith 17f_zj_2_p25_tr.fits / 17f_zj_2_p25_trblk1.fits 17f_zj_2_p25_trn.fits
ccdmask 17f_zj_2_p25_trn.fits 17f_zj_2_p25_trn
imcopy 17f_zj_2_p25_trn.fits 17f_zj_2_p25_trnm.fits
fixpix 17f_zj_2_p25_trnm.fits 17f_zj_2_p25_trn.pl << flat for day2
--------------
MASK for OBJ
--------------
day1:
ccdmask 17o_zJ_76-78_tr.fits 17o_zJ_76-78_tr
ccdmask 17o_zJ_80-82_tr.fits 17o_zJ_80-82_tr
ccdmask 17o_zJ_84-86_tr.fits 17o_zJ_84-86_tr
ccdmask 17o_zJ_90-92_tr.fits 17o_zJ_90-92_tr
imcopy 17o_zJ_76-78_tr.fits 17o_zJ_76-78_trm.fits
imcopy 17o_zJ_80-82_tr.fits 17o_zJ_80-82_trm.fits
imcopy 17o_zJ_84-86_tr.fits 17o_zJ_84-86_trm.fits
imcopy 17o_zJ_90-92_tr.fits 17o_zJ_90-92_trm.fits
fixpix 17o_zJ_76-78_trm.fits 17o_zJ_76-78_tr.pl
fixpix 17o_zJ_80-82_trm.fits 17o_zJ_80-82_tr.pl
fixpix 17o_zJ_84-86_trm.fits 17o_zJ_84-86_tr.pl
fixpix 17o_zJ_90-92_trm.fits 17o_zJ_90-92_tr.pl
day2:
ccdmask 17o_zJ_66-68_tr.fits 17o_zJ_66-68_tr
ccdmask 17o_zJ_70-72_tr.fits 17o_zJ_70-72_tr
imcopy 17o_zJ_66-68_tr.fits 17o_zJ_66-68_trm.fits
imcopy 17o_zJ_70-72_tr.fits 17o_zJ_70-72_trm.fits
fixpix 17o_zJ_66-68_trm.fits 17o_zJ_66-68_tr.pl
fixpix 17o_zJ_70-72_trm.fits 17o_zJ_70-72_tr.pl
memo:
ここまでの200512のSXDS解析、(slit-12,6,17)は以下の設定のままやっている。デフォルト値のまま。
ccdmask ncmed=nlmed=7
ncsig=nlsig=15
lsigma=hsigma=6
ngood=5
ちょっとむらむらした感じになる。
------------
OBJ / FLAT
------------
for day1:
imarith 17o_zJ_76-78_trm.fits / 17f_zj_1_p25_trnm.fits 17o_zJ_76-78_trmf.fits
imarith 17o_zJ_80-82_trm.fits / 17f_zj_1_p25_trnm.fits 17o_zJ_80-82_trmf.fits
imarith 17o_zJ_84-86_trm.fits / 17f_zj_1_p25_trnm.fits 17o_zJ_84-86_trmf.fits
imarith 17o_zJ_90-92_trm.fits / 17f_zj_1_p25_trnm.fits 17o_zJ_90-92_trmf.fits
for day2:
imarith 17o_zJ_66-68_trm.fits / 17f_zj_2_p25_trnm.fits 17o_zJ_66-68_trmf.fits
imarith 17o_zJ_70-72_trm.fits / 17f_zj_2_p25_trnm.fits 17o_zJ_70-72_trmf.fits
-----------
重ね合わせ
-----------
day1:
imarith 17o_zJ_76-78_trmf.fits * -1.0 17o_zJ_76-78_trmf-1.fits
imshift 17o_zJ_76-78_trmf-1.fits 17o_zJ_76-78_trmf-1sh.fits 0 44
imarith 17o_zJ_76-78_trmf.fits + 17o_zJ_76-78_trmf-1sh.fits 17o_zJ_76-78_trmf_comb.fits
imarith 17o_zJ_80-82_trmf.fits * -1.0 17o_zJ_80-82_trmf-1.fits
imshift 17o_zJ_80-82_trmf-1.fits 17o_zJ_80-82_trmf-1sh.fits 0 44
imarith 17o_zJ_80-82_trmf.fits + 17o_zJ_80-82_trmf-1sh.fits 17o_zJ_80-82_trmf_comb.fits
imarith 17o_zJ_84-86_trmf.fits * -1.0 17o_zJ_84-86_trmf-1.fits
imshift 17o_zJ_84-86_trmf-1.fits 17o_zJ_84-86_trmf-1sh.fits 0 44
imarith 17o_zJ_84-86_trmf.fits + 17o_zJ_84-86_trmf-1sh.fits 17o_zJ_84-86_trmf_comb.fits
imarith 17o_zJ_90-92_trmf.fits * -1.0 17o_zJ_90-92_trmf-1.fits
imshift 17o_zJ_90-92_trmf-1.fits 17o_zJ_90-92_trmf-1sh.fits 0 44
imarith 17o_zJ_90-92_trmf.fits + 17o_zJ_90-92_trmf-1sh.fits 17o_zJ_90-92_trmf_comb.fits
imcombine 17o_zJ_76-78_trmf_comb.fits,17o_zJ_80-82_trmf_comb.fits,17o_zJ_84-86_trmf_comb.fits,17o_zJ_90-92_trmf_comb.fits 17o_zJ_day1.fits combine="sum" reject="none"
imcopy 17o_zJ_day1.fits[*,70:140] 17o_zJ_day1s.fits
day2:
imarith 17o_zJ_66-68_trmf.fits * -1.0 17o_zJ_66-68_trmf-1.fits
imshift 17o_zJ_66-68_trmf-1.fits 17o_zJ_66-68_trmf-1sh.fits 0 44
imarith 17o_zJ_66-68_trmf.fits + 17o_zJ_66-68_trmf-1sh.fits 17o_zJ_66-68_trmf_comb.fits
imarith 17o_zJ_70-72_trmf.fits * -1.0 17o_zJ_70-72_trmf-1.fits
imshift 17o_zJ_70-72_trmf-1.fits 17o_zJ_70-72_trmf-1sh.fits 0 44
imarith 17o_zJ_70-72_trmf.fits + 17o_zJ_70-72_trmf-1sh.fits 17o_zJ_70-72_trmf_comb.fits
imcombine 17o_zJ_66-68_trmf_comb.fits,17o_zJ_70-72_trmf_comb.fits 17o_zJ_day2.fits combine="sum" reject="none"
#imcopy 17o_zJ_day2.fits[*,70:130] 17o_zJ_day2s.fits
１日目も２日目も何もみえない。
----------------------------------------------------------------------------------------
SXDS HK
imarith ../data/MCSA00014988.fits - ../data/MCSA00014990.fits 14988-14990.fits
-----------
Slit 6 HK total = 5400s = 100min
-----------
imcopy ../data/MCSA00014988.fits[1:1750,1245:1370] 17o_HK_14988.fits HK500 600sec
imcopy ../data/MCSA00014990.fits[1:1750,1245:1370] 17o_HK_14990.fits
imcopy ../data/MCSA00014996.fits[1:1750,1245:1370] 17o_HK_14996.fits
imcopy ../data/MCSA00014998.fits[1:1750,1245:1370] 17o_HK_14998.fits
imcopy ../data/MCSA00015000.fits[1:1750,1245:1370] 17o_HK_15000.fits
imcopy ../data/MCSA00015002.fits[1:1750,1245:1370] 17o_HK_15002.fits
imcopy ../data/MCSA00015004.fits[1:1750,1245:1370] 17o_HK_15004.fits
imcopy ../data/MCSA00015006.fits[1:1750,1245:1370] 17o_HK_15006.fits
imcopy ../data/MCSA00014984.fits[1:1750,1245:1370] 17o_HK_14984.fits HK500 300sec
imcopy ../data/MCSA00014986.fits[1:1750,1245:1370] 17o_HK_14986.fits
imarith 17o_HK_14988.fits - 17o_HK_14990.fits 17o_HK_88-90.fits
imarith 17o_HK_14996.fits - 17o_HK_14998.fits 17o_HK_96-98.fits
imarith 17o_HK_15000.fits - 17o_HK_15002.fits 17o_HK_00-02.fits
imarith 17o_HK_15004.fits - 17o_HK_15006.fits 17o_HK_04-06.fits
imarith 17o_HK_14984.fits - 17o_HK_14986.fits 17o_HK_84-86.fits 300s
identify 17o_HK_14988.fits section="middle line" coordli="OHlines_HK500.dat" nsum=5 function="cheby"
9+both order=3 0.9268
reidentify 17o_HK_14988.fits 17o_HK_14988.fits section="middle line" step=10 nsum=5 overrid+ interac+ nlost=3
いつものことだがフィットがきたない
fitcoord 17o_HK_14988 function="cheby" xorder=1 yorder=1 interac+
2.798
X order = 3
Y order = 3
(1, 1) = 10806.02 (1750, 1) = 24791.6
(1, 126) = 10831.81 (1750, 126) = 24791.54
transform 17o_HK_88-90.fits 17o_HK_88-90_tr.fits 17o_HK_14988 interpt="linear"
transform 17o_HK_96-98.fits 17o_HK_96-98_tr.fits 17o_HK_14988 interpt="linear"
transform 17o_HK_00-02.fits 17o_HK_00-02_tr.fits 17o_HK_14988 interpt="linear"
transform 17o_HK_04-06.fits 17o_HK_04-06_tr.fits 17o_HK_14988 interpt="linear"
transform 17o_HK_84-86.fits 17o_HK_84-86_tr.fits 17o_HK_14988 interpt="linear"
----------
FLAT p25mm
----------
imcopy ../data/MCSA00014882.fits[1:1750,1245:1370] 17f_HK_14882.fits p25mm
imcopy ../data/MCSA00014884.fits[1:1750,1245:1370] 17f_HK_14884.fits
imcopy ../data/MCSA00014886.fits[1:1750,1245:1370] 17f_HK_14886.fits
imcombine 17f_HK_14882.fits,17f_HK_14884.fits,17f_HK_14886.fits 17f_HK_p25.fits combine="median" reject="none"
hedit 17f_HK_p25.fits DISPAXIS 1 add+ ver-
transform 17f_HK_p25.fits 17f_HK_p25_tr.fits 17o_HK_14988 interpt="linear"
blkavg 17f_HK_p25_tr.fits 17f_HK_p25_trblk.fits 1 125 option="average"
imcopy 17f_HK_p25_trblk.fits[*,1] 17f_HK_p25_trblk1.fits
imarith 17f_HK_p25_tr.fits / 17f_HK_p25_trblk1.fits 17f_HK_p25_trn.fits
ccdmask 17f_HK_p25_trn.fits 17f_HK_p25_trn
imcopy 17f_HK_p25_trn.fits 17f_HK_p25_trnm.fits
fixpix 17f_HK_p25_trnm.fits 17f_HK_p25_trn.pl << flat
-------------
MASK for OBJ
-------------
ccdmask 17o_HK_88-90_tr.fits 17o_HK_88-90_tr
ccdmask 17o_HK_96-98_tr.fits 17o_HK_96-98_tr
ccdmask 17o_HK_00-02_tr.fits 17o_HK_00-02_tr
ccdmask 17o_HK_04-06_tr.fits 17o_HK_04-06_tr
ccdmask 17o_HK_84-86_tr.fits 17o_HK_84-86_tr
imcopy 17o_HK_88-90_tr.fits 17o_HK_88-90_trm.fits
imcopy 17o_HK_96-98_tr.fits 17o_HK_96-98_trm.fits
imcopy 17o_HK_00-02_tr.fits 17o_HK_00-02_trm.fits
imcopy 17o_HK_04-06_tr.fits 17o_HK_04-06_trm.fits
imcopy 17o_HK_84-86_tr.fits 17o_HK_84-86_trm.fits
fixpix 17o_HK_88-90_trm.fits 17o_HK_88-90_tr.pl
fixpix 17o_HK_96-98_trm.fits 17o_HK_96-98_tr.pl
fixpix 17o_HK_00-02_trm.fits 17o_HK_00-02_tr.pl
fixpix 17o_HK_04-06_trm.fits 17o_HK_04-06_tr.pl
fixpix 17o_HK_84-86_trm.fits 17o_HK_84-86_tr.pl
------------
OBJ / FLAT
------------
imarith 17o_HK_88-90_trm.fits / 17f_HK_p25_trnm.fits 17o_HK_88-90_trmf.fits
imarith 17o_HK_96-98_trm.fits / 17f_HK_p25_trnm.fits 17o_HK_96-98_trmf.fits
imarith 17o_HK_00-02_trm.fits / 17f_HK_p25_trnm.fits 17o_HK_00-02_trmf.fits
imarith 17o_HK_04-06_trm.fits / 17f_HK_p25_trnm.fits 17o_HK_04-06_trmf.fits
imarith 17o_HK_84-86_trm.fits / 17f_HK_p25_trnm.fits 17o_HK_84-86_trmf.fits
-----------
重ね合わせ
-----------
imarith 17o_HK_88-90_trmf.fits * -1.0 17o_HK_88-90_trmf-1.fits
imshift 17o_HK_88-90_trmf-1.fits 17o_HK_88-90_trmf-1sh.fits 0 44
imarith 17o_HK_88-90_trmf.fits + 17o_HK_88-90_trmf-1sh.fits 17o_HK_88-90_trmf_comb.fits
imarith 17o_HK_96-98_trmf.fits * -1.0 17o_HK_96-98_trmf-1.fits
imshift 17o_HK_96-98_trmf-1.fits 17o_HK_96-98_trmf-1sh.fits 0 44
imarith 17o_HK_96-98_trmf.fits + 17o_HK_96-98_trmf-1sh.fits 17o_HK_96-98_trmf_comb.fits
imarith 17o_HK_00-02_trmf.fits * -1.0 17o_HK_00-02_trmf-1.fits
imshift 17o_HK_00-02_trmf-1.fits 17o_HK_00-02_trmf-1sh.fits 0 44
imarith 17o_HK_00-02_trmf.fits + 17o_HK_00-02_trmf-1sh.fits 17o_HK_00-02_trmf_comb.fits
imarith 17o_HK_04-06_trmf.fits * -1.0 17o_HK_04-06_trmf-1.fits
imshift 17o_HK_04-06_trmf-1.fits 17o_HK_04-06_trmf-1sh.fits 0 44
imarith 17o_HK_04-06_trmf.fits + 17o_HK_04-06_trmf-1sh.fits 17o_HK_04-06_trmf_comb.fits
imarith 17o_HK_84-86_trmf.fits * -1.0 17o_HK_84-86_trmf-1.fits
imshift 17o_HK_84-86_trmf-1.fits 17o_HK_84-86_trmf-1sh.fits 0 44
imarith 17o_HK_84-86_trmf.fits + 17o_HK_84-86_trmf-1sh.fits 17o_HK_84-86_trmf_comb.fits
imcombine 17o_HK_88-90_trmf_comb.fits,17o_HK_96-98_trmf_comb.fits,17o_HK_00-02_trmf_comb.fits,17o_HK_04-06_trmf_comb.fits,17o_HK_84-86_trmf_comb.fits 17o_HK_trmf_all.fits combine="sum" reject="none"
imcopy 17o_HK_trmf_all.fits[*,65:130] 17o_HK_trmf_alls.fits