Here we demonstrate how to extract an image over various energy bands. The PC data from the first observation of GRB 091029 are used as an example.
>xselect
** XSELECT V2.2a **
> Enter session name >[xsel]
xsel:ASCA >read event
> Enter the Event file dir >[./]
> Enter Event file list >[ ] sw00374210000xpcw3po_cl.evt
Got new mission: SWIFT
> Reset the mission ? >[yes]
Notes: XSELECT set up for SWIFT
Time keyword is TIME in units of s
Default timing binsize = 5.0000
Setting...
Image keywords = X Y with binning = 1
WMAP keywords = X Y with binning = 1
Energy keywords = PI with binning = 1
Getting Min and Max for Energy Column...
Got min and max for PI: 0 1023
Got the minimum time resolution of the read data: 2.5073
MJDREF = 5.1910000742870E+04 with TIMESYS = TT
Number of files read in: 1
******************** Observation Catalogue ********************
Data Directory is: /home/work/kpa/2009/grb091029/00374210000-xrt/
HK Directory is: /home/work/kpa/2009/grb091029/00374210000-xrt/
OBJECT OBS_ID DATE-OBS DATAMODE
1 GRB091029 ( 00374210991 2009-10-29T PHOTON
> extract image
extractor v5.11 11 May 2009
Getting FITS WCS Keywords
Doing file: /home/work/kpa/2009/grb091029/00374210000-xrt/sw00374210000xpcw3po_cl.evt
100% completed
Total Good Bad: Time Phase Grade Cut
5858 5858 0 0 0 0
===============================================================================
Grand Total Good Bad: Time Phase Grade Cut
5858 5858 0 0 0 0
in 15790. seconds
Image has 5858 counts for 0.3710 counts/sec
> plot image
The image above was extracted for the full energy range of the XRT: about 0.2-10 keV. However, the data can be filtered to produce images over any energy band within this range, using the command filter pha_cutoff, where the range must be given in terms of channels. For the XRT, 1 channel = 10 eV.
> filter pha_cutoff 30 100 > extract image > save image 0.3-1keV.im > clear pha_cutoff > filter pha_cutoff 100 1000 > extract image > save im 1-10keV.im
The commands above extract and save images over soft (0.3-1 keV) and hard (1-10 keV) bands for the full field of view. Filtering on regions and time can also be performed, as explained in the Introduction to XSELECT thread.