(Single) station issues

The behaviour of the individual stations is relevant for all baselines, not only international ones. Because long baselines require the maximum achievable S/N, it is of particular importance to optimise the stations' performance. This section of the wiki covers station issues, including single-station observations.

Effelsberg HBA sensitivity problem

Olaf Wucknitz 31-Aug-2012 11:19

It has been observed for quite a while that the Effelsberg station DE601 often has a rather low sensitivity, particularly in the HBA. Bad sensitivity in general can be explained by bad station calibration, so it was not clear if this is a particular problem of Effelsberg. However, theory predicts that the sensitivity in the nominal phase centre of the beam should not depend on the position on the sky, no matter how good or bad the calibration is. The only position-dependent effect in the beam centre should be due to the element beams, which are not only relatively well behaved but also have certain symmetry properties. In contrast to this expectation, several observations shows serious variations with sky position.

Interferometry

In long-baseline interferometric observations, we have noticed earlier that the sensitivity is not as symmetric around the zenith as expected. This can be seen in AIPS calibration solution plots of L45786, a rather short observation. Different colours are for different frequency bands within the HBA low range. Page 8 shows the amplitude solutions of DE601 in comparison with other stations. Remember that high values mean low sensitivity. The Effelsberg gains vary over short durations much more than expected.

Another and better example is L35012, a 12h long track observation. Page 9 shows how the sensitivity of DE601 drops to almost zero around 23h-00h in a frequency-dependent way. No other station shows a similar behaviour.

Single-station pulsar observations

Masaya Kuniyoshi performed a number of tests observing the bright pulsar B0329+54 and comparing the sensitivity for different parts of its track around the sky. He found that the sensitivity shows strong variations (by more than an order of magnitude), and that these variations are not symmetric around zenith. Depending on the version of the calibration table, the sensitivity was best either east or west of the zenith.

Mapping the beam with single-station observations

I thought about several ways to map out the beam shape of the station in order obtain more information about its behaviour. Just the amplitude in the centre does not tell us much, not even if we should blame the analog tile beam or the digital station beam. In the end I decided to do single-station observations using CygA and CasA as bright test sources. I mapped the beam by forming >200 beamlets in a grid around either of these sources. At any given time, all beamlets used the same subband (to reduce the influence of RFI etc), but the subbands were switched between a selection of relatively clean frequencies. Only the beamlet statistics were recorded, but no voltages. As additional diagnostic, I also stored local station correlations and subband statistics. In the night 21-22 Aug 2012 I did this for Effelsberg using the subbands 75,310,420 (114.6,160.5,182.0 MHz).

Beam maps in different coordinate systems are collected here:

sb 075 sb 310 sb 420
CygA radec lm radec lm radec lm
CasA radec lm radec lm radec lm

Each PDF file has one page per snapshot (30sec). Here 'radec' refers to relative equatorial coordinates. In this system the beam rotates and gets distorted due to projection effects. To get rid of these effects, I also present projections to the horizontal components of the local station system, labelled as 'lm'. The latter system has the advantage that the elements are horizontally aligned, which results in translation-independent beams and an automatic compensation for any projection effects. We thus expect that in this system the beams should be stationary.

However, this is not what we find. The beams moves around significantly, by about 1-2 degrees.

Because the resolution of the beamlet plots is limited, I also include (lower two panels) images derived from the local correlations after applying the calibration table. These maps are entirely consistent with the upper panels, which confirms my correct understanding of coordinates and phase sign conventions. Having this option available also allows us to redo the maps with different calibration parameters or element position files without actually repeating the observations. This diagnostic turned out to be very helpful.

As comparison I produced similar images for DE603 (Tautenburg) and DE605 (Jülich), derived from observations of 27-28 Aug 2012. Many thanks to science support at ASTRON and to Eva Juette (Jülich) and Matthias Hoeft (Tautenburg) for making this possible! This was a great help.

Here is the corresponding table for DE603:

sb 075 sb 310 sb 420
CygA radec lm radec lm radec lm
CasA radec lm radec lm radec lm

… and for DE605:

sb 075 sb 210 sb 270 sb 310 sb 375 sb 420
CygA radec lm radec lm radec lm radec lm radec lm radec lm
CasA radec lm radec lm radec lm radec lm radec lm radec lm

Note that the observing times are not all the same.

