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

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

The effect of beam time structure on counting detectors in SRS experiments

aRutherford Appleton Laboratory, Chilton, Didcot, Oxon OX11 0QX, UK
*Correspondence e-mail: j.e.bateman@rl.ac.uk

(Received 10 October 1999; accepted 22 May 2000)

Counting detector systems are increasingly used in X-ray experiments because of their attractive properties as regards linearity, large dynamic range and simple noise properties. In synchrotron radiation source (SRS) applications of X-ray detectors, counting rates are generally high enough to require dead-time correction. The time structure of an SRS beam interacts with the dead-time characteristics of the detector in a way that the simple stochastic dead-time models cannot always handle. This report generates analytical and Monte Carlo mathematical models which describe the rate performance of any given detector system when used with the typical beam structures encountered in an SRS.

1. Introduction

Counting detector systems are attractive for X-ray detection applications because of their inherent linearity, large potential dynamic range and well defined (statistical) noise characteristics. Their main limitation is the inevitable dead-time losses which appear at data rates comparable with the detector dead time. These dead-time losses are well understood when (as in the case of the detection of radioactive decay particles) the arrival of events is random in time (Knoll, 1989[Knoll, G. F. (1989). Radiation Detection and Measurement, ch. 17. New York: Wiley.]). Experiments on synchrotron radiation sources, on the other hand, combine very high data rates with a highly time-structured beam which complicates the simple stochastic model. The bunch structure of the electron beam in a synchrotron radiation source (SRS) produces an X-ray beam with the same time modulation. This time structure will (in general) interact with the intrinsic timing properties of any counting system used to detect the X-­rays and can often limit the rate performance of the system.

Aspects of the problem have been discussed by Arndt (1978[Arndt, U. W. (1978). J. Phys. E, 11, 671-673.]), Batterman (1980[Batterman, B. W. (1980). CHESS Technical Memorandum No. 6. CHESS, Ithaca, NY, USA.]), Lee & Mills (1992[Lee, W. K. & Mills, D. M. (1992). Rev. Sci. Instrum. 63(1), 1225-1228.]), Cousins et al. (1993[Cousins, C. S. G., Gerward, L., Laundy, D., Meads, R. E., Staun Olsen, J. & Sheldon, B. J. (1993). Nucl. Instrum. Methods Phys. Res. A, 324, 598-608.]) and Cousins (1994[Cousins, C. S. G. (1994). J. Appl. Cryst. 27, 159-163.]). Using mathematical models, this report investigates more general formulations for assessing the performance of any given type of detector system when used with some typical SRS beam structures. The modelling approach required to describe the inter­action between the beam structure and the detector system varies according to the relation between the counter dead time and the beam structure. It is convenient to consider the three distinct types of beam structure separately.

In the following the symbols Nin and Nout refer to the time-averaged rate of events (i.e. over a period long compared with any beam structure) which a detector system would deliver if it had no dead time (Nin) and the rate that it actually does deliver (Nout) in the context of a given beam structure. The two distinct types of dead-time model are denoted by the terms counting above a simple discriminator in which a single dead-time period is operative for each event and is not re-triggerable during that period, and pulse height analysis in which a clean spectrum demands rejection of pile-up events.

2. SRS beam structures

For the purposes of developing the models, we use the time structures used on the Daresbury Laboratory SRS. These are:

  • (i) Flat fill. The ring is full of evenly spaced bunches of electrons which give continuously an X-ray pulse of width 0.2 ns every 2 ns.

  • (ii) Single bunch. There is only one bunch in the ring which produces an X-ray pulse of width 0.2 ns every circuit of the ring, i.e. every 320 ns.

  • (iii) Gapped beam. In this case there is a gap in the bunch structure so that X-rays are delivered as in the flat-fill case for 200 ns of each beam circuit and are absent for the remaining 120 ns of the 320 ns period.

Different SRS machines will have a cycle time determined by the machine diameter (3 µs for the ESRF and ∼1 µs for Diamond) with a similar fine structure and various gap patterns.

