Nuclear Instuments and Methods, Vol. 161 pp. 137-140
On Digital Filtration in Correlation Time-of-Flight Spectrometry
Laszlo J. Naszodi
Abstract: A new method of filtering the correlation time-of-flight spectrum is proposed, which fits well to the direct evaluation of the spectrum. It makes additional mechanical and electronical devices unnecessary and the loss in intensity is less than that of the usual filtering technique.
1. Introduction
In the correlation technique a pseudo-random binary sequence (PRBS) is used for chopping the incident beam 1, 2). As its definition, a PRBS {aj } containing m ones and N-m zeros has the autocorrelation function
Caa(v) = Σ aj aj+v = | { | m if v=0 (mod N) h = mc if v≠0 (mod N) |
(1) |
where c=(m-1)/(N-1) is called the duty cycle of the PRBS. This technique yields the time-of-flight (TOF) spectrum with the counts
Zk = Σ aki Si + bk , k=0,1,2,…,N-1, (2)
where
aki = ak-i(mod N), the (k, i) element of the circular matrix A;
Si = the ith element of the conventional TOF spectrum;
bk = the background in the kth time-channel, its expected value <bk>=b.
(For the sake of simplicity both the number of the elements of the PRBS and the number of the time-channels are chosen as N, i.e. the indices k, i = 0, 1, …, N-1.)
Using eq. (1), the spectrum {Si } can be reconstructed from eq. (2):
Si = [m(1-c)] -1 Σ aki Zk – c[m(1-c)] -1 Σ Zk – b/m (3)
Expression (3) means that the inverse transformation of the correlation (2) is its transposition, regardless of an additive constant and of a factor. A great advantage of the correlation method using PRBS is that the solution of the linear system (2) can be obtained without a difficult inversion procedure. This is true until we do not want a further elaboration of the data. From the mathematical-statistical viewpoint PRBSs yield the minimal total error for the spectrum {Si} 3, 4).
In the theoretical research of the correlation TOF method the gain factor arises as a fundamental category. Its square
gi2 = <Δ2Siconv>/<Δ2Sicorr>, (4)
the ratio of the variances of the ith point of the spectrum obtained by the conventional and the correlation methods and it expresses the superiority of the correlation technique over the conventional one. With the help of the gain factor one can optimize the correlation measurement through its controllable parameters. It has been performed in one- and two-dimensional cases5-8) and for filtration9, 10).
In the use of the gain factor (4) all the elements Si are supposed to be reconstructable by the transformation (3). It has been shown 5, 10) that the gain for Si increases due to the increase of the relative signal height. This means that the identification of a small peak is difficult, if the spectrum contains a large nuisance peak increasing the “background of ignorance” 6), even in the case of optimally chosen parameters of the PRBS. Here a successful filtering of the large peak increases the relative signal height of the low intensity portions, improving the gain for the points of the interesting small peak.
The most widespread method deviced for filtration is the use of an additional chopper and of gating electronics following the filtering chopper 9, 10). In the generally accepted opinion a chopper that has the same PRBS as the correlation one should not be used for filtration 11). One goal of this paper is to prove that one chopper can fulfill both the task of correlation and that of filtration.
2. Filtration and correlation with one chopper
In an arrangement an electronic gate is set up along with the correlation chopper so that the count of an event, if it may be involved in the unwanted Si, say in Se, is prevented by the gate. Mathematically, the counts will be the following:
Zk =(1- aki )Σ aki Si + bk , k=0,1,2,…,N-1 (5)
Amadori and Hossfeld tried to construct the Si‘s by a procedure similar to the transformation (3) and they established that the construction would be possible only if the third correlation function
Caaa(u,v) = Σ aj aj+u aj+v (6)
of the PRBS had only three levels10, 11). They declared that among the known PRBS, provided N < 1024, they did not find such a sequence. Incidentally, I found one in their earlier cited paper 5) that meets the above-mentioned criterion. However, it is easy to prove that every PRBS with h = 1 meets this criterion.
Despite of this finding, I suggest a new method which still keeps the arrangement of one chopper but changes the mathematical elaboration of the experimental data.
How can system (5) be solved, taking into account that the inverse of its matrix does not exist? (5) is an m-fold underdetermined system for the Si‘s, the rank of its matrix is N – m.
An evident way is to make m preliminary assumptions or constraints. For example, if we know that in a given part of the spectrum there appears no peak, we can assume that these Si‘s are constant.
There is another way, which seems to be the most adequate one in the mathematical-statistical sense and it is straightforwardly directed to the goal of the investigations. As the final result of the experiments one wants to obtain estimations for some parameters of the spectrum, such as for maximum height, the position and the standard deviation of a peak, rather than for the spectrum itself. Without filtration these few parameters can be estimated either through the usual two steps: decorrelation (3) and a fitting procedure, or through a direct fitting procedure12). In the case of the two steps filtration the decorrelation has to overcome the problem of underdetermination, which does not occur in the direct fitting. For the latter let us have an example with the assumption that the spectrum is well described by the sum of J Gaussians:
Si(P) = Σ Ar exp[-(i-Mr)2/Dr2], i=0,1,2,…,N-1; r=1,2, …, J (7)
i.e. Si (P) is the expected value of Si. P is the vector of the 3J + 1 parameters; as the (3J + 1)st parameter, the expected value b of the background is taken.
The expected value of eq. (5):
Zk = { | Σ aki Si (P)+ bk 0 |
if ake=0 if ake=1 |
(8) |
for which the parameters can be estimated by the weighted least-squares criterion:
Σ wk {Zk – Σ aki [Si (P)+ b/m]}2 => min (9)
where weighting by the reciprocals of the counts yields the best estimation for P in the sense of having the (3J + 1) dimensional dispersion ellipsoid of minimum volume. On the other hand, the filtration is also performed by the mathematical weighting, and electronic gating is unnecessary:
wk = { | 1/Zk 0 |
if ake=0 if ake=1 |
(10) |
where e is the index of the spectrum element to filter out. The solution of eq. (9) with (10) can be obtained by the iterative solution of the matrix equation
GTATWAGΔP = GTATWΔZ,
where
G is a N x (3J + 1) matrix, its elements
Giu = δSi(P)/δPu P = Pf i = 0, 1, …, N-1; u = 1, 2, …, 3J + 1;
ΔP is the vector of the changes, its elements
ΔPu = Pu-Pfu u = 1, 2, …, 3J + 1;
Pf is the vector of the parameters, determined in the previous step of iteration;
W is the diagonal matrix of the weights, its elements
Wjk = { | wk if j=k 0 if j ≠ k |
ΔZ is the vector of differences between the measured and the computed counts, and the superscript T denotes transposition.
Here I mention without proof that filtration by omission of measured data generally implies a loss of information theoretically in the sense that the determinant of the information matrix M = GTATWAG decreases due to leaving out symmetric dyad(s) from the positive definite matrix. Consequently, the volume of the dispersion ellipsoid of the estimated parameters increases. It is approximately equal to the expression (detM)-1/2 . Notwithstanding this theorem, the filtration can turn out to be useful in practice. Although the correlation in eq. (2) is supposed to take place by a rectangular chopping function, with the only values of zero and one, the real modulating functions are not of this form. This inperfection results in the so-called input-transducer error5). This effect can cause a stronger disturbance than the background of ignorance which also decreases by filtering the large uninteresting peak.
Now let us calculate the losses of intensity, which result from the two types of filtration. The arrangement of the additional filtering chopper with parameters (Nf, mf, cf) and of the following gate transmits the portion (m / N) (mf / Nf) (1 – mf / Nf) ≈ ccf (1 – cf) of the incident intensity. Its maximum is c/4. In the new method of one chopper the remaining part is: (m / N)(1- m / N) ≈ c(1 – c) times the incident intensity. (Of course, the optimum is influenced not only by the maximum intensity, but by the signal-to-background rate, too.)
If the ignorable peak is so wide that it is worth to cut out more than one channel, say s channels, the remaining portion is approximately (m / N)(1 – sm / N) ≈ c(1 – sc) for the new method, which has the maximum of order ¼ s, if c = s/2. The filtration of the channels e1, e2, …, es is performed by the weights
Wjk = { | 1/Zk if ake1= ake2 = … akes=0 0 otherwise, |
instead of eq. (10). For two filtered channels the exact value of the remaining portion is (m / N) [1 – (2m – h) / N]. For more channels the parameters of the autocorrelation functions of higher order should also be considered for the exact calculation of loss in intensity. For more channels one must consider that by increasing the number of the zero weights, the information matrix M can become singular. To avoid this, the duty cycle c should be decreased. The optimization is not too difficult by computer trials and it can also be performed analytically for fixed models [e.g. for spectrums following eq. (7)]. The use of a smaller duty cycle does not mean any disadvantage for the reason that filtering is effective only for low background, and in these cases a smaller than 1/2 duty cycle can be optimal6) . If c is small, the posterior filtration outlined here does not erase too many Zk, so it can be extended to more channels without making the system underdetermined.
3. Dropped-out points
The posterior filtration can be applied for the appearance of dropped-out points, too. Because of electronic error some Zk can become false (e.g. contents of some time-channels are erased). These channels can be filtered by the weights wk = 0, if Zk is dropped out. Here the decorrelation (3) is of no use. This is shown by computer-simulated measurements. We produced a TOF spectrum by eq. (2), adding Poisson-distributed errors to the Zk ‘s. The parameters of the used PRBS: N = 31, m = 6, c = 1/6. The model function: Si = A exp[– (i – M)2 / D2] with A = 100, M = 6, D2 = 3 and with background b = 100. Half of the measured points Zk were simulated as rejected points. These randomly chosen channels got the mean content Z . (One could not do better.) The spectrum {Si } was reconstructed by eq. (3), the traditional decorrelation (Fig. 1), and by my method according to eq. (9) (Fig. 2).
|
||||||||
The former resulted in an unevaluable spectrum. The latter yields nearly the exact parameters. The mean square deviations of the parameters are about √2 time the MSDs given by the simulation without dropped-out points, but still they are sufficiently small (table 1). It means that the proposed evaluation of experiments with lost counts provides a good result that reflects the loss in accuracy proportionally to the loss of information. As opposed, the formerly used method results in total loss of accuracy due to some loss of information.
|
||||||||||||||||||||||||
Table 1. Estimation for the parameters of TOF spectrum, with posterior filtration of the dropped-out points. In the simulation 50% of the measured points were dropped out. A is the amplitude, M is the abscissa, D2 is the square width of the peak, b is the background level.
4. Summary
A new method is introduced for filtering posterior the correlation TOF spectrum. It has the following advantages over the usual decorrelation technique:
- The loss of intensity is significantly less.
- It is no use to add another chopper and gating electronics to the usual TOF correlation arrangement.
- Filtering is performed after the measurement, so it is adjustable, the channels to be filtered can be chosen a posteriori, and the unfiltered spectrum can also be reconstructed.
- It fits well to the data evaluation process, yielding the best estimation for the parameters of the spectrum.
- The dropped-out points can also be filtered by this technique contrary to the usual decorrelation method.
References
- L. Pal, N. Kroo, P. Pellionisz, F. Szlavik and I. Vizi, Symp. on Neutron inelastic scattering, Copenhagen, (IAEA, Vienna, 1968) vol. 2, p. 407.
- F. Gompf, W. Reichardt, W. Glaser and K. Beckurts, ibid: p. 417.
- G. Wilchelmi and F. Gompf, Nucl. Instr. and Meth. 81 (1970) 36.
- K. Weise, Nucl. Instr. and Meth. 98 (1972) 119.
- F. Hossfeld, R. Amadori and R. Scherm, Proc. Symp. on Instrumentation of neutron inelastic scattering (IAEA, Vienna, 1970) p. 117.
- R. von Jan and R. Scherm, Nucl. Instr. and Meth. 80 (1970) 69.
- R. Amadori and F. Hossfeld, Symp. on Neutron inelastic scattering, Grenoble, IAEA/SM-155/F-5 (1972).
- F. Hossfeld and R. Amadori, Berichte der Kernforschunganlage Julich (1972).
- P. Pellionisz, N. Kroo and F. Mezei, Symp. on Neutron inelastic scattering, Grenoble, IAEA/SM-155/F-8 (1972).
- R. Amadori and F. Hossfeld, Nucl. Instr. and Meth. 124 (1975) 429.
- R. Amadori, Berichte der Kernforschungsanlage Julich, 1050-FF (1974).
- L. Cser and L. J. Naszodi, following article.