Skip banner and navigation tools.
The Swift-XRT data products generator
Build XRT products.
Skip contents.
Contents
Overview
This facility allows the creation of X-ray light curves, spectra and positions of any object in the Swift XRT
field of view. Images of the field of view can also be created.
Details of how the light curves, spectra and enhanced positions are produced are given in Evans
et al. (2009), (MNRAS, 397, 1177).
The tools which create these products are based on those which run
automatically for GRBs (which in turn
are based on the Swift
software and FTOOLs), with the following modifications:
- Light curves can be binned by fixing the number of counts per
bin and allowing the bin duration to vary, (as for GRBs) or
with fixed-duration bins, or one bin per snapshot or
obsID
- If no UVOT observations exist in the white, v or b filters, the UV filters can be used to determine an enhanced position (note:
the systematic error is larger in this circumstance).
- Up to four time-resolved spectra can be created, instead of a single time-averaged spectrum.
- Unenhanced positions are also produced (using a PSF fit rather).
- For spectra there is only a single absorber in the automatic fit, and this is not fixed at the Galactic value.
- The observations used to create each product can be specified.
- Images in user-defined energy bands can be created.
Although spectra are automatically fitted, the spectral files are
available for download as a tar file should users wish to fit
different models or interact with the fit. (The RMF file is not
included in this file, but is part of the CALDB.)
To create products, simply fill out the form. The object details
section is mandatory and the information provided here will be used to identify the source
for which products should be built. You also use this section to specify which products you
require. For each of the products you choose, an extra form will appear
asking for details of the desired product. The label for each form field is a link; following it opens a help
box in the top left of the window, which gives summary information about the field. Full details
for each option follow here. Note that any control characters (& ;
etc) will be removed from your input.
Back to contents | Back to product generator
To build products you need to supply basic details of the object and
observations, and choose which products you want. The following
sections describe each of the forms in detail. Once you submit the
form you will be taken to a page telling you the status of your job.
The products will appear on this page when they are built, and will
be kept for 7 days. Note that it is safe to leave the page or close
the browser while the products are building and the job will not be
terminated. However we recommend that you bookmark the page first so
you can find your products later. Alternatively, if you supply your
e-mail address on the main form you will be notified by e-mail when
each of your products is built. That e-mail will also contain a link
to the page in case you did not bookmark it. Note that every time
you submit the products form a new job is executed, so re-entering
the product details in the form will start the product building
again from scratch, it will not take you to the status page of the
existing job.
- Name
- The name of the object. Primarily a cosmetic detail, however the "Find"
button will use the name to determine the target ID, RA, Dec and start time if possible. The
RA and Dec will be taken from SIMBAD, if the name can be resolved. If not, the values will be
taken from the spacecraft pointing information during the first observation of the source.
This may differ from the object position, so the default search radius is increased to 3' in
this case. If you change the name but do not update the co-ordinates, a warning will be
issued. Note -- when you submit the form, the co-ordinates supplied are not checked against
any catalogued position for the source name. Although your products will be labelled with the
name you supplied, they will be produced for the position you supplied.
- Target ID
- The target ID of the object. If your object
has multiple target IDs which you wish to include in the light
curve, these should be separated by commas, with no spaces, e.g.
00050100,00058990. This field is filled in
automatically if the name is resolved successfully.
- Start time
- The zero-point of the observations. This is
used to scale the x-axis in the plots. If this is left blank,
the TSTART value of the first event list will be used. Accepted time formats are given
below. If the name is successfully resolved this field will be automatically set to
contain the Swift Mission Elapsed Time at which the first observation of your target began.
- Coordinates
- The position of the object of interest. This is
filled in automatically if the name is resolved succesfully. Note
that, if the name could not be resolved by SIMBAD, but matches a
target name in the Swift
observation database, the RA and Dec are taken from this database. Coordinate
entry is free-format (see the acceptable coordinate formats section below).
- Try to centroid?
- Whether the software should attempt to
centroid on the source, to get the best position in the XRT
astrometric frame. This only affects the light curves and
spectra. It is normally recommended that this be set to
"yes" since the position in the XRT coordinate frame may
slightly from the object position in other frames (the XRT has an
astrometric uncertainty of 3.5" to 90% confidence).
However, the centroid algorithm may not find the source if it is
very faint. Also, if multiple sources lie in the search region,
the brightest is assumed to be the target of interest. In such
cases, it is recommended that the user examines the XRT data
themselves to obtain the best XRT position of the object, and set
this option to "no".
- Search radius
- If you set "Try to centroid?" to "yes" you must also
provide an error radius. The code will centroid on the brightest object within this
distance from the RA/Dec provided.
- Build light curve/spectrum/position/image
- Specify which products to create.
- E-mail address
- If a valid address is supplied, an e-mail will be dispatched when each of the requested products
has been produced. We recommend that you use this facility, as jobs can take a long time to
complete. The e-mail address will include the URL of your finished products.
Back to contents | Back to product generator
Choosing observations or time intervals
For each product you can choose which data get used to build the products.
The options (not all present for all products) are to use all available data,
only observations starting within some period of the first one, or to define
which observations to use. The first two cases are self-explanetory however the
latter needs some elucidation.
To define which observations are included you can either supply a list of ObsIDs
if you know them, a time range, or a mixture of the two. ObsIDs should be entered
as an 11-digit number and time ranges as start-stop. You can enter several
time intervals, separated by commas. For example, if your object had been observed all year, but you
did not want to include the summer observations, you could state 2008/01/01-2008/05/31, 2008/09/01-2008/12/31.
Valid date/time formats are detailed below. Note that, when
entering dates, only observations which both start and finish within the specified time
will be included.
Back to contents | Back to product generator
There are 4 binning options available:
- Time
- Counts
- Snapshot
- Observation
The latter two settings have no further option associated with them
(although see Common Options below), and
simply produce a single light curve and hardness ratio bin per Swift
snapshot or observation (see the definitions of
observations and snapshots). Note that the "per observation" option
assumes that the observations do not overlap in time. If this is not true, this binning
method should not be used.
For the "Time" or "Counts" binning options there are various
customisable settings, as described below.
Binning by time
This is the binning method used for most applications, where bins
are defined as being of a fixed duration. Please note that, in this
case, no checks are made to guarantee that the source is detected and
that Gaussian statistics are valid in a given bin. Please read the caveats before using this for the first time.
- Bin length (s)
- The duration of each bin. Can be different
for WT and PC mode.
- Hardness ratio bin length
- The length of each bin in the hardness
ratio. By selecting "Same as main curve" you can force this to follow
the main light curve settings.
- Min WT Counts
- Sometimes the XRT enters Windowed Timing mode for reasons other than
source brightness, such as earth limb contamination. This can result in
spurious WT data points with non-Gaussian errors which ruin the scaling of
the light curve plots. By setting this field, any WT mode bin with fewer
than the specified number of counts is assumed to be a spurious bin, and is
rejected.
Back to contents | Back to product generator
Binning by counts
This is the method used for binning GRBs and is described in detail in Evans et
al. (2007).
- Counts per bin
- The minimum number of counts necessary for a bin
to be filled.
- Min bin length
- The shortest duration a bin can have.
- Dynamic binning
-
In the standard products, the data are binned dynamically; that
is, the binning criteria vary with count rate. In this case the
"Counts per bin" field is valid at 1 count per second.
At any other count-rate the counts-per-bin value (X) is
determined dynamically thus:
When the count rate changes from 1, by some factor (the rate
factor), X changes by some other factor (the
bin factor). The default values for these factors are 10
and 1.5 respectively, thus an order-of-magnitude change in count
rate produces a factor of 1.5 change in bin size. This is done
discretely rather than continuously, so if X=20 then
where the count rate lies in the range 1-9.99 cps a bin must
contain at least 20 counts. However, when the count rate is
10-99.99 cps, there must be 30 counts for a bin to be full. At
lower rates of 0.1-0.99 cps, only 13 (=20/1.5) counts are needed
to complete a bin, although this is rounded to 15: no matter how
low the count rate goes, a bin will always need at least 15
counts, to ensure that using Gaussian methods to propagate errors
statistics is valid.
- Rate/Bin factor
- The settings to control the dynamic binning (see above).
- Minimum counts per bin
- Some minimum below which the dynamic binning
cannot drop. We suggest that you leave this at 15 (or higher) so that the assumption
that the errors are Gaussian is valid.
- Minimum sigma
- For all but the final bin, the bin is not full until
the bin's sigma value (i.e. src_cts/ sqrt(bg_counts_in_src-reg)) is
equal to or exceeds the value specified here.
Back to contents | Back to product generator
Common options
- Energy and grade selection
- The default setting is to
extract light curves using events of grade 0-12 (or 0-2 for WT
mode), between 0.3 and 10 keV. For the hardness ratio, the soft band
covers 0.3-1.5 keV and the hard band 1.5-10 keV. If this option is
set to "User specified" these numbers can be
changed. For grades there are only two options: to use the default
grades (0-12 for PC and 0-2 for WT mode), which is recommended; or only
grade 0 events. If you do not know what this means, use the default.
- Specify observations
- See Choosing observations or time intervals, above.
Back to contents | Back to product generator
- Use redshift?
- If this is set to "No" (default)
the automatic spectral fit is simply an absorbed power-law, with no components fixed.
However if you choose to supply a redshift the absorption is split into two components.
One is at redshift 0 and frozen to the Galactic column, the second is frozen at
the redshift you supplied, and the column density if free.
- Redshift
- If you said "Yes" to the above you must
supply the redshift to be used.
- Use which observations
- See Choosing
observations or time intervals, above.
- Grade range
- Whether to use events with the default grade selection
criteria (0-12 for PC, 0-2 for WT mode) or only grade 0 events. Unless you know what this
means, use the default. Note that the counts-to-flux conversion factor
deduced from this spectrum is appropriate to the grade selection used, so
if you intend to convert a count-rate light curve into flux, the light
curve must be created with the same grade selection.
- Time for spectrum
- By default, all data contained within
the selected observations are included in the spectrum. If you wish
to define time ranges, select "User defined" from the
drop-down menu. You can then enter up to four time regions. Each
region must have a name and a time region. The latter is of the
format start-stop (see accepted time
formats) and can include multiple intervals separated by commas.
Back to contents | Back to product generator
Two types of position are produced: enhanced and unenhanced. Each of
these uses a PSF fit, with bad column and pileup correction built in,
rather than the barycentric approach in the xrtcentroid
tool. For the enhanced positions the absolute astrometry is
corrected using field stars in the UVOT, this yields a systematic
uncertainty of 1.4" (90% confidence), compared to the 3.5"
systematic for the unenhanced positions, where the astrometry is
determined from the star trackers.
As well as correcting for the bad columns and any pile up, the PSF
fit has a lower statistical error than the barycentric one. In the
latter case the statistical error is calculated as
Err=22*C-0.48, where C is the number of
counts in the image. For PSF fitting the error is determined for
each fit, however simulations show it scales approximately as
Err=14.6*C-0.48.
The fields in the position form are details thus:
- Inclusion radius
- For each XRT image created as part of position determination,
source detection and centroiding is performed. The source nearest to the object position
given in the Object Details form is taken as the object of interest, provided it lies within this inclusion radius.
- Use which observations
- See Choosing
observations or time intervals, above.
- PSF Method
-
There are 3 methods which can be used to determine the position:
- Sum all data
- First snapshot only
- Multiple snapshots (default)
In most cases the third method is advised. Details of each method, and when to use
which, follow:
- Sum all data
- This is necessary for sources which are
too faint to be found in a single snapshot. A single image and
exposure map is created from all available data, and then
fitted. Note that this can be slow (as many exposure maps are
created and summed). This method is not recommended for highly
variable sources (unless they are never piled up), as the PSF of the aggregate
image in such cases will not correspond to any individual calibrated profile
which may reduce the position accuracy.
- First snapshot only
- In this method a single image
and exposure map are created from only the first snapshot of PC
mode data in which a source can be detected. This is the fastest
method. This is only recommended for bright sources (with more than
about 300 counts in the first snapshot) where the statistical error
from a single snapshot is low.
- Multiple snapshots
- This is the recommended approach
for most cases. The data are divided into snapshots, and a
centroid is performed for each snapshot. The PSF profile used is
chosen on a snapshot-by-snapshot basis accounting for
variability. The final position is formed from the weighted mean
of these positions, with the XRT systematic error added in
quadrature as a final step. If this method is chosen, the
results page will give the centroid (and image) from each
snapshot for diagnostic purposes.
A sub-option of this
method allows you to set the statistical error threshold,
discussed below.
- Statistical error threshold
- If the "Multiple
snapshots" method is chosen, you can specify the statistical
error threshold. A running measure of the the statistical error in
the weighted mean position is kept, and once it drops below the
threshold value, the process stops and the position is reported.
This is an efficiency setting: since the systematic error is
3.5" it is inefficient to use all (for example) 50 snapshots of
data for a source if the statistical error is 0.2" after three
snapshots. By default however, this value is set to 0, so all snapshots will be used.
Back to contents | Back to product generator
Images are produced in gif and FITS format.
- Energy Bands
- Here you specify the energy band(s) for your image(s), as a comma separated list
of
lowen-highen values. Energies should be in keV. The default setting is
to create 3 images, 0.3-10 keV, 0.3-1.5 keV and 1.5-10 keV.
- Use which observations
- See Choosing
observations or time intervals, above.
Back to contents | Back to product generator
While efforts have been made to allow users to enter data in any format they desire,
it is not possible to transparently handle any conceivable input. Details of the supported
input formats are below.
Although none of the time fields is mandatory, there are several
places you may wish to enter a time: the "Start time" box,
to specify which data should be included in your products, and to
generate time-resolved spectra. The same formats are acceptable in
each box, although see the important caveat for calendar dates.
There are 3 general time formats permitted: calendar
dates, [Modified] Julian Dates and Swift Mission Elapsed Time (MET). For time-resolved spectra
you can also enter the time in seconds since the "Start time".
Note that dates before the Swift epoch of 2001 January 01 are not accepted.
Calendar dates
The general form for calendar dates is year month day,
delimited with spaces, slashes (/) or hyphens (-). A time may also
be appended as hours:minutes:seconds. This should be separated from
the date either by a space, a T, @ or the word "at".
The year must be a 4 digit number, month can be either the full
month name in English, the first 3 digits of the month, or the month
number (January=1). Some examples are given below.
Important caveat: The use of hyphens as delimiters is only possible
in the "Start time" box. Any box accepting a time range uses hyphens to separate
the start and end times, so, for example 2009-09 would be interpreted as "from 2009 until 9"
which in turn would generate an error as these are not valid inputs.
- 2009-10-27
- 20091027
- 2009/10/27 at 12:32
- 2009-10-27T12:32:48.7
Back to contents | Back to product generator
Julian Dates
Julian Dates or Modified Julian Dates can simply be entered as they
are, optionally prefixed with JD or MJD.
For example 2454950.5, or 54950 are acceptable Julian Dates.
Back to contents | Back to product generator
Mission Elapsed Time (MET)
MET is the time used internally by Swift; any times you enter are ultimately converted to MET.
If you know the MET of the time in question, you can enter it directly.
Note that any purely numeric (or exponential) number entered in a time box which cannot be parsed
as a calendar date of Julian date will be assumed to be an MET value (for time-resolved spectra,
if the number is less than 108 it is interpreted as the time since the start time).
Back to contents | Back to product generator
Coordinates must be entered as RA followed by Dec, separated either by a comma or space.
Each co-ordinate can be in decimal degrees or sexagesimal. In the latter case RA is interpreted as
being in hours, not degrees. Coordinates are assumed to be in J2000.0 and will not be precessed.
Sexagesimal input must have at least minute accuracy, but apart from
this any valid format should be correctly parsed. For RA, hours,
minutes and seconds can be separated by spaces, colons or h, m or s.
Sexagesimal declination entries can be separated by colons, spaces,
d, m, ', s or ". Note that meaningless use of these delimiters
will result in an error. For example legal values for RA include
12h34m12s, 12 34.2 or 188.55,
but 12d34m would generate an error, d
cannot be used to separate RA hours from RA minutes!
Some examples of valid inputs are below. This is not exhaustive.
- 03 31 11.82 +43 54 16.8
- 03h 31m 11.82s, 43.90467
- 03:31:11.8s, 43d54'16.8''
- 52.79925, 43.90467
- 52.79925, +43:54:16.8
Back to contents | Back to product generator
Usage and acknowledgement policy
Users may publish any products created using these web tools,
provided that this facility is acknowledged. We request this be done
both by citing Evans
et al. 2009, (MNRAS, 397, 1177) when the data are introduced (if
you use an enhanced position, please also cite Goad
et al. 2007, A&A, 476, 1401), and including the following
text in the acknowledgements section of the paper (not necessary for
ATELs or circulars): This work made use of data supplied by the
UK Swift Science Data Centre at the University of Leicester.
Back to contents | Back to product generator
Caveats
General
The software used to create XRT products have been designed to work for point sources.
Source localisation, pile-up detection and correction, bad column correction etc. are
directly affected by this assumption. These tools should not be used to create
products for extended objects.
The time taken to create your products can take anything from a few minutes to a few hours,
depending on the current number of jobs running, the number of observations included in your
request, and the number of source events (for light curves). It is thus recommended that you
bookmark the "working" page to which you are taken after submitting the form. We also
suggest you supply your e-mail address, so that you are notified once your products are built,
and given the link to them. In extreme cases, for very bright sources with many observations, we
recommend that users submit multiple jobs (ideally sequentially), requesting only a few
observations in each one.
Light curves
The code used here assumes that errors can be propagated according to
Gaussian statistics -- except for the final bin of the data (a holdover from
the GRB code, which is useful and so kept). Thus, if the bin size is such
that some bins contain fewer than 15 counts, the error values may well
be wrong. The onus is on the user to specify sufficient counts that
the statistics are Gaussian. A warning will be given on the light curve
results page if any bins contain fewer than 15 counts.
The maximum number of bins permitted in a single light curve is 500,000. When a light curve
is selected with constant-duration bins, the total exposure time in each mode is divided by the binsize.
If either mode exceeds half a million bins, the light curve will not be built. In such cases users can either select
larger time bins, or create several light curves, each using less than the full set of observations, and then combine
these manually to give the complete light curve.
Back to contents | Back to product generator
Change log
- 2009 December 14 - Unenhanced positions and images
added to the list of available products, and the user interface redesigned.
- 2009 November 4 - The spectrum and light curve
generators were upgraded. Now (if "centroid" is set to
"yes") a new centroid is determined for each snapshot with
more than 40 source counts. If the source is fainter than this the
exposure-weighted mean position determined so far is used. This
improves the products in the rare cases where the attitude
information is of low accuracy. A second change was made at the same
time, whereby the definition of the PSF wings used for pileup
determination is brightness dependent. This improves pile-up
correction in very bright sources (i.e. where XRT would normally be
in WT mode but has been put in PC mode at the user's request).
- 2009 May 29 - Fixed a bug in the spectrum generator that affects pile-up correction.
For objects which are piled up at the start of a snapshot, the piled-up data were included twice, once
corrected for pileup and once not. For GRBs this makes very little difference to the spectral shape
(hence it was not detected when testing this software). For non-GRB objects we urge users to check their spectra
with the corrected software. We apologise for this problem.
Back to contents | Back to product generator