At these stations, no funny effects are observed, besides the occasional RFI producing additional peaks. The beams look much nicer, and they are solid as a rock. This is how it should be!

Incorrect element coordinates

What could cause the motion of the beam relative to its nominal position? A bug in the beam former is possible, but not very likely given the relatively small effect. Another option are incorrect element coordinates. However, they cannot be entirely wrong, because that would distort the beam much more. The gentle shift does instead point towards minor problems with the local coordinate system.

If the coordinate system is rotated a bit relative to reality, we can describe the observed shift as follows: As a result of the rotation of the sky (or the Earth for the sceptics amoung you), the source moves around in the local (real) lmn coordinate system in a predictable way. The beamformer compensates for this, but it may use its own (incorrect) lmn system. Because the rotation of the sky is the same in both cases, the transformation from real to incorrect coordinates is simply the small rotation between the systems. What we observe is the apparent shift, which is the difference between the real and incorrect coordinates.

This means the shift is given by the rotation angle multiplied by the generating matrix of the rotation. It consists of three free parameters that can be fitted analytically in the first-order approximation of small angles. The procedure basically consists of fitting a linear (including constant offset) model for the shift as function of lmn coordinates. The best fit is consistent for both sources and all observed frequencies and corresponds to no rotation around the l and m axes but a rotation of 1.5 to 1.6 degrees around the (vertical) n axis. In other words the station coordinates are rotated by 1.55 degrees in the horizontal plane.

Test with corrected coordinates

Instead of repeating the beamlet observations with modified element coordinates, I reused the local correlations but used the rotated coordinates for the imaging. Results in the lm system are given here:

sb 075 sb 310 sb 420
CygA lm lm lm
CasA lm lm lm

Note that the upper two panels are not changed. They are just included as comparison.

We find that the rotated coordinates seem to solve the problem: The beams do not move around anymore but are stationary as they should. The beam shape is still bad and offset from the centre, because the station calibration has not been redone (yet).

Calibration tables

I can also apply different (or no) calibration tables before imaging the correlations. Maps are available and can be added to the wiki on request.

Effect on sensitivity

Here are plots on the sensitivity at the beam peak and nominal beam centre for the original data:

sb 075 sb 310 sb 420
CygA orig rotated orig rotated orig rotated
CasA orig rotated orig rotated orig rotated

'orig' is for the original version, 'rotated' has the coordinates rotated for the imaging results. Each file shows as function of time: min/max of all beamlets, value in the nominal centre, and the min/max(2) from the imaging data, and the centre(2) of the imaging. We generally see that the centre of the beam is much lower than the max and that the value at the centre shows the observed variations as function of time. The maximum itself, even though it may be lower than for perfect calibration, shows about the expected variation with elevation and azimuth. Elevation calculated for the local station coordinate system is shown for comparison.

For 'rotated', only the centre2 curves change, because only those beam images have been corrected. We see that the level as such does not always get better (because the station calibration is still bad), but the changes with time now follow the behaviour of the maximum that also resembles the expectations from the element response.

Curves for the other stations can be added on request.

Conclusions

The coordinates of the Effelsberg HBA tiles are incorrect. Rotating them by 1.5 to 1.6 degrees around the vertical axis seems to fix the problem. Some questions remain:

  • how did this happen?
  • when did it happen?
  • should we also rotate the local coordinate system or only the element positions?
  • what about the LBA positions?

We have 'new' measurements of the element positions available and should put them in the AntennaField.conf file as soon as possible. After that the station calibration should be repeated. For this we can even reuse the previous observations.

Olaf Wucknitz 31-Aug-2012 12:54

Addendum

The original determination of the rotation angle used the beginning of each 30 sec snapshot to calculate the appropriate delays, not the middle. Correcting for this increases the angle slightly, to about 1.6 deg.

Here are results from individual fits:

