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 = MeqMaker.MeqMaker(solvable=do_solve and run_option == 'calibrate')
# 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])
# E - beam
# add a fixed primary beam first
from Calico.OMS.wsrt_cos3_beam import WSRTCos3Beam #,wsrt_beams_zernike
meqmaker.add_sky_jones('E', 'WSRT cos^3 beam',
[WSRTCos3Beam("E", solvable=False)])
meqmaker.add_sky_jones('Es', 'solvable WSRT cos^3 beam',
[WSRTCos3Beam("Es", solvable=True)])
## add solvable refraction
# from Calico.OMS import solvable_position_shifts
# meqmaker.add_sky_jones('R','position shifts',solvable_position_shifts);
# add differential gains
from Calico.OMS import solvable_sky_jones
meqmaker.add_sky_jones('dE', 'differential gains', [
solvable_sky_jones.DiagRealImag('dE'),
solvable_sky_jones.FullRealImag('dE'),
solvable_sky_jones.DiagAmplPhase('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.DiagRealImag("B"),
solvable_jones.FullRealImag("B"),
solvable_jones.DiagAmplPhase("B")
],
flaggable=True)