Self-calibration

WSClean can be used to perform self-cal. There are two main approaches, supporting a range of scenarios:

  • Use WSClean to fill the MODEL_DATA with predicted visibilities from an imaging run. When using -mgain with a value less than 1, WSClean will fill the column. This column can then be used during calibration (e.g. with DP3). The advantage of this is that it is relatively fast and easy.

  • Make WSClean output a source component list (see Component lists) and use this list during the calibration (DP3 supports this format). The advantage of this approach is that it is easier to apply the beam (and other effects) on a source list, it is very accurate, and it is possible to prune/edit the source list. A disadvantage can be that if the number of sources is very large, calibration will be very slow.

Image-based self-calibration

One scenario where image-based self-cal can be useful is when combined with the ability to do (almost-)full-sky imaging, since in some situations MWA’s field of view (for example) might require including highly off-axis sources in the calibration model.

To perform self-cal, you need to use major iterations, as these will fill/update the MODEL_DATA column. When using Casa’s tasks for calibration, the calibration tasks will use the MODEL_DATA and calibrate the CORRECTED_DATA using this column. Other calibration utilities, like MWA’s mitchcal or LOFAR’s DP3 can similarly be instructed to calibrate using the MODEL_DATA column.

Filling the MODEL_DATA requires a setting of mgain < 1, e.g. an mgain of 0.9. As long as mgain is not 1, WSClean will end with a major iteration, and the MODEL_DATA column will be set to the “best” model from the cleaning model.

Self-calibration from existing image

The ‘-predict’ option can be used to fill the MODEL_DATA column with a prediction from a pre-existing image (see prediction). After having predicted model visibilities, these visibilities can be used to calibrate the data (e.g. with DP3).

Polarized imaging & calibration

You can self-cal on Stokes I or on multiple polarizations, where the latter is more accurate (at least in the case of the MWA), but takes more time. If you run WSClean on the desired polarizations one by one, e.g. on XX and then on YY, or joinedly clean them, the MODEL_DATA column will have all imaged polarizations correctly filled in. The first run will create the MODEL_DATA column, set all values to zero and then fill the XX column, the second run will notice the MODEL_DATA already exists, and only update the YY column. B. McKinley has used this method and did a few self-cal loops to create very deep and well-calibrated Fornax A images. You can run all polarizations at once with -pol xx,xy,yx,yy, -pol iquv, or -pol rr,rl,lr,ll. For info on polarimetric deconvolution settings, see polarimetry deconvolution.

CASA on-the-fly mode

Certain CASA commands (e.g. ft) will put keywords in a measurement set that turn on the “on-the-fly” (otf) mode. In OTF mode, CASA will ignore the MODEL_DATA column and use other keywords to determine the model data. To make use of the MODEL_DATA afterwards, you can use the delmod CASA command to disable OTF mode:

delmod(vis=’myobs.ms’,otf=True,scr=False)

WSClean will never use or change OTF keywords in the measurement set.

Next chapter: Image weighting