PLOTS_temp/DE601HBA_casa_sb075_peak.txt
ipol 0  peak2  read 31 good   rotations: -0.155159 +0.016478 +1.567199  [deg]    mean err 0.054048   no fit: 2.347833   max 0.163467 4.107604  [deg]
ipol 1  peak2  read 31 good   rotations: -0.144448 +0.114021 +1.566374  [deg]    mean err 0.051291   no fit: 2.591965   max 0.141686 4.469822  [deg]
PLOTS_temp/DE601HBA_casa_sb310_peak.txt
ipol 0  peak2  read 31 good   rotations: -0.068465 -0.054442 +1.652892  [deg]    mean err 0.035403   no fit: 1.047701   max 0.086019 2.235176  [deg]
ipol 1  peak2  read 31 good   rotations: -0.146425 -0.121926 +1.678144  [deg]    mean err 0.045300   no fit: 0.954634   max 0.128048 2.106454  [deg]
PLOTS_temp/DE601HBA_casa_sb420_peak.txt
ipol 0  peak2  read 31 good   rotations: -0.019954 -0.012736 +1.663448  [deg]    mean err 0.038866   no fit: 0.872822   max 0.110787 1.956218  [deg]
ipol 1  peak2  read 31 good   rotations: -0.037443 -0.086530 +1.681440  [deg]    mean err 0.040619   no fit: 0.817438   max 0.100355 1.768338  [deg]
PLOTS_temp/DE601HBA_cyga_sb075_peak.txt
ipol 0  peak2  read 24 good   rotations: -9.768076 +5.951909 -1.678177  [deg]    mean err 3.580895   no fit: 4.181896   max 8.600379 9.021761  [deg]
ipol 1  peak2  read 24 good   rotations: -1.136360 -0.468115 +1.510575  [deg]    mean err 0.058256   no fit: 3.033756   max 0.255960 5.529044  [deg]
PLOTS_temp/DE601HBA_cyga_sb310_peak.txt
ipol 0  peak2  read 24 good   rotations: -0.286233 -0.544840 +1.638029  [deg]    mean err 0.039315   no fit: 1.454305   max 0.203168 2.803888  [deg]
ipol 1  peak2  read 24 good   rotations: -0.072902 -0.481818 +1.679534  [deg]    mean err 0.037520   no fit: 1.267094   max 0.159909 2.484623  [deg]
PLOTS_temp/DE601HBA_cyga_sb420_peak.txt
ipol 0  peak2  read 23 good   rotations: -0.118661 -0.380058 +1.645125  [deg]    mean err 0.016417   no fit: 1.227655   max 0.048035 2.375695  [deg]
ipol 1  peak2  read 23 good   rotations: -0.236542 -0.465762 +1.651933  [deg]    mean err 0.017252   no fit: 0.964057   max 0.059468 2.115911  [deg]

The results are the rotations values. We see that most of them agree quite well, with the exception of CygA SB075 for ipol=0. But here the residuals after the fit (mean position error of the beam 3.58 deg) is also much larger than for all other fits. Looking at the movies of the beam shapes shows that some frames are affected by RFI or sidelobes. We can conclude that all other fits show a rotation around the vertical axis of 1.5 to 1.7 deg. 1.6 seems to be a good average. Rotations around the other axes are negligible.

I also ran similar fits for a 24h LBA observations. The measurement is much less clean, also due to strong intermittent RFI, but the results seems to be consistent with no rotation at all.

Olaf Wucknitz 03-Sep-2012 15:21

Today James Anderson took recent (May 2012) measurements of the element positions and produced a table that can be compare with AntennaField.conf. It turns out that the new coordinates are indeed rotated relative to the old ones by 1.6 degrees. The deviations from this rotation matrix are less than 0.05 deg.

The LBA positions are rotated by max. 0.1 deg.

Olaf Wucknitz 03-Sep-2012 21:42

Effelsberg cable lengths (HBA)

Initial analysis

In station calibration tests with local correlations (CasA, HBA low, but other data also available) taken 5/6 Sep 2012, I noticed that some HBA tiles have delays that are very different from the majority. Most are very close to zero with a scatter of ca. 2-3 nsec, but there are delays up to +-125 nsec, corresponding to almost 40 m in vacuum.

Here are numbers from some fits. They do not have the highest possible accuracy, but are good enough for this purpose:

el  3  ipol 0  phase   19.868 deg  delay   64.150 nsec
el  4  ipol 0  phase    7.273 deg  delay   63.666 nsec
el  5  ipol 0  phase   70.305 deg  delay  124.354 nsec
el 11  ipol 0  phase   22.098 deg  delay   63.429 nsec
el 20  ipol 0  phase   30.057 deg  delay   64.047 nsec
el 37  ipol 0  phase -109.703 deg  delay  -51.903 nsec
el 38  ipol 0  phase  -90.417 deg  delay  -51.864 nsec
el 46  ipol 0  phase -100.576 deg  delay  -52.030 nsec
el 47  ipol 0  phase  -82.251 deg  delay  -51.493 nsec
el 82  ipol 0  phase   53.040 deg  delay -126.612 nsec

