next up previous contents
Next: Consequences Up: class-weighting Previous: Fixing folded frequency switching   Contents

Fixing the AVERAGE output system temperature

The commands AVERAGE, ACCUMULATE, and STITCH combine the spectra available in the current index according to their relative weights, and produce a single output spectrum with a proper description in terms of observing time, resolution, and system temperature.

Assuming the simplest case of a set of $N$ spectra with identical observing time, resolution, and system temperature, the output spectrum has a resolution identical to the input one, and an observing time which is the sum of all the inputs one. The output system temperature is then derived by reverting the theoritical time weight formula

\begin{displaymath}
w_{\ensuremath{\mathrm{out}}} = \sum_N w_{\ensuremath{\mathrm{in}}}
\end{displaymath} (3)


\begin{displaymath}
T_{\ensuremath{\mathrm{sys,out}}} = w_T^{-1}(w_{\ensuremath...
...{\ensuremath{\mathrm{res,out}}},t_{\ensuremath{\mathrm{out}}})
\end{displaymath} (4)

When averaging folded fsw with TIME weighting, this gives:
\begin{displaymath}
w_{\ensuremath{\mathrm{T,out}}} =
\sum_N w_{\ensuremath{\...
...ensuremath{\mathrm{res}}} / T_{\ensuremath{\mathrm{sys,in}}}^2
\end{displaymath} (5)

If we reverse the generic $w_{\ensuremath{\mathrm{T}}}$ formula, this gives an unexpected change of system temperature:
\begin{displaymath}
T_{\ensuremath{\mathrm{sys,out}}} = \frac{T_{\ensuremath{\mathrm{sys,in}}}}{\sqrt{2}}
\end{displaymath} (6)

In other words, when the factor 2 was introduced in 2019, it was not added symetrically in the $w_T^{-1}$ reverse operation done by AVERAGE. The result was that when all the averaged spectra were folded fsw, the output Tsys was underestimated by a factor $\sqrt {2}$. As of 24-feb-2021 and release mar21, this is fixed now:

\begin{displaymath}
w_{\ensuremath{\mathrm{T,ffsw}}}^{-1} = \frac{1}{2} \times w_{\ensuremath{\mathrm{T}}}^{-1}
\end{displaymath} (7)

This reverse function is used when mixing folded fsw spectra only. When mixing psw or wsw spectra only, the generic $w_{T}^{-1}$ function is used. Finally, when mixing folded fsw with psw or wsw, the output Tsys is more complicated to evaluate, as there is no unique weight formula to revert. The chosen solution is to use also the generic $w_{\ensuremath{\mathrm{T}}}^{-1}$ function, leading to an approximate evaluation of the output system temperature. In this latter case, we also introduce a mix switching mode (reflecting the combination of spectra observed in different modes), and the associated section in the spectrum header is emptied except for this code.

For example, when mixing a one psw or wsw spectrum with one fsw spectrum of same system temperature, integration time, and frequency resolution, this gives:

$\displaystyle w_{\ensuremath{\mathrm{T,psw}}}$ $\textstyle =$ $\displaystyle t_{\ensuremath{\mathrm{in}}} \times f_{\ensuremath{\mathrm{res}}} / T_{\ensuremath{\mathrm{sys,in}}}^2$ (8)
$\displaystyle w_{\ensuremath{\mathrm{T,ffsw}}}$ $\textstyle =$ $\displaystyle 2 \times t_{\ensuremath{\mathrm{in}}} \times f_{\ensuremath{\mathrm{res}}} / T_{\ensuremath{\mathrm{sys,in}}}^2$ (9)
$\displaystyle w_{\ensuremath{\mathrm{T,mix}}}$ $\textstyle =$ $\displaystyle t_{\ensuremath{\mathrm{out}}} \times f_{\ensuremath{\mathrm{res}}} / T_{\ensuremath{\mathrm{sys,mix}}}^2$ (10)

with $t_{\ensuremath{\mathrm{out}}} = 2 \times t_{\ensuremath{\mathrm{in}}}$ (sum of input integration times). Since $w_{\ensuremath{\mathrm{T,mix}}} = w_{\ensuremath{\mathrm{T,psw}}}+w_{\ensuremath{\mathrm{T,ffsw}}}$, this yields:
\begin{displaymath}
T_{\ensuremath{\mathrm{sys,mix}}} = \sqrt{\frac{t_{\ensurem...
... =
T_{\ensuremath{\mathrm{sys,in}}} \times \sqrt{\frac{2}{3}}
\end{displaymath} (11)

We see that the resulting system temperature is underestimated because the reverse weight formula is not ideal in the mixed case.


next up previous contents
Next: Consequences Up: class-weighting Previous: Fixing folded frequency switching   Contents
Gildas manager 2023-06-01