CN107272049B - Digital n- γ discriminating method based on pulse width - Google Patents
Digital n- γ discriminating method based on pulse width Download PDFInfo
- Publication number
- CN107272049B CN107272049B CN201710569098.6A CN201710569098A CN107272049B CN 107272049 B CN107272049 B CN 107272049B CN 201710569098 A CN201710569098 A CN 201710569098A CN 107272049 B CN107272049 B CN 107272049B
- Authority
- CN
- China
- Prior art keywords
- point
- digital waveform
- amplitude
- pulse
- pulse width
- 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
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
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
Abstract
The invention discloses a kind of digital n- γ discriminating method based on pulse width, comprising: carried out using detector and capture card from trigger-type digital waveform S in mixed radiation fieldnIt samples and stores；Digital waveform S will be samplednCubic spline interpolation processing is carried out, the digital waveform S of 4 sampling points is obtained4n；Digital waveform S4nIt is handled through First derivative spectrograply, obtains digital waveform S4nPeak point P and peak value Vp；Calculate maximum white noise VnWith activation threshold value VTRatio R, and also be pulse start-stop point amplitude VsWith peak value VpRatio；Using ratio R, counter push away obtains start-stop point amplitude Vs, 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, obtains n- γ using pulse width W and screens result.The discriminating method, which has many advantages, such as that calculating is easy, excludes subjective parameters setting, to be influenced, it is accurate to screen, in radiation detection technology field with very high practical value and promotional value.
Description
Technical field
The present invention relates to radiation detection technology fields, and in particular to the digital n- γ discriminating method based on pulse width.
Background technique
Since neutron discovery, it has been widely used in radiation detection field.In mixed radiation field, gamma-rays and neutron
It is also most important two kinds of rays that n, which is most common in radiation field, also the same in the occasion of neutron since two kinds of rays are not charged
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- γ is screened
Technology is to improve single radionetric survey precision in mixed radiation field, reduces one of the key of neutron (or gamma-rays) background.
N- γ examination is when generating different charged particles in detector using neutron and gamma-rays to cause different decaying
Between speed component ratio, to the difference of impulse waveform occur.Currently, the most common number n- γ discriminating method is that charge compares
Method screens incoming particle by fast component (or slow component) charge of impulse waveform and the proportional difference of total electrical charge, and this method is not
Foot place be, need specify 3 parameter calculated charge ratios, i.e., total electrical charge integral starting point, total electrical charge integral 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- γ is screened
The factor.Due to the specified not stringent reference standard of these three parameters, needs repeatedly to optimize determination according to quality factor, lead
Cause calculating process cumbersome, and the quality factor that n- γ is screened is larger by subjective impact.
Patent publication No. be CN201620788584.8 Chinese patent disclose it is a kind of digitize neutron spectrometer neutron,
γ discrimination system, 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 are artificially set by setting the one-to-one correspondence mapping relations of maximum value X% and maximum value Y% time difference
Fixed maximum value X% and maximum value Y% is changing value, and optimal value could be obtained by needing to carry out repeatedly calculation optimization.This method
The calculating process for the quality factor that optimization n- γ is screened is simplified to a certain extent, but two parameters set are without stringent ginseng
Standard is examined, still needs to optimize calculating, it is cumbersome there is no solving the problems, such as to calculate.
To sum up, it is necessary to n- γ screening techniques be proposed to improve, to obtain easy discriminating method.
Summary of the invention
The purpose of the present invention is to provide a kind of digital n- γ discriminating method based on pulse width, mainly solves existing skill
Examination algorithm present in art specifies parameter more, and calculating is complex, and quality factor optimization process is biggish by subjective impact
Problem.
To achieve the goals above, The technical solution adopted by the invention is as follows:
Digital n- γ discriminating method based on pulse width, is mainly resulted in using the die-away time difference of different particles
Corresponding pulses decay to required time difference when identical amplitude proportion, and it is wide that this difference can intuitively be embodied in pulse
On degree.Therefore n- γ examination can effectively be carried out according to pulse width, the specific steps are as follows:
Step 1: being carried out using detector and capture card from trigger-type digital waveform S in mixed radiation fieldnSampling is simultaneously
Storage, digital waveform SnFor activation threshold value VTPreceding m point and activation threshold value VTN-m point afterwards, wherein m is greater than 0, n and is greater than m.
Step 2: digital waveform S will be samplednCubic spline interpolation processing is carried out, the digital waveform of 4 sampling points is obtained
S4n。
Step 3: digital waveform S4nIt is handled through First derivative spectrograply, obtains digital waveform S4nPeak point P and peak value Vp。
Step 4: by digital waveform S4nPreceding 2m point maximum amplitude as maximum white noise Vn。
Step 5: calculating maximum white noise VnWith activation threshold value VTRatio R, and R is also pulse start-stop point amplitude VsWith peak
Value VpRatio.
Step 6: counter push away obtains pulse start-stop point amplitude V using ratio Rs, and determine the starting point K and cut off of pulse
Q, the specific steps are as follows:
(1) from digital waveform S4nPeak point P be moved to the left point by point, and successively with pulse start-stop point amplitude VsCompare, if
Starting point amplitude VK≤ pulse start-stop point amplitude VsWhen, then K point is starting point, otherwise continues to compare to the left；
(2) from digital waveform S4nPeak point P move right point by point, and successively with pulse start-stop point amplitude VsCompare, if
Cut off amplitude VQ≤ pulse start-stop point amplitude VsWhen, then Q point is cut off, otherwise continues to compare to the right.
Wherein, the width between starting point K and cut off Q is pulse width W.
Finally, being distributed using the pulse width W, obtains n- γ and screen result.
Compared with prior art, the invention has the following advantages:
(1) present invention is in such a way that detector and capture card trigger digit pulse wave, by the digit pulse wave of generation into
The processing of row cubic spline interpolation, only sets a parameter on waveform, i.e., maximum white noise VnWith activation threshold value VTRatio R,
In, activation threshold value VTTo determine parameter when impulse waveform use, ratio R only with maximum white noise VnIt is related, arbitrary point amplitude and arteries and veins
Rush start-stop point amplitude VsComparison obtains pulse width W, using pulse width to realize that n- γ is screened, simplifies cumbersome parameter
The calculating process of optimization.
(2) moreover, pulse start-stop point amplitude VsWith peak VpValue ratio R only by the maximum white noise of sample waveform and touching
It sends out threshold value and determines that quality factor are poor caused by can effectively avoid because of subjective factor, n- γ examination is impacted.
(3) present invention is using easy discriminating method, the result screened compared with the Charge comparison method of existing maturation,
The peak value of two kinds of particles can be more differentiated, and screens that process is more simple, and quality factor are not influenced by subjective setup parameter, the party
Method, which has many advantages, such as that calculating is easy, excludes subjective parameters setting, to be influenced, it is accurate to screen, in radiation detection technology field with very high
Practical value and promotional value.
Detailed description of the invention
Fig. 1 is the schematic diagram of the digital n- γ discriminating method based on pulse width.
Fig. 2 is the flow chart of the digital n- γ discriminating method based on pulse width.
Fig. 3 is that the hangover integral charge of Charge comparison method and total mark charge compare distribution map.
Fig. 4 is the charge that is obtained using Charge comparison method than the 2D distribution map with pulse amplitude
Fig. 5 is the pulse width distribution map of the digital n- γ discriminating method based on pulse width.
Fig. 6 is the pulse width of the digital n- γ discriminating method based on pulse width and the 2D distribution map of pulse amplitude.
Specific embodiment
Present invention will be further explained below with reference to the attached 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 when 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
The difference of pulse width W can effectively carry out n- γ examination.
As shown in Fig. 2, the digital n- γ discriminating method based on pulse width, discriminating method include the following steps:
Step 1: being carried out using detector and capture card from trigger-type digital waveform S in mixed radiation fieldnSampling is simultaneously
Storage, digital waveform SnFor activation threshold value VTPreceding m point and activation threshold value VTN-m point afterwards, wherein m is greater than 0, n and is greater than m.
Step 2: digital waveform S will be samplednCubic spline interpolation processing is carried out, the digital waveform of 4 sampling points is obtained
S4n。
Step 3: digital waveform S4nIt is handled through First derivative spectrograply, obtains digital waveform S4nPeak point P and peak value Vp。
Step 4: by digital waveform S4nPreceding 2m point maximum amplitude as maximum white noise Vn。
Step 5: calculating maximum white noise VnWith activation threshold value VTRatio R, and R is also pulse start-stop point amplitude VsWith peak
Value VpRatio.
Step 6: counter push away obtains pulse start-stop point amplitude V using ratio Rs, and determine the starting point K and cut off of pulse
Q, the specific steps are as follows:
(1) from digital waveform S4nPeak point P be moved to the left point by point, and successively with pulse start-stop point amplitude VsCompare, if
Starting point amplitude VK≤ pulse start-stop point amplitude VsWhen, then K point is starting point, otherwise continues to compare to the left；
(2) from digital waveform S4nPeak point P move right point by point, and successively with pulse start-stop point amplitude VsCompare, if
Cut off amplitude VQ≤ pulse start-stop point amplitude VsWhen, then Q point is cut off, otherwise continues to compare to the right.
Wherein, the width between starting point K and cut off Q is pulse width W.
Finally, being distributed using the pulse width W, obtains n- γ and screen result.
Experimental data use EJ299-33A detector electron energy on the phot-neutron source of 7MeV by 2.5Gs/s's
Digital oscilloscope sampling measures.Wherein, the hangover integral charge and total mark charge that Fig. 3 is Charge comparison method are than distribution map, Fig. 5
For pulse width distribution map of the invention, Fig. 4 is the charge that is obtained using Charge comparison method than the 2D distribution map with pulse amplitude,
Fig. 6 is the 2D distribution map using pulse width and pulse amplitude of the invention, can be seen that the method for the present invention more from comparison diagram
The peak value of two kinds of particles can be differentiated, n- γ examination quality factor are higher, it is more easy to screen process, and quality factor are not by subjectivity
Setup parameter influences.
Above-described embodiment is merely a preferred embodiment of the present invention, and it is not intended to limit the protection scope of the present invention, as long as using
Design principle of the invention, and the non-creative variation worked and made is carried out on this basis, it should belong to of the invention
Within protection scope.
Claims (1)
1. the digital n- γ discriminating method based on pulse width, which comprises the steps of:
(I) it in mixed radiation field, is carried out using detector and capture card from trigger-type digital waveform SnIt samples and stores；It is described
Digital waveform SnFor activation threshold value VTPreceding m point and activation threshold value VTN-m point afterwards, wherein m is greater than 0, n and is greater 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 S4nIt is handled through First derivative spectrograply, obtains digital waveform S4nPeak point P and peak value Vp；
(IV) by digital waveform S4nPreceding 2m point maximum amplitude as maximum white noise Vn；
(V) setup algorithm maximum white noise VnWith activation threshold value VTRatio R, and R is also pulse start-stop point amplitude VsWith peak value Vp
Ratio；
(VI) ratio R is utilized, counter push away obtains pulse start-stop point amplitude Vs, and determine the starting point K and cut off Q of pulse, from
Digital waveform S4nPeak point P be moved to the left point by point, and successively with pulse start-stop point amplitude VsCompare, if starting point amplitude VK≤
Pulse start-stop point amplitude VsWhen, then K point is starting point, otherwise continues to compare to the left, from digital waveform S4nPeak point P point by point to
Move right, and successively with start-stop point amplitude VsCompare, if cut off amplitude VQ≤ pulse start-stop point amplitude VsWhen, then Q point is cut-off
Otherwise point continues to compare to the right；Width between the starting point K and the cut off Q is pulse width W；
(VII) it is distributed using the pulse width W, obtains n- γ and screen result.
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 CN107272049A (en) | 2017-10-20 |
| CN107272049B true 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) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108627867A (en) * | 2018-03-14 | 2018-10-09 | 中国人民解放军火箭军工程大学 | A kind of real-time discrimination systems of digitlization n- γ |
| CN108490477A (en) * | 2018-03-14 | 2018-09-04 | 中国人民解放军火箭军工程大学 | A kind of n- γ pulses discriminating methods and system based on improvement centroid algorithm |
| CN110781906B (en) * | 2018-07-31 | 2023-05-30 | 中国工程物理研究院电子工程研究所 | N-gamma screening method based on semi-supervised support vector machine |
| CN110231505B (en) * | 2019-06-06 | 2020-12-08 | 西安交通大学 | Oscilloscope waveform jitter correction method based on cubic spline interpolation |
| CN112180421A (en) * | 2020-04-21 | 2021-01-05 | 宁波甬东核辐射监测有限公司 | Alpha and beta pulse discrimination method and device |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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 |
| CN106443760A (en) * | 2016-09-14 | 2017-02-22 | 湖北锐世数字医学影像科技有限公司 | Scintillation pulse model-simplifying, reconstruction and energy-obtaining methods |
Family Cites Families (7)
| 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 |
| CN101034120A (en) * | 2007-03-09 | 2007-09-12 | 中国科学院上海光学精密机械研究所 | Pulse shape measuring device and measuring method |
| CN101303377B (en) * | 2008-06-11 | 2011-08-10 | 成都南迪测量科技有限公司 | 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 |
| CN105093261B (en) * | 2015-08-20 | 2018-08-24 | 中国原子能科学研究院 | A kind of neutron, the γ discrimination systems of digitlization neutron spectrometer |
| CN105789001B (en) * | 2016-03-18 | 2018-05-01 | 南京瑞派宁信息科技有限公司 | The method and apparatus that a kind of ion beam produces |
-
2017
- 2017-07-13 CN CN201710569098.6A patent/CN107272049B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| 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 |
| 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 |
| CN106443760A (en) * | 2016-09-14 | 2017-02-22 | 湖北锐世数字医学影像科技有限公司 | Scintillation pulse model-simplifying, reconstruction and energy-obtaining methods |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107272049A (en) | 2017-10-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107272049B (en) | Digital n- γ discriminating method based on pulse width | |
| CN102753963B (en) | Use the transmission function of material and the method for identification material and device | |
| US20180113203A1 (en) | Interpolation Measurement of the Arrival Time and/or Amplitude of a Digitized Electronic Pulse | |
| CN102916683B (en) | A kind of Parameter-adjustablenuclear nuclear pulse simulation method | |
| Haugen et al. | Optimization of data‐acquistion rates in time‐correlated single‐photon fluorimetry | |
| Belli et al. | Conceptual design, development and preliminary tests of a compact neutron spectrometer for the JET experiment | |
| Yang et al. | Digital discrimination of neutron and γ ray using an organic scintillation detector based on wavelet transform modulus maximum | |
| US20100308231A1 (en) | method and system for discrimination pulse shape | |
| Liu et al. | Counting-loss correction method based on dual-exponential impulse shaping | |
| EP3529628A1 (en) | Ratio-reference measurement of the arrival time and/or amplitude of a digitized electronic pulse | |
| CN106841238A (en) | The method and system for positron annihilation dopplerbroadening spectrum based on anticoincidence | |
| Katabuchi et al. | Pulse-width analysis for neutron capture cross-section measurement using an NaI (Tl) detector | |
| CN105911078A (en) | Measurement method and measurement system for annihilation distribution of positron beam bunches in material | |
| CN113126140B (en) | System and method for real-time discrimination of high-speed digital n/gamma waveform | |
| Tardini et al. | Characterisation of a BC501A compact neutron spectrometer for fusion research | |
| CN109239765A (en) | A kind of imitative core pulse generation device and method | |
| EP3676640B1 (en) | Methods and systems for calibration of particle detectors | |
| WO2018081166A1 (en) | Ratio-reference measurement of the arrival time and/or amplitude of a digitized electronic pulse | |
| Ruben et al. | A new four channel pulse shape discriminator | |
| Simpson et al. | Performance of four digital algorithms for γ− γ timing with LaBr3 (Ce) scintillators | |
| Reva et al. | Computer methods of recognition and analysis of X-ray and gamma radiation parameters | |
| CN1145809C (en) | Waveform discriminating method for gamma-ray detector | |
| Nabavi et al. | Modified phase-only correlator for pileup correction for gamma-ray spectrometry | |
| CN114254495B (en) | Monte Carlo method-based simulated energy spectrum generation method | |
| Paepen et al. | Use of a CAEN digitiser for nuclear safeguards and security applications with a scintillator detector |
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 |