 |
|
|
|
|
|
V adresári s dátami sa zvyčajne nachádza adresár odf.
/home/pato/AP3-2/0111570101/odf
Príklad pre zobrazenie všetkých pozorovaní ("evt" súborov), skript v pythone:
import glob, pyfits
print("Pocet eventov:", sum([pyfits.open(f)[1].header["NAXIS2"] for f in glob.glob("/home/data/sandbox/3C273_all/*/odf/*/*ingEvts.ds")]))
Základné nastavenia v adresári odf:
heainit xmmsasinit2
export SAS_ODF=$PWD
export SAS_CCF=${SAS_ODF}/ccf.cif
Vytvorenie svetelnej krivky:
evselect table=2408_0693781301_EPN_S003_ImagingEvts.ds rateset='sel_rate3.lc' withrateset=yes timebinsize=50 maketimecolumn=yes timecolumn=TIME expression='#XMMEA_EP && (PI in [2000:12000]) && (PATTERN<=12)'
tabgtigen table=sel_rate3.lc gtiset=sel_reg.gti expression="COUNTS.gt.10.0.and.COUNTS.lt.500.0"
http://xmm.esac.esa.int/sas/current/doc/tabgtigen/node3.html
Vytvorenie snímky:
evselect table=2408_0693781301_EPN_S003_ImagingEvts.ds withimageset=true imageset=ima2_view.fits xcolumn=DETX ycolumn=DETY ximagebinsize=10 yimagebinsize=10 expression=" #XMMEA_EP && (PI in [300:12000]) && (PATTERN<=12) && GTI(sel_reg.gti,TIME)"
Vyextrahovanie spektra objektu:
evselect table=2408_0693781301_EPN_S003_ImagingEvts.ds withspectrumset=true withspecranges=true energycolumn=PI specchannelmin=0 specchannelmax=11999 spectralbinsize=15 updateexposure=yes spectrumset=agn_spec.pha expression=" #XMMEA_EP && (PI in [300:12000]) && (PATTERN<=12) && GTI(sel_reg.gti,TIME) && ((DETX,DETY) IN circle(72,-220,400))"
Vyextrahovanie spektra pozadia:
evselect table=2408_0693781301_EPN_S003_ImagingEvts.ds withspectrumset=true withspecranges=true energycolumn=PI specchannelmin=0 specchannelmax=11999 spectralbinsize=15 updateexposure=yes spectrumset=bkg_spec.pha expression=" #XMMEA_EP && (PI in [300:12000]) && (PATTERN<=12) && GTI(sel_reg.gti,TIME) && ((DETX,DETY) IN circle(2232,3230,800))"
Vygenerovanie matíc rmf a arf:
rmfgen spectrumset=agn_spec.pha rmfset=spec.rmf extendedsource=no badpixlocation=evtlist.ds detmaptype=psf
arfgen spectrumset=agn_spec.pha arfset=spec.arf withrmfset=yes rmfset=spec.rmf extendedsource=no badpixlocation=evtlist.ds detmaptype=psf?
Spojenie spektier a matíc do jedného súboru:
specgroup spectrumset=agn_spec.pha mincounts=25 oversample=3 rmfset=spec.rmf arfset=spec.arf backgndset=bkg_spec.pha groupedset=PN_spectrum_AGN.fits
 |
|
|
|