hetg_ss433

As a first step we are trying to plot an older HETG observation of SS 433, viz. obsid 1019.

We download the PHA2 file along with associated ARFs and RMFs from TGCat. Here are the contents of the tarball downloaded from tgcat:

xliu@ast006:~/ss433_old$ tar -zxvf tgcat.tar.gz 
tgcat/
tgcat/obs_1019_tgid_3617/
tgcat/obs_1019_tgid_3617/heg_-1.rmf.gz
tgcat/obs_1019_tgid_3617/meg_1.arf.gz
tgcat/obs_1019_tgid_3617/heg_1.arf.gz
tgcat/obs_1019_tgid_3617/heg_1.rmf.gz
tgcat/obs_1019_tgid_3617/meg_-1.rmf.gz
tgcat/obs_1019_tgid_3617/meg_1.rmf.gz
tgcat/obs_1019_tgid_3617/pha2.gz
tgcat/obs_1019_tgid_3617/heg_-1.arf.gz
tgcat/obs_1019_tgid_3617/vv_report.txt
tgcat/obs_1019_tgid_3617/meg_-1.arf.gz
tgcat/checksums.dat

Then we fire up isis and load the PHA file

isis> load_data("pha2.gz");
Reading: ............
Integer_Type[12]
isis> list_data;
Current Spectrum List:
 id    instrument part/m  src    use/nbins   A   R     totcts   exp(ksec)  target
  1     HETG-ACIS  heg-3   1    8192/ 8192   -   -  2.2400e+02    23.352  SS 433
file:  pha2.gz
  2     HETG-ACIS  heg-2   1    8192/ 8192   -   -  5.6600e+02    23.352  SS 433
file:  pha2.gz
  3     HETG-ACIS  heg-1   1    8192/ 8192   -   -  7.3020e+03    23.352  SS 433
file:  pha2.gz
  4     HETG-ACIS  heg+1   1    8192/ 8192   -   -  6.8670e+03    23.352  SS 433
file:  pha2.gz
  5     HETG-ACIS  heg+2   1    8192/ 8192   -   -  4.3500e+02    23.352  SS 433
file:  pha2.gz
  6     HETG-ACIS  heg+3   1    8192/ 8192   -   -  1.8700e+02    23.352  SS 433
file:  pha2.gz
  7     HETG-ACIS  meg-3   1    8192/ 8192   -   -  1.0730e+03    23.352  SS 433
file:  pha2.gz
  8     HETG-ACIS  meg-2   1    8192/ 8192   -   -  2.6900e+02    23.352  SS 433
file:  pha2.gz
  9     HETG-ACIS  meg-1   1    8192/ 8192   -   -  1.0598e+04    23.352  SS 433
file:  pha2.gz
 10     HETG-ACIS  meg+1   1    8192/ 8192   -   -  1.2199e+04    23.352  SS 433
file:  pha2.gz
 11     HETG-ACIS  meg+2   1    8192/ 8192   -   -  2.2900e+02    23.352  SS 433
file:  pha2.gz
 12     HETG-ACIS  meg+3   1    8192/ 8192   -   -  8.5200e+02    23.352  SS 433
file:  pha2.gz
isis>

Next we load the response matrices (ARFs and RMFs)

heg_m1_arf = load_arf("heg_-1.arf.gz");
heg_p1_arf = load_arf("heg_1.arf.gz");
meg_p1_arf = load_arf("meg_1.arf.gz");
meg_m1_arf = load_arf("meg_-1.arf.gz");

heg_m1_rmf = load_rmf("heg_-1.rmf.gz");
heg_p1_rmf = load_rmf("heg_1.rmf.gz");
meg_p1_rmf = load_rmf("meg_1.rmf.gz");
meg_m1_rmf = load_rmf("meg_-1.rmf.gz");

Then we assign the ARFs and RMFs to the appropriate data sets, and compute the flux correction [Are we doing this correctly?]

assign_arf(heg_m1_arf, 3);  % Dataset  3 = HEG-1
assign_arf(heg_p1_arf, 4);  % Dataset  4 = HEG+1
assign_arf(meg_m1_arf, 9);  % Dataset  9 = MEG-1
assign_arf(meg_p1_arf,10);  % Dataset 10 = MEG+1

assign_rmf(heg_m1_rmf, 3);  % Dataset  3 = HEG-1
assign_rmf(heg_p1_rmf, 4);  % Dataset  4 = HEG+1
assign_rmf(meg_m1_rmf, 9);  % Dataset  9 = MEG-1
assign_rmf(meg_p1_rmf,10);  % Dataset 10 = MEG+1

flux_corr([3,4,9,10]);

Next step is to group the data and ignore bad wavelength ranges. But we're not sure which of the various grouping/rebinning commands is best to use in this case. Also, what energy range to 'xnotice' for the HEG and the MEG?

Unless otherwise stated, the content of this page is licensed under Creative Commons Attribution-ShareAlike 3.0 License