# Detector Dead Time Determination and Optimal Counting Rate for a Detector Near a Spallation Source or a Subcritical Multiplying System

Authors: V. Bécares; J. Blázquez
Journal: Science and Technology of Nuclear Installations
 (2012)

---

## Abstract

The operation of accelerator-driven systems or spallation sources requires the monitoring of intense neutron fluxes, which may be billions-fold more intense than the fluxes obtained with usual radioactive sources. If a neutron detector is placed near a very intense source, it can become saturated because of detector dead time. On the contrary, if it is placed far away from the source, it will lose counting statistics. For this reason, there must exist an optimal position for placing the detector. The optimal position is defined as the one with the minimal relative uncertainty in the counting rate. In this work, we review the techniques to determine the detector dead time that can be applied with an accelerator-driven subcritical system or a spallation source. For the case of a spallation source, counting rates do not follow Poisson's statistics because of the multiplicity of the number of neutrons emitted by incident proton. It has been found a simple expression that relates the optimal counting rate with the source multiplicity and the uncertainty in the determination of the dead time.


---

## Body

## 1. Introduction

A strong interest in intense accelerator-driven neutron sources for different applications has grown in the last decades. Although a number of reactions is available for neutron production with accelerated particles (photoneutrons, D-Be, and others), preferred technology for these neutron sources is spallation, because it provides the highest neutron yield per unit of energy of the incident particles. Neutron spallation sources consist of a high-energy proton accelerator coupled to a high Z target. Typically, a 1000 MeV incident proton beam could produce ~50 nucleons per incident proton, approximately half of  them neutrons [1]. In this way, very intense neutron fluxes can be obtained (over 1015 n/cm2 s). Although current mode detectors are generally preferred for monitoring such high fluxes, pulse mode detectors are also of interest because of their ability to record information from the individual detector signals. However, if a pulsed mode detector is used to monitor such high fluxes, it can become saturated because of dead time effects.An important application of spallation sources is to drive accelerator-driven systems (ADSs). An ADS consists of a subcritical reactor driven by an external spallation neutron source. The interest in these systems has risen since the 1990s [2–6] as it has been recognized their ability to reduce the volume and radiotoxicity of high-level nuclear waste. Operation of ADSs requires the monitoring of the neutron source intensity and the neutron flux within the reactor, for instance, to apply reactivity monitoring techniques such as the current-to-flux technique or pulsed neutron source (PNS) techniques [7–9]. Again, for this purpose, it is required to measure neutron fluxes high enough to likely cause saturation of the detectors used.Therefore, the operation of both spallation sources and ADSs requires the monitoring of intense neutron fluxes that can cause detector saturation problems. Hence, a good characterization of the detector dead times becomes mandatory.The measurement of dead times is usually accomplished using two constant-intensity neutron sources (the so-called two-source method). However, in Section2 we propose alternative techniques that are suitable to be directly applied to measurements in spallation sources or ADS, without the need of dedicated calibration experiments and that can provide a better characterization of the dead time of the detector system than the two-source method. These techniques have been applied to the experimental results obtained at the Yalina-Booster subcritical facility [10, 11] during the EUROTRANS [12] experimental campaign carried out in this facility.Finally, the combination of the saturation effects at high counting rates and the poor statistics at low counting rates cause that there must be an optimal counting rate with a minimum relative uncertainty. The determination of this optimal counting rate for the case of an spallation source is complicated by the non-Poissonian nature of the spallation process. This problem will be addressed in Section3.


## 2. Measurement of Dead Time

The dead time of a detector is defined as the minimum time interval that two consecutive counts must be separated in order to be recorded as two different events. The effect of having a dead time in a detector used to monitor counting rates is that the measured counting rates will be lower than the real ones. However, the real counting rate can be determined from the measured one if the dead time of the detector is known.To obtain a formula for relating the real and the measured counting rates in a detector system affected by dead time, it must be taken into account that there are two fundamental behaviors for the dead time. They are referred asparalyzable and non-paralyzable [13, 14] (also called extendable and non-extendable, resp.). In a detector affected by paralyzable dead time, events that occur during the dead time of a previous one, and consequently not recorded, also produce a dead time. In other words, the effect of this last event can be regarded as extending the dead time of the previous one. On the contrary, in a non-paralyzable detector, events occurred during the dead time of a previous event are neither detected nor cause an additional dead time. For instance, fission chambers are affected by paralyzable dead time because a fission occurring during the time that the detector is ionized after a previous event also causes additional ionization, thus extending the dead time. On the contrary, some elements of the electronic chains can introduce a non-paralyzable dead time if they ignore the counts arrived during the time they are processing or recording a previous one.The relationship of the real counting rates with the measured counting rates is well known for these two basic models. If we denote byN the real counting rate, by M the measured counting rate, and by τ the dead time and considering that the non-dead time disturbed distribution is Poissonian,   M and N are related by
(1)N=M1-Mτ,
for the case of a non-paralyzable dead time, and by(2)M=Ne-Nτ,
for the case of a paralyzable dead time.These formulae are derived for an initially Poissonian process. This is the case of radioactive decay, but not the case of either fission or spallation. In the case of nuclear fission, the probability to have a fission in the immediate instants after a previous one is higher than that  for the Poissonian process because of the secondary neutrons emitted in the fission process. In fact, the non-Poissonian nature of the fission process, constitutes the basis of the neutron noise techniques. Nevertheless, for dead time correction purposes the deviation from the Poisson distribution is small enough and the formulae derived for the Poisson process can be applied for fission. The spallation process is not a Poissonian process. This fact will be further discussed in Section3.In practice, however, these two models of dead time are idealized models, and actual detectors systems have neither a paralyzable nor a non-paralyzable dead time, but a combination of both. Furthermore, detector systems comprise additional elements (amplifiers, etc.) that add additional dead times to the dead time of the detector itself. Hence, it is common to study the behavior ofseries arrangements of dead times. Quite obviously, a series arrangement of dead times is only relevant when the first dead time is shorter than the following ones, because otherwise only the first dead time is to be taken into account. See, for instance, [14] for further details. Notice that for small dead time corrections, that is, if τ≪1/N, both correction formula (1) and (2) can be approximated by the same first-order Taylor term:(3)N=M(1+Mτ).And therefore, first-order dead time corrections are independent on the model of dead time. The relative value of the dead time correction considering each one of (1), (2), and (3) for different counting rates is presented in Figure 1.Figure 1
Relative value of the dead time correction for different counting rates considering (1), (2), and (3).Several techniques are available for the experimental determination of dead times [13, 14]. The most simple one is possibly the two-source method, described in Section 2.1. Beside these techniques, in this study we suggest another two techniques for dead time determination. The first technique (Section 2.2) is applicable with variable intensity sources such as spallation sources or ADSs. With the second technique (Section 2.3), the dead time is determined from the distribution of time intervals between consecutive detector counts.


### 2.1. The Two-Source Technique

The most usual technique to measure the dead timeτ is the two-source technique [13, 15]. Let us consider that we have two sources S1 and S2 placed at distances d1 and d2 from the detector, as it is shown schematically in Figure 2.Figure 2
Two-source method.In the absence of sources, the detector gets only the background radiation levelB c/s. When the source S2 is removed and only S1 is present, the detector gets N1+B c/s and detects M1 c/s. Therefore, from (3),(4)N1+B=M1(1+M1τ).Similarly, in absence ofS1 and with only S2 present, we will have that,(5)N2+B=M2(1+M2τ).And with both sources,(6)N1+N2+B=M12(1+M12τ).Working out the dead time from these three equations,(7)τ=M1+M2-M12-BM122-M12-M22.From the uncertainties in the right-hand side (Poisson statistics)’s it is immediate to estimate the uncertainty in the dead time. The backgroundB has been included in case that sources are not intense enough. Usually, it can be neglected, but this simplification introduces a systematic error in the measurement.


### 2.2. Dead Time Determination Varying the Source Intensity

In an ADS or, more generally, in an accelerator-driven neutron source, the ability to change the accelerator intensity can be applied for detector dead time determination. Typically, the neutron flux measured at any point of the assembly is proportional to the source intensity. Therefore, if the detector used to measure the flux within the assembly is not affected by dead time, the measured counting rate in this detector will be linear with the accelerator intensity; sayN=kIA. If it is affected by dead time, the measured counting rate in the detector will be related with the source intensity by an equation similar to (3):(8)kIA=M(1+Mτ)⟹IAM=1k+τkM.Hence, from the linear fit ofIA/M to M, the dead time τ of the detector can be determined. This technique can be also applied using the counting rates in a second detector with a much lower sensitivity so that it is not affected by dead time, instead of the accelerator current IA.A variant of this technique using a single detector can be also applied in an ADS if the source-jerk technique is used for reactivity monitoring. The source-jerk technique [16–18] is based on the kinetic response of the neutron flux in the system to a sudden removal of the external source. In this case, the counting rate in the flux monitors experiences a fast decay due to the decay of the prompt neutrons followed by a much slower decay due to the delayed neutrons, in such a way that it can be considered that the counting rate ends in a constant level in the millisecond scale. This situation is shown in Figure 3, where the neutron level before the source removal is denoted by N0 and the neutron level after the source removal (due to the slowly decaying delayed neutrons) is denoted by N1.Figure 3
Graphical scheme of the source-jerk technique.It can be obtained in the point kinetics model that the reactivity of the reactor (in units of dollars) is related with the ratioN0/N1 by the expression(9)ρ($)=1-N0N1⟹N0=(1-ρ($))N1.Equation (9) constitutes the basis of the source-jerk technique. If we consider that at the counting rate N0 the detector is affected by dead time (at therefore it measures a counting rate M0) but the counting rate N1 is low enough not to be affected by dead time, we can write an expression similar to (8) replacing IA by N1, M by M0, and k by 1-ρ($):(10)N1M0=11-ρ($)+τ1-ρ($)M0.Hence, if it is possible to operate the ADS with different values of the source intensity (and therefore ofN1 and M0), we can use source-jerk experiments to determine both the reactivity of the assembly and the dead time of the detector. An example of the application of this technique to the experimental results in Yalina-Booster is presented in Figure 4.Figure 4
Example of the application of (10) to the source-jerk experiments at the Yalina-Booster subcritical facility.However, the application of (8) or (10) for the determination of the detector dead time, either with two detectors with different sensibilities or with a single detector and the source-jerk technique, has been found to be limited in practice by a number of causes, including the presence of spatial and energy effects that make the point-kinetics model lose validity, source stability requirements before the source removal, and dynamical effects that may cause the reactivity to vary with the power level and thus with the source intensity.


### 2.3. Dead Time Determination from the Distribution of Time Intervals between Consecutive Counts

