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
[TiggerSkyModel(), central_point_source, fitsimage_sky])
# now add optional Jones terms
# these will show up in the menu automatically
# E - beam
# add a fixed primary beam first
from Calico.OMS import wsrt_beams
from Calico.OMS import solvable_pointing_errors
meqmaker.add_sky_jones('E',
'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")])
# specify available sky models
# these will show up in the menu automatically
from Siamese.OMS import fitsimage_sky
from Lions import gridded_sky
import Meow.LSM
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])