def meqmaker_compile_options():
meqmaker = MeqMaker.MeqMaker(solvable=True,
use_jones_inspectors=True,
use_skyjones_visualizers=False,
use_decomposition=False)
if TiggerSkyModel is not None:
tsm = [TiggerSkyModel()]
else:
tsm = []
meqmaker.add_sky_models(tsm +
[central_point_source, fitsimage_sky, gridded_sky])
meqmaker.add_sky_jones('dE', 'differential gains', [
solvable_sky_jones.DiagRealImag('dE'),
solvable_sky_jones.FullRealImag('dE'),
solvable_sky_jones.DiagAmplPhase('dE')
])
meqmaker.add_uv_jones('P', 'feed orientation', [feed_angle])
meqmaker.add_uv_jones('B', 'bandpass', [
solvable_jones.DiagRealImag('B'),
solvable_jones.FullRealImag('B'),
solvable_jones.DiagAmplPhase('B')
])
meqmaker.add_uv_jones('G', 'receiver gains/phases', [
solvable_jones.DiagRealImag('G'),
solvable_jones.FullRealImag('G'),
solvable_jones.DiagAmplPhase('G')
])
meqmaker.add_vis_proc_module('IG', 'multiplicative IFR errors',
[ifr_based_errors.IfrGains()])
meqmaker.add_vis_proc_module('IC', 'additive IFR errors',
[ifr_based_errors.IfrBiases()])
TDLCompileOptions(*meqmaker.compile_options())
return meqmaker
meqmaker.add_uv_jones('P', 'feed orientation', [feed_angle])
# B - bandpass, G - gain
from Calico.OMS import solvable_jones
meqmaker.add_uv_jones(
'B', 'bandpass',
[solvable_jones.DiagAmplPhase("B"),
solvable_jones.FullRealImag("B")])
meqmaker.add_uv_jones(
'G', 'receiver gains/phases',
[solvable_jones.DiagAmplPhase("G"),
solvable_jones.FullRealImag("G")])
from Calico.OMS import ifr_based_errors
meqmaker.add_vis_proc_module('IG', 'interferometer gains',
[ifr_based_errors.IfrGains()])
meqmaker.add_vis_proc_module('IC', 'interferometer biases',
[ifr_based_errors.IfrBiases()])
# very important -- insert meqmaker's options properly
TDLCompileOptions(*meqmaker.compile_options())
import Purr.Pipe
def _define_forest(ns, parent=None, **kw):
if run_purr:
Timba.TDL.GUI.purr(mssel.msname + ".purrlog", [mssel.msname, '.'])
# create Purr pipe
global purrpipe
purrpipe = Purr.Pipe.Pipe(mssel.msname)