We have considered another technique to determine the detector dead time that requires a single detector and a single counting rate. This technique arises from the distribution of time intervals between consecutive counts. Modern data acquisition systems such as the one used during the EUROTRANS experiments at Yalina-Booster can register the times of individual detector signals instead of average counting rates. It is well known that for a Poissonian process the distribution of time intervals between consecutive counts (let us denote it byI1(t)) takes the shape of a simple exponential. However, the presence of a dead time alters this distribution. Analytical formulae for this distribution are known for a long time for the cases of a simple paralyzable or non-paralyzable dead time [14]. For the non-paralyzable case this formula reads(11)I1(t)=Ne-N(t-τ),fort>τ;
while, for the paralyzable case (12)I1(t)=N∑j=1J1(j-1)![-N(t-jτ)]j-1e-jNτ,
where J is the largest integer below t/τ. For the case of a series arrangement of dead times, the determination of the analytical shape of I1(t) becomes cumbersome and the reader is referred to the bibliography [19, 20]. Nevertheless, the shape of I1(t) for complex arrangements of dead times can also be determined with Monte Carlo simulations. For this, we have implemented a program that generates events simulating an initial Poisson random process and then applies successively dead times of different types. Some results obtained with this program are presented in Figure 5. They are the cases of a purely paralyzable and non-paralyzable dead time, as well as the cases of series arrangements of two dead times: paralyzable/non-paralyzable and non-paralyzable/paralyzable. Normalization has been chosen so that ∫0∞I1(t)dt=1 for the unperturbed Poisson distribution; with this normalization I1(0) is equal to the counting rate of the unperturbed Poisson distribution. Notice the clear differences in the shape of I1(t) among these cases. This causes that the shape of I1(t) can be used to determine not only the value of the dead time of the system but also its nature (paralyzable or non-paralyzable) even in the presence of complex series of dead times.Distribution functionI1(t) for intervals between adjacent random events (counting rate 106 counts/s, 109 counts in the simulation).
(a)
Unperturbed distribution of time intervals and distribution of time intervals with paralyzable and non-paralyzable dead time (400 ns in any case)(b)
Distribution of time intervals for a series arrangement of a paralyzable and a non-paralyzable dead time (200 ns paralyzable and 400 ns non-paralyzable). The cases of a single paralyzable or non-paralyzable dead time (both of 400 ns) are also shown for comparison(c)
Distribution of time intervals for a series arrangement of a non-paralyzable and a paralyzable dead time (200 ns non-paralyzable and 400 ns paralyzable). The cases of a single paralyzable or non-paralyzable dead time (both of 400 ns) are also shown for comparisonAn example of the application of this technique obtained in Yalina-Booster is presented in Figure6. The determination of the value of the dead time was straightforward, but the determination of the type of dead time or if several types were present was impossible because of the presence of some ringing effects that are visible in the figure. Notice as well that (11) and (12) are derived for a strictly constant counting rate; variations of the source intensity in actual systems also alter these results.Figure 6
Example of the distribution of intervals between consecutive counts obtained from the Yalina-Booster experiments. Gray line represents the fit to a exponential (over a large range than the plotted interval).


## 2.1. The Two-Source Technique

The most usual technique to measure the dead timeτ is the two-source technique [13, 15]. Let us consider that we have two sources S1 and S2 placed at distances d1 and d2 from the detector, as it is shown schematically in Figure 2.Figure 2
Two-source method.In the absence of sources, the detector gets only the background radiation levelB c/s. When the source S2 is removed and only S1 is present, the detector gets N1+B c/s and detects M1 c/s. Therefore, from (3),(4)N1+B=M1(1+M1τ).Similarly, in absence ofS1 and with only S2 present, we will have that,(5)N2+B=M2(1+M2τ).And with both sources,(6)N1+N2+B=M12(1+M12τ).Working out the dead time from these three equations,(7)τ=M1+M2-M12-BM122-M12-M22.From the uncertainties in the right-hand side (Poisson statistics)’s it is immediate to estimate the uncertainty in the dead time. The backgroundB has been included in case that sources are not intense enough. Usually, it can be neglected, but this simplification introduces a systematic error in the measurement.


## 2.2. Dead Time Determination Varying the Source Intensity

In an ADS or, more generally, in an accelerator-driven neutron source, the ability to change the accelerator intensity can be applied for detector dead time determination. Typically, the neutron flux measured at any point of the assembly is proportional to the source intensity. Therefore, if the detector used to measure the flux within the assembly is not affected by dead time, the measured counting rate in this detector will be linear with the accelerator intensity; sayN=kIA. If it is affected by dead time, the measured counting rate in the detector will be related with the source intensity by an equation similar to (3):(8)kIA=M(1+Mτ)⟹IAM=1k+τkM.Hence, from the linear fit ofIA/M to M, the dead time τ of the detector can be determined. This technique can be also applied using the counting rates in a second detector with a much lower sensitivity so that it is not affected by dead time, instead of the accelerator current IA.A variant of this technique using a single detector can be also applied in an ADS if the source-jerk technique is used for reactivity monitoring. The source-jerk technique [16–18] is based on the kinetic response of the neutron flux in the system to a sudden removal of the external source. In this case, the counting rate in the flux monitors experiences a fast decay due to the decay of the prompt neutrons followed by a much slower decay due to the delayed neutrons, in such a way that it can be considered that the counting rate ends in a constant level in the millisecond scale. This situation is shown in Figure 3, where the neutron level before the source removal is denoted by N0 and the neutron level after the source removal (due to the slowly decaying delayed neutrons) is denoted by N1.Figure 3
Graphical scheme of the source-jerk technique.It can be obtained in the point kinetics model that the reactivity of the reactor (in units of dollars) is related with the ratioN0/N1 by the expression(9)ρ($)=1-N0N1⟹N0=(1-ρ($))N1.Equation (9) constitutes the basis of the source-jerk technique. If we consider that at the counting rate N0 the detector is affected by dead time (at therefore it measures a counting rate M0) but the counting rate N1 is low enough not to be affected by dead time, we can write an expression similar to (8) replacing IA by N1, M by M0, and k by 1-ρ($):(10)N1M0=11-ρ($)+τ1-ρ($)M0.Hence, if it is possible to operate the ADS with different values of the source intensity (and therefore ofN1 and M0), we can use source-jerk experiments to determine both the reactivity of the assembly and the dead time of the detector. An example of the application of this technique to the experimental results in Yalina-Booster is presented in Figure 4.Figure 4
Example of the application of (10) to the source-jerk experiments at the Yalina-Booster subcritical facility.However, the application of (8) or (10) for the determination of the detector dead time, either with two detectors with different sensibilities or with a single detector and the source-jerk technique, has been found to be limited in practice by a number of causes, including the presence of spatial and energy effects that make the point-kinetics model lose validity, source stability requirements before the source removal, and dynamical effects that may cause the reactivity to vary with the power level and thus with the source intensity.


## 2.3. Dead Time Determination from the Distribution of Time Intervals between Consecutive Counts

We have considered another technique to determine the detector dead time that requires a single detector and a single counting rate. This technique arises from the distribution of time intervals between consecutive counts. Modern data acquisition systems such as the one used during the EUROTRANS experiments at Yalina-Booster can register the times of individual detector signals instead of average counting rates. It is well known that for a Poissonian process the distribution of time intervals between consecutive counts (let us denote it byI1(t)) takes the shape of a simple exponential. However, the presence of a dead time alters this distribution. Analytical formulae for this distribution are known for a long time for the cases of a simple paralyzable or non-paralyzable dead time [14]. For the non-paralyzable case this formula reads(11)I1(t)=Ne-N(t-τ),fort>τ;
while, for the paralyzable case (12)I1(t)=N∑j=1J1(j-1)![-N(t-jτ)]j-1e-jNτ,
where J is the largest integer below t/τ. For the case of a series arrangement of dead times, the determination of the analytical shape of I1(t) becomes cumbersome and the reader is referred to the bibliography [19, 20]. Nevertheless, the shape of I1(t) for complex arrangements of dead times can also be determined with Monte Carlo simulations. For this, we have implemented a program that generates events simulating an initial Poisson random process and then applies successively dead times of different types. Some results obtained with this program are presented in Figure 5. They are the cases of a purely paralyzable and non-paralyzable dead time, as well as the cases of series arrangements of two dead times: paralyzable/non-paralyzable and non-paralyzable/paralyzable. Normalization has been chosen so that ∫0∞I1(t)dt=1 for the unperturbed Poisson distribution; with this normalization I1(0) is equal to the counting rate of the unperturbed Poisson distribution. Notice the clear differences in the shape of I1(t) among these cases. This causes that the shape of I1(t) can be used to determine not only the value of the dead time of the system but also its nature (paralyzable or non-paralyzable) even in the presence of complex series of dead times.Distribution functionI1(t) for intervals between adjacent random events (counting rate 106 counts/s, 109 counts in the simulation).
(a)
Unperturbed distribution of time intervals and distribution of time intervals with paralyzable and non-paralyzable dead time (400 ns in any case)(b)
Distribution of time intervals for a series arrangement of a paralyzable and a non-paralyzable dead time (200 ns paralyzable and 400 ns non-paralyzable). The cases of a single paralyzable or non-paralyzable dead time (both of 400 ns) are also shown for comparison(c)
Distribution of time intervals for a series arrangement of a non-paralyzable and a paralyzable dead time (200 ns non-paralyzable and 400 ns paralyzable). The cases of a single paralyzable or non-paralyzable dead time (both of 400 ns) are also shown for comparisonAn example of the application of this technique obtained in Yalina-Booster is presented in Figure6. The determination of the value of the dead time was straightforward, but the determination of the type of dead time or if several types were present was impossible because of the presence of some ringing effects that are visible in the figure. Notice as well that (11) and (12) are derived for a strictly constant counting rate; variations of the source intensity in actual systems also alter these results.Figure 6
Example of the distribution of intervals between consecutive counts obtained from the Yalina-Booster experiments. Gray line represents the fit to a exponential (over a large range than the plotted interval).


## 3. Optimal Counting Rate with a Spallation Neutron Source

If the detector is placed very close to a spallation neutron source, it can become saturated; if it is placed too far away, statistical resolution is very small. From here it can be inferred that there must exist an optimal position in which the relative uncertainty of the counting rate is minimal. Let us consider the case of a non-paralyzable dead time. In this case, if we denote bym=MΔt the number of counts in the detector in the interval Δt, we have from (1) that(13)N=m/Δt1-τm/Δt.The uncertainty inN can be estimated from the variances:(14)σN2=(∂N∂m)2σm2+(∂N∂τ)2στ2.Partial derivatives can be computed from (13) as follows(15)∂N∂m=Nm(1+τN),∂N∂τ=N2.For a radioactive source, which follows a Poisson-like distribution of detection probabilities, we have thatσm2=m, but a spallation source, because of the neutron multiplicity μ, the distribution of detection probabilities differs from the Poisson distribution. In a first approximation, we have that [21](16)σm2=m(1+Dμ),
where D is the Diven factor, which is close to the unity. Notice that if the neutron multiplicity is zero, the variance corresponds to that of a Poisson process. With these conditions,(17)σN2=N2m2(1+τN)2m(1+Dμ)+N4στ2.Replacingm by the first-order approximation m≃NΔt(1-Nτ) and approximating the first term in the right-hand side to the first order in Nτ, we find the following expression for the relative error in σN:(18)εN2=(σNN)2=(1+3τN)N(1+Dμ)Δt+N2στ2.The optimal counting rate is the one that minimizes the relative errorεN:(19)∂εN2∂N|N=Nop=0⟹-1Nop2(1+DμΔt)+2Nopστ2=0.Working outNop in this equation, we finally find that(20)Nop=1+Dμ2στ2Δt3.Typical values for the problem are given in Table1. With these values, optimal counting rate is 1.05×105(c/s). Correction by dead with these parameters is about 20%, in the limit of the first-order approximation (3). The only parameter one can change in (20) is Δt so if we choose 10 s instead of 1 s, the new optimal counting rate becomes 0.48×105(c/s), and the correction is now 10%, low enough for using (3). Therefore, the detector must be placed at points with counting rates about 4×104(c/s). Diven’s factor can be taken as the unity without a large error in the computation of the optimal counting rate.Table 1
Multiplicity, Diven’s factor dead time, and measurement time for a typical problem.
μDτ (μs)Δt (s)201.12±0.11


## 4. Conclusions

In addition to traditional techniques for dead time measurement, we propose in this work two alternative methods to determine detector dead times. The first of them can be applied with variable intensity sources, such as accelerator-driven neutron sources. This method consists in studying the relationship between the detector counting rate and the source intensity that results of the presence of dead time in the detector system. A variation of this technique can be applied measuring the counting rate in a second detector with a lower counting rate (i.e., with negligible dead time effects) instead of the source intensity. Another variation of this technique, which does not require either a second detector or knowing the source intensity, can be applied in an ADS making most of the results of the source-jerk technique, used to determine the reactivity of the system.The second method we propose to determine the dead time of a detector system can be applied if the data acquisition system registers the times of individual detector events instead of counting rates. In this way the distribution of time intervals between consecutive counts can be obtained, from which the detector dead time can be determined. This method has the advantage of requiring no additional detector and only one counting rate. In addition, this method is in principle capable of determining the type of dead time (paralyzable or non-paralyzable) even in the case of complex arrangements of dead times.Finally, we remark that neutrons produced in a spallation source do not follow the Poisson statistics because of the multiplicity in the neutron production. An expression has been found that relates the multiplicity, the measurement interval, and the variance of the detector dead time. This expression allows placing neutron monitors where the relative error is minimal.



---

*Source: 240693-2012-04-29.xml*