 |
|
|
|
|
|
Umiestnenie potrebných dátových súborov v metacentre:
/storage/brno2/home/mkocka/data/chandra/m82
Zoznam všetkých pozorovaní ("evt" súborov), uloží sa do súboru "zoznam.txt":
find /storage/brno2/home/mkocka/data/chandra/m82/ -name "*_evt2.fits.gz" > zoznam.txt
Výber pozorovaní, ktoré sa budú spracovávať:
cat zoznam.txt | grep -e "/105" > zoznam2.txt
Vytvorenie reprojektovaných dát:
reproject_obs @zoznam2.txt reproj/
cd reproj;
Vytvorenie farebnej snímky:
dmcopy infile="merged_evt.fits[energy=200:1500][bin x=3900:4500,y=3800:4400]" outfile=soft_img.fits;
dmcopy infile="merged_evt.fits[energy=1500:2500][bin x=3900:4500,y=3800:4400]" outfile=medium_img.fits;
dmcopy infile="merged_evt.fits[energy=2500:8000][bin x=3900:4500,y=3800:4400]" outfile=hard_img.fits;
Spustenie ds9 a vytvorenie regiónov:
ds9 -rgb -red soft_img.fits -green medium_img.fits -blue hard_img.fits
POZN.: Galaxia M82 nemá aktívne galaktické jadro a preto región označený ako "jadro" vôbec nemusí odpovedať galaktickému jadru.
cd .;
specextract "reproj/10542_reproj_evt.fits[sky=region(okolie.reg)]" bkgfile="reproj/10542_reproj_evt.fits[sky=region(bg.reg)]" outroot="spec/10542" grouptype=NUM_CTS binspec=10 bkg_grouptype=NUM_CTS bkg_binspec=10;
specextract "reproj/10543_reproj_evt.fits[sky=region(okolie.reg)]" bkgfile="reproj/10543_reproj_evt.fits[sky=region(bg.reg)]" outroot="spec/10543" grouptype=NUM_CTS binspec=10 bkg_grouptype=NUM_CTS bkg_binspec=10;
specextract "reproj/10544_reproj_evt.fits[sky=region(okolie.reg)]" bkgfile="reproj/10544_reproj_evt.fits[sky=region(bg.reg)]" outroot="spec/10544" grouptype=NUM_CTS binspec=10 bkg_grouptype=NUM_CTS bkg_binspec=10;
specextract "reproj/10545_reproj_evt.fits[sky=region(okolie.reg)]" bkgfile="reproj/10545_reproj_evt.fits[sky=region(bg.reg)]" outroot="spec/10545" grouptype=NUM_CTS binspec=10 bkg_grouptype=NUM_CTS bkg_binspec=10;
specextract "reproj/10542_reproj_evt.fits[sky=region(jadro.reg)]" bkgfile="reproj/10542_reproj_evt.fits[sky=region(bg.reg)]" outroot="spec/10542j" grouptype=NUM_CTS binspec=10 bkg_grouptype=NUM_CTS bkg_binspec=10;
specextract "reproj/10543_reproj_evt.fits[sky=region(jadro.reg)]" bkgfile="reproj/10543_reproj_evt.fits[sky=region(bg.reg)]" outroot="spec/10543j" grouptype=NUM_CTS binspec=10 bkg_grouptype=NUM_CTS bkg_binspec=10;
specextract "reproj/10544_reproj_evt.fits[sky=region(jadro.reg)]" bkgfile="reproj/10544_reproj_evt.fits[sky=region(bg.reg)]" outroot="spec/10544j" grouptype=NUM_CTS binspec=10 bkg_grouptype=NUM_CTS bkg_binspec=10;
specextract "reproj/10545_reproj_evt.fits[sky=region(jadro.reg)]" bkgfile="reproj/10545_reproj_evt.fits[sky=region(bg.reg)]" outroot="spec/10545j" grouptype=NUM_CTS binspec=10 bkg_grouptype=NUM_CTS bkg_binspec=10;
cd spec;
combine_spectra src_spectra="10542.pi,10543.pi,10544.pi,10545.pi" outroot=okolie;
combine_spectra src_spectra="10542j.pi,10543j.pi,10544j.pi,10545j.pi" outroot=jadro;
XSPEC
XSPEC12>cpd /xw
XSPEC12>data jadro_src.pi
XSPEC12>setplot energy
XSPEC12>plot ld
XSPEC12>ignore **-0.6
XSPEC12>ignore 8.0-**
XSPEC12>model phabs(bknpow)
XSPEC12>fit
======================================================================== Model phabs<1>*bknpower<2> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 phabs nH 10^22 0.915131 +/- 1.41051E-02 2 2 bknpower PhoIndx1 0.994078 +/- 1.66577E-02 3 2 bknpower BreakE keV 5.43710 +/- 3.31739E-02 4 2 bknpower PhoIndx2 -2.39898 +/- 8.94488E-02 5 2 bknpower norm 3.61719E-04 +/- 7.82159E-06 ________________________________________________________________________ Fit statistic : Chi-Squared = 696.10 using 505 PHA bins. Test statistic : Chi-Squared = 696.10 using 505 PHA bins. Reduced chi-squared = 1.3922 for 500 degrees of freedom Null hypothesis probability = 1.369670e-08POZN.: Spektrum z oblasti "jadro" pravdepodobne neobsahuje žiadne emisné čiary.
XSPEC12>data okolie_src.pi
XSPEC12>ignore **-0.6
XSPEC12>ignore 8.0-**
Po hodine hrania sa s rôznymi modelmi a zlepšovaním fitovania:
======================================================================== Model phabs<1>(apec<2> + bknpower<3>) + lorentz<4> + lorentz<5> + lorentz<6> Source No.: 1 Active/On Model Model Component Parameter Unit Value par comp 1 1 phabs nH 10^22 1.47342 +/- 1.93732E-02 2 2 apec kT keV 0.926549 +/- 9.85555E-03 3 2 apec Abundanc 1.46203 +/- 0.445641 4 2 apec Redshift 0.0 frozen 5 2 apec norm 1.11871E-03 +/- 3.35374E-04 6 3 bknpower PhoIndx1 1.62447 +/- 0.179203 7 3 bknpower BreakE keV 5.27910 +/- 6.41401E-02 8 3 bknpower PhoIndx2 9.50000 +/- 1.05144 9 3 bknpower norm 1.41252E-04 +/- 3.86825E-05 10 4 lorentz LineE keV 0.821034 +/- 5.53119E-03 11 4 lorentz Width keV 0.150807 +/- 1.51569E-02 12 4 lorentz norm 1.95282E-05 +/- 1.34140E-06 13 5 lorentz LineE keV 6.68900 +/- 9.65514E-03 14 5 lorentz Width keV 3.67708E-02 +/- 2.68482E-02 15 5 lorentz norm 3.88768E-06 +/- 5.07861E-07 16 6 lorentz LineE keV 6.22648 +/- 2.84352E-02 17 6 lorentz Width keV 0.292283 +/- 0.107545 18 6 lorentz norm 3.39931E-06 +/- 8.83293E-07 ________________________________________________________________________ Fit statistic : Chi-Squared = 782.18 using 505 PHA bins. Test statistic : Chi-Squared = 782.18 using 505 PHA bins. Reduced chi-squared = 1.6028 for 488 degrees of freedom Null hypothesis probability = 5.389734e-16
Čiary nafitované modelom apec sa nachádzajú v súbore ciary.txt
Výber niektorých čiar s vysokou emisivitou:
Wavelength [keV] ion upper/lower level emisivita [ph cm³/s]
# 1.0215 Ne X 3 1 5.72E-017
# 1.022 Ne X 4 1 1.14E-016
# 1.211 Ne X 7 1 1.52E-017
# 1.3522 Mg XI 7 1 3.08E-017
# 1.3311 Mg XI 2 1 1.39E-017
# 1.4717 Mg XII 3 1 3.40E-017
# 1.4726 Mg XII 4 1 6.70E-017
# 1.8538 Si XIII 5 1 1.12E-017
# 1.865 Si XIII 7 1 7.61E-017
# 1.8394 Si XIII 2 1 3.64E-017
# 2.0043 Si XIV 3 1 2.32E-017
# 2.0061 Si XIV 4 1 4.57E-017
# 2.4606 S XV 7 1 3.09E-017
# 2.4303 S XV 2 1 1.47E-017
# 1.0093 Fe XXI 40 1 5.38E-016
# 1.0534 Fe XXII 21 1 2.57E-016
# 1.0564 Fe XXIII 20 5 1.74E-016
Manuálne vložené čiary s lorentzovským profilom:
Čiary 6.7 a 6.2 keV podľa tohto zdroja patria železu.
Pre určenie zastúpenia čiary Si XIII boli vybraté energie v intervale od 1.80 keV do 1.90 keV.
Galaxia M82 v takto zvolenom intervale vyzerá takto:
Pre lepšie ukázanie zastúpenia čiary kremíka v galaxii M82 je na naslednovnej snímke zobrazený celý interval 0.2-8.0 keV modrou farbou a čiara kremíka červenou farbou. Obe farby boli logaritmicky naškálované podľa intenzít. Čím je teda oblasť červenšia, tým je väčšia pravdepodobnosť zdroja emisnej kremíkovej čiary. Zaujímavosťou je modrá oblasť blízko centra, kde je buď menší výskyt takto ionizovaného kremíka alebo tam dochádza k absorpcii.
Pre porovnanie je tu ešte čiara železa 6.7 keV. Intenzita čiary je avšak veľmi nízka vzhľadom na šum v tejto oblasti energií a preto nasledovná snímka nemôže byť dôveryhodná ohľadom zastúpenia železa.
A možno je všetko úplne inak ...
specextract "reproj/1546_reproj_evt.fits[sky=region(ds9.reg)]" bkgfile="reproj/1546_reproj_evt.fits[sky=region(bg.reg)]" outroot="spec/1546" grouptype=NUM_CTS binspec=10 bkg_grouptype=NUM_CTS bkg_binspec=10;
 |
|
|
|