[go: up one dir, main page]

CN115793032B - Peak drift correction system, method, equipment and storage medium for proportional counter - Google Patents

Peak drift correction system, method, equipment and storage medium for proportional counter Download PDF

Info

Publication number
CN115793032B
CN115793032B CN202211585791.XA CN202211585791A CN115793032B CN 115793032 B CN115793032 B CN 115793032B CN 202211585791 A CN202211585791 A CN 202211585791A CN 115793032 B CN115793032 B CN 115793032B
Authority
CN
China
Prior art keywords
correction
spectrogram
convolution
proportional counter
peak drift
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211585791.XA
Other languages
Chinese (zh)
Other versions
CN115793032A (en
Inventor
王环
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Suzhou Wuxiang Scientific Instrument Co ltd
Original Assignee
Suzhou Wuxiang Scientific Instrument Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Suzhou Wuxiang Scientific Instrument Co ltd filed Critical Suzhou Wuxiang Scientific Instrument Co ltd
Priority to CN202211585791.XA priority Critical patent/CN115793032B/en
Publication of CN115793032A publication Critical patent/CN115793032A/en
Application granted granted Critical
Publication of CN115793032B publication Critical patent/CN115793032B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Measurement Of Radiation (AREA)

Abstract

The invention belongs to the technical field of proportional counter tube correction, and particularly relates to a peak drift correction system, a peak drift correction method, peak drift correction equipment and a peak drift correction storage medium for proportional counter tubes, wherein the peak drift correction system, the peak drift correction method, the peak drift correction equipment and the peak drift storage medium comprise a spectrogram convolution calculation unit and a convolution correction unit which is in communication connection with the spectrogram convolution calculation unit; the spectrogram convolution calculation unit is used for calculating a spectrogram convolution sum; the convolution correction unit is used for calculating correction factors after convolution and correction of the spectrogram, and correcting the spectrogram based on the correction factors. According to the method, the light intensity of the excitation source is not required to be adjusted, the amplification factor is not required to be changed, calculation is carried out on the convolution of the original spectrogram, then the correction is carried out on the convolution of the spectrogram and the discarding time, and then the signal gain correction is carried out, the correction factor is calculated, and the channel of the original spectrogram is translated based on the correction factor, so that the drift correction is completed.

Description

Peak drift correction system, method, equipment and storage medium for proportional counter
Technical Field
The invention belongs to the technical field of proportional counter tube correction, and particularly relates to a peak drift correction system, method and equipment of a proportional counter tube and a storage medium.
Background
The proportional counter tube is an X-ray detector based on the photoionization effect, can measure the intensity and energy of incident X-rays at the same time, is mainly used for energy calibration and intensity measurement of X-rays, and has wide application in various fields such as X-ray astronomy, X-ray spectrum analysis, plasma diagnosis, surface physics, stress analysis and the like, in particular to the XRF fields such as plating thickness measurement, noble metal component analysis and the like.
The proportional counter tube detector is an X-ray detector based on the photoionization effect, and can measure the intensity and energy of incident X-rays at the same time. The center is an anode filament, the shell of the cylindrical barrel is a cathode, and incident particles collide with gas atoms in the barrel to ionize the atoms and generate electrons and positive ions. Under the action of the electric field, the electrons move to the central anode wire, and a pulse electric signal with the amplitude proportional to photon energy is generated.
As shown in fig. 2, when the energy spectrum analysis is performed, the element channel of the conventional proportional counter detector is shifted along with the change of the counting rate, and the energy channel of the conventional proportional counter detector is shifted left along with the increase of the counting rate. This creates interference with the analysis of the spectra we follow. At present, the common practice is to adjust the light intensity of an excitation source to ensure the consistency of sample counting rates, or respectively calibrate, and change the amplification factor to keep the consistency of peak channels so as to avoid calculation deviation caused by peak drift. The common practice is complicated in actual operation, and the implementation conditions cannot be met in many application scenes, such as the conditions that an excitation source is not adjustable or the adjustment precision cannot be achieved.
Disclosure of Invention
In view of the above-mentioned shortcomings, the present invention aims to provide a peak drift correction system, method, device and storage medium for proportional counter.
The invention provides the following technical scheme:
a peak drift correction system of a proportional counter tube comprises a spectrogram convolution calculation unit and a convolution correction unit which is in communication connection with the spectrogram convolution calculation unit; the spectrogram convolution calculation unit is used for calculating a spectrogram convolution sum; the convolution correction unit is used for calculating correction factors after convolution and correction of the spectrogram, and correcting the spectrogram based on the correction factors.
The device also comprises a spectrogram generating unit; the spectrogram generating unit is used for generating an uncorrected spectrogram after the X-rays enter the proportional counter tube.
The spectrogram generating unit comprises a proportional counter tube, a pre-amplifying circuit, a signal amplifying and ADC circuit, a multichannel pulse amplitude analyzer and a spectrogram storage unit before correction, wherein the multichannel pulse amplitude analyzer comprises a digital forming filtering unit and a pulse screening logic unit. After entering the proportional counter tube, the X-ray photons are processed by a pre-amplifying circuit, each photon is converted into an electric signal pulse, the pulse amplitude is in direct proportion to the photon energy, the signals are processed by a signal amplifying and ADC circuit and then are converted into a spectrogram by a multi-channel pulse amplitude analyzer, the horizontal axis of the spectrogram is photon energy, the vertical axis of the spectrogram is photon number, and the spectrogram is stored in a spectrogram storage unit before correction.
The system also comprises a corrected spectrogram storage unit which is in communication connection with the convolution correction unit. The corrected spectrogram storage unit is used for storing the corrected spectrogram.
A peak drift correction method of a proportional counter tube is applied to a peak drift correction system of the proportional counter tube, and comprises the following steps:
s1, calculating a spectral graph convolution sum;
s2, carrying out discarding time correction on spectrogram convolution sum, and then carrying out signal gain correction;
s3, calculating a correction factor;
s4, correcting the spectrogram based on the correction factors.
In S4, the calculation formula of the spectrogram correction is: n (N) t1 =n×γ; wherein N is the number of channels before correction, N t1 For the corrected number of channels, γ is a correction factor.
The calculation formula of the correction factor is: gamma=g (N) t2 * h, performing H; wherein h is a correction constant, g (N) t2 The volume is the corrected spectral plot.
The calculation formula for the spectrogram convolution and the discard time correction is as follows: g (N) t1 =g (N)/readtime; wherein g (N) t1 For the compensated spectral graph convolution sum, deadtime is the dead time of the processor, g (N) is the spectral graph convolution sum;
readtime=ocr/ICR, OCR being the output effective count rate and ICR being the input count rate;
the calculation formula for the spectrogram convolution sum and the signal gain correction is as follows: g (N) t2 =g(N) t1 A0; wherein A0 is the total amplification factor of the multichannel pulse amplitude analyzer in the signal processing process.
The calculation formula of the volume of the spectrum chart is as follows:
where In is the photon count for that channel and Nmax is the maximum channel value.
The peak drift correction equipment for the proportional counter comprises a memory, a processor and a computer program which is stored in the memory and can run on the processor, wherein the processor realizes the peak drift correction method for the proportional counter when executing the program.
A computer readable storage medium having stored thereon a computer program which when executed by a processor implements the peak drift correction method of proportional counter as described above.
The beneficial effects of the invention are as follows: according to the method, the light intensity of the excitation source is not required to be adjusted, the amplification factor is not required to be changed, calculation is carried out on the convolution of the original spectrogram, then the correction is carried out on the convolution of the spectrogram and the discarding time, and then the signal gain correction is carried out, the correction factor is calculated, and the channel of the original spectrogram is translated based on the correction factor, so that the drift correction is completed. The method can be implemented in a digital multichannel pulse amplitude analyzer (in a lower computer FPGA or a DSP) or in upper computer software, the obtained final spectrogram is corrected, the correction precision is high, the influence of the counting rate on the peak channel can be eliminated, and the method is simple and effective on the application operation level.
Drawings
FIG. 1 is a schematic diagram of the system of the present application;
FIG. 2 is a schematic diagram before spectrogram correction;
fig. 3 is a schematic diagram after spectrogram correction.
Marked in the figure as: the proportional counter and pre-amplifier circuit 101, the signal amplifying and ADC circuit 102, the digital shaping filter unit 103, the pulse screening logic unit 104, the spectrogram before correction storage unit 105, the spectrogram convolution calculation unit 106, the convolution correction unit 107, and the spectrogram after correction storage unit 108.
Detailed Description
Example 1
As shown in fig. 1, the peak shift correction system of the proportional counter includes a spectrogram convolution calculating unit 106, a convolution correcting unit 107 communicatively connected to the spectrogram convolution calculating unit 106, a spectrogram generating unit communicatively connected to the spectrogram convolution calculating unit 106, and a corrected spectrogram storage unit 108 communicatively connected to the convolution correcting unit 107.
Specifically, the spectrogram generating unit is used for generating an uncorrected spectrogram after the X-rays enter the proportional counter tube. The spectrogram generating unit comprises a proportional counter tube and pre-amplifying circuit 101, a signal amplifying and ADC circuit 102, a multichannel pulse amplitude analyzer and a spectrogram storage unit before correction 105, wherein the multichannel pulse amplitude analyzer comprises a digital shaping filtering unit 103 and a pulse screening logic unit 104. After entering the proportional counter tube, the X-ray photons are processed by a pre-amplifying circuit, each photon is converted into an electric signal pulse, the pulse amplitude is in direct proportion to the photon energy, the signals are processed by the signal amplifying and ADC circuit 102 and then converted into a spectrogram by a multi-channel pulse amplitude analyzer, the horizontal axis of the spectrogram is photon energy, the vertical axis of the spectrogram is photon number, and the spectrogram is stored in the spectrogram storage unit 105 before correction.
The spectrogram convolution calculation unit 106 is configured to calculate a spectral image convolution sum. The convolution correction unit 107 is configured to calculate a correction factor after convolving and correcting the spectrogram, and correct the spectrogram based on the correction factor. The corrected spectrogram storage unit 108 is used for storing the corrected spectrogram.
The correction can be completed based on FPGA or DSP of the digital multichannel pulse amplitude analyzer, the digital multichannel pulse amplitude analyzer stores the spectrogram before correction in the RAM block, then completes calculation through logic or DSP of the digital multichannel pulse amplitude analyzer, stores the spectrogram after correction in the other RAM block, and waits for reading by an upper computer.
The correction can also be based on the upper computer software, the data of the digital multichannel pulse amplitude analyzer is extracted into the memory, and the correction is completed according to the same method.
Example two
The embodiment provides a peak drift correction method of a proportional counter, which is applied to a peak drift correction system of the proportional counter of the first embodiment, and comprises the following steps:
s1, calculating a spectral graph convolution sum;
s2, carrying out discarding time correction on spectrogram convolution sum, and then carrying out signal gain correction;
s3, calculating a correction factor;
s4, correcting the spectrogram based on the correction factors.
Specifically, in S1, the calculation formula of the spectral graph convolution sum is:
where In is the photon count for that channel and Nmax is the maximum channel value.
In S2, the calculation formula for carrying out the convolution sum of the spectrogram and the discarding time correction is as follows: g (N) t1 =g (N)/readtime; wherein g (N) t1 For the compensated spectral graph convolution sum, deadtime is the dead time of the processor and g (N) is the spectral graph convolution sum.
readtime=ocr/ICR, OCR being the output effective count rate and ICR being the input count rate;
the calculation formula for the spectrogram convolution sum and the signal gain correction is as follows: g (N) t2 =g(N) t1 A0; wherein A0 is the total amplification factor of the multichannel pulse amplitude analyzer in the signal processing process.
A0 contains all signal gain amplification from the signal output from the proportional counter to the spectrogram output. When the multichannel pulse amplitude analyzer adopts a digital multichannel pulse amplitude analyzer, the digital amplification is also in the calculation range.
In S3, the calculation formula of the correction factor is: gamma=g (N) t2 * h, performing H; wherein h is a correction constant, g (N) t2 The volume is the corrected spectral plot.
The correction constant h is measured through actual experimental data, in practical application, h can be used as a settable quantity, and the final corrected channel offset is ensured to be zero by adjusting the value of h, so that the optimal effect is realized.
In S4, the calculation formula of the spectrogram correction is: n (N) t1 =n×γ; wherein N is the number of channels before correction, N t1 For the corrected number of channels, γ is a correction factor.
As shown in FIG. 3, the count of channel N of the original spectrogram is translated to channel N t1 After the completion of the translation, the drift correction is completed.
Example III
The present embodiment provides a peak drift correction device for a proportional counter, including a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor implements the peak drift correction method for a proportional counter according to the second embodiment when executing the program.
Example IV
The present embodiment provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the peak shift correction method of the proportional counter as described in the second embodiment.
The foregoing description is only a preferred embodiment of the present invention, and the present invention is not limited thereto, but it is to be understood that modifications and equivalents of some of the technical features described in the foregoing embodiments may be made by those skilled in the art, although the present invention has been described in detail with reference to the foregoing embodiments. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (8)

1. A peak drift correction system of a proportional counter tube is characterized in that: the device comprises a spectrogram convolution calculation unit and a convolution correction unit which is in communication connection with the spectrogram convolution calculation unit; the spectrogram convolution calculation unit is used for calculating a spectrogram convolution sum, and a calculation formula of the spectrogram convolution sum is as follows:in is photon count of the energy channel, and Nmax is the maximum channel value; the convolution correction unit is used for calculating correction factors after carrying out discarding time correction and signal gain correction on spectrogram convolution sum, and correcting the spectrogram based on the correction factors;
the calculation formula for the spectrogram convolution and the discard time correction is as follows: g (N) t1 =g (N)/readtime; wherein g (N) t1 For the compensated spectral graph convolution sum, deadtime is the dead time of the processor, g (N) is the spectral graph convolution sum; readtime=ocr/ICR, OCR being the output effective count rate and ICR being the input count rate;
convolving the spectrogramThe calculation formula for signal gain correction is: g (N) t2 = g(N) t1 A0; wherein A0 is the total amplification factor of the multichannel pulse amplitude analyzer in the signal processing process.
2. The peak shift correction system of proportional counter according to claim 1, wherein: the device also comprises a spectrogram generating unit; the spectrogram generating unit is used for generating an uncorrected spectrogram after the X-rays enter the proportional counter tube.
3. The peak shift correction system of proportional counter according to claim 1, wherein: the system also comprises a corrected spectrogram storage unit; the corrected spectrogram storage unit is used for storing the corrected spectrogram.
4. A peak drift correction method for a proportional counter, applied to the peak drift correction system for a proportional counter according to claim 1, comprising the following steps:
s1, calculating a spectral graph convolution sum;
s2, carrying out discarding time correction on spectrogram convolution sum, and then carrying out signal gain correction;
s3, calculating a correction factor;
s4, correcting the spectrogram based on the correction factor;
the calculation formula of the volume of the spectrum chart is as follows:in is photon count of the energy channel, and Nmax is the maximum channel value;
the calculation formula for the spectrogram convolution and the discard time correction is as follows: g (N) t1 =g (N)/readtime; wherein g (N) t1 For the compensated spectral graph convolution sum, deadtime is the dead time of the processor, g (N) is the spectral graph convolution sum;
readtime=ocr/ICR, OCR being the output effective count rate and ICR being the input count rate;
the calculation formula for the spectrogram convolution sum and the signal gain correction is as follows: g (N) t2 = g(N) t1 A0; wherein A0 is the total amplification factor of the multichannel pulse amplitude analyzer in the signal processing process.
5. The method for correcting peak drift of proportional counter according to claim 4, wherein: in S4, the calculation formula of the spectrogram correction is: n (N) t1 =n×γ; wherein N is the number of channels before correction, N t1 For the corrected number of channels, γ is a correction factor.
6. The method for correcting peak drift of proportional counter according to claim 5, wherein: the calculation formula of the correction factor is: gamma=g (N) t2 * h, performing H; wherein h is a correction constant, g (N) t2 The volume is the corrected spectral plot.
7. A peak drift correction device for a proportional counter, comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, characterized by: the processor, when executing the program, implements the peak drift correction method for proportional counter according to any one of claims 4 to 6.
8. A computer-readable storage medium having stored thereon a computer program, characterized by: the program when executed by a processor implements the peak drift correction method for proportional counter according to any one of claims 4-6.
CN202211585791.XA 2022-12-09 2022-12-09 Peak drift correction system, method, equipment and storage medium for proportional counter Active CN115793032B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211585791.XA CN115793032B (en) 2022-12-09 2022-12-09 Peak drift correction system, method, equipment and storage medium for proportional counter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211585791.XA CN115793032B (en) 2022-12-09 2022-12-09 Peak drift correction system, method, equipment and storage medium for proportional counter

Publications (2)

Publication Number Publication Date
CN115793032A CN115793032A (en) 2023-03-14
CN115793032B true CN115793032B (en) 2024-01-30

Family

ID=85418497

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211585791.XA Active CN115793032B (en) 2022-12-09 2022-12-09 Peak drift correction system, method, equipment and storage medium for proportional counter

Country Status (1)

Country Link
CN (1) CN115793032B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63298087A (en) * 1987-05-29 1988-12-05 Nippon Atom Ind Group Co Ltd Correcting method for counting rate
CN108042148A (en) * 2017-11-30 2018-05-18 江苏赛诺格兰医疗科技有限公司 The method of the real time correction of pet detector spectrum drift and spectrum drift correction system
WO2019025479A1 (en) * 2017-08-02 2019-02-07 The European Atomic Energy Community (Euratom), Represented By The European Commission Method and system for pulse multiplicity counting with dead time correction
CN111025373A (en) * 2019-12-26 2020-04-17 中广核久源(成都)科技有限公司 Method for digitally correcting decay time of sodium iodide crystal in real time
JP2021148513A (en) * 2020-03-17 2021-09-27 福島SiC応用技研株式会社 Spectrum display device, spectrum display method, program for spectrum display device, and spectrum measurement system
CN114252901A (en) * 2021-11-18 2022-03-29 清华大学 Dead time correction method and system

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2953298B1 (en) * 2009-11-30 2014-10-31 Commissariat Energie Atomique METHOD OF CORRECTING THE STACK PHENOMENON APPLIED TO ACQUIRED X-RAY SPECTRA USING A SPECTROMETRIC SENSOR
US8937275B2 (en) * 2012-10-12 2015-01-20 Thermo Fisher Scientific Inc. Method and apparatus to monitor gain of a proportional counter including correcting the counting threshold of a pulse height spectrum
US8927920B2 (en) * 2013-02-07 2015-01-06 Schlumberger Technology Corporation Correcting gamma-ray energy spectra for pileup degradation
GB2514418A (en) * 2013-05-24 2014-11-26 Oxford Instr Nanotechnology Tools Ltd Method of processing a particle spectrum
FR3069066B1 (en) * 2017-07-17 2019-08-16 Commissariat A L'energie Atomique Et Aux Energies Alternatives METHOD FOR TREATING AN IMPULSE GENERATED BY A IONIZING RADIATION DETECTOR
JP6808693B2 (en) * 2018-08-16 2021-01-06 日本電子株式会社 X-ray analyzer and counting rate correction method
US11344266B2 (en) * 2019-09-16 2022-05-31 Redlen Technologies, Inc. Calibration methods for improving uniformity in X-ray photon counting detectors

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS63298087A (en) * 1987-05-29 1988-12-05 Nippon Atom Ind Group Co Ltd Correcting method for counting rate
WO2019025479A1 (en) * 2017-08-02 2019-02-07 The European Atomic Energy Community (Euratom), Represented By The European Commission Method and system for pulse multiplicity counting with dead time correction
CN108042148A (en) * 2017-11-30 2018-05-18 江苏赛诺格兰医疗科技有限公司 The method of the real time correction of pet detector spectrum drift and spectrum drift correction system
CN111025373A (en) * 2019-12-26 2020-04-17 中广核久源(成都)科技有限公司 Method for digitally correcting decay time of sodium iodide crystal in real time
JP2021148513A (en) * 2020-03-17 2021-09-27 福島SiC応用技研株式会社 Spectrum display device, spectrum display method, program for spectrum display device, and spectrum measurement system
CN114252901A (en) * 2021-11-18 2022-03-29 清华大学 Dead time correction method and system

Also Published As

Publication number Publication date
CN115793032A (en) 2023-03-14

Similar Documents

Publication Publication Date Title
Day et al. Photoelectric quantum efficiencies and filter window absorption coefficients from 20 eV to 10 keV
US5463219A (en) Mass spectrometer system and method using simultaneous mode detector and signal region flags
US5677536A (en) Gamma camera with on the fly calibration for PMT drift
WO2010150301A1 (en) Mass spectrometer
CN109799541B (en) A Spectral Drift and Integral Nonlinear Correction Method for Gamma Spectrum Measurement
JP2013512443A5 (en)
CN104898159B (en) A kind of calibration method of energy resolution CT detectors
CN107390253B (en) Single photon emission tomographic imaging equipment Quality Control scaling method and system
CN115793032B (en) Peak drift correction system, method, equipment and storage medium for proportional counter
JP5945245B2 (en) Signal pulse detection device, mass spectrometer, and signal pulse detection method
WO2001031552A1 (en) Differential correction method and apparatus
CN114740515B (en) A method for correcting trajectory loss with amplitude ratio
Carramate et al. THCOBRA X-ray imaging detector operating in pure Kr
CN206725471U (en) A kind of quantitatively calibrating system of X-ray imaging plate
CN106950233B (en) Quantitative calibration system of X-ray imaging plate
US7312457B2 (en) Spectrometer with charge-carrier-trapping correction
CN114089409B (en) Detector correction method and system
JP2921854B2 (en) Method for eliminating count loss effect of X-ray counter and linear amplifier
CN112946726B (en) Nonlinear correction method and system based on output signal of FPGA-ADC module
US10510522B2 (en) Mass spectrometry using plasma ion source
CN114779313A (en) Ray energy spectrum monitoring system and energy spectrum correction method
Blaj et al. Robust pixel gain calibration with limited statistics
Moran et al. SD-WDS: Bremsstrahlung Determination by Theoretical Calculation and Subsequent Estimation of Mass Attenuation Coefficients
MacGregor et al. The third forbidden beta spectrum of rubidium-87
JPH0262993A (en) Position detection type radiation detecting device and its adjusting method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant