EP2656111A1 - Procede de mesure du taux de comptage de coïncidences, utilisant une conversion temps-numerique et une methode de temps mort reconductible avec mesure du temps actif - Google Patents

Procede de mesure du taux de comptage de coïncidences, utilisant une conversion temps-numerique et une methode de temps mort reconductible avec mesure du temps actif

Info

Publication number
EP2656111A1
EP2656111A1 EP11804670.5A EP11804670A EP2656111A1 EP 2656111 A1 EP2656111 A1 EP 2656111A1 EP 11804670 A EP11804670 A EP 11804670A EP 2656111 A1 EP2656111 A1 EP 2656111A1
Authority
EP
European Patent Office
Prior art keywords
time
measurement
detectors
detector
events
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Withdrawn
Application number
EP11804670.5A
Other languages
German (de)
English (en)
Inventor
Christophe BOBIN
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Original Assignee
Commissariat a lEnergie Atomique CEA
Commissariat a lEnergie Atomique et aux Energies Alternatives CEA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Commissariat a lEnergie Atomique CEA, Commissariat a lEnergie Atomique et aux Energies Alternatives CEA filed Critical Commissariat a lEnergie Atomique CEA
Publication of EP2656111A1 publication Critical patent/EP2656111A1/fr
Withdrawn legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R29/00—Arrangements for measuring or indicating electric quantities not covered by groups G01R19/00 - G01R27/00
    • G01R29/02—Measuring characteristics of individual pulses, e.g. deviation from pulse flatness, rise time or duration
    • G01R29/027—Indicating that a pulse characteristic is either above or below a predetermined value or within or beyond a predetermined range of values
    • G01R29/0273—Indicating that a pulse characteristic is either above or below a predetermined value or within or beyond a predetermined range of values the pulse characteristic being duration, i.e. width (indicating that frequency of pulses is above or below a certain limit)
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/17—Circuit arrangements not adapted to a particular type of detector
    • G01T1/171—Compensation of dead-time counting losses
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/17—Circuit arrangements not adapted to a particular type of detector
    • G01T1/172—Circuit arrangements not adapted to a particular type of detector with coincidence circuit arrangements

Definitions

  • the present invention relates to a method for measuring the counting rate of coincident events between a plurality of radiation detectors that operate in parallel.
  • radiation is meant photons or particles.
  • the measurement is made by recording the distribution of time intervals which are defined by the duration between a trigger signal, or "start” signal, provided by a measurement channel, stop signal, or “stop” signal, from another measuring channel.
  • a trigger signal or "start” signal
  • stop signal or “stop” signal
  • a time-to-digital conversion is combined, more precisely a conversion in digital form of the recorded time intervals, with a protection against the "start" and "stop” signal trains that can cause time-biased measurements with known measurement techniques.
  • the above-mentioned protection uses an extendible dead- time method; this dead time is common to the different ways; it is built from the signals "start” and “stop”; this method of renewable dead time is combined with a real-time measurement of active time (in English, live-time).
  • the application of the active time method makes it possible to measure the effective counting periods to determine the coincidence rate between the different channels of the radiation detection system.
  • the present invention applies in particular to nuclear instruments that implement the measurement of coincidences between radiation detectors.
  • the invention makes it possible, for example, to determine the activity of a radionuclide.
  • the term "dead time” refers to a period of paralysis of a measurement system. This period follows the detection of a pulse in the system. During this period, any new pulse can not be processed correctly for the acquisition of information such as a count or an amplitude for example.
  • the origin of the paralysis depends on the system used; it may be a correlated pulse saturation or post-pulse (in English, afterpulses).
  • the expression "method of the renewable dead time” designates a method for managing the paralysis of the detection system under consideration. It consists in preventing, during a predefined duration, any processing of a pulse to extract information, following the detection of this pulse.
  • any new signal arriving during this period is only used to extend it by the same predefined duration.
  • the system becomes free or active again when no new pulse has been detected.
  • the active time method consists of measuring the cumulative active time between the dead time periods by counting the pulses of a clock. We can then consider that the real time of a measurement is sampled. In the field of nuclear instrumentation, it is common to implement a counting method of measuring the coincidence rate between two or more radiation detectors that operate in parallel.
  • the first is based on the use of a logic circuit that produces a coincidence resolution time for each channel.
  • Coincidence counting is performed by constructing a recovery period between logical signals.
  • the second consists in measuring time intervals between the channels to obtain the recording of a temporal distribution.
  • the time slot measurement sequence begins with a "start” signal that is provided by one channel. And the time intervals are measured by the "stop" signals that appear on the other channels.
  • This second method makes it possible to preserve the information given by the temporal fluctuations of the detection system for counting coincidences. From the recording of histograms of the durations of occurrence between the channels, an off-line processing allows to adjust the delay and the length of the coincidence resolution time so as to they are adapted to the detection system.
  • the recording of a time interval requires two steps: the time interval between the signals "Start” and “stop” are first converted into an amplitude, and then the amplitude is exploited by a multichannel analyzer to record the time distribution.
  • the disadvantage of such a method is the introduction of a period of paralysis. This is due in particular to the amplitude conversion process which requires a charge period of a capacitor.
  • the multichannel analyzer can also contribute significantly to the paralysis of the measurement system.
  • a time-to-digital converter performs the time measurement from the counting of the pulses of a clock whose frequency defines the optimal temporal resolution of the device. More sophisticated systems can improve this resolution through the use of interpolation methods or the Vernier method.
  • the known systems for measuring time intervals have a drawback: for these systems, the problem of paralysis does not take into account the detection means associated with them.
  • This problem is particularly significant when, in the radiation detector, the same event generates correlated post-pulses at the input of the same channel.
  • post-pulses may also be due to metastable states that are characteristic of certain radionuclides.
  • An object of the invention is to overcome this disadvantage.
  • This algorithm combines a time-to-digital conversion, intended for the measurement of time intervals, with a transposition of a renewable time-out method. This dead time is common to all channels.
  • the algorithm is implemented in a FPGA-type programmable component, or in-field programmable programmable gate array, for real-time processing.
  • the actual measurement time is determined by a transposition of the active time method.
  • sampling is performed by means of the clock included in the programmable component.
  • a transposition in digital form of the module described in document [4] in the case of a coincidence measurement using a recovery time, is known from the following document to which we refer. report to :
  • a particularity of the invention lies in the fact that it does not include any particular channel that would be dedicated to triggering the measurement of time intervals: the beginning of a measurement is triggered by any one of the detectors. Consequently, the invention can be adapted to a symmetrical detection system, for example a system intended for the application of the RCTD method (English, TDCR) or the Triple-to-Double Coincidence Method (Triple, Triple). to Double Coincidence Ratio).
  • the algorithm manages the arrival times of the pulses in the other channels to establish temporal histograms of the multiple coincidences between the channels (double coincidences, triple coincidences, etc.). .).
  • the correlated post-pulses in the different channels are used for the renewal of the dead time.
  • the acquired information namely the histograms corresponding to the multiple coincidences between the channels, is recorded in an acquisition computer; the measured active time is recorded there.
  • an algorithm is applied in real time to a detection system comprising at least two channels, whereas only one channel is considered in the document. And the measurement of the active time is applied directly by sampling the periods outside the dead time using a clock.
  • an object of the present invention is the measurement of a counting rate of coincidences between several detectors.
  • the present invention essentially solves a problem of processing correlated post-pulses, namely post-pulses which are generated in the detectors used or which may result from metastable states of certain radionuclides. This is the reason why the dead time is common to the different channels of the measuring system. These are not just periods of paralysis caused by the duration of discrimination.
  • the invention is closer to the almost direct digital transposition of the analog module MAC3 which has already been done (see documents [4] and [5]).
  • the subject of the present invention is a method for measuring the counting rate of coincident events between N radiation detectors which operate in parallel and are respectively associated with N detection channels, where N is an integer at least equal to to 2, each detector being able to send an electrical pulse to the detection channel with which it is associated when an event occurs in this detector, in which:
  • the counting rate of the coincident events is measured using the recorded temporal distributions.
  • the active time is measured in real time.
  • the dead time is preferably common to the N detection channels.
  • the measurement of the count rate is triggered by one of the N detection channels when an event occurs in the detector associated with it, and
  • an algorithm is implemented which establishes temporal histograms corresponding to the multiple coincidences between the N detection channels, that is to say the coincident events between P detectors among the N radiation detectors, where P traverses the set of integers ranging from 2 to N, from the arrival times of the pulses in the N1 other detection channels.
  • the detectors may be identical to each other.
  • N is equal to 3
  • the detectors are photomultipliers and the method is used to implement the triple-to-double coincidence ratio method.
  • the detectors may also not be identical to each other.
  • N is equal to 2
  • the detectors are respectively a detector of gamma photons and an electron detector and the method is used to implement the beta-gamma coincidence method.
  • FIG 1 is a schematic view of a measurement chain which is intended to determine the activity of a radionuclide by the RCTD method and in which is implemented an example of the method, object of the invention,
  • FIG. 2 is a flow chart of an algorithm which is used in this example and which deals in parallel with time measurements and the management of the renewable dead time, and
  • FIG. 3 is a timing diagram relating to this algorithm.
  • the example that will be described relates to an application of the invention to the RCTD method which is commonly used for measuring an activity by liquid scintillation counting (in English, liquid scintillation counting).
  • the implementation of this method makes it possible to highlight the advantages of the invention in the case where the arrival of the signals "start” and “stop” is random because of the existence of correlated post-pulses which can lead to at confusing times in the measurement of time intervals.
  • the measurement chain which will be discussed later, comprises three counting channels, or detection channels, which are identical; and each of these pathways begins with a photomultiplier.
  • the detection is done symmetrically so that each of the channels is capable of triggering a measurement sequence and renewable timeout. There is therefore no specific path dedicated to triggering a time slot measurement process.
  • the other two channels are used to implement histograms corresponding to second and third pulses (for counting double and triple coincidences).
  • a measuring device using the RCTD method was made using a commercially available digital card, namely an Altré ® development kit, equipped with a Stratix ® III FPGA chip.
  • the time per channel is defined as a multiple of the minimum temporal resolution.
  • the total measurement time and the minimum dead time are defined by the user according to a number of clock ticks.
  • connection between the digital part of the measurement chain and the acquisition computer is achieved through an Ethernet connection.
  • Figure 1 is a schematic view of the measurement chain in which is implemented an example of the method, object of the invention.
  • a bottle 2 is introduced containing a mixture of a liquid scintillator and a radioactive source whose activity is to be determined by the RCTD method.
  • three photomultipliers 4, 6, 8 are used which are symmetrically arranged around the vial 2, at 120 ° from each other. Following radioactive decay, ionizing radiation is emitted. This results in an energy deposit in the liquid scintillator. This leads to the emission of light photons which are distributed randomly between the three photomultipliers.
  • each photomultiplier converts an incident light photon into a photoelectron. This conversion is performed by means of a photocathode located at the entrance of the photomultiplier and depends on the quantum yield of the photocathode.
  • the photomultiplier comprises a chain of dynodes.
  • dynodes In the latter is a process of multiplication of photoelectrons that are produced at the photocathode. This gives a current large enough to be transformed into a voltage pulse that is exploitable by a fast amplifier.
  • This pulse is sent to a detection path that connects the photomultiplier to the fast amplifier.
  • the detection channels 10, 12, and 14 which respectively connect the photomultipliers 4, 6, 8 to the fast amplifiers 16, 18, 20 have been represented.
  • These amplifiers 16, 18, 20 are respectively associated with analog devices 22, 24, 26, namely DFCs (in English, CFD), that is to say constant fraction discriminators (in English, constant fraction discriminators). .
  • DFCs in English, CFD
  • constant fraction discriminators in English, constant fraction discriminators
  • the measurement system also comprises a digital device 28 which consists of an FPGA in the example described.
  • a digital device 28 which consists of an FPGA in the example described.
  • the time-to-digital conversion and the protection discussed above have been programmed and will be discussed later.
  • the set of DFCs makes it possible to switch from the analog part of the detection system to the digital time-to-digital conversion device 28.
  • Each DFC produces a logic signal that is directly exploitable by this device 28 and has the advantage of reducing the temporal fluctuations. that is generally encountered in the case of conventional threshold discrimination.
  • the three DFCs which respectively equip the detection channels, provide logic impulses. These pulses materialize the "start" and "stop” signals that arrive at the input of the digital device 28.
  • the measurement chain also includes an acquisition computer 30 which processes the time histograms provided by the digital device 28 and determines the counting rates sought.
  • this computer is connected to the digital device 28 via an Ethernet link.
  • the computer 30 is provided with a device 32 for displaying the results of the measurements.
  • the active time is incremented with each clock stroke.
  • the active time is incremented with each clock stroke.
  • the active time is incremented with each clock stroke.
  • the dead time counter is set to a predefined value.
  • the channel number of the double coincidence histogram is 1 and the channel number of the triple coincidence histogram is 1. In 48, the channel number of the triple coincidence histogram is 1.
  • the dead time counter is decremented by 1.
  • the channel number of the histogram of the double coincidences is increased by 1 and the channel number of the histogram of triple coincidences is increased by 1.
  • the channel number of the histogram of triple coincidences is increased by 1 and there is a new clock tick.
  • the channel No. 1 is the first channel on which a signal is detected from the corresponding PMT, namely one of the PMTs A, B and C with the notations of FIG. 3;
  • PMT A is one of the three PMTs of Figure 1
  • PMT B is one of the other two PMTs of Figure 1
  • PMT C is the last of the three PMTs of Figure 1;
  • the channel No. 2 (respectively No. 3) is the second (respectively third) channel on which a signal from the PMT corresponding to this second (respectively third) channel is detected.
  • one wonders whether or not the channel of the histogram of the triple coincidences is 2048 (for the lane 3), where 2048 is the largest channel in the example under consideration.
  • the first three lines of this timeline are relative to PMTs A, B and C which have already been mentioned above.
  • the arrows Fh denote dotted lines representing the clock signals whose period is 8 ns in the example considered.
  • the arrows Fp represent the prolongation of the dead time.
  • the line marked Tmt corresponds to the total dead time.
  • Tmm is the minimum dead time.
  • P is the measurement period of the active time. The crenals that we see on the Tmt line materialize the periods of dead time obtained by successive renewals of the minimum dead time Tmm.
  • Zone ZI corresponds to:
  • Zone Z2 corresponds to:
  • Zone Z3 corresponds to:
  • Zone Z4 corresponds to:
  • the implementation of this algorithm in the FPGA circuit 28 is based on a synchronous management of the logic pulses delivered by the DFC modules 16, 18 and 20.
  • the temporal resolution of the histograms is thus limited by the sampling frequency in the FPGA (which constitutes a digital device).
  • this implementation does not have disadvantages for the accuracy of the measurement of the coincidences between the channels.
  • the FPGA can trigger a time measurement phase only if it is in active time beforehand, that is to say outside the dead time.
  • the temporal measurement sequence and the dead time are described in the flowchart of FIG. 2 and the timing diagram of FIG. 3. They are triggered when a clock pulse is synchronous with at least one logic signal coming from the three DFC modules .
  • the minimum dead time is predefined by the user. It must always be larger than the temporal dynamics of the histograms, which defines the maximum measurable time interval.
  • the algorithm scans the arrival of a logical signal in the channels that did not initiate the phase of the temporal measurements.
  • the duration between the arrivals of the first and second pulses is expressed in number of clock ticks. This number is used to increment in real time the corresponding channel in the time histogram of double coincidences.
  • the measurement phase is triggered and the first channel of the double coincidence histogram is incremented.
  • the last channel is incremented in the purpose of keeping information from a timeout period that has been triggered.
  • the algorithm scans the arrival of a logic signal in the third way.
  • the duration between the arrivals of the first and third pulses is expressed in number of clock strokes. This number is used to increment in real time the corresponding channel in the time histogram of the triple coincidences.
  • this part is a transposition of the analog module MAC3.
  • This transposition has already been implemented in a digital coincidence measurement system by the recovery time method (see documents [4] and [5]).
  • this technique does not preserve temporal information.
  • the triggering of the dead time is performed synchronously with the start of the time measurement phase for a minimum duration which is expressed in number of clock ticks.
  • a minimum duration which is expressed in number of clock ticks.
  • Active time is the effective duration of the measurement. It is measured in real time, sampling the periods outside the idle time with the FPGA clock.
  • the timing diagram of Figure 3 gives a representation of the flow of the algorithm of Figure 2 (which is used in an example of the present invention). It should be noted that the duration of discrimination of the logic signals is also taken into account in the renewal of the dead time.
  • the management of the dead time is carried out starting from a synchronization of the signals with the clock pulses, as for the management of the temporal measurements.
  • the first channels corresponding to an optimal time of coincidence resolution are chosen, in order to avoid any loss of information that is contained in histograms of double and triple coincidences.
  • the present invention has various industrial applications.
  • the present invention overcomes a very specific drawback, due to radiation detectors, namely the random nature of the arrival of impulses with the possible presence of post-impulses of which it is difficult to characterize the temporal distribution.
  • this case corresponds for example to the 4 ⁇ - ⁇ coincidence method.
  • This method can associate two different detectors, for example a proportional counter and a scintillator type detector.
  • a time-to-digital converter is useful in the case of measuring the lifetime of a metastable state in the decay scheme of particular radionuclides.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Measurement Of Radiation (AREA)

Abstract

Procédé de mesure du taux de comptage de coïncidences, utilisant une conversion temps-numérique et une méthode de temps mort reconductible avec mesure du temps actif. On mesure le taux de comptage d' événements coïncidents entre des détecteurs de rayonnement (4, 6, 8) fonctionnant en parallèle, on convertit sous forme numérique les fluctuations temporelles des événements coïncidents et l'on utilise la méthode du temps mort reconductible avec mesure du temps actif pour éliminer tous les autres événements corrélés, susceptibles de se produire dans un même détecteur. Selon l'invention, on enregistre les distributions temporelles des intervalles de temps qui séparent les impulsions et l'on mesure le taux de comptage des événements coïncidents à l'aide des distributions temporelles enregistrées. L'invention s'applique notamment à la méthode du rapport des coïncidences triples à doubles.

Description

PROCEDE DE MESURE DU TAUX DE COMPTAGE DE COÏNCIDENCES, UTILISANT UNE CONVERSION TEMPS-NUMERIQUE ET UNE METHODE DE TEMPS MORT RECONDUCTIBLE AVEC MESURE DU TEMPS ACTIF
DESCRIPTION
DOMAINE TECHNIQUE
La présente invention concerne un procédé de mesure du taux de comptage d' événements coïncidents entre une pluralité de détecteurs de rayonnement qui fonctionnent en parallèle.
Par « rayonnement », on entend des photons ou des particules .
La mesure est effectuée grâce à l'enregistrement de la distribution d' intervalles de temps qui sont définis par la durée entre un signal de déclenchement, ou signal « start », fourni par une voie (en anglais, channel) de mesure, et un signal d'arrêt, ou signal « stop », provenant d'une autre voie de mesure. Selon le nombre de voies utilisées, on aura des signaux « stop» multiples, conduisant à une mesure de coïncidences multiples, par exemple des coïncidences doubles ou des coïncidences triples.
Indiquons dès à présent que, selon un aspect de l'invention, on combine une conversion temps-numérique (en anglais, time-to-digital conversion) , plus précisément une conversion sous forme numérique des intervalles de temps enregistrés, avec une protection contre les trains de signaux « start » et « stop » qui peuvent être à l'origine de mesures temporelles biaisées avec les techniques de mesure connues. En particulier, il s'agit de disposer d'une conversion temps-numérique adaptée aux distributions aléatoires des temps d'arrivée des signaux délivrés par les détecteurs de rayonnement.
Selon d'autres aspects de l'invention, la protection mentionnée ci-dessus utilise une méthode de temps mort reconductible (en anglais, extendible dead- time) ; ce temps mort est commun aux différentes voies ; il est construit à partir des signaux « start » et « stop » ; cette méthode de temps mort reconductible est combinée à une mesure en temps réel du temps actif (en anglais, live-time) .
L'application de la méthode du temps actif permet de mesurer les périodes effectives de comptage pour déterminer le taux des coïncidences entre les différentes voies du système de détection du rayonnement .
Au sujet de la méthode du temps mort reconductible avec mesure du temps actif, on se reportera aux documents suivants :
[1] J. Bouchard, « MTR2 : a discriminator and dead-time module used in counting Systems », Applied Radiation and Isotopes 52 (2000) 441-446 ;
[2] EP 0 574 287, « Circuit de temps mort de type reconductible », invention de J. Bouchard.
La présente invention s'applique notamment à des instruments nucléaires qui mettent en œuvre la mesure de coïncidences entre des détecteurs de rayonnement.
Dans le cas particulier de la métrologie de la radioactivité, l'invention permet par exemple de déterminer l'activité d'un radionucléide . É TAT DE LA TECHNIQUE ANTÉRIEURE
Par définition, l'expression « temps mort » désigne une période de paralysie d'un système de mesure. Cette période suit la détection d'une impulsion dans le système. Durant cette période, toute nouvelle impulsion ne peut être traitée correctement pour l'acquisition d'informations telles qu'un comptage ou une amplitude par exemple.
L'origine de la paralysie dépend du système utilisé ; il peut s'agir d'une saturation des impulsions ou de post-impulsions (en anglais, after- pulses) corrélées.
Dans ce qui suit, l'expression « méthode du temps mort reconductible » désigne une méthode permettant de gérer les paralysies du système de détection considéré. Elle consiste à empêcher, pendant une durée prédéfinie, tout traitement d'une impulsion pour en extraire une information, à la suite de la détection de cette impulsion.
Contrairement à une méthode de temps mort non reconductible (en anglais, non extendible dead- time) , tout nouveau signal arrivant durant cette période est seulement utilisé pour la prolonger de la même durée prédéfinie. Le système redevient libre ou actif lorsqu' aucune nouvelle impulsion n'a été détectée .
La méthode du temps actif consiste à mesurer le cumul des temps actifs entre les périodes de temps mort en comptant les impulsions d'une horloge. On peut alors considérer que le temps réel d'une mesure est échantillonné. Dans le domaine de l'instrumentation nucléaire, il est courant de mettre en œuvre une méthode de comptage consistant à mesurer le taux des coïncidences entre deux ou plus de deux détecteurs de rayonnement qui fonctionnent en parallèle.
Deux techniques sont généralement appliquées pour le comptage des coïncidences.
La première est fondée sur l'usage d'un circuit logique qui produit un temps de résolution de coïncidence pour chaque voie. Le comptage des coïncidences est effectué par construction d'une période de recouvrement entre les signaux logiques.
La seconde consiste à mesurer des intervalles de temps entre les voies pour obtenir l'enregistrement d'une distribution temporelle. La séquence de mesure des intervalles de temps débute par un signal « start » qui est fourni par une voie. Et les intervalles de temps sont mesurés grâce aux signaux « stop » qui apparaissent sur les autres voies.
Cette seconde méthode permet de conserver l'information donnée par les fluctuations temporelles du système de détection pour le comptage des coïncidences. A partir de l'enregistrement d'histogrammes des durées d'occurrence entre les voies, un traitement en temps différé (en anglais, off-line processing) permet d'ajuster le retard et la longueur du temps de résolution de coïncidence de façon à ce qu'ils soient adaptés au système de détection.
Dans un dispositif analogique, l'enregistrement d'un intervalle de temps nécessite deux étapes : l'intervalle de temps entre les signaux « start » et « stop » est d'abord converti sous la forme d'une amplitude puis cette dernière est exploitée par un analyseur multicanaux pour l'enregistrement de la distribution temporelle.
L'inconvénient d'une telle méthode est l'introduction d'une période de paralysie. Celle-ci est due en particulier au procédé de conversion en amplitude qui nécessite une période de charge d'un condensateur. L'analyseur multicanaux peut aussi contribuer de façon significative à la paralysie du système de mesure.
Selon les taux de comptage et le système de détection utilisé, cette paralysie du dispositif conduit à un biais dans la mesure des taux de coïncidences. A ce sujet, on se référera par exemple au document suivant :
[3] Time-to-amplitude converters and time calibrator, ORTEC®, 17 décembre 2009.
Le passage à la technologie numérique permet d'enregistrer la distribution temporelle directement à partir des mesures des durées entre les signaux « start » et « stop ». La paralysie rencontrée dans le dispositif analogique est par conséquent réduite de façon significative.
Dans une version simple, un convertisseur temps-numérique effectue la mesure temporelle à partir du comptage des impulsions d'une horloge dont la fréquence définit la résolution temporelle optimale du dispositif. Des systèmes plus sophistiqués permettent d'améliorer cette résolution grâce à l'utilisation de méthodes d'interpolation ou de la méthode du Vernier. Les systèmes connus, destinés à la mesure d'intervalles de temps, présentent un inconvénient : pour ces systèmes, le problème des paralysies ne prend pas en compte les moyens de détection qui leur sont associés.
A titre d'exemple, il existe des approches analytiques permettant de corriger les pertes de comptage dans le cas des dispositifs analogiques. Cependant, les corrections sont alors fondées sur l'utilisation d'une méthode de temps mort non reconductible. Cette méthode est mal adaptée aux paralysies engendrées par des trains d' impulsions aléatoires, conduisant à des mélanges de signaux « start » et « stop ».
Ce problème est particulièrement significatif lorsque, dans le détecteur de rayonnement, un même événement génère des post-impulsions corrélées à l'entrée d'une même voie.
Cet inconvénient ne peut être ignoré dans le domaine de la métrologie de la radioactivité. On précise que les post-impulsions peuvent d'ailleurs être dues à des états métastables qui sont caractéristiques de certains radionucléides .
EXPOSÉ DE L' INVENTION
Un but de l'invention est de remédier à cet inconvénient .
Pour ce faire, selon un aspect de l'invention, on propose de mesurer le taux de comptage des coïncidences entre différentes voies de détection qui fonctionnent en parallèle, au moyen d un système numérique fonctionnant en temps réel.
Un algorithme est mis en œuvre dans ce système. Cet algorithme combine une conversion temps- numérique, destinée à la mesure d'intervalles de temps, avec une transposition d'une méthode de temps mort reconductible. Ce temps mort est commun à l'ensemble des voies. A ce sujet, on se reportera au document suivant :
[4] J. Bouchard et al., "MAC3 : an electronic module for the processing of puises delivered by a three photomultiplier liquid scintillation counting system", Applied radiation and isotopes, 52 (2000), pp. 669-672.
De préférence, l'algorithme est mis en œuvre dans un composant programmable de type FPGA, ou réseau de portes programmables in situ (en anglais, field-programmable gâte array) , en vue d'un traitement en temps réel.
Le temps effectif de mesure est déterminé grâce à une transposition de la méthode du temps actif. A cet effet, un échantillonnage est effectué au moyen de l'horloge que comporte le composant programmable. Il convient de noter qu'une transposition sous forme numérique du module décrit dans le document [4], dans le cas d'une mesure de coïncidences à l'aide d'un temps de recouvrement, est connue par le document suivant auquel on se reportera :
[5] C. Bobin et al., "First results in the development of an on-line digital counting platform dedicated to primary measurements", Applied radiation and isotopes, 68 (2010), pp. 1519-1522.
Une particularité de l'invention réside dans le fait qu'elle ne comprend aucune voie particulière qui serait dédiée au déclenchement de la mesure des intervalles de temps : le début d'une mesure est déclenché par l'un quelconque des détecteurs. En conséquence, l'invention peut s'adapter à un système de détection symétrique, par exemple un système destiné à l'application de la méthode RCTD (en anglais, TDCR) ou méthode du Rapport des Coïncidences Triples à Doubles (en anglais, Triple to Double Coïncidence Ratio) .
Une fois la mesure déclenchée par l'une des voies, l'algorithme assure la gestion des temps d'arrivée des impulsions dans les autres voies pour établir des histogrammes temporels des coïncidences multiples entre les voies (coïncidences doubles, coïncidences triples, etc ...) .
Selon la méthode du temps mort reconductible, les post-impulsions corrélées dans les différentes voies sont utilisées pour la reconduction du temps mort.
L'information acquise, à savoir les histogrammes correspondant aux coïncidences multiples entre les voies, est enregistrée dans un ordinateur d'acquisition ; on y enregistre le temps actif mesuré.
On se reportera au document suivant :
[6] WO 2010/125062, « Procédé de mesure du taux de comptage d'impulsions, utilisant une méthode du type des temps morts reconductibles avec mesure du temps actif », invention de B.
Censier . Il existe plusieurs différences entre les divers aspects de la présente invention et le procédé décrit dans ce document [6] .
Tout d' abord, dans la présente invention, on applique en temps réel un algorithme à un système de détection comprenant au moins deux voies, alors qu'une seule voie est considérée dans le document. Et la mesure du temps actif s'applique directement par échantillonnage des périodes en dehors du temps mort à l'aide d'une horloge.
On n'enregistre donc pas, dans un ordinateur d'acquisition, les dates d'occurrence des impulsions en vue d'un post-traitement du comptage et du temps actif, alors que le procédé décrit dans le document utilise un traitement en temps différé.
De plus, un but de la présente invention est la mesure d'un taux de comptage de coïncidences entre plusieurs détecteurs.
En outre, en termes de paralysie du système de détection, la présente invention résout essentiellement un problème de traitement des post-impulsions corrélées, à savoir des post-impulsions qui sont engendrées dans les détecteurs utilisés ou qui peuvent résulter d'états métastables de certains radionucléides . C'est la raison pour laquelle le temps mort est commun aux différentes voies du système de mesure. Il ne s'agit pas seulement de périodes de paralysie engendrées par la durée de discrimination.
Du point de vue de l'algorithme, l'invention est plus proche de la transposition numérique quasiment directe du module analogique MAC3 qui a déjà été réalisée (voir les documents [4] et [5] ) .
Mais une différence majeure entre la présente invention et ce qui a déjà été effectué réside dans la mesure des coïncidences : dans l'invention, elle est réalisée grâce à l'enregistrement des distributions temporelles.
De façon précise, la présente invention a pour objet un procédé de mesure du taux de comptage d'événements coïncidents entre N détecteurs de rayonnement qui fonctionnent en parallèle et sont respectivement associés à N voies de détection, où N est un nombre entier au moins égal à 2, chaque détecteur étant apte à envoyer une impulsion électrique sur la voie de détection à laquelle il est associé lorsqu'un événement se produit dans ce détecteur, dans lequel :
-on convertit sous forme numérique les fluctuations temporelles des événements coïncidents, et
-on utilise la méthode du temps mort reconductible avec mesure du temps actif pour éliminer tous les autres événements corrélés qui sont susceptibles de se produire dans un même détecteur,
caractérisé en ce que :
-on enregistre les distributions temporelles des intervalles de temps qui séparent les impulsions, et
-on mesure le taux de comptage des événements coïncidents à l'aide des distributions temporelles enregistrées. De préférence, le temps actif est mesuré en temps réel.
De plus, le temps mort est de préférence commun aux N voies de détection.
Selon un mode de réalisation préféré de
1 ' invention :
-la mesure du taux de comptage est déclenchée par l'une des N voies de détection lors de la survenue d'un événement dans le détecteur qui lui est associé, et
-lorsque la mesure est déclenchée, on met en œuvre un algorithme qui établit des histogrammes temporels correspondant aux coïncidences multiples entre les N voies de détection, c'est-à-dire aux événements coïncidents entre P détecteurs parmi les N détecteurs de rayonnement, où P parcourt l'ensemble des nombres entiers allant de 2 à N, à partir des temps d'arrivée des impulsions dans les N-l autres voies de détection .
Dans la présente invention, les détecteurs peuvent être identiques les uns aux autres.
Dans ce cas, selon un mode de réalisation particulier de l'invention, N est égal à 3, les détecteurs sont des photomultiplicateurs et l'on utilise le procédé pour mettre en œuvre la méthode du rapport des coïncidences triples à doubles.
Mais les détecteurs peuvent aussi ne pas être identiques les uns aux autres.
Dans ce cas, selon un autre mode de réalisation particulier de l'invention, N est égal à 2, les détecteurs sont respectivement un détecteur de photons gamma et un détecteur d'électrons et l'on utilise le procédé pour mettre en œuvre la méthode des coïncidences bêta-gamma (en anglais, beta-gamma coïncidence method) .
BRÈVE DESCRIPTION DES DESSINS
La présente invention sera mieux comprise à la lecture de la description d'exemples de réalisation donnés ci-après, à titre purement indicatif et nullement limitatif, en faisant référence aux dessins annexés sur lesquels :
-la figure 1 est une vue schématique d'une chaîne de mesure qui est destinée à déterminer l'activité d'un radionucléide par la méthode RCTD et dans laquelle est mis en œuvre un exemple du procédé, objet de l'invention,
-la figure 2 est un organigramme d'un algorithme qui est utilisé dans cet exemple et qui traite en parallèle des mesures temporelles et la gestion du temps mort reconductible, et
-la figure 3 est un chronogramme relatif à cet algorithme.
EXPOSE DE TAILLE DE MODES DE REALISATION PARTICULIERS
L'exemple que l'on va décrire est relatif à une application de l'invention à la méthode RCTD qui est couramment utilisée pour la mesure d'une activité par comptage à scintillation liquide (en anglais, liquid scintillation counting) . La mise en œuvre de cette méthode permet de mettre en évidence les avantages de l'invention dans le cas où l'arrivée des signaux « start » et « stop » est aléatoire en raison de l'existence de post-impulsions corrélées qui peuvent conduire à des périodes confuses dans la mesure des intervalles de temps.
Les deux phénomènes les plus importants à l'origine des post-impulsions sont :
-l'émission lumineuse en scintillation liquide qui résulte de différents mécanismes physico-chimiques conduisant à un étalement de la distribution temporelle des photons émis, et
-l'ionisation du gaz résiduel qui est présent dans les photomultiplicateurs utilisés pour la mise en œuvre de la méthode RCTD.
La chaîne de mesure, sur laquelle on reviendra par la suite, comprend trois voies de comptage, ou voies de détection, qui sont identiques ; et chacune de ces voies commence par un photomultiplicateur.
La détection est faite de façon symétrique de sorte que chacune des voies est capable de déclencher une séquence de mesure et de temps mort reconductible. Il n'existe donc aucune voie spécifique, dédiée au déclenchement d'un processus de mesure d'intervalles de temps.
Une fois la période de mesure initiée par une première impulsion électrique qui a été envoyée à l'une des trois voies, les deux autres voies sont utilisées pour mettre en œuvre des histogrammes temporels correspondant à des deuxième et troisième impulsions (pour le comptage des coïncidences doubles et triples) .
Dans le cas d'une mesure par la méthode RCTD, des données sont enregistrées dans un ordinateur d'acquisition sous la forme de deux histogrammes temporels, correspondant aux temps d'arrivée des deuxièmes et troisièmes impulsions. L'enregistrement du temps actif total (cumul des périodes en dehors du temps mort) permet de calculer les taux des coïncidences doubles et triples entre les voies.
Rappelons que, conformément à l'invention, on mesure des événements coïncidents mais seulement entre différents détecteurs. Dans un même détecteur, il peut y avoir de tels événements mais ils n'apportent aucune information sur la mesure. On se protège de ces événements corrélés dans un même détecteur grâce à la méthode du temps mort reconductible avec mesure du temps actif : on ne tient compte que du premier événement dans un détecteur mais pas d'autres événements ultérieurs dans ce détecteur pendant la période du temps mort.
A titre purement indicatif et nullement limitatif, on a réalisé un dispositif de mesure par la méthode RCTD à l'aide d'une carte numérique commercialement disponible, à savoir un kit de développement Altéra®, équipé d'un circuit FPGA Stratix® III.
Et, dans ce circuit logique, on a implémenté un programme permettant d'enregistrer des histogrammes temporels avec une profondeur d'échantillonnage qui a été fixée à 2048 canaux. De plus, on a choisi une fréquence d'échantillonnage égale à 125 MHz, ce qui conduit à une résolution temporelle minimale de 8 ns .
Pour 1 ' implémentation des histogrammes, le temps par canal est défini selon un multiple de la résolution temporelle minimale.
Cette augmentation de la dynamique temporelle se fait toutefois au prix d'une perte de résolution ; mais la souplesse de programmation du circuit FPGA permet de remédier à cette perte en augmentant la profondeur d' échantillonnage .
Le temps total de mesure et le temps mort minimum sont définis par l'utilisateur selon un nombre de coups horloge .
Et la liaison entre la partie numérique de la chaîne de mesure et l'ordinateur d'acquisition est réalisée grâce à une connexion Ethernet.
La figure 1 est une vue schématique de la chaîne de mesure dans laquelle est mis en œuvre un exemple du procédé, objet de l'invention.
Dans cette chaîne de mesure, la détection de la lumière a lieu dans une chambre optique (non représentée) . Dans cette chambre optique, on introduit un flacon 2 contenant un mélange d'un scintillateur liquide et d'une source radioactive dont on veut déterminer l'activité par la méthode RCTD.
Pour détecter la lumière de scintillation, on utilise trois photomultiplicateurs 4, 6, 8 qui sont disposés de façon symétrique autour du flacon 2, à 120° les uns des autres. A la suite d'une désintégration radioactive, un rayonnement ionisant est émis. Il en résulte un dépôt d'énergie dans le scintillateur liquide. Cela conduit à l'émission de photons lumineux qui se distribuent aléatoirement entre les trois photomultiplicateurs .
Dans la chaîne de mesure, chaque photomultiplicateur assure la conversion d'un photon lumineux incident en un photoélectron. Cette conversion est effectuée par l'intermédiaire d'une photocathode située à l'entrée du photomultiplicateur et dépend du rendement quantique de la photocathode.
De plus, le photomultiplicateur comprend une chaîne de dynodes. Dans cette dernière a lieu un processus de multiplication des photoélectrons qui sont produits au niveau de la photocathode. On obtient ainsi un courant suffisamment important pour pouvoir être transformé en une impulsion de tension qui soit exploitable par un amplificateur rapide.
Cette impulsion est envoyée sur une voie de détection qui relie le photomultiplicateur à l'amplificateur rapide. Sur la figure 1, on a d'ailleurs représenté les voies de détection 10, 12, 14 qui relient respectivement les photomultiplicateurs 4, 6, 8 aux amplificateurs rapides 16, 18, 20.
Ces amplificateurs 16, 18, 20 sont respectivement associés à des dispositifs analogiques 22, 24, 26, à savoir des DFC (en anglais, CFD) c'est-à- dire des discriminateurs à fraction constante (en anglais, constant fraction discriminators) . Le signal délivré par chaque amplificateur alimente le DFC qui lui est associé.
La chaîne de mesure comprend aussi un dispositif numérique 28 qui est constitué par un FPGA dans l'exemple décrit. Dans ce FPGA, on a programmé la conversion temps-numérique et la protection dont il a été question plus haut et sur lesquelles on reviendra par la suite.
L'ensemble des DFC permet de passer de la partie analogique du système de détection au dispositif numérique de conversion temps-numérique 28. Chaque DFC produit un signal logique qui est directement exploitable par ce dispositif 28 et présente l'avantage de réduire les fluctuations temporelles que l'on rencontre généralement dans le cas d'une discrimination par seuil classique.
Plus précisément, les trois DFC, qui équipent respectivement les voies de détection, fournissent des impulsions logiques. Ces impulsions matérialisent les signaux « start » et « stop » qui arrivent à l'entrée du dispositif numérique 28.
La chaîne de mesure comprend également un ordinateur d'acquisition 30 qui traite les histogrammes temporels fournis par le dispositif numérique 28 et détermine les taux de comptage recherchés. Dans l'exemple décrit, cet ordinateur est relié au dispositif numérique 28 par une liaison Ethernet. De plus, l'ordinateur 30 est muni d'un dispositif 32 d'affichage des résultats des mesures.
On va maintenant considérer l'algorithme qui est représenté sur la figure 2. Cet algorithme est mis en œuvre dans le FPGA 28 faisant partie de la chaîne de mesure représentée sur la figure 1 et traite en parallèle la conversion temps-numérique et la gestion du temps mort reconductible.
On va également considérer le chronogramme qui est représenté sur la figure 3 et qui est relatif à cet algorithme.
L'algorithme est tel que représenté sur la figure 2. Ci-après, on en précise simplement quelques éléments pour clarifier certaines abréviations.
En 34, on part d'un état où la chaîne de mesure n'est pas paralysée ; le temps actif est incrémenté à chaque coup d'horloge. En 36, on se demande si oui ou non on détecte un signal de l'un quelconque des photomultiplicateurs (en abrégé : PMT) . En 38, on se demande si oui ou non on a trois signaux synchrones provenant respectivement des trois PMT. En 40, le compteur de temps mort est mis à une valeur prédéfinie .
En 42, on se demande si oui ou non deux PMT parmi les trois PMT fournissent des signaux synchrones. En 44, on se demande si oui ou non l'un des PMT fournit un signal. En 46, le numéro du canal de l'histogramme des coïncidences doubles vaut 1 et le numéro du canal de l'histogramme des coïncidences triples vaut 1. En 48, le numéro du canal de l'histogramme des coïncidences triples vaut 1.
En 50, lors d'un nouveau coup d'horloge, le compteur de temps mort est décrémenté de 1. En 52, le numéro du canal de l'histogramme des coïncidences doubles est augmenté de 1 et le numéro du canal de l'histogramme des coïncidences triples est augmenté de 1. En 54, le numéro du canal de l'histogramme des coïncidences triples est augmenté de 1 et il y a un nouveau coup d'horloge. En 56, on se demande si oui ou non le compteur de temps mort est mis à 0.
En 58, on se demande si oui ou non on a un signal de PMT sur la voie n°3. A ce propos, on apporte les précisions suivantes :
-la voie n°l est la première voie sur laquelle on détecte un signal provenant du PMT correspondant, à savoir l'un des PMT A, B et C avec les notations de la figure 3 ;
-le PMT A est l'un des trois PMT de la figure 1, le PMT B est l'un des deux autres PMT de la figure 1 et le PMT C est le dernier des trois PMT de la figure 1 ; et
-la voie n°2 (respectivement n°3) est la deuxième (respectivement troisième) voie sur laquelle on détecte un signal provenant du PMT correspondant à cette deuxième (respectivement troisième) voie.
En 60, on se demande si oui ou non on détecte un signal de PMT sur chacune des voies n°2 et 3. En 62, on se demande si oui ou non le canal de l'histogramme des coïncidences triples vaut 2048 (pour la voie n°3), où 2048 est le canal le plus grand dans l'exemple considéré.
En 64, on se demande si oui ou non le canal de l'histogramme des coïncidences doubles vaut 2048 et l'on se demande si oui ou non le canal de l'histogramme des coïncidences triples vaut 2048 (pour les voies n°2 et 3) . En 66, on se demande si oui ou non les signaux qui proviennent des PMT correspondant respectivement aux voies n°2 et 3 sont synchrones.
Considérons maintenant le chronogramme de la figure 3.
Les trois premières lignes de ce chronogramme sont relatives aux PMT A, B et C qui ont déjà été mentionnés plus haut. Les flèches Fh désignent des pointillés représentant les signaux d'horloge dont la période vaut 8 ns dans l'exemple considéré. Les flèches Fp représentent la prolongation du temps mort. La ligne notée Tmt correspond au temps mort total. Tmm désigne le temps mort minimal. P désigne la période de mesure du temps actif. Les crénaux que l'on voit sur la ligne Tmt matérialisent les périodes de temps mort obtenues par des reconductions successives du temps mort minimal Tmm.
La zone ZI correspond :
-au déclenchement de la mesure par le PMT A qui fournit alors une première impulsion,
-à une deuxième impulsion, fournie par le PMT B, avec un retard de 1 coup d'horloge, et
-à l'incrémentation du 2ieme canal de l'histogramme des coïncidences doubles.
La zone Z2 correspond :
-à une troisième impulsion, fournie par le PMT C, avec un retard de 7 coups d'horloge, et
-à l'incrémentation du 8ieme canal de l'histogramme des coïncidences triples.
La zone Z3 correspond :
-au déclenchement de la mesure par le PMT A et le PMT B, et -à l'incrémentation du 1er canal de l'histogramme des coïncidences doubles.
La zone Z4 correspond :
-à une troisième impulsion, fournie par le PMT C, avec un retard de 1 coup d'horloge, et
-à l'incrémentation du 2ieme canal de l'histogramme des coïncidences triples.
On revient maintenant sur l'algorithme de la figure 2.
L' implémentation de cet algorithme dans le circuit FPGA 28 (figure 1) est fondée sur une gestion synchrone des impulsions logiques délivrées par les modules DFC 16, 18 et 20.
La résolution temporelle des histogrammes est ainsi limitée par la fréquence d'échantillonnage dans le FPGA (qui constitue un dispositif numérique) . Toutefois, cette implémentation ne présente pas d'inconvénients pour la justesse de la mesure des coïncidences entre les voies.
Durant le temps d'acquisition défini par l'utilisateur, le FPGA ne peut déclencher une phase de mesures temporelles que s'il est au préalable en temps actif c'est-à-dire en dehors du temps mort.
La séquence de mesure temporelle et le temps mort sont décrits dans l'organigramme de la figure 2 et le chronogramme de la figure 3. Ils sont déclenchés lorsqu'une impulsion d'horloge est synchrone avec au moins un signal logique provenant des trois modules DFC.
II convient de noter que le traitement des histogrammes des intervalles de temps et du temps mort de type reconductible est implémenté selon deux processus qui s'exécutent en parallèle à partir des mêmes signaux logiques.
Le temps mort minimum est prédéfini par l'utilisateur. Il doit être toujours plus grand que la dynamique temporelle des histogrammes, dynamique qui définit l'intervalle de temps maximum mesurable.
Une fois la séquence déclenchée, l'algorithme scrute l'arrivée d'un signal logique dans les voies qui n'ont pas initié la phase des mesures temporelles. La durée entre les arrivées des première et deuxième impulsions est exprimée en nombre de coups d'horloge. Ce nombre est utilisé pour incrémenter en temps réel le canal correspondant dans l'histogramme temporel des coïncidences doubles.
Lorsque deux signaux logiques provenant de deux voies différentes sont détectés de façon synchrone par une impulsion d'horloge, la phase de mesure est déclenchée et le premier canal de l'histogramme des coïncidences doubles est incrémenté.
Dans le cas particulier où aucune seconde impulsion n'est détectée durant la période correspondant à la profondeur de codage de l'histogramme temporel (valeur minimale dans l'exemple décrit : 2047 X 8 ns = 16,376 με), le dernier canal est incrémenté dans le but de conserver l'information d'une période de temps mort qui a été déclenchée.
Lorsque deux signaux logiques provenant de deux voies différentes ont déjà été détectés durant la séquence en cours (simultanément ou non), l'algorithme scrute l'arrivée d'un signal logique dans la troisième voie. La durée entre les arrivées des première et troisième impulsions est exprimée en nombre de coups d'horloge. Ce nombre est utilisé pour incrémenter en temps réel le canal correspondant dans l'histogramme temporel des coïncidences triples.
Lorsque les deux signaux logiques provenant des deux autres voies (2eme et 3eme voies) sont détectés de façon synchrone par une impulsion d'horloge, le même canal, correspondant à la durée entre les arrivées des première et deuxième impulsions (exprimée en nombre de coups d'horloge), est incrémenté dans les histogrammes des coïncidences doubles et triples.
Dans le cas particulier où aucune troisième impulsion n'est détectée durant la période correspondant à profondeur de codage, le dernier canal de l'histogramme des coïncidences triples est incrémenté .
En ce qui concerne la partie de l'algorithme qui est dédiée à la gestion du temps mort de type reconductible et à la mesure du temps actif, cette partie est une transposition du module analogique MAC3. Cette transposition a déjà été mise en œuvre dans un système numérique de mesure de coïncidences par la méthode du temps de recouvrement (voir les documents [4] et [5]) . On rappelle que cette technique ne conserve pas l'information temporelle.
Le déclenchement du temps mort est effectué de façon synchrone avec le début de la phase de mesure temporelle pour une durée minimale qui est exprimée en nombre de coups d'horloge. Selon le principe du temps mort de type reconductible, tout nouveau signal logique provenant des trois modules DFC et arrivant durant une période de temps mort prolonge cette période de la même durée minimale (qui est définie par l'utilisateur) .
Le temps actif représente la durée effective de la mesure. Il est mesuré en temps réel, en échantillonnant les périodes en dehors du temps mort avec l'horloge du FPGA.
Le chronogramme de la figure 3 donne une représentation du déroulement de l'algorithme de la figure 2 (qui est utilisé dans un exemple de la présente invention) . Il convient de noter que la durée de discrimination des signaux logiques est aussi prise en compte dans la reconduction du temps mort.
Il convient également de noter que la gestion du temps mort est réalisée à partir d'une synchronisation des signaux avec les impulsions d'horloge, comme pour la gestion des mesures temporelles.
Après correction des histogrammes en divisant le contenu des canaux par la mesure du temps actif, on peut procéder au calcul des taux de comptage des coïncidences multiples. Les valeurs de ces taux sont données par le cumul du contenu des canaux correspondant à la région temporelle qui est définie par l'utilisateur.
Dans le cas de la méthode RCTD, on choisit les premiers canaux correspondant à un temps optimal de résolution de coïncidence, dans le but d'éviter toute perte des informations qui sont contenues dans les histogrammes des coïncidences doubles et triples.
La présente invention a diverses applications industrielles.
En effet, dans l'industrie, de nombreux domaines nécessitent des mesures faisant intervenir un convertisseur temps-numérique, par exemple les mesures dans le domaine des temps-fréquences ou les mesures dans le domaine de la photonique.
En particulier, dans le domaine de la métrologie de la radioactivité et, plus généralement, de la physique nucléaire et des particules, la présente invention remédie à un inconvénient bien spécifique, dû aux détecteurs de rayonnement, à savoir le caractère aléatoire des temps d'arrivée des impulsions avec la présence éventuelle de post-impulsions dont il est difficile de caractériser la distribution temporelle.
Par ailleurs, on a donné des exemples de l'invention dans lesquels on utilise trois détecteurs de rayonnement. Mais la présente invention n'est pas limitée à ce cas : l'homme du métier peut adapter les exemples donnés à un cas où l'on utilise plus de trois détecteurs de rayonnement.
L'homme du métier peut même adapter ces exemples au cas où l'on utilise seulement deux détecteurs de rayonnement. Dans le domaine de la métrologie de la radioactivité, ce cas correspond par exemple à la méthode des coïncidences 4πβ-γ. Cette méthode peut associer deux détecteurs différents, par exemple un compteur proportionnel et un détecteur de type scintillateur . En outre, un convertisseur temps-numérique est utile dans le cas de la mesure de la durée vie d'un état métastable dans le schéma de désintégration de radionucléides particuliers.

Claims

REVENDICATIONS
1. Procédé de mesure du taux de comptage d' événements coïncidents entre N détecteurs de rayonnement (4, 6, 8) qui fonctionnent en parallèle et sont respectivement associés à N voies de détection (10, 12, 14), où N est un nombre entier au moins égal à 2, chaque détecteur étant apte à envoyer une impulsion électrique sur la voie de détection à laquelle il est associé lorsqu'un événement se produit dans ce détecteur, dans lequel :
-on convertit sous forme numérique les fluctuations temporelles des événements coïncidents, et
-on utilise la méthode du temps mort reconductible avec mesure du temps actif pour éliminer tous les autres événements corrélés qui sont susceptibles de se produire dans un même détecteur,
caractérisé en ce que :
-on enregistre les distributions temporelles des intervalles de temps qui séparent les impulsions, et
-on mesure le taux de comptage des événements coïncidents à l'aide des distributions temporelles enregistrées.
2. Procédé selon la revendication 1, dans lequel le temps actif est mesuré en temps réel.
3. Procédé selon l'une quelconque des revendications 1 et 2, dans lequel le temps mort est commun aux N voies de détection (10, 12, 14) .
4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel :
-la mesure du taux de comptage est déclenchée par l'une des N voies de détection (10, 12, 14) lors de la survenue d'un événement dans le détecteur qui lui est associé, et
-lorsque la mesure est déclenchée, on met en œuvre un algorithme qui établit des histogrammes temporels correspondant aux coïncidences multiples entre les N voies de détection, c'est-à-dire aux événements coïncidents entre P détecteurs parmi les N détecteurs de rayonnement, où P parcourt l'ensemble des nombres entiers allant de 2 à N, à partir des temps d'arrivée des impulsions dans les N-l autres voies de détection .
5. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel les détecteurs (4, 6, 8) sont identiques les uns aux autres.
6. Procédé selon la revendication 5, dans lequel N est égal à 3, les détecteurs sont des photomultiplicateurs (4, 6, 8) et l'on utilise le procédé pour mettre en œuvre la méthode du rapport des coïncidences triples à doubles.
7. Procédé selon l'une quelconque des revendications 1 à 4, dans lequel les détecteurs ne sont pas identiques les uns aux autres.
8. Procédé selon la revendication 7, dans lequel N est égal à 2, les détecteurs sont
respectivement un détecteur de photons gamma et un détecteur d'électrons et l'on utilise le procédé pour mettre en œuvre la méthode des coïncidences bêta-gamma.
EP11804670.5A 2010-12-20 2011-12-19 Procede de mesure du taux de comptage de coïncidences, utilisant une conversion temps-numerique et une methode de temps mort reconductible avec mesure du temps actif Withdrawn EP2656111A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1060771A FR2969310B1 (fr) 2010-12-20 2010-12-20 Procede de mesure du taux de comptage de coincidences, utilisant une conversion temps-numerique et une methode de temps mort reconductible avec mesure du temps actif
PCT/EP2011/073213 WO2012084802A1 (fr) 2010-12-20 2011-12-19 Procede de mesure du taux de comptage de coïncidences, utilisant une conversion temps-numerique et une methode de temps mort reconductible avec mesure du temps actif

