2010-07-27

製作 window file

在進行 X-ray binary 的 long-term modulation 分析的時候,常常需要先把 eclipse 的部份由光變曲線中先扣除,再進行 binning 的動作。當然,扣除 eclipse 的部份可以先把 raw data 先用90分鐘做一次 binning 輸出 light curve,再自己寫程式扣除之後做 rebin。但是這樣實在是太累了,而 heasoft 其實就有提供一個可以做 window file 的程式,叫做 xronwin。
##ReadMore##

進入 xronwin 之後,所有該有的工具都會顯示在終端機裡面:
  • [T] Change TIME Windows
  • [P] Change PHASE Windows
  • [F] Change FLUX Windows
  • [E] Change EXPOSURE Windows
  • [R] READ an Input file
  • [S] SHOW All Defined Windows
  • [C] CLEAR All Defined Windows
  • [W] WRITE Window file
  • [Q] QUIT Window Program
  • [DS] Define Series Number (Currently 1)
  • [DB] Define Bin Type (Currently Orig. Bins)
  • [DE] Define Epoch for Phase Windows (Currently Undefined)
  • [DP] Define Period for Phase Windows (Currently Undefined)
  • [DW] Copy Defined Windows from Series to Series
不過並不是所有的指令都用得到。如果要製作一個 remove eclipse 的 window file,首先要用的是 DE 和 DP 這兩個指令。DE 可以輸入 eclipse ephemeris 的 phase zero epoch,而 DP 則是定義 orbital period。DE 的內定單位是 TJD,DP 的單位則是秒。

接著,就可以按 P,進入產生 phase window 的畫面:
  • [A] Add more windows
  • [#] delete window number: #
  • to main menu
此時按 A 就可以定義在哪些 phase 範圍內的資料才是你要的。我通常習慣把 eclipse 中點附近放在 phase 0 的地方,這樣只要做一個 window 就夠了,例如 phase 0.1-0.9 這樣 (視 eclipse 時間長度而定)。

最後,回到主選單之後,再按 W,就可以輸出 window file 了。這個檔案在 lcurve 和 efold 之中都可以用到。不過在做 lcurve 之前,我都習慣先用 efold 來看看扣得乾不乾淨。


沒有留言: