Patent Grant
7345285
This application claims priority to French Application No. 04 50627, filed on Mar. 31, 2004, entitled: “Spectra Acquisition System with Threshold Adaptation” by Guillaume MONTEMONT, Patrice OUVRIE-BUFFET, and Loick VERGER and was not published in English.
The invention concerns the field of spectrometry, particularly X or gamma, by means of a semi-conductor detector.
Such detectors often have charge carrier collection problems.
An incident radiation spectrometry makes it possible to determine the quantity of charge carriers created by the interaction of a photon, for example X or gamma.
To this end, one applies an electric field to the detector material by means of electrodes in such a way as to drain these charges and induce an electric signal the amplitude of which is proportional to their quantity. Unfortunately, the measurement of this quantity of carriers is made difficult by the imperfect transport properties of the semi-conductor: indeed, part of the charges does not reach the electrodes. This is the problem of incomplete collection.
Furthermore, particularly in the case of CdZnTe, transport properties of the holes are particularly poor whereas those of the electrons are good.
One then measures a signal that corresponds to the transit of the electrons in the material.
One of the methods proposed to correct this incomplete collection, which is due to poor transport of the holes, is to carry out an electronic processing of the measured signal.
In addition to the amplitude measurement measuring the total collected charge, one or several other parameters of the form of the obtained electric signal are measured (for example, its rise time). This makes it in particular possible to go back to the site of interaction of the photon in the semi-conductor medium.
Thanks to a calibration of the collection efficiency as a function of the interaction site, one may then determine the charge actually deposited by the photon.
This method is the bi-parametric “electrons” method disclosed in document FR-2738919.
The rise time may be determined in different ways.
A first method consists in measuring the time spent between the charge signal going through a triggering threshold (performed by a comparator) determining the start of the pulse and its arrival at its maximum amplitude (measured by a peak detector).
Another method consists in working on the current signal (obtained by means of a differentiator circuit). One then measures the time spent between the current passing above a threshold value and its return below said threshold. This method has the advantage of being precise and easy to carry out since it uses a comparator and a time to amplitude converter.
However, said methods have limitations, particularly with regard to the linearity of the measurements on an important signal dynamic. Moreover, one can configure them so as to minimise the noise uniquely around certain signal values.
In this figure, Te is the transit time of the electrons for a given interaction site: T2 is the time measured for a high energy and T1 is the time measured in the case of a weaker energy pulse.
One therefore obtains measured times which are different to each other and different to the transit time of the electrons, which is normally unique.
Moreover, in terms of noise, there will be an optimum threshold position (in general, at the point where the slope is at a maximum). However, this position can only be met for a certain amplitude value.
Indeed, there is a compromise to be made between dynamic and noise: filtering limits noise but reduces stability of performances over a large dynamic.
A second limitation of conventional methods is the impossibility of processing complex pulses, due for example to multiple interactions.
A double interaction ID is represented in
Two incident photons of different energies give the measured amplitude/time pairing.
The conventional bi-parametric method therefore fails to identify the energies.
It would be interesting to identify these cases of multiple pulses in order to avoid degrading the spectra.
Moreover, the possibility of separating the contributions of the 2 electron clouds is interesting for recognising simple interactions from multiple interactions and to more easily identify the peaks on a Compton continuum.
Indeed, in the energy range 100 keV -1000 keV, the 2 types of important interactions are the photoelectric effect (the photon is completely absorbed) and the Compton Effect (diffusion with deposition of a fraction of energy).
The simple interactions make it possible to obtain a better energy resolution whereas the double interactions have high probabilities of corresponding to reabsorbed Compton photons and therefore belonging to a complete deposition peak.
The combination of these two items of information offers the possibility of reducing the Compton continuum while at the same time improving resolution.
Finally, a third limitation is as follows. Not all detectors provide a current having a “flat” shape. For the geometries commonly used, the current increases when the electrons approach the anode. It is therefore necessary to apply an end of pulse threshold higher than the start of pulse threshold in order to optimise the performance. The measurement is therefore not facilitated.
The invention first concerns a method for exploiting or analyzing or processing a signal, particularly a current, provided by a detector, in which:
According to a specific embodiment, a method according to the invention further comprises the following steps:
The invention further concerns a method for exploiting or analyzing or processing a signal, particularly a current, provided by a detector, in which:
The invention further concerns a method for exploiting or analyzing or processing a signal, particularly a current, provided by a detector, comprising a first component superimposed on a second component, in which:
Whatever the envisaged embodiment, a method according to the invention involves a re-adaptation, by means of the previous estimation of the charge and the transit time, of the threshold value used for measuring the time and the amplitude.
The invention therefore further concerns a spectrometry method with adaptive bi-parametric correction, preferably by means of digital means.
The new threshold values may be calculated from previously memorised reference data.
Reiteration may take place until at least one difference between two successive values of at least one of the charges determined and/or of the time spent above at least one of the thresholds and/or one of the time intervals is less than a predetermined value.
The steps of determining charges may be performed through integration.
The signal is advantageously filtered before integration.
The new threshold values may be determined by an analytical method, or by an experimental calibration or by modelling.
The signal is generated for example by a semi-conductor detector, for example a CdTe or CdZnTe semi-conductor detector.
The invention further concerns a device for exploiting or processing a signal provided by a detector, comprising:
Such a device may further comprise:
According to another aspect, the invention further concerns a device for exploiting or processing or analyzing a signal provided by a detector, comprising:
It further concerns a device for exploiting or processing or analyzing a signal provided by a detector and comprising a first component superimposed on a second component, said device comprising:
The means for determining the charge(s) may comprise one or several integrators.
Such a device according to one of the above embodiments may further comprise means for filtering the signal or the signals before integration.
Said means for determining new threshold values can be programmed to implement an analytical method, or an experimental calibration or a modelling.
A detection device according to the invention comprises a detector and a device such as described here-above.
A first embodiment of the invention is illustrated in
The system comprises a semi-conductor type detector 1, for example in the field of X or gamma spectrometry.
A charge preamplifier 2 provides an electrical charge Q (t) to a differentiator 4, from which a current pulse I (t) is obtained on the basis of charge Q (t).
It further comprises comparator type thresholding means 16.
Said thresholding means thus make it possible to identify the instants or the times between which the signal is above a given threshold value. They produce a logic signal, for example a pulse or a gate signal of the type illustrated in
This device makes it possible to measure the time or the duration the current remains above a certain threshold S1, or is superior to said threshold, which enables a measurement or an estimation of the transit time T (
Time to amplitude converter means 14 perform this measurement or this estimation of time from the logic signal provided by the thresholding means 16.
An analogue embodiment of means 14 integrate a constant current during the active state of the logic signal.
A digital embodiment of means 14 may for example comprise a clock pulse counter, the frequency of which is compatible with the desired precision.
A commutated integrator 8 integrates the current over this period and thus provides a signal relative to the corresponding electrical charge. This charge is the integral of the signal over the time T.
The informations relating to charge and time (provided by means 8 and 14) are provided to calculation means 10, to calculate or determine a new threshold S2 which will replace threshold S1 in means 16, and which will be better suited to the pulse to be analysed.
Digital threshold calculation means 10 may for example comprise a microprocessor or a programmable arithmetic unit.
One uses, to measure the amplitude, a derivation of the charge signal Q (t), then re-integrated (by the means 8) over the measured time T of the pulse, after thresholding.
The system iterates the process a sufficient number of times to attain a desired precision: at each iteration, one uses a new threshold value. The signal has been previously memorised in the memorisation means 7.
Preferably, the system further comprises a low pass shaping filter 6.
The characteristics of said low pass filter 6 are determined according to the measurement noise spectral density (generated for the most part by detector 1 and pre-amplification electronics 2). Said noise spectral density may either be measured or estimated theoretically.
The impulse response of this filter 6 is, according to one embodiment, calculated so as to compensate the noise colouration.
A finite impulse response filter 6 of the digital type using, for example, multi-rate techniques may also be suitable and is easier to achieve.
Moreover, in a feedback structure as that of
One may cite different methods for determining how means 10 can calculate the new threshold.
A first method is known as the analytical method. One estimates the charge Q of the pulse and its time or duration Assuming the current pulse has a square shape, the average current I=Q/T gives a good estimation of the upper level of the pulse. For a filter having a symmetric pulse response, placing the threshold at I/2=Q/2T therefore allows one to lock on the maximum slope. One thus refines the measurement.
A second method is known as the “Experimental calibration method”.
By means of a generator, one sends into the input of the system pulses simulating a detector, with a known charge and time or duration (Q,T) pairing. One then scans the domain (Q,T) and, for each adjustment, one determines the threshold value that enables a correct measurement to be obtained (the Q and T values measured using said threshold are in conformance with the real Q and T of the pulse). One is then able to “tabulate” the function threshold=f(Q,T) over the whole domain. A measurement of Qm and Tm on a real signal, provided by means 8,14, then enables a threshold S=f(Qm, Tm) to be determined.
A third method is known as the “Modelling Method”. This approach is similar but uses a noised pulse model that one derivates and filters (the impulse response of the filter is known).
One then determines, by calculation, the threshold position that enables the best estimation of Q and T, for example by means of the MAP (Maximum A Posteriori) criterion. It involves a probabilistic approach that consists in placing the threshold at a level such that one obtains an estimation of Q and of T with a maximum probability.
Whatever method is used, one can digitize the signal either just before the filtering, or as of the output of the preamplifier 2.
As criterion for convergence, and therefore for stopping the iterations, one may choose a difference, less than a predetermined value, between two successive thresholds Si and Si+1 and/or two successive quantities of charges Qi and Qi+1.
The invention further concerns a device and a spectrometry method with multi-parametric correction for recognising multiple interactions.
Such a device is illustrated in
It comprises elements similar to those of
But here several thresholdings are applied in order to separate the different components of the pulses.
This system makes it possible, in particular, to process the case of two components (
To break down the pulse, one may for example apply 3 thresholds: a first threshold (S1) to detect the instant of interaction (t0), corresponding to the start or the beginning of the pulse; a second threshold (S2) to detect the end of component n° I (t1) and a third threshold (S3) for the end of component n° II (t2).
The signal is then integrated over two periods (from the start to the end of pulse n° I and from the start to the end of pulse n° II) so as to evaluate the respective contributions of Q1 and Q2.
The integration over T2 therefore gives for example Q1+Q2 whereas the integration over T1 will give Q1+Q2.T1/T2.
By means of these 2 data or equations, one obtains the 2 unknown values Q1 and Q2.
Said three thresholds S1, S2, S3 are respectively associated with thresholding means 16, 26, 36.
The first two thresholds make it possible to determine the start or the beginning and the end of the first component of the signal.
Thresholds S1 and S3 make it possible to determine the start or the beginning and the end of the second component of the signal.
The 3 thresholds S1, S2 and S3 are linked and adjusted by means 10 after calculation of charges Q1 and Q2 and intervals T1 and T2 associated with these components.
Here again, means 10 may be programmed according to one of the three above-mentioned methods (analytical method, or by experimental calibration or by modelling).
It is also possible to process a signal such as that of
In this figure, reference 46 designates thresholding means, reference 48 an integrator, and reference 44 time to amplitude converter means.
First and second thresholds are readjusted by means 10 from the moment data relative to charges and times have been calculated, and these readjusted values are re-injected in the thresholding means 16, 46.
Iterations are repeated until a predetermined convergence threshold is attained. Here again, one can retain, as convergence criterion and therefore for stopping the iterations, a difference, less than a predetermined value, between two successive thresholds Si and Si+1 and/or two successive quantities of charges Qi and Qi+1.
The invention may also be adapted to the different geometries of detectors, particularly detectors of non-planar geometry, giving signals of more complex form.
In this case, the start and the end of the pulse are processed differently, with a start threshold value and an end threshold value that are different to each other.
Preferentially the pulse is first filtered in order to obtain a signal with steeper sides, so as to optimise the precision of the time measurement.
As illustrated in
Thresholds S1 and S2 define a time T between start or beginning and end of the pulse.
The charge Q corresponding to the difference between the integrals of the filtered curve (from 0 to T2) and the real curve (from 0 to T1) is calculated (see formula on the following page).
The values of Q and T are entered in calculation means 110 to determine new start and end pulse threshold S′1 and S′2. There ensues a new calculation of charge Q′, a new time interval T′, which in turn are entered in calculation means 10.
The iterations are repeated until a predetermined convergence threshold is attained.
A device for the implementation of this method is illustrated in
In concrete terms, the filtering may be determined in the following manner.
One first establishes a model of the shape of the current detector at the moment of the transition (which may be achieved either in digital form or in analytical form, by a differential equation for example).
One then determines a response of the linear filter transforming this shape into a pulse or a gate signal, which consists in carrying out a mathematical operation of deconvolution.
A commutated differential integrator 118 realises the difference of the integrals of the 2 input signals between the two start and end dates.
This may be carried out analogically in the manner illustrated in
The calculation made corresponds to the formula:
Whatever the envisaged embodiment, the invention results in the determination of a pairing (charge Q, time T) that makes it possible to determine the energy of the incident photon thanks to a calibration curve established by means of an energy calibrated source. The resulting value of Q is in principle proportional to the energy for a given time T.
A method involving the determination of the energy of the photons as a function of the measured charge and rise time is described in FR-2 738 919.
Calculation means 10, 110 may memorise a correspondence table that makes it possible, during subsequent measurements, to know the threshold to apply to return, by iteration, to the charge representative of the interaction phenomenon.
For a charge Q measured and a time T obtained, for example by comparison of the signal with a threshold S, one searches in the table or in means 10, 110 the pairing (Q, T) and one obtains the corresponding threshold S′, itself re-injected as threshold value in the comparator etc.
| Number | Date | Country | Kind |
|---|---|---|---|
| 04 50627 | Mar 2004 | FR | national |
| Number | Name | Date | Kind |
|---|---|---|---|
| 4543530 | Kovach | Sep 1985 | A |
| 5821538 | De Antoni et al. | Oct 1998 | A |
| 5854489 | Verger et al. | Dec 1998 | A |
| 6211664 | Bonnefoy et al. | Apr 2001 | B1 |
| 20030226973 | Beusch | Dec 2003 | A1 |
| Number | Date | Country |
|---|---|---|
| 2 738 919 | Sep 1995 | FR |
| Number | Date | Country | |
|---|---|---|---|
| 20050230632 A1 | Oct 2005 | US |