Publications (1)

Publication Number Publication Date
EP2656111A1 true EP2656111A1 (fr) 2013-10-30

Family

ID=44275654

Family Applications (1)

Application Number Title Priority Date Filing Date
EP11804670.5A Withdrawn EP2656111A1 (fr) 2010-12-20 2011-12-19 Procede de mesure du taux de comptage de coïncidences, utilisant une conversion temps-numerique et une methode de temps mort reconductible avec mesure du temps actif

Country Status (4)

Country Link
US (1) US20130268224A1 (fr)
EP (1) EP2656111A1 (fr)
FR (1) FR2969310B1 (fr)
WO (1) WO2012084802A1 (fr)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6326347B2 (ja) * 2014-10-20 2018-05-16 日本電子株式会社 ライブタイム比演算回路、ライブタイム比演算方法、放射線検出装置、および試料分析装置
CN112068182B (zh) * 2020-09-15 2022-05-17 中国核动力研究设计院 基于多丝正比室的4πβ-γ符合测量装置及测量方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3786256A (en) * 1971-11-18 1974-01-15 Nat Nuclear Corp Method and apparatus for nuclear fuel assay with a neutron source and coincident fission neutron detectors
FR2692088B1 (fr) * 1992-06-03 1994-07-22 Commissariat Energie Atomique Circuit de temps mort de type reconductible.
US5953421A (en) * 1995-08-16 1999-09-14 British Telecommunications Public Limited Company Quantum cryptography
US7818047B2 (en) * 2001-11-09 2010-10-19 Nova R&D, Inc. X-ray and gamma ray detector readout system
US8194813B2 (en) * 2004-10-19 2012-06-05 Lawrence Livermore National Security, Llc Absolute nuclear material assay using count distribution (LAMBDA) space
FR2945129B1 (fr) 2009-04-29 2011-06-03 Commissariat Energie Atomique Procede de mesure du taux de comptage d'implusions, utilisant une methode du type des temps morts reconductibles avec mesure du temps actif

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO2012084802A1 *

