The new channel-averaged data base is often referred to as the channel 0 data base, and I will refer to the AVSPC output in this way henceforth. AVSPC will copy across the relevant extension files that are attached to the input file. It is wise not to include all channels in this average. The first and last few are likely to contain rubbish as the band response falls off fairly fast at the band edges. There is not much point to averaging in noise with your signal. The default for AVSPC is to take the central 75% of the band, and this should be adequate; experience with ATCA data indicates that a good channel range for 128 MHz data is 7 to 27. For other bandwidths, you must make your own assessment (see below), as the band edges roll off at different rates from bandwidth to bandwidth.
Setting POSSM in motion at this point will produce an amp-scalar averaged (in time) plot from all the visibilities. This means working out the amplitudes for each visibility and then averaging them, rather than averaging the real and imaginary parts of all the visibilities which is called a vector average. The phases are always computed with a vector average because of the circular nature of phase.
The scalar averaged spectrum contains contributions from emission in the whole primary beam. If the source is resolved, then different baselines will have different amplitudes, and the scalar average will relect that. Note that channels which contain only noise will appear biased after the scalar average at the mean noise level of the individual visibilities. This is because you are averaging a positive definite quantity; the amplitude of a noise signal. Any interference will be obvious in this plot.
On the other hand, the vector average aparm(1)=1 can be very useful for supressing interference. If the phase stability was good and as is usual, you did the basic on-line calibration with a strong calibrator (on-line programs CACAL or DELCOR), then the phases on all baselines would have been set to zero and the vector sum would be meaningful. Any interference generally has semi-random phase so it tends to average out. The vector average will make a spectrum at one location on the sky (it is basically a spectral-line cube with one spatial pixel) so make sure you set offset to point at the piece of sky you are interested in. Alternatively you could restrict the average to the short baselines so that the exact values of offset are less critical.
You can display this plot directly on the TV (set dotv=1 in POSSM), or display the resultant plot file (created if dotv=-1) with TKPL (plot on Tektronix or SUN Tektronix emulator), MIXPL (display on terminals or laser printers), or TVPL (display on a TV), and then you should have a fairly reasonable idea of which channels to keep and which channels to throw away. Note also that displaying on the SUN's Tektronix emulator with TKPL is much faster than displaying on the SUN's TV with TVPL or directly on the TV with dotv =1. However, for small plots, such as those produced by POSSM, it makes little difference, most of the time going into starting the task up rather than executing it.
The `divide by channel 0' (POSSM will just use the central 75% of the band for this) option (bparm(1)=1) might also prove handy. This option divides the frequency dependent visibilities by the channel 0 visibilities. This removes most of the calibration errors as well as most of the source structure and turns the source into a pseudo point source (you could then use something other than a point source for these plots). The remaining response across the band is because of the filters and electronics. Experiment with this option, and see how much it affects the outcome. Make sure there is no severely corrupted data in the channel 0 data channels if you use this option.
Note that channel 1 in ATCA spectra contains the zero frequency or DC term. This is why there are 33 channels in 128 MHz bandwidth spectra, rather than 32 as you would expect. This is very useful in determining if your data contain DC offsets, a problem discussed more below.
POSSM | |
source='1934-638',' ' | Select source by name or perhaps |
calcode='p' | select the primary with its calcode |
if you set it in SETJY | |
freqid=1 | select desired frequency |
antennas=0 | average over all antennas |
docalib=-1 | no calibration to apply |
aparm=0 | scalar average |
aparm(1)=1 | vector average |
bparm=0 | |
solint=0 | average all data |
solint=-1 | scan averages |
solint=# | average data over this many minutes |
dotv=1 | plot directly on TV or |
dotv=-1 | make a plot file for display later |
MIXPL | |
outfile=' ' | Don't write plot to disk file |
invers=0 | Plot highest version plot file |
device=0 | Laser printer |
device=10 | Visual 603 terminal |
If the lag spectrum is corrupted by a DC offset, its Fourier transform is modifed by the addition of a zero frequency delta function. The convolution of this delta function (which might be very strong) and the SinC squared function can thus produce a strong ringing (i.e., propogation of the SinC squared sidelobes) throughout the spectrum. You can see any such DC offset in the spectrum that POSSM produced (see above) by examining channel 1 of the spectrum (see above). If the offset just occurs in one baseline you would need to plot individual baselines with POSSM to see this clearly.
However, there is an easy cure. The SinC squared function has zeros every other odd channel. Therefore, you need only select the odd numbered channels in the good channel range that you determined above with POSSM. In addition, because adjacent channels are not independent, you sacrifice very little with regards to the signal-to-noise ratio in doing this. AVSPC has an easy option for doing just this. It may be worth doing regardless of whether you can see any DC term or not. However, note that this discussion is for the 128 MHz bandwidth only. The situation is more complex at other bandwidths and this simple procedure is not a cure for them.
AVSPC | |
flagver=0 | Apply highest FG table if any |
chansel=7,27,2 | Select odd channels in range 7 to 27 |
chansel=7,27,1 | Select all channels in range 7 to 27 |