el  3  ipol 1  phase    6.786 deg  delay   63.891 nsec
el  4  ipol 1  phase    7.176 deg  delay   64.081 nsec
el  5  ipol 1  phase   75.115 deg  delay  122.093 nsec
el 11  ipol 1  phase    7.172 deg  delay   63.773 nsec
el 20  ipol 1  phase   21.977 deg  delay   64.309 nsec
el 37  ipol 1  phase -111.961 deg  delay  -51.549 nsec
el 38  ipol 1  phase -112.679 deg  delay  -52.373 nsec
el 46  ipol 1  phase -102.704 deg  delay  -52.253 nsec
el 47  ipol 1  phase -105.893 deg  delay  -50.844 nsec
el 82  ipol 1  phase   69.122 deg  delay -125.398 nsec

Both polarisation are affected in the same way. All others are well below 10 nsec. Because these offsets are so close to the cable length differences, and because there are equal numbers of positive and negative offsets for 50-64 nsec and 120-126 nsec, a plausible cause could be swapped cable lengths in the CableDelays.conf file. Note that the RCU numbers in that file are 2*(element-number)+ipol.

The current calibration table corrects for some of these delays, but not for all. This is not too surprising, because offsets off this magnitude are not really expected, and the calibration software may not be suited for them.

This is from the CableDelays.conf file:

# Lenghts are in meters, delays are in ns.
#
# Note: The first order values are:
#	50m		199.2573
#	80m		326.9640
#	85m		342.5133
#	115m	465.5254
#	130m	530.6981

Estimated from these numbers, my 50-64 and 120-126 nsec offsets correspond to 15 m and 30 m cable length, respectively. The only way to match the corrected cable lengths with the available lengths of 50/80/85/115/130m is to apply the corrections to these numbers with negative sign.

Here are the values from the file for the affected elements together with the corrected cable lengths:

                     HBA
#RCUnr	        len	delay           corrected
#----------------------------------------------------
6		130     530.6981           115
7		130     530.6981           115
8		130     530.6981           115
9		130     530.6981           115
10		115     465.4000            85
11		115     465.6000            85
22		130     530.6981           115
23		130     530.6981           115
40		130     530.6981           115
41		130     530.6981           115
74		115     465.5254           130
75		115     465.5254           130
76		115     465.5254           130
77		115     465.5254           130
92		115     465.5000           130
93		115     465.5000           130
94		115     465.5000           130
95		115     465.5000           130
164		 85     342.5133           115
165		 85     342.5133           115

Or in summary:

115m instead of 130m: RCUs   6,  7,  8,  9, 22, 23, 40, 41  antennas  3, 4,11,20
130m instead of 115m: RCUs  74, 75, 76, 77, 92, 93, 94, 95  antennas 37,38,46,47
 85m instead of 115m: RCUs  10, 11                          antennas  5
115m instead of  85m: RCUs 164,165                          antennas 82

This does indeed look as if four antennas were swapped between 115 and 130m and one antenna between 85 and 115m.

Note that this still has to be confirmed.

Olaf Wucknitz 13-Sep-2012 14:22

Old CableDelays.conf files

An old version of the cable length file de601-cabledelays.conf-rev14249.svn000.tmp.conf actually seems to be more accurate than the current version. It is unclear why the files has been changed.

This old version and the differences have been found by Menno Norden. Thanks for that!

Here is a diff output of the two files:

diff DE601-CableDelays.conf-rev14249.svn000.tmp.conf DE601-CableDelays.conf
1a2,4
> # CableDelays.conf for DE601
> #
> #
33,38c36,41
< 6             130     530.6981        130     530.6981        115     465.5000
< 7             130     530.6981        130     530.6981        115     465.5000
< 8             130     530.6981        130     530.6981        115     465.5000
< 9             130     530.6981        130     530.6981        115     465.4000
< 10            130     530.6981        130     530.6981         85     342.4000
< 11            130     530.6981        130     530.6981         85     342.6000
---
> 6             130     530.6981        130     530.6981        130     530.6981
> 7             130     530.6981        130     530.6981        130     530.6981
> 8             130     530.6981        130     530.6981        130     530.6981
> 9             130     530.6981        130     530.6981        130     530.6981
> 10            130     530.6981        130     530.6981        115     465.4000
> 11            130     530.6981        130     530.6981        115     465.6000
49,50c52,53
< 22            130     530.6981        130     530.6981        115     465.6000
< 23            130     530.6981        130     530.6981        115     465.4000
---
> 22            130     530.6981        130     530.6981        130     530.6981
> 23            130     530.6981        130     530.6981        130     530.6981
67,68c70,71
< 40            130     530.6981        130     530.6981        115     465.5000
< 41            130     530.6981        130     530.6981        115     465.5000
---
> 40            130     530.6981        130     530.6981        130     530.6981
> 41            130     530.6981        130     530.6981        130     530.6981
85,86c88,89
< 58            130     530.6981        130     530.6981        130     519.7000
< 59            130     530.6981        130     530.6981        130     519.7000
---
> 58            130     530.6981        130     530.6981        115     465.5254
> 59            130     530.6981        130     530.6981        115     465.5254
101,104c104,107
< 74            130     530.6981        130     530.6981        130     519.8000
< 75            130     530.6981        130     530.6981        130     519.7000
< 76            130     530.6981        130     530.6981        130     519.6000
< 77            130     530.6981        130     530.6981        130     519.8000
---
> 74            130     530.6981        130     530.6981        115     465.5254
> 75            130     530.6981        130     530.6981        115     465.5254
> 76            130     530.6981        130     530.6981        115     465.5254
> 77            130     530.6981        130     530.6981        115     465.5254
113,114c116,117
< 86            130     530.6981        130     530.6981         85     342.5000
< 87            130     530.6981        130     530.6981         85     342.5000
---
> 86            130     530.6981        130     530.6981        115     465.5254
> 87            130     530.6981        130     530.6981        115     465.5254
191,192c194,195
< 164           130     530.6981        130     530.6981        115     465.5000
< 165           130     530.6981        130     530.6981        115     465.4000
---
> 164           130     530.6981        130     530.6981         85     342.5133
> 165           130     530.6981        130     530.6981         85     342.5133

It is easy to confuse things now, but I think this is the correct summary:

6-11,22-23,40-41,74-77,164-165 have been correct in rev14249 and are wrong now
58-59,86-87 have been wrong in rev14249 and are correct now
92-95 have been wrong and are still wrong

This seems to be consistent with Menno's analysis. In other words neither of the files is correct and we should create a new version, but maybe only after confirming the conclusions at the station.

Olaf Wucknitz 14-Sep-2012 09:47

Swapped tiles?

It cannot be completely ruled out that some of the tiles are connected incorrectly, which would also produce incorrect cable lengths per rcu number. I only had a quick look at this possibility. In such a case we would expect the delays to change with sky position, because the geometric delays would not be taken care of entirely. For the following plot I first produced a calibration table for a certain time (integration 16 in this case), then calibrated with this table and determines residual delays relative to it as function of time and rcu number:

residual delays

This is only for polarisation Y (X has more RFI). The elevation varies between about 80 deg and something like 50 deg or so. If there are any systematic delay changes, they should not be more than ca. 4 nsec, corresponding to 1.2 m. Even with the projection effects, no tile seems to be off by more than 2.5 m. Swapped neighbouring tiles are still a (unlikely) possibility, but swaps of more separated tiles don't seem to exist.

Olaf Wucknitz 14-Sep-2012 12:16

Is that it?

Maybe not. James mentions that he did not find any large delays in HBA observations he took with the old version of the table. This would be evidence against incorrect cable lengths for 58-59,86-87,92-95 in rev14249.

Menno Norden provided a new file with corrections applied. There are also comments in the file explaining what has been changed:

de601-cabledelays.conf-modified-20120914-norden.conf

Olaf Wucknitz 14-Sep-2012 17:46

Check of cable lengths at station

On 16th October Menno Norden, Henri Meulman, Ian Stewart, Masaya Kuniyoshi, Andreas Horneffer and Olaf Wucknitz were at the station and (amongst other things) check the cable lengths in question. Reading the cable labels in the container is virtually impossible, so we did that at the tiles:

RCU antenna current length correct length according to analysis read at tile
58-59 29 130 115 115
86-87 43 85 115 115
92-93 46 115 130 130
94-95 47 115 130 130

This confirms the analysis but is not consistent with no large delays in James' earlier experiments. A new version of the table has been installed now: de601-cabledelays.conf-modified-20121016-norden.conf

A diff of old versus new (without the comment lines) produces the following:

7,12c7,12
< 6		130	530.6981	130	530.6981	130     530.6981
< 7		130	530.6981	130	530.6981	130     530.6981
< 8		130	530.6981	130	530.6981	130     530.6981
< 9		130	530.6981	130	530.6981	130     530.6981
< 10		130	530.6981	130	530.6981	115     465.4000
< 11		130	530.6981	130	530.6981	115     465.6000
---
> 6		130	530.6981	130	530.6981	115     465.5000
> 7		130	530.6981	130	530.6981	115     465.5000
> 8		130	530.6981	130	530.6981	115     465.5000
> 9		130	530.6981	130	530.6981	115     465.4000
> 10		130	530.6981	130	530.6981	 85     342.4000
> 11		130	530.6981	130	530.6981	 85     342.6000
23,24c23,24
< 22		130	530.6981	130	530.6981	130     530.6981
< 23		130	530.6981	130	530.6981	130     530.6981
---
> 22		130	530.6981	130	530.6981	115     465.6000
> 23		130	530.6981	130	530.6981	115     465.4000
41,42c41,42
< 40		130	530.6981	130	530.6981	130     530.6981
< 41		130	530.6981	130	530.6981	130     530.6981
---
> 40		130	530.6981	130	530.6981	115     465.5000
> 41		130	530.6981	130	530.6981	115     465.5000
75,78c75,78
< 74		130	530.6981	130	530.6981	115     465.5254
< 75		130	530.6981	130	530.6981	115     465.5254
< 76		130	530.6981	130	530.6981	115     465.5254
< 77		130	530.6981	130	530.6981	115     465.5254
---
> 74		130	530.6981	130	530.6981	130     519.8000
> 75		130	530.6981	130	530.6981	130     519.7000
> 76		130	530.6981	130	530.6981	130     519.6000
> 77		130	530.6981	130	530.6981	130     519.8000
93,96c93,96
< 92		130	530.6981	130	530.6981	115     465.5000
< 93		130	530.6981	130	530.6981	115     465.5000
< 94		130	530.6981	130	530.6981	115     465.5000
< 95		130	530.6981	130	530.6981	115     465.5000
---
> 92		130	530.6981	130	530.6981	130     519.6000
> 93		130	530.6981	130	530.6981	130     519.6000
> 94		130	530.6981	130	530.6981	130     519.6000
> 95		130	530.6981	130	530.6981	130     519.6000
165,166c165,166
< 164		130	530.6981	130	530.6981	 85     342.5133
< 165		130	530.6981	130	530.6981	 85     342.5133
---
> 164		130	530.6981	130	530.6981	115     465.5000
> 165		130	530.6981	130	530.6981	115     465.4000

This is entirely consistent with the differences found in my analysis earlier. The CableDelays.conf file should be correct now.

Olaf Wucknitz 16-Oct-2012 15:32

Effelsberg post-maintenance tests

One week after the maintenance visit by Menno Norden and Henri Meulman (plus locals), I ran some test observations to check the cable length files and other aspects.

No X-Y swaps in HBA

I ran an additional test to check for any X/Y swaps of HBA elements, using astronomical data.

I first found a relatively clean subband (205 in this case). Then I modified the BeamServer.conf to update the analog beamformer only once an hour (and thus not at all in this short test), to allow to first form a beam and then switch off some elements.

I then set the analog beamformer to a bright source (CasA in this case). I select a number of reference tiles (0,1,2 in one case and 80,81,82 in the second). For these reference tiles I keep all elements (to have a better S/N). For the other tiles I switch off all but one element (looping through all 16). Then I record cross-correlations (rspctl –xcstatistics) for 1min each with 1sec integrations. After the observation I check for each element number if the parallel-hand correlations to the reference tiles are stronger or weaker than the cross-hands. For each tile/element to test I print the fraction of integrations/reference tiles (out of 60×3) for which the parallel-hand is stronger than the cross-hand. This fraction should be almost 1 if everything is okay, and this is indeed the case for all tiles and all elements. The fractions are all > > 0.9, which means that no X/Y dipoles are swapped.

Of course one could use a better statistics than just counting “good” ratios of parallal/cross, but for this case the simple approach seems to be sufficient.

Remaining LBA problems

I still find some problems with a couple of LBA:

rcu 39 sometimes high
rcu 68 low (but stable)
rcu 74 low (but stable)
rcu 119 often high

HBA element positions and calibration

The following figure shows images and phases in the direction of CasA using 1 sec of data per subband. I initially included subbands 50-462 but then (automatically) flagged those with strong RFI. For the calibration I fitted phases and delays for all HBA tiles. The plotted images and phases are averaged over the good subbands. The angular scale is in units of rad for lambda=2m. Subbands are averaged rescaling the angular scale to have the same resolution. In this way all subbands should have the same station beam. The coordinates are relative 2d station coordinates (called lm before). The numbers above the images are diagnostics to judge the calibration quality. centre/peak/coher stand for the value in the centre, the value at the peak, and the value in the centre divided by the sum of visibility amplitudes. The latter would be unity for perfect calibration and a signal dominated by a point-source in the centre. Applying the calibration increases this from about 0.1 to almost 1.

Note that the colour scale is different for the calibrated phases. Calibration errors would show up as horizontal and vertical bars with systematic deviations, where horizontal and vertical have different signs. We see on the right that these errors disappear after calibration.

comparison of uncalibrated and calibrated data (click to enlarge)

I also produced a movie (also available as PDF) that spans 10 h and applies an independent calibration for each of the 68 snapshots. Particularly interesting is the uncalibrated left-hand side. Even though the beam does not look great, it now stays constant and does not move around with time. This tells us that the problem with the rotated element positions is ironed out eventually.

The resulting calibration tables are pretty constant so that (as expected) one calibration should be good for the entire observation. I applied the calibration from the snapshot shown above over the full run, which results in a calibrated movie (also as PDF). Not surprisingly the calibration gets worse near the horizon.

I applied the same calibration also to observations of CygA, resulting in a movie and a corresponding PDF file. The coherence is only slightly worse than for the CasA data themselves.

The delays from the fits are mostly very small, around 1-2 nsec. The largest one is 3.4 nsec, corresponding to less than 1m of cable delay. This means that all entries in CableDelays.conf are now correct, which indirectly implies that both of the previously installed versions were incorrect. I cannot tell why the large delays have not shown up much earlier. Recall that we had cable length errors of up to 40m before!

Carmen Toribio produced a calibration table from observations taken a few days earlier. I went through the same exercise as above but initially applied that calibration table (which is now installed on the LCU).

Here is the snapshot:

official and own calibration (click to enlarge)

We see that most antennas already have the correct phases. The additional delays needed now are all very small (< 1 nsec) with the exception of the following:

 7  ipol 0  phase  -85.120 deg  delay    4.742 nsec
19  ipol 0  phase  -82.224 deg  delay    4.659 nsec
23  ipol 0  phase  -87.879 deg  delay    4.680 nsec
26  ipol 0  phase  -85.406 deg  delay    4.620 nsec
27  ipol 0  phase  -83.324 deg  delay    4.664 nsec
31  ipol 0  phase  -86.116 deg  delay    4.648 nsec
 7  ipol 1  phase  -84.459 deg  delay    4.559 nsec
19  ipol 1  phase  -83.685 deg  delay    4.800 nsec
23  ipol 1  phase  -87.458 deg  delay    4.651 nsec
26  ipol 1  phase  -84.579 deg  delay    4.704 nsec
27  ipol 1  phase  -83.671 deg  delay    4.799 nsec
31  ipol 1  phase  -86.001 deg  delay    4.685 nsec

All other delays shifted to about -0.3 nsec, because the mean is fixed to 0. This means the additional delays are almost exactly 5.0 nsec. It is thus plausible that they are resulting from the infamous 5 nsec problem that will be fixed with a hardware upgrade. In this particular occasion we are losing about 15% sensitivity. Note perfect, but much better than what we had before!

Olaf Wucknitz 29-Oct-2012 18:48

New clock and PPS synchronisation

On 15th November 2012, Andreas Horneffer (with help from Masaya Kuniyoshi and Olaf Wucknitz watching) installed the new SyncOptics boards that hopefully improve the relative timing and solve the 5 nsec problem.

Click here for details of performance tests done later.

 
lbg/single/start.txt · Last modified: 15-Mar-2013 12:03 by Olaf Wucknitz
[unknown button type]
 
Recent changes RSS feed Donate Powered by PHP Valid XHTML 1.0 Valid CSS Driven by DokuWiki