بسم الله الرحمن الرحيم

1. بسم الله الرحمن الرحيم ”وقل ربي زدني علماً“

2. Developing a DSP Core using an FPGA Prototype for Scintillation Detector Signals Eng. Mahmoud Kamel Prepared by: Submitted to: Communication & Electronics Dept., Al Azhar University Supervised by: Prof. Dr. Ahmed Safwat Prof. Dr. Mahmoud Ashour Dr. Ashraf Aboshosha

3. Outline This core gives us all important features of the scintillation detector signals such as shaping, counting, pulse height and multichannel analyzing. The main purpose of this research work is to de-noise, compress and reconstruct the scintillation signals by which the processing speed, storage and precision will be improved.

4. Outline This core is implemented to apply the forward wavelet transform and interpolation technique. A new contribution of this framework arises from employing the interpolation techniques to reconstruct the signals where the mother wavelet and details are not required. Building a Multi-Channel Analyzer of the scintillation detector signals

5. Index of Content Scintillation detectors Importance of scintillation detectors Data Acquisition System Proposed digital processing algorithm Wavelets – Interpolation Technique Comparative study with the previous techniques Single channel and multi channel analyzer Conclusions and future work

6. Scintillation Detector Figure 1 : Schematic diagram of a scintillation detector

7. Scintillation Detector A scintillator is a material that emits light, scintillates, when absorbing radiation. The energy can be determined by measuring the pulse height spectrum. This is called spectroscopy. A scintillation detector is obtained when a scintillator is coupled to an electronic light sensor such as PMT or photodiode.

8. Importance of Scintillation Detectors Detection of mixed ionizing fluxes near nuclear objects. Radionuclide control of samples and radiation pollution. Determination of the type and energy of high-energy particles and products of their reactions with targets. Nuclear medicine (Gamma Camera, PET Tomography, …)

9. Data Acquisition System Figure 2 : The practical data acquisition system of scintillation detector Signals. (1) Scope , (2) high voltage source, (3) scintillator , (4) power supply

10. Why FPGA ? FPGA incorporates thousands of logic cells linked by programmable switches Highly parallel configurable digital signal processor A many channel signal processing was required in these detector to obtain a precise signals Availability of high-level design entry method FPGA designs easily changed, recompiled and low cost

11. FPGA Design Flow of the Solution Design Specification Design Entry/RTL Coding Behavioral or Structural Description of Design • RTL Simulation • Functional Simulation • Verify Logic Model & Data Flow • (No Timing Delays) LE • Synthesis • Translate Design into Device Specific Primitives • Optimization to Meet Required Area & Performance Constraints MEM I/O • Place & Route • Map Primitives to Specific Locations inside • Target Technology with Reference to Area & • Performance Constraints • Specify Routing Resources to Be Used

12. FPGA Design Flow of the solution Timing Analysis - Verify Performance Specifications Were Met - Static Timing Analysis tclk Gate Level Simulation -Timing Simulation - Verify Design Will Work in Target Technology Program & Test - Program & Test Device on Board

13. Pre-processing Phase 1-Wavelet based Decomposition 2- Interpolation based Reconstruction Pulse Shaping & Counting Multichannel analyzer Store & Show data Figure 3: The overall proposed solution

14. Pre Amplifier A SCA Counter Main Amplifier B MCA Figure 4: The proposed solution

15. Hardware System Figure 5: The FPGA XSC50k-Spartan II and the PC-based parallel interface

16. Pre-processing Phase • De-noising • Compression • Reconstruction

17. Effect of Noise onPulse Shaping & Counting Figure 6: Effect of noise on pulse shaping

18. Wavelets The wavelet analysis procedure is to adopt a wavelet prototype function, called an analyzing wavelet or mother wavelet. Wavelet transform decompose the original signal into different scales of resolution; these called the approximation and detail coefficients .

19. Wavelet Decomposition Levels d1 2 2 2 2 2 2 H X0 d2 H G d3 H G A3 G Figure 7: Three wavelet decomposition levels

20. Wavelet Families Haar Daubechies Biorthogonal Coifelt Symelet Myer

21. Haar Wavelet Design Pro: Allows good approximation with a subset of coefficients. It can be computed quickly and easily. Implemented easily by FPGA.

22. Design Block Diagram Figure 8: Design block diagram

23. Selecting the best Decomposition Level The quality of the compressed signals is the main criterion to select the best decomposition level in terms of Peak Signal to Noise Ratio (PSNR). The other similarity measure are Euclidean Distance (ED) , Cross Correlation coefficient (CC) and Mean Square Error(MSE).

24. Decomposition Levels Figure 9: Four approximation coefficients of Haar wavelet transform

25. Statistics of Four Decomposition Levels Table 1: Statistics of four levels Haar transform

26. Comparison of Different Mother Wavelets Table 2: Similarity measure of constructed and original signals of the different mother wavelets

27. Interpolation The Interpolation is a method of constructing new data points within the range of a discrete set of known data points. Interpolation is performed by fitting the supplied data with polynomial functions between data points and evaluating the appropriate function at the desired points.

28. Reconstruct Signals Using Interpolation Figure 10: a) Original signals. b) Transformed signal. c) Reconstructed signals

29. Interpolation Algorithms Nearest neighbor interpolation Linear interpolation Cubic Hermit Interpolation Cubic spline interpolation

30. Comparison of Applying Different Interpolation Techniques Table 3: Statistics of different interpolation techniques

31. Previous Pre-processing Techniques Accumulation Technique Median filter

32. Accumulation Technique Figure 11: Digital processing algorithm of scintillation detector signals

33. Median Filter Figure 12: Reconstructed signals using Median filter The value of an output sample is determined by the median of the neighborhood signals.

34. Comparison of the Preprocessing Results Table 4: Statistics of the preprocessing techniques

35. Figure 13: Pulse shaping after denoising

36. Figure 14: Pulse counting

37. Multi Channel Analyzer The MCA system is used to measure the height of each output pulse and the number of each output pulses simultaneously. By performing this operation for all detector events in a given interval the MCA generates a spectrum of the distribution of energy for a measured events with the y axis representing counts and the x axis representing channel value.

38. Multi Channel Analyzer Figure 15: Divided original signals into 16 channels

39. Multi Channel Analyzer Figure 16: Energy spectrum with 16 channels

40. Channel Calibration Energy channel values are converted into kilo electron volts with a channel-to-kilo electron volt conversion factor which is determined from a comparison of photo peak energies and channel location close to the energy of interest.

41. Conclusions One of the most important advantages of this system is the high compression rate (12.5%) using the interpolated wavelets Compared with the accumulation technique and median filtering, the proposed design achieved the best precision Capability of constructing MCA from SCA Coiflet is the best mother wavelet and Cubic spline is the best interpolation technique. Combining both of them for down and up sapling in wavelets is a new theoretical contribution of this framework

42. Future work Applying more complex wavelet filters. Modifying the proposed architecture to process more scintillator detectors. Employing the presented results as a base to identify radiation type and isotopes.

43. Thank You