[an error occurred while processing this directive]
Introduction This thread describes how to process EPIC mosaic-mode observations and create merged background subtracted and exposure corrected images combining data from all subpointings.This thread uses data from the observations of Fermi Bubble: ObsID 0672401301. These data are also used as the Mosaic Mode Analysis example where a script and output files are provided. Expected Outcome The final outcome of this thread are images of the merged data of the Fermi Bubble observation.SAS Tasks to be Used
Prerequisites Useful Links
|
For this example we will assume the directory structure:
/fermi/0672401301/
for the directory where the initial processing of the data will take place
with the ODF data in the directory:
/fermi/0672401301/odf
epchain
epchain withoutoftime=true
emchain
pn-filter
mos-filter
Note that mos-filter will indicate which CCDs are operating in an anomalous mode (the line with the corresponding CCD would be marked by **** in that case. In the the example below all CCDs are fine). These CCDs should be excluded in downstream processing.
1S001
Limit=1.0 CCD=2 Hardness=6.235 Uncertainty=0.815
Limit=1.0 CCD=3 Hardness=6.679 Uncertainty=0.957
Limit=0.90 CCD=4 Hardness=3.807 Uncertainty=0.410
Limit=1.84 CCD=5 Hardness=4.843 Uncertainty=0.496
Limit=1.0 CCD=6 Hardness=0.000 Uncertainty=0.000 ****
Limit=1.0 CCD=7 Hardness=7.094 Uncertainty=0.947
2S002
Limit=1.0 CCD=2 Hardness=3.116 Uncertainty=0.296
Limit=1.0 CCD=3 Hardness=5.197 Uncertainty=0.651
Limit=1.0 CCD=4 Hardness=3.966 Uncertainty=0.410
Limit=1.32 CCD=5 Hardness=5.416 Uncertainty=0.625
Limit=1.0 CCD=6 Hardness=5.297 Uncertainty=0.722
Limit=1.0 CCD=7 Hardness=4.739 Uncertainty=0.544
In this case the only CCD not useful is MOS1 CCD#6, which was lost to a micrometeorite hit.
The diagnostic plots created by the filtering tasks should also be examined as they provide an indication of the quality of the data. These have names like mos1S001-hist.qdp, mos2S002-hist.qdp, and pnS003-hist.qdp and can be plotted using the command, for example, qdp mos1S001-hist.qdp. Figure 1 shows the light-curve screening for the MOS1 instrument.
Figure 1: MOS1 light curve for the Fermi Bubble mosaic mode observation. This was a particularly clean observation.
However, before running emosaic_prep the SAS_ODF environment variable needs to be set so it points directly at the *.SAS file created by odfingest:
setenv SAS_ODF odf/2154_0672401301_SCX00000SUM.SAS
emosaic_prep mos1evtfile='mos1S001-clean.fits' mos2evtfile='mos2S002-clean.fits' \
pnevtfile='pnS003-clean.fits' atthkfile='atthk.fits'
The data from the subpointings are placed in directories with the nomenclature:
/fermi/0672401301/prep_mosaic_001
/fermi/0672401301/prep_mosaic_002
...
To avoid conflicts, delete the MOS files created in the first run. They will be recreated in the second run.
rm -f prep_mosaic*/mos*-clean*.ds
Run the task emosaic_prep the second time, this time with the PN OOT event file.
emosaic_prep mos1evtfile='mos1S001-clean.fits' mos2evtfile='mos2S002-clean.fits' \
pnevtfile='pnS003-clean-oot.fits' atthkfile='atthk.fits'
Run the task esprep_mosaic_time to extract the time intervals for the individual pointings (written to the file mosaic-times.txt), which are used in esas-mosaic to correct all the pointing information.
esprep_mosaic_time atsfile=../odf/2154_0672401301_SCX00000ATS.FIT \
sasfile=../odf/2154_0672401301_SCX00000SUM.SAS
Run the task esas-mosaic to correct the header information of the subpointings. esas-mosaic will step through the subpointing directories and run the task esprep_mosaic to make the corrections.
esas-mosaic sublist=mosaic-times.txt m1prefix=1S001 m2prefix=2S002 pnprefix=S003
/fermi/0672401301/merge
/fermi/0672401301/odf
/fermi/0672401301/prep_mosaic_001
/fermi/0672401301/prep_mosaic_002
/fermi/0672401301/prep_mosaic_003
/fermi/0672401301/prep_mosaic_004
/fermi/0672401301/prep_mosaic_005
/fermi/0672401301/prep_mosaic_006
/fermi/0672401301/prep_mosaic_007
/fermi/0672401301/prep_mosaic_008
/fermi/0672401301/prep_mosaic_009
/fermi/0672401301/prep_mosaic_010
/fermi/0672401301/prep_mosaic_011
/fermi/0672401301/prep_mosaic_012
where /fermi/0672401301/merge is the directory for the merged image processing, the prep_mosaic_* directories contain the individual subpointing image processed data. The ODF data are in the /fermi/0672401301/odf directory. For the processing point the environmental variable $SAS_CCF to the file created in the initial processing: /fermi/0672401301/proc/ccf.cif
Figure 2: Merged EPIC adaptively smoothed images of
the mosaic-mode observation of the Fermi Bubble.
Caveats
For information about known issues, refer to the ESAS Warnings and Watchouts pages. |