3. Dead-time modelling with a flat-fill beam

No counting detector in routine use can respond to the fine structure of the beam. The fastest counter in general use is probably a scintillator/PMT device which can operate with a dead time of a few tens of nanoseconds. Recently, avalanche photodiodes have become available with the potential to yield a detector dead time of a few nano­seconds, which is, however, still well beyond the single-bunch width [see Kishimoto (1997[Kishimoto, S. (1997). Nucl. Instrum. Methods Phys. Res. A, 397, 343-353.]) for a discussion of this case]. The flat-fill case therefore reduces, as far as the counter system is concerned, to uniform delivery of X-rays in time. For this we can use the standard dead-time models. At any data capture rate accessible to the detector system there will be very much less than one detectable event per bunch so that the arrival of events will be effectively random in time.

3.1. Counting above a simple discriminator

In this case the response of the counting system (as a function of the input rate) is

[N_{\rm out}=N_{\rm in}/(1+N_{\rm in}\tau),\eqno(1)]

where Nout is the data capture rate, Nin is the random X-ray rate driving the system and τ is the counting-system dead time which is usually a combination of intrinsic detector dead time and electronic dead time (Knoll, 1989[Knoll, G. F. (1989). Radiation Detection and Measurement, ch. 17. New York: Wiley.]).

In this case the data capture rate is 50% of Nin at Nin = [1/\tau] and the data capture rate Nout asymptotes to 1/τ as Nin tends to infinity. Fig. 1[link] illustrates the behaviour of equation (1)[link] in terms of the natural parameters Ninτ and Noutτ.

[Figure 1]
Figure 1
The characteristic throughput curve predicted by formula (1)[link] for counting above a simple discriminator and by formula (2)[link] for the uncorrupted rate in a pulse-height spectrum. The natural parameters Ninτ and Noutτ are used where Nin is the event rate hitting the detector and Nout is the data capture rate. τ represents the detector/electronic dead time in the first case and the event occupancy time in the second case. This is the case of random incidence with no detectable beam structure.

3.2. Pulse height analysis

In the case of pulse height analysis we require a clean signal free of pile-up distortion. This requires that a dead-time period must always be free of another event. Poisson statistics show that in this case

[N_{\rm out}=N_{\rm in}\exp(-N_{\rm in}\tau).]

This function peaks at a value of Nin/e at Nin = 1/τ and declines to zero as Nin increases further. It is important to understand the significance of τ in this equation: it represents the sum of the occupancy times (the period in which the analogue pulse of the detector is significantly above threshold) of the pulse under consideration and any preceding pulse. If we wish to be consistent with the definition of τ in (1)[link], it is necessary to write this function as

[N_{\rm out}=N_{\rm in}\exp(-2N_{\rm in}\tau).\eqno(2)]

This is the actual counting rate that will be observed if a pile-up rejector is fitted to the system. If one is not fitted then the detected rate will be given by (1)[link], and (2)[link] will give the rate of `good' events in the PHA spectrum with Nout(1) − Nout(2) events in the `pile-up' spectrum. The word `significantly' in the above definition of the occupancy time sets the limits to the acceptable distortion in the pulse-height spectrum; usually the occupancy time is defined as the period that the pulse is >1% of its amplitude above the base line. Fig. 1[link] illustrates the form of equation (2)[link] graphically in terms of the natural parameters Ninτ and Noutτ.

In conclusion we note that the rate performance of the detector systems is entirely governed by their own dead time and is not affected by the beam structure.

4. Dead-time modelling with a single-bunch beam

In this case the situation is dominated by the fact that all the X-rays emitted during a bunch crossing are inseparable (in time) by any normal detector system. The important statistic is therefore the number of X-rays (M) incident on the detector per beam pulse of width 0.2 ns. The following two distinct cases must be considered.

4.1. Case 1: the detector dead time is shorter than the circulation period (τB)

In this situation the detector registers one count only, irrespective of the value of M, but is ready to trigger on the arrival of the next beam pulse. In order to evaluate the number (m) of counts recorded we must consider the Poisson distribution of the number of X-ray counts (x) actually produced for a mean number of M,

[P(x,M)=\exp(-M)M^x/x\,!\eqno(3)]

where

[\textstyle\sum\limits_{x}P(x,M)=1.]

The number of events at each occurrence number (x) is xP(x,M), i.e. exp(-M)Mx/(x-1)!, where

[\textstyle\sum\limits_{x}xP(x,M)=M.]

Summing the valid events we see that we can count the occurrences of single events, but double, treble etc. events only count as one,

[m=1P(1,M)+1P(2,M)+1P(3,M) +.\,.\,..]

Using equation (3)[link] we see that

[\eqalign{m&=P(1,M)+\textstyle\sum\limits_{x2}P(2,M)\cr&=P(1,M)+[1-P(0,M)-P(1,M)]\cr&=M\exp(-M)+1-\exp(-M)-M\exp(-M)\cr&=1-\exp(-M).}]

Since we obtain M X-ray counts, of which we detect m at every beam pulse, our observed counting rate is

[N_{\rm out}=[1-\exp(-M)]/\tau_B.\eqno(4)]

If Nin is the rate of events hitting the detector then

[N_{\rm in}=M/\tau_B]

and we can reformulate (4)[link] as

[N_{\rm out}=[1-\exp(-N_{\rm in}\tau_B)]/\tau_B.\eqno(5)]

As shown in Fig. 2[link], the data capture rate asymptotes to 1/τB as the beam intensity (Nin) is increased. This is the case in which we are simply counting above a simple discriminator.

[Figure 2]
Figure 2
The corresponding curves generated by equations (5)[link] and (5a)[link] for the case of a single-bunch beam combined with a detector dead time which is less than the beam circulation period τB. It will be noted that the detector dead time is irrelevant if it is less than τB. In this figure, and in most of the subsequent figures, the natural parameters NinτB and NoutτB are used.

When uncorrupted data (i.e. no pile-ups) is required then we are allowed to count only the term P(1, M) in the Poisson distribution. This gives a detected rate of

[\eqalignno{N_{\rm out}&=M\exp(-M)/\tau_B,\cr&=N_{\rm in}\tau_B\exp(-N_{\rm in}\tau_B)/\tau_B\cr&=N_{\rm in}\exp(-N_{\rm in}\tau_B).&(5a)}]

This is the same expression as (2)[link] above with the difference that the dead time is the beam circulation period and, because the beam is synchronous in its temporal pattern, the factor of two is not present. This curve is also shown in Fig. 2[link].

In conclusion we note that the system throughput curves are very similar to the flat-fill case except that the governing time constant is now the beam circulation time. Thus counting above a simple discriminator, the ultimate rate of any counter is just the beam circulation frequency. The `no pile-up' case is academic in this situation since it is impossible to make a pile-up rejector which can function within the 0.2 ns of the fine structure. However, the `pile up' curve (Fig. 2[link]) is of interest in that it shows the rate of uncorrupted events with the difference between the two curves showing the rate of `piled up' events.

4.2. Case 2: the detector dead time is longer than the beam circulation period

In this case any bunches which follow one with a detected event will not be seen by the detector until the dead time of the detector has elapsed. We can apply the models derived above as follows.

Counting above a simple discriminator, the counting system sees a data capture rate governed by relation (1)[link]. In this case the input rate is controlled by the statistics of the countable events (m), i.e. relation (4)[link]. Thus,

[N_{\rm out}={\left[1-\exp(-M)\right]/\tau_B\over1+\left[1-\exp(-M)\right]\tau/\tau_B}.]

By substituting Nin = M/τB, this becomes

[N_{\rm out}={{[1-\exp(-N_{\rm in}\tau_B)]/\tau_B}\over{1+[1-\exp(-N_{\rm in}\tau_B)]\tau/\tau_B}}.]

To make this formula represent the true behaviour of the counting system we must observe that in the single-bunch mode the beam loss in the counter dead time is digitized to the number of bunch crossings covered by τ. Thus the final formula for this case becomes

