FPGA-based α/γ pulse shape discrimination for BaF₂ detector using 2 Gsps fast waveform sampling

FPGA-based α/γ pulse shape discrimination for BaF 2 detector using 2 Gsps fast waveform sampling 1 2 first-author 1 2 corresp Corresponding author. E-mail address: fengcq@ustc.edu.cn Corresponding author. E-mail address: fengcq@ustc.edu.cn fengcq@ustc.edu.cn 1 2 1 2 State Key Laboratory of Particle Detection and Electronics, University of Science and Technology of China, Hefei 230026, China State Key Laboratory of Particle Detection and Electronics, University of Science and Technology of China, Hefei 230026, China Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China Department of Modern Physics, University of Science and Technology of China, Hefei 230026, China Four FPGA-based α/γ pulse shape discrimination algorithms for BaF 2 detector are investigated and compared in this paper. A 2 Gsps fast waveform sampling board based on DRS4 chip is employed to sample the pulses. The test results with a 22 Na γ-source and the natural radioactivity of BaF 2 show good discrimination performance of the algorithms, with false rates around 1%. Small logical resource occupancy and short dead time are achieved. About 4400 slices are used in FPGA for pulse sampling and real-time discrimination altogether. BaF 2 Pulse extraction α/γ discrimination DRS4 FPGA 1 Introduction 1 Since the discovery of the intense, fast component of light pulse from BaF 2 crystal in 1980s [ 1 1 ] , BaF 2 scintillators have been extensively used in nuclear physics and elementary-particle physics for detecting γ-rays and light charged particles [ 2 2 , 3 3 ] , especially in medium and high energy regions. Commonly, BaF 2 crystals contain radium impurity, leading to background constituted mainly by α-particles, and in minor extent by γ-rays and electrons [ 4 4 ] . The possibility of discriminating particles from γ-rays and rejecting the background is an important issue in γ-spectroscopy. The radioluminescence spectrum of BaF 2 crystal contains a fast component with an extremely short decay time of 0.6 ns and a slow component with decay time of 620 ns.The fast/slow component ratio is large for γ-rays and small for heavy particles [ 5 5 - 7 7 5-7 ] . This allows one to discriminate heavy particles (such as α-particles) from γ-rays. Traditionally, analog methods were used for discrimination, with complex circuits and long dead times [ 8 8 , 9 9 ] . Over the past decade, with the advancement of electronics and digital signal processing technology, digital pulse shape discrimination techniques are extensively used [ 2 2 , 9 9 ] . To date, most works used ADCs for signal digitalization, and executed the digital algorithms by software, which would limit the processing speed. In this paper, we discuss four digital α/γ pulse shape discrimination (PSD) algorithms for BaF 2 scintillator, using 2 Gsps (Giga samples per second) fast waveform sampling and FPGA-based pulse information extraction. The DRS4, a switched capacitor arrays chip, is used for waveform sampling. The PSD algorithms are compared utilizing the natural radioactivity of BaF 2 crystal and a 22 Na γ-ray source. 2 Digital signal processing of pulses 1 Caused by the light quenching effect in BaF2 crystal, the fast/slow component ratio is greater for γ-induced pulse than α-induced pulse. Fig. 1 Fig. 1 shows two waveforms from a BaF2 crystal, the green waveform from α-particle, and the red from γ-ray. The α/γ PSD algorithms include the Charge Comparison Method (CCM)[ 10 10 ], Pulse Peak Analysis (PPA), Pulse Gradient Analysis (PGA)[ 11 11 ] and Finite impulse response filter Comparison Method (FCM). 2.1. Charge comparison method 2 The CCM is based on comparing the long and short integrals over different scales of one pulse. The long integral is a summation of the entire pulse samples, while the short integral is a summation of several samples near the maximum sample. In this paper, the optimum short integral is the summation of 16 samples: 5 samples before and 10 samples after the maximum sample. The regions for the long and short integrals are illustrated in Fig. 1. Fig. 1. The short integral corresponds to the region of the fast component, whilst the long integral represents the sum of the fast and slow components. The γ-induced events, with larger fast/slow component ratios than the α-induced events, should have larger short integrals for long integrals in the same length as those of the α-induced events. 2.2. Pulse peak analysis 2 The PPA is based on comparing the long integral and amplitude of the maximum sample, i.e. peak of the pulse. As shown in Fig.1 Fig.1 , for the same peak amplitudes, γ-induced pulses have smaller long integrals than α-induced pulses, due to larger fast/slow component ratios in BaF 2 scintillation. 2.3. Pulse gradient analysis 2 The PGA is based on comparing the long integral and the amplitude difference between the peak and the sample occurring at a defined time interval after the peak. The amplitude difference is often referred as the discrimination gradient, and the specific sample’s amplitude is known as the discrimination amplitude. The time interval was 5 ns in this work, considering different decay times of the fast and slow components, and timing property of the photomultiplier tube (PMT). The decay time of the digitized fast component was 3–4 ns, therefore the discrimination amplitude mainly depends on the amplitude of the slow component, while the peak mainly depends on amplitude of the fast component. Due to the larger fast/slow component ratios, the γ-induced events have bigger discrimination gradients for long integrals in same length as those of α-induced events. 2.4. Finite impulse response filter comparison method (FCM) 2 Through frequency spectrum analysis, the fast components in BaF2 pulses correspond to higher frequencies, at about 100 MHz, whereas the slow components are mostly lower than 2 MHz in the frequency domain. The pulses from BaF2 crystal are processed by a FIR low-pass filter. A proper low-pass filter reduces the amplitude of fast component much more than that of the slow component. In this study, the passband frequency of 3 MHz was chosen for the filter after tests from 2 MHz to 100 MHz. The filter-output-amplitude/filter-input-amplitude ratio of a γ-induced event, owing to its greater fast/slow component ratio, should be smaller than that of an α-induced event. This ratio is called the FIR-Peak-Ratio. 2.5. Pulse information extraction 2 When sampling BaF2 scintillation signal, the data rate is extremely high in consideration of the high sampling rate and high event rate (~50–200 Hz, depending on the source activity). Instead of sending the raw pulse shapes (a 2080 bytes package for each pulse) and analyzing the pulses in a computer as in most prior works, bad pulse elimination and pulse information extraction are executed in field programmable gate array (FPGA), and the pulse information (a 32 bytes package for each pulse) are sent to a computer. This reduces the data transfer by 64 times. At the USB transfer rate of 20 MBps, the transmission time decrease to 1.6 µs from 100 µs, which tremendously reduces the dead time correspondingly. After digitization of the pulse, a simple baseline restorer algorithm is applied. The baseline value is calculated as mean value of the first 16 samples, as shown in Fig.1 Fig.1 , and subtracted from the whole pulse form. Then all the information needed for the four PSD methods are extracted from each pulse form, namely long integral, short integral, peak amplitude, discrimination gradient and filter-output-amplitude. Fig. 2 Fig. 2 shows the flow of digital processing. Finally, together with packet flags, information of each pulse is sent to the computer through USB port as a 32 bytes package. 3 Experimental set-up and results 1 3.1. Experimental set-up 2 In this work, a cylindrical BaF 2 crystal of Φ 5 cm × 3 cm was instrumented with a XP2020Q PMT, and a collimated 22 Na source was placed 5 cm apart from it. The PMT was biased at −1.8 kV, which was considerably lower than the maximum rating value (−2.6 kV). To verify the PSD algorithms, a NaI detector was placed opposite to BaF 2 detector for coincidence. The PMT anodes were connected to inputs of a 2 Gsps fast sampling board based on a DRS4 chip [ 12 12 ] , the fourth version of Domino Ring Sampler (DRS) from Paul Scherrer Institute (PSI), Switzerland. The sampling board was capable of sampling 4 independent input channels of 1024-cell sampling depth, and with 300 MHz sine wave input, the signal-to-noise and distortion ratio (SINAD) of this sampling board reached 39.7dB, which corresponded to an effective number of bits (ENOB) of about 7 bits [ 13 13 ] . Using this sampling board, we achieved 7.5% energy resolution of the full-energy peak of 22 Na at 1.275 MeV. A low-end FPGA from Xilinx Spartan-3 family (XC3S5000) [ 14 14 ] was employed for digital signal processing on the sampling board. The experimental arrangement is shown schematically in Fig. 3. Fig. 3. 3.2. Results: Natural background 2 As previously mentioned, a characteristic of BaF2 crystals is the background caused by radium impurities constituted by α-particles, γ-rays and electrons. Thus, we used the internal natural radioactivity to calibrate digital pulse sampling and α/γ discrimination algorithms, and to determine the discrimination threshold. Only the BaF2 detector was used in this part. After extracting the pulse information, the results are presented in a 2-D scatter plot of the variables corresponded to each PSD algorithm, as shown in Fig. 4 Fig. 4 . Fig.4 (a) Fig.4 (a) shows the scatter plot of the CCM algorithm. The short integrals for 170000 events plotted against the long integrals: γ-rays are identified as points above the discrimination threshold (the black solid line), while α-particles are identified as points lay below it. Figs.4(b) and 4(c) Figs.4(b) and 4(c) are scatter plots of PPA and PGA. Particularly, in Fig.4 (d) Fig.4 (d) , for the FCM algorithm, with the FIR-peak-ratio being in the ordinate, and the ratio of short integral/long integral in the abscissa, γ-rays are identified as points below the discrimination threshold, while α-particles are identified as points above it. Furthermore, to evaluate the separation of γ-ray and α-particle regions, in Fig.5 Fig.5 the identification spectra are obtained by calculating the distance to the discrimination threshold line in the scatter plots for each event. The spectra are fitted with Gauss function shown as the red curves, and the full width at half maximum (FWHM) of the γ-peak as well as the α-peak in the spectra are compared in Table 1 Table 1 . More intuitively, figure-of-merit (FOM) values for the four methods can be calculated by Eq. (1), as shown in Table 1 Table 1 : 1 1 Methods 1 1 FWHM γ 1 1 FWHM α 1 1 |Peak γ −Peak α | 1 1 FOM 1 1 FAR(%) 1 1 CCM 1 1 1.39×10 4 1 1 1.95×10 4 1 1 37180 1 1 1.113 1 1 1.359 1 1 PPA 1 1 1.43×10 3 1 1 1.53×10 3 1 1 3568 1 1 1.205 1 1 1.339 1 1 PGA 1 1 1.46×10 3 1 1 2.00×10 3 1 1 4201 1 1 1.216 1 1 1.324 1 1 FCM 1 1 4.76×10 −2 1 1 9.89×10 −2 1 1 0.1651 1 1 1.127 1 1 1.526 <math display="block" id="M1"><mrow><mtext>FOM=</mtext><mrow><mo>|</mo><mrow><mfrac><mrow><msub><mrow><mtext>Peak</mtext></mrow><mtext>γ</mtext></msub><mo>−</mo><msub><mrow><mtext>Peak</mtext></mrow><mtext>α</mtext></msub></mrow><mrow><msub><mrow><mtext>FWHM</mtext></mrow><mtext>γ</mtext></msub><mo>−</mo><msub><mrow><mtext>FWHM</mtext></mrow><mtext>α</mtext></msub></mrow></mfrac></mrow><mo>|</mo></mrow></mrow></math> 3.3. Results: 22 Na γ-source 2 Due to the two-photon radiation of 22Na, the coincidence events of BaF2 detector and NaI detector can be identified as γ-ray events except for α-induced random coincidence events. Fig.6 Fig.6 shows distribution histograms of the distances to the discrimination threshold line of 20000 coincidence events from BaF2 detector, correspond to each of the four PSD algorithms. Most of the coincidence events are identified as γ-rays for all algorithms. False Alarm Rate (FAR) is defined as the proportion of coincidence events identified as α-particles, including α-induced random coincidences and wrong results of corresponding algorithms. The FARs of each algorithm are shown in Table 1 Table 1 . 3.4. Method comparison 2 Performance of the four PSD algorithms is compared in Table 1 Table 1 . Little variations are seen in the FOMs and FARs of the algorithms. FOM indicates the performance of separation, and the PPA and PGA algorithm are slightly better than the other two. According to the consistency of FARs for different algorithms, most of the misidentified events are likely belong to α-induced random coincidence events, being a fixed proportion among the whole events. FPGA resource occupancy is also a considerable issue. With the Xilinx Spartan-3 family (XC3S5000) FPGA, for CCM, PPA and PGA algorithms, 10% LUTs (6831 of 66560), and 13% Slices (4389 of 33280) were used. And for FCM algorithm, due to the use of the FIR filter, 59% LUTs (39545 of 66560), and 85% Slices (28618 of 33280) were used. Except for FCM algorithm, other algorithms show small FPGA resource utilization, considering that Spartan-3 family FPGA is low-cost and has small chip size. In summary, FCM achieves similar performance but costs much more FPGA logical resources than the other algorithms. Comprehensively consider the performance and FPGA usage, PGA turns out the best of all the four algorithms. 4 Discussion 1 Pulse information extraction is the key technique of the PSD algorithms. The good separations of α/γ events in identification spectra of each algorithm are mainly based on the information extracted in FPGA, well presenting notable features of different pulses. Besides, the extraction executed in FPGA greatly compresses the packet size, hence the reduced transmission time, shortened dead time, and reduced usage of FPGA logical resources. The events type can be judged in FPGA after obtaining the discrimination thresholds, instead of after transmitting to computer. Using several multiplier-adders and comparators after pulse information extraction, and the real-time PSD can be achieved. With the good performance, very short dead time and low logical resource usage, it’s prospective to use PSD algorithms like PGA in future BaF 2 detector experiments, such as CSNS-WNS (white neutron source at China Spallation Neutron Source) [ 15 15 , 16 16 ] , to discriminate α/γ events and eliminate α background. 5 Conclusion 1 Four FPGA-based digital PSD algorithms were investigated and compared in this work, similar PSD performance were observed in BaF 2 detector. Besides, basing on the pulse information extraction technique, the dead time and logical resource cost are significantly reduced, which makes the PSD algorithms suitable for real-time pulse processing in high-rate field experiments. References M. Laval , M. Moszyński , R. Allemand , et al ., Barium fluoride—inorganic scintillator for subnanosecond timing . J. Nuclear Instruments and Methods in Physics Research . 206 , 169 - 176 ( 1983 ). doi: 10.1016/0167-5087(83)91254-1 http://doi.org/10.1016/0167-5087(83)91254-1 Barium fluoride—inorganic scintillator for subnanosecond timing F. Amorini , E. De Filippo , P. Guazzoni , et al ., Digital pulse shape acquisition from BaF2: preliminary results . Paper presented at the 2006 IEEE Nuclear Science Symposium Conference Record , IEEE , San Diego , 29 Oct - 1 Nov 2006 Digital pulse shape acquisition from BaF2: preliminary results R. Novotny , Performance of the BaF2-calorimeter TAPS . J. Nuclear Physics B-Proceedings Supplements . 61 , 137 - 142 ( 1998 ). doi: 10.1016/S0920-5632(97)00552-5 http://doi.org/10.1016/S0920-5632(97)00552-5 Performance of the BaF2-calorimeter TAPS K. Wisshak K , F. Käppeler , Large barium fluoride detectors . J. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment . 227 , 91 - 96 ( 1984 ). doi: 10.1016/0168-9002(84)90105-0 http://doi.org/10.1016/0168-9002(84)90105-0 Large barium fluoride detectors E. Dafni , A note on the application of BaF2 scintillators to γ-ray and charged particle detection . J. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment . 254 , 54 - 60 ( 1987 ). doi: 10.1016/0168-9002(87)90481-5 http://doi.org/10.1016/0168-9002(87)90481-5 A note on the application of BaF2 scintillators to γ-ray and charged particle detection E. De Filippo , G. Lanzanó , A. Pagano , et al ., A study of light quenching in BaF2 crystals for heavy ions at intermediate energies . J. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment . 342 , 527 - 533 ( 1994 ). doi: 10.1016/0168-9002(94)90282-8 http://doi.org/10.1016/0168-9002(94)90282-8 A study of light quenching in BaF2 crystals for heavy ions at intermediate energies G. Lanzanó , A. Pagano , S. Urso , et al ., Using BaF2 crystals as detectors of light charged particles at intermediate energies . J. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment . 312 , 515 - 520 ( 1992 ). doi: 10.1016/0168-9002(92)90199-E http://doi.org/10.1016/0168-9002(92)90199-E Using BaF2 crystals as detectors of light charged particles at intermediate energies J. Gál , G. Kalinka , B.M. Nyakó , et al ., Particle discriminator for the identification of light charged particles with CsI (Tl) scintillator + PIN photodiode detector . J. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment . 366 , 120 - 128 ( 1995 ). doi: 10.1016/0168-9002(95)00560-9 http://doi.org/10.1016/0168-9002(95)00560-9 Particle discriminator for the identification of light charged particles with CsI (Tl) scintillator + PIN photodiode detector R. Aryaeinejad , J.K. Hartwell , D.F. Spencer , Comparison between digital and analog pulse shape discrimination techniques for neutron and gamma ray separation . Paper presented at the 2005 IEEE Nuclear Science Symposium Conference Record , IEEE , Fajardo , 23-29 Oct 2005 Comparison between digital and analog pulse shape discrimination techniques for neutron and gamma ray separation J.M. Adams , G. White , A versatile pulse shape discriminator for charged particle separation and its application to fast neutron time-of-flight spectroscopy . J. Nuclear Instruments and Methods . 156 , 459 - 476 ( 1978 ). doi: 10.1016/0029-554X(78)90746-2 http://doi.org/10.1016/0029-554X(78)90746-2 A versatile pulse shape discriminator for charged particle separation and its application to fast neutron time-of-flight spectroscopy B. D’Mellow , M.D. Aspinall , R.O. Mackin , et al ., Digital discrimination of neutrons and γ-rays in liquid scintillators using pulse gradient analysis . J. Nuclear Instruments and Methods in Physics Research Section A: Accelerators, Spectrometers, Detectors and Associated Equipment . 578 , 191 - 197 ( 2007 ). doi: 10.1016/j.nima.2007.04.174 http://doi.org/10.1016/j.nima.2007.04.174 Digital discrimination of neutrons and γ-rays in liquid scintillators using pulse gradient analysis S. Ritt , Design and performance of the 6 GHz waveform digitizing chip DRS4 . Paper presented at the 2008 IEEE Nuclear Science Symposium Conference Record , IEEE , Dresden , 19-25 Oct 2008 Design and performance of the 6 GHz waveform digitizing chip DRS4 J.H. Wang , L. Zhao , C.Q. Feng , et al ., Evaluation of a fast pulse sampling module with switched-capacitor arrays . J. IEEE Transactions on Nuclear Science . 59 , 2435 - 2443 ( 2012 ). doi: 10.1109/TNS.2012.2208656 http://doi.org/10.1109/TNS.2012.2208656 Evaluation of a fast pulse sampling module with switched-capacitor arrays Xilinx, Spartan-3 FPGA Family Data Sheet. (XILINX INC, 2013) , http://www.xilinx.com/support/documentation/data_sheets/ds099.pdf http://www.xilinx.com/support/documentation/data_sheets/ds099.pdf . Accessed 15 Aug 2016 Xilinx, Spartan-3 FPGA Family Data Sheet. (XILINX INC, 2013) J.Y. Tang , S.N. Fu , H.T. Jing , et al ., J. Chinese physics C 34 , 121 ( 2010 ) B. He , P. Cao , D.L. Zhang , et al ., Clock Distributing for BaF2 Readout Electronics at CSNS-WNS . arXiv: 1602.06635 [physics.ins-det] . Clock Distributing for BaF2 Readout Electronics at CSNS-WNS The authors would also thank Prof. Stefan Ritt and Dr. Jinhong Wang for their help. 10.1007/s41365-016-0173-8 This work was supported by the National Natural Science Fundation of China (No. 11205154). 2016-05-28 2016-08-16 2016-08-22 2016-12-24 2017-02-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2017 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-016-0173-8 http://dx.doi.org/10.1007/s41365-016-0173-8 . 2017 Nuclear Science and Techniques 02 28