Brevet FR3086764A1 — publié le 20200403

SYSTEME DE SPECTROMETRIE, PROCEDE DE SPECTROMETRIE ET PRODUIT PROGRAMME D'ORDINATEUR ASSOCIES DESCRIPTION DOMAINE TECHNIQUE La présente invention concerne un système de spectrométrie comprenant un étage de détection comportant un détecteur principal et au moins un détecteur auxiliaire, le détecteur principal étant configuré pour détecter au moins des photons présentant une énergie appartenant à une plage d'énergie principale dans le domaine X et/ou Y, et pour délivrer un signal de détection principal comportant des impulsions principales, chaque impulsion principale étant représentative d'une interaction du détecteur principal avec un photon de la plage d'énergie principale, chaque détecteur auxiliaire étant configuré pour détecter au moins des photons présentant une énergie appartenant à une plage d'énergie auxiliaire dans le domaine X et/ou Y, et pour délivrer un signal de détection auxiliaire comportant des impulsions auxiliaires, chaque impulsion auxiliaire d'un signal de détection auxiliaire donné étant représentative d'une interaction du détecteur auxiliaire correspondant avec un photon de la plage d'énergie auxiliaire associée. L'invention concerne également un procédé de spectrométrie et un produit programme d'ordinateur. L'invention s'applique notamment à la détection et la mesure de faibles traces d'isotopes radioactifs, c'est-à-dire la détection d'isotopes radioactifs présentant une activité faible au regard d'un bruit de fond imposé par des conditions expérimentales données (par exemple du même ordre de grandeur), et la mesure de ladite activité. ÉTAT DE LA TECHNIQUE ANTÉRIEURE De façon connue, le processus d'atténuation des rayonnements X et Y dans la matière, par exemple dans un détecteur, est principalement régi par deux effets : 2 - l'effet photoélectrique, par lequel un photon incident dépose l'intégralité de son énergie dans le détecteur, ce qui se traduit, dans le signal électrique de détection issu du détecteur, par une impulsion dite « impulsion d'absorption totale », participant à un « pic d'absorption totale » dans le spectre d'un signal de détection issu du détecteur, exploitable dans le cadre de la caractérisation d'une source ayant émis ledit photon ; et - l'effet Compton, par lequel un photon incident dépose une partie seulement de son énergie dans le détecteur en le traversant, ce qui se traduit, dans le spectre du signal électrique de détection, par un continuum qui est assimilable à du bruit 10 et qui est susceptible de masquer, dans ledit spectre, des pics d'absorption totale de faible amplitude générés par effet photoélectrique. Il en résulte que l'effet Compton est indésirable, dans la mesure où il rend plus difficile la caractérisation de la source. Il est également connu d'utiliser des systèmes de spectrométrie dans 15 lesquels un détecteur auxiliaire est associé à un détecteur principal, et de construire un signal hybride à partir du signal de détection issu de chacun parmi le détecteur principal et le détecteur auxiliaire. Le signal hybride correspond au signal de détection du détecteur principal, duquel a été supprimée chaque impulsion telle que, à la même date, il existe une 20 impulsion coïncidente dans le signal de détection issu du détecteur auxiliaire. Ainsi est partiellement filtré le bruit de diffusion Compton imputable à l'activité de la source, c'est-à-dire imputable aux photons qui sont issus de la source et qui déposent une partie de leur énergie dans le détecteur principal et une autre partie de leur énergie dans le détecteur auxiliaire. 25 Néanmoins, un tel système de spectrométrie ne donne pas entière satisfaction. En effet, le signal de détection issu d'un système de spectrométrie X et/ou Y, dans lequel un détecteur principal et un détecteur auxiliaire sont associés, peut être décomposé en plusieurs composantes : 3 - le signal utile, qui est la composante produite par le dépôt d'énergie totale d'un photon émis par la source dans le détecteur principal (formant donc une impulsion d'absorption totale) ; - le bruit de diffusion de type I, qui comprend : - une composante I.a, représentative d'une première interaction d'un photon issu de la source avec le détecteur principal par diffusion Compton, puis d'une interaction secondaire avec le détecteur auxiliaire ; - une composante I.b, représentative d'une première 10 interaction d'un photon issu de la source avec le détecteur auxiliaire par diffusion Compton, puis d'une interaction secondaire avec le détecteur principal ; - une composante I.c, représentative d'une première interaction d'un photon issu de l'environnement avec le détecteur 15 principal par diffusion Compton, puis d'une interaction secondaire avec le détecteur auxiliaire ; - une composante I.d, représentative d'une première interaction d'un photon issu de l'environnement avec le détecteur auxiliaire par diffusion Compton, puis d'une interaction secondaire avec 20 le détecteur principal ; - le bruit de diffusion de type II, qui comprend : - une composante II.a, représentative d'une interaction d'un photon issu de la source avec le détecteur principal par diffusion Compton, sans interaction secondaire avec le détecteur auxiliaire ; 25 - une composante II.b, représentative d'une première interaction d'un photon issu de la source avec l'environnement par diffusion Compton, puis d'une interaction secondaire avec le détecteur principal ; 4 - une composante II.c, représentative d'une interaction d'un photon issu de l'environnement avec le détecteur principal, sans interaction secondaire avec le détecteur auxiliaire. En ce qui concerne la composante II.b, les photons sont considérés comme provenant de l'environnement, car ayant perdu leur corrélation à la source. Bien que le traitement fondé sur la construction du signal hybride permette de rejeter, c'est-à-dire filtrer, le bruit de diffusion de type I, les systèmes de spectrométrie connus ne sont pas en mesure de rejeter le bruit de diffusion de type II. Un but de l'invention est donc de proposer un système de spectrométrie 10 opérant dans le domaine X et/ou Y, et qui permette une meilleure réjection du bruit de diffusion, en particulier du bruit de diffusion de type II, que les systèmes de spectrométrie connus. EXPOSÉ DE L'INVENTION 15 A cet effet, l'invention a pour objet un système de spectrométrie du type précité, comportant, en outre, un étage de traitement comprenant : - un organe d'acquisition configuré pour acquérir le signal de détection principal en un signal de détection principal acquis, et chaque signal de détection auxiliaire en un signal de détection auxiliaire acquis correspondant ; 20 - un organe de calcul configuré pour : - détecter les impulsions principales dans le signal de détection principal acquis et les impulsions auxiliaires dans chaque signal de détection auxiliaire acquis ; - déterminer, pour chaque impulsion principale, une 25 énergie principale correspondante et, pour chaque impulsion auxiliaire, une énergie auxiliaire correspondante ; - calculer un premier spectre en fonction du signal de détection principal acquis et, éventuellement, de chaque signal de détection auxiliaire acquis ; 5 - calculer au moins un deuxième spectre en fonction d'au moins un signal de détection auxiliaire acquis et, éventuellement, du signal de détection principal acquis ; - une mémoire configurée pour stocker un modèle numérique associant 5 une énergie donnée d'un photon incident parvenant à l'étage de détection, à : - une énergie principale déposée par le photon incident dans le détecteur principal, dont la valeur dépend de ce que le photon incident est issu de la source à caractériser ou d'un environnement distinct de la source ; et 10 - une énergie auxiliaire déposée par le photon incident dans chaque détecteur auxiliaire, dont la valeur dépend de ce que le photon incident est issu de la source à caractériser ou de l'environnement ; l'organe de calcul étant, en outre, configuré pour calculer un vecteur 15 isotopique et/ou un flux de bruit comme étant chacun une solution d'un jeu d'équations comprenant : - une première équation dans laquelle le premier spectre s'exprime comme la somme du produit du vecteur isotopique par un premier paramètre du modèle numérique d'une part, et du produit du 20 flux de bruit par un deuxième paramètre du modèle numérique d'autre part ; et - au moins une deuxième équation dans laquelle chaque deuxième spectre s'exprime comme la somme du produit du vecteur isotopique par un troisième paramètre du modèle numérique 25 d'une part, et du produit du flux de bruit par un quatrième paramètre du modèle numérique d'autre part ; le vecteur isotopique étant un vecteur dont chaque composante d'indice i est représentative d'un nombre de photons émis par la source à caractériser, par unité de temps, dans un intervalle d'énergie E, associé à ladite composante d'indice i, lesdits 6 photons ayant interagi exclusivement avec le détecteur principal et ayant déposé entièrement leur énergie dans le détecteur principal, le flux de bruit étant un vecteur dont chaque composante d'indice i est représentative d'un nombre de photons, par unité de temps, provenant de l'environnement, dans un intervalle d'énergie Ei associé à ladite composante d'indice i, et incidents sur l'étage de détection. En effet, dans la mesure où le calcul du vecteur isotopique et du flux de bruit met en oeuvre des paramètres du modèle numérique, le système de spectrométrie selon l'invention autorise une réjection du bruit qui s'appuie sur un modèle complet dudit 10 système de spectrométrie, qui incorpore notamment les réponses du détecteur principal et de ou des détecteur(s) auxiliaire(s). La modélisation conjointe du détecteur principal et de ou des détecteur(s) auxiliaire(s) et de leur comportement autorise, en particulier, la réjection des bruits issus des diffusions par l'environnement du système de spectrométrie (composantes 15 II.b et II.c), réjection que les systèmes de spectrométrie de l'état de la technique ne sont pas aptes à réaliser. De façon avantageuse, le calcul du premier spectre en fonction du signal de détection principal acquis et de chaque signal de détection auxiliaire acquis autorise une rejection des composantes de bruit de diffusion de type I. 20 Ainsi, grâce au système de spectrométrie selon l'invention, une spectrométrie avec un niveau de bruit plus faible que celui autorisé par les systèmes de spectrométrie de l'état de la technique est possible, notamment grâce à un filtrage des composantes de bruit qui ne pouvaient être discriminées par les approches connues. Suivant d'autres aspects avantageux de l'invention, le système de 25 spectrométrie comporte une ou plusieurs des caractéristiques suivantes, prise(s) isolément ou suivant toutes les combinaisons techniquement possibles : - l'équation relative au premier signal s'écrit : no S1 = + WteBe o=1 7 où Si est le premier spectre, égal au spectre du signal de détection principal acquis ; Fi est une matrice dont un terme général -Fiji; est égal à la probabilité qu'un photon incident d'énergie E et issu de la source à caractériser induise un dépôt d'énergie E1 dans le détecteur principal ; A est le vecteur isotopique ; no est le nombre de détecteurs auxiliaires, chaque détecteur auxiliaire d'indice e étant associé à un secteur d'angles géométriques d'incidence propre noté AS ; Wise est une matrice dont un terme général est égal à la probabilité 10 qu'un photon provenant de l'environnement, incident sur l'étage de détection selon le secteur d'angles géométriques d'incidence AS, et présentant une énergie E induise un dépôt d'énergie E1 dans le détecteur principal ; et Be est la partie du flux de bruit correspondant à des photons incidents sur l'étage de détection selon le secteur d'angles géométriques d'incidence AS ; 15 - le calculateur est configuré pour calculer un signal hybride, le signal hybride étant obtenu à partir du signal de détection principal, duquel a été supprimée chaque impulsion principale telle que, à la même date, il existe une impulsion auxiliaire coïncidente dans au moins un signal de détection auxiliaire acquis, un deuxième spectre parmi l'au moins un deuxième spectre étant le spectre du signal hybride ; 20 - une deuxième équation parmi l'au moins une deuxième équation s'écrit : ne S'i = F'1A + Wri,oBa 0=1 où S'i est un deuxième spectre parmi l'au moins un deuxième spectre, égal au spectre du signal hybride ; 25 ri est une matrice dont un terme général fui est égal à la probabilité qu'un photon incident d'énergie E et issu de la source induise un dépôt d'énergie E1 dans le détecteur principal sachant qu'auront été rejetées les impulsions principales telles que, à la même date, il existe au moins une impulsion auxiliaire coïncidente ; A est le vecteur isotopique ; 8 no est le nombre de détecteurs auxiliaires, chaque détecteur auxiliaire d'indice 6 étant associé à un secteur d'angles géométriques d'incidence propre noté A6 ; W'i,e est une matrice dont un terme général est égal à la probabilité qu'un photon provenant de l'environnement, incident sur l'étage de détection selon le secteur d'angles géométriques d'incidence A,I3, et présentant une énergie E induise un dépôt d'énergie E dans le détecteur principal sachant qu'auront été rejetées les impulsions principales telles que, à la même date, il existe au moins une impulsion auxiliaire coïncidente ; et Be est la partie du flux de bruit correspondant à des photons incidents sur 10 l'étage de détection selon le secteur d'angles géométriques d'incidence Mi ; - une deuxième équation parmi l'au moins une deuxième équation s'écrit : 52B = F2A W2,eBe où 6 est un indice associé à un détecteur auxiliaire donné, le détecteur 15 auxiliaire d'indice 6 étant associé à un secteur d'angles géométriques d'incidence propre noté As ; S2,0 est le deuxième spectre, égal au spectre du signal auxiliaire acquis correspondant au détecteur auxiliaire d'indice 61; F2 est une matrice dont un terme général f2,,,, est égal à la probabilité 20 qu'un photon incident d'énergie E et issu de la source à caractériser induise un dépôt d'énergie E, dans le détecteur auxiliaire d'indice 6 ; A est le vecteur isotopique ; W2,13 est une matrice dont un terme général w2,,,,,0 est égal à la probabilité qu'un photon provenant de l'environnement, incident sur l'étage de détection selon le 25 secteur d'angles géométriques d'incidence AS, et présentant une énergie E, induise un dépôt d'énergie E dans le détecteur auxiliaire d'indice 6 ; et Be est la partie du flux de bruit correspondant à des photons incidents sur l'étage de détection selon le secteur d'angles géométriques d'incidence AS. En outre, l'invention a pour objet un procédé de spectrométrie comprenant les étapes suivantes : 9 - détection d'impulsions principales dans un signal de détection principal acquis, le signal de détection principal acquis étant le résultat de l'acquisition d'un signal de détection principal délivré par un détecteur principal, le détecteur principal étant configuré pour détecter au moins des photons présentant une énergie appartenant à une plage d'énergie principale dans le domaine X et/ou Y, et pour délivrer le signal de détection principal, le signal de détection principal comportant les impulsions principales, chaque impulsion principale étant représentative d'une interaction du détecteur principal avec un photon de la plage d'énergie principale ; - détection d'impulsions auxiliaires dans au moins un signal de détection 10 auxiliaire acquis, chaque signal de détection auxiliaire acquis étant le résultat d'un signal de détection auxiliaire correspondant délivré par un détecteur auxiliaire correspondant, chaque détecteur auxiliaire étant configuré pour détecter au moins des photons présentant une énergie appartenant à une plage d'énergie auxiliaire dans le domaine X et/ou Y, et pour délivrer le signal de détection auxiliaire correspondant, chaque signal de 15 détection auxiliaire comportant des impulsions auxiliaires, chaque impulsion auxiliaire d'un signal de détection auxiliaire donné étant représentative d'une interaction du détecteur auxiliaire correspondant avec un photon de la plage d'énergie auxiliaire associée ; - détermination, pour chaque impulsion principale, d'une énergie principale correspondante et, pour chaque impulsion auxiliaire, d'une énergie auxiliaire 20 correspondante ; - calcul d'un premier spectre en fonction du signal de détection principal acquis et, éventuellement, de chaque signal de détection auxiliaire acquis ; - calcul d'au moins un deuxième spectre en fonction d'au moins un signal de détection auxiliaire acquis et, éventuellement, du signal de détection principal acquis ; 25 - calculer d'un vecteur isotopique et d'un flux de bruit comme étant les solutions d'un jeu d'équations comprenant : - une première équation dans laquelle le premier spectre s'exprime comme la somme du produit du vecteur isotopique par un premier paramètre d'un modèle numérique prédéterminé d'une part, 10 et du produit du flux de bruit par un deuxième paramètre du modèle numérique d'autre part ; et - au moins une deuxième équation dans laquelle chaque deuxième spectre s'exprime comme la somme du produit du vecteur isotopique par un troisième paramètre du modèle numérique d'une part, et du produit du flux de bruit par un quatrième paramètre du modèle numérique d'autre part ; le vecteur isotopique étant un vecteur dont chaque composante d'indice i est représentative d'un nombre de photons émis par la source à caractériser, par unité de 10 temps, dans un intervalle d'énergie E, associé à ladite composante d'indice i, lesdits photons ayant interagi exclusivement avec le détecteur principal et ayant déposé entièrement leur énergie dans le détecteur principal, le flux de bruit étant un vecteur dont chaque composante d'indice i est représentative d'un nombre de photons, par unité de temps, provenant de 15 l'environnement, dans un intervalle d'énergie E, associé à ladite composante d'indice i, et incidents sur l'étage de détection, le modèle numérique associant une énergie donnée d'un photon incident à: - une énergie principale déposée par le photon incident 20 dans le détecteur principal, dont la valeur dépend de ce que le photon incident est issu d'une source à caractériser ou d'un environnement distinct de la source ; et - une énergie auxiliaire déposée par le photon incident dans chaque détecteur auxiliaire, dont la valeur dépend de ce que le 25 photon incident est issu de la source à caractériser ou de l'environnement. En outre, l'invention a pour objet un produit programme d'ordinateur comprenant des instructions de code de programme qui, lorsqu'elles sont exécutées par un ordinateur, mettent en oeuvre le procédé tel que défini ci-dessus. 30 11 BRÈVE DESCRIPTION DES DESSINS L'invention sera mieux comprise à l'aide de la description qui va suivre, donnée uniquement à titre d'exemple non limitatif et faite en se référant aux dessins annexés sur lesquels : 5 - la figure 1 est une représentation schématique d'un mode de réalisation d'un système de spectrométrie selon l'invention ; - la figure 2 est une représentation schématique d'une variante du système de spectrométrie de la figure 1. 10 EXPOSÉ DÉTAILLÉ DE MODES DE RÉALISATION PARTICULIERS Un système de spectrométrie 2 selon l'invention est représenté sur la figure 1. Le système de spectrométrie 2 est destiné à établir un spectre de rayonnement d'un échantillon 3, dans le domaine X et/ou Y. Plus précisément, le système 15 de spectrométrie 2 est configuré pour établir un spectre de rayonnement de l'échantillon 3, dans lequel le bruit correspondant à un fond de diffusion Compton a été rejeté, c'est-à-dire filtré. Le système de spectrométrie 2 est particulièrement destiné à l'établissement du spectre de rayonnement d'échantillons comprenant des isotopes 20 radioactifs d'activité faible. Dans ce cas, chacun parmi un signal de détection principal et un signal de détection auxiliaire, qui seront décrits ultérieurement, comporte des impulsions successives au cours du temps, majoritairement distinctes les unes des autres (régime dit de « temps mort faible »). Le système de spectrométrie 2 comporte un étage de détection 4 et un 25 étage de traitement 6 connecté en sortie de l'étage de détection 4. L'étage de détection 4 est configuré pour détecter des photons, et pour générer, en sortie, un signal électrique de détection principal et au moins un signal électrique de détection auxiliaire. En particulier, l'étage de détection 4 est configuré pour détecter des photons appartenant au domaine X et/ou Y. 12 Le signal de détection principal est représentatif de la détection de photons par un détecteur principal. Chaque signal de détection auxiliaire est représentatif de la détection de photons par un détecteur auxiliaire correspondant. Les détecteurs principal et auxiliaire(s) seront décrits ultérieurement. L'étage de traitement 6 est configuré pour traiter le signal de détection principal et chaque signal de détection auxiliaire, afin de déterminer le spectre de rayonnement de l'échantillon 3. Par « spectre de rayonnement », il est entendu, au sens de la présente invention, un vecteur dont chaque composante d'indice i (i entier naturel) est 10 représentative d'un nombre de photons émis par l'échantillon 3, par unité de temps, dans un intervalle d'énergie E associé à ladite composante d'indice i, lesdits photons ayant été générés par l'échantillon 3 et ayant interagi uniquement avec le détecteur principal et ayant déposé entièrement leur énergie dans celui-ci. Le spectre de rayonnement est également appelé « activité isotopique » de l'échantillon 3. 15 L'étage de détection 4 comprend le détecteur principal 8 et l'au moins un détecteur auxiliaire 10. L'étage de détection 4 comprend, en outre, un organe d'alimentation 12 pour l'alimentation du détecteur principal 8 et de chaque détecteur auxiliaire 10 en énergie électrique. Le détecteur auxiliaire 10 est disposé autour du détecteur principal 8. En 20 d'autres termes, le détecteur auxiliaire 10 recouvre au moins partiellement le détecteur principal 8. De préférence, le détecteur auxiliaire 10 recouvre le détecteur principal 8 de façon à intercepter un secteur angulaire valant au moins 2n sr (stéradian). Dans le cas au l'étage de détection 4 comporte une pluralité de détecteurs auxiliaires 10, les détecteurs auxiliaires 10 sont disposés autour du détecteur 25 principal 8, de façon à recouvrir au moins partiellement le détecteur principal 8. De préférence, les détecteurs auxiliaires 10 sont agencés de façon à ce qu'un secteur angulaire valant au moins 2n sr (stéradian) autour du détecteur principal 8 soit intercepté par les détecteurs auxiliaires 10. Le détecteur principal 8 est configuré pour détecter des photons dans une 30 plage d'énergie principale appartenant au domaine X et/ou Y, et pour générer un signal de 13 détection principal représentatif d'une telle détection. Plus précisément, le signal de détection principal comporte des impulsions, encore appelées « impulsions principales », chacune représentative d'une interaction du détecteur principal 8 avec un photon de la plage d'énergie principale. La plage d'énergie principale est continue ou discontinue. Le détecteur principal 8 est, par exemple, une diode germanium hyperpur (HPGe). Chaque détecteur auxiliaire 10 est configuré pour détecter des photons dans une plage d'énergie auxiliaire appartenant au domaine X et/ou Y, et pour générer un 10 signal de détection auxiliaire représentatif d'une telle détection. Plus précisément, chaque signal de détection auxiliaire comporte des impulsions, encore appelées « impulsions auxiliaire », chacune représentative d'une interaction du détecteur auxiliaire 10 correspondant avec un photon de la plage d'énergie auxiliaire associée. Chaque plage d'énergie auxiliaire est continue ou discontinue. 15 L'intersection entre la plage d'énergie principale et la plage d'énergie auxiliaire est non nulle. De préférence, le système de spectrométrie 2 est configuré pour établir le spectre de rayonnement de l'échantillon 3 dans la plage d'énergie correspondant à l'intersection de la plage d'énergie principale et de la plage d'énergie auxiliaire. Chaque détecteur auxiliaire 10 est, de préférence, un scintillateur, par 20 exemple un scintillateur au germanate de bismuth BGO, un scintillateur à l'iodure de sodium dopé au thallium Nal(TI), un scintillateur au bromure de lanthane dopé au cerium LaBr3(Ce), ou encore un scintillateur plastique, tel qu'un scintillateur commercialisé sous la référence « E.1-200 » par la société Eljen TechnologyTM. Avantageusement, l'étage de détection 4 comporte, en outre, un 25 blindage 14, configuré pour absorber (voire réfléchir) des particules dont l'énergie appartient à la plage d'énergie principale et/ou à la plage d'énergie auxiliaire. Le blindage 14 est disposé autour du (des) détecteur(s) auxiliaire(s) 10. En d'autres termes, le blindage 14 recouvre, de préférence intégralement, le(s) détecteur(s) auxiliaire(s) 10. Il en résulte que le blindage 14 recouvre le détecteur principal 14 8, de préférence de façon à intercepter un secteur angulaire valant au moins 4n sr autour du détecteur principal 8. Le blindage 14 est, par exemple, réalisé en plomb ou en tungstène, ou en tout autre matériau présentant une haute densité et/ou un numéro atomique élevé. Un tel blindage 14 est destiné à réduire la quantité de photons provenant de l'environnement (c'est-à-dire ne provenant pas de l'échantillon 3) qui parviennent au détecteur principal 8 et au(x) détecteur(s) auxiliaire(s) 10. Comme illustré par la figure 1, un volume 16 est ménagé entre le(s) détecteur(s) auxiliaire(s) 10 et le blindage 14. Le volume 16 est destiné à accueillir 10 l'échantillon 3. Par exemple, l'échantillon 3 est un échantillon solide. Selon un autre exemple, l'échantillon 3 est un fluide circulant dans une tuyauterie dont au moins une section est disposée dans le volume 14. Est maintenant décrit l'étage de traitement 6 dans le cas où l'étage de détection 4 comporte un unique détecteur auxiliaire 10. L'architecture de l'étage de 15 traitement 6 dans le cas où l'étage de détection 4 comporte plus d'un détecteur auxiliaire 10 s'en déduit directement. L'étage de traitement 6 comporte un convertisseur analogique-numérique 17, une mémoire 18 et un organe de calcul 20. Le convertisseur analogique- numérique 17 et l'organe de calcul 20 sont tous deux connectés à la mémoire 18. 20 Le convertisseur analogique numérique 17 est configuré pour acquérir, c'est-à-dire numériser, chacun du signal de détection principal et du signal de détection auxiliaire en provenance de l'étage de détection 4, en un signal de détection principal acquis et un signal de détection auxiliaire acquis, respectivement. La mémoire 18 est configurée pour stocker chacun du signal de détection 25 principal acquis et du signal de détection auxiliaire acquis. La mémoire 18 est, en outre, configurée pour stocker un modèle numérique du fonctionnement de l'étage de détection 4. Un tel modèle numérique est, par exemple, obtenu par en mettant en oeuvre une simulation numérique du transport de particules ionisantes par méthode de 30 Monte-Carlo (par exemple, au moyen des codes MCNP, TRIPOLI, GEANT, etc.). De 15 préférence, le modèle numérique ainsi obtenu est validé expérimentalement avec des sources-étalons. Un tel modèle comporte six matrices, respectivement notée Fi, ri, Fz, Wi, W'i et W2. Chacune des matrices Fi, ri, F2, Wi, W'i et W2 comporte jmax lignes et imax colonnes. Le terme général fi de la matrice Fl est égal à la probabilité qu'un photon incident d'énergie E et issu de l'échantillon 3 induise un dépôt d'énergie E1 dans le détecteur principal 8. Le terme général wij,i de la matrice Wi est égal à la probabilité qu'un 10 photon correspondant à un bruit de diffusion de type I.c, I.d, II.b, II.c et présentant une énergie E induise un dépôt d'énergie E1 dans le détecteur principal 8. Le terme général f2,0 de la matrice Fz est égal à la probabilité qu'un photon incident d'énergie E et issu de l'échantillon 3 induise un dépôt d'énergie E1 dans le détecteur auxiliaire 10. 15 Le terme général w2j,i de la matrice W2 est égal à la probabilité qu'un photon correspondant à un bruit de diffusion de type I.c, I.d, II.b, II.c et présentant une énergie E induise un dépôt d'énergie E1 dans le détecteur auxiliaire 10. Le terme général de la matrice ri est égal à la probabilité qu'un photon incident d'énergie E et issu de l'échantillon 3 induise un dépôt d'énergie E1 dans le 20 détecteur principal 8 sachant qu'auront été rejetés les événements où un dépôt d'énergie simultané non nul aura été détecté dans le détecteur auxiliaire 10. Par « impulsions simultanées », il est entendu, au sens de la présente invention, une impulsion principale et une impulsion auxiliaire qui sont séparées temporellement par une durée inférieure à une durée maximale prédéterminée, par 25 exemple inférieure à quelques microsecondes. Le terme général de la matrice W'i est égal à la probabilité qu'un photon correspondant à un bruit de diffusion de type I.c, I.d,11.b,11.c (c'est-à-dire provenant de l'environnement) et présentant une énergie E induise un dépôt d'énergie E1 dans le détecteur principal 8 sachant qu'auront été rejetés les événements où un dépôt d'énergie 30 simultané non nul aura été détecté dans le détecteur auxiliaire 10. 16 Il résulte de ce qui précède que l'indice in., est l'indice servant à repérer l'énergie maximale Eimax de l'échantillon 3, et imax l'indice servant à repérer l'énergie maximale Elmax déposée dans les détecteurs 8, 10. Par conséquent, l'inégalité suivante est vérifiée : 5 En outre, il ressort de ce qui précède que les matrices F1, ri, F2 sont des matrices de transfert du signal photonique issu de l'échantillon 3. En outre, les matrices Wi, W'i et W2 sont des matrices de transfert du signal photonique issu de l'environnement autour de l'étage de détection 4, assimilé à du bruit. En d'autres termes, les matrices Wi, et W2 correspondent chacune à un modèle du bruit, c'est-à-dire des photons ne 10 provenant pas de l'échantillon 3. Les bruits de diffusion de type I.a, I.b et II.a ne sont pas considérés car ils dépendent directement de l'échantillon 3. Le bruit de diffusion de type II.b est, quant à lui, pris en compte dans les matrices et W2 car, dans ce cas, le photon a perdu sa cohérence avec l'échantillon 3. L'organe de calcul 20 est configuré pour traiter le signal de détection 15 principal acquis et le signal de détection auxiliaire acquis stockés dans la mémoire 18. De préférence, l'organe de calcul 20 est configuré pour réaliser un filtrage du signal de détection principal acquis et du signal de détection auxiliaire acquis. Un tel filtrage consiste, par exemple, à affecter, dans chaque signal de détection acquis, une valeur nulle aux échantillons dont l'amplitude est inférieure ou égale à un seuil 20 prédéterminé. Ceci conduit à la génération d'un signal de détection principal filtré et un signal de détection auxiliaire filtré. En variante, ou en complément, l'organe de calcul 20 est configuré pour réaliser un filtrage du signal de détection principal acquis et du signal de détection auxiliaire acquis en affectant une valeur nulle aux impulsions dont l'amplitude est inférieure ou égale 25 à une amplitude minimale prédéterminée. L'organe de calcul 20 est, en outre, configuré pour détecter des impulsions coïncidentes dans le signal de détection principal acquis et le signal de détection auxiliaire acquis (ou le cas échéant, dans le signal de détection principal filtré et le signal de détection auxiliaire filtré), c'est-à-dire des impulsions qui, dans le signal de détection 30 principal acquis et le signal de détection auxiliaire acquis (ou, le cas échéant, dans le signal 17 de détection principal filtré et le signal de détection auxiliaire filtré), surviennent à une même date. L'organe de calcul 20 est également configuré pour construire un signal hybride, correspondant au signal de détection principal acquis (ou, le cas échéant, au signal de détection principal filtré) dont ont été supprimées les impulsions coïncidant avec des impulsions du signal de détection auxiliaire acquis (ou, le cas échéant, du signal de détection auxiliaire filtré). L'organe de calcul 20 est également configuré pour déterminer, pour chaque impulsion du signal de détection principal, du signal hybride et du signal de 10 détection auxiliaire, l'énergie qui a été déposée dans le détecteur 8, 10 correspondant par le photon ayant engendré ladite impulsion. Par exemple, pour une impulsion donnée, l'énergie déposée dans le détecteur 8, 10 par le photon correspondant est obtenue à partir de l'amplitude maximale de l'impulsion, ou encore à partir de l'aire de l'impulsion. L'organe de calcul 20 est, en outre, configuré pour calculer chacun parmi : 15 - un spectre Si du signal de détection principal acquis (ou, le cas échéant, du signal de détection principal filtré) ; - un spectre S'i du signal hybride ; et - un spectre Sz du signal de détection auxiliaire acquis (ou, le cas échéant, du signal de détection auxiliaire filtré). 20 Chaque spectre Si, S'i et 52 est un vecteur à jrnax lignes et une colonne. Chaque spectre Si, S'i et Sz est la réponse de l'étage de détection 4 aux radiations issues de l'échantillon 3 et de l'environnement, de sorte que, par exemple, Si vérifie la relation (1) suivante : 51 = FlA W113 (1) 25 où A, appelé « vecteur isotopique », est un vecteur à in,ax lignes et une colonne, dont chaque composante d'indice i est représentative d'un nombre de photons émis par l'échantillon 3, par unité de temps, dans un intervalle d'énergie E associé à ladite composante d'indice i, lesdits photons ayant interagi exclusivement avec le détecteur principal 8 et ayant déposé entièrement leur énergie dans celui-ci ; et 18 B, appelé « flux de l'ambiance radioactive », ou encore « flux de bruit », est un vecteur à imax lignes et une colonne, dont chaque composante d'indice i est représentative d'un nombre de photons, par unité de temps, provenant de l'environnement autour de l'étage de détection 4, dans un intervalle d'énergie E, associé à 5 ladite composante d'indice i, et incidents sur l'étage de détection 4. La relation (1) s'écrit également : {511 = fia,iAl + ... + fiilAi + ... + fi,i,'',tAimax + w11,1131 + - -- + w 1,i, 1 Bi + ... + w1,in,a',113imax Sm = ft,i,iAt + ... + fi,i,jAi + - + fi,i,'axjAimax + (2) wt,i,iBt + ... + wt,i,iBi + ... + wLimax,iBimax Stimax = fm,imax Al + - + fi,i,imax Ai + ... + (1,i,'a',imaxAimax + wuj,'ax Bi ± ..- + wi,i,imax Bi + ... + wumax,;,,,a,,Bi,',' De la même manière, S'I. et 52 vérifient les relations (3) et (4) suivantes : S'i = F'1A + W'1B (3) S2 = F2A + W2B (4) 10 L'organe de calcul 20 est également configuré pour calculer une estimée E). du flux de bruit B (ou « estimée du flux de bruit É1. »). Plus précisément, le système d'équations (1), (3) et (4) étant surdéterminé, l'étage de traitement 6 est configuré pour calculer l'estimée du flux de bruit .1; selon : = [VV2- l(Sz (S) 15 ou encore selon : B [w2 - 1(52- F2F11.1511) (6) ou encore selon : = P [W2- F2FilW1i 1(52 - (1-P)[W2- F2Frille1j l(Sz- F2F11Sr1) (7) avec p un réel tel que prédéterminé ou non. 20 Le symbole « » représente l'opérateur « inverse », l'organe de calcul 20 étant configuré pour mettre en oeuvre l'une quelconque des méthodes numériques ou algébriques d'inversion de matrice connues de l'homme du métier. En variante, pour réaliser une telle inversion, l'organe de calcul 20 est configuré pour mettre en oeuvre toute autre méthode n'impliquant par une inversion de matrice stricto sensu, mais plutôt une 19 estimation de la meilleure fonction de transfert entre un terme de droite et un terme de gauche sous une forme déterminée (par exemple Wz-F2F11W1 dans l'équation (5)). Ceci permet, notamment, de tenir compte, dans la recherche des estimées de la fonction de transfert, des incertitudes inhérentes à la nature stochastique des objets mathématiques que sont les matrices de transfert Fi, ri, F2 et W2. Une telle méthode est, par exemple, une méthode de type « Maximum Likelihood Expectation-Maximization » (pour méthode espérance-maximisation du maximum de vraisemblance) dédiées aux signaux nucléaires. L'étage de traitement 6 est également configuré pour calculer une 10 estimée Â de l'activité isotopique A (ou « estimée de l'activité isotopique Â »). Plus précisément, l'étage de traitement 6 est configuré pour calculer l'estimée de l'activité isotopique Â, en particulier si la relation (5) a été utilisée auparavant, selon : Â = W113) (8) ou encore, en particulier si la relation (6) a été utilisée auparavant, selon : 15 Â IN'e (9) ou encore selon : Â = W16) + (1-0F'11(5'1- W'11.3) (10) avec r un réel tel que prédéterminé ou non. Toutefois, le recours aux formules (6) et (9) est préféré, ces dernières 20 fournissant a priori une meilleure estimée de l'activité isotopique Â grâce au recours à Sit car, dans S'i, le bruit de diffusion de type I a été préalablement filtré. L'estimée de l'activité isotopique Â constitue une estimée du spectre d'émission de l'échantillon 3. Le fonctionnement du système de détection 2 va maintenant être décrit, 25 dans le cas où le système de détection 2 comporte un détecteur principal 8 et un détecteur auxiliaire 10. Au cours d'une phase d'initialisation, un modèle numérique du fonctionnement de l'étage de détection 4 est établi et stocké dans la mémoire 18 de l'étage de traitement 6. En particulier, les matrices de transfert F1, FIL F2, W1, W'i et W2 sont 30 stockées dans la mémoire 18. 20 Puis, le système de détection 2, et en particulier l'étage de détection 4, est exposé à un échantillon 3. Le détecteur principal 8 et le détecteur auxiliaire 10, qui sont alimentés en énergie électrique par l'organe d'alimentation 12, délivrent en sortie, un signal de détection principal et un signal de détection auxiliaire respectivement. Puis, le convertisseur analogique-numérique 17 de l'étage de traitement 6 acquiert le signal de détection principal et le signal de détection auxiliaire reçus depuis l'étage de détection 4, en un signal de détection principal acquis et un signal de détection auxiliaire acquis, respectivement. 10 Le signal de détection principal acquis et le signal de détection auxiliaire acquis sont stockés dans la mémoire 18. Puis, de préférence, l'organe de calcul 20 de l'étage de traitement 6 filtre le signal de détection principal acquis et le signal de détection auxiliaire acquis en un signal de détection principal filtré et un signal de détection auxiliaire filtré. 15 Puis, l'organe de calcul 20 détecte les impulsions coïncidentes dans le signal de détection principal acquis et le signal de détection auxiliaire acquis (ou, le cas échéant, dans le signal de détection principal filtré et le signal de détection auxiliaire filtré). L'organe de calcul 20 construit, en outre, le signal hybride, correspondant au signal de détection principal acquis dont ont été supprimées les impulsions coïncidant 20 avec des impulsions du signal de détection auxiliaire acquis (ou, le cas échéant, correspondant au signal de détection principal filtré dont ont été supprimées les impulsions coïncidant avec des impulsions du signal de détection auxiliaire filtré). Puis, pour chaque impulsion du signal de détection principal acquis et du signal de détection auxiliaire acquis (ou, le cas échéant, du signal de détection principal 25 filtré et du signal de détection auxiliaire filtré), et du signal hybride, l'organe de calcul 20 calcule l'énergie qui a été déposée dans le détecteur 8, 10 correspondant, par le photon ayant engendré ladite impulsion. Puis, l'organe de calcul 20 calcule chacun parmi : - le spectre Si du signal de détection principal acquis ; 30 - le spectre du signal hybride ; et 21 - le spectre S2 du signal de détection auxiliaire acquis. Puis, l'organe de calcul 20 calcule l'estimée du flux de bruit 13, en mettant en oeuvre la formule (5), (6) ou (7). Puis, l'organe de calcul 20 calcule l'estimée de l'activité isotopique Â, en 5 mettant en oeuvre la formule (8), (9) ou (10). De préférence, pour le calcul de l'estimée de l'activité isotopique Â, la formule (8) est utilisée si l'estimée du flux de bruit â est préalablement calculée au moyen de la formule (5), et la formule (9) est utilisée si l'estimée du flux de bruit â est préalablement calculée au moyen de la formule (6). Un deuxième mode de réalisation de l'invention diffère du système de 10 spectrométrie 2 décrit précédemment en ce que l'étage de détection 4 comporte no détecteurs auxiliaires (no étant un entier naturel supérieur ou égal à 2). Dans ce cas, l'indice e, compris entre 1 et no, permet d'établir une bijection entre un détecteur auxiliaire 10 et un secteur d'angles géométriques d'incidence correspondant, noté A6. Avantageusement, la somme de l'ensemble des secteurs d'angles géométriques d'incidence vaut 2n. 15 Dans ce cas, l'étage de détection 4 délivre no signaux de détection auxiliaires à l'étage de traitement 6. Dans ce cas, le modèle numérique du fonctionnement de l'étage de détection 4 comporte, outre les matrices de transfert Fi, F'1, F2, no matrices Wi,e, no matrices vtri,e et no matrices Wz,e. 20 Chacune des matrices Wi,e, W'i,e ou Wz,e est assimilable à une matrice de transfert ou Wz du signal photonique issu de l'environnement, pour le détecteur auxiliaire 10 d'indice 0 exclusivement. Dans ce cas, le convertisseur analogique-numérique 17 est configuré pour acquérir chaque signal de détection auxiliaire en un signal de détection auxiliaire acquis. 25 En outre, dans ce cas, l'organe de calcul 20 est configuré pour calculer, pour chaque signal de détection auxiliaire acquis, éventuellement filtré, issu d'un détecteur auxiliaire 10 d'indice 0, le spectre un spectre 52,0 correspondant. Chaque spectre Si, S'i et 52,0 est la réponse de l'étage de détection 4 aux radiations issues de l'échantillon 3 et de l'environnement. En particulier, les relations (1') 30 suivantes sont vérifiées : 22 Si = FIA + 4'21 Wi,e Be { s'i = FriA + 419_1 WiteBe S2,1 = F2A + W2,1Bi Be, étant un vecteur à imax lignes et une colonne, dont chaque composante d'indice i est représentative d'un nombre de photons, par unité de temps, provenant de l'environnement, dans un intervalle d'énergie Ei associé à ladite composante d'indice i, et 5 incidents sur l'étage de détection 4 selon le secteur d'angles géométriques d'incidence M. A titre d'exemple, et par analogie avec ce qui précède, la décomposition du spectre Si est donnée par les relations suivantes : 51,1 = --- --- fi,:'',1A:md' ne Xwt,t,i,eBi,e+---+wt,,,i,eB,,e+-wi,,,riax,t,eBimax,e 8=1 - = fiajAi fi,imexjAimax ne XW1,1,j,0 B1,0+---+Wl,i,j,0Bi,0+---Wl,iniax,10Bimax,el (2') 0=1 - Al+ + f A- + + f A + limax = umax max Lin,',jnia' max ne Xwt,ti,'',oBio+---+wi,,,i,,,, Bi +...1A/1,imaxiimay0 Biniaxm 0=1 Dans ce cas, l'organe de calcul est configuré pour calculer l'estimée de 10 l'activité isotopique Â d'une façon similaire à ce qui a été exposé précédemment, en mettant notamment en oeuvre les estimées du flux de bruit Be dans chaque secteur d'angles géométriques d'incidence M. En particulier, et comme cela ressort des relations (2'), l'ensemble des relations relatives aux spectres S2,1 à 52,n9 sont mises en oeuvre pour le calcul de l'estimée du flux de bruit â et de l'estimée de l'activité isotopique Â. 15 Selon une variante du premier ou du deuxième mode de réalisation, illustrée par la figure 2, le blindage 16 comporte une ouverture 22 formant collimateur. 5.2,0 = F2A WzoBe Szno= F2A W2,n0Bno 23 Dans ce cas, l'échantillon 3 est, de préférence, disposé d'un côté opposé du blindage 16 par rapport aux détecteur 8, 10, de préférence dans le champ de vue du collimateur. Alternativement, les matrices F11, F12-1, F2-1, W11, Wi-1 et W2-1 sont préalablement stockées dans la mémoire 18. Dans ce cas, l'organe de calcul 20 est 5 configuré pour déterminer l'estimée de l'activité isotopique P, et l'estimée du flux de bruit B par la mise en oeuvre directe des relations (5), (6), (7) et (8) sans que le calcul de l'inverse de chacune des matrices Fi, F'1, Fz, W1, Wi et W2 respectivement ne soit requis. Ceci s'applique également dans le cas où l'étage de détection 4 comporte no détecteurs auxiliaires. 10 Selon un mode de réalisation de l'invention, l'organe de calcul est configuré pour supprimer, dans chaque signal de détection auxiliaire, chaque impulsion telle que, à la même date, il existe une impulsion coïncidente dans le signal de détection issu du détecteur principal. Dans ce cas, les matrices F2 et W2 sont modifiées en conséquence, d'une façon similaire à celle permettant l'obtention REVENDICATIONS1. Système de spectrométrie (2) comprenant un étage de détection (4) comportant un détecteur principal (8) et au moins un détecteur auxiliaire (10), le détecteur principal (8) étant configuré pour détecter au moins des photons présentant une énergie appartenant à une plage d'énergie principale dans le domaine X et/ou Y, et pour délivrer un signal de détection principal comportant des impulsions principales, chaque impulsion principale étant représentative d'une interaction du détecteur principal (8) avec un photon de la plage d'énergie principale, chaque détecteur auxiliaire (10) étant configuré pour détecter au moins des photons présentant une énergie appartenant à une plage d'énergie auxiliaire dans le domaine X et/ou Y, et pour délivrer un signal de détection auxiliaire comportant des impulsions auxiliaires, chaque impulsion auxiliaire d'un signal de détection auxiliaire donné étant représentative d'une interaction du détecteur auxiliaire (10) correspondant avec un photon de la plage d'énergie auxiliaire associée, le système de spectrométrie (2) étant caractérisé en ce qu'il comporte, en outre, un étage de traitement (6) comprenant : - un organe d'acquisition (17) configuré pour acquérir le signal de détection principal en un signal de détection principal acquis, et chaque signal de détection auxiliaire en un signal de détection auxiliaire acquis correspondant ; - un organe de calcul (20) configuré pour : - détecter les impulsions principales dans le signal de détection principal acquis et les impulsions auxiliaires dans chaque signal de détection auxiliaire acquis ; - déterminer, pour chaque impulsion principale, une énergie principale correspondante et, pour chaque impulsion auxiliaire, une énergie auxiliaire correspondante ; 25 - calculer un premier spectre en fonction du signal de détection principal acquis et, éventuellement, de chaque signal de détection auxiliaire acquis ; - calculer au moins un deuxième spectre en fonction d'au moins un signal de détection auxiliaire acquis et, éventuellement, du signal de détection principal acquis ; - une mémoire (18) configurée pour stocker un modèle numérique associant une énergie donnée d'un photon incident parvenant à l'étage de détection (4), à: 10 - une énergie principale déposée par le photon incident dans le détecteur principal (8), dont la valeur dépend de ce que le photon incident est issu de la source à caractériser ou d'un environnement distinct de la source ; et - une énergie auxiliaire déposée par le photon incident 15 dans chaque détecteur auxiliaire (10), dont la valeur dépend de ce que le photon incident est issu de la source à caractériser ou de l'environnement ; l'organe de calcul (20) étant, en outre, configuré pour calculer un vecteur isotopique et/ou un flux de bruit comme étant chacun une solution d'un jeu d'équations 20 comprenant : - une première équation dans laquelle le premier spectre s'exprime comme la somme du produit du vecteur isotopique par un premier paramètre du modèle numérique d'une part, et du produit du flux de bruit par un deuxième paramètre du modèle numérique d'autre 25 part ; et - au moins une deuxième équation dans laquelle chaque deuxième spectre s'exprime comme la somme du produit du vecteur isotopique par un troisième paramètre du modèle numérique d'une part, et du produit du flux de bruit par un quatrième paramètre du modèle numérique d'autre part ; 26 le vecteur isotopique étant un vecteur dont chaque composante d'indice i est représentative d'un nombre de photons émis par la source à caractériser, par unité de temps, dans un intervalle d'énergie E, associé à ladite composante d'indice i, lesdits photons ayant interagi exclusivement avec le détecteur principal (8) et ayant déposé entièrement leur énergie dans le détecteur principal (8), le flux de bruit étant un vecteur dont chaque composante d'indice i est représentative d'un nombre de photons, par unité de temps, provenant de l'environnement, dans un intervalle d'énergie E, associé à ladite composante d'indice i, et incidents sur l'étage de détection (4). 10 2. Système selon la revendication 1, dans lequel l'équation relative au premier signal s'écrit : no FIA + Wi,eBe 0=1 où Si est le premier spectre, égal au spectre du signal de détection 15 principal acquis ; Fi est une matrice dont un terme général fi,;' est égal à la probabilité qu'un photon incident d'énergie E et issu de la source à caractériser induise un dépôt d'énergie Ej dans le détecteur principal ; A est le vecteur isotopique ; 20 no est le nombre de détecteurs auxiliaires, chaque détecteur auxiliaire d'indice 0 étant associé à un secteur d'angles géométriques d'incidence propre noté AS ; Wise est une matrice dont un terme général est égal à la probabilité qu'un photon provenant de l'environnement, incident sur l'étage de détection selon le secteur d'angles géométriques d'incidence AS, et présentant une énergie E, induise un 25 dépôt d'énergie E1 dans le détecteur principal ; et Be est la partie du flux de bruit correspondant à des photons incidents sur l'étage de détection selon le secteur d'angles géométriques d'incidence AS. 27 3. Système selon la revendication 1 ou 2, dans lequel le calculateur est configuré pour calculer un signal hybride, le signal hybride étant obtenu à partir du signal de détection principal, duquel a été supprimée chaque impulsion principale telle que, à la même date, il existe une impulsion auxiliaire coïncidente dans au moins un signal de détection auxiliaire acquis, un deuxième spectre parmi l'au moins un deuxième spectre étant le spectre du signal hybride. 4. Système selon la revendication 3, dans lequel une deuxième équation 10 parmi l'au moins une deuxième équation s'écrit : ne = Filé% WrteBe où est un deuxième spectre parmi l'au moins un deuxième spectre, égal au spectre du signal hybride ; ri est une matrice dont un terme général fil, est égal à la probabilité 15 qu'un photon incident d'énergie E, et issu de la source induise un dépôt d'énergie Ej dans le détecteur principal (8) sachant qu'auront été rejetées les impulsions principales telles que, à la même date, il existe au moins une impulsion auxiliaire coïncidente ; A est le vecteur isotopique ; no est le nombre de détecteurs auxiliaires, chaque détecteur auxiliaire 20 d'indice e étant associé à un secteur d'angles géométriques d'incidence propre noté A(3 ; est une matrice dont un terme général est égal à la probabilité qu'un photon provenant de l'environnement, incident sur l'étage de détection (4) selon le secteur d'angles géométriques d'incidence AS, et présentant une énergie E, induise un dépôt d'énergie E1 dans le détecteur principal (8) sachant qu'auront été rejetées 25 les impulsions principales telles que, à la même date, il existe au moins une impulsion auxiliaire coïncidente ; et Be est la partie du flux de bruit correspondant à des photons incidents sur l'étage de détection selon le secteur d'angles géométriques d'incidence M. 28 5. Système selon l'une quelconque des revendications 1 à 4, dans lequel une deuxième équation parmi l'au moins une deuxième équation s'écrit : Sze = F2A W2,0130 Où 0 est un indice associé à un détecteur auxiliaire donné, le détecteur 5 auxiliaire d'indice 0 étant associé à un secteur d'angles géométriques d'incidence propre noté M ; 52,e est le deuxième spectre, égal au spectre du signal auxiliaire acquis correspondant au détecteur auxiliaire d'indice 0 ; F2 est une matrice dont un terme général f2,j,, est égal à la probabilité 10 qu'un photon incident d'énergie E, et issu de la source à caractériser induise un dépôt d'énergie E1 dans le détecteur auxiliaire d'indice 6 ; A est le vecteur isotopique ; Wze est une matrice dont un terme général w2,j,,,o est égal à la probabilité qu'un photon provenant de l'environnement, incident sur l'étage de détection selon le 15 secteur d'angles géométriques d'incidence M, et présentant une énergie E, induise un dépôt d'énergie E1 dans le détecteur auxiliaire d'indice 0 ; et Be est la partie du flux de bruit correspondant à des photons incidents sur l'étage de détection selon le secteur d'angles géométriques d'incidence MI. 20 6. Procédé de spectrométrie comprenant les étapes suivantes : - détection d'impulsions principales dans un signal de détection principal acquis, le signal de détection principal acquis étant le résultat de l'acquisition d'un signal de détection principal délivré par un détecteur principal (8), le détecteur principal (8) étant configuré pour détecter au moins des photons présentant une énergie appartenant à une 25 plage d'énergie principale dans le domaine X et/ou Y, et pour délivrer le signal de détection principal, le signal de détection principal comportant les impulsions principales, chaque impulsion principale étant représentative d'une interaction du détecteur principal (8) avec un photon de la plage d'énergie principale ; - détection d'impulsions auxiliaires dans au moins un signal de détection 30 auxiliaire acquis, chaque signal de détection auxiliaire acquis étant le résultat d'un signal 29 de détection auxiliaire correspondant délivré par un détecteur auxiliaire (10) correspondant, chaque détecteur auxiliaire (10) étant configuré pour détecter au moins des photons présentant une énergie appartenant à une plage d'énergie auxiliaire dans le domaine X et/ou Y, et pour délivrer le signal de détection auxiliaire correspondant, chaque signal de détection auxiliaire comportant des impulsions auxiliaires, chaque impulsion auxiliaire d'un signal de détection auxiliaire donné étant représentative d'une interaction du détecteur auxiliaire (10) correspondant avec un photon de la plage d'énergie auxiliaire associée ; - détermination, pour chaque impulsion principale, d'une énergie 10 principale correspondante et, pour chaque impulsion auxiliaire, d'une énergie auxiliaire correspondante ; - calcul d'un premier spectre en fonction du signal de détection principal acquis et, éventuellement, de chaque signal de détection auxiliaire acquis ; - calcul d'au moins un deuxième spectre en fonction d'au moins un signal 15 de détection auxiliaire acquis et, éventuellement, du signal de détection principal acquis ; - calculer d'un vecteur isotopique et d'un flux de bruit comme étant les solutions d'un jeu d'équations comprenant : - une première équation dans laquelle le premier spectre s'exprime comme la somme du produit du vecteur isotopique par 20 un premier paramètre d'un modèle numérique prédéterminé d'une part, et du produit du flux de bruit par un deuxième paramètre du modèle numérique d'autre part ; et - au moins une deuxième équation dans laquelle chaque deuxième spectre s'exprime comme la somme du produit du 25 vecteur isotopique par un troisième paramètre du modèle numérique d'une part, et du produit du flux de bruit par un quatrième paramètre du modèle numérique d'autre part ; le vecteur isotopique étant un vecteur dont chaque composante d'indice i est représentative d'un nombre de photons émis par la source à caractériser, par unité de 30 temps, dans un intervalle d'énergie E, associé à ladite composante d'indice i, lesdits 30 photons ayant interagi exclusivement avec le détecteur principal (8) et ayant déposé entièrement leur énergie dans le détecteur principal (8), le flux de bruit étant un vecteur dont chaque composante d'indice i est représentative d'un nombre de photons, par unité de temps, provenant de l'environnement, dans un intervalle d'énergie E: associé à ladite composante d'indice i, et incidents sur l'étage de détection (4), le modèle numérique associant une énergie donnée d'un photon incident à: - une énergie principale déposée par le photon incident 10 dans le détecteur principal (8), dont la valeur dépend de ce que le photon incident est issu d'une source (3) à caractériser ou d'un environnement distinct de la source (3); et - une énergie auxiliaire déposée par le photon incident dans chaque détecteur auxiliaire (10), dont la valeur dépend de ce que le 15 photon incident est issu de la source (3) à caractériser ou de l'environnement. 7. Produit programme d'ordinateur comprenant des instructions de code de programme qui, lorsqu'elles sont exécutées par un ordinateur, mettent en oeuvre le 20 procédé selon la revendication 6.