[N_{\rm out}={{[1-\exp(-N_{\rm in}\tau_B)]/\tau_B}\over{1+[1-\exp(-N_{\rm in}\tau_B)]n}},\eqno(6)]

where n is an integer defined by n = Int(τ/τB).

Checking the limiting cases of (6)[link] we see that if τ/τB < 1 (n = 0) then (6)[link] becomes (5)[link] as required. If τ >> τB and rates are high, then (6)[link] limits to 1/τ, and if rates are low then (6)[link] limits to (1)[link].

Fig. 3[link] shows the input/output curves given by (6)[link] for a few values of n. A rough rule of thumb emerges from Fig. 3[link] for very high rates (>1/τB), namely that Nout[\tau_B^{-1}(n+1)^{-1}]. As noted above, at low rates Nout tends to the standard formula (1)[link] for n >> 1.

[Figure 3]
Figure 3
The predictions of equation (6)[link] are plotted for various values of n [= Int(τ/τB)] for counting above a simple discriminator with dead time longer than the beam circulation time in the case of a single-bunch beam.

Similar reasoning can be applied to the pulse-height analysis case (i.e. no pile-ups allowed). Pile-up in this case arises in two distinct settings: pile-up within a bunch, the fraction of non-piled up events being

[N_{\rm out}/N_{\rm in}=\exp(-N_{\rm in}\tau_B),]

and pile-up in the detector system, where

[N_{\rm out}/N_{\rm in}=\exp(-2N_{\rm in}\tau).]

Combining these two fractions the overall rate of `good' events is

[\eqalignno{N_{\rm out}&=N_{\rm in}\exp(-N_{\rm in}\tau_B)\exp(-2N_{\rm out}\tau_Bn)\cr&=N_{\rm in}\exp[-N_{\rm in}\tau_B(2n + 1)],&(7)}]

where, as usual, n = Int(τ/τB). In the limiting case of n = 0 (no counter dead time) we see the formula for the bunch pile-up as expected, and at large detector dead times (2n >> 1) equation (7)[link] reduces to the pile-up formula for the detector system alone, as expected.

Fig. 4[link] shows the behaviour of equation (7)[link] for a range of values of n. In this figure the fraction Nout/Nin denotes the fraction of `clean' pulses in a pulse-height spectrum. Pile-up protection on the counter circuit can only remove the fraction generated by the counter time constant and the unresolvable pile-up within the bunch remains. This contribution to the pile-up will only become small at rates such that NinτB << 1.

[Figure 4]
Figure 4
The equivalent predictions from equation (7)[link] (pulse-height analysis) when the event occupancy is greater than the beam circulation time.

5. Dead-time modelling with a gapped beam

In the case of a gapped beam there are three interacting parameters describing the situation: the detector dead time (τ), the basic beam cycle time (τB) and the length of the `on' time of the beam (τon). This leads to a situation which is analytically untractable, and to which the simplest modelling approach is a Monte Carlo computer program. The Monte Carlo model also provides the opportunity to check the analytical formulae developed for the simpler cases.

In order to compute the Nout versus Nin curves for the case of counting above a simple discriminator, we focus our attention on a large number of beam cycles (∼1000) and generate events randomly in this interval building in the `on' period. We express τ and τon as fractions of τB and work in units of τB. In order to apply the dead time in the analogous fashion to an electronic system the chosen number of `input' pulses (a few thousand) are sorted into ascending order and only counted into the `output' total if they are more than one dead-time period from the preceding event. The input and output rates are expressed naturally as NinτB and NoutτB.

Fig. 5[link] shows the plots obtained from the program when τon is set to the value appropriate for the SRS (τon = 0.625τB) and the detector dead time is varied from τ = 0.1τB to 10τB (i.e. 32 nS to 3.2 µs corresponding to a scintillation detector at one end to a germanium detector at the other, for example). The input range of the counting rate is 0–6.25 MHz and the output range is 0–4.7 MHz.

[Figure 5]
Figure 5
The throughput characteristics predicted by the Monte Carlo model for a gapped beam with the `on' time (τon) set at the SRS value (0.625τB) for counting above a simple discriminator with a range of detector dead times τ expressed as fractions of τB.

The program can be modified in a simple manner to generate the corresponding curves for the case in which we require uncorrupted events (i.e. a clean pulse-height spectrum). The same logic is applied as above except that the dead-time period is extended forwards and backwards in time from the event under consideration.

Fig. 6[link] shows the Nout versus Nin curves generated in this way with the same parameter ranges as in Fig. 5[link]. As expected, these curves do not asymptote to a maximum rate but peak at values of Nin and Nout inversely related to τ, but clearly not in any simple way.

[Figure 6]
Figure 6
The corresponding model predictions for the `clean' pulse-height analysis rates with a range of occupancy times.

Computation of the system throughput curves in terms of these generalized parameters means that the curves in Figs. 5[link] and 6[link] can be used for any detector used with the given pattern of gapped beam. If a different τon is required then the curves must be recalculated. Any pattern of beam and gap can be modelled if it is required.

6. Comparison of the formulae with the Monte Carlo model

6.1. Counting above a simple discriminator

The simplest case of a flat-fill beam can be represented in the model by setting τon = 1.0. If τ = 1 (using units of τB) then the plot of NoutτB (y) versus NinτB (x) is just y = 1/(1+x) [equation (1)[link]]. As Fig. 7[link] shows, the agreement between the model and the function is perfect within the statistical noise of the Monte Carlo model.

[Figure 7]
Figure 7
A comparison of the predictions of the Monte Carlo model in a flat-fill beam with equation (1)[link] for the case of counting above a threshold. A range of dead times is shown.

The model can represent a single-bunch beam by setting τon = 0.001 (0.32 ns). Fig. 8[link] compares the output of the model with formulae (5)[link] and (6)[link]. Again the agreement is perfect within the statistical error. Because the bunch is very narrow, the dead-time curves switch abruptly with the value of n. In reality, with detector time jitter the transitions will be less precise.

[Figure 8]
Figure 8
A comparison of the predictions of the Monte Carlo model in a single-bunch beam for the output rate (counting above a simple discriminator) with equation (6)[link]. A range of dead times is shown.

6.2. Pulse-height analysis

Fig. 9[link] compares the rate of `clean' events in the case of a flat-fill beam with τ in the range 0.5–2, as predicted by equation (2)[link] and the Monte Carlo model.

