Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://xmm.vilspa.esa.es/sas/7.0.0/documentation/threads/MOS_blank_fields.html
Äàòà èçìåíåíèÿ: Tue Aug 22 17:50:41 2006
Äàòà èíäåêñèðîâàíèÿ: Sat Dec 22 14:50:31 2007
Êîäèðîâêà:
Ïîèñêîâûå ñëîâà: arp 220

How to extract MOS background spectra from blank fields event lists

MOS blank field background event lists

EPIC background spectra are normally extracted from regions of the observation field-of-view, once "contamination" due to individual discrete celestial sources has been removed. This technique is fully appropriate when analyzing spectra of point-like sources in full window modes. However, there are cases where the determination of the background from the same observation is difficult, either due to the intrinsic extended nature of a celestial object a user is interested to analyze, or because a bright source was observed in a windows mode. In the latter case, it is sometimes impossible to identify field-of-view regions, where the target source does not substantially contribute to the total X-ray counts.

In order to cope with this situation, blank sky event lists can be downloaded from the XMM-Newton Science Operations Center web site. This thread briefly explain how to use them.

This thread does not deal with the rather complex problem of a careful determination of the background normalization appropriate for a given observation. Users are referred to the EPIC background analysis web pages for more information.

set up your SAS environment (following the SAS start-up thread)

extract an image of the observation in detector coordinates

evselect table=MOSclean.fits imagebinning=binSize imageset=MOSimage.fits withimageset=yes \
  xcolumn=DETX ycolumn=DETY ximagebinsize=80 yimagebinsize=80

display the image

imgdisplay withimagefile=true imagefile=MOSimage.fits

identify the center coordinates and radius of the source spectrum extraction region. Let's assume that they are:
- DETX(center) = 100
- DETY(center) = -100
- Radius = 8 x 80 = 640

extract the source spectrum from the identified region:

evselect table=MOSclean.fits withspectrumset=yes spectrumset=MOSsource_spectrum.fits \
  energycolumn=PI spectralbinsize=15 withspecranges=yes specchannelmin=0 specchannelmax=11999 \
  expression='#XMMEA_EM && (PATTERN<=12) && ((DETX,DETY) IN circle(100,-100,640))'

download the appropriate blank field event list from the blank sky event lists repository (see here for an explanation of the naming convention). Let it be: MOSbackground.ds.

extract the background spectrum from the same region in detector coordinates as the source, using the same event selection expression:

evselect table=MOSbackground.ds withspectrumset=yes spectrumset=MOSbackground_spectrum.fits \
  energycolumn=PI spectralbinsize=15 withspecranges=yes specchannelmin=0 specchannelmax=11999 \
  expression='#XMMEA_EM && (PATTERN<=12) && ((DETX,DETY) IN circle(100,100,640))'

back to the MOS spectrum extraction thread, step n.13

It is in principle possible as well to extract background spectra from blank field event lists, if the spectrum accumulation regions are defined in sky coordinates. In this case, one needs to recast the blank field events to mimic the nominal pointing direction of the observer's field, using, for instance, the SAS task attcalc:

attcalc eventset=MOSbackground.ds attitudelabel=fixed withatthkset=N refpointlabl=user nominalra=... nominaldec=...