連続波ポインティング測定のデータ解析

1. ファイル群のコピー

アンテナログファイルをaste-mx2にコピー

  1. アンテナログファイルを aste-1c → aste-2c → aste-at3i → aste-mx2 の順番でコピー。(2011年度版)
  % ssh aste10@aste-mx2
  % ssh aste10@aste-2c

解析に使うファイルリストを作成する

  1. 指示書と観測データファイル(*.LCK)のリストを作成する。
    • 下記の作業ディレクトリで指示書名と観測データファイル名(*.LCK)のリストを作成する。
      • ディレクトリ: /home/astebeam/pointing/observations
      • リスト名: C01015jup.list
 asteobs@aste-2c% grep "Apr[ ]* 21" /log/6/mess | grep START | grep -v Cosmos | grep -v SEND > tmp.txt
 asteobs@aste-at3i% sftp asteobs@aste-2c
 sftp>get tmp.txt
 asteobs@aste-mt2%sftp asteobs@aste-at3i
 sftp>get tmp.txt
 astebeam@aste-mx2 % sftp asteobs@aste-mt2
 sftp> get tmp.txt
 astebeam@aste-mx2 % cat tmp.txt | grep sbrf11 | awk '{print $9,$8}' | sed -e "s/aste11.sbrf11.//;s/'[ ]*11/.start C1/;s/00.000/.LCK/" > tmp2.txt
 astebeam@aste-mx2 % cp tmp2.txt ~/pointing/observations
  1. アンテナログファイルのリストを作成する。
  1. tmp2.txtとtmp3.txtを整形して「指示書名」「lock-inファイル名」「アンテナログファイル名」「dAZ, dEL」を列挙したリスト C00519sat.list を作成。
    • 列を並べるにはpasteコマンドが便利。
    • 【重要】リストの書き方
      • 同時に行っている光ポインティングの指示書はすっ飛ばしてください。
      • LCK と real の時刻が対応していることを確認する。
      • LCKファイルの拡張子 (.LCK) は不要です。
      • 測定の最初に行った jpx5cc3.start を用いたポインティングについては、記載しないでください。
      • 列の順序は下記の通りにして下さい。列の区切りはスペース・タブ、両方可能。
         # 指示書名      LCKファイル名 アンテナログファイル名  dAZ  dEL   コメント
         ijpx5c3.start C010150719 real20101015071759 -7.0 +3.0 
         ijpx5c3.start C010150811 real20101015081026 -7.0 +3.0 
         ijpx5c3.start C010150843 real20101015084232 -7.0 +18.  # EL<10

ファイルリスト作成時のtips

2. iOPTテーブルの転送

  1. iOPT の結果が書き込まれたテーブルの名前を調べる。asteobs@aste-mx2 → asteobs@aste-2c → astepon@aste-c1c (pw: kollama!) の順番でログイン。
astepon@aste-c1c % cd ~/optpon/
astepon@aste-c1c % ls 20101015*
201010150123.tab                                  # ← ここでは、このファイルを例にとります。
  1. aste-2cにログイン
     astebeam@aste-mx2 % ssh aste-2c -l asteobs
     asteobs@aste-2c % cd subref/iopttab
  2. 測定時に起動していたiOPTが吐いた巨大なテーブル201006020123.tabに含まれる必要な情報(ここでは172-24Dという恒星)をreduce_iopttab で抜き出して、201010150123_red.tabという名前で保存。
    # USAGE: reduce_iopttab (iOPTが出力したテーブル) (恒星の名前)
    asteobs@aste-2c % ssh astepon@aste-c1c '~/optpon/reduce_iopttab ~/optpon/201010150123.tab 172-24D' > 201010150123_red.tab
  3. 三鷹に転送 (aste-2c → aste-mx2)
    astebeam@aste-mx2 % scp asteobs@aste-2c:~/subref/iopttab/201010150123_red.tab ~/pointing/observations

3. 連続波クロススキャンデータの一次解析

  1. ログイン。(astebeam@aste-mx2)
  2. 一次解析。(update on 07-Oct-10)
     % source ~/.cshrc.linux
     % cd ~/pointing/observations/
     % point_redux -l C01015jup.list

C01015jup.listに記述されている測定時間ごとの .dazdel ファイル (e.g., C010150447.dazdel)と、それらをまとめたファイル(C01015jup.dazdel)が作られる。

4. 電波軸と光学望遠鏡軸のオフセットの測定

  1. オフセットを測定。
     astebeam@aste-mx2 % cd ~/pointing/observations
     astebeam@aste-mx2 % iopt_track C01015jup.list 201010150123_red.tab
    iopt_trackは、主鏡成分を引き去った副鏡指向誤差成分をファイル C01015jup.r2o.dazdel に書き出す(r2oはradio-to-opticalの意味)。
  2. 画像 (PNG, PS) が生成されているので、これをwikiに貼付ける。
     astebeam@aste-mx2 % ls *.png
     C01015jup.r2o.png

#ref(): File not found: "iopt_track.png" at page "Observation and Analysis/SubRef/DataAnalysis_Pointing"

5. ポイティングのAz, EL依存性

(※) PS画像を書き出すには、~/pointing/observations/ ディレクトリに必要ファイルをコピーするなどして作業する必要がある。

6. 器差ファイルの更新と再測定による確認

以上。


Front page   New Page list Search Recent changes   Help   RSS of recent changes