CIAO的简单使用，对观测数据的时变分析.

---

前言

为了记录正在学习使用钱德拉望远镜的天文软件包CIAO的一些简单的使用，写下这篇文章.

提示：以下是本篇文章正文内容，下面案例可供参考

---

一、质心修正

tar -zxvf package_12899_201115215157.tar.gz
id="953"
source activate ciao-4.12
chandra_repro ${id} outdir=""
cd ${id}/primary/
pwd
gunzip *.gz
cd ../repro/
pwd
dmkeypar acisf00${id}_repro_evt2.fits TSTART echo+
dmkeypar acisf00${id}_repro_evt2.fits TSTOP echo+
dmlist acisf00${id}_repro_evt2.fits header | egrep 'RA_TARG' > RA.log 2>&1
dmlist acisf00${id}_repro_evt2.fits header | egrep 'DEC_TARG' > DEC.log 2>&1 
ra=`cat RA.log | grep -P -o '[\d]{1,3}\.\d{5}|\-[\d]{1,3}\.\d{5}'`
dec=`cat DEC.log | grep -P -o '[\d]{1,3}\.\d{5}|\-[\d]{1,3}\.\d{5}'`
punlearn axbary
pset axbary infile=acisf00${id}_repro_evt2.fits
pset axbary orbitfile=../primary/orbitf069509100N001_eph1.fits
pset axbary outfile=acis_${id}_bary_evt2.fits 
pset axbary ra=${ra} dec=${dec}
axbary
punlearn axbary      
pset axbary infile=pcadf069583361N004_asol1.fits
pset axbary orbitfile=../primary/orbitf069509100N001_eph1.fits
pset axbary outfile=pcadf_${id}_bary_evt2.fits 
pset axbary ra=${ra} dec=${dec}
axbary mode=h
dmhedit acis_${id}_bary_evt2.fits file= op=add key=ASOLFILE value="pcad_${id}_bary_asol1.fits"
dmlist acis_${id}_bary_evt2.fits"[cols time]" data rows=1:6 
dmlist acisf00953_repro_evt2.fits"[cols time]" data rows=1:6 

---

二、时变分析

#创建两个源和一个背景
ds9 acis_${id}_bary_evt2.fits -cmap b -zoom 2 -scale log
# 确定正在使用的芯片
punlearn dmstat
dmstat "acis_${id}_bary_evt2.fits[sky=region(src1.reg)][cols ccd_id]" > ccd_id.log 2>1&
ccd="3"
dmstat "acis_${id}_bary_evt2.fits[sky=region(src2.reg)][cols ccd_id]"
dmstat "acis_${id}_bary_evt2.fits[sky=region(bkg.reg)][cols ccd_id]"
#                创建扣除背景光照曲线
#          src1
punlearn dmextract
pset dmextract infile="acis_${id}_bary_evt2.fits[ccd_id=${ccd},sky=region(src1.reg)][bin time=::2000]"
pset dmextract outfile="src1_sub_lc.fits"
pset dmextract bkg="acis_${id}_bary_evt2.fits[ccd_id=${ccd},sky=region(bkg.reg)]"
pset dmextract opt="ltc1"
pset dmextract clobber=yes
pset dmextarct verbose=0
dmextract   
#          src2
punlearn dmextract        
pset dmextract infile="acis_${id}_bary_evt2.fits[ccd_id=${ccd},sky=region(src2.reg)][bin time=::2000]"
pset dmextract outfile="src2_sub_lc.fits"
pset dmextract bkg="acis_${id}_bary_evt2.fits[ccd_id=${ccd},sky=region(bkg.reg)]"
pset dmextract opt="ltc1"                                                
pset dmextract clobber=yes
pset dmextract verbose=0                                   
dmextract                 
#          bkg
punlearn dmextract        
pset dmextract infile="acis_${id}_bary_evt2.fits[ccd_id=${ccd},sky=region(bkg.reg)][bin time=::2000]"
pset dmextract outfile="bkg_sub_lc.fits"
pset dmextract opt="ltc1"                                                   
pset dmextract clobber=yes
pset dmextract verbose=0                                                    
dmextract                 

---

1.使用Python画第一个源的光变曲线

                光变曲线可以使用matplotlib绘制    
#              src1_lc
python <<< "
from pycrates import read_file
import matplotlib.pylab as plt 
tab = read_file('src1_sub_lc.fits')
dt = tab.get_column('dt').values
rate = tab.get_column('count_rate').values
erate = tab.get_column('count_rate_err').values 
plt.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('src1_sub_lc.fits')
plt.savefig('src1_lc.pdf')
#plt.show()
"

结果：

---

2.使用Python画第二个源的光变曲线

#              src2_lc
python <<< "             
from pycrates import read_file
import matplotlib.pylab as plt 
tab = read_file('src2_sub_lc.fits')
dt = tab.get_column('dt').values
rate = tab.get_column('count_rate').values
erate = tab.get_column('count_rate_err').values 
plt.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.fits')
plt.savefig('src2-1_lc.pdf')   
#plt.show()               
"

结果：

---

3.使用Python画背景的光变曲线

#              bkg

python <<< "             
from pycrates import read_file
import matplotlib.pylab as plt 
tab = read_file('bkg_sub_lc.fits')
dt = tab.get_column('dt').values
rate = tab.get_column('count_rate').values
erate = tab.get_column('count_rate_err').values 
plt.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('bkg_sub_lc.fits')
plt.savefig('bkg_lc.pdf')   
#plt.show()               
"

结果：

---

总结

一、质心修正.
二、时变分析.
1.使用Python画第一个源的光变曲线.
2.使用Python画第二个源的光变曲线.
3.使用Python画背景的光变曲线.