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
'primary beam', [wsrt_beams],
pointing=solvable_pointing_errors)
# then add differential gains
from Calico.OMS import solvable_sky_jones
meqmaker.add_sky_jones('dE', 'differential gains',
[solvable_sky_jones.FullRealImag('dE')])
# P - feed angle
from Siamese.OMS import feed_angle
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())
lsm = Meow.LSM.MeowLSM(include_options=False)
meqmaker.add_sky_models([lsm, fitsimage_sky, gridded_sky])
from Lions import ZJones
meqmaker.add_sky_jones('Z', 'iono', [ZJones.ZJones()])
from Calico.OMS import solvable_sky_jones
meqmaker.add_sky_jones('GD', 'directional receiver gains/phases', [
solvable_sky_jones.DiagAmplPhase("GD"),
solvable_sky_jones.FullRealImag("GD")
])
from Calico.OMS import solvable_jones
meqmaker.add_uv_jones(
'G', 'receiver gains/phases',
[solvable_jones.DiagAmplPhase("G"),
solvable_jones.FullRealImag("G")])
meqmaker.add_uv_jones('B', 'bandpass', [solvable_jones.DiagAmplPhase("B")])
# simulate G Jones
from Siamese.OMS import oms_gain_models
meqmaker.add_uv_jones('G_sim', 'simulate receiver gains/phases',
[oms_gain_models])
# very important -- insert meqmaker's options properly
TDLCompileOptions(*meqmaker.compile_options())
def _define_forest(ns):
# make pynodes, xyzcomponent for sources
