Brevet FR2945129A1 — publié le 20101105

PROCEDE DE MESURE DU TAUX DE COMPTAGE D'IMPULSIONS, UTILISANT UNE METHODE DU TYPE DES TEMPS MORTS RECONDUCTIBLES AVEC MESURE DU TEMPS ACTIF DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne un procédé de mesure du taux de comptage d'impulsions (en anglais, pulse count rate), utilisant une méthode du type des temps morts reconductibles avec mesure du temps actif (en anglais, extendible dead-times with live-time measurement). Elle s'applique notamment aux mesures d'activités de radionucléides, plus particulièrement en métrologie de la radioactivité. ÉTAT DE LA TECHNIQUE ANTÉRIEURE La mesure de l'activité d'un radionucléide consiste à déterminer le nombre moyen de désintégrations subies par ce radionucléide par unité de temps. Dans la plupart des systèmes de mesure, la détection d'une radiation émise par suite d'une désintégration se traduit par la formation d'une impulsion de courant à la sortie des détecteurs que comportent de tels systèmes ; et cette impulsion est amplifiée et mise en forme électroniquement si nécessaire. La détection d'un évènement revient donc à la détection d'une impulsion dépassant un seuil fixé par les utilisateurs. 2 Comme tout système de détection a un temps de réponse non nul, il présente des périodes d' � aveuglement � qu'on appellera des temps morts : la détection d'un évènement rend le système inapte à détecter les dépassements de seuil pendant un certain laps de temps. De plus, les caractéristiques de l'instrumentation utilisée et de l'électronique associée peuvent introduire des artefacts de mesure et des temps de traitement du signal, qui sont aussi des sources de temps mort. Ces phénomènes sont d'autant plus pénalisants que le taux de comptage mesuré est grand, ce qui peut causer des pertes de comptage et une mauvaise estimation de l'activité réelle de la source radioactive étudiée. Il est alors nécessaire d'estimer le temps actif de la mesure, c'est-à-dire le temps pendant lequel le système de mesure est considéré comme apte à détecter le rayonnement émis par la source. Indiquons dès à présent que l'invention utilise une transposition, sous la forme d'un algorithme, de la méthode dite "des temps morts reconductibles avec mesure du temps actif". Cette dernière consiste (1) à appliquer un temps mort reconductible, de durée donnée, après chaque évènement détecté, et (2) à mesurer la durée totale des périodes en dehors des temps morts, cette durée constituant le temps actif de la mesure. Cette méthode permet la gestion des temps mort de manière robuste, jusqu'à de grands taux de comptage et pour une grande diversité de systèmes de mesures. Elle a déjà été mise en oeuvre à l'aide de 3 modules électroniques dédiés, que l'on peut insérer dans une chaîne d'acquisition analogique. Au sujet de cette méthode, 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. Avec le développement des numériseurs rapides et du traitement des informations en temps réel, il est utile de transposer cette méthode dans le domaine numérique. En effet, l'un des problèmes des chaînes d'acquisition numérique, utilisant des 15 numériseurs rapides, est la gestion des importants flux de données qui sont générés, par exemple lors du traitement et du stockage des données. L'utilisation de composants programmables tels que les FPGA, c'est-à-dire les réseaux de portes programmables in 20 situ (en anglais, field-programmable gate arrays), permet de réduire ces flux de données pour n'en conserver que les informations nécessaires. L'algorithme utilisé dans l'invention permet d'appliquer la méthode des temps morts 25 reconductibles avec mesure du temps actif à partir des données de base que sont la durée de chaque impulsion détectée et la durée de chaque intervalle de temps séparant des impulsions. Il est donc applicable à tout système de mesure fournissant ces informations à un 30 ordinateur d'acquisition. 10 4 L'idée de base pour le traitement des temps morts est l'application volontaire d'une période de temps mort, de durée donnée, après chaque impulsion. Cela permet, dans une certaine mesure, de s'affranchir des problèmes liés aux durées variables des impulsions, à l'existence de temps de traitement dans la chaîne d'acquisition, ou à l'existence d'artefacts résultant des impulsions. Cette méthode peut être mise en oeuvre en suivant principalement deux stratégies différentes qui utilisent respectivement un temps mort fixe et un temps mort reconductible. Dans le cas du temps mort fixe, on applique une durée de temps mort donnée après chaque impulsion. Pendant la durée de ce temps mort, le système de déclenchement sur seuil est inhibé. Dans le cas du temps mort reconductible, le temps mort est reconduit à chaque déclenchement sur seuil, y compris lors d'une période de temps mort. La technique du temps mort fixe est largement utilisée car elle est relativement simple à mettre en oeuvre. Elle présente cependant des inconvénients qui sont d'autant plus importants que les taux de comptage sont élevés. En effet, si une impulsion N+1 commence alors que le système de déclenchement se trouve dans un temps mort imposé par l'impulsion précédente N, elle ne sera pas détectée par le système inhibé. Si, du fait de l'impulsion N+1, le retour du signal en dessous du seuil s'effectue en partie en dehors de la période de temps mort imposé, le temps mort est sous-estimé. A ce titre, le terme de temps mort fixe est critiquable : il concerne la technique utilisée mais ne représente pas fidèlement la réalité expérimentale dans laquelle le temps mort véritable n'est jamais vraiment 5 fixe et dépend notamment du taux de comptage, de la mise en forme des impulsions et des caractéristiques du détecteur utilisé. C'est pourquoi, l'utilisation de temps morts reconductibles est nécessaire. Se pose ensuite la question de l'évaluation du temps actif de la mesure. Dans le cas d'un temps mort fixe, l'utilisation de méthodes statistiques a posteriori est largement répandue. Ces méthodes sont fondées sur une modélisation plus ou moins fine du système de mesure et de sa chaîne d'acquisition et dépendent en ce sens des caractéristiques du système utilisé. Elles deviennent rapidement complexes à mesure que l'on tente de rendre compte de caractéristiques de plus en plus précises. Cependant, malgré cette complexité, elles ne peuvent jamais rendre compte de la totalité des possibilités expérimentales (telles que les caractéristiques de l'électronique d'amplification et de mise en forme, les caractéristiques de l'électronique de déclenchement et les artefacts de mesure) et s'avèrent souvent sensibles aux détails de schémas de désintégration mal connus. Elles présentent donc un défaut de robustesse. Ces méthodes sont en outre difficilement applicables aux temps morts reconductibles, la modélisation devenant trop complexe. 6 En ce qui concerne plus particulièrement l'instrumentation numérique, on peut citer l'exemple de la technique DCC (pour Digital Coincidence Counting). A ce sujet, on se reportera aux documents suivants : [3] J.D. Keightley, G.C. Watt, "Digital coincidence counting (DCC) and its use in the corrections for out-of-channel gamma events in 4i4-y coincidence counting", Applied Radiation and Isotopes 56 (2002) 205-210 ; [4] John Keightley and Tae Soon Park, "Digital coincidence counting for radionuclide standardization", Metrologia 44 (2007) S32-S35. Dans les systèmes qui y sont décrits, une carte d'acquisition numérise le signal à traiter ; puis un FPGA est chargé du déclenchement sur seuil et de la datation des évènements détectés. L'ensemble des échantillons numérisés constituant une impulsion détectée est alors envoyé à un ordinateur d'acquisition. Et le traitement des temps morts est de type � temps mort fixe � avec estimation statistique du temps actif. Ces systèmes illustrent les inconvénients associés à la structure du traitement de données et à la méthode utilisée pour le traitement des temps morts. D'une part, la numérisation du signal et le stockage de l'intégralité des échantillons numérisés formant les impulsions engendrent d'importants flux de données entre la carte et l'ordinateur d'acquisition. Pour un taux de comptage élevé, typiquement supérieur à quelques dizaines de coups par secondes, le système d'acquisition introduit une nouvelle source de temps 7 mort difficilement caractérisable, qui est liée au débit irrégulier et limité des données entre la carte et l'ordinateur. D'autre part, l'utilisation du temps mort fixe limite également les taux de comptage accessibles sans perte et complique l'analyse du signal car, pour ce faire, elle nécessite le développement de méthodes complexes et spécifiques. EXPOSÉ DE L'INVENTION La présente invention a pour but de remédier aux inconvénients précédents. Elle a pour objet un 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, caractérisé en ce que : - on choisit une valeur Tm de temps mort reconductible, - on acquiert N impulsions successives I1r I2r IN d'un signal électrique impulsionnel, on détermine la durée Di de chaque impulsion acquise Ti, ainsi que la durée Ti de l'intervalle de temps entre l'impulsion précédente Ii-1 et l'impulsion Ti, Ti étant choisi égal à l'intervalle de temps entre le début de l'acquisition et le début de la première impulsion, on effectue une itération en donnant successivement à l'indice i les valeurs 1 à N, et : • pour une valeur i de l'indice, on met à jour la durée totale Ttot de l'acquisition, en ajoutant 8 Ti et Di à la valeur que Ttot avait auparavant, T t étant initialement égal à 0, • on compare Ti à Tm, • et si Ti est inférieur ou égal à Tm, on passe à l'indice suivant i+1, • ou si Ti est supérieur à Tm, on met à jour la durée du temps actif Ta_tot en ajoutant la différence Ti-Tm à la valeur que Ta_tot avait auparavant, Ta_tot étant initialement égal à 0, puis on met à jour une variable C égale au nombre d'impulsions détectées en dehors du temps mort, en ajoutant 1 à la valeur que C avait auparavant, C étant initialement égal à 0, et l'on passe à l'indice i+1, et - on détermine le taux de comptage en 15 divisant la valeur finalement obtenue pour C par la valeur finalement obtenue pour Ta_tot. Selon un mode de réalisation particulier du procédé, objet de l'invention : - on choisit un seuil d'amplitude S, 20 - on détermine en outre l'amplitude Ai de chaque impulsion acquise Ii, 1 En variante, l'itération peut être effectuée après avoir acquis la totalité des impulsions I1r I2, IN et déterminé les valeurs D1, D2, ..., DN, T2, TN, et éventuellement les valeurs Al, A2, ..., AN des amplitudes des impulsions acquises. Et lorsque le signal électrique impulsionnel est un signal analogique, les valeurs Dl, D2, ..., DN, T2, ..., TN, et éventuellement les valeurs Al, A2, ..., AN des amplitudes des impulsions acquises, peuvent être déterminées par des moyens électroniques analogiques, après quoi elles sont numérisées en vue du traitement des impulsions acquises, ou bien le signal 10 électrique impulsionnel peut être numérisé, après quoi les valeurs Dlr D2, ..., DN, T2, ..., TN, et éventuellement les valeurs Al, A2, ..., AN des amplitudes des impulsions acquises, sont déterminées par des moyens électroniques numériques en vue du traitement des impulsions acquises. Selon un mode de réalisation particulier de l'invention, on acquiert le signal électrique impulsionnel à la sortie d'un dispositif capable d'engendrer ce signal électrique impulsionnel sous l'effet d'un phénomène physique. Le phénomène physique peut être un rayonnement, auquel cas le dispositif est un détecteur de ce rayonnement. En particulier, le phénomène physique peut être un rayonnement émis par un radionucléide, auquel cas le dispositif est un détecteur de ce rayonnement, la mesure du taux de comptage des impulsions permettant ainsi la mesure de l'activité du radionucléide. 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 et dans laquelle les informations nécessaires pour mettre en oeuvre un 11 procédé conforme à l'invention sont déterminées de façon analogique ; - la figure 2 est une vue schématique d'une autre chaîne de mesure qui est destinée à déterminer l'activité d'un radionucléide et dans laquelle les informations nécessaires pour mettre en oeuvre un procédé conforme à l'invention sont déterminées de façon numérique ; - la figure 3 est un organigramme qui illustre un exemple du procédé, objet de l'invention ; - la figure 4 illustre schématiquement une portion d'un signal impulsionnel que l'on peut traiter conformément à l'invention ; et - la figure 5 montre : en A, un signal comportant quatre impulsions dont les deux dernières forment un empilement ; en B, un signal de discrimination associé au signal montré en A ; en C, un signal représentant un temps mort imposé reconductible ; et en D, un signal réprésentant un temps mort total. EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS La chaîne de mesure, qui est schématiquement représentée sur la figure 1, est destinée à étudier une source radioactive 2 constituée par un radionucléide qui émet un rayonnement 4. La chaîne de mesure permet de mesurer l'activité de ce radionucléide et comprend successivement : un détecteur 6 pour détecter le rayonnement 4, 12 - un pré-amplificateur 8 pour amplifier le signal électrique analogique impulsionnel, fourni par le détecteur 6 lorsque ce dernier reçoit le rayonnement, - un dispositif électronique analogique 10 pour traiter (sous forme analogique) le signal (électrique) amplifié, fourni par le pré-amplificateur 8, afin de déterminer les informations de temps (durée de chaque impulsion sélectionnée du signal, et temps séparant deux impulsions successives) et éventuellement d'amplitude qui sont nécessaires à la mise en oeuvre d'un procédé conforme à l'invention, un numériseur 12 pour numériser les informations fournies par le dispositif 10 sous forme analogique, - un ordinateur 14 pour calculer l'activité du radionucléide à partir des informations ainsi numérisées et à l'aide d'un procédé conforme à l'invention, et - des moyens 16 pour afficher les résultats obtenus par l'ordinateur 14. Le dispositif 10 est réalisable par l'homme du métier à partir de modules électroniques commercialement disponibles. En particulier, on utilise par exemple des TAC c'est-à-dire des convertisseurs temps-amplitude (en anglais, time-to-amplitude converters) pour la mesure des temps. Une variante de la chaîne de mesure représentée sur la figure 1 est schématiquement illustrée par la figure 2. La chaîne que l'on voit sur cette figure 2 comprend un détecteur 18 et un pré- 13 amplificateur 20 qui ont les mêmes fonctions que le détecteur 6 et le pré-amplificateur 8 de la figure 1. La chaîne de la figure 2 comprend aussi : - un numériseur rapide 22 pour numériser le signal électrique impulsionnel total, fourni par le pré-amplificateur 20, et - un dispositif électronique numérique 24 pour traiter le signal numérique fourni par le numériseur 22, afin de déterminer les informations de temps (durée de chaque impulsion sélectionnée du signal et temps séparant deux impulsions successives) et éventuellement d'amplitude qui sont nécessaires à la mise en oeuvre d'un procédé conforme à l'invention. Le dispositif 24 est suivi d'un ordinateur 26 pour calculer l'activité du radionucléide à partir des informations fournies par le dispositif 24 et à l'aide d'un procédé conforme à l'invention. Et des moyens d'affichage 28 sont prévus pour afficher les résultats obtenus par l'ordinateur 26. On précise que le procédé, objet de l'invention, ne dépend pas des caractéristiques du détecteur 6 (figure 1) ou 18 (figure 2), qui se compose par exemple d'un scintillateur suivi d'un photomultiplicateur. En revenant à la figure 2, on précise aussi que le numériseur rapide 22 est adapté aux durées (typiquement 1 ps) des impulsions qu'il reçoit du pré- amplificateur 20 : il doit avoir une fréquence d'échantillonnage suffisante pour permettre la résolution temporelle du signal qu'il reçoit de ce pré- amplificateur. 14 A titre purement indicatif et nullement limitatif, la fréquence d'échantillonnage du numériseur 22 vaut 200 MHz. En d'autres termes, le signal qu'il reçoit est numérisé toutes les 5 ns par lui. Le dispositif 24 est également réalisable par l'homme du métier, à partir de composants programmables tels que des processeurs ou des FPGA c'est-à-dire des réseaux de portes programmables in situ (en anglais, field-programmable gate arrays). A titre d'exemple, le dispositif 24 de la figure 2 est un FPGA ; ce dernier est programmé de manière à prendre en charge la détection des dépassements de seuil du signal numérisé qui lui parvient du numériseur 22, ainsi que la mesure des durées, et éventuellement des amplitudes, qui ont été mentionnées plus haut et sur lesquelles on reviendra par la suite. L'algorithme qui est mis en oeuvre dans le procédé, objet de l'invention, nécessite les informations relatives à : - la durée de chaque impulsion, cette durée constituant le � temps de discrimination �, c'est-à-dire le temps pendant lequel le signal impulsionnel reste au dessus d'un seuil de déclenchement prédéfini, et - la durée de l'intervalle de temps entre chaque impulsion et l'impulsion précédente. Comme on l'a vu, ces informations peuvent être indifféremment fournies par un système analogique, utilisant par exemple un TAC, ou par un système numérique, utilisant par exemple un FPGA pour le 15 traitement, dès lors que ces informations sont disponibles sur un ordinateur ou, plus généralement, tout moyen de calcul capable d'implémenter l'algorithme. On peut également utiliser l'information relative à l'amplitude de chacune des impulsions pour faire une mesure d'activité au dessus d'un seuil d'énergie prédéfini. Les données temporelles (durées), mentionnées plus haut, constituent les informations minimales qui sont nécessaires pour effectuer une mesure d'activité avec traitement des temps morts conformément à l'invention. En travaillant sur cette base, on autorise le système en amont à compresser au maximum le flux de données initial qui est généré par le signal numérisé, ce qui assouplit les limitations qui sont liées aux transferts de données entre l'ensemble carte d'acquisition-système de mesure et l'ordinateur d'acquisition. On a mis en oeuvre le procédé en utilisant un détecteur de type cristal-puits NaI, avec des sources de 37Cs ; et l'on n'a alors observé aucune perte de comptage significative jusqu'à un taux de comptage de 105 coups par seconde. Le système utilisé pour cette mise en oeuvre associe un numériseur rapide 12 bits/200 MHz à un FPGA qui effectue la détection des dépassements de seuil ainsi que les mesures de temps et d'amplitude. Ces informations sont traitées sur l'ordinateur d'acquisition, par exemple de type PC, pour appliquer le traitement de temps mort et le calcul du taux de comptage. 16 Les principaux avantages procurés par le procédé, objet de la présente invention, sont l'amélioration des taux de comptage mesurables ainsi que la robustesse et l'adaptabilité à divers types de systèmes de détection-acquisition : ce procédé est efficace, la précision étant conservée jusqu'aux forts taux de comptage, et il est robuste car applicable à une large gamme de systèmes de mesure, quelque soit le temps de réponse de ces systèmes (temps mort minimal) ou les artefacts générés par ces derniers. L'algorithme utilisé dans un exemple du procédé, objet de l'invention, est schématiquement représenté sur la figure 3. Sur cette figure 3, i désigne une variable d'indentation pour la boucle utilisée dans l'algorithme, sur les différents évènements détectés ; C (comptage) désigne une variable contenant le nombre d'impulsions détectées en dehors du temps mort ; Di représente la durée de discrimination de l'impulsion Ii considérée (1 S ?) Si l'amplitude Ai dépasse ce seuil, elle est prise en compte dans l'analyse. Pour l'impulsion Ii, l'algorithme compare alors la valeur Ti (durée entre la fin de l'impulsion Ii_1 et le début de l'impulsion Ii) à la valeur du temps mort imposé Tm. Si Ti est supérieur à Tm, l'impulsion Ii n'est pas dans le temps mort imposé du fait de l'impulsion I. La variable de comptage est donc incrémentée (C=C+1) et la durée totale du temps actif est mise à jour, en ajoutant à la variable Ta_tot la 19 durée entre la fin du temps mort imposé, précédant l'impulsion Ii, et le début de l'impulsion Ii (Ta_tot=Ta_tot+Ti-Tm) . Et l'on passe à l'impulsion suivante (i=i+1). Et, si Ti est inférieur ou égal à Tm, on passe directement à l'impulsion suivante. Si l'amplitude Ai de l'impulsion Ii ne dépasse pas le seuil S (Ai Il convient en effet de noter qu'en général, dans toute chaîne de mesure utilisant des moyens électroniques analogiques, la partie de la chaîne destinée à détecter les impulsions comporte un discriminateur qui détecte le dépassement d'un seuil. Il n'est donc pas nécessaire d'utiliser un seuil supplémentaire. Toutefois, on peut souhaiter ajouter un tel seuil, supérieur à celui qui existe déjà, par exemple pour éliminer des impulsions parasites (bruit). Comme on l'a vu, la dernière étape de l'algorithme est l'incrémentation de l'indice i, afin 20 de passer à la prise en compte de l'impulsion suivante 11+1. On peut facilement ajouter un test supplémentaire après incrémentation de i, par exemple pour fixer un nombre total d'impulsions, détecter la fin d'un fichier rassemblant les données d'entrée, ou encore définir un temps actif maximum Tam ou un temps total maximum Ttm : par exemple, dans ce dernier cas, cela signifie que l'on ne � reboucle � pas (arrêt de l'incrémentation de i) si le temps total devient supérieur ou égal à Ttm. Le taux de comptage est finalement estimé (après que i a parcouru toutes les valeurs de 1 à N) en calculant le rapport C/Ta_tot, c'est-à-dire le nombre d'impulsions détectées hors temps mort, dont l'amplitude est supérieure à la variable de seuil S, divisé par la durée totale du temps actif. La figure 4 est un exemple-type d'une portion de signal impulsionnel, comprenant deux impulsions In et In+1. On voit les paramètres correspondants Dn, Dn+1 et Tn+l. Cette figure est relative à un traitement numérique du signal. Le temps t est porté en abscisse et l'amplitude a du signal est portée en ordonnée (unités arbitraires). La figure 5 est également relative à un traitement numérique du signal. Le temps t est porté en abscisse et l'amplitude a en ordonnée. Sur la partie A de cette figure, on voit une portion de signal-type, comprenant quatre impulsions, dont les deux dernières forment un � empilement �. 21 Sur la partie B, on voit un signal de discrimination, qui vaut 1 lorsque le signal représenté dans la partie A est au dessus du seuil de déclenchement S, et qui vaut 0 sinon ; ce signal de discrimination représente le temps mort minimal qui est imposé par le système. Sur la partie C, on voit le signal représentant le temps mort imposé reconductible Tm. Et, sur la partie D, on voit le signal représentant le temps mort total, c'est-à-dire le cumul du temps mort minimal du système et du temps mort imposé ; lorsque l'amplitude de ce signal vaut 1, il s'agit d'une période de temps mort ; et lorsque l'amplitude du signal vaut 0, il s'agit d'une période de temps actif. La présente invention n'est pas limitée à la mesure de l'activité d'un radionucléide. Elle s'applique aussi à la mesure du taux de comptage de tout rayonnement (photons ou particules) détecté par des moyens appropriés, aptes à fournir un signal électrique impulsionnel par interaction avec ce rayonnement. Plus généralement, l'invention s'applique à la mesure du taux de comptage d'impulsions électriques issues d'un dispositif apte à fournir de telles impulsions sous l'effet d'un phénomène physique. A titre d'exemple, on peut citer un détecteur de mouvements. L'invention s'applique même à la mesure du 30 taux de comptage des potentiels d'actions, c'est-à-dire des impulsions fournies par les neurones. REVENDICATIONS1. 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, caractérisé en ce que : - on choisit une valeur Tm de temps mort reconductible, - on acquiert N impulsions successives I1r 10 I2r IN d'un signal électrique impulsionnel, on détermine la durée Di de chaque impulsion acquise Ii, ainsi que la durée Ti de l'intervalle de temps entre l'impulsion précédente Ii_1 et l'impulsion Ii, Ti étant choisi égal à l'intervalle 15 de temps entre le début de l'acquisition et le début de la première impulsion, - on effectue une itération, en donnant successivement à l'indice i les valeurs 1 à N, et : • pour une valeur i de l'indice, on met à 20 jour la durée totale Ttot de l'acquisition, en ajoutant Ti et Di à la valeur que Ttot avait auparavant, Ttot étant initialement égal à 0, • on compare Ti à Tm, • et si Ti est inférieur ou égal à Tm, on 25 passe à l'indice suivant i+1, • ou si Ti est supérieur à Tm, on met à jour la durée du temps actif Ta_tot en ajoutant la différence Ti-Tm à la valeur que Ta_tot avait auparavant, Ta_tot étant initialement égal à 0, on met à jour une variable C 30 égale au nombre d'impulsions détectées en dehors du temps mort, en ajoutant 1 à la valeur que C avaitauparavant, C étant initialement égal à 0, et l'on passe à l'indice suivant i+1, et on détermine le taux de comptage en divisant la valeur finalement obtenue pour C par la valeur finalement obtenue pour Ta_tot. 2. Procédé selon la revendication 1, dans - on choisit un seuil d'amplitude S, - on détermine en outre l'amplitude Ai de chaque impulsion acquise Ti, 1 3. Procédé selon l'une quelconque des revendications 1 et 2, dans lequel on effectue un test supplémentaire après être passé à l'indice suivant, pour savoir si oui ou non : 25 on a traité un nombre prédéfini d'impulsions, - ou l'on est arrivé à la fin d'un fichier contenant les valeurs Dlr D2, ..., DN et T1, T2, ..., TN, et éventuellement les valeurs Al, A2, ..., AN des amplitudes 30 des impulsions acquises, 24 ou l'on a atteint un temps actif prédéfini, ou l'on a atteint un temps total d'acquisition prédéfini. 4. Procédé selon l'une quelconque des revendications 1 à 3, dans lequel on effectue l'itération au fur et à mesure de l'acquisition des impulsions I2r I2, IN et de la détermination des valeurs Dlr D2, ..., DN, T2, ..., TN, et éventuellement des valeurs Al, A2, AN des amplitudes des impulsions acquises. 5. Procédé selon l'une quelconque des 15 revendications 1 à 3, dans lequel on effectue l'itération après avoir acquis la totalité des impulsions I2r I2, IN et déterminé les valeurs D1, D2, DN, T2r ..., TN, et éventuellement les valeurs Al, A2, ..., AN des amplitudes des impulsions acquises. 20 6. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel le signal électrique impulsionnel est un signal analogique, les valeurs Dlr D2, ..., DN, T2, ..., TN, et éventuellement les valeurs Al, 25 A2, ..., AN des amplitudes des impulsions acquises, sont déterminées par des moyens électroniques analogiques puis sont numérisées en vue du traitement des impulsions acquises. 30 7. Procédé selon l'une quelconque des revendications 1 à 5, dans lequel le signal électrique 25 impulsionnel est un signal analogique, ce signal électrique impulsionnel est numérisé puis les valeurs Dlr D2, ..., DN, T2, ..., TN, et éventuellement les valeurs Al, A2, ..., AN des amplitudes des impulsions acquises, sont déterminées par des moyens électroniques numériques en vue du traitement des impulsions acquises. 8. Procédé selon l'une quelconque des revendications 1 à 7, dans lequel on acquiert le signal électrique impulsionnel à la sortie d'un dispositif (6, 18) capable d'engendrer ce signal électrique impulsionnel sous l'effet d'un phénomène physique (4). 9. Procédé selon la revendication 8, dans lequel le phénomène physique est un rayonnement(4) et le dispositif est un détecteur (6, 18) de ce rayonnement. 10. Procédé selon la revendication 8, dans lequel le phénomène physique est un rayonnement (4) émis par un radionucléide (2) et le dispositif est un détecteur (6, 18) de ce rayonnement, la mesure du taux de comptage des impulsions permettant ainsi la mesure de l'activité du radionucléide.