前言：
為了記錄正在學習使用錢德拉望遠鏡的天文軟件包ciao，寫下這篇博客.
如有錯誤之處及建議，歡迎到評論區(qū)討論.
heasoft的安裝請參看：https://blog.csdn.net/qq_46753404/article/details/117304662?utm_source=app&app_version=4.7.1&code=app_1562916241&uLinkId=usr1mkqgl919blen

---

一.下載樣本數(shù)據(jù)

下載有專屬的命令，本文以下載觀測號13858為例.當然，你也可以官網(wǎng)先下載樣本數(shù)據(jù)，之后解壓就行.

1.下載命令：

download_chandra_obsid 13858

注意：一定要初始化ciao，否則只會顯示以下內(nèi)容：

download_chandra_obsid: command not found

下載完以后會在你的當前目錄下生成一個名為13858的目錄.

2.下載數(shù)據(jù)后，首先進行預(yù)處理.這是每次處理數(shù)據(jù)前一定要做的事情！！！以保證軟件包使用的是望遠鏡最新的相應(yīng)文件.

預(yù)處理命令：

chandra_repro 13858 outdir=""

效果顯示：

Running chandra_repro version: 03 April 2020Processing input directory '/home/wang/Downloads/13858'Resetting afterglow status bits in evt1.fits file...Running the destreak tool on the evt1.fits file...Running acis_build_badpix and acis_find_afterglow to create a new bad pixel file...Running acis_process_events to reprocess the evt1.fits file... Output from acis_process_events: # acis_process_events (CIAO 4.12): The following error occurred 2 times:dsAPEPULSEHEIGHTERR -- WARNING: pulse height is less than split threshold when performing serial CTI adjustment. Filtering the evt1.fits file by grade and status and time... Applying the good time intervals from the flt1.fits file... The new evt2.fits file is: /home/wang/Downloads/13858/repro/acisf13858_repro_evt2.fitsUpdating the event file header with chandra_repro HISTORY record Creating FOV file... Setting observation-specific bad pixel file in local ardlib.par.Cleaning up intermediate filesWARNING: Observation-specific bad pixel file set for session use:/home/wang/Downloads/13858/repro/acisf13858_repro_bpix1.fitsRun 'punlearn ardlib' when analysis of this dataset completed.The data have been reprocessed. Start your analysis with the new products in /home/wang/Downloads/13858/repro

之后就會在你之前下載的觀測號文件中生成一個名為repro的目錄

---

二.建立源區(qū)域和背景區(qū)域

可以先查看repro目錄下都有哪些文件.
1.切換目錄到repro：

cd 13858/repro

2.使用ls命令查看.

ls

效果顯示如下：

acisf13858_000N001_bpix1.fits acisf13858_repro_bpix1.fits acisf13858_000N001_fov1.fits acisf13858_repro_evt2.fits acisf13858_000N001_msk1.fits acisf13858_repro_flt2.fits acisf13858_000N001_mtl1.fits acisf13858_repro_fov1.fits acisf13858_000N001_stat1.fits acisf456520720N001_pbk0.fits acisf13858_asol1.lis pcadf456520092N001_asol1.fits

我們需要定義兩個區(qū)域，一個用于源，另一個用于背景.要做到這一點，首先用ds9顯示圖像:

ds9 acisf13858_repro_evt2.fits

效果如下：
3.換顏色
點擊color， 選b

4.放大
點擊zoom，再選擇zoom 4

5.選源（src_region）.
點擊edit，后選擇region，然后在源附近長按左鍵拖動，就會出現(xiàn)圓圈.這是src圖效果
6.選背景區(qū)域（bkg_region）
遠離亮點，選一個較大的區(qū)域作為背景.步驟同第5步

三、能譜擬合.

1.仍然在ciao初始化的環(huán)境下.依次輸入以下命令.

punlearn specextract pset specextract infile="acisf13858_repro_evt2.fits[sky=region(src.reg)]" pset specextract bkgfile="acisf13858_repro_evt2.fits[sky=region(bkg.reg)]" pset specextract outroot=13858 pset specextract clobber=yes pset specextract grouptype=NUM_CTS pset specextract binspec=30 pset specextract weight=no specextract

2.進入sherpa.

sherpa

進入界面如下：
依次輸入以下內(nèi)容：

load_pha('13858_grp.pi') ignore(':0.8,8.0:') subtract() plot_data() plt.savefig('notice.jpg') set_source(xsphabs.abs1*powlaw1d.p1) abs1.nH = 0.07 freeze(abs1.nH) guess(p1) fit() plot_fit_delchi() plt.savefig('fit.jpg')

我的擬合結(jié)果如下：

---

總結(jié)

一、下載需要的數(shù)據(jù).
二、在ds9選區(qū)域.
三、能譜擬合.

總結(jié)

以上是生活随笔為你收集整理的记录钱德拉望远镜天文软件ciao的一些简单使用，本文能谱拟合基于sherpa的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。