1.0 Download the data

观测号为953和461

官网下载

下载数据的网址：https://cda.harvard.edu/chaser/

在Observation ID输入观测ID953,461
点击Search
Select all
Add Products to Retrieval List
Retrieve Products

点击网址

下载压缩包
解压

ciao下载

~/chandra/test5 » ciao  #初始化ciao
(ciao-4.12) --------------------------------------------------------------------
~/chandra/test5 » download_chandra_obsid 953,461 #下载数据
Downloading files for ObsId 11823, total size is 366 Mb.Type     Format      Size  0........H.........1  Download Time Average Rate...
Total download size for ObsId 11823 = 366 Mb
Total download time for ObsId 11823 = 7 m 19 s

1.1 数据预处理

~/chandra » chandra_repro 953，461 outdir="" verbose=5   Running chandra_repro
version: 03 April 2020

进入预处理后的repro目录

1.2 质心修正

别忘了初始化ciao

ciao
punlearn axbary
pset axbary infile=acisf00953_repro_evt2.fits
pset axbary orbitfile=../primary/orbitf069509100N001_eph1.fits
pset axbary outfile=acis_953_bary_evt2.fits
pset axbary ra=6.02430 dec=-72.08090
axbary

dmhedit acisf00953_repro_evt2.fits file= op=add key=ASOLFILE value="pcad_133_bary_asol1.fits"
dmhedit acisf00953_repro_evt2.fits"[cols time]" data rows=1:6

用ds9打开acis_953_bary_evt2.fits文件

ds9 acis_953_bary_evt2.fits -cmap b -scale log &

1.3 Edit -> Region

分别做出如下三个区域，并用ds9打开

ds9 acis_953_bary_evt2.fits -cmap b -scale log -region src1.reg -region src2.reg -region bkg.reg

a: For src1

punlearn dmextract
dmextract infile="acis_953_bary_evt2.fits [ccd_id=3,sky=region(src1.reg)][bin time=::2000]" outfile="src1_sub_lc.fits" bkg="acis_953_bary_evt2.fits[ccd_id=3,sky=region(bg.reg)]" opt="ltc1" clobber=yes verbose=2

b: For src2

punlearn dmextract
dmextract infile="acis_953_bary_evt2.fits[ccd_id=3,sky=region(src2.reg)][bin time=::2000]" outfile="src2_sub_lc.fits" bkg="acis_953_bary_evt2.fits[ccd_id=3,sky=region(bg.reg)]" opt="ltc1" clobber=yes verbose=0

1.4 Plot the light curves

a: For src1 by fv in Heasoft

初始化heasoft,用fv打开fv src2_sub_lc.fits

h
fv src2_sub_lc.fits &

LIGHTCURVE>>Plot
设置坐标

然后Print,保存

b: For src1 by python

写python脚本

vi pl.py

脚本内容如下


```python
from pycrates import read_file
import matplotlib.pylab as plttab = read_file("src2_sub_lc_400.fits")
dt = tab.get_column("dt").values
rate = tab.get_column("net_rate").values
erate = tab.get_column("err_rate").valuesplt.errorbar(dt, rate, yerr=erate,marker="o",color="red", mfc="black", mec="black", ecolor="grey")
plt.xlabel("$\Delta$ T (sec)")
plt.ylabel("Net Count Rate (counts/sec)")
plt.title("src2_sub_lc_400.fits")
plt.show()

用python运行pl.py

python pl.py

如图

2.0 Surce2 lightcurve with 200s time-bin(分为200bins,上面操作为2000bins)

改为200bins

punlearn dmextract
dmextract infile="acis_953_bary_evt2.fits[ccd_id=3,sky=region(src2.reg)][bin time=::200]" outfile="src2_sub_lc_200.fits" bkg="acis_953_bary_evt2.fits[ccd_id=3,sky=region(bg.reg)]" opt="ltc1" clobber=yes verbose=0

画图

3.0 Light curve for HRC(高分辨率CDD461)

Update two files

punlearn axbary
pset axbary infile=pcadf064938693N003_asol1.fits
pset axbary orbitfile=../primary/orbitf064886464N002_eph1.fits
pset axbary outfile=pcad_461_bary_asol1.fits
pset axbary ra=187.2779170 dec=2.0523890
pset axbary clobber=yes
axbary mode=h

X-ray数据处理笔记：Chandra Timing Analysis相关推荐

python离散积分_科学网—python数据处理笔记（三）通道积分图 - 钱磊的博文
通道积分是除了p-v图外又一件分析分子云气体运动的利器.在不同的通道查看,可以看到一些速度不同寻常的气体,这些气体就是恒星形成研究中十分感兴趣的高速气体.这些高速气体对于分子云中的能量平衡以及恒星的形 ...
《Static Timing Analysis for Nanometer Designs》汇总传送门
文章目录 [第一章绪论](https://blog.csdn.net/yueguangmuyu/article/details/125721646) [第二章 STA相关概念](https://bl ...
python笔记图片_科学网—python数据处理笔记（二）p-v图 - 钱磊的博文
python数据处理笔记(二)p-v图已有 9903 次阅读 2012-5-24 17:46 |个人分类:知识|系统分类:科研笔记| Python, 分子云, 常用工具, p-v图 p-v图是分析分 ...
Operation condition And timing analysis
Operation condition 主要指芯片运行的环境的不同,在这里加上Interconnect net 模型. Operating condition Description Process ...
python数值比较器_python笔记16（数据处理笔记1）
学习线性代数的时候,想要可视化,但是仅仅是三维坐标系,我都感觉到手画图很不直观.所以这次的python学习,一方面是把之前没学完的学完,另一方面,也是想了解下图形绘制的技术.以下内容和6.00笔记10 ...
自监督论文阅读笔记 Urban feature analysis from aerial remote sensing imagery using self-supervised and semi-s
"Urban feature analysis from aerial remote sensing imagery using self-supervised and semi-super ...
python如何提取数据中的年月_Python数据处理笔记——Pandas时间数据提取处理
做数据分析时,对于有时间数据的数据来源,在时间维度上的剖析必不可少比如: 在一天的销售时间内,哪些时间段是高峰是否与星期有关 RFM模型中的R怎么快速计算距离天数如何找出某个时间点或时间段的数据 ...
读书笔记《Outlier Analysis》第八章分类、文本和混合属性中的异常检测
1.引言前面讨论的都是数值数据,然而,在现实生活生成中,还有很多其他类型的数据,如性别.种族.邮编.人员和实体的名称.IP地址等.通常处理这些数据更有挑战,因为难以以均匀和一致的方式来处理各种类型的 ...
cs224w（图机器学习）2021冬季课程学习笔记4 Link Analysis: PageRank (Graph as Matrix)
诸神缄默不语-个人CSDN博文目录 cs224w(图机器学习)2021冬季课程学习笔记集合文章目录 1. Graph as Matrix 2. PageRank / the Google Algor ...

X-ray数据处理笔记：Chandra Timing Analysis