`mkexpmap' generates an exposure map which may be used
to convert a counts image of a source to an image in
flux units. An image in sky-coordinates, the exposure
map is essentially an image of the effective area
associated with each sky position as a result of the
aspect motion of the telescope over the course of an
observation. It accounts for the effects of dither
motion which are especially important near the edges of
the detector and also near bad pixels and bad columns.
The default units of the exposure map (M) are
M = [cm**(2)*counts/photon]
and are the same units as the effective area (ARF).
Optionally, the exposure time factor (t) may be
included, yielding a map
M' = M*t = [sec*cm**(2)*counts/photon].
To include the exposure time factor in the output map,
set normalize="no".
To compute an image in photon flux
units, divide the counts image (N) [counts/pixel] by
the product of the exposure map (M) and the exposure time (t):
photon flux [photons/cm**2/s/pixel] = N / (M*t) = N / M'
Dmimgcalc may be used to carry out this division.
As input, 'mkexpmap' requires an aspect histogram file
(computed by asphist) and an instrument map file
(computed by mkinstmap). The output exposure map is
generated for a specified sky-coordinate grid.
The input aspect histogram gives the total dwell time
for each sky pointing direction accumulated by the
telescope over the course of the observation.
The input instrument map is an image in detector
coordinates which, for a given incident spectrum, is
essentially the product of the mirror effective area
projected onto the detector surface with the QE of the
detector. It also incorporates bad pixels and ACIS
windows (if any) by setting the efficiency to zero at
the appropriate locations on the detector.
The exposure map is a dwell-time-weighted sum of
overlapping instrument maps, with one term in the sum
for each pointing orientation specified
in the aspect histogram. Because this repeated
shifting and adding of instrument map images can
require an extremely large number of floating-point
operations, the exposure map calculation may take a
long time. The running time scales roughly as the
product of the number of aspect histogram records and
the number of sky pixels. In cases where computation
speed is more important than a detailed treatment of
the pointing motion, one can greatly speed up the
computations by computing an exposure map for the
average pointing direction; to select this option,
set useavgaspect=yes.
For a detailed, quantitative definition of the exposure
map, see Davis, John E. 2001, ApJ, 548, 1010.
See also
http://space.mit.edu/CXC/docs/expmap_intro.ps.gz
.