Safer default facetselfcal config for delay calibration
For complicated sources the current default settings for facetselfcal are sub-optimal. I've almost never not used it, so I think it is safe to always go for the "reset" approach where phases are first solved on the international stations and the remote stations are subsequently solved for on longer timescales or with a higher smoothness constraint. The screenshots below compare the 9th iteration on a couple of LBCS candidates between the current approach on the left and the reset approach on the right, all starting from identical point source models. It seems to consistently match or outperform the current default. ![Screenshot_from_2024-11-28_15-23-28](/uploads/5e7f612c3d3b42ce736071d782baf9a7/Screenshot_from_2024-11-28_15-23-28.png) ![Screenshot_from_2024-11-28_15-24-18](/uploads/d50e4ee8b90b0fd04575328ee0cd934f/Screenshot_from_2024-11-28_15-24-18.png) ![image](/uploads/eafe604ba0236459574aefd473891054/image.png) ![image](/uploads/89b36183f45c3624643d60d5833033da/image.png) ![image](/uploads/ecb2ad2c7f686ce5bb82603c2d01f651/image.png) **Default config** ``` imsize = 1600 pixelscale = 0.075 robust = -0.5 uvmin = 40000 maskthreshold = [7.0] soltype_list = ['scalarphasediff','scalarphase','scalarcomplexgain'] soltypecycles_list = [0,0,2] solint_list = [4,1,100] nchan_list = [1,1,1] smoothnessconstraint_list = [10.0,2.0,10.0] antennaconstraint_list = ['alldutch',None,None] docircular = True forwidefield = True stop = 10 beamcor = "auto" phaseupstations = "core" ``` **Suggested "reset" config** ``` imsize = 1600 pixelscale = 0.075 robust = -0.5 uvmin = 40000 maskthreshold = [7.0] soltype_list = ['scalarphasediff','scalarphase','scalarphase','scalarcomplexgain'] soltypecycles_list = [0,0,0,2] solint_list = [4,1,1,100] nchan_list = [1,1,1,1,] smoothnessconstraint_list = [10.0,2.0,10.0,10.0] antennaconstraint_list = ['alldutch',None,None,None] resetsols_list = [None,'alldutch',None,None] docircular = True forwidefield = True stop = 10 beamcor = "auto" phaseupstations = "core" ``` I can open a MR for further discussion or testing.
issue