FR 3065814 A1 20181102 FR 1753811 A 20170428 1 PROCÃDÃ ET SYSTÃME DE CARACTÃRISATION D'UNE SOURCE DE RAYONNEMENT IONISANT Ã CARACTÃRISER DESCRIPTION DOMAINE TECHNIQUE La prÃ©sente invention concerne un procÃ©dÃ© et un systÃ¨me de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser, ainsi qu'un support de donnÃ©es lisible par un ordinateur sur lequel est enregistrÃ© un programme incluant les instructions pour exÃ©cuter des Ã©tapes du procÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser. L'invention s'applique Ã la dÃ©tection et Ã la caractÃ©risation de sources de rayonnement ionisant. La terminologie de Â« rayonnement ionisant Â» doit s'entendre ci-dessus et dans le reste de ce document dans son entendement le plus large en incluant ainsi aussi bien les rayonnements dits Â« indirectement ionisants Â», tels que les photons gamma y et les neutrons, que les rayonnement dits Â« ionisants Â» en tant que tels, tels que les particules bÃ©ta p (Ã©lectrons et positrons) et les particules alpha a (protons). ÃTAT DE LA TECHNIQUE ANTÃRIEURE Avec le contexte du terrorisme international et les risques liÃ©s Ã l'utilisation de matÃ©riaux radioactifs pour la fabrication de bombes dites Â« sales Â», le dÃ©veloppement d'outils permettant de dÃ©tecter et d'identifier ces matÃ©riaux radioactifs est devenu une prioritÃ©. Ces matÃ©riaux radioactifs sont gÃ©nÃ©ralement identifiables du fait qu'ils constituent des sources de rayonnement ionisant Ã©mettant des neutrons. Les mÃ©thodes usuelles pour dÃ©tecter un rayonnement ionisant comportant des neutrons requiÃ¨rent gÃ©nÃ©ralement l'utilisation de l'isotope 3 de l'hÃ©lium 3He et de ses capacitÃ©s Ã absorber les neutrons. NÃ©anmoins, avec la multiplication des besoins en dÃ©tection de matÃ©riau radioactif, il y a pÃ©nurie d'isotopes 3 de l'hÃ©lium3He. Il 3065814 2 est donc nÃ©cessaire de dÃ©velopper des mÃ©thodes alternatives pour la dÃ©tection de rayonnement ionisant n'utilisant pas l'isotope 3 de l'hÃ©lium 3He. Contrairement aux mÃ©thodes utilisant l'isotope 3 de l'hÃ©lium 3He, les mÃ©thodes alternatives prÃ©sentent l'inconvÃ©nient d'Ãªtre aussi bien sensibles aux neutrons 5 qu'aux rayonnements gamma. Or les photons gamma sont gÃ©nÃ©ralement prÃ©pondÃ©rants sur les neutrons et peuvent Ãªtre Ã©mis par des sources parfaitement lÃ©gitime, telles que par exemple les bananes. Parmi ces mÃ©thodes alternatives Ã l'utilisation de l'isotope 3 de l'hÃ©lium, celles basÃ©es sur les scintillateurs organiques, telles que dÃ©crites par G. H. V. Bertrand et 10 ses coauteurs dans leurs travaux publiÃ©s en 2015 dans la revue scientifique Â« Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment Â», volume 776 pages. 114 Ã 128, sont particuliÃ¨rement prometteuses. En effet, les scintillateurs organiques sont relativement peu coÃ»teux et ils permettent de dÃ©tecter l'ensemble des neutrons, sans se limiter aux 15 seuls neutrons thermiques. NÃ©anmoins, avec de tels scintillateurs, la discrimination entre les signaux issus de neutrons et ceux issus de photons gamma reste particuliÃ¨rement problÃ©matique. En effet, cette discrimination, si elle basÃ©e le temps de vol, nÃ©cessite l'utilisation deux scintillateurs distincts. Des rÃ©cents travaux, tels que ceux de A.Buffler et 20 ses coauteurs publiÃ©s en 2016 dans la revue scientifique Â« International Journal of Modern Physics Conference Series Â» volume 44 pages 1660228, proposent, afin de supprimer cette nÃ©cessitÃ© Ã faire appel Ã deux scintillateurs, d'utiliser une discrimination par forme d'impulsion assistÃ©e par ordinateur. Une telle discrimination requiÃ¨re un Ã©talonnage Ã partir de signaux produits par l'interaction uniquement de neutrons qui sont 25 relativement difficiles, voire impossibles, Ã obtenir, si ce n'est par une prÃ©-sÃ©paration avec une technique basÃ©e sur le temps de vol. Ainsi l'utilisation de scintillateur organique pour caractÃ©riser des sources de rayonnement ionisant, notamment pour discriminer les signaux produits par l'interaction des neutrons, requiÃ¨re une mÃ©thode de discrimination simple, c'est-Ã -dire ne 3065814 3 requÃ©rant pas l'utilisation de deux scintillateurs en parallÃ¨le ou un Ã©talonnage Ã partir de signaux produits par l'interaction uniquement de neutrons. On notera qu'une telle problÃ©matique ne se limite pas aux seuls scintillateurs organiques, celle-ci s'appliquant Ã©galement aux scintillateur inorganique, et 5 est Ã©galement observÃ©e pour la discrimination d'autres types de rayonnement ionisant, tels que la discrimination entre des sources de rayonnements ionisants dont le rayonnement ionisant comporte des photons gamma y et des particules beta F3 ou des particules alpha a. EXPOSÃ DE L'INVENTION 10 L'invention vise Ã rÃ©soudre au moins partiellement les inconvÃ©nients ci- dessus et a ainsi notamment pour but de fournir un procÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser et dÃ©terminer ou non la prÃ©sence de neutrons Ã©mis par ladite source, le procÃ©dÃ© ne requÃ©rant qu'un seul scintillateur et ne nÃ©cessitant pas un accÃ¨s Ã un rayonnement ionisant comportant uniquement des 15 neutrons, ou d'effectuer au prÃ©alable une sÃ©paration des neutrons. A cet effet, l'invention a pour objet un procÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser dont le rayonnement ionisant comporte un premier type de particules et est susceptible de comporter un deuxiÃ¨me type de particules, le procÃ©dÃ© comportant les Ã©tapes suivantes : 20 - fourniture de N signaux temporels d'impulsions de rÃ©fÃ©rence, chaque signal temporel d'impulsion de rÃ©fÃ©rence correspondant Ã la dÃ©tection du premier type de particules par un scintillateur, N Ã©tant un entier naturel non nul ; - choix d'au moins une grandeur discriminante ; - pour chaque signal temporel d'impulsion de rÃ©fÃ©rence, extraction de la 25 valeur de chaque grandeur discriminante choisie pour former une valeur de rÃ©fÃ©rence correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de rÃ©fÃ©rence ; - acquisition de M signaux temporels d'impulsions acquis, chaque signal temporel d'impulsion acquis correspondant Ã la dÃ©tection, par le scintillateur, d'un 3065814 4 photon issu de la source de rayonnement ionisant Ã caractÃ©riser, M Ã©tant un entier naturel non nul ; - pour chaque signal temporel d'impulsion acquis, extraction de la valeur de chaque grandeur discriminante choisie pour former une valeur caractÃ©ristique 5 correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion acquis ; - calcul d'une premiÃ¨re valeur caractÃ©ristique de dissemblance Ã partir d'un premier calcul de dissemblance entre une matrice de rÃ©fÃ©rence et une matrice caractÃ©ristique, la matrice de rÃ©fÃ©rence Ã©tant formÃ©e par les valeurs de rÃ©fÃ©rence, la 10 matrice caractÃ©ristique Ã©tant formÃ©e par les valeurs caractÃ©ristiques ; - comparaison de la premiÃ¨re valeur caractÃ©ristique de dissemblance calculÃ©e Ã un premier seuil prÃ©dÃ©terminÃ© ; - caractÃ©risation de la prÃ©sence du deuxiÃ¨me type de particules s dans un rayonnement Ã©mis par la source de rayonnement ionisant Ã caractÃ©riser en fonction 15 du rÃ©sultat de la comparaison entre la premiÃ¨re valeur caractÃ©ristique de dissemblance calculÃ©e et le premier seuil prÃ©dÃ©terminÃ©, dans lequel l'Ã©tape de choix d'au moins une grandeur discriminante comporte le choix de R grandeurs discriminantes, R Ã©tant un entier naturel non nul, et le premier de dissemblance entre la matrice de rÃ©fÃ©rence et la matrice caractÃ©ristique 20 comporte les sous-Ã©tapes suivantes : - pour chaque grandeur discriminante identifiÃ©e par un indice i, subdivision d'un intervalle reprÃ©sentatif des valeurs prises par la grandeur discriminante en k, sous-intervalles, k, Ã©tant un entier naturel non nul ; - calcul, pour chacune de la matrice de rÃ©fÃ©rence et la matrice 25 caractÃ©ristique, de la valeur de Mi k, coefficients 0, chaque coefficient 0 Ã©tant associÃ© Ã R sous-intervalles, chacun desdits R sous-intervalles Ã©tant associÃ© Ã une grandeur discriminante respective et Ã©tant Ã©gal Ã un sous-intervalle de ladite grandeur discriminante respective, un premier coefficient 01 et un deuxiÃ¨me coefficient 02 distincts Ã©tant tels qu'il existe au moins une grandeur caractÃ©ristique pour laquelle le sous- intervalle correspondant du premier coefficient 01 est diffÃ©rent du sous-intervalle 3065814 5 correspondant du deuxiÃ¨me coefficient 02, la valeur de chacun des coefficients 0 Ã©tant calculÃ©e comme le nombre de signaux temporels d'impulsions de rÃ©fÃ©rence, respectivement de signaux temporels d'impulsions acquis, dont les grandeurs discriminantes sont dans les P sous-intervalles associÃ©s audit coefficient, les valeur de 5 coefficients 0 calculÃ©es pour l'une de la matrice de rÃ©fÃ©rence et de la matrice caractÃ©ristique Ã©tant multipliÃ©es par un facteur correctif correspondant au ration du nombre signaux temporels d'impulsion entre la matrice de rÃ©fÃ©rence et la matrice caractÃ©ristique ; - calcul de la premiÃ¨re valeur caractÃ©ristique de dissemblance Ã partir 10 d'une diffÃ©rence en valeur absolue mettant en oeuvre les coefficients 0 calculÃ©s Ã partir de la matrice de rÃ©fÃ©rence et de la matrice caractÃ©ristique et de R sommes successives, chacune des R sommes Ã©tant associÃ©e Ã une grandeur discriminante respective, et ladite diffÃ©rence se faisant directement entre chaque coefficient 0 calculÃ© sur la base de la matrice de rÃ©fÃ©rence et le coefficient 0 calculÃ© sur la base de la matrice caractÃ©ristique 15 correspondant aux mÃªmes R sous-intervalles, la premiÃ¨re somme se faisant sur la base du rÃ©sultat de ladite diffÃ©rence et les R-1 sommes suivantes sur la base du rÃ©sultat de la somme la prÃ©cÃ©dent directement. En effet, les inventeurs ont dÃ©couvert qu'en utilisant au moins une grandeur de discriminante et le calcul d'une valeur caractÃ©ristique de dissemblance entre 20 une matrice de valeurs de grandeur de discriminante calculÃ©es pour une source de rayonnement ionisant constituÃ© de photons gamma et une matrice de valeurs de grandeur de discriminante calculÃ©es pour la source de rayonnement ionisant Ã caractÃ©riser, il Ã©tait possible de discriminer la prÃ©sence de signaux temporels d'impulsions issues de neutrons. 25 De plus l'utilisation d'un calcul dissemblance basÃ© sur le dÃ©coupage en sous-intervalles de chacune grandeur discriminante suivit de N sommes permet d'identifier l'ensemble des dissemblances entre la population de signaux temporels d'impulsions de rÃ©fÃ©rence et la population de signaux temporels d'impulsions acquis. La discrimination est donc particuliÃ¨rement amÃ©liorÃ©e vis-Ã -vis des procÃ©dÃ©s de 30 caractÃ©risation selon l'art antÃ©rieur. 3065814 6 On peut noter que ainsi, avec un tel procÃ©dÃ©, dans le cadre d'une application prÃ©fÃ©rÃ©e de l'invention dans lequel le premier type de particules est celui des photons gamma y et le deuxiÃ¨me type de particules est celui des neutrons, il est possible de caractÃ©riser une source de rayonnement ionisant Ã caractÃ©riser afin de dÃ©terminer si 5 cette source de rayonnement ionisante Ã©met des neutrons. Le procÃ©dÃ© selon l'invention ne nÃ©cessite donc qu'un seul scintillateur et est rÃ©alisÃ© Ã partir d'un Ã©talonnage basÃ© sur une source ionisante de rÃ©fÃ©rence n'Ã©mettant que des photons gamma. On peut noter que par valeur caractÃ©ristique de dissemblance, il doit Ãªtre entendu une valeur caractÃ©risant la diffÃ©rence entre les N signaux temporels 10 d'impulsion de rÃ©fÃ©rence et les M signaux temporels d'impulsion acquis comportant les Ã©tapes de calcul prÃ©citÃ©, cette valeur caractÃ©ristique pouvant Ãªtre une valeur caractÃ©ristique de dissemblance en tant que telle, telle qu'une diffÃ©rence statistique entre les valeurs de rÃ©fÃ©rence et les valeurs caractÃ©ristique, ou une valeur de vraisemblance, telle que l'inverse d'une diffÃ©rence statistique entre les valeurs de 15 rÃ©fÃ©rence et les valeurs caractÃ©ristique. Bien entendu, l'homme du mÃ©tier est Ã mÃªme d'utiliser aussi bien une valeur caractÃ©ristique de dissemblance qu'une valeur de vraisemblance, ces valeurs Ã©tant gÃ©nÃ©ralement liÃ©es l'une Ã l'autre par un simple calcul d'inversion. Dans le cas oÃ¹ la valeur caractÃ©ristique de dissemblance est une valeur 20 de vraissemblance, le calcul de la valeur caractÃ©ristique de dissemblance peut Ã©galement comporter une opÃ©ration d'inversion du rÃ©sultat obtenu aprÃ¨s la NIÃ¨me somme. On notera qu'un le terme de matrice de rÃ©fÃ©rence et de matrice caractÃ©ristique doivent s'entendre de maniÃ¨re large et en tant que tableau de valeurs numÃ©riques quel qu'en soit la reprÃ©sentation rÃ©ellement utilisÃ© pour le traitement des 25 donnÃ©es pour mettre en oeuvre le procÃ©dÃ© selon l'invention. Ainsi, il est parfaitement envisageable, sans que l'on sorte du cadre de l'invention, que les valeurs de rÃ©fÃ©rence, respectivement les valeurs caractÃ©ristiques, soient manipulÃ©es sous la forme de N vecteurs, respectivement M vecteurs, chaque vecteur comportant l'ensemble des valeurs de rÃ©fÃ©rence, respectivement des valeurs caractÃ©ristiques, extraites pour un signal 30 temporel d'impulsion donnÃ©. 3065814 7 Les Ã©tapes d'extraction de la valeur de chaque grandeur discriminante, que ce soit pour les signaux temporels d'impulsions de rÃ©fÃ©rence et pour les signaux temporels d'impulsion acquis, peuvent comprendre au moins une sous-Ã©tape choisie dans les sous-Ã©tapes suivantes : 5 normalisation de chaque signal temporel d'impulsion ; rÃ©duction du bruit Ã©lectronique de chaque signal temporel d'impulsion. De telles sous-Ã©tapes permettent d'obtenir une comparaison aisÃ©e entre les signaux temporels d'impulsions et donc un calcul d'une valeur 10 caractÃ©ristique de dissemblance parfaitement reprÃ©sentatif de la dissemblance entre les signaux temporels d'impulsion de rÃ©fÃ©rence et ceux acquis. L'au moins une grandeur discriminante peut Ãªtre sÃ©lectionnÃ©e parmi : - une intÃ©grale d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis; 15 - un temps caractÃ©ristique du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis; - une valeur, en au moins un instant, du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis; - une fonction de la variance d'au moins une portion du signal temporel 20 d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis; - une valeur de dÃ©rivÃ©e du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis. De telles grandeurs discriminantes peuvent permettre de mettre en Ã©vidence les caractÃ©ristiques distinctives entre les signaux temporels d'impulsion 25 signature d'une interaction d'une particule du premier type et ceux d'une interaction d'une particule du deuxiÃ¨me type. En effet, l'intÃ©grale d'au moins une portion peut Ãªtre par exemple reprÃ©sentative d'un pic secondaire prÃ©sent uniquement dans les signaux temporels d'impulsion obtenus pour les particules du deuxiÃ¨me type alors que le temps caractÃ©ristique peut Ãªtre reprÃ©sentative d'un temps d'interaction plus important pour ces 30 mÃªme particules du deuxiÃ¨me type. 3065814 8 La terminologie de Â« valeur de dÃ©rivÃ©e d'un signal temporel d'impulsion Â» doit s'entendre de maniÃ¨re large en englobant ainsi, par exemple, les valeurs de dÃ©rivÃ©e moyenne sur une portion du signal temporel d'impulsion, les valeurs de dÃ©rivÃ©e maximale du signal temporel d'impulsion, et les valeurs de dÃ©rivÃ©e du signal 5 temporel d'impulsion Ã un instant donnÃ©. L'Ã©tape de choix d'au moins une grandeur discriminante peut Ãªtre postÃ©rieure Ã l'Ã©tape de fourniture des N signaux temporels d'impulsions de rÃ©fÃ©rence et Ã l'Ã©tape de normalisation de chaque signal temporel d'impulsion de rÃ©fÃ©rence, l'Ã©tape de choix d'au moins une grandeur discriminante comportant les 10 sous-Ã©tapes suivantes : - fourniture de P signaux temporels d'impulsions de comparaison, chaque signal temporel d'impulsion de comparaison correspondant Ã la dÃ©tection, par le scintillateur, d'une particule du premier et ou deuxiÃ¨me type de particules Ã©mis par une source de rayonnement ionisant Ã©mettant un rayonnement ionisant comportant des 15 particules du premier type et des particules du deuxiÃ¨me type, P Ã©tant un entier naturel non nul ; - construction d'une pluralitÃ© d'ensembles de grandeurs discriminantes ; - pour chaque signal temporel d'impulsion de rÃ©fÃ©rence, et pour chaque ensemble de grandeurs discriminantes, extraction de la valeur de chaque grandeur 20 discriminante dudit ensemble de grandeurs discriminantes pour former une premiÃ¨re valeur de comparaison correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de rÃ©fÃ©rence ; - pour chaque signal temporel d'impulsion de comparaison, et pour chaque ensemble de grandeurs discriminantes, calcul de la valeur de chaque grandeur 25 discriminante dudit ensemble de grandeurs discriminantes pour former une deuxiÃ¨me valeur de comparaison ; - pour chaque ensemble de grandeurs discriminantes, calcul d'une deuxiÃ¨me valeur caractÃ©ristique de dissemblance Ã partir d'un deuxiÃ¨me calcul de dissemblance entre une premiÃ¨re matrice constituÃ©e des premiÃ¨res valeurs de 30 comparaison extraites pour ledit ensemble de grandeurs discriminantes et une deuxiÃ¨me 3065814 9 matrice constituÃ©e des deuxiÃ¨mes valeurs de comparaison extraites pour ledit ensemble de grandeurs discriminantes ; - choix d'au moins une grandeur discriminante correspondant Ã l'un des ensembles de grandeurs discriminantes, le choix Ã©tant rÃ©alisÃ© sur la base des deuxiÃ¨mes 5 valeurs caractÃ©ristique de dissemblance calculÃ©e pour chacun des ensembles de grandeurs discriminantes. Lors de la sous Ã©tape de choix d'au moins une grandeur discriminante correspondant Ã l'un des ensembles de grandeurs discriminantes, l'au moins une grandeur discriminante correspondant choisie peut Ãªtre celle qui correspond Ã l'ensemble 10 de grandeurs discriminantes pour lequel la deuxiÃ¨me valeur caractÃ©ristique de dissemblance prÃ©sente un extrÃ©mum caractÃ©ristique d'une dissemblance maximale. Avec de telles sous Ã©tape de choix d'au moins une grandeur discriminante, il est possible d'optimiser la discrimination, puisqu'il est possible de baser le choix de cette ou ces grandeurs discriminantes sur une maximisation de la grandeur 15 discriminante. L'Ã©tape de choix d'au moins une grandeur discriminante peut comporter en outre les sous-Ã©tapes suivantes : - fourniture de Q signaux temporels d'impulsions de test, chaque signal temporel d'impulsion de test correspondant Ã la dÃ©tection, par le scintillateur, d'une 20 particule du premier type Ã©mis par une source de rayonnement ionisant Ã©mettant un rayonnement ionisant comportant des particules du premier type, Q Ã©tant un entier naturel non nul ; - pour chaque signal temporel d'impulsion de test, et pour chaque ensemble de grandeurs discriminantes, extraction de la valeur de chaque grandeur 25 discriminante dudit ensemble de grandeurs discriminantes pour former une troisiÃ¨me valeur de comparaison correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de test ; - pour chaque ensemble de grandeurs discriminantes, calcul d'une troisiÃ¨me valeur caractÃ©ristique de dissemblance Ã partir du deuxiÃ¨me calcul de 30 dissemblance entre une premiÃ¨re matrice constituÃ©e des premiÃ¨res valeurs de 3065814 10 comparaison extraites pour ledit ensemble de grandeurs discriminantes et une deuxiÃ¨me matrice constituÃ©e des troisiÃ¨me valeurs de comparaison extraites pour ledit ensemble de grandeurs discriminantes ; - comparaison de la troisiÃ¨me valeur caractÃ©ristique de dissemblance 5 calculÃ©e pour chacun des ensembles de grandeurs discriminantes Ã un deuxiÃ¨me seuil prÃ©dÃ©terminÃ© ; - sÃ©lection parmi les ensembles de grandeurs discriminantes de ceux qui prÃ©sentent une dissemblance minimale Ã partir de ladite comparaison entre la troisiÃ¨me valeur caractÃ©ristique de dissemblance calculÃ©e pour chacun des ensembles de grandeurs 10 discriminantes Ã un deuxiÃ¨me seuil prÃ©dÃ©terminÃ© ; - choix d'au moins une grandeur discriminante correspondant Ã l'un des ensembles de grandeurs discriminantes, l'au moins une grandeur discriminante correspondant choisie Ã©tant celle correspondant Ã l'ensemble de grandeurs discriminantes parmi les ensembles de grandeurs discriminantes sÃ©lectionnÃ©s dont la 15 deuxiÃ¨me valeur caractÃ©ristique de dissemblance prÃ©sente un extrÃ©mum caractÃ©ristique d'une dissemblance maximale. Avec de telles sous Ã©tape de choix d'au moins une grandeur discriminante, il est possible d'optimiser la discrimination sans risque de dÃ©tection erronÃ©e, puisqu'il est possible de baser le choix de cette ou ces grandeurs discriminantes 20 sur une maximisation de la grandeur discriminante tout en limitant les risques de faux positif grÃ¢ce au calcul de la troisiÃ¨me valeur caractÃ©ristique de dissemblance. Au moins une grandeur discriminante peut Ãªtre notÃ©e (lot et calculÃ©e selon : ti Qt0t= ft.s(t) dt oÃ¹ s est le signal temporel d'impulsion de rÃ©fÃ©rence ou le signal 25 temporel d', respectivement le signal temporel d'impulsion acquis; tmax est la coordonnÃ©e temporelle du maximum du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis; et t1 est une coordonnÃ©e temporelle calculÃ©e selon : 3065814 11 t1= tmax+p.tHm oÃ¹ p est un rÃ©el strictement positif ; et tHM est la coordonnÃ©e temporelle, ultÃ©rieure Ã tmax, pour laquelle la valeur du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, vaut 0,5. 5 Une telle valeur (lot permet de discriminer les signaux temporels d'impulsion prÃ©sentant une interaction plus forte, ou plus faible, vis-Ã -vis de ceux des signaux temporels d'impulsions. Au moins une grandeur discriminante peut Ãªtre notÃ©e Qtad et calculÃ©e selon : Qtall= f dt 10 oÃ¹ s est le signal temporel d'impulsion de rÃ©fÃ©rence, respectivement le signal temporel d'impulsion acquis; to est une coordonnÃ©e temporelle calculÃ©e selon : tmax+a.tHm oÃ¹ a est un rÃ©el strictement positif ; tmax est la coordonnÃ©e temporelle du maximum du 15 signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis; et tHM est la coordonnÃ©e temporelle, ultÃ©rieure Ã tmax, pour laquelle la valeur du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, vaut 0,5 ; 20 et t1 est une coordonnÃ©e temporelle calculÃ©e selon : t1= tmax+13.tHm oÃ¹ p est un rÃ©el strictement positif. Une telle valeur Qtad permet de discriminer les signaux temporels d'impulsion prÃ©sentant une interaction retardÃ©e vis-Ã -vis de ceux des signaux temporels d'impulsions. 3065814 12 Chaque ensemble de grandeurs discriminantes comporte l'un au moins de (lot et Qtail, les ensembles se distinguant deux Ã deux par au moins l'une des valeurs a et (3. De cette maniÃ¨re, il est possible de dÃ©finir de maniÃ¨re optimisÃ©e les 5 bornes des sommes mis en oeuvre pour le calcul de la valeur des grandeurs discriminante et ainsi la valeur caractÃ©ristique de dissemblance qui sera calculÃ©e. L'au moins une grandeur discriminante peut Ãªtre notÃ©e D et est la variance de la portion du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement le signal temporel d'impulsion acquis, comprise entre tic, et t1, oÃ¹ : 10 tic, est une coordonnÃ©e temporelle calculÃ©e selon : -4= tmax+a.tHm oÃ¹ a est un rÃ©el strictement positif ; tmax est la coordonnÃ©e temporelle du maximum du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis; et 15 tHM est la coordonnÃ©e temporelle, ultÃ©rieure Ã tmax, pour laquelle la valeur du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, vaut 0,5 ; t1 est une coordonnÃ©e temporelle calculÃ©e selon : t1= tmax+13.tHm oÃ¹ p est un rÃ©el strictement positif. 20 Une telle grandeur de discriminante permet de caractÃ©riser la prÃ©sence de signaux temporels d'impulsions se diffÃ©renciant des signaux temporels d'impulsions de rÃ©fÃ©rence. Cette grandeur discriminante est donc particuliÃ¨rement adaptÃ©e pour mettre en Ã©vidence une population de deux types. Chaque ensemble de grandeurs discriminantes peut comporter D, les 25 ensembles se distinguant deux Ã deux par au moins l'une des valeurs a et R. Ainsi, il est possible de sÃ©lectionnÃ© la grandeur discriminante D permettant de fournir une grandeur de dissemblance parfaitement caractÃ©ristique de la prÃ©sence de particules du deuxiÃ¨me type. 3065814 13 Lors de l'Ã©tape de construction de la pluralitÃ© d'ensembles de grandeurs discriminantes, un ensemble donnÃ© de grandeurs discriminantes peut comporter S grandeurs discriminantes, S Ã©tant un entier naturel non nul, et le deuxiÃ¨me calcul de dissemblance entre une premiÃ¨re matrice et une deuxiÃ¨me matrice, pour l'ensemble 5 donnÃ© de S grandeurs discriminantes, peut comporter les sous-Ã©tapes suivantes : - pour chaque grandeur discriminante identifiÃ©e par un indice j, subdivision d'un intervalle reprÃ©sentatif des valeurs prises par la grandeur discriminante en k, sous-intervalles, k, Ã©tant un entier naturel non nul ; - calcul, pour chacune de la premiÃ¨re matrice et de la deuxiÃ¨me matrice, 10 de la valeur de Tls1kj coefficients 6, chaque coefficient 6 Ã©tant associÃ© Ã S sous- intervalles, chacun desdits S sous-intervalles Ã©tant associÃ© Ã une grandeur discriminante respective et Ã©tant Ã©gal Ã un sous-intervalle de ladite grandeur discriminante respective, un premier coefficient 61 et un deuxiÃ¨me coefficient 62 distincts Ã©tant tels qu'il existe au moins une grandeur caractÃ©ristique pour laquelle le sous-intervalle correspondant du 15 premier coefficient 61 est diffÃ©rent du sous-intervalle correspondant du deuxiÃ¨me coefficient 62, la valeur de chacun des coefficients 6 Ã©tant calculÃ©e comme le nombre de signaux temporels d'impulsions de rÃ©fÃ©rence, respectivement de signaux temporels d'impulsions acquis, dont les grandeurs discriminantes sont dans les S sous-intervalles associÃ©s audit coefficient, les valeurs de coefficients 6 calculÃ©s pour l'une de la premiÃ¨re 20 matrice et de la deuxiÃ¨me matrice Ã©tant multipliÃ©es par un facteur correctif correspondant au ration du nombre de signaux temporels d'impulsions entre la premiÃ¨re matrice et la deuxiÃ¨me matrice ; - calcul d'une valeur caractÃ©ristique de dissemblance Ã partir d'une diffÃ©rence en valeur absolue mettant en oeuvre les coefficients 6 calculÃ©s Ã partir de la 25 premiÃ¨re matrice et de la deuxiÃ¨me matrice et de S sommes successives, chacune des S sommes Ã©tant associÃ©e Ã une grandeur discriminante respective, et ladite diffÃ©rence se faisant directement entre chaque coefficient 6 calculÃ© sur la base de la premiÃ¨re matrice et le coefficient 6 calculÃ© sur la base de la deuxiÃ¨me matrice correspondant au mÃªme R sous-intervalles, la premiÃ¨re somme se faisant sur la base du rÃ©sultat de ladite diffÃ©rence 3065814 14 et les S-1 sommes suivantes sur la base du rÃ©sultat de la somme la prÃ©cÃ©dent directement. Un tel deuxiÃ¨me calcul d'une valeur caractÃ©ristique de dissemblance permet d'obtenir une somme de diffÃ©rences correspondant Ã chaque sous-intervalle des 5 grandeurs discriminantes. Il est donc aisÃ© avec un tel calcul de repÃ©rer la prÃ©sence d'une deuxiÃ¨me population de type de particules, telle que les signaux temporels d'impulsion des particules du deuxiÃ¨me type, vis-Ã -vis de la premiÃ¨re population homogÃ¨ne que constitue les signaux temporels d'impulsion de rÃ©fÃ©rence qui ne corresponde qu'Ã un seul type de particules, Ã savoir les particules du premier type. 10 Il peut Ãªtre en outre prÃ©vu une Ã©tape d'estimation du nombre et/ou d'un taux d'impulsions correspondant au deuxiÃ¨me type de particules. Le premier type de particules et le deuxiÃ¨me type de particules peuvent Ãªtre sÃ©lectionnÃ©s parmi le groupe comprenant : - les photons gamma en tant que premier type de particules et les 15 neutrons en tant que deuxiÃ¨me type de particules, - les photons gamma en tant que premier type de particules et les particules beta en tant que deuxiÃ¨me type de particules, - les photons gamma en tant que premier type de particules et les particules alpha en tant que deuxiÃ¨me type de particules, 20 - les particules beta en tant que premier type de particules et les particules alpha en tant que deuxiÃ¨me type de particules. L'invention concerne en outre un systÃ¨me de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser comportant : - un scintillateur, 25 - un photomultiplicateur couplÃ© au scintillateur, et - une unitÃ© de traitement adaptÃ©e pour traiter les signaux temporels d'impulsions fournis par le photomultiplicateur et obtenus lors de l'absorption de rayonnements ionisants par le scintillateur, le systÃ¨me de caractÃ©risation Ã©tant en ce que l'unitÃ© de traitement est 30 configurÃ©e pour mettre en oeuvre un procÃ©dÃ© de caractÃ©risation selon l'invention. 3065814 15 Un tel systÃ¨me de caractÃ©risation prÃ©sente les avantages liÃ©s Ã la mise en oeuvre du procÃ©dÃ© selon l'invention. L'invention concerne en outre un support de donnÃ©es lisible par un ordinateur sur lequel est enregistrÃ© un programme incluant les instructions pour exÃ©cuter 5 des Ã©tapes du procÃ©dÃ© selon l'invention. Un tel support de donnÃ©es prÃ©sente les avantages liÃ©s Ã la mise en oeuvre du procÃ©dÃ© selon l'invention. BRÃVE DESCRIPTION DES DESSINS L'invention sera mieux comprise Ã l'aide de la description qui va suivre, 10 donnÃ©e uniquement Ã titre d'exemple non limitatif et faite en se rÃ©fÃ©rant aux dessins annexÃ©s sur lesquels : - la figure 1 illustre une vue schÃ©matique d'un systÃ¨me de caractÃ©risation selon l'invention, - la figure 2 est un graphique de deux signaux temporels d'impulsions 15 moyennÃ©s obtenues avec un scintillateur organique, le premier signal temporel d'impulsion moyennÃ© ayant Ã©tÃ© obtenu avec une source de rayonnement ionisant n'Ã©mettant que des photons gamma, le deuxiÃ¨me signal temporel d'impulsion moyennÃ© ayant Ã©tÃ© obtenu avec une source de rayonnement ionisant Ã©mettant des photons gamma et des neutrons, 20 - la figure 3 est un graphique d'un signal temporel d'impulsion obtenue pour un scintillateur, le graphique figurant des valeurs utilisÃ©es pour le calcul de grandeurs discriminantes selon l'invention, - les figures 4A et 4B sont deux graphiques figurant la rÃ©partition selon deux grandeurs discriminantes de respectivement un ensemble d'impulsions issu d'une 25 source de rayonnement ionisant comportant uniquement des photons gamma et un ensemble d'impulsions issu d'une source de rayonnement ionisant comportant des neutrons et des photons gamma, - les figures 5A et 5B illustrent graphiquement les diffÃ©rences sur lesquelles se base le calcul intermÃ©diaire d'une valeur caractÃ©ristique de dissemblance 3065814 16 selon l'invention ceci pour respectivement une grandeur discriminante (lot et une grandeur discriminante Qtai, - les figures 6A et 6B sont des graphiques illustrant la variation d'une valeur caractÃ©ristique de dissemblance pour deux grandeurs discriminantes entre une 5 matrice de signaux temporels d'impulsions de rÃ©fÃ©rence produits par l'interaction de photons gamma et respectivement une matrice de signaux temporels d'impulsions de comparaison produits par l'interaction de neutrons et de photons gamma et un matrice de signaux temporels d'impulsions de comparaison produits par l'interaction uniquement de photons gamma, 10 - la figure 7 illustre le rÃ©sultat d'un calcul intermÃ©diaire pour le calcul d'une valeur caractÃ©ristique de discrimination Ã partir d'une grandeur discriminante de variance. Les diffÃ©rentes possibilitÃ©s (variantes et modes de rÃ©alisation) doivent Ãªtre comprises comme n'Ã©tant pas exclusives les unes des autres et peuvent se combiner 15 entre elles. EXPOSÃ DÃTAILLÃ DE MODES DE RÃALISATION PARTICULIERS La figure 1 illustre un systÃ¨me de caractÃ©risation 10 d'une source de rayonnement ionisant 20 Ã caractÃ©riser dont le rayonnement ionisant comporte un premier type de particules et est susceptible de comporter un deuxiÃ¨me type de 20 particules selon l'invention. Si dans le cadre de cet exposÃ© dÃ©taillÃ© seul le cas d'un premier type de particules qui est les photons gamma et d'un deuxiÃ¨me type de particules qui est les neutrons est dÃ©crit, l'invention ne se limite bien entendu pas Ã cette seule possibilitÃ©. En effet l'homme du mÃ©tier est parfaitement Ã mÃªme Ã partir de la prÃ©sente divulgation 25 d'appliquer l'invention Ã la discrimination d'autres types de particules prÃ©sents dans les rayonnements ionisants. L'invention vise ainsi Ã©galement plus particuliÃ¨rement la discrimination de : - sources de rayonnement ionisant comprenant des photons gamma y et des particules beta p vis-Ã -vis d'une source de rayonnement ionisant de rÃ©fÃ©rence 3065814 17 dont le rayonnement ne comporte que des photons gamma y, le premier type de particules Ã©tant alors les photons gamma y et le deuxiÃ¨me type de particules Ã©tant les particules beta p ; - sources de rayonnement ionisant comprenant des photons gamma 5 y et des particules alpha a vis-Ã -vis d'une source de rayonnement ionisant de rÃ©fÃ©rence dont le rayonnement ne comporte que des photons gamma y, le premier type de particules Ã©tant alors les photons gamma y et le deuxiÃ¨me type de particules Ã©tant les particules alpha a ; - sources de rayonnement ionisant comprenant des particules beta p et 10 des particules alpha a vis-Ã -vis d'une source de rayonnement ionisant de rÃ©fÃ©rence dont le rayonnement ne comporte que des particules beta 13; , le premier type de particules Ã©tant alors les particules beta p et le deuxiÃ¨me type de particules Ã©tant les particules alpha a. Un tel systÃ¨me de caractÃ©risation 10 comporte : 15 un scintillateur organique 11, un photomultiplicateur 12 couplÃ© au scintillateur organiques 11, et une unitÃ© de traitement 13 adaptÃ©e pour traiter les signaux temporels d'impulsions fourni par le photomultiplicateur 12 et obtenus lors de l'absorption de rayonnements ionisants par le scintillateur organique 11. 20 On peut noter que si l'invention vise plus particuliÃ¨rement les systÃ¨mes de caractÃ©risation 1 Ã scintillateur organique, elle ne se limite pas Ã ces seuls systÃ¨mes de caractÃ©risation et concerne Ã©galement les systÃ¨mes de caractÃ©risation Ã scintillateur inorganique tels que les scintillateurs Ã semiconducteur. Ainsi selon une telle possibilitÃ© non illustrÃ©, le systÃ¨me de caractÃ©risation peut ne pas comporter de photomultiplicateur 25 et peut comporter un systÃ¨me d'amplification du signal issu du scintillateur Ã semiconducteur autre qu'un photomultiplicateur. Ainsi d'une maniÃ¨re gÃ©nÃ©rale le systÃ¨me de caractÃ©risation comporte un scintillateur, tel que le scintillateur organique 11, un systÃ¨me d'amplification, tel que le photomultiplicateur 12, et une unitÃ© de traitement adaptÃ©e pour traiter les signaux 3065814 18 temporels d'impulsions fourni par le systÃ¨me d'amplification et obtenus lors de l'absorption de rayonnements ionisants par le scintillateur L'unitÃ© de traitement 13 est notamment adaptÃ©e pour mettre en oeuvre un procÃ©dÃ© de caractÃ©risation de la source de rayonnement ionisant 20 Ã caractÃ©riser. 5 L'unitÃ© de traitement 13 est sÃ©lectionnÃ©e parmi : - un circuit de traitement dÃ©diÃ© dont la configuration, ceci aussi bien en termes de composants et de programmation, est adaptÃ©e pour la mise en oeuvre du procÃ©dÃ© de caractÃ©risation de la source de rayonnement ionisant 20 Ã caractÃ©riser, - un ordinateur comportant une carte d'acquisition par laquelle 10 l'ordinateur est reliÃ© au photomultiplicateur 12 et un support de donnÃ©es lisible par un ordinateur sur lequel est enregistrÃ© un programme incluant les instructions pour exÃ©cuter des Ã©tapes du procÃ©dÃ© de caractÃ©risation de la source de rayonnement ionisant 20 Ã caractÃ©riser, - une configuration hybride comprenant un circuit de traitement et un 15 ordinateur configurÃ©s pour la mise en oeuvre du procÃ©dÃ© de caractÃ©risation de la source de rayonnement ionisant 20 Ã caractÃ©riser. Un tel procÃ©dÃ© de caractÃ©risation comporte les Ã©tapes suivantes : - fourniture de N signaux temporels d'impulsions de rÃ©fÃ©rence, chaque signal temporel d'impulsion de rÃ©fÃ©rence correspondant Ã la dÃ©tection d'un photon y par 20 un scintillateur, N Ã©tant un entier naturel non nul ; - choix d'au moins une grandeur discriminante ; - pour chaque signal temporel d'impulsion de rÃ©fÃ©rence, extraction de la valeur de chaque grandeur discriminante choisie pour former une valeur de rÃ©fÃ©rence correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de 25 rÃ©fÃ©rence ; - acquisition de M signaux temporels d'impulsions acquis, chaque signal temporel d'impulsion acquis correspondant Ã la dÃ©tection, par le scintillateur, d'un photon issu de la source de rayonnement ionisant Ã caractÃ©riser, M Ã©tant un entier naturel non nul ; 3065814 19 - pour chaque signal temporel d'impulsion acquis, extraction de la valeur de chaque grandeur discriminante choisie pour former une valeur caractÃ©ristique correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion acquis ; 5 - calcul d'une premiÃ¨re valeur caractÃ©ristique de dissemblance Ã partir d'un premier calcul de dissemblance entre une matrice de rÃ©fÃ©rence et une matrice caractÃ©ristique, la matrice de rÃ©fÃ©rence Ã©tant formÃ©e par les valeurs de rÃ©fÃ©rence, la matrice caractÃ©ristique Ã©tant formÃ©e par les valeurs caractÃ©ristiques ; - comparaison de la premiÃ¨re valeur caractÃ©ristique de dissemblance 10 calculÃ©e Ã un premier seuil prÃ©dÃ©terminÃ© ; - caractÃ©risation de la prÃ©sence de neutrons dans un rayonnement Ã©mis par la source de rayonnement ionisant Ã caractÃ©riser en fonction du rÃ©sultat de la comparaison entre la premiÃ¨re valeur caractÃ©ristique de dissemblance calculÃ©e et le premier seuil prÃ©dÃ©terminÃ©. 15 Avec un tel procÃ©dÃ©, il est possible de dÃ©terminer si la source de rayonnement caractÃ©ristique 20 Ã caractÃ©riser prÃ©sente une Ã©mission comprenant des neutrons. L'Ã©tape de fourniture des N signaux temporels d'impulsions de rÃ©fÃ©rence, peut consister en : 20 soit en une caractÃ©risation d'une source de rayonnement ionisant de rÃ©fÃ©rence dont l'Ã©mission de rayonnement ionisant est connue pour ne comporter que des photons gamma y, cette caractÃ©risation Ã©tant prÃ©fÃ©rentiellement rÃ©alisÃ©e avec le systÃ¨me de caractÃ©risation 10, - soit en une caractÃ©risation de plusieurs sources de rayonnement 25 ionisant de rÃ©fÃ©rence dont les Ã©missions de rayonnement ionisant sont connues pour ne comporter que des photons gamma y, ces caractÃ©risations Ã©tant prÃ©fÃ©rentiellement rÃ©alisÃ©es avec le systÃ¨me de caractÃ©risation 10, - soit la rÃ©cupÃ©ration de signaux temporels d'impulsion reprÃ©sentatif de simulÃ©s Ã partir de prÃ©cÃ©dentes caractÃ©risations de source de rayonnement ionisant de 3065814 20 rÃ©fÃ©rence dont les Ã©missions de rayonnement ionisant sont connues pour ne comporter que des photons gamma y. Afin d'illustrer le choix de l'au moins une grandeur discriminante, la figure 2 met en parallÃ¨le ces signaux temporels d'impulsions moyennÃ©s 301, 302 obtenus 5 avec respectivement une source de rayonnement ionisant Ã©mettant uniquement des photons gamma y et une source de rayonnement ionisant Ã©mettant des photons y et des neutrons. On peut ainsi voir que si ces signaux temporels d'impulsions moyennÃ©s 301, 302 sont relativement similaires, ceux-ci se distinguent principalement par deux caractÃ©ristiques qui sont : 10 - le signal temporel moyennÃ© de la source de rayonnement ionisant Ã©mettant uniquement des photons gamma y prÃ©sente une dÃ©croissance plus rapide, et donc une largeur Ã mi-hauteur plus faible que le signal temporel moyennÃ© de la source de rayonnement ionisant Ã©mettant des photons y et des neutrons, - le signal temporel moyennÃ© de la source de rayonnement ionisant 15 Ã©mettant des photons y et des neutrons prÃ©sente un pic secondaire prÃ©pondÃ©rant vis-Ã - vis de celui du signal temporel moyennÃ© de la source de rayonnement ionisant Ã©mettant uniquement des photons gamma y. Ainsi, lors du choix de l'au moins une grandeur discriminante, la grandeur discriminante peut Ãªtre sÃ©lectionnÃ©e parmi : 20 - une intÃ©grale d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis, tel que par exemple l'intÃ©grale de la portion correspondant au pic secondaire prÃ©pondÃ©rant pour la source de rayonnement ionisant Ã©mettant des photons y et des neutrons ; - un temps caractÃ©ristique du signal temporel d'impulsion de rÃ©fÃ©rence 25 ou signal temporel d'impulsion acquis, telle qu'un temps de descente ; - une valeur, en au moins un instant, du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis, tel qu'une amplitude correspondant au pic secondaire ; 3065814 21 - une fonction de la variance d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis, une telle valeur permettant d'illustrer la prÃ©pondÃ©rance ou non du pic secondaire ; - une valeur de dÃ©rivÃ©e du signal temporel d'impulsion de rÃ©fÃ©rence ou 5 signal temporel d'impulsion acquis, une telle valeur pouvant Ãªtre reprÃ©sentative de la dÃ©croissance plus rapide observÃ©e pour le signal temporel moyennÃ© de la source de rayonnement ionisant Ã©mettant uniquement des photons gamma y. Dans une configuration classique de l'invention, chacune des Ã©tapes d'extraction de la valeur de chaque grandeur discriminante, que ce soit pour les signaux 10 temporels d'impulsions de rÃ©fÃ©rence ou pour les signaux temporels d'impulsion acquis, comprend les sous-Ã©tapes suivantes : - normalisation de chaque signal temporel d'impulsion de rÃ©fÃ©rence ou acquis ; rÃ©duction du bruit Ã©lectronique de chaque signal temporel 15 d'impulsion de rÃ©fÃ©rence ou acquis n, - pour chaque signal temporel d'impulsion de rÃ©fÃ©rence ou acquis calcul de la valeur de chaque grandeur discriminante choisie pour former la valeur caractÃ©ristique correspondante. Bien entendu, en fonction de configuration du systÃ¨me de 20 caractÃ©risation 1, chacune des Ã©tapes d'extraction de la valeur de chaque grandeur discriminante peut ne pas comporter l'une ou les deux sous-Ã©tapes parmi la sous-Ã©tape de normalisation et la sous-Ã©tape rÃ©duction du bruit Ã©lectronique. On peut noter que dans le cadre de l'invention, lors de l'Ã©tape de choix de l'au moins une grandeur discriminante, il est prÃ©fÃ©rentiellement choisie au moins deux 25 grandeurs discriminantes, ces grandeurs discriminantes Ã©tant prÃ©fÃ©rentiellement choisies parmi une intÃ©grale d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis et une fonction de la variance d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis. 3065814 22 Lors de cette Ã©tape de choix de l'au moins une grandeur discriminante, il est donc choisi R grandeurs discriminantes, R Ã©tant un entier naturel non nul. Lors de l'Ã©tape de calcul de la premiÃ¨re valeur caractÃ©ristique de dissemblance, la matrice de rÃ©fÃ©rence peut Ãªtre une matrice comportant R colonnes et N 5 lignes, chaque colonne correspondant Ã une grandeur discriminante respective et contenant les valeurs de rÃ©fÃ©rence extraites pour ladite grandeur discriminante. De la mÃªme faÃ§on, lors de de cette mÃªme Ã©tape la matrice caractÃ©ristique peut Ãªtre une matrice comportant R colonnes et M lignes, chaque colonne correspondant Ã une grandeur discriminante respective et contenant les valeurs caractÃ©ristiques extraites pour 10 ladite grandeur discriminante. Bien entendu, il est Ã©galement envisageable, sans que l'on sorte du cadre de l'invention, que la matrice de rÃ©fÃ©rence soit une matrice comportant N colonnes et R lignes ou encore une matrice comportant une unique colonne et NxR lignes, voire NxR colonnes et une unique ligne, tant que ladite matrice est formÃ©e par l'ensemble des 15 valeurs de rÃ©fÃ©rence extraites. Afin de faciliter le calcul de la valeur caractÃ©ristique de dissemblance, la matrice caractÃ©ristique prÃ©sente, bien entendu, un dimensionnement similaire Ã celui de la matrice de rÃ©fÃ©rence. L'Ã©tape de calcul de la valeur caractÃ©ristique de dissemblance comporte les sous-Ã©tapes suivantes : 20 - pour chaque grandeur discriminante identifiÃ©e par un indice i, subdivision d'un intervalle reprÃ©sentatif des valeurs prises par la grandeur discriminante en k, sous-intervalles, k, Ã©tant un entier naturel non nul ; - calcul, pour chacune de la matrice de rÃ©fÃ©rence et la matrice caractÃ©ristique, de la valeur de Mi k, coefficients 0, chaque coefficient 0 Ã©tant associÃ© Ã 25 R sous-intervalles, chacun desdits R sous-intervalles Ã©tant associÃ© Ã une grandeur discriminante respective et Ã©tant Ã©gal Ã un sous-intervalle de ladite grandeur discriminante respective, un premier coefficient 01 et un deuxiÃ¨me coefficient 02 distincts Ã©tant tels qu'il existe au moins une grandeur caractÃ©ristique pour laquelle le sous-intervalle correspondant du premier coefficient 01 est diffÃ©rent du sous-intervalle 30 correspondant du deuxiÃ¨me coefficient 02, la valeur de chacun des coefficients 0 Ã©tant 3065814 23 calculÃ©e comme le nombre de signaux temporels d'impulsions de rÃ©fÃ©rence, respectivement de signaux temporels d'impulsions acquis, dont les valeurs de grandeurs discriminantes extraites sont dans les P sous-intervalles associÃ©s audit coefficient, les valeur de coefficients 0 calculÃ©es pour l'une de la matrice de rÃ©fÃ©rence et de la matrice 5 caractÃ©ristique Ã©tant multipliÃ©es par un facteur correctif correspondant au ration du nombre de signaux temporels d'impulsions entre la matrice de rÃ©fÃ©rence et la matrice caractÃ©ristique ; - calcul de la premiÃ¨re valeur caractÃ©ristique de dissemblance Ã partir d'une diffÃ©rence en valeur absolue mettant en oeuvre les coefficients 0 calculÃ©s Ã partir 10 de la matrice de rÃ©fÃ©rence et de la matrice caractÃ©ristique, et de R sommes successives, chacune des R sommes Ã©tant associÃ©e Ã une grandeur discriminante respective, et ladite diffÃ©rence se faisant directement entre les coefficients 0 calculÃ©s sur la base de la matrice de rÃ©fÃ©rence et les coefficients 0 calculÃ©s sur la base de la matrice caractÃ©ristique correspondant au mÃªme R sous-intervalles, la premiÃ¨re somme se faisant sur la base du 15 rÃ©sultat de ladite diffÃ©rence et les R-1 sommes suivantes sur la base du rÃ©sultat de la somme la prÃ©cÃ©dent directement . On notera bien entendu que par soucis de concision, la terminologie Â« une diffÃ©rence entre chaque coefficient 0 calculÃ© sur la base de la premiÃ¨re matrice et le coefficient 0 calculÃ© sur la base de la deuxiÃ¨me matrice correspondant au mÃªme R 20 sous-intervalles Â», s'entend Trh ki diffÃ©rences, chacune correspondant Ã un coefficient 0 respectif calculÃ© sur la base de la premiÃ¨re matrice. De la mÃªme faÃ§on, une somme donnÃ©e parmi les R sommes, mis Ã part, la en' somme, correspond Ã une pluralitÃ© de sommes. Ainsi la TIÃ¨me somme, T Ã©tant compris entre 1 et R, R exclut, correspond Ã T-1 k, sommes, ki correspondant aux sous- 25 intervalles des grandeurs pour la ou lesquelles la somme correspondante n'a pas encore Ã©tÃ© calculÃ©e. De cette maniÃ¨re, la valeur caractÃ©ristique de dissemblance correspond Ã une somme des diffÃ©rences en valeur absolue entre la rÃ©partition des valeurs de grandeurs discriminantes obtenues pour les signaux temporels d'impulsions de rÃ©fÃ©rence 30 et les signaux temporels d'impulsions caractÃ©ristiques. Il est Ã©galement possible de 3065814 24 mettre en exergue les diffÃ©rences de rÃ©partition des signaux temporels d'impulsion liÃ©es Ã la prÃ©sence d'Ã©ventuels neutrons dans le rayonnement ionisant 21 Ã©mis par la source de rayonnement ionisant 20 Ã caractÃ©riser. On peut noter qu'il est Ã©galement envisageable dans le cadre de 5 l'invention, comme cela sera illustrÃ© en lien avec les figures 6A et 6B, d'effectuer une ou plusieurs sommes prÃ©alable Ã la diffÃ©rence, cette ou ces sommes Ã©tant chacune associÃ©e Ã une grandeur discriminante respective. En effet, de telles sommes effectuÃ©es avant la diffÃ©rence est parfaitement Ã©quivalent Ã avoir choisi un nombre de grandeurs discriminantes rÃ©duites du nombre de sommes rÃ©alisÃ©es prÃ©alablement Ã ladite 10 diffÃ©rence. Ainsi, si une telle possibilitÃ© n'est clairement pas prÃ©fÃ©rable, elle est parfaitement compatible avec l'invention, puisqu'elle induit juste, en plus des Ã©tapes revendiquÃ©es dans l'invention, des calculs supplÃ©mentaires rÃ©duisant le nombre de grandeurs discriminantes rÃ©ellement mise en oeuvre. Afin d'optimiser le choix de l'au moins une grandeur discriminante, 15 celle-ci Ã©tant choisie aprÃ¨s la fourniture de N signaux temporels d'impulsions de rÃ©fÃ©rence, l'Ã©tape de choix de ladite au moins une grandeur discriminante peut comporter les Ã©tapes suivantes : - fourniture de P signaux temporels d'impulsions de comparaison, chaque signal temporel d'impulsion de comparaison correspondant Ã la dÃ©tection, par le 20 scintillateur, d'un photon y et/ou d'un neutron Ã©mis par une source Ã©mettant un rayonnement ionisant comportant des photons y et des neutrons, P Ã©tant un entier naturel non nul ; - construction d'une pluralitÃ© d'ensembles de grandeurs discriminantes ; - pour chaque signal temporel d'impulsion de rÃ©fÃ©rence, et pour chaque 25 ensemble de grandeurs discriminantes, extraction de la valeur de chaque grandeur discriminante dudit ensemble de grandeurs discriminantes pour former une premiÃ¨re valeur de comparaison correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de rÃ©fÃ©rence; - pour chaque signal temporel d'impulsion de comparaison, et pour 30 chaque ensemble de grandeurs discriminantes, extraction de la valeur de chaque 3065814 25 grandeur discriminante dudit ensemble de grandeurs discriminantes pour former une deuxiÃ¨me valeur de comparaison correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de comparaison ; - pour chaque ensemble de grandeurs discriminantes, calcul d'une 5 deuxiÃ¨me valeur caractÃ©ristique de dissemblance Ã partir d'un deuxiÃ¨me calcul de dissemblance entre une premiÃ¨re matrice constituÃ©e des premiÃ¨res valeurs de comparaison extraites pour ledit ensemble de grandeurs discriminantes et une deuxiÃ¨me matrice constituÃ©e des deuxiÃ¨mes valeurs de comparaison extraites pour ledit ensemble de grandeurs discriminantes ; 10 - choix d'au moins une grandeur discriminante correspondant Ã l'un des ensembles de grandeurs discriminantes, le choix Ã©tant rÃ©alisÃ© sur la base des deuxiÃ¨mes valeurs caractÃ©ristique de dissemblance calculÃ©e pour chacun des ensembles de grandeurs discriminantes. Dans une configuration classique de l'invention, chacune des Ã©tapes 15 d'extraction de la valeur de chaque grandeur discriminante, que ce soit pour les signaux temporels d'impulsions de rÃ©fÃ©rence ou pour les signaux temporels d'impulsion de comparaison, comprend les sous-Ã©tapes suivantes : - normalisation de chaque signal temporel d'impulsion de rÃ©fÃ©rence ou de comparaison ; 20 rÃ©duction du bruit Ã©lectronique de chaque signal temporel d'impulsion de rÃ©fÃ©rence ou de comparaison, - pour chaque signal temporel d'impulsion de rÃ©fÃ©rence ou de comparaison calcul de la valeur de chaque grandeur discriminante choisie pour former la valeur caractÃ©ristique correspondante. 25 Bien entendu, en fonction de configuration du systÃ¨me de caractÃ©risation 1, chacune des Ã©tapes d'extraction de la valeur de chaque grandeur discriminante peut ne pas comporter l'une ou les deux sous-Ã©tapes parmi la sous-Ã©tape de normalisation et la sous-Ã©tape rÃ©duction du bruit Ã©lectronique. La sous-Ã©tape de construction d'une pluralitÃ© d'ensembles de 30 grandeurs discriminantes, peut consister en une sÃ©lection de plusieurs ensembles de 3065814 26 grandeurs discriminantes considÃ©rÃ©es comme reprÃ©sentative de la diffÃ©rence entre le signal temporel d'impulsion d'un photon gamma y et le signal temporel d'impulsion d'un neutron aprÃ¨s que ces deux signaux aient Ã©tÃ©. Ainsi lors de la construction de chacun des ensembles de grandeurs discriminantes, chacune des grandeurs discriminantes peut Ãªtre 5 sÃ©lectionnÃ©e parmi : - une intÃ©grale d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis, tel que par exemple l'intÃ©grale de la portion correspondant au pic secondaire prÃ©pondÃ©rant pour la source de rayonnement ionisant Ã©mettant des photons y et des neutrons ; 10 - un temps caractÃ©ristique du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis, telle qu'un temps de descente ; - une valeur, en au moins un instant, du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis, tel qu'une amplitude correspondant au pic secondaire ; 15 - une fonction de la variance d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis, une telle valeur permettant d'illustrer la prÃ©pondÃ©rance ou non du pic secondaire ; - une valeur de dÃ©rivÃ©e du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis, une telle valeur pouvant Ãªtre reprÃ©sentative de la 20 dÃ©croissance plus rapide observÃ©e pour le signal temporel moyennÃ© de la source de rayonnement ionisant Ã©mettant uniquement des photons gamma y. On notera que cette mÃªme sous Ã©tapes de construction d'une pluralitÃ© d'ensemble de grandeurs discriminantes peut consister en un choix d'au moins une fonction de grandeur discriminante qui est fonction d'une variable ou plusieurs variables, 25 chacun des ensembles de grandeurs discriminantes correspondant Ã une valeur de la variable ou des variables donnÃ©es et comportant une grandeur discriminante rÃ©sultant de la fonction de grandeur discriminante pour ladite ou lesdites valeurs de la ou des variables donnÃ©es. Ainsi, une telle fonction de grandeur de discrimination peut Ãªtre, par 30 exemple, une intÃ©grale d'une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou 3065814 27 signal temporel d'impulsion acquis, une ou les limites de ladite portion Ã©tant dÃ©finies par la ou les variables donnÃ©es, ou encore une valeur de dÃ©rivÃ©e du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis Ã une position dudit signal temporel d'impulsion dÃ©finie Ã partir de ladite ou lesdites variables donnÃ©es. 5 Selon une possibilitÃ© avantageuse de l'invention, au moins une fonction de grandeur discriminante de chacun des ensembles de grandeurs discriminantes sont sÃ©lectionnÃ©s dans le groupe comportant : une grandeur discriminante notÃ©e (-lot correspondant Ã une intÃ©grale d'une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel 10 d'impulsion acquis, une grandeur discriminante notÃ©e (lad, correspondant Ã une intÃ©grale d'une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis, une grandeur discriminante est notÃ©e D correspondant Ã une 15 fonction de la variance d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis. La grandeur discriminante (lot est une grandeur discriminante calculÃ©e selon : Qt0t= ft.s(t) dt oÃ¹ s est le signal temporel d'impulsion de rÃ©fÃ©rence, respectivement le 20 signal temporel d'impulsion acquis, tmax est la coordonnÃ©e temporelle du maximum du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, et t1 est une coordonnÃ©e temporelle calculÃ©e selon : t1= tmax+13.tHm oÃ¹ p est un rÃ©el strictement positif, et tHM est la coordonnÃ©e temporelle, ultÃ©rieure Ã tmax, pour laquelle la valeur du signal temporel d'impulsion de rÃ©fÃ©rence, 25 respectivement du signal temporel d'impulsion acquis, vaut 0,5. La grandeur discriminante Qtad est une grandeur discriminante calculÃ©e selon : 3065814 28 Qtall= f dt to oÃ¹ s est le signal temporel d'impulsion de rÃ©fÃ©rence, respectivement le signal temporel d'impulsion acquis, to est une coordonnÃ©e temporelle calculÃ©e selon : tmax+a.tHm oÃ¹ a est un rÃ©el strictement positif, tmax est la coordonnÃ©e temporelle du maximum du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal 5 temporel d'impulsion acquis, et tHM est la coordonnÃ©e temporelle, ultÃ©rieure Ã tmax, pour laquelle la valeur du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, vaut 0, et t1 est une coordonnÃ©e temporelle calculÃ©e selon : t1= tmax+13.tHm oÃ¹ p est un rÃ©el strictement positif. La grandeur discriminante D est la variance de l'intÃ©gral de la portion du 10 signal temporel d'impulsion de rÃ©fÃ©rence, respectivement le signal temporel d'impulsion acquis, comprise entre to et t1, oÃ¹ to est une coordonnÃ©e temporelle calculÃ©e selon : tmax+a.tHm oÃ¹ a est un rÃ©el strictement positif, tmax est la coordonnÃ©e temporelle du maximum du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, et tHM est la coordonnÃ©e temporelle, ultÃ©rieure Ã tmax, pour 15 laquelle la valeur du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, vaut 0,5, t1 est une coordonnÃ©e temporelle calculÃ©e selon : t1= tmax+13.tHm oÃ¹ p est un rÃ©el strictement positif. Autrement dit, D est la variance de la granduer discriminante Qtad du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement le signal temporel d'impulsion 20 acquis. La figure 3 illustre le calcul des valeurs Qtot, Qtail et D en figurant sur un exemple de signal temporel d'impulsion, outre les portions de signal temporel d'impulsion correspondant Ã Qtot et Qtail, les temps correspondant Ã to, t1, tmax -max et THM. Les paramÃ¨tres a et p sont des valeurs rÃ©elles qui sont utilisÃ©es pour prendre en compte 3065814 29 le temps de rÃ©ponse du scintillateur 11 et du photomultiplicateur 12. S'il est possible de les dÃ©terminer Ã partir des donnÃ©es caractÃ©ristiques de l'ensemble scintillateur 11 et photomultiplicateur 12, il est nÃ©anmoins Ã©galement possible, comme cela est illustrÃ© Ã la suite de ce document, de les dÃ©terminer par la mise en oeuvre d'une Ã©tape de choix de 5 l'au moins une grandeur discriminante. Ainsi, avec de telles fonctions de grandeurs discriminantes, il est possible de former une pluralitÃ© d'ensembles de grandeurs discriminantes basÃ©s se distinguant deux Ã deux par au moins l'un des paramÃ¨tres a et R. Lors de la sous Ã©tape de choix d'au moins une grandeur discriminante 10 correspondant Ã l'un des ensembles de grandeurs discriminantes, l'au moins une grandeur discriminante correspondant choisie peut Ãªtre celle qui correspond Ã l'ensemble de grandeurs discriminantes pour lequel la deuxiÃ¨me valeur caractÃ©ristique de dissemblance prÃ©sente un extrÃ©mum caractÃ©ristique d'une dissemblance maximale. Ainsi, dans le cas oÃ¹ la valeur caractÃ©ristique de dissemblance est une 15 valeur caractÃ©ristique de dissemblance en tant que telle qui permet de caractÃ©riser les diffÃ©rences entre les signaux d temporels d'impulsions de rÃ©fÃ©rence et les signaux temporels d'impulsions de comparaison, l'extremum caractÃ©ristique d'une dissemblance est un maximum. Dans le cas oÃ¹ la valeur caractÃ©ristique de dissemblance est une valeur de vraisemblance qui permet de caractÃ©riser les similaritÃ©s entre les signaux temporels 20 d'impulsions de rÃ©fÃ©rence et les signaux temporels d'impulsions de comparaison l'extremum caractÃ©ristique d'une dissemblance maximale est un minimum. En variante Ã une telle sous Ã©tapes de choix d'au moins une grandeur discriminante correspondant Ã l'un des ensembles de grandeurs discriminantes et ceci afin de permettre de limiter les risque, lors de la caractÃ©risation de la source de 25 rayonnement ionisant 20 Ã caractÃ©riser, l'Ã©tape de choix d'au moins une grandeur discriminante comporte en outre les sous-Ã©tapes suivantes : - fourniture de Q signaux temporels d'impulsions de test, chaque signal temporel d'impulsion de test correspondant Ã la dÃ©tection, par le scintillateur, d'un photon y Ã©mis par une source Ã©mettant un rayonnement ionisant comportant des 30 photons y, Q Ã©tant un entier naturel non nul, 3065814 30 - pour chaque signal temporel d'impulsion de test, et pour chaque ensemble de grandeurs discriminantes, extraction de la valeur de chaque grandeur discriminante dudit ensemble de grandeurs discriminantes pour former une troisiÃ¨me valeur de comparaison correspondant Ã ladite grandeur discriminante et audit signal 5 temporel d'impulsion de test ; - pour chaque ensemble de grandeurs discriminantes, calcul d'une troisiÃ¨me valeur caractÃ©ristique de dissemblance Ã partir du deuxiÃ¨me calcul de dissemblance entre une premiÃ¨re matrice constituÃ©e des premiÃ¨res valeurs de comparaison extraites pour ledit ensemble de grandeurs discriminantes et une deuxiÃ¨me 10 matrice constituÃ©e des troisiÃ¨me valeurs de comparaison extraites pour ledit ensemble de grandeurs discriminantes ; - comparaison de la troisiÃ¨me valeur caractÃ©ristique de dissemblance calculÃ©e pour chacun des ensembles de grandeurs discriminantes Ã un deuxiÃ¨me seuil prÃ©dÃ©terminÃ© ; 15 - sÃ©lection parmi les ensembles de grandeurs discriminantes ceux qui prÃ©sentent une dissemblance minimale Ã partir de ladite comparaison entre la troisiÃ¨me valeur caractÃ©ristique de dissemblance calculÃ©e pour chacun des ensembles de grandeurs discriminantes Ã un deuxiÃ¨me seuil prÃ©dÃ©terminÃ© ; - choix d'au moins une grandeur discriminante correspondant Ã l'un des 20 ensembles de grandeurs discriminantes, l'au moins une grandeur discriminante correspondant choisie est celle correspondant Ã l'ensemble de grandeurs discriminantes parmi les ensembles de grandeur discriminante sÃ©lectionnÃ© dont la deuxiÃ¨me valeur caractÃ©ristique de dissemblance prÃ©sente un extrÃ©mum caractÃ©ristique d'une dissemblance maximale. 25 Comme indiquÃ© prÃ©cÃ©demment, l'Ã©tape d'extraction de la valeur de chaque grandeur discriminante comprend les sous-Ã©tapes suivantes : normalisation de chaque signal temporel d'impulsion de test ; rÃ©duction du bruit Ã©lectronique de chaque signal temporel d'impulsion de test, 3065814 31 - pour chaque signal temporel d'impulsion de test calcul de la valeur de chaque grandeur discriminante choisie pour former la valeur caractÃ©ristique correspondante. Bien entendu, en fonction de configuration du systÃ¨me de 5 caractÃ©risation 1, chacune des Ã©tapes d'extraction de la valeur de chaque grandeur discriminante peut ne pas comporter l'une ou les deux sous-Ã©tapes parmi la sous-Ã©tape de normalisation et la sous-Ã©tape rÃ©duction du bruit Ã©lectronique. Lors de cette Ã©tape de choix de l'au moins une grandeur discriminante, un ensemble donnÃ© de grandeurs discriminantes comporte S grandeurs discriminantes, S 10 Ã©tant un entier naturel non nul. La deuxiÃ¨me et la troisiÃ¨me valeur caractÃ©ristique de dissemblance peuvent Ãªtre calculÃ©es Ã partir d'un calcul de dissemblance similaire Ã celui de la premiÃ¨re valeur caractÃ©ristique de dissemblance. Ainsi les calculs de dissemblance de la deuxiÃ¨me et la troisiÃ¨me valeur caractÃ©ristique de dissemblance peuvent comporter les sous-Ã©tapes suivantes : 15 - pour chaque grandeur discriminante identifiÃ©e par un indice j, subdivision d'un intervalle reprÃ©sentatif des valeurs prises par la grandeur discriminante en k, sous-intervalles, k, Ã©tant un entier naturel non nul ; - calcul, pour chacune de la premiÃ¨re matrice et de la deuxiÃ¨me matrice, de la valeur de Tls1kj coefficients 6, chaque coefficient 6 Ã©tant associÃ© Ã S sous- 20 intervalles, chacun desdits S sous-intervalles Ã©tant associÃ© Ã une grandeur discriminante respective et Ã©tant Ã©gal Ã un sous-intervalle de ladite grandeur discriminante respective, un premier coefficient 61 et un deuxiÃ¨me coefficient 62 distincts Ã©tant tels qu'il existe au moins une grandeur caractÃ©ristique pour laquelle le sous-intervalle correspondant du premier coefficient 61 est diffÃ©rent du sous-intervalle correspondant du deuxiÃ¨me 25 coefficient 62, la valeur de chacun des coefficients 6 Ã©tant calculÃ©e comme le nombre de signaux temporels d'impulsions de la premiÃ¨re matrice, respectivement de signaux temporels d'impulsions de la deuxiÃ¨me matrice, dont les valeurs de grandeurs discriminantes extraites sont dans les S sous-intervalles associÃ©s audit coefficient, les valeurs de coefficients 6 calculÃ©s pour l'une de la premiÃ¨re matrice et de la deuxiÃ¨me 3065814 32 matrice Ã©tant multipliÃ©es par un facteur correctif correspondant au ration du nombre signaux temporels d'impulsions entre la premiÃ¨re matrice et la deuxiÃ¨me matrice ; - calcul d'une valeur caractÃ©ristique de dissemblance Ã partir d'une diffÃ©rence en valeur absolue mettant en oeuvre les coefficients 6 calculÃ©s Ã partir de la 5 premiÃ¨re matrice et de la deuxiÃ¨me matrice et de S sommes successives, chacune des S sommes Ã©tant associÃ©e Ã une grandeur discriminante respective, et ladite diffÃ©rence se faisant directement entre chaque coefficient 6 calculÃ© sur la base de la premiÃ¨re matrice et le coefficient 6 calculÃ© sur la base de la deuxiÃ¨me matrice correspondant au mÃªme R sous-intervalles, Ila premiÃ¨re somme se faisant sur la base du rÃ©sultat de ladite 10 diffÃ©rence et les S-1 sommes suivantes sur la base du rÃ©sultat de la somme la prÃ©cÃ©dent directement. On notera bien entendu que par soucis de concision, la terminologie Â« une diffÃ©rence entre chaque coefficient 6 calculÃ© sur la base de la premiÃ¨re matrice et le coefficient 6 calculÃ© sur la base de la deuxiÃ¨me matrice correspondant au mÃªme S 15 sous-intervalles Â», s'entend re_i kj diffÃ©rences, chacune correspondant Ã un coefficient 6 respectif calculÃ© sur la base de la premiÃ¨re matrice. De la mÃªme faÃ§on, une somme donnÃ©e parmi les S sommes, mis Ã part, la SIÃ¨me somme, correspond Ã une pluralitÃ© de sommes. Ainsi la UlÃ¨me somme, U Ã©tant compris entre 1 et S, S exclut, correspond Ã ri,_1k, sommes, kj correspondant aux sous- 20 intervalles des grandeurs discriminante pour la ou lesquelles la somme correspondante n'a pas encore Ã©tÃ© calculÃ©e. Ainsi, avec un tel calcul, les calculs de la deuxiÃ¨me et de la troisiÃ¨me valeur caractÃ©ristique de dissemblance se distinguent uniquement en ce que la deuxiÃ¨me matrice est constituÃ©e des deuxiÃ¨mes valeurs de comparaison calculÃ©es pour ledit 25 ensemble de grandeurs discriminantes pour le calcul de la deuxiÃ¨me valeur caractÃ©ristique de dissemblance alors que la deuxiÃ¨me matrice est constituÃ©e des troisiÃ¨me valeurs de comparaison calculÃ©es pour ledit ensemble de grandeurs discriminantes pour le calcul de la deuxiÃ¨me valeur caractÃ©ristique de dissemblance. On peut noter que, selon une possibilitÃ© avantageuse de l'invention, le 30 premier seuil prÃ©dÃ©terminÃ© peut Ãªtre dÃ©terminÃ© Ã partir du rÃ©sultat de l'Ã©tape de choix 3065814 33 de l'au moins une grandeur discriminante. Ainsi dans le cas oÃ¹ le choix de l'au moins une grandeur discriminante fait intervenir les calculs de la deuxiÃ¨me et la troisiÃ¨me valeur caractÃ©ristique de dissemblance, le premier seuil prÃ©dÃ©terminÃ© peut Ãªtre sÃ©lectionnÃ© dans un intervalle dÃ©limitÃ© par les deuxiÃ¨me et troisiÃ¨me valeurs caractÃ©ristique de 5 dissemblance pour l'ensemble de grandeurs discriminantes sÃ©lectionnÃ©es. Selon une autre possibilitÃ© de l'invention, le procÃ©dÃ© de caractÃ©risation peut Ã©galement comporter une Ã©tape supplÃ©mentaire d'estimation du nombre et/ou d'un taux de signaux temporels d'impulsions correspondant Ã un neutron. Une telle estimation peut notamment Ãªtre calculÃ©e Ã partir de la premiÃ¨re valeur caractÃ©ristique de 10 dissemblance puisque cette valeur de caractÃ©ristique est, lorsque cette derniÃ¨re est une valeur caractÃ©ristique de dissemblance en tant que tel, directement proportionnelle aux taux de signaux temporels d'impulsions correspondant Ã un neutron. Les figures 4A Ã 6B illustrent un premier exemple de mise en oeuvre d'un tel procÃ©dÃ© dans lequel le choix de l'au moins une grandeur discriminante est rÃ©alisÃ©e Ã 15 partir d'un ensemble de grandeurs discriminantes comportant (lot et Qtad et pour lequel le choix de l'au moins une grandeur discriminante fait appel Ã un calcul des troisiÃ¨me valeurs caractÃ©ristiques de dissemblance. Dans le cadre de ce premier exemple de mise en oeuvre, quatre sources de rayonnement ionisant ont sÃ©lectionnÃ©, trois sources de rayonnement ionisant dont 20 l'Ã©mission de rayonnement ionisant est constituÃ©e uniquement de photons gamma y et une source dont l'Ã©mission de rayonnement ionisant est constituÃ©e de photons gamma y et de neutrons. Les trois sources de rayonnement ionisant dont l'Ã©mission de rayonnement ionisant est constituÃ©e uniquement de photons gamma sont des sources constituÃ©es respectivement de l'isotope 22 du sodium 22Na, de l'isotope 137 du cÃ©sium 25 137Cs et de l'isotope 60 du cobalt 60Co. La source dont l'Ã©mission de rayonnement ionisant est constituÃ©e de photons gamma y et de neutrons est une source constituÃ©e d'isotopes 252 du californium 252Cf. L'ensemble de ces sources de rayonnement ionisant a Ã©tÃ© caractÃ©risÃ© Ã partir d'un scintillateur liquide standard, le scintillateur BC-501-AÂ® commercialisÃ© par 30 Saint-GobainÂ® afin de dÃ©finir les N signaux temporels d'impulsions de rÃ©fÃ©rence et les M 3065814 34 signaux temporels d'impulsion acquis. Ainsi, avec N Ã©gal Ã 16000, les signaux temporels d'impulsions de rÃ©fÃ©rence sont constituÃ©s de signaux temporels de rÃ©fÃ©rence de chacune des trois sources de rayonnement ionisant dont l'Ã©mission de rayonnement ionisant est constituÃ©e uniquement de photons gamma. De la mÃªme faÃ§on, M est Ã©galement Ã©gal Ã 5 16000 et les signaux temporels d'impulsions acquis sont constituÃ©s de 16000 signaux temporels acquis de la source de l'isotope 252 du californium 252CF. ConformÃ©ment au procÃ©dÃ© de caractÃ©risation selon l'invention, chacun des signaux temporels d'impulsions de rÃ©fÃ©rence et acquis a Ã©tÃ© ceci en les divisant par la valeur maximale de de leur amplitude. La figure 2 illustre ainsi les moyennes respectives 10 des signaux temporels d'impulsion de rÃ©fÃ©rence et des signaux temporels d'impulsion acquis. Une fois les signaux temporels d'impulsions, les Ã©tapes de calcul des valeurs de rÃ©fÃ©rence et des valeurs caractÃ©ristiques ont Ã©tÃ© mis en oeuvre. Pour se faire, les grandeurs discriminantes Qtot et Qtail ont Ã©tÃ© calculÃ©s pour chacun des signaux 15 temporels d'impulsion de rÃ©fÃ©rence et acquis. Pour se faire, et selon le rÃ©sultat du choix de l'au moins une grandeur discriminante prÃ©sentÃ© Ã la suite, les paramÃ¨tres a et p ont Ã©tÃ© fixÃ©s Ã respectivement 4 et 20. Les 32000 valeurs de rÃ©fÃ©rence et de ces 32000 valeurs caractÃ©ristiques extraites ont ensuite permis de constituer la matrice de rÃ©fÃ©rence et la matrice 20 caractÃ©ristique. La matrice de rÃ©fÃ©rence est ainsi une matrice comportant deux colonnes, l'une correspondant Ã (lot et l'autre Ã Qtail et 16000 lignes correspondant aux signaux temporels d'impulsions. Afin de permettre un calcul de la valeur caractÃ©ristique de dissemblance, l'intervalle reprÃ©sentatif des valeurs prises par Qtot a Ã©tÃ© subdivisÃ© en 300 25 sous-intervalles tandis que l'intervalle reprÃ©sentatif des valeurs prises par Qtail a Ã©tÃ© subdivisÃ© en 150 sous-intervalles. Pour chacune de la matrice de rÃ©fÃ©rence et de la matrice caractÃ©ristique, 45000 coefficients 0 correspondant chacun Ã un sous-intervalle donnÃ© de (lot et un sous-intervalle donnÃ© de Qtail ont Ã©tÃ© calculÃ©s. Chacun des coefficients 0 est Ã©gal 30 au nombre de signaux temporels d'impulsions de ladite matrice dont les valeurs de 3065814 rÃ©fÃ©rence, respectivement les valeurs caractÃ©ristiques, dont les valeurs de grandeurs discriminantes sont dans les sous-intervalles de (lot et de Qtad correspondant. Les 45000 coefficients 0 calculÃ©s Ã partir de la matrice de rÃ©fÃ©rence permettent ainsi d'obtenir une matrice de 300 par 150. De maniÃ¨re identique, les 45000 coefficients 0 calculÃ©s Ã partir 5 de la matrice caractÃ©ristique permettent ainsi d'obtenir une matrice de 300 par 150. La figure 4A illustre graphiquement les valeurs de coefficients 0 obtenues pour la matrice de rÃ©fÃ©rence, avec en abscisse les valeurs des sous-intervalles de Qtot et en ordonnÃ©e les valeurs des sous-intervalles de Qta,i, les valeurs des coefficients 0 Ã©tant figurÃ©es au moyen d'une Ã©chelle de couleur. On peut ainsi voir sur cette figure 10 que dans le cas des valeurs de rÃ©fÃ©rence, les grandeurs discriminantes prÃ©sentent une variance relativement faible puisqu'elles sont rÃ©parties en un ensemble 321 relativement dense. Ainsi, les valeurs de la grandeur discriminante (lot sont comprises entre 60 et 90 et celles de la grandeur discriminante Qtad sont comprises entre 0 et 13. Concernant les valeurs de coefficients 0 obtenues pour la matrice de 15 caractÃ©ristique, qui sont figurÃ©es sur la figure 4B selon la mÃªme reprÃ©sentation que la figure 4A, on note une variance beaucoup plus forte des grandeurs discriminantes. En effet, on remarque que ces grandeurs discriminantes sont rÃ©parties en un premier ensemble 322 relativement dense et un deuxiÃ¨me ensemble 323 plus Ã©tendu. Ainsi, les valeurs de la grandeur discriminante (lot sont comprises entre 55 et 140 et celles de la 20 grandeur discriminante Qtad sont comprises entre 2 et 47. Une fois les coefficients 0 calculÃ©s Ã la fois pour la matrice de rÃ©fÃ©rence et pour la matrice caractÃ©ristique, le calcul de la premiÃ¨re valeur caractÃ©ristique de dissemblance est rÃ©alisÃ© Ã partir de deux sommes et d'une diffÃ©rence. Le calcul de la premiÃ¨re valeur caractÃ©ristique de dissemblance peut 25 consister Ã rÃ©aliser pour chacun des intervalles des grandeurs (lot et Qta,i, une premiÃ¨re diffÃ©rence en valeur absolue entre les coefficients 0 de la matrice de rÃ©fÃ©rence et de la matrice caractÃ©ristique, suivit de deux sommes successives du rÃ©sultat de cette diffÃ©rence sur respectivement l'une des grandeurs discriminantes parmi (lot et Qtad et l'autre des grandeurs discriminantes parmi Cltot et Qtail. 3065814 36 Afin d'illustrÃ© le rÃ©sultat d'un tel calcul de la premiÃ¨re valeur caractÃ©ristique de dissemblance, un calcul intermÃ©diaire a Ã©tÃ© rÃ©alisÃ© pour illustrer les dissemblances observable Ã partir de la grandeur discriminante (lot. Ce calcul consiste Ã faire successivement : 5 - une premiÃ¨re somme des coefficients 0 de la matrice de rÃ©fÃ©rence, respectivement de la matrice caractÃ©risant, en sommant, pour chacun des sous-intervalles de la grandeur discriminante Otot, le nombre de signaux temporels d'impulsions de chaque sous intervalles de la grandeur discriminante Qtad correspondant audit intervalle, le rÃ©sultat de chacune des premiÃ¨res sommes pour la matrice de 10 rÃ©fÃ©rence 331 et pour la matrice caractÃ©ristique 332 Ã©tant figurÃ© sur la figure 5A- pour chacun des sous-intervalles de la grandeur discriminante Otot, une diffÃ©rence en valeur absolue entre le rÃ©sultat de la premiÃ¨re somme obtenue pour la matrice de rÃ©fÃ©rence et le rÃ©sultat de la premiÃ¨re somme obtenue pour la matrice caractÃ©ristique, le rÃ©sultat de cette diffÃ©rence Ã©tant figurÃ© sous la forme de portions hachurÃ©es 333 sur la figure 5A. 15 De la mÃªme faÃ§on, un calcul intermÃ©diaire similaire a Ã©tÃ© rÃ©alisÃ© pour illustrer les dissemblances observables Ã partir de la grandeur discriminante Qtad. Ce calcul consiste Ã faire successivement : - une premiÃ¨re somme des coefficients 0 de la matrice de rÃ©fÃ©rence, respectivement de la matrice caractÃ©risant, en sommant, pour chacun des sous- 20 intervalles de la grandeur discriminante Qta,i, le nombre de signaux temporels d'impulsions de chaque sous intervalles de la grandeur discriminante Qtot correspondant audit intervalle, le rÃ©sultat de chacune des premiÃ¨res sommes pour la matrice de rÃ©fÃ©rence 341 et pour la matrice caractÃ©ristique 342 Ã©tant figurÃ© sur la figure 5B, - pour chacun des sous-intervalles de la grandeur discriminante Qta,i, une 25 diffÃ©rence en valeur absolue entre le rÃ©sultat de la premiÃ¨re somme obtenue pour la matrice de rÃ©fÃ©rence et le rÃ©sultat de la premiÃ¨re somme obtenue pour la matrice caractÃ©ristique, le rÃ©sultat de cette diffÃ©rence Ã©tant figurÃ© sous la forme de portions hachurÃ©es 343 sur la figure 5B. Ainsi avec un tel calcul d'une valeur caractÃ©ristique de dissemblance, il 30 est possible de sommer les diffÃ©rences observables sur la figure 5A et figurÃ©es par les 3065814 37 portions hachurÃ©es 333 avec celles observables sur la figure 5A et figurÃ©es par les portions hachurÃ©es 343. A partir du calcul de la valeur caractÃ©ristique de dissemblance, il est possible d'estimer le taux de signaux temporels d'impulsions produits par l'interaction de 5 neutron. En effet, le taux de signaux temporels d'impulsion produits par l'interaction de neutrons est directement proportionnel Ã la valeur caractÃ©ristique de dissemblance calculÃ©e. Dans la configuration utilisÃ©e pour cette mise en oeuvre, c'est-Ã -dire avec une source de rayonnement ionisant dont le rayonnement ionisant comporte des photons gamma y et des neutrons qui est une source de rayonnement ionisant comprenant 10 l'isotope 252 du Californium 252Cf, il a Ã©tÃ© estimÃ© que 29,9% des 16000 signaux temporels d'impulsion Ã©taient des signaux temporels d'impulsion produits par l'interaction de neutrons. Un tel procÃ©dÃ© ne se limite, bien entendu, pas aux seuls scintillateurs organiques liquides et est applicable Ã tout scintillateur organique quel que soit sont type. 15 Afin d'illustrer cette possibilitÃ©, le procÃ©dÃ© a Ã©tÃ© mis en oeuvre selon ce mÃªme premier exemple de mise en oeuvre sur trois autre scintillateurs organiques non liquides, un scintillateur E32Â® conÃ§u et rÃ©alisÃ© par le CEA et des scintillateurs EJ-299-33Â® et EJ-200Â® commercialisÃ©s tous deux par Ejen TechnologyÂ®. Les rÃ©sultats sont synthÃ©tisÃ©s dans le tableau 1 suivant dans lequel est indiquÃ© pour chacun des scintillateurs l'estimation du 20 taux de signaux temporels d'impulsions produits par l'interaction de neutron obtenue pour la source de rayonnement ionisant comprenant l'isotope 252 du californium 252Cf en parallÃ¨le avec un taux de signaux temporels d'impulsion produits par l'interaction de neutrons pour une source de rayonnement ionisant dont le rayonnement ionisant ne comporte que des photons gamma y, cette source ayant Ã©tÃ© simulÃ©e Ã partir d'un sous 25 Ã©chantillonnage alÃ©atoire Ã 8000 signaux temporels d'impulsion des 16000 signaux temporels d'impulsions de rÃ©fÃ©rence. Scintillateur Estimation du taux de signaux temporels d'impulsion produits par l'interaction de neutrons 252Cf Source y 3065814 38 BC-501-AÂ® 29,9 4 E32Â® 26,0 4 EJ-299-33Â® 26,3 4 EJ-200Â® 23,2 4 Tableau 1: Estimation du taux de signaux temporels d'impulsion de neutrons pour respectivement une source de rayonnement ionisant comprenant l'isotope 252 du Californium et pour un sous Ã©chantillon des signaux temporels d'impulsions de rÃ©fÃ©rence ceci pour quatre scintillateurs organiques. 5 Bien entendu, pour chacun de ces scintillateurs les paramÃ¨tres a et p ont Ã©tÃ© adaptÃ©s afin de prendre en compte les caractÃ©ristiques temporelles de chacun des scintillateurs. La procÃ©dure mis en place pour une telle adaptation est prÃ©cisÃ©e par la suite en lien avec le scintillateur BC-501-A. 10 On peut ainsi noter que si les taux de signaux temporels d'impulsions prÃ©sentent une variation relativement importante, le taux estimÃ© est Ã chaque fois significatif vis-Ã -vis de celui, reprÃ©sentatif du taux d'erreur, obtenu pour des signaux temporels d'impulsions obtenus pour une source de rayonnement ionisant dont le rayonnement ionisant ne comporte que des photons y. 15 Dans le cadre de ce mÃªme premier exemple de mise en oeuvre, il a Ã©galement pu Ãªtre montrÃ© qu'un nombre de 2000 signaux temporels d'impulsions Ã©tant suffisant pour estimer un taux de signaux temporels d'impulsions produits par l'interaction de neutrons avec une marge d'erreur de 4% et que cette marge d'erreur devient infÃ©rieure Ã 1% pour un nombre de 8000 signaux temporels d'impulsions. De plus 20 mÃªme avec un nombre rÃ©duit de signaux temporels d'impulsions, Ã savoir 100 signaux temporels d'impulsions, il a Ã©tÃ© possible d'obtenir une valeur caractÃ©ristique de dissemblance significative dÃ©montrant la prÃ©sence de neutrons dans le rayonnement ionisant Ã©mis pour la source de rayonnement ionisant comportant l'isotope 252 du californium 252Cf. 25 On peut Ã©galement noter que les rÃ©sultats de taux temporels d'impulsions obtenues pour la source de rayonnement ionisant comportant l'isotope 252 du californium 252Cf, et donc la valeur caractÃ©ristique de dissemblance correspondante, et 3065814 39 celle obtenu pour le sous-Ã©chantillonnage des signaux temporels d'impulsion de rÃ©fÃ©rence, et donc la valeur caractÃ©ristique de dissemblance correspondante, peuvent Ãªtre utilisÃ©s pour dÃ©finir le premier seuil prÃ©dÃ©terminÃ©. Ainsi dans le cadre de ce premier exemple de mise en oeuvre, le premier seuil prÃ©dÃ©terminÃ© peut Ãªtre choisi compris entre 5 la valeur caractÃ©ristique de dissemblance obtenue pour la source de rayonnement ionisant comportant l'isotope 252 du californium 252Cf et celle obtenue pour le sous-Ã©chantillonnage des signaux temporels d'impulsion de rÃ©fÃ©rence. Les figures 6A et 6B illustrent l'Ã©tape de choix de l'ensemble de grandeurs caractÃ©ristiques dans le cadre de ce premier exemple de mise en oeuvre avec 10 en figure 6A le calcul de la deuxiÃ¨me valeur caractÃ©ristique de dissemblance et en figure 6B le calcul de la troisiÃ¨me valeur caractÃ©ristique de dissemblance. Dans le cadre de ce premier exemple de mise en oeuvre, chaque ensemble de grandeurs discriminantes est constituÃ© de la grandeur discriminante (lot et de la grandeur discriminante Qtad et la distinction entre les ensembles de grandeurs 15 discriminantes est fournie par une variation des paramÃ¨tres a et R. Ainsi, comme le montre les figures 6A et 6B, le paramÃ¨tre a a Ã©tÃ© variÃ© entre 0,5 est 12 et le paramÃ¨tre (3 a Ã©tÃ© variÃ© entre 4 et 100. Une fois les ensembles de grandeurs discriminantes ainsi construits, le choix de l'au moins une grandeur discriminante parmi ces ensembles a Ã©tÃ© rÃ©alisÃ© au 20 moyens de 16000 signaux temporels d'impulsions obtenus Ã partir de la source de rayonnement ionisant comportant l'isotope 252 du californium 252Cf, et de 8000 signaux temporels d'impulsions obtenus par un sous-Ã©chantillonnage alÃ©atoire des signaux de rÃ©fÃ©rence. Ainsi, lesdits 16000 signaux temporels d'impulsions forment les signaux temporels d'impulsion de comparaison et lesdits 8000 signaux temporels d'impulsions 25 forment des signaux temporels d'impulsions de test. A partir de ces signaux temporels d'impulsions de rÃ©fÃ©rence, en application du procÃ©dÃ© selon l'invention et de maniÃ¨re similaire au calcul de la valeur caractÃ©ristique de dissemblance explicitÃ©e en rÃ©fÃ©rence aux figure 4A Ã 5B, des premiÃ¨res valeurs de comparaison sont calculÃ©es pour chacune des grandeurs discriminantes et 30 pour chacun des ensembles de grandeurs discriminantes. De la mÃªme faÃ§on, Ã partir des 3065814 signaux temporels d'impulsion de comparaison et de respectivement les signaux temporels d'impulsion de test, sont calculÃ©es des deuxiÃ¨mes valeurs de comparaison, respectivement des troisiÃ¨mes valeurs de comparaison, pour chacune des grandeurs discriminantes et pour chacun des ensembles de grandeurs discriminantes. 5 Afin d'obtenir les deuxiÃ¨mes valeurs caractÃ©ristique de dissemblance, un calcul de dissemblance, selon le principe dÃ©jÃ dÃ©crit, a Ã©tÃ© rÃ©alisÃ© Ã partir des premiÃ¨res et deuxiÃ¨mes valeurs de comparaison ceci pour chacun des ensembles de grandeurs discriminantes. Ces valeurs caractÃ©ristiques de dissemblance permettent d'estimer un taux de signaux temporels d'impulsion produits par l'interaction de neutrons 10 parmi les signaux temporels d'impulsions de comparaison pour chacun des ensembles de grandeurs discriminantes. Les taux de signaux temporels d'impulsion produits par l'interaction de neutrons ainsi estimÃ©s sont figurÃ©s graphiquement sur la figure 6A ceci en fonction des valeurs des paramÃ¨tres a et F3 pour lesdits ensembles de grandeurs discriminantes correspondant. Bien entendu, il est prÃ©fÃ©rable d'obtenir un taux maximal, 15 celui-ci indiquant qu'un nombre maximal de signaux temporels d'impulsions produits par l'interaction de neutrons ont pu Ãªtre dÃ©tectÃ©s et permettant une discrimination optimisÃ©e des sources de rayonnement ionisant dont le rayonnement ionisant comporte des neutrons et des photons y. De maniÃ¨re identique et afin d'obtenir les troisiÃ¨mes valeurs 20 caractÃ©ristiques de dissemblance, un calcul de dissemblance a Ã©tÃ© rÃ©alisÃ© Ã partir des premiÃ¨res et troisiÃ¨mes valeurs de comparaison ceci pour chacun des ensembles de grandeurs discriminantes, tout en prenant en compte le ratio de signaux temporels d'impulsions entre les signaux temporels d'impulsions de rÃ©fÃ©rence et les signaux temporels d'impulsions de test. Ces valeurs caractÃ©ristiques de dissemblance permettent 25 d'estimer un taux de signaux temporels d'impulsion produits par l'interaction de neutrons parmi les signaux temporels d'impulsions de test pour chacun des ensembles de grandeurs discriminantes. Les taux de signaux temporels d'impulsion produits par l'interaction de neutrons ainsi estimÃ©s sont figurÃ©s graphiquement sur la figure 6B ceci en fonction des valeurs des paramÃ¨tres a et F3 pour lesdits ensembles de grandeurs 30 discriminantes correspondant. Ici, puisque les signaux temporels d'impulsions de test ne 3065814 41 comportent que des photons y, il est nÃ©cessaire de minimiser le taux de signaux temporels d'impulsions produits par l'interaction de neutrons. Ainsi, le choix de l'ensemble de grandeurs discriminantes consiste donc Ã choisir les paramÃ¨tres a et p permettant de maximiser le taux de signaux temporels 5 d'impulsions produits par l'interaction de neutrons pour les signaux temporels d'impulsions de comparaison tout en assurant un taux de signaux temporels d'impulsions produits par l'interaction de neutrons pour les signaux temporels d'impulsions de test ne dÃ©passe un deuxiÃ¨me seuil prÃ©dÃ©terminÃ©, Ã savoir 3%. Dans l'exemple illustrÃ© sur les figures 6A et 6B, le couples de paramÃ¨tres a et (3 (4, 20) permet notamment de rÃ©pondre 10 Ã ces contraintes. On notera que ce graphique montre que d'autres couples des valeurs sont Ã©galement envisageables. Ainsi, avec un tels choix de l'au moins une grandeur discriminante constituÃ© des grandeurs discriminantes (-lot et Qtail pour le couple de paramÃ¨tre a et p (4,20), on obtient un systÃ¨me de caractÃ©risation offrant une discrimination optimisÃ©e 15 entre les sources de rayonnement ionisant Ã©mettant des rayonnements ionisants comportant des photons y et des neutrons vis-Ã -vis de celles dont le rayonnement ionisant ne comporte que des photons y. La figure 7 illustre le rÃ©sultat d'un calcul intermÃ©diaire selon un deuxiÃ¨me exemple de mise en oeuvre dans lequel la grandeur discriminante est la 20 grandeur discriminante D, les paramÃ¨tres a et p Ã©tant prÃ©sentant les valeurs 4, 20 dÃ©terminÃ©s pour le premier exemple de mise en oeuvre. Selon le mÃªme principe dÃ©crit dans le cadre du premier exemple de mise en oeuvre et Ã partir des mÃªmes signaux temporels d'impulsion de rÃ©fÃ©rence et acquis, il est possible de calculer les valeurs de rÃ©fÃ©rence et les valeurs caractÃ©ristiques de 25 la grandeur discriminante D sur la base du calcul dÃ©jÃ dÃ©crit. Ces 16000 valeurs de rÃ©fÃ©rence et ces 16000 valeurs de caractÃ©ristiques sont ensuite respectivement utilisÃ©es pour constituer une matrice de rÃ©fÃ©rence et une matrice caractÃ©ristique. La matrice de rÃ©fÃ©rence est ainsi une matrice comportant une colonne correspondant Ã D et 16000 lignes correspondant aux signaux temporels d'impulsions. 3065814 42 Afin de permettre un calcul de la valeur caractÃ©ristique de dissemblance, l'intervalle reprÃ©sentatif des valeurs prises par D a Ã©tÃ© subdivisÃ© en 200 sous-intervalles. Pour chacune de la matrice de rÃ©fÃ©rence et de la matrice 5 caractÃ©ristique, 200 coefficients 0 correspondant chacun Ã un sous-intervalle donnÃ© de D ont Ã©tÃ© calculÃ©s. Chacun des coefficients 0 est Ã©gal au nombre de signaux temporels d'impulsions de ladite matrice sont dans le sous-intervalle de D correspondant. Les 200 coefficients 0 calculÃ©s Ã partir de la matrice de rÃ©fÃ©rence permettent ainsi d'obtenir une matrice de 200 par 1 De maniÃ¨re identique, les 200 coefficients 0 calculÃ©s Ã partir de la 10 matrice caractÃ©ristique permettent ainsi d'obtenir une matrice de 200 par 1. La figure 7 illustre graphiquement les valeurs de coefficients 0 obtenues pour la matrice de rÃ©fÃ©rence 362 et pour la matrice caractÃ©ristique 361, avec en abscisse les valeurs des sous-intervalles de D et en ordonnÃ©e les valeurs des coefficients O. On peut ainsi voir sur cette figure que dans le cas des valeurs de rÃ©fÃ©rence, les grandeurs 15 discriminantes se prÃ©sentent sous la forme d'un unique pic avec une forte amplitude de 4000. Pour la matrice caractÃ©ristique, on note une variance plus forte des grandeurs discriminante. En effet, on remarque que ces grandeurs discriminantes sont rÃ©parties sous la forme d'un premier pic d'amplitude relativement importante de 2700, et un pic secondaire d'une amplitude plus faible atteignant 250. 20 Une fois les coefficients 0 calculÃ©s Ã la fois pour la matrice de rÃ©fÃ©rence et pour la matrice caractÃ©ristique, le calcul de la premiÃ¨re valeur caractÃ©ristique de dissemblance est rÃ©alisÃ© Ã partir d'une diffÃ©rence et d'une somme. Ce calcul de la premiÃ¨re valeur caractÃ©ristique, dont le rÃ©sultat d'une Ã©tape intermÃ©diaire est figurÃ© par la figure 7, consiste Ã faire successivement : 25 - pour chacun des sous-intervalles de la grandeur discriminante D, une diffÃ©rence en valeur absolue entre les coefficients 0 calculÃ©s pour la matrice de rÃ©fÃ©rence et pour la matrice caractÃ©ristique, le rÃ©sultat de cette diffÃ©rence Ã©tant figurÃ© sous la forme de portion hachurÃ©e 363 sur la figure 7, - une somme du rÃ©sultat des diffÃ©rences obtenues pour chacun des 30 intervalles de la grandeur discriminante D. 3065814 43 A partir de ce calcul de la valeur caractÃ©ristique de dissemblance, de la mÃªme maniÃ¨re que pour le premier exemple de mise en oeuvre, il est possible d'estimer le taux de signaux temporels d'impulsions produits par l'interaction de neutron. En effet, le taux de signaux temporels d'impulsion produits par l'interaction de neutrons est 5 directement proportionnel Ã la valeur caractÃ©ristique de dissemblance calculÃ©e. Dans la configuration utilisÃ©e pour cette mise en oeuvre, c'est-Ã -dire avec une source de rayonnement ionisant dont le rayonnement ionisant comporte des photons gamma y et des neutrons qui est une source de rayonnement ionisant comprenant l'isotope 252 du Californium 252Cf, il a Ã©tÃ© estimÃ© que 27,2% des 16000 signaux temporels d'impulsion 10 Ã©taient des signaux temporels d'impulsion produits par l'interaction de neutrons. Ce rÃ©sultat est donc trÃ¨s proche de celui obtenu dans le cadre de la premiÃ¨re du premier exemple de mise en oeuvre alors que ce deuxiÃ¨me exemple de mise en oeuvre n'est basÃ© que sur une seule grandeur discriminante. REVENDICATIONS1. ProcÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser dont le rayonnement ionisant comporte un premier type de particules et est susceptible de comporter un deuxiÃ¨me type de particules, le procÃ©dÃ© comportant les Ã©tapes suivantes : - fourniture de N signaux temporels d'impulsions de rÃ©fÃ©rence, chaque signal temporel d'impulsion de rÃ©fÃ©rence correspondant Ã la dÃ©tection du premier type de particules par un scintillateur, N Ã©tant un entier naturel non nul ; - choix d'au moins une grandeur discriminante ; - pour chaque signal temporel d'impulsion de rÃ©fÃ©rence, extraction de la valeur de chaque grandeur discriminante choisie pour former une valeur de rÃ©fÃ©rence correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de rÃ©fÃ©rence ; - acquisition de M signaux temporels d'impulsions acquis, chaque signal temporel d'impulsion acquis correspondant Ã la dÃ©tection, par le scintillateur, d'un photon issu de la source de rayonnement ionisant Ã caractÃ©riser, M Ã©tant un entier naturel non nul ; - pour chaque signal temporel d'impulsion acquis' extraction de la valeur de chaque grandeur discriminante choisie pour former une valeur caractÃ©ristique correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion acquis ; - calcul d'une premiÃ¨re valeur caractÃ©ristique de dissemblance Ã partir d'un premier calcul de dissemblance entre une matrice de rÃ©fÃ©rence et une matrice caractÃ©ristique, la matrice de rÃ©fÃ©rence Ã©tant formÃ©e par les valeurs de rÃ©fÃ©rence, la matrice caractÃ©ristique Ã©tant formÃ©e par les valeurs caractÃ©ristiques ; - comparaison de la premiÃ¨re valeur caractÃ©ristique de dissemblance calculÃ©e Ã un premier seuil prÃ©dÃ©terminÃ© ; 3065814 45 - caractÃ©risation de la prÃ©sence du deuxiÃ¨me type de particules dans un rayonnement Ã©mis par la source de rayonnement ionisant Ã caractÃ©riser en fonction du rÃ©sultat de la comparaison entre la premiÃ¨re valeur caractÃ©ristique de dissemblance calculÃ©e et le premier seuil prÃ©dÃ©terminÃ©, 5 dans lequel l'Ã©tape de choix d'au moins une grandeur discriminante comporte le choix de R grandeurs discriminantes, R Ã©tant un entier naturel non nul, et dans lequel le premier de dissemblance entre la matrice de rÃ©fÃ©rence et la matrice caractÃ©ristique comporte les sous-Ã©tapes suivantes : - pour chaque grandeur discriminante identifiÃ©e par un indice i, 10 subdivision d'un intervalle reprÃ©sentatif des valeurs prises par la grandeur discriminante en k, sous-intervalles, k, Ã©tant un entier naturel non nul ; - calcul, pour chacune de la matrice de rÃ©fÃ©rence et la matrice caractÃ©ristique, de la valeur de Trh ki coefficients 0, chaque coefficient 0 Ã©tant associÃ© Ã R sous-intervalles, chacun desdits R sous-intervalles Ã©tant associÃ© Ã une grandeur 15 discriminante respective et Ã©tant Ã©gal Ã un sous-intervalle de ladite grandeur discriminante respective, un premier coefficient 01 et un deuxiÃ¨me coefficient 02 distincts Ã©tant tels qu'il existe au moins une grandeur caractÃ©ristique pour laquelle le sous-intervalle correspondant du premier coefficient 01 est diffÃ©rent du sous-intervalle correspondant du deuxiÃ¨me coefficient 02, la valeur de chacun des coefficients 0 Ã©tant 20 calculÃ©e comme le nombre de signaux temporels d'impulsions de rÃ©fÃ©rence, respectivement de signaux temporels d'impulsions acquis, dont les valeurs de grandeurs discriminantes extraites sont dans les P sous-intervalles associÃ©s audit coefficient, les valeur de coefficients 0 calculÃ©es pour l'une de la matrice de rÃ©fÃ©rence et de la matrice caractÃ©ristique Ã©tant multipliÃ©es par un facteur correctif correspondant au ration du 25 nombre de signaux temporels d'impulsions entre la matrice de rÃ©fÃ©rence et la matrice caractÃ©ristique ; - calcul de la premiÃ¨re valeur caractÃ©ristique de dissemblance Ã partir d'une diffÃ©rence en valeur absolue mettant en oeuvre les coefficients 0 calculÃ©s Ã partir de la matrice de rÃ©fÃ©rence et de la matrice caractÃ©ristique et de R sommes successives, chacune des R sommes Ã©tant associÃ©e Ã une grandeur discriminante respective, et ladite 3065814 46 diffÃ©rence se faisant directement entre chaque coefficient 0 calculÃ© sur la base de la matrice de rÃ©fÃ©rence et le coefficient 0 calculÃ© sur la base de la matrice caractÃ©ristique correspondant aux mÃªmes R sous-intervalles, la premiÃ¨re somme se faisant sur la base du rÃ©sultat de ladite diffÃ©rence et les R-1 sommes suivantes sur la base du rÃ©sultat de la 5 somme la prÃ©cÃ©dent directement. 2. ProcÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser selon la revendication 1, dans lequel au moins une grandeur discriminante est sÃ©lectionnÃ©e parmi : 10 - une intÃ©grale d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis; - un temps caractÃ©ristique du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis; - une valeur, en au moins un instant, du signal temporel d'impulsion de 15 rÃ©fÃ©rence ou signal temporel d'impulsion acquis; - une fonction de la variance d'au moins une portion du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis; - une valeur de dÃ©rivÃ©e du signal temporel d'impulsion de rÃ©fÃ©rence ou signal temporel d'impulsion acquis. 20 3. ProcÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser selon la revendication 1 ou 2, dans lequel l'Ã©tape de choix d'au moins une grandeur discriminante est postÃ©rieure Ã l'Ã©tape de fourniture des N signaux temporels d'impulsions de rÃ©fÃ©rence et Ã l'Ã©tape de normalisation de chaque signal temporel 25 d'impulsion de rÃ©fÃ©rence, l'Ã©tape de choix d'au moins une grandeur discriminante comportant les sous-Ã©tapes suivantes : - fourniture de P signaux temporels d'impulsions de comparaison, chaque signal temporel d'impulsion de comparaison correspondant Ã la dÃ©tection, par le 30 scintillateur, d'une particule du premier type ou du deuxiÃ¨me type Ã©mise par une source 3065814 47 de rayonnement ionisant Ã©mettant un rayonnement ionisant comportant des particules du premier type et des particules du deuxiÃ¨me type, P Ã©tant un entier naturel non nul ; - construction d'une pluralitÃ© d'ensembles de grandeurs discriminantes ; - pour chaque signal temporel d'impulsion de rÃ©fÃ©rence, et pour chaque 5 ensemble de grandeurs discriminantes, extraction de la valeur de chaque grandeur discriminante dudit ensemble de grandeurs discriminantes pour former une premiÃ¨re valeur de comparaison correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de rÃ©fÃ©rence ; - pour chaque signal temporel d'impulsion de comparaison, et pour 10 chaque ensemble de grandeurs discriminantes, extraction de la valeur de chaque grandeur discriminante dudit ensemble de grandeurs discriminantes pour former une deuxiÃ¨me valeur de comparaison correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de comparaison ; - pour chaque ensemble de grandeurs discriminantes, calcul d'une 15 deuxiÃ¨me valeur caractÃ©ristique de dissemblance Ã partir d'un deuxiÃ¨me calcul de dissemblance entre une premiÃ¨re matrice constituÃ©e des premiÃ¨res valeurs de comparaison extraites pour ledit ensemble de grandeurs discriminantes et une deuxiÃ¨me matrice constituÃ©e des deuxiÃ¨mes valeurs de comparaison extraites pour ledit ensemble de grandeurs discriminantes ; 20 - choix d'au moins une grandeur discriminante correspondant Ã l'un des ensembles de grandeurs discriminantes, le choix Ã©tant rÃ©alisÃ© sur la base des deuxiÃ¨mes valeurs caractÃ©ristique de dissemblance calculÃ©e pour chacun des ensembles de grandeurs discriminantes. 25 4. ProcÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser selon la revendication 3, dans lequel lors de la sous Ã©tape de choix d'au moins une grandeur discriminante correspondant Ã l'un des ensembles de grandeurs discriminantes, l'au moins une grandeur discriminante correspondant choisie est celle qui correspond Ã l'ensemble de grandeurs discriminantes pour lequel la deuxiÃ¨me valeur 3065814 48 caractÃ©ristique de dissemblance prÃ©sente un extrÃ©mum caractÃ©ristique d'une dissemblance maximale. 5. ProcÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã 5 caractÃ©riser selon la revendication 3, dans lequel l'Ã©tape de choix d'au moins une grandeur discriminante comporte en outre les sous-Ã©tapes suivantes : - fourniture de Q signaux temporels d'impulsions de test, chaque signal temporel d'impulsion de test correspondant Ã la dÃ©tection, par le scintillateur, d'une particule du premier type Ã©mis par une source de rayonnement ionisant Ã©mettant un 10 rayonnement ionisant comportant des particules du premier type, Q Ã©tant un entier naturel non nul ; - pour chaque signal temporel d'impulsion de test, et pour chaque ensemble de grandeurs discriminantes, extraction de la valeur de chaque grandeur discriminante dudit ensemble de grandeurs discriminantes pour former une troisiÃ¨me 15 valeur de comparaison correspondant Ã ladite grandeur discriminante et audit signal temporel d'impulsion de test ; - pour chaque ensemble de grandeurs discriminantes, calcul d'une troisiÃ¨me valeur caractÃ©ristique de dissemblance Ã partir du deuxiÃ¨me calcul de dissemblance entre une premiÃ¨re matrice constituÃ©e des premiÃ¨res valeurs de 20 comparaison extraites pour ledit ensemble de grandeurs discriminantes et une deuxiÃ¨me matrice constituÃ©e des troisiÃ¨me valeurs de comparaison extraites pour ledit ensemble de grandeurs discriminantes ; - comparaison de la troisiÃ¨me valeur caractÃ©ristique de dissemblance calculÃ©e pour chacun des ensembles de grandeurs discriminantes Ã un deuxiÃ¨me seuil 25 prÃ©dÃ©terminÃ© ; - sÃ©lection parmi les ensembles de grandeurs discriminantes de ceux qui prÃ©sentent une dissemblance minimale Ã partir de ladite comparaison entre la troisiÃ¨me valeur caractÃ©ristique de dissemblance calculÃ©e pour chacun des ensembles de grandeurs discriminantes Ã un deuxiÃ¨me seuil prÃ©dÃ©terminÃ© ; 3065814 49 - choix d'au moins une grandeur discriminante correspondant Ã l'un des ensembles de grandeurs discriminantes, l'au moins une grandeur discriminante correspondant choisie Ã©tant celle correspondant Ã l'ensemble de grandeurs discriminantes parmi les ensembles de grandeurs discriminantes sÃ©lectionnÃ©s dont la 5 deuxiÃ¨me valeur caractÃ©ristique de dissemblance prÃ©sente un extrÃ©mum caractÃ©ristique d'une dissemblance maximale. 6. ProcÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser selon l'une quelconque des revendications 1 Ã 5, dans lequel au moins une 10 grandeur discriminante est notÃ©e (lot et est calculÃ©e selon : ti Qt0t = ftmaxS dt oÃ¹ s est le signal temporel d'impulsion de rÃ©fÃ©rence ou le signal temporel d', respectivement le signal temporel d'impulsion acquis; tmax est la coordonnÃ©e temporelle du maximum du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis; et 15 t1 est une coordonnÃ©e temporelle calculÃ©e selon : t1= tmax+13.tHm oÃ¹ p est un rÃ©el strictement positif ; et tHM est la coordonnÃ©e temporelle, ultÃ©rieure Ã tmax, pour laquelle la valeur du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, vaut 0,5. 20 7. ProcÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser selon l'une quelconque des revendications 1 Ã 6, dans lequel au moins une grandeur discriminante est notÃ©e Qtad et calculÃ©e selon : Qtall = f S(t) dt to oÃ¹ s est le signal temporel d'impulsion de rÃ©fÃ©rence, respectivement le 25 signal temporel d'impulsion acquis; 3065814 50 tic, est une coordonnÃ©e temporelle calculÃ©e selon : -4= tmax+a.tHm oÃ¹ a est un rÃ©el strictement positif ; tmax est la coordonnÃ©e temporelle du maximum du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal 5 temporel d'impulsion acquis; et tHM est la coordonnÃ©e temporelle, ultÃ©rieure Ã tmax, pour laquelle la valeur du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, vaut 0,5 ; et t1 est une coordonnÃ©e temporelle calculÃ©e selon : t1= tmax+13.tHm 10 oÃ¹ p est un rÃ©el strictement positif. 8. ProcÃ©dÃ© selon les revendications 3, 6 et 7 seules ou en combinaison avec la revendication 4 ou 5, dans lequel chaque ensemble de grandeurs discriminantes comporte l'un au moins de (-lot et Qtail, les ensembles se distinguant deux Ã deux par au 15 moins l'une des valeurs a et R. 9. ProcÃ©dÃ© selon l'une quelconque des revendications 1 Ã 8, dans lequel au moins une grandeur discriminante est notÃ©e D et est la variance de la portion du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement le signal temporel d'impulsion acquis, 20 comprise entre tic, et t1, oÃ¹ : tic, est une coordonnÃ©e temporelle calculÃ©e selon : -4= tmax+a.tHm oÃ¹ a est un rÃ©el strictement positif ; tmax est la coordonnÃ©e temporelle du maximum du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal 25 temporel d'impulsion acquis; et 3065814 51 tHM est la coordonnÃ©e temporelle, ultÃ©rieure Ã tmax, pour laquelle la valeur du signal temporel d'impulsion de rÃ©fÃ©rence, respectivement du signal temporel d'impulsion acquis, vaut 0,5 ; t1 est une coordonnÃ©e temporelle calculÃ©e selon : t1= tmax+13.tHm 5 oÃ¹ p est un rÃ©el strictement positif. 10. ProcÃ©dÃ© selon les revendications 3 et 9 seules ou en combinaison avec l'une quelconque des revendications 4 Ã 8, dans lequel chaque ensemble de grandeurs discriminantes comporte D, les ensembles se distinguant deux Ã deux par au 10 moins l'une des valeurs a et R. 11. ProcÃ©dÃ© selon la revendication 3 seule ou prise en combinaison avec l'une quelconque des revendications 4 Ã 10, dans lequel, lors de l'Ã©tape de construction de la pluralitÃ© d'ensembles de grandeurs discriminantes, un ensemble donnÃ© de 15 grandeurs discriminantes comporte S grandeurs discriminantes, S Ã©tant un entier naturel non nul, et dans lequel le deuxiÃ¨me calcul de dissemblance entre une premiÃ¨re matrice et une deuxiÃ¨me matrice, pour l'ensemble donnÃ© de S grandeurs discriminantes, comporte les sous-Ã©tapes suivantes : - pour chaque grandeur discriminante identifiÃ©e par un indice j, 20 subdivision d'un intervalle reprÃ©sentatif des valeurs prises par la grandeur discriminante en k, sous-intervalles, k, Ã©tant un entier naturel non nul ; - calcul, pour chacune de la premiÃ¨re matrice et de la deuxiÃ¨me matrice, de la valeur de FIs1kj coefficients 6, chaque coefficient 6 Ã©tant associÃ© Ã S sous- intervalles, chacun desdits S sous-intervalles Ã©tant associÃ© Ã une grandeur discriminante 25 respective et Ã©tant Ã©gal Ã un sous-intervalle de ladite grandeur discriminante respective, un premier coefficient 61 et un deuxiÃ¨me coefficient 62 distincts Ã©tant tels qu'il existe au moins une grandeur caractÃ©ristique pour laquelle le sous-intervalle correspondant du premier coefficient 61 est diffÃ©rent du sous-intervalle correspondant du deuxiÃ¨me coefficient 62, la valeur de chacun des coefficients 6 Ã©tant calculÃ©e comme le nombre de 3065814 52 signaux temporels d'impulsions de la premiÃ¨re matrice, respectivement de signaux temporels d'impulsions de la deuxiÃ¨me matrice, dont les valeurs de grandeurs discriminantes extraites sont dans les S sous-intervalles associÃ©s audit coefficient, les valeurs de coefficients 6 calculÃ©s pour l'une de la premiÃ¨re matrice et de la deuxiÃ¨me 5 matrice Ã©tant multipliÃ©es par un facteur correctif correspondant au ration du nombre signaux temporels d'impulsions entre la premiÃ¨re matrice et la deuxiÃ¨me matrice ; - calcul d'une valeur caractÃ©ristique de dissemblance Ã partir de S sommes et d'une diffÃ©rence en valeur absolue mettant en oeuvre les coefficients 6 calculÃ©s Ã partir de la premiÃ¨re matrice et de la deuxiÃ¨me matrice, chacune des S sommes 10 Ã©tant associÃ©e Ã une grandeur discriminante respective, et ladite diffÃ©rence se faisant Ã partir de l'ensemble des rÃ©sultats de l'une desdites S sommes, entre les rÃ©sultats obtenus Ã partir des coefficients 6 calculÃ©s sur la base de la premiÃ¨re matrice et les rÃ©sultats obtenus Ã partir des coefficients 6 calculÃ©s sur la base de la deuxiÃ¨me matrice, le reste des S sommes rÃ©alisÃ©es aprÃ¨s le calcul de la diffÃ©rence se faisant sur la base du rÃ©sultat de 15 ladite diffÃ©rence. - calcul d'une valeur caractÃ©ristique de dissemblance Ã partir d'une diffÃ©rence en valeur absolue mettant en oeuvre les coefficients 6 calculÃ©s Ã partir de la premiÃ¨re matrice et de la deuxiÃ¨me matrice et de S sommes successives, chacune des S sommes Ã©tant associÃ©e Ã une grandeur discriminante respective, et ladite diffÃ©rence se 20 faisant directement entre chaque coefficient 6 calculÃ© sur la base de la premiÃ¨re matrice et le coefficient 6 calculÃ© sur la base de la deuxiÃ¨me matrice correspondant au mÃªme R sous-intervalles, la premiÃ¨re somme se faisant sur la base du rÃ©sultat de ladite diffÃ©rence et les S-1 sommes suivantes sur la base du rÃ©sultat de la somme la prÃ©cÃ©dent directement. 25 12. ProcÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser selon l'une quelconque des revendication 1 Ã 11, dans lequel il est en outre prÃ©vu une Ã©tape d'estimation du nombre et/ou d'un taux d'impulsions correspondant au deuxiÃ¨me type de particules. 30 3065814 53 13. procÃ©dÃ© de caractÃ©risation d'une source de rayonnement ionisant Ã caractÃ©riser selon l'une quelconque des revendications 1 Ã 12, dans lequel le premier type de particules et le deuxiÃ¨me type de particules sont sÃ©lectionnÃ©s parmi le groupe comprenant : 5 - les photons gamma en tant que premier type de particules et les neutrons en tant que deuxiÃ¨me type de particules, - les photons gamma en tant que premier type de particules et les particules beta en tant que deuxiÃ¨me type de particules, - les photons gamma en tant que premier type de particules et les 10 particules alpha en tant que deuxiÃ¨me type de particules, - les particules beta en tant que premier type de particules et les particules alpha en tant que deuxiÃ¨me type de particules. 14. SystÃ¨me de caractÃ©risation (10) d'une source de rayonnement 15 ionisant Ã caractÃ©riser comportant : - un scintillateur (11), - un photomultiplicateur (12) couplÃ© au scintillateur (11), et - une unitÃ© de traitement (13) adaptÃ©e pour traiter les signaux temporels d'impulsions fournis par le photomultiplicateur (12) et obtenus lors de 20 l'absorption de rayonnements ionisants par le scintillateur (11), le systÃ¨me de caractÃ©risation (10) Ã©tant en ce que l'unitÃ© de traitement (13) est configurÃ©e pour mettre en oeuvre un procÃ©dÃ© de caractÃ©risation selon l'une quelconque des revendication 1 Ã 13. 25 15. Support de donnÃ©es lisible par un ordinateur sur lequel est enregistrÃ© un programme incluant les instructions pour exÃ©cuter des Ã©tapes du procÃ©dÃ© selon l'une quelconque des revendications 1 Ã 13.