Most information can be found in the cookbooks. The list below is only for tools specific for long baselines / for the Long Baseline Week.
This script will fit phase and delay offsets between superterp stations. It does not use the intra-superterp baselines, because they are generally dominated by RFI and interfering sources that would lead to incorrect offsets. Instead it uses baselines to (one or several) external stations and determines the superterp offsets using the closure phase relation. The user should select reference stations so that baselines to them have a nice and clean signal that is not affected by interfering sources (baselines not too short) but that is still strong enough (baselines not too short).
Please try the script lbw_fit_phases_superterp.py. In order to set the parameters, you have to edit the script itself, see the comments there for more details. Or ask for instructions…
This script needs owutil.py.
This script can add up visibilities to all superterp stations in order to create a virtual superterp after correlation. It also combines the flags (leaving the data unflagged only if all superterp stations are okay!) and the weights. For the weights it can take correlations between superterp stations into account, which generally reduces the weights.
The script can also correct for pre-determined phase (or delay) offsets between the superterp stations. The solution file for that (produced by the script discussed above and edited manually) has to be provided. An example file for calibration offsets can be found here. These numbers have been determined for L57169 (TauA) using DE602 as external station.
The file can include solutions for phase offsets alone or for phase+delay
offsets. If both exist, the latter will be used. If the files has several
solutions in it, the last one will be used! It is recommended to have a look
on the file instead of using the output from lbw_fit_phases_superterp.py
blindly.
Please try the script lbw_phase_up.py. In order to set the parameters, you have to edit the script itself, see the comments there for more details. Or ask for instructions…
This script needs owutil.py.
Tobia's converter can be found on cep1 at ~carozzi/scripts/mscorpol/, the documentation is in mscorpol_doc.pdf. You can also download a tar file here.
This converter uses a simple geometrical dipole beam model.
A very simple taql version without element beam model can be found here. Another script doing the same with plain pyrap was written by George Heald and is available on cep1 at ~carozzi/scripts/mslin2circ.py. Or download it here.
You can also try Olaf's lbw_dipole_correction.py instead. This is actually a multi-purpose script that can convert to circular with or without X-Y offset correction (phases only, no delays!), determine X-Y phase offsets, or only correct for X-Y offsets without conversion or …
As all of Olaf's scripts presented here, you have to set your paramters in the script. You cannot just call it with options.
This script needs owutil.py.
The beam model is the same as used in earlier work and also the same as applied by Tobia. The outputs of both converters seem to be consistent with each other, modulo an opposite sign convention for RL and LR.
Currently we can only fit and correct phase offsets, no delays yet.
Please try lbw_dipole_correction.py. This is actually a multi-purpose script that can convert to circular with or without X-Y offset correction (phases only, no delays!), determine X-Y phase offsets, or only correct for X-Y offsets without conversion or …
As all of Olaf's scripts presented here, you have to set your paramters in the script. You cannot just call it with options.
This script needs owutil.py.
The script fits for the X-Y offsets by minimising the leakage into RL and LR (which should be zero for unpolarised sources). In principle there is additional information in RR and LL, which should have the same amplitude but maybe different phases as result of DFR, but this information is not used yet.
Please also have a look at the example conversion to circular with X-Y offset determination.
ms2uvfits (casacore tool)
Parameters like the following seem to work:
ms2uvfits in=<whatever.ms> out=<whatever.fits> multisource=<False or True> writesyscal=False
difmap needs multisource=False (I think).
This script plots amplitudes as function of delay and time for all four polarisation products (in XY or RL). In addition it also plots phase combinations that can be used as estimates for X-Y offsets (if in XY) or DFR (if in RL).
You can download the script lbw_plot_pol_leakage.py here. In order to set the parameters, you have to edit the script itself, see the comments there for more details. Or ask for instructions…
This script needs owutil.py.