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
# these will show up in the menu automatically
from Siamese.OMS import gridded_sky
from Siamese.AGW import azel_sky
from Siamese.OMS import transient_sky
from Siamese.OMS import fitsimage_sky
## OMS: time to retire this one
#import Meow.LSM
#lsm = Meow.LSM.MeowLSM(include_options=False);
models = [gridded_sky, azel_sky, transient_sky, fitsimage_sky]
# ,lsm ];
try:
from Siamese.OMS.tigger_lsm import TiggerSkyModel
models.insert(0, TiggerSkyModel())
except:
print 'Failure to import TiggerSkyModel module'
print 'Is the location of Tigger defined in your PYTHONPATH environment variable?'
pass
meqmaker.add_sky_models(models)
# now add optional Jones terms
# these will show up in the menu automatically
# Ncorr - correct for N
from Siamese.OMS import oms_n_inverse
meqmaker.add_sky_jones('Ncorr', 'n-term correction', oms_n_inverse)
# Z - ionosphere
def enable_calibration(enable):
meqmaker.use_decomposition_opt.show(not enable)
if enable:
meqmaker.use_decomposition_opt.set(False)
cal_toggle.when_changed(enable_calibration)
# specify available sky models
# these will show up in the menu automatically
from Calico.OMS import central_point_source
from Siamese.OMS import fitsimage_sky
from Siamese.OMS.tigger_lsm import TiggerSkyModel
meqmaker.add_sky_models(
[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')])
# now load optional modules for the ME maker
from Meow import MeqMaker
meqmaker = MeqMaker.MeqMaker()
# specify available sky models
# these will show up in the menu automatically
from Siamese.OMS import gridded_sky
from Siamese.OMS import transient_sky
from Siamese.OMS import fitsimage_sky
import Meow.LSM
from Siamese.OMS.tigger_lsm import TiggerSkyModel
lsm = Meow.LSM.MeowLSM(include_options=False)
meqmaker.add_sky_models([
gridded_sky, transient_sky, fitsimage_sky, lsm,
TiggerSkyModel(verbose=2)
])
# now add optional Jones terms
# these will show up in the menu automatically
# Ncorr - correct for N
from Siamese.OMS import oms_n_inverse
meqmaker.add_sky_jones('Ncorr', 'n-term correction', oms_n_inverse)
# Z - ionosphere
from Siamese.OMS import oms_ionosphere
meqmaker.add_sky_jones('Z', 'ionosphere', oms_ionosphere)
# L - dipole projection
from Siamese.OMS import oms_dipole_projection