Документ взят из кэша поисковой машины. Адрес оригинального документа : http://xmm.vilspa.esa.es/sas/7.0.0/doc/rgscombine.ps.gz Дата изменения: Thu Jun 29 00:36:22 2006 Дата индексирования: Tue Oct 2 08:29:56 2012 Кодировка: Поисковые слова: crab nebula

XMM-Newton Science Analysis System Page: 1
rgscombine
June 28, 2006
Abstract
Combines a list of spectrum les to create one combined spectrum le, and the
corresponding list of response matrices and background spectrum, to create one nal
response matrix and background spectrum le, for a paticular RGS instrument and
spectral order.
1 Instruments/Modes
Instrument Mode
RGS Spectroscopy
2 Use
pipeline processing no
interactive analysis yes
3 Description
This task adds spectra, background spectra and response matrices of the same instrument and spectral
order, to produce combined les that can be used for quantitative analysis in a spectral tting package
such as Xspec. The component spectra should have been produced by rgsproc, using the attsyle=user
parameters so that all spectra and responses are aligned in wavelength. For example, and for each
observation of the set, the call to rgsproc could be:
rgsproc withsrc=yes srclabel=Capella srcstyle=radec srcra=79.172330 srcdec=+45.997991 attstyle=user
attra=79.172330 attdec=+45.997991 attapos=259.81
Right ascension and declination are the object coordinates and are used twice for the parameters srcra
and srcdec, and also attra and attdec. The position angle is a feature of an individual obersvation and can
be found in the POS details for each observation in XMM-Newton archive, XSA, or in the Observation
Log Browser on the XMM-Newton website.
The task diers from rgs uxer, in that it is able to combine response matrices and produce a single
response matrix for spectral analysis. The spectra attributes and response values, bar counts, are weighted
xmmsas 20060628 1801-7.0.0

XMM-Newton Science Analysis System Page: 2
prior to addition, in terms of the exposure time and areascal of the individual observation. In this way
the spectra (NET or TOTAL), including background spectra, can be combined.
4 Parameters
This section documents the parameters recognized by this task (if any).
Parameter Mand Type Default Constraints
pha yes list of
les
A list of one or more spectrum les, all produced by rgsspectrum. The task does not enforce any
constraint on the combining of dissimilar spectra.
rmf no string
A list of one or more response matrix les, all produced by rgsrmfgen.
bkg no string
A list of one or more background spectrum les, if spectra to be combined are of type TOTAL.
lepha no le
The le name for the output combined spectrum. The default name varies according to the format:
combinedspectra.ds.
lebkg no le
The le name for the output combined background spectrum if applicable. The default name varies
according to the format: combinedbkgspec.ds
lermf no le
The le name for the output combined response matrix. The default name varies according to the format:
combinedrespmat.ds
rmfgrid no int 4000
The number of energy bins in the response matrix. The default is 4000, as is the default from rgsrmfgen.
All matrices in the observation set must have the same number of bins to be combined.
min no real 4
The minimum wavelength of output for creating the grid over which to combine spectra.
xmmsas 20060628 1801-7.0.0

XMM-Newton Science Analysis System Page: 3
max no real 39
The maximum wavelength of output for creating the grid over which to combine spectra.
bins no int 3400
The number of pha channels in the spectra les.
5 Errors
This section documents warnings and errors generated by this task (if any). Note that warnings and
errors can also be generated in the SAS infrastructure libraries, in which case they would not be docu-
mented here. Refer to the index of all errors and warnings available in the HTML version of the SAS
documentation.
badResponseGrid (error)
The binning of the nal response matrix with the parameter rmfgrid, is dierent from that
found as the number for rows, in one of the input individual response matrices
pairing (error)
The number of response matrix les does not match the number of spectrum les; there
must be a one-to-one pairing.
wrongOrder (error)
Either: The RFLORDER attribute (re ection order) of the response matrix and background
spectra (if applicable) must match the RFLORDER of the corresponding spectrum.
wrongInstrument (error)
Either: The RGS instrument of the background spectra le does not match that of the pha
spectra le, if applicable and adding background spectra, or the order of the response matrix
does not match that of the spectra.
badPhaName (error)
The name of a response matrix le can be inferred from the name of the corresponding
spectrum le only if the nal 14 characters of the spectrum le name begin with \SRSPEC".
illegalBkgName (error)
The rst background spectra le name cannot be \."; it must either be \@" or else be given
explicitly.
illegalSpectraHDUCLAS (error)
A HDUCLAS attribute (OGIP), HDUCLAS2 attribute (TOTAL or NET spectrum) or HDUCLAS3
attribute (COUNT) is not consistent with a previous le in the pha spectra le list.
Alignment of spectra (error)
The pha les are not aligned in terms of wavelength. The task rgsproc will have to be
called for each pha le with the appropriate attitude parameters (attastyle, attra, attdec,
attapos).
xmmsas 20060628 1801-7.0.0

XMM-Newton Science Analysis System Page: 4
Filename parameter (warning)
No lename given for results les for either spectrum or response matrix
corrective action: No action required. Defaults used.
6 Input Files
 A list of spectrum les as produced by rgsspectrum or rgsproc.
 A list of response matrix les, as produced by rgsrmfgen or rgsproc.
 A list of background spectrum les as produced by rgsspectrum or rgsproc for TOTAL
spectra.
7 Output Files
The output spectrum le is a FITS table with the same format as the original spectrum les.
FITS A FITS-compliant DAL dataset with the table placed in an extension named \SPECTRUM". The
ve columns for a TOTAL spectra le are:
CHANNEL real32 bin center
COUNTS real32 total counts for each channel
QUALITY real32 quality ag as assessed
AREASCAL int16 weighted areascal
BACKSCAL int16 weighted backscal
8 Algorithm
FOR EACH spectrum and response file:
totalexposure = totalexposure + EXPOSURE
FOR EACH spectrum channel:
IF ( acceptable QUALITY )
channelexp = EXPOSURE * AREASCAL
backscalexp = EXPOSURE * BACKSCAL
totalcounts = totalcounts + COUNTS
Then:
Accumulate counts, backscalexp and channelexp or staterrquad
ELSE
do not add counts
FOR EACH response bin (rows):
FOR EACH response channel (columns):
IF ( corresponding spectrum channel not discarded )
newresponse = newresponse + MATRIX * EXPOSURE * AREASCAL
xmmsas\_20060628\_1801-7.0.0

XMM-Newton Science Analysis System Page: 5
FOR EACH spectrum channel:
IF (all quality bad over all spectra files)
Set quality flag to 1
Finalareascal = channelexp / totalexposure
Finalbackscal = backscalexp / totalexposure
Finalresponse = newresponse / Finalareascal
9 Comments

References
xmmsas 20060628 1801-7.0.0