research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Empirical dead-time corrections for synchrotron sources

CROSSMARK_Color_square_no_text.svg

aAdvanced Photon Source, Argonne National Laboratory, Argonne, IL 60439, USA
*Correspondence e-mail: d-walko@anl.gov

(Received 3 April 2008; accepted 16 July 2008; online 3 October 2008)

An experimental comparison of models for performing dead-time corrections of photon-counting detectors at synchrotron sources is presented. The performance of several detectors in the three operating modes of the Advanced Photon Source is systematically compared, with particular emphasis on asymmetric fill patterns. Several simple and well known correction formulas are evaluated. The results demonstrate the critical importance of detector speed and synchrotron fill pattern in selecting the proper dead-time correction.

1. Introduction

The dead-time of a photon-counting detector is the time needed for a detector to recover after counting a photon, before another photon can be counted. Dead-time corrections extend the usable range of a detector into a high-count-rate regime in which X-rays arrive too fast for all to be individually counted. Models for such corrections were originally based on Poisson statistics; that is, the intervals between X-rays arriving at the detector are assumed to be random and independent. In such cases, as a tube source, the dead-time constant τ is determined simply by how quickly the detector can recover after a counting event. At synchrotron X-ray sources, the time intervals between X-ray arrivals are not random but are determined by the fill pattern of electron bunches in the accelerator ring. For the most efficient use of a photon-counting detector, dead-time corrections which take into account the time structure of the X-ray source must be used. This article assesses various dead-time corrections for the main operating modes of the Advanced Photon Source (APS), comparing measurements made with two scintillator detectors and an avalanche photodiode (APD) detector. To our knowledge, this is the first systematic measurement of dead-time effects which studies variations in both the time structure of the source and the speed of the detector response. This is also, we believe, the first experimental investigation of dead-time effects in an asymmetric fill pattern. The results can be generalized to other synchrotron fill patterns and other detectors and could be used to optimize the data collection rate of an experiment.

The purpose of this paper is to determine the simplest and best dead-time correction for a given experimental condition, that is, for a given detector and the given source's pulse structure. We report tests on detectors with significantly different characteristics (speed and use of pulse-height analysis) in the various operating modes of the APS. §2[link] describes various dead-time models appropriate for synchrotron sources. In §3[link] and §4[link] we describe the fill patterns of the APS and our experimental set-up, respectively. In §5[link] we present the results and compare the accuracy of the various dead-time models.

2. Dead-time models

The model used in a dead-time correction will depend on the characteristics of the detector and the source. Several authors (e.g. Arndt, 1978[Arndt, U. W. (1978). J. Phys. E, 11, 671-673.]; Quintana, 1991[Quintana, J. P. (1991). J. Appl. Cryst. 24, 261-262.]; Cousins, 1994[Cousins, C. S. G. (1994). J. Appl. Cryst. 27, 159-163.]; Kishimoto, 1997[Kishimoto, S. (1997). Nucl. Instrum. Methods Phys. Res. A, 397, 343-353.]; Knoll, 2000[Knoll, G. F. (2000). Radiation Detection and Measurement, 3rd ed. New York: John Wiley and Sons.]; Bateman, 2000[Bateman, J. E. (2000). J. Synchrotron Rad. 7, 307-312.]; Ida & Iwata, 2005[Ida, T. & Iwata, Y. (2005). J. Appl. Cryst. 38, 426-432.]) have developed a variety of dead-time models; in this section we will not derive them in detail but will review the models most appropriate for synchrotron sources. We will also present the effective dead-time τ for a given model, which may depend on the intrinsic detector speed, the extrinsic photon pulse rate, and whether pulse-height analysis (PHA) is used to discriminate between single-photon and multiphoton events. The ultimate goal is to be able to efficiently determine NT, the true count rate, from NO, the observed count rate, for a given experimental situation.

As a rule of thumb, at a source with bunch separation T, NO is the product of three terms: the arrival frequency 1/T, the probability of a countable event, and the fraction of time that the detector is rendered inoperative (`dead') owing to dead-time effects. By `countable event' we mean a photon arriving at the detector which will be counted by the pulse-processing electronics, which depends on whether PHA is used. PHA is often implemented for scintillator detectors at synchrotrons to avoid contamination of the signal by higher harmonic photons; with PHA, a countable event is when one single photon arrives from a given bunch in the storage ring. (More advanced PHA techniques, which would determine the number of photons arriving from a given bunch, may lead to a complete solution for dead-time issues but will not be considered here.) The photons from a single bunch can be considered to have a Poisson distribution (Knoll, 2000[Knoll, G. F. (2000). Radiation Detection and Measurement, 3rd ed. New York: John Wiley and Sons.]). We use P(j) to represent the probability of j photons arriving from one bunch, calculating it using Poisson statistics. The probability of arrival of one photon is

[P(1)=N_{\rm{T}}T\exp\left(-N_{\rm{T}}T\right),\eqno(1)]

with the mean of the Poisson distribution being NTT. However, if PHA is not used, as is often the case with APDs, then a countable event occurs whenever one or more photons from a given bunch arrives at the detector. The probability of any number of photons arriving is the probability that zero photons do not arrive, that is,

[\textstyle\sum\limits_{j\,=\,1}^\infty P(\,j)=1-P(0)=1-\exp\left(-N_{\rm{T}}T\right).\eqno(2)]

There are several ways to model the dead-time-induced inefficiency of a detector. In the non-extended dead-time model, the detector is considered inoperative for a fixed period of time τ after each observed photon; the dead-time-related efficiency of the detector is thus (1 − τNO). In the simplest and perhaps most common formulation of non-extended dead-time, the observed count rate is simply the product of the true count rate and the efficiency factor, which rearranges to give

[N_{\rm{T}}={{N_{\rm{O}}}\over{1-\tau N_{\rm{O}}}}.\eqno(3)]

Several authors (Arndt, 1978[Arndt, U. W. (1978). J. Phys. E, 11, 671-673.]; Cousins, 1994[Cousins, C. S. G. (1994). J. Appl. Cryst. 27, 159-163.]; Kishimoto, 1997[Kishimoto, S. (1997). Nucl. Instrum. Methods Phys. Res. A, 397, 343-353.]) have presented modifications of this formula to take into account the time structure of a synchrotron source.

However, most X-ray photon-counting detectors exhibit extended dead-time (Quintana, 1991[Quintana, J. P. (1991). J. Appl. Cryst. 24, 261-262.]; Bateman, 2000[Bateman, J. E. (2000). J. Synchrotron Rad. 7, 307-312.]). In this case the detector is freshly rendered inoperative after each arriving photon. The detector recovery is delayed by the arrival of any photons during the detector's intrinsic recovery time τd. Together, τd (which is due to the physical response of the photon-sensitive element plus the speed of the electronics) and T, the time separation of the synchrotron pulses, make up the two basic timescales of interest in dead-time corrections. We assume that the duration of a single X-ray bunch, of the order of 100 ps at the APS, is by far short enough to be neglected. According to the rule of thumb discussed above, the observed count rate for the extended dead-time model using PHA will be (1/T)P(1)P(0)n−1. The integer n, which reflects the discrete nature of the source, is the number of bunches that occur in time τd. If a photon in a given bunch is to be counted, no photons must have arrived in the previous amount of time τd, that is, in the previous n − 1 bunches. The extended dead-time model, with PHA, then simplifies to

[N_{\rm{O}}=N_{\rm{T}}\exp\left(-{\tau}N_{\rm{T}}\right),\eqno(4)]

where the effective dead-time is nT, or

[\tau=T\left[{\rm{Int}}\left(\tau_{\rm{d}}/T\right)+1\right].\eqno(5)]

Here, the Int(x) function is the integer less than or equal to x. Equations (4)[link] and (5)[link] should generally be applicable to scintillators at synchrotron sources. The solution of (4)[link] for NT is called the Lambert W function (see Corless et al., 1996[Corless, R. M., Gonnet, G. H., Hare, D. E. G., Jeffrey, D. J. & Knuth, D. E. (1996). Adv. Comput. Math. 5, 329-359.]); strategies for solution are discussed in Appendix A[link]. Note that when the detector's intrinsic dead-time is less than the bunch separation, equation (5)[link] simplifies to τ = T and the bunches are well separated. In the opposite limit, when τd [\gg] T, equation (5)[link] yields ττd as the source effectively appears continuous and the bunch structure becomes irrelevant.

Without PHA, single- and multiple-photon events are not distinguished and therefore the extended dead-time model gives an observed count rate of (1/T)[1 − P(0)]P(0)n−1, or

[N_{\rm{O}}=(1/T)\left[1-\exp\left(-TN_{\rm{T}}\right)\right]\exp\left[-(n-1)TN_{\rm{T}}\right].\eqno(6)]

Although equation (6)[link] does not have an analytic solution for NT, a recursive solution similar to that presented in Appendix A[link] can be derived. An important limiting case of equation (6)[link] occurs when the intrinsic detector response is significantly shorter than the separation of the synchrotron bunches (τd < T). This `isolated' model, with n = 1, should represent APDs under most operating modes of a synchrotron, and simplifies to

[N_{\rm{T}}=-(1/T)\ln\left(1-N_{\rm{O}}T\right),\eqno(7)]

where the effective dead-time is simply the bunch period T. Equation (7)[link] remains applicable for situations in which only photons from specific bunches are counted, such as an isolated bunch in a synchrotron's asymmetric fill pattern (see §2[link]). Specific examples include the incident beam being mechanically chopped to periodically block some X-ray bunches, or the detector output being electronically gated to match the pump frequency in a pump–probe experiment. In these cases, T in (7)[link] becomes the period of the bunch(es) of interest. Fig. 1[link] compares the non-extended, extended and isolated models described here. The observed count rates for the non-extended and isolated models eventually saturate at 1/τ for high count rates, with the isolated model saturating quicker. The observed count rate for the extended model will reach a maximum and then decrease when NT ≥ 1/τ.

[Figure 1]
Figure 1
Comparison of observed and true count rates for the non-extended (solid line), extended (dashed line) and isolated (dash-dot line) dead-time models. The axes are normalized to the dead-time τ.

The dead-time correction of an asymmetric fill pattern is more difficult to model, since T varies between bunches. Lee & Mills (1992[Lee, W. K. & Mills, D. M. (1992). Rev. Sci. Instrum. 63, 1225-1228.]) derived a model in which they divided the fill pattern of charged particles in the accelerator into m groups, where a `group' could be a collection of closely spaced bunches, not necessarily individual bunches. In this model, bunches within a time τd can be grouped together, and are separated from other groups by at least τd. With αi representing the fraction of charge in the ith group, and Tperiod representing the round-trip time of a bunch in the accelerator, then

[N_{\rm{O}}=N_{\rm{T}}\textstyle\sum\limits_{i\,=\,1}^m\alpha_i\exp\left(-\alpha_iN_{\rm{T}}T_{\rm{period}}\right)\eqno(8)]

when PHA is used, and

[N_{\rm{O}}=1/T_{\rm{period}}\textstyle\sum\limits_{i\,=\,1}^m\left[1-\exp\left(-\alpha_iN_{\rm{T}}T_{\rm{period}}\right)\right]\eqno(9)]

without PHA. In effect, different bunches in an asymmetric fill pattern will experience different dead-time effects. Again, this model requires τd to be longer than the group duration and shorter than the group separation. Monte Carlo simulations have been used to predict the dead-time behavior of asymmetric modes (Bateman, 2000[Bateman, J. E. (2000). J. Synchrotron Rad. 7, 307-312.]), but to our knowledge these formulas have not been tested empirically.

3. Time structure of the APS ring

The APS accelerator is 1104 m in circumference, yielding a round-trip time of Tperiod = 3.682 µs; it normally operates with a current of 102 mA. The circulating electrons are grouped into bunches which can, in principle, fill the ring in a wide variety of patterns in any of its 1296 buckets. In practice, the APS in recent years has operated with three fill patterns. The standard operating mode, accounting for about 75% of operating time, fills the ring with 24 equally spaced bunches. The nominal current in this `24-bunch mode' is thus 4.2 mA per bunch, with a spacing of Tperiod/24 = 153 ns between bunches. In `324-bunch mode', every fourth bucket is filled with 0.31 mA per bunch for a separation of Tperiod/324 = 11.4 ns between bunches.

The final mode is an asymmetrically filled mode optimized for timing experiments. In this `hybrid singlet mode', one singlet bunch is isolated from the others; this bunch is filled up to 16 mA. On the opposite side of the ring, eight sets of seven bunches (`septets') are filled equally with the remaining current, for 10.8 mA per septet. The singlet is separated from the nearest septets by a 1.59 µs gap on either side. The septets arrive at intervals of 68 ns, including the 17 ns width of each septet.

4. Experiment description

This series of experiments was performed at beamline 7ID of the APS. A double-crystal monochromator selected 10 keV X-rays from the undulator source. The second crystal of the monochromator was detuned to reduce spectral contamination from higher harmonics. In the 7ID-C experimental enclosure, motorized slits determined the incident X-ray flux, which was monitored by an ion chamber. The amplifier on the ion chamber was always set to keep its operation in the linear regime. Following the ion chamber, an amorphous material (Kapton or nylon) scattered X-rays into one of the photon-counting detectors used in this study, i.e. two scintillator detectors and one avalanche photodiode.

The two scintillator detectors are commonly used, commercially available detectors (Cyberstar). The first is a Tl-doped NaI scintillator crystal and an integrated electronics unit (model X1000, which includes an analog shaping amplifier, single-channel analyzer and high-voltage power supply). The peaking time of the amplifier on this detector, which we refer to as the `NaI' detector, was 0.3 µs for all measurements reported here, which is well matched to the intrinsic decay time constant of NaI of about 230 ns. The second, referred to as the `YAP' detector, is a Ce-doped YAlO3 scintillator crystal with X2000 electronics. While YAP has poorer energy resolution and less efficiency at lower energies compared with NaI, its decay time constant is much shorter at 27 ns; a 50 ns peaking time was used for this detector. For both detectors, PHA was used by setting lower and upper discriminators in single-channel analyzers.

The APD detector used in this experiment was a model EG&G C30707 10 mm × 10 mm silicon APD built into a home-made amplifier circuit and aluminium housing that provided sufficient gain for single-photon counting with a constant fraction discriminator (CFD). Before the CFD, the amplifier produced a ∼10 ns duration electrical pulse; the CFD output duration was ∼30 ns.

In each measurement the count rates of the photon-counting detector and the monitoring ion chamber were recorded as the monochromatic X-ray flux was scanned from low to high. These measurements were performed for all detectors over the course of several months for each of the operating modes described above. The data collected in each operating mode were then fit to equations (3)[link] and (4)[link] for the scintillators, and equations (3)[link] and (6)[link] or (7)[link], as appropriate, for the APD; the two fitting parameters were the effective dead-time τ and an overall scaling factor. This scale factor was first estimated in the low-count-rate regime, as the ratio of the count rates of the photon-counting detector and the ion chamber. The exact value of the scale factor depended on parameters such as the detector geometry, thickness of the amorphous material, and gain of the ion chamber's pre­amplifier; as such it varied from measurement to measurement. In the following section we present the true count rate NT incident upon the scintillator detectors as the count rate of the ion chamber multiplied by the scale factor.

5. Results

Comparisons of the observed to the true count rates for both detectors in each mode are shown in Figs. 2[link], 3[link] and 4[link]. Note the differences in scales on the three figures as the faster detectors progressively operate to higher count rates. The relatively slow NaI detector has almost the same behavior in the 324- and 24-bunch modes, since the dead-time correction is dominated by the detector response time. On the other hand, the YAP and APD detectors are sensitive to differences between these symmetric modes. Meanwhile, the dynamic range of all detectors is relatively hampered in hybrid singlet mode. Also shown are best fits of the data to several models: the non-extended model [equation (3)[link]] for all three detectors; the extended model with PHA [equation (4)[link]] for the scintillators, and the extended model without PHA [equation (6)[link]] or the isolated model [equation (7)[link]] for the APD. Although both the YAP and APD detectors are fast enough to isolate bunches in 24-bunch mode, the use of PHA for the YAP but not the APD requires different models. Certainly, the non-extended dead-time model follows the general trend for all three detectors and could be used for a detector being operated somewhat beyond its linear response regime. However, the non-extended dead-time model tends to underestimate the correction in the middle of the range, typically by 2–5%, and overestimate it in the high-count-rate end of the range by 5% or more. The R-factor for fitting with the non-extended model was generally two to four times greater than for the more appropriate model.

[Figure 2]
Figure 2
Response of the NaI detector in the three operating modes of the APS. The count rate of the ion chamber is scaled to that of the scintillator in the low-count-rate regime, and as such is treated as the `true' count rate of the scintillator. For each operating mode, models of non-extended dead-time (dotted lines) and extended dead-time (solid lines) are shown. The dashed line for hybrid singlet mode is the model of Lee & Mills (1992[Lee, W. K. & Mills, D. M. (1992). Rev. Sci. Instrum. 63, 1225-1228.]), as explained in the text.
[Figure 3]
Figure 3
Response of the YAP detector in the three operating modes of the APS. Symbols are the same as in Fig. 2[link], but note the axes extend to count rates about five times higher.
[Figure 4]
Figure 4
Response of the APD detector in the three operating modes of the APS. The dotted lines represent the non-extended dead-time model. The solid line corresponds to the isolated-bunch model in 24-bunch and hybrid singlet modes. The dash-dot line in 324-bunch mode is the extended dead-time model without PHA. The dashed line for hybrid singlet mode is the model of Lee & Mills (1992[Lee, W. K. & Mills, D. M. (1992). Rev. Sci. Instrum. 63, 1225-1228.]), which essentially overlaps the fit using the isolated-bunch model.

The Lee & Mills (1992[Lee, W. K. & Mills, D. M. (1992). Rev. Sci. Instrum. 63, 1225-1228.]) calculations for the hybrid singlet case are also shown in the three figures. Differences in peaking times and use of PHA leads to three distinct dead-time models for this fill pattern. The NaI detector is too slow to resolve the septets, so there are two `groups' of charge, with α1 = 16/102 (the singlet) and α2 = 86/102 (all the septets grouped together). Since the YAP and APD detectors can resolve the singlets, there are nine `groups' of charge in this case; the singlet is still α1 = 16/102, but the septets are αi = 10.75/102 for i = 2 to 9. The detectors also vary in their PHA properties; equation (8)[link] is appropriate for the NaI and YAP detectors with PHA, as is equation (9)[link] for the APD without PHA.

In Table 1[link] we present the measured dead-times and compare them with the most appropriate model as described in §2[link]. Uncertainties in the measured values of τ are typically 1–2%. In these calculations, τd is assumed to be the shaping time, i.e. twice the peaking time, for the scintillator detectors and the CFD output pulse duration for the APD. Comparing the measured and calculated values of τ for the scintillator detectors, it appears that the peaking time of 100 ns is an accurate value of τd for the YAP detector but 600 ns is an underestimate for the NaI detector. Using equation (5)[link] to compare the measurements in 24- and 324-bunch modes, a more accurate value of τd for the NaI detector is 720 ns. The isolated model does not fit the APD data in 324-bunch mode as well as it does in 24-bunch mode. Since the output of the CFD is longer than the bunch separation in 324-bunch mode, the bunches are not truly isolated. Thus, equation (6)[link] should be used rather than equation (7)[link], with n = 3. Ideally one would use faster electronics with an APD in a mode with such closely spaced bunches. When the output of the APD is gated to an individual bunch in 24-bunch mode (not shown), the isolated-bunch model applies just as well as for the ungated APD; the best-fit time constant was 3.61 µs, which compares with Tperiod as expected.

Table 1
Comparison of measured and calculated dead-times in the applicable model

Also shown are assumed values of τd for each detector. See text for details of the calculations.

    Effective τ (ns)  
Detector/τd (ns) Operating mode Measured Calculated Model used
NaI/600 324 731 602 Extended, with PHA
24 749 614 Extended, with PHA
Hybrid singlet 2186 2500 Lee & Mills, with PHA
YAP/100 324 104 102 Extended, with PHA
24 140 153 Extended, with PHA
Hybrid singlet 493 415 Lee & Mills, with PHA
APD/30 324 10.9 11.3 Extended, without PHA (n = 3)
24 142 153 Isolated, without PHA
Hybrid singlet 432 415 Lee & Mills, without PHA
†Data and model fit to equation (4)[link].
‡Data and model fit to equation (7)[link].

The equations for hybrid singlet mode are not very convenient to use in practice, since neither equation (8)[link] or equation (9)[link] is analytically solvable for NT. Therefore, we fit the results of these calculations as described above to equation (4)[link] and equation (7)[link], respectively, to obtain the calculated values of τ listed in Table 1[link]. The experimental data in this mode also have good fits to equations (4)[link] and (7)[link], as seen in Figs. 2[link]–4[link]. However, there are differences of ∼20% between the values of τ fit from the scintillator data and from the fits to the models, as listed in Table 1[link]. The discrepancies are probably due to a breakdown in one of the assumptions of the model. For the YAP detector, whose calculated τ underestimates the measured value, the limiting assumption is probably that the septets can be completely separated into non-interacting groups; τd is, in fact, slightly longer than the septet separation. For the NaI detector, in which the calculation overestimates τ, the error may arise from the assumption that bunches arriving within a group can be treated as arriving simultaneously. In this case, the relevant `group' is the whole train of septets, 500 ns long in all, whose duration is comparable with τd.

6. Conclusion

We have demonstrated the need for using the appropriate dead-time model with a proper value of the effective dead-time parameter in order to accurately correct high-count-rate data at a pulsed synchrotron source. The extended dead-time model with PHA was best for scintillator detectors, and the isolated-bunch model without PHA was usually appropriate for APDs. Specific models were, within their ranges of validity, fairly accurate at predicting the effective dead-time τ to use, but for best results τ should be measured whenever possible and its empirical value used.

The non-extended model is, perhaps, the more popular of dead-time corrections, with the advantage of a simple mathematical formula. Since all of the corrections discussed here are similar to first order, the non-extended model can be used for a detector operating slightly beyond its linear-response regime, but it is not generally trustworthy over a detector's full dynamic range. For example, if the non-extended model [equation (3)[link]] is used to correct data from a detector that obeys the non-extended model, the resulting systematic error will be under 1% as long as the true count rate is 0.15/τ or below.

These results can provide guidance in the development of future storage-ring operating modes. Certainly, detector dead-time is not the only factor to consider when designing a fill pattern, but it is one issue to address. For example, doubling the number of bunches in the standard mode of the APS from 24 to 48 would not improve the performance of the NaI or the YAP detector [as calculated by equation (5)[link]], although it should cut the APD's dead-time by one half. Fewer bunches will worsen τ for all three detectors and would be an additional disadvantage to those experiments which benefit from a quasi-continuous source. These results can also be used in designing asymmetric fill patterns, which must balance the need of experiments requiring special timing patterns with the need of reasonable, predictable dead-times for other experiments.

APPENDIX A

Solutions of the extended dead-time model

The equation for the extended dead-time model with pulse-height analysis, equation (4)[link], is called the Lambert W function (Corless et al., 1996[Corless, R. M., Gonnet, G. H., Hare, D. E. G., Jeffrey, D. J. & Knuth, D. E. (1996). Adv. Comput. Math. 5, 329-359.]). Various approximations exist for finding NT, the true count rate. One solution is found in a Taylor expansion (Corless et al., 1996[Corless, R. M., Gonnet, G. H., Hare, D. E. G., Jeffrey, D. J. & Knuth, D. E. (1996). Adv. Comput. Math. 5, 329-359.]),

[N_{\rm{T}}=N_{\rm{O}}{\sum\limits_{n\,=\,1}^\infty}{{\left(n\tau{N_{\rm{O}}}\right)^{n-1}}\over{n!}},\eqno(10)]

which will converge over the range of NO for which equation (4)[link] is applicable, that is, for NTτ < 1. Alternatively, one can use the following recursive algorithm based on the Newton–Raphson method to solve NO = NTexp(−τNT) for the true count rate NT. In this algorithm, NG is a guess for the true count rate at each step, NR is the resulting value of the observed count rate, and is the required accuracy, which can easily be 1 count s−1.

(i) Check observed rate. The maximum possible value for NO is (τe)−1, so the algorithm should stop with an error if NO > (τe)−1.

(ii) Initialize guess. Set NGNO. This is reasonable as an initial guess as NO is never greater than NT.

(iii) Calculate resultant rate. Compute NRNGexp(−τNG).

(iv) Finished? If (NONR) < , the guess is sufficiently close; set NTNG and exit the procedure.

(v) Improved guess. Use

[N_{\rm{G}}\,\,\leftarrow\,\,N_{\rm{G}}+\left[{{\exp\left(\tau{N_{\rm{G}}}\right)}\over{1-\tau{N_{\rm{G}}}}}\right]\left(N_{\rm{O}}-N_{\rm{R}}\right).]

The guess is increased by the resulting discrepancy, scaled by the slope of the function. Go back to step (iii).

This recursion relation converges much faster than equation (10)[link]. NT is approached from below, which ensures the algorithm does not arrive at the NTτ > 1 solution for a given NO (see Fig. 1[link]). As an example, for τ = 150 ns and = 1 count s−1, this recursion relation converges within nine iterations up to the maximum true count rate of 1/τ, whereas equation (10)[link] requires hundreds of terms for similar accuracy.

Footnotes

Current address: Department of Physics, DePaul University, Chicago, IL 60614, USA.

Acknowledgements

We thank John Quintana for helpful discussions, Steve Ross for supplying APDs, the APS Detector Pool for equipment loans, and the staff of beamline XOR/7ID for providing the beam time used for these measurements. Use of the APS was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357.

References

First citationArndt, U. W. (1978). J. Phys. E, 11, 671–673.  CrossRef CAS Web of Science
First citationBateman, J. E. (2000). J. Synchrotron Rad. 7, 307–312.  Web of Science CrossRef CAS IUCr Journals
First citationCorless, R. M., Gonnet, G. H., Hare, D. E. G., Jeffrey, D. J. & Knuth, D. E. (1996). Adv. Comput. Math. 5, 329–359.  CrossRef Web of Science
First citationCousins, C. S. G. (1994). J. Appl. Cryst. 27, 159–163.  CrossRef CAS Web of Science IUCr Journals
First citationIda, T. & Iwata, Y. (2005). J. Appl. Cryst. 38, 426–432.  Web of Science CrossRef CAS IUCr Journals
First citationKishimoto, S. (1997). Nucl. Instrum. Methods Phys. Res. A, 397, 343–353.  CrossRef CAS Web of Science
First citationKnoll, G. F. (2000). Radiation Detection and Measurement, 3rd ed. New York: John Wiley and Sons.
First citationLee, W. K. & Mills, D. M. (1992). Rev. Sci. Instrum. 63, 1225–1228.  CrossRef Web of Science
First citationQuintana, J. P. (1991). J. Appl. Cryst. 24, 261–262.  CrossRef Web of Science IUCr Journals

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds