예제 #1
0
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
예제 #2
0
                       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())

import Purr.Pipe