Also Published As

Publication number Publication date
WO2012084802A1 (fr) 2012-06-28
FR2969310B1 (fr) 2014-11-28
US20130268224A1 (en) 2013-10-10
FR2969310A1 (fr) 2012-06-22

Similar Documents

Publication Publication Date Title
JP6657398B2 (ja) パルスパイルアップイベントのリアルタイム処理の方法及びシステム
Powolny et al. Time-based readout of a silicon photomultiplier (SiPM) for time of flight positron emission tomography (TOF-PET)
EP2257835B1 (fr) Détecteur de rayonnement de photon unique
US9354332B2 (en) Detector array with time-to-digital conversion having improved temporal accuracy
EP1550884B1 (fr) Système de détection de rayonnements comportant une pluralité d'éléments détecteurs ayant chacun un circuit d'acquisition et de comptage
Tyndall et al. A 100Mphoton/s time-resolved mini-silicon photomultiplier with on-chip fluorescence lifetime estimation in 0.13 μm CMOS imaging technology
US8969814B2 (en) System and method of determining timing triggers for detecting gamma events for nuclear imaging
EP2771716B1 (fr) Procede de determination de la dose d'irradiation deposee dans un scintillateur par un rayonnement ionisant et dispositif associe
EP2507652A1 (fr) Procédé d'étalonnage d'un détecteur de rayonnement x
EP0131478B1 (fr) Procédé de traitement des impulsions de localisation délivrées par une gamma caméra et gamma caméra mettant en oeuvre ce procédé
EP0937265A1 (fr) Dispositif et procede de determination de la position presumee d'un evenement par rapport a un ensemble de photodetecteurs, et application aux gamma-c
CA2765336C (fr) Procede apte a discriminer une composante gamma et une composante neutronique dans un signal electronique
WO2012084802A1 (fr) Procede de mesure du taux de comptage de coïncidences, utilisant une conversion temps-numerique et une methode de temps mort reconductible avec mesure du temps actif
FR2615961A1 (fr) Procede de prise en compte des impulsions de localisation delivrees par une gamma camera
JP4330904B2 (ja) 放射線検出方法および装置
EP3835831B1 (fr) Procédé de détection de radiations et système associé mettant en oeuvre une discrimination des coïncidences neutron-gamma
WO2010125062A1 (fr) Procédé de mesure du taux de comptage d'impulsions, utilisant une méthode du type des temps morts reconductibles avec mesure du temps actif
EP3051316B1 (fr) Circuit électronique configuré pour contrôler un module d'acquisition d'un dispositif de détection de photons, module d'acquisition, procédé de contrôle, et programme d'ordinateur correspondant
EP2668521A1 (fr) Procede et dispositif de correction d'empilement d'impulsions pour un detecteur de rayonnement
Xu et al. Comparison of digital and analog silicon photomultiplier for positron emission tomography application
Rech et al. Toward Constraintless Time-Correlated Single-Photon Counting Measurements: A New Method to Remove Pile-Up Distortion
Gasparini et al. Characterizing single-and multiple-timestamp time of arrival estimators with digital SiPM PET detectors
US20230204799A1 (en) Radiation detector
FR1465698A (fr) Discriminateur de particules nucléaires
US20190170882A1 (en) System and method to unpile overlapping pulses

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20130620

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AL AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HR HU IE IS IT LI LT LU LV MC MK MT NL NO PL PT RO RS SE SI SK SM TR

DAX Request for extension of the european patent (deleted)
17Q First examination report despatched

Effective date: 20140610

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE APPLICATION IS DEEMED TO BE WITHDRAWN

18D Application deemed to be withdrawn

Effective date: 20160701