[Figure 9]
Figure 9
A comparison of the Monte Carlo model in a flat-fill beam (`clean' pulse-height analysis rate) with equation (2)[link]. A range of occupancy times is shown.

Fig. 10[link] compares the predictions of formula (7)[link] and the model in the case of a single-bunch beam when the detector occupancy time increases from less than τB to greater than τB.

[Figure 10]
Figure 10
A comparison of the Monte Carlo model in a single-bunch beam (`clean' pulse-height analysis rate) with equation (7)[link]. A range of occupancy times is shown.

In all cases the agreement is perfect within the statistical noise of the Monte Carlo process. Since the Monte Carlo model simply mimics the corresponding electronic function, these results give confidence in the mathematical analysis (and vice versa in the case of the simplest formulae).

7. Conclusions

The models developed in this report permit detailed exploration of the performance of counting systems in any of the beam time structures generally met. As expected, counting above a simple discriminator leads to an asymptotic data capture rate and the rate of unpiled-up events peaks and thereafter declines. The main point to emerge is that in the flat-fill case the detector dead time (τ) is the controlling parameter and in the single-bunch case the beam circulation time (τB) is the controlling parameter. The result of the latter effect is that a detector with a dead time significantly less than τB will be limited to a rate of [\sim1/\tau_B]. For such a detector the flat-fill case is more advantageous. For very slow detectors (τ >> τB) the difference is small.

There is one fast detector system which can significantly reduce the losses associated with single-bunch operation. Using a gas microstrip detector (GMD) with the electron drift direction normal to the beam direction (Bateman et al., 1998[Bateman, J. E., Connolly, J. F., Derbyshire, G. E., Marsh, A. S., Stephenson, R., Farrow, R. C., Helsby, W. I., Dobson, B. R, Mutikainen, R. & Suni, I. (1998). Report RAL-TR-1998-073. Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, UK. (http://www-dienst.rl.ac.uk/library/1998/tr/raltr-1998073.pdf )], 2000[Bateman, J. E., Connolly, J. F., Derbyshire, G. E., Marsh, A. S., Stephenson, R., Farrow, R. C., Helsby, W. I., Dobson, B. R., Mutikainen, R. & Suni, I. (2000). Nucl. Instrum. Methods Phys. Res. A. To be published.]) one can use the drift time of the X-ray-generated electron clouds to randomize the times of arrival at the detector anodes. The 320 ns period of the SRS corresponds to a drift distance of 15 mm, which is quite comfortable. For longer periods of up to several microseconds it is possible to choose slower drift conditions without compromising the basic speed of the anode pulse. Operating in this way the GMD can randomize the events in the beam spike and detect them with a limit set by its own dead time (∼100 ns). On the SRS the gain in limiting rate is of the order of ×3 but, on the larger machines, factors of >10 could be achieved in single-bunch mode.

The performance of any detector system in gapped-beam mode is obviously intermediate between the single-bunch and flat-fill cases. The interaction between the parameters is complex and the practical method of solving it is to use the Monte Carlo model with the parameters of the situation appropriately inserted. The simple model used is approximate but adequate to elucidate the key features of the situation.

The models show clearly the great difficulty in obtaining clean pulse-height spectra (and clean discrimination) at high rates. As Figs. 6[link] and 10[link] show, if the occupancy time is larger than τB in single-bunch or gapped beams, the occupancy time of the detector must be less than τB to achieve a useful fraction of `clean' data. Similarly, in a flat-fill beam (Fig. 9[link]) the value of τ is critical to achieving `clean' data at high rates.

References

First citationArndt, U. W. (1978). J. Phys. E, 11, 671–673.  CrossRef CAS Web of Science Google Scholar
First citationBateman, J. E., Connolly, J. F., Derbyshire, G. E., Marsh, A. S., Stephenson, R., Farrow, R. C., Helsby, W. I., Dobson, B. R, Mutikainen, R. & Suni, I. (1998). Report RAL-TR-1998–073. Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, UK. (http://www-dienst.rl.ac.uk/library/1998/tr/raltr-1998073.pdfGoogle Scholar
First citationBateman, J. E., Connolly, J. F., Derbyshire, G. E., Marsh, A. S., Stephenson, R., Farrow, R. C., Helsby, W. I., Dobson, B. R., Mutikainen, R. & Suni, I. (2000). Nucl. Instrum. Methods Phys. Res. A. To be published.  Google Scholar
First citationBatterman, B. W. (1980). CHESS Technical Memorandum No. 6. CHESS, Ithaca, NY, USA.  Google Scholar
First citationCousins, C. S. G. (1994). J. Appl. Cryst. 27, 159–163.  CrossRef CAS Web of Science IUCr Journals Google Scholar
First citationCousins, C. S. G., Gerward, L., Laundy, D., Meads, R. E., Staun Olsen, J. & Sheldon, B. J. (1993). Nucl. Instrum. Methods Phys. Res. A, 324, 598–608.  CrossRef Web of Science Google Scholar
First citationKishimoto, S. (1997). Nucl. Instrum. Methods Phys. Res. A, 397, 343–353.  CrossRef CAS Web of Science Google Scholar
First citationKnoll, G. F. (1989). Radiation Detection and Measurement, ch. 17. New York: Wiley.  Google Scholar
First citationLee, W. K. & Mills, D. M. (1992). Rev. Sci. Instrum. 63(1), 1225–1228.  CrossRef Web of Science Google Scholar

© 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