NAME j_ima_mosaic v7.0.2 - Produce weighted mosaics of JEM-X images USAGE Call with at least the parameter "inObsGrp" defined to the current observation group DOL. j_ima_mosaic inObsGrp=og_jmx1.fits+1 DESCRIPTION Version 7.0.2/20070824/JC 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". 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: 1) [RECONSTRUCTED] (default) or [RECTIFIED] or [RAWINTENSITY] (depending on the input intensity IMATYPE; see above): it is the weighted intensity map (counts/cm2/s) produced by use of the parameter "viewIntens", 2) [VARIANCE]: the weighted variance map (use "viewVar"), 3) [SIGNIFICANCE]: the significance map (use "viewSig") defined as the intensity divided by the square root of the variance at each pixel. 4) [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). In order to combine mosaics, e.g. from JEM-X1 and JEM-X2, it is possible to make a mosaic of mosaics by use of the parameter "moscomb". Caveat: the input mosaics must have the same dimensions, which is the case if they have been produced by using the same values for the parameters "diameter" and "cdelt". In the current version the output data structure of the combined mosaic is always JMX1-MOSA-IMA even if JEM-X2 data have been involved (SCREW 2008). The list of input mosaics is attached as an index in a new group read by the parameter inObsGrp. To get the largest possible output mosaic the parameters "radiusSelect" and "diameter" must be set to -1. 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. PARAMETERS inObsGrp, s,h, , ,, "DOL of observation group" mapSelect, s,h, "RECON", ,, "Default for Reconstruced (Residal+Sources) Intensity maps, else RECTI for Rectified Intensity maps or RAWIN for Raw Intensity maps" dolBPL, s,h, "", ,, "DOL of the vignetting maps, e.g., ../../ic/jmx1/rsp/jmx1_bpl_grp_0002.fits[JMX1-DMAP-BPL]" radiusSelect, r,h, 5., 0.,6.5, "Input images selection radius in deg." eminSelect, r,h, 0., 0.,80., "Minimum energy (keV) to select energy interval" emaxSelect, r,h, 80., 0.,80., "Maximum energy (keV) to select energy interval" diameter, r,h, 0, -1,90., "Diameter in degrees of the mosaic image; 0 allows to fit the size from the angle between input images, <0: may be used to force larger mosaics" cdelt, r,h, 0.03, 0.02,1., "Pixel size in degrees at mosaic center" RAcenter, r,h, -1, -1,360., "Center of the mosaic image Right Ascension in degrees (if <0 use the computed middle position of the mosaic)" DECcenter, r,h, 0., -90.,90., "Center of the mosaic image Declination in degrees (not used if RAcenter <0)" outfile, s,h, "", ,, "Mosaic FITS file prefix name (optional)" viewTime, b,h, Y, ,, "Create total exposure Time map? Y|N (if Y: produce both raw exposure time map and energy dependent effective exposure maps)" viewIntens, b,h, Y, ,, "Create raw Intensity map? Y|N" viewVar, b,h, N, ,, "Create Variance map? Y|N" viewSig, b,h, Y, ,, "Create Significance map? Y|N" print_ScWs, b,h, N, ,, "List input ScWs in mosaic headers? Y|N" moscomb, b,h, N, ,, "Combine mosaics? Y|N" chatter, i,h, 2, 0,5, "Level of test output" EXAMPLES Default usage : Combines all j_ima_iros intensity images reconstructed with sources (IMATYPE: RECONSTRUCTED) to produce, for each energy band found, weighted intensity and significance mosaic maps as well as exposure time maps in the file jmx1_mosa_ima.fits (default name in the case of JEM-X1): j_ima_mosaic inObsGrp=og_jmx1.fits+1 dolBPL="ic/jmx1/rsp/jmx1_bpl_grp_0004.fits[JMX1-DMAP-BPL]" To select IROS Intensity images of type (IMATYPE) RECTIFIED: j_ima_mosaic inObsGrp=og_jmx1.fits+1 mapSelect=RECTI dolBPL="ic/jmx1/rsp/jmx1_bpl_grp_0004.fits[JMX1-DMAP-BPL]" To select only a specific energy band inside a given interval (in keV): j_ima_mosaic inObsGrp=og_jmx1.fits+1 eminSelect=5. emaxSelect=11. dolBPL="ic/jmx1/rsp/jmx1_bpl_grp_0004.fits[JMX1-DMAP-BPL]" To give the mosaic output file another name (e.g. "JEMX_MOSA_IMA.fits"): j_ima_mosaic inObsGrp=og_jmx1.fits+1 outfile=JEMX_MOSA_IMA dolBPL="ic/jmx1/rsp/jmx1_bpl_grp_0004.fits[JMX1-DMAP-BPL]" To produce arbitrary sized mosaics (can lead to very large files): j_ima_mosaic inObsGrp=og_jmx1.fits+1 diameter=-1 dolBPL="ic/jmx1/rsp/jmx1_bpl_grp_0004.fits[JMX1-DMAP-BPL]" To make a zoom on a given position (and obtain an ordered list of the ScWs weights at that position): j_ima_mosaic inObsGrp=og_jmx1.fits+1 radiusSelect=3. diameter=5. RAcenter=123.45 DECcenter=67.89 dolBPL="ic/jmx1/rsp/jmx1_bpl_grp_0004.fits[JMX1-DMAP-BPL]" To produce the weighted variance map but avoid the significance map: j_ima_mosaic inObsGrp=og_jmx1.fits+1 viewVar=Y viewSign=N dolBPL="ic/jmx1/rsp/jmx1_bpl_grp_0004.fits[JMX1-DMAP-BPL]" To combine JEM-X1 mosaics with JEM-X2 mosaics: 1) dal_create mosa_grp.fits dal_attach mosa_grp.fits obs/J0039_1/jmx1_mosa_ima.fits obs/J0039_2/jmx2_mosa_ima.fits obs/J0102_1/jmx1_mosa_ima.fits obs/J0102_2/jmx2_mosa_ima.fits 2) j_ima_mosaic inObsGrp=mosa_grp.fits radiusSelect=-1 diameter=-1 moscomb=y outfile=jmx_1p2_mosa.fits NB: For backward compatibility imperatives it is still possible to combine input images produced by older OSA versions as well as to use the old values of mapSelect: "RAW_RECT", "RES+SRC", or "RAW_INT". ERROR CODES -2002 : inObsGrp missing. -2005 : Cannot create output FITS file if it already exits. BUGS Please report problems to author. AUTHOR Jerome Chenevez, Danish National Space Center jerome@dnsc.dk SEE ALSO j_ima_iros This software is frequently updated. Please check latest updates available at: ftp://ftp.spacecenter.dk/pub/jemx/jerome/J_IMA_MOSAIC