lofar-astron / factor

Facet calibration for LOFAR
http://www.astron.nl/citt/facet-doc
GNU General Public License v2.0
19 stars 12 forks source link

Verifying facet run #199

Closed duyhoang-astro closed 7 years ago

duyhoang-astro commented 7 years ago

Hi, Factor said that it has successfully finished one of the facets. The max residuals in the widefield images (*.wsclean_pre-image.fits and .wsclean_post-image.fits) actually decrease (from 0.03216 to 0.03197). However, the max residuals in the facet region significantly increase (from 0.01497 to 0.03197), as we can see in the images below.

pre-residual facet region image post-residual facet region image

I guess Factor checks for the max residuals in the whole widefield residual image (as the old facet scripts do). But in this case it seems not to work properly. An additional check is probably needed to compare max residuals within both the facet and the widefield images?

P/S: perhaps, a not related question to the thread: how can I undo this direction (as we do in the facet scripts)? I can leave this facet here, let factor run on other facets, then come back to this facet later. The option -r RESET, --reset=RESET seems to both undo and rerun the facet. Thanks! Duy

AHorneffer commented 7 years ago

To me it looks like some weaker sources were not subtracted, probably because they were not within the masked region. If that is indeed the case, then you could supply a region file to Factor with these sources in there.

If you just want to rest one direction, but not selfcal it again, then you can do it the following way:

  1. re-run factor with the -r option
  2. stop it after the facetsubrest operation is done and before the following faceselfcal operation finishes. (Because a faceselfcal usually takes quite some time this shouldn't be too hard.)
  3. change the order of the directions in your directions file and move the direction, that you don't want to selfcal further down.
  4. re-start factor
duyhoang-astro commented 7 years ago

@AHorneffer: Thanks for the tips. I reseted the facet and will try this facet later when the residuals are improved.