CN107272049A - Digital n γ discriminating methods based on pulse width - Google Patents
Digital n γ discriminating methods based on pulse width Download PDFInfo
- Publication number
- CN107272049A CN107272049A CN201710569098.6A CN201710569098A CN107272049A CN 107272049 A CN107272049 A CN 107272049A CN 201710569098 A CN201710569098 A CN 201710569098A CN 107272049 A CN107272049 A CN 107272049A
- Authority
- CN
- China
- Prior art keywords
- digital waveform
- pulse
- pulse width
- amplitude
- digital
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T5/00—Recording of movements or tracks of particles; Processing or analysis of such tracks
- G01T5/02—Processing of tracks; Analysis of tracks
Abstract
The invention discloses a kind of digital n γ discriminating methods based on pulse width, including：Carried out in mixed radiation field using detector and capture card from trigger-type digital waveform SnSample and store；Will sampling digital waveform SnCubic spline interpolation processing is carried out, the digital waveform S of 4 sampling points is obtained4n；Digital waveform S4nHandled through First derivative spectrograply, obtain digital waveform S4nPeak point P and peak value Vp；Calculate maximum white noise VnWith activation threshold value VTRatio R, and also be pulse terminal amplitude VsWith peak value VpRatio；It is counter to push away acquisition terminal amplitude V using ratio Rs, and determine the starting point K and cut off Q of pulse；Width between starting point K and cut off Q is pulse width W, and obtaining n γ using pulse width W screens result.The discriminating method has the advantages that to calculate easy, exclusion subjective parameters setting influence, screens accurate, has very high practical value and promotional value in radiation detection technology field.
Description
Technical field
The present invention relates to radiation detection technology field, and in particular to the digital n- γ discriminating methods based on pulse width.
Background technology
Since neutron is found, it has been widely used in radiation detection field.In mixed radiation field, gamma-rays and neutron
N be it is most common in radiation field be also most important two kinds of rays, because two kinds of rays are not charged, neutron occasion similarly
With the presence of gamma-rays, neutron detector is also sensitive to gamma-rays, and gamma ray detector can be to neutron-sensitive.Therefore, n- γ are screened
Technology is to improve single radionetric survey precision in mixed radiation field, reduces one of the key of neutron (or gamma-rays) background.
It is when producing different charged particles in detector using neutron and gamma-rays to cause different decay that n- γ, which are screened,
Between speed component ratio, so as to the difference of impulse waveform occur.At present, most common digital n- γ discriminating methods be charge ratio compared with
Method, screens incoming particle, this method is not by fast component (or slow component) electric charge of impulse waveform and the proportional difference of total electrical charge
Foot part be, it is necessary to specify 3 parameter calculated charge ratios, i.e., total electrical charge integration starting point, total electrical charge integration cut off and
Fast component charge integration cut off (or slow composition charge integration starting point), these three parameters directly determine the quality that n- γ are screened
The factor.Because the specified not strict normative reference of these three parameters is, it is necessary to according to many suboptimization determinations of quality factor, lead
Cause calculating process cumbersome, and the quality factor that n- γ are screened is larger by subjective impact.
Patent publication No. for CN201620788584.8 Chinese patent disclose it is a kind of digitize neutron spectrometer neutron,
γ discrimination systems, it uses pulse fall time to screen algorithm, and pulse signal peak atenuation is discriminated to different proportion time difference Δ t
Other neutron and gamma-rays, by setting maximum X% and maximum Y% time differences one-to-one corresponding mapping relations, it artificially sets
Fixed maximum X% and maximum Y% is that changing value could obtain optimal value, it is necessary to carry out repeatedly calculation optimization.This method
The calculating process for the quality factor that optimization n- γ are screened is simplified to a certain extent, but two parameters set are without strict ginseng
Standard is examined, still needs optimization to calculate, the problem of not solving to calculate cumbersome.
To sum up, it is necessary to n- γ screening techniques are proposed to improve, to obtain the discriminating method of simplicity.
The content of the invention
It is an object of the invention to provide a kind of digital n- γ discriminating methods based on pulse width, existing skill is mainly solved
Algorithm is screened present in art specifies parameter more, calculates complex, quality factor optimization process is larger by subjective impact
Problem.
To achieve these goals, the technical solution adopted by the present invention is as follows：
Digital n- γ discriminating methods based on pulse width, it mainly result in using the die-away time difference of different particles
Corresponding pulses decay to time difference required during identical amplitude proportion, and it is wide that this difference can intuitively be embodied in pulse
On degree.Therefore n- γ examinations can effectively be carried out according to pulse width, comprised the following steps that：
The first step：In mixed radiation field, carried out using detector and capture card from trigger-type digital waveform SnSampling is simultaneously
Storage, digital waveform SnFor activation threshold value VTPreceding m point and activation threshold value VTN-m point afterwards, wherein, m is more than 0, n and is more than m.
Second step：Will sampling digital waveform SnCubic spline interpolation processing is carried out, the digital waveform of 4 sampling points is obtained
S4n。
3rd step：Digital waveform S4nHandled through First derivative spectrograply, obtain digital waveform S4nPeak point P and peak value Vp。
4th step：By digital waveform S4nThe maximum amplitude of preceding 2m point be used as maximum white noise Vn。
5th step：Calculate maximum white noise VnWith activation threshold value VTRatio R, and R also be pulse terminal amplitude VsWith peak
Value VpRatio.
6th step：It is counter to push away acquisition pulse terminal amplitude V using ratio Rs, and determine the starting point K and cut off of pulse
Q, is comprised the following steps that：
(1) from digital waveform S4nPeak point P pointwises be moved to the left, and successively with pulse terminal amplitude VsCompare, if
Starting point amplitude VK≤ pulse terminal amplitude VsWhen, then K points are starting point, otherwise continue to compare to the left；
(2) from digital waveform S4nPeak point P pointwises move right, and successively with pulse terminal amplitude VsCompare, if
Cut off amplitude VQ≤ pulse terminal amplitude VsWhen, then Q points are cut off, otherwise continue to compare to the right.
Wherein, the width between starting point K and cut off Q is pulse width W.
Finally, it is distributed using the pulse width W, obtains n- γ and screen result.
Compared with prior art, the invention has the advantages that：
(1) present invention enters the digit pulse ripple of generation by way of detector and capture card trigger digit pulse ripple
The processing of row cubic spline interpolation, only sets a parameter on waveform, i.e., maximum white noise VnWith activation threshold value VTRatio R, its
In, activation threshold value VTDetermine parameter when being used for impulse waveform, ratio R only with maximum white noise VnIt is relevant, arbitrfary point amplitude and arteries and veins
Rush terminal amplitude VsContrast, obtains pulse width W, using pulse width to realize that n- γ are screened, simplifies cumbersome parameter
The calculating process of optimization.
(2) moreover, pulse terminal amplitude VsWith peak VpValue than R only by sample waveform maximum white noise and touch
Send out threshold value to determine, it is poor to be prevented effectively from the quality factor caused by subjective factor, n- γ examinations are impacted.
(3) present invention is using easy discriminating method, the result that it is screened compared with existing ripe Charge comparison method,
It can more differentiate the peak value of two kinds of particles, and screen that process is more simple, quality factor are not influenceed by subjective setup parameter, the party
Method has the advantages that to calculate easy, exclusion subjective parameters setting influence, screens accurate, has in radiation detection technology field very high
Practical value and promotional value.
Brief description of the drawings
Fig. 1 is the schematic diagram of the digital n- γ discriminating methods based on pulse width.
Fig. 2 is the flow chart of the digital n- γ discriminating methods based on pulse width.
Fig. 3 is the hangover integral charge and total mark charge ratio distribution map of Charge comparison method.
Fig. 4 is the 2D distribution maps of the charge ratio and pulse amplitude obtained using Charge comparison method
Fig. 5 is the pulse width distribution map of the digital n- γ discriminating methods based on pulse width.
Fig. 6 is the pulse width and the 2D distribution maps of pulse amplitude of the digital n- γ discriminating methods based on pulse width.
Embodiment
The invention will be further described with reference to the accompanying drawings and examples, and embodiments of the present invention include but is not limited to
The following example.
Embodiment
As shown in figure 1, neutron pulse waveform S1 and gamma-ray burst waveform S2 decay to it is corresponding during identical amplitude proportion
Time point is different, i.e. time t' and time t, the difference of different particle die-away times are intuitively embodied on pulse width W, according to
Pulse width W difference can effectively carry out n- γ examinations.
As shown in Fig. 2 the digital n- γ discriminating methods based on pulse width, discriminating method comprises the following steps：
The first step：In mixed radiation field, carried out using detector and capture card from trigger-type digital waveform SnSampling is simultaneously
Storage, digital waveform SnFor activation threshold value VTPreceding m point and activation threshold value VTN-m point afterwards, wherein, m is more than 0, n and is more than m.
Second step：Will sampling digital waveform SnCubic spline interpolation processing is carried out, the digital waveform of 4 sampling points is obtained
S4n。
3rd step：Digital waveform S4nHandled through First derivative spectrograply, obtain digital waveform S4nPeak point P and peak value Vp。
4th step：By digital waveform S4nThe maximum amplitude of preceding 2m point be used as maximum white noise Vn。
5th step：Calculate maximum white noise VnWith activation threshold value VTRatio R, and R also be pulse terminal amplitude VsWith peak
Value VpRatio.
6th step：It is counter to push away acquisition pulse terminal amplitude V using ratio Rs, and determine the starting point K and cut off of pulse
Q, is comprised the following steps that：
(1) from digital waveform S4nPeak point P pointwises be moved to the left, and successively with pulse terminal amplitude VsCompare, if
Starting point amplitude VK≤ pulse terminal amplitude VsWhen, then K points are starting point, otherwise continue to compare to the left；
(2) from digital waveform S4nPeak point P pointwises move right, and successively with pulse terminal amplitude VsCompare, if
Cut off amplitude VQ≤ pulse terminal amplitude VsWhen, then Q points are cut off, otherwise continue to compare to the right.
Wherein, the width between starting point K and cut off Q is pulse width W.
Finally, it is distributed using the pulse width W, obtains n- γ and screen result.
Experimental data use EJ299-33A detectors electron energy on 7MeV phot-neutron source by 2.5Gs/s's
Digital oscilloscope sampling is measured.Wherein, Fig. 3 is the hangover integral charge and total mark charge ratio distribution map, Fig. 5 of Charge comparison method
For the pulse width distribution map of the present invention, Fig. 4 is the 2D distribution maps of the charge ratio and pulse amplitude obtained using Charge comparison method,
Fig. 6 is the 2D distribution maps using pulse width and pulse amplitude of the invention, and the inventive method is can be seen that more from comparison diagram
The peak value of two kinds of particles can be differentiated, n- γ examinations quality factor are higher, and examination process is more easy, and quality factor are not by subjectivity
Setup parameter influences.
Above-described embodiment is only the preferred embodiments of the present invention, not limiting the scope of the invention, as long as using
The design principle of the present invention, and the change for carrying out non-creativeness work on this basis and making, all should belong to the present invention's
Within protection domain.
Claims (2)
1. the digital n- γ discriminating methods based on pulse width, it is characterised in that comprise the following steps：
(I) in mixed radiation field, carried out using detector and capture card from trigger-type digital waveform SnSample and store；It is described
Digital waveform SnFor activation threshold value VTPreceding m point and activation threshold value VTN-m point afterwards, wherein, m is more than 0, n and is more than m；
(II) by the sampling digital waveform SnCubic spline interpolation processing is carried out, the digital waveform S of 4 sampling points is obtained4n；
(III) the digital waveform S4nHandled through First derivative spectrograply, obtain digital waveform S4nPeak point P and peak value Vp；
(IV) by digital waveform S4nThe maximum amplitude of preceding 2m point be used as maximum white noise Vn；
(V) maximum white noise V is calculatednWith activation threshold value VTRatio R, and R also be pulse terminal amplitude VsWith peak value VpRatio
Value；
(VI) ratio R is utilized, it is counter to push away acquisition pulse terminal amplitude Vs, and determine the starting point K and cut off Q of pulse；Institute
It is pulse width W to state the width between starting point K and the cut off Q；
(VII) it is distributed using the pulse width W, obtains n- γ and screen result.
2. the digital n- γ discriminating methods according to claim 1 based on pulse width, it is characterised in that the step
(VI) in, comprising the following steps that for starting point K and cut off Q is determined：
(1) from digital waveform S4nPeak point P pointwises be moved to the left, and successively with pulse terminal amplitude VsCompare, if starting
Point amplitude VK≤ pulse terminal amplitude VsWhen, then K points are starting point, otherwise continue to compare to the left；
(2) from digital waveform S4nPeak point P pointwises move right, and successively with terminal amplitude VsCompare, if cut off width
Value VQ≤ pulse terminal amplitude VsWhen, then Q points are cut off, otherwise continue to compare to the right.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710569098.6A CN107272049B (en) | 2017-07-13 | 2017-07-13 | Digital n- γ discriminating method based on pulse width |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710569098.6A CN107272049B (en) | 2017-07-13 | 2017-07-13 | Digital n- γ discriminating method based on pulse width |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107272049A true CN107272049A (en) | 2017-10-20 |
| CN107272049B CN107272049B (en) | 2019-01-08 |
Family
ID=60072901
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710569098.6A Active CN107272049B (en) | 2017-07-13 | 2017-07-13 | Digital n- γ discriminating method based on pulse width |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107272049B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108490477A (en) * | 2018-03-14 | 2018-09-04 | 中国人民解放军火箭军工程大学 | A kind of n- γ pulses discriminating methods and system based on improvement centroid algorithm |
| CN108627867A (en) * | 2018-03-14 | 2018-10-09 | 中国人民解放军火箭军工程大学 | A kind of real-time discrimination systems of digitlization n- γ |
| CN110231505A (en) * | 2019-06-06 | 2019-09-13 | 西安交通大学 | A kind of waveform shake modification method based on cubic spline interpolation |
| CN110781906A (en) * | 2018-07-31 | 2020-02-11 | 中国工程物理研究院电子工程研究所 | N-gamma discrimination method based on semi-supervised support vector machine |
| CN112180421A (en) * | 2020-04-21 | 2021-01-05 | 宁波甬东核辐射监测有限公司 | Alpha and beta pulse discrimination method and device |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS58167988A (en) * | 1982-03-30 | 1983-10-04 | Solar Optical Japan Kk | Solid flying track detector material for measuring thermal neutron |
| US5880469A (en) * | 1995-01-31 | 1999-03-09 | Miller; Thomas Gill | Method and apparatus for a directional neutron detector which discriminates neutrons from gamma rays |
| CN101034120A (en) * | 2007-03-09 | 2007-09-12 | 中国科学院上海光学精密机械研究所 | Pulse shape measuring device and measuring method thereof |
| CN101303377A (en) * | 2008-06-11 | 2008-11-12 | 成都南迪测量科技有限公司 | Parameter acquisition method of discharging radiation pulse |
| CN101793514A (en) * | 2002-11-08 | 2010-08-04 | 皮克托米特里国际公司 | Method and apparatus for capturing, geolocating and measuring oblique images |
| CN102004260A (en) * | 2009-09-02 | 2011-04-06 | 成都理工大学 | Method for discriminating pulse width |
| EP2504721A1 (en) * | 2009-11-25 | 2012-10-03 | UT-Battelle, LLC | Shifting scintillator neutron detector |
| US8399849B1 (en) * | 2009-08-08 | 2013-03-19 | Redpine Signals, Inc | Fast neutron detector |
| CN103913766A (en) * | 2012-12-28 | 2014-07-09 | 株式会社东芝 | Digital rate meter and radiation monitoring system using digital rate meter |
| CN105093261A (en) * | 2015-08-20 | 2015-11-25 | 中国原子能科学研究院 | Neutron and gamma ray discrimination system of digital neutron spectrometer |
| CN105789001A (en) * | 2016-03-18 | 2016-07-20 | 南京瑞派宁信息科技有限公司 | Ion beam generation method and device |
| CN106443760A (en) * | 2016-09-14 | 2017-02-22 | 湖北锐世数字医学影像科技有限公司 | Scintillation pulse model-simplifying, reconstruction and energy-obtaining methods |
-
2017
- 2017-07-13 CN CN201710569098.6A patent/CN107272049B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS58167988A (en) * | 1982-03-30 | 1983-10-04 | Solar Optical Japan Kk | Solid flying track detector material for measuring thermal neutron |
| US5880469A (en) * | 1995-01-31 | 1999-03-09 | Miller; Thomas Gill | Method and apparatus for a directional neutron detector which discriminates neutrons from gamma rays |
| CN101793514A (en) * | 2002-11-08 | 2010-08-04 | 皮克托米特里国际公司 | Method and apparatus for capturing, geolocating and measuring oblique images |
| CN101034120A (en) * | 2007-03-09 | 2007-09-12 | 中国科学院上海光学精密机械研究所 | Pulse shape measuring device and measuring method thereof |
| CN101303377A (en) * | 2008-06-11 | 2008-11-12 | 成都南迪测量科技有限公司 | Parameter acquisition method of discharging radiation pulse |
| US8399849B1 (en) * | 2009-08-08 | 2013-03-19 | Redpine Signals, Inc | Fast neutron detector |
| CN102004260A (en) * | 2009-09-02 | 2011-04-06 | 成都理工大学 | Method for discriminating pulse width |
| EP2504721A1 (en) * | 2009-11-25 | 2012-10-03 | UT-Battelle, LLC | Shifting scintillator neutron detector |
| CN103913766A (en) * | 2012-12-28 | 2014-07-09 | 株式会社东芝 | Digital rate meter and radiation monitoring system using digital rate meter |
| CN105093261A (en) * | 2015-08-20 | 2015-11-25 | 中国原子能科学研究院 | Neutron and gamma ray discrimination system of digital neutron spectrometer |
| CN105789001A (en) * | 2016-03-18 | 2016-07-20 | 南京瑞派宁信息科技有限公司 | Ion beam generation method and device |
| CN106443760A (en) * | 2016-09-14 | 2017-02-22 | 湖北锐世数字医学影像科技有限公司 | Scintillation pulse model-simplifying, reconstruction and energy-obtaining methods |
Non-Patent Citations (1)
| Title |
|---|
| 左卓: "《一种单探头中子/伽马甄别能谱仪研究》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108490477A (en) * | 2018-03-14 | 2018-09-04 | 中国人民解放军火箭军工程大学 | A kind of n- γ pulses discriminating methods and system based on improvement centroid algorithm |
| CN108627867A (en) * | 2018-03-14 | 2018-10-09 | 中国人民解放军火箭军工程大学 | A kind of real-time discrimination systems of digitlization n- γ |
| CN110781906A (en) * | 2018-07-31 | 2020-02-11 | 中国工程物理研究院电子工程研究所 | N-gamma discrimination method based on semi-supervised support vector machine |
| CN110231505A (en) * | 2019-06-06 | 2019-09-13 | 西安交通大学 | A kind of waveform shake modification method based on cubic spline interpolation |
| CN112180421A (en) * | 2020-04-21 | 2021-01-05 | 宁波甬东核辐射监测有限公司 | Alpha and beta pulse discrimination method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107272049B (en) | 2019-01-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107272049A (en) | Digital n γ discriminating methods based on pulse width | |
| CN102753963B (en) | Use the transmission function of material and the method for identification material and device | |
| CN105319574A (en) | An alpha and beta pulse discriminating method | |
| CN107526096A (en) | It is a kind of to correct method of the temperature to SiPM gain effects | |
| CN102916683B (en) | A kind of Parameter-adjustablenuclear nuclear pulse simulation method | |
| CN110308476A (en) | A kind of particle radiation detection method and detection device | |
| Belli et al. | Conceptual design, development and preliminary tests of a compact neutron spectrometer for the JET experiment | |
| US20100308231A1 (en) | method and system for discrimination pulse shape | |
| JP2004108796A (en) | Radiation measurement device | |
| Pechousek et al. | Virtual instrumentation technique used in the nuclear digital signal processing system design: Energy and time measurement tests | |
| CN106842277A (en) | A kind of pile-up events processing method and processing device | |
| CN106841238A (en) | The method and system for positron annihilation dopplerbroadening spectrum based on anticoincidence | |
| JP7239921B2 (en) | Positron annihilation spectrometer | |
| Spieler | Pulse processing and analysis | |
| CN115980815A (en) | Method and device for acquiring energy of scintillation pulse, electronic equipment and storage medium | |
| CN108490477A (en) | A kind of n- γ pulses discriminating methods and system based on improvement centroid algorithm | |
| CN108535764A (en) | A kind of waveform to a variety of pulses screens algorithm | |
| EP3676640B1 (en) | Methods and systems for calibration of particle detectors | |
| Flaska et al. | Requirements and quantitative comparison of fast waveform digitizers for data-acquisition systems designed for nuclear nonproliferation applications | |
| Simpson et al. | Performance of four digital algorithms for γ− γ timing with LaBr3 (Ce) scintillators | |
| JPH06331750A (en) | Radiation energy spectrum measuring equipment | |
| CN1145809C (en) | Waveform discriminating method for gamma-ray detector | |
| Mitra et al. | Concurrent measurement of individual gamma-ray spectra during and between fast neutron pulses | |
| CN113126140A (en) | System and method for real-time discrimination of high-speed digital n/gamma waveforms | |
| CN107422362A (en) | A kind of noise-reduction method of energy spectrum and time spectrum |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20180129 Address after: Three road 610059 Sichuan city of Chengdu province Chenghua District Erxian Qiaodong No. 1 Applicant after: Chengdu University of Technology Applicant after: Sichuan University of Science & Engineering Applicant after: Southwest University of Science and Technology Address before: Three road 610059 Sichuan city of Chengdu province Chenghua District Erxian Qiaodong No. 1 Applicant before: Chengdu University of Technology Applicant before: Southwest University of Science and Technology |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |