Tools and scripts

Cookbooks

Most information can be found in the cookbooks. The list below is only for tools specific for long baselines / for the Long Baseline Week.

Fitting offsets between superterp stations

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.

Phasing up the superterp post-correlation

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.

Conversion to circular

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.

Fitting and correcting X-Y offsets

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.

Convert to uvfits for AIPS and difmap

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).

Plotting tools

Delay and polarisation leakage plots

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.

 
lbg/lbw1/tools.txt · Last modified: 22-May-2012 18:44 by Olaf Wucknitz
[unknown button type]
 
Recent changes RSS feed Donate Powered by PHP Valid XHTML 1.0 Valid CSS Driven by DokuWiki