Next: src_collect
Up: Observation group level analysis
Previous: Observation group level analysis
Contents
j_ima_mosaic
This executable combines input sky maps produced by j_ima_iros to make JEM-X
mosaic images in a single FITS file (named by default jmxi_mosa_ima.fits)
comprising each energy band found in the input observation group or selected
by the user.
The input intensity IMATYPE (RECONSTRUCTED or RECTIFIED or RAWINTENSITY) is
selected by use of the parameter "mapSelect". Note that if the user wants to extract fluxes from the image, only "RECONSTRUCTED"
(recommended) or "RECTIFIED" make sense. There are four parameters
defining the mosaic skymaps produced in the output file so it can have, for
each energy band, up to four image extensions (depending on the chosen
options) showing in the following respective order:
- RECONSTRUCTED (default) or RECTIFIED or RAWINTENSITY (depending on the
input intensity IMATYPE; see above): it is the weighted intensity map
(counts/dm2/s) produced by use of the parameter "viewIntens",
- VARIANCE: the weighted variance map (use "viewVar"),
- SIGNIFICANCE: the
significance map (use "viewSig") defined as the intensity divided by the
square root of the variance at each pixel.
- EXPOSURE: the effective
exposure times (use "viewTime") expressed in seconds taking account of the dead
time, grey filter effects and vignetting of the JEM-X instruments (Notice: the
vignetting used is normalised to 1 at the center of JEM-X FoV so as this
corresponds to an effective detector aera of about 70 cm2 in average,
depending on the actual energy range). With the latter option an additional map
is also produced in the last extension OBS_TIME of the output file that
shows the accumulated raw observation times (in seconds).
If sky coordinates (RAcenter, DECcenter) are given as parameters to define the
center of the mosaic map, an ordered list of the input images as a function of
their respective weights at the given coordinates is displayed for each energy
band in the log file.
Since OSA10 it is possible (switching the new parameter view_nb to ``yes'') to
output in a separate extension the number of SCWs contributing to each pixel of the mosaic.
This new image will be of type NB_SING, and will be the last extension of the mosaic fits file.
j_ima_mosaic can combine any set of JEM-X mosaics,
even if they have been obtained with different units. It is however
mandatory that the user must attach all the requested
mosaics to a group. For instance, if the user wants combine mosaics
from the three OSA analyses: mysrc_1_jmx1, mysrc_2_jmx1, mysrc_3_jmx2,
the user shall type:
dal_create mos_group.fits ""
dal_attach mos_group.fits obs/mysrc_1_jmx1/jmx1_mosa_ima.fits
obs/mysrc_2_jmx1/jmx1_mosa_ima.fits
obs/mysrc_3_jmx2/jmx2_mosa_ima.fits "" ""
("" in the first line indicates that we use the default template,
which is a group; the final "" "" represents two empty child
parameters of dal_attach; indeed dal_attach can attach up to five
children in a single invocation).
Now the group is ready, and one can call j_ima_mosaic:
j_ima_mosaic inObsGrp=mos_group.fits outfile=combined_mosaic.fits moscomb=y
The moscomb parameter is very important, as it tells j_ima_mosaic
that we are trying to combine mosaics. The user can use either JMX1 or JMX2.
j_ima_mosaic can now produce mosaics either in the default tangential projection
or using AIToff-Hammer projection (useful in particular when mosaicking large parts of the sky, such as for
the Galactic Plane Scans).
AIToff-Hammer projection can be selected with the parameter IMA2_AITproj="yes"
.
Table 23:
j_ima_mosaic parameters
|
|
|
Name |
Type |
Description |
inObsGrp |
string |
DOL of observation group
default `` '' |
mapSelect |
string |
Default for Reconstructed (Residual+Sources)
Intensity maps, else RECTI for Rectified Intensity maps or
RAWIN for Raw Intensity maps
default ``RECON'' |
dolBPL |
string |
DOL of the vignetting maps,
e.g., ../../ic/jmx1/rsp/jmx1_bpl_grp_0002.fits[JMX1-DMAP-BPL]
default `` '' |
radiusSelect |
real |
Input images selection radius in deg
default ``5'' |
eminSelect |
real |
Minimum energy (keV) to select energy
interval
default ``0'' |
emaxSelect |
real |
Maximum energy (keV) to select energy
interval
default ``80'' |
diameter |
real |
Diameter in degrees of the mosaic image;
0 allows to fit the size from the angle between input images,
: may be used to force larger mosaics
default ``0'' |
cdelt |
real |
Pixel size in degrees at mosaic center
default ``0.026'' |
RAcenter |
real |
Center of the mosaic image Right Ascension in degrees
(if
use the computed middle position of the mosaic)
default ``-1'' |
DECcenter |
real |
Center of the mosaic image Declination
in degrees (not used if RAcenter
)
default ``0'' |
outfile |
string |
Mosaic FITS file prefix name (optional)
default `` '' |
viewTime |
boolean |
Create total exposure Time map?
(if Y: produce both raw exposure time map and energy dependent
effective exposure maps)
default ``Y'' |
viewIntens |
boolean |
Create raw Intensity map?
default ``Y'' |
viewVar |
boolean |
Create Variance map?
default ``Y'' |
viewSig |
boolean |
Create Significance map?
default ``Y'' |
view_nb |
boolean |
Create a map that shows the number of
input images contributing to each pixel?
default ``N'' |
AITproj |
boolean |
Produce map in galactic coordinates using
AIToff-Hammer projection?
default ``N'' |
print_ScWs |
boolean |
List input Science Windows in mosaic headers?
default ``N'' |
chatter |
integer |
Level of test output
default ``2'' |
Next: src_collect
Up: Observation group level analysis
Previous: Observation group level analysis
Contents
inthelp@sciops.esa.int