Experimental validation of material discrimination ability of muon scattering tomography at the TUMUTY facility

Experimental validation of material discrimination ability of muon scattering tomography at the TUMUTY facility 1 first-author 2 3 1 corresp Corresponding author, zengzhi@tsinghua.edu.cn Corresponding author, zengzhi@tsinghua.edu.cn zengzhi@tsinghua.edu.cn 1 4 Key Laboratory of Particle and Radiation Imaging (Ministry of Education) and Department of Engineering Physics, Tsinghua University, Beijing 100084, China Key Laboratory of Particle and Radiation Imaging (Ministry of Education) and Department of Engineering Physics, Tsinghua University, Beijing 100084, China Department of Nuclear Engineering, University of California Berkeley, CA 94709, USA Department of Nuclear Engineering, University of California Berkeley, CA 94709, USA Department of Radiology and Biomedical Imaging, University of California, San Francisco, CA 94143, USA Department of Radiology and Biomedical Imaging, University of California, San Francisco, CA 94143, USA College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China College of Nuclear Science and Technology, Beijing Normal University, Beijing 100875, China Muon scattering tomography is believed to be a promising technique for cargo container inspection owing to the ability of natural muons to penetrate into dense materials and the absence of artificial radiation. In this work, the material discrimination ability of muon scattering tomography is evaluated based on experiments at the Tsinghua University cosmic ray muon tomography (TUMUTY) facility, with four materials: flour (as drugs substitute), aluminum, steel, and lead. The features of the different materials could be discriminated with cluster analysis and classifiers based on support vector machine (SVM). The overall discrimination precisions for these four materials could reach 70, 95, and 99% with 1-, 5-, and 10-min-long measurement, respectively. Muon tomography Cargo container inspection Material discrimination SVM classifier 1. Introduction 1 Cosmic rays, mainly originating from outside the solar system, produce showers of secondary particles, which are able to reach the sea-level. At the Earth’s surface, muons are the most numerous charged particles and have an intensity of approximately 1 cm -2 min -1 [ 1 1 ] . The muon scattering tomography technique [ 2 2 ] , which exploits the multiple Coulomb scattering of muons, has been applied in various fields, including cargo container inspection [ 3 3 - 6 6 3-6 ] and spent nuclear fuel monitoring [ 7 7 , 8 8 ] . This technique shows potential advantages owing to the ability of natural muons to penetrate dense materials and the absence of artificial radiation. Previous studies have shown that muon scattering tomography is more sensitive to high- Z materials, including uranium, tungsten, and lead. The performance of this technique with light materials is limited by the deep penetration of muons. However, during cargo container inspection, attention should be paid to both high- Z materials and light materials, such as drugs and explosives. If the identification ability of muon scattering tomography for light materials could be improved, it would play a more important role in cargo container inspection. In our previous work, we proposed the application of machine learning to identify drugs and explosives based on simulated datasets [ 9 9 ] . Evidence shows that the statistical results of the scattering densities of different materials are distinguishable. To perform material classification utilizing statistical results, a machine learning method based on support vector machine (SVM) was applied. In machine learning, SVMs are supervised learning models that analyze data for classification and regression analysis. The basic idea is trying to divide examples into different categories by a gap as wide as possible, which agrees with the problem of material classification. The SVM method is fully capable of performing materials identification if the features of different materials are derived. According to our previous simulation, drugs and explosives can be discriminated from air and metals with acceptable measurement durations using SVM classifiers. However, the spatial resolution of the detector was ignored and the simulated classifier did not distinguish the actual materials perfectly. In this work, a series of experiments were designed and conducted at the TUMUTY facility to validate the material discrimination ability of muon scattering tomography with four different materials [ 10 10 ] . 2. Methodology 1 2.1 Principle 2 Muon scattering tomography utilizes the multiple Coulomb scattering interactions between muons and materials. Muons traversing a medium could be deflected by many small-angle scatters. The net 2D scattering angles approximately follow a Gaussian distribution with zero mean and variance, given by [ 1 1 ] : <math display="block" id="M1"><mrow><msub><mi>σ</mi><mi>θ</mi></msub><mo>≅</mo><mfrac><mrow><mn>13.6</mn><mo> </mo><mtext>MeV</mtext></mrow><mrow><mi>β</mi><mi>c</mi><mi>p</mi></mrow></mfrac><msqrt><mrow><mi>X</mi><mo>/</mo><msub><mi>X</mi><mn>0</mn></msub></mrow></msqrt><mo>,</mo></mrow></math> where <math id="M2"><mi>c</mi></math> is the speed of light, <math id="M3"><mi>β</mi></math> is the ratio of the muons’ speed to the speed of light, <math id="M4"><mi>p</mi></math> is the muons’ momentum, <math id="M5"><mi>X</mi></math> is the thickness of the medium, and <math id="M6"><mrow><msub><mi>X</mi><mn>0</mn></msub></mrow></math> is the radiation length, which can be described by [ 1 1 ] : <math display="block" id="M7"><mrow><msub><mi>X</mi><mn>0</mn></msub><mo>≅</mo><mfrac><mrow><mi>A</mi><mo>⋅</mo><mn>716.4</mn><mtext> g</mtext><mo>/</mo><msup><mrow><mtext>cm</mtext></mrow><mn>2</mn></msup></mrow><mrow><mi>ρ</mi><mo>⋅</mo><mi>Z</mi><mrow><mo>(</mo><mrow><mi>Z</mi><mo>+</mo><mn>1</mn></mrow><mo>)</mo></mrow><mtext>ln</mtext><mrow><mo>(</mo><mrow><mn>287</mn><mo>/</mo><msqrt><mi>Z</mi></msqrt></mrow><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow></math> where <math id="M8"><mi>A</mi></math> is the atomic mass, <math id="M9"><mi>Z</mi></math> is the atomic number, and <math id="M10"><mi>ρ</mi></math> is the density. It can be seen that the variance of the scattering angles is dependent on the properties of the materials for muons with the same energies. Generally, materials with higher densities and larger atomic numbers always give larger scattering angle variances. Based on Eq.(1), the scattering density, which is dependent on the radiation length and muon momentum, can be defined as: <math display="block" id="M11"><mrow><mtext>λ</mtext><mrow><mo>(</mo><mrow><mi>p</mi><mo>,</mo><msub><mi>X</mi><mn>0</mn></msub></mrow><mo>)</mo></mrow><mo>=</mo><msubsup><mi>σ</mi><mi>θ</mi><mn>2</mn></msubsup><mo>/</mo><mi>X</mi><mo>≅</mo><msup><mrow><mrow><mo>(</mo><mrow><mfrac><mrow><mn>13.6</mn><mo> </mo><mtext>MeV</mtext></mrow><mi>p</mi></mfrac></mrow><mo>)</mo></mrow></mrow><mn>2</mn></msup><mfrac><mn>1</mn><mrow><msub><mi>X</mi><mn>0</mn></msub></mrow></mfrac><mo>,</mo></mrow></math> where <math id="M12"><mtext>λ</mtext></math> represents the scattering density and has units of rad 2 /cm. The approximation <math id="M13"><mrow><mi>β</mi><mi>c</mi><mo>≅</mo><mn>1</mn></mrow></math> for muons is used in this equation. Thus, the properties of the target materials could be found, if precise scattering densities are given. A number of algorithms have been developed to reconstruct the scattering densities, including the PoCA, MLSD, and MAP [ 11 11 , 12 12 ] . Among all the algorithms, the PoCA algorithm is the one with the lower computation and memory usage, and it performs well with short measurement durations. However, it should be noted that the energy of natural muons is not constant, which leads non-constant scattering densities for certain materials. The statistical fluctuation of reconstructed scattering densities is the main factor limiting the material discrimination ability, especially for low- Z materials. To decrease the effect of the statistical fluctuation, a discrimination method based on cluster analysis has been developed. The basic assumption of this method is that materials in neighbor voxels are very likely to be the same. Because containers are filled with the same materials during transportation and given that the voxels in this work have a size of 1 cm, which is much lower than the size of cargoes, it is reasonable to assume that voxels within a certain region (for example a region of 20×20×20 cm 3 ) are the same. By combining all the reconstructed scattering densities in the region, it is very likely that the effects of the statistical fluctuation are reduced. Our previous simulation has confirmed that the mean value and standard deviation of the scattering densities in a region can be utilized as features to perform materials discrimination [ 9 9 ] . 2.2 Experiment Setup 2 To evaluate the material discrimination ability of muon scattering tomography, a series of experiments were designed and conducted at the TUMUTY facility ( Fig. 1 Fig. 1 ). The facility consists of 6 layers of large scale position sensitive MRPC detectors, forming an efficient detection area with a size of 70×70×70 cm 3 . The detectors show a good spatial resolution of approximately 0.6 mm. However, the detection efficiency of each detector is approximately 80% at present and the total efficiency of the system is only approximately 26%. Three experiments with different designs were carried out with 4 different materials: flour, aluminum, steel, and lead ( Fig. 2 Fig. 2 ). Flour was chosen as the substitute of drugs, which are materials of great concerns for cargo container inspection. Flour was placed in a 20×20×20 cm 3 cubic plastic box. Aluminum, steel, and lead were selected to represent low- Z , medium- Z , and high- Z metals, respectively. The first experiment was designed to simulate a container of different cubic materials. The flour box and a 20×20×20 cm 3 cubic aluminum box were placed side by side with a central distance of approximately 30 cm. Two lead bricks were stacked together, forming a 20×10×10 cm 3 lead box. The three materials were placed on a platform in the middle of the TUMUTY facility. For the first experiment, measurements were taken over 10 days. The second experiment was carried out to show the discrimination ability for hidden materials. The flour box was surrounded by several steel cylinders to simulate a situation in which a low- Z non-metallic material is hidden inside medium- Z metals. The size of each steel cylinder was 10 cm in diameter and 20 cm in height. For the second experiment, measurements were taken over 5 days. The third experiment was similar to the second experiment in some aspects. The flour box was placed in the middle of a hollow steel box with a length of 40 cm, a height of 5 cm, and a thickness of 2 cm. For the third experiment, measurements were taken over 4 days. 3. Results and Discussion 1 3.1 Image Reconstruction 2 In this work, the PoCA algorithm was applied to reconstruct the scattering densities of spatial voxels. Sometimes a false triggering of a single detector may occur, which result in wrong muon tracks, leading to poor results. To avoid this, the trajectories of muons were fitted with the principle component analysis (PCA) method before reconstruction and the trajectories with poor linearity were rejected. It should be noted that the effective measuring time is not equal to the actual data collection time because of the detector efficiency and trajectory selection. The effective time can be calculated from the number of muons impinging on the detection area. A Monte Carlo simulation based on the actual geometry of the TUMUTY using the Geant4 toolkit was performed to calculate the muon acceptance of the system. Considering that the muon flux at sea level is approximately 1 cm -2 min -1 , the simulation showed that approximately 1600 muons per min impinged on the detection area. During the first experiment, 433560 valid trajectories were collected after data cut off, and the equivalent time was approximately 270 min. The equivalent times of experiments 2 and 3 were 220 and 50 min, respectively. The detection efficiency of the system was ignored temporarily because we wanted to investigate the inherent material discrimination ability of the muon scattering tomography system. The scattering densities in the three experiments were reconstructed using the PoCA algorithm. The size of the voxels was set to 1×1×1 cm 3 . The densities of each experiment were summed together along the Z and Y axes, as shown in Fig. 3 Fig. 3 , for a more effective demonstration. The regions of different materials were marked with rectangles in X-Y projection images. The color bars in Figure 3 Figure 3 represent the absolute values of the scattering densities. The shapes and sizes of lead and steel materials can be clearly identified in the projection images of all the experiments. The signal from the aluminum box in experiment 1 is low but still distinguishable. However, the scattering densities of the flour box in all the experiments are similar to that of the background, whose material is air. The inherent scattering density of flour is significantly low due to its small atomic number and low density. Thus, the statistical fluctuations caused by the detectors’ spatial resolution and PoCA algorithm make it more difficult to distinguish flour from the background. To evaluate the performance of the short time scanning, the datasets from the first 10 min of all the experiments were reconstructed separately, as shown in Fig. 4 Fig. 4 . In this situation, slight differences can be observed in the region of dense materials. However, it is still difficult to distinguish different materials, because of large fluctuations caused by the small quantity of incident muons. For datasets from shorter measuring times, the results of different materials are more difficult to be distinguished. 3.2 Material Discrimination Ability 2 To evaluate the statistic features of different materials with short measuring times, the dataset of each experiment was divided into subsets with results taken over 1, 5, and 10 min. For example, the dataset of experiment 1 can be divided into 270 subsets of 1-min measurements, 54 subsets of 5-min measurements, and 27 subsets of 10-min measurements. Each subset can be treated as an individual measurement. The following analysis is based on the subsets. Some regions were selected as regions of interest (ROIs) for further analysis depending on the actual positions of the materials, as shown in Fig. 3 Fig. 3 . Two cubic regions were selected for the steel cylinders in experiment 2 within the side of the flour box, marked with red rectangles in Fig. 3(b) Fig. 3(b) . Besides the measured materials, regions of the background were selected as well to represent the results from air, marked with blue rectangles in Figure 3 Figure 3 . 3.2.1 Discrimination of Materials from the Background 3 Owing to the PoCA algorithm principle, when the quantity of muons is not large enough, which is a common for short measuring times, the reconstructed density of many voxels is zero, because no PoCA points arise in these voxels. For voxels with high atomic numbers, there is a higher probability that PoCA points are reconstructed inside these voxels. The percentages of non-zero voxels in different material regions were calculated based on subsets with different measuring times, as shown in Fig. 5 Fig. 5 . The percentages of non-zero voxels is expected to follow a normal distribution based on the central limit theorem, which was fitted in Fig. 5 Fig. 5 . It can be seen that the percentages for the same material with different measuring times are nearly proportional to the measuring times. Furthermore, the denser is the material, the smaller is the percentage of non-zero voxels. The percentage of air is clearly lower than that of other materials. However, the difference between other materials is not so obvious, especially between flour and aluminum. Thus, air can be effectively distinguished from other materials. Among all the materials, flour has the biggest effect on air identification because of its low atomic number. The overlap probability of air and flour can be obtained using fitted normal distributions, which is approximately 2.31% for 1-min datasets. For 5- and 10-min datasets, the false probability is lower than 0.01%, which can be ignored. Moreover, the percentages of non-zero voxels are significantly low for all the materials. Even for dense materials, such as Fe, the value is only 0.5% per min, leading to the poor discrimination ability using the scattering density for short time measurement, as shown in Fig. 4 Fig. 4 . 3.2.2 Materials Discrimination 3 To analyze the differences between the measured materials, the reconstructed scattering densities in ROIs with different materials were extracted; the distributions are shown in Fig. 6 Fig. 6 . The values of air, flour, aluminum, and lead come from experiment 1 and the values of steel come from experiment 2. It can be seen that even in the region of low- Z material, large scattering densities are obtained. The deviations of the scattering densities for all the materials are large compared to the mean values, which makes it difficult to identify the material in a single voxel. However, the difference in the distributions of the scattering densities for the materials are still present. Obviously, the mean values and standard deviations for dense materials are larger than those for low- Z materials. Although it is difficult to discriminate a material from a single voxel, it is still possible to perform discrimination in a region of the same materials with multiple voxels according to the distribution of the scattering densities. The mean value and standard deviation are two of the most important features for the distribution, which could be used as features for material discrimination. The mean values and standard deviations of the scattering densities in all ROIs of the three experiments were calculated, as shown in Fig. 7 Fig. 7 . Each data point represents the mean value and standard deviation of the scattering densities in a certain ROI, obtained using sub-datasets with different measuring times. A difference between all materials is visible, especially for the 10-min dataset. With a 10-min measurement, the four materials (flour, aluminum, steel, and lead) could be clearly distinguished. For 1- and 5-min measurements, the difference is not as clear but it provides a possible method to perform material discrimination. However, the result of air are not as satisfying because it overlaps with the results of flour and aluminum, even with 10-min measurements. Nevertheless, this is not particularly important because the discrimination between air and other materials can be achieved by the percentage of non-zero voxels, as previously described. The comparison of the results of the sub-datasets from different measuring times show that different mean values and standard deviations from the same material are not the same. For example, the mean values of lead are distributed around 60 mrad 2 /cm in 1-min measurements, but those in the 5- and 10-min sub-datasets are 30 and 20 mrad 2 /cm, respectively. A decrease of the standard deviations with increasing measuring time was also observed. The same problem occurred in the previous simulation [ 9 9 ] . The simulated results show that the decreasing tendency will become less pronounced, i.e., they will become more stable, for measuring times longer than 20 min. This may be because of the limitation of the PoCA algorithm. When the size of the dataset is not large enough, the statistical fluctuations of the reconstructed scattering densities are very large. Despite the inconsistency of the results with different measuring times, the results for different materials with the same measuring time seem to be relatively stable. Thus, different material discrimination models could be trained according to the measuring time. The mean values and standard deviations of the scattering densities could be used as two features for further training of the classifiers. In this study, SVM-based classifiers were used because of their good performance for classified separation. Before training the SVM models, the features were normalized by: <math display="block" id="M14"><mrow><msup><mi>x</mi><mo>′</mo></msup><mo>=</mo><mfrac><mrow><mi>x</mi><mo>−</mo><mover accent="true"><mi>x</mi><mo>¯</mo></mover></mrow><mi>σ</mi></mfrac><mo>,</mo></mrow></math> where <math id="M15"><mi>x</mi></math> represents the raw value of the feature, <math id="M16"><mover accent="true"><mi>x</mi><mo>¯</mo></mover></math> represents its mean value, <math id="M17"><mi>σ</mi></math> represents its standard deviation, and <math id="M18"><msup><mi>x</mi><mo>′</mo></msup></math> represents the normalized value of the feature. Normalization is a common preprocessing method in the field of machine learning. In this study, only half of the datasets was used to train the SVM classifiers, and the other half was used as test datasets to evaluate the performance of the classification. Due to the inconsistency of features with different measuring times, different models were trained for 1-, 5-, and 10-min sub-datasets, as shown in Figure 8 Figure 8 . The region in Figure 8 Figure 8 was subdivided into four regions with different colors, representing the results of the SVM classification. The data points were actual results derived from sub-datasets of the three experiments. The performance of the classifiers was evaluated with test-datasets; the detailed precisions of different materials are listed in Table 1 Table 1 . The overall precisions of the four materials for 1-, 5-, and 10-min measurement are 70, 95, and 99%, respectively. It could be concluded that the longer is the measuring time, the more precise is the classification. Even when the measuring time is only 1 min, the precisions for most materials are higher than 70%, which is acceptable for short time measurements. With this method, material discrimination with muon scattering within short measuring times was achieved. 1 1 Measuring time (min) 5 1 Discrimination Precision 1 1 1 1 Flour 1 1 Al 1 1 Fe 1 1 Pb 1 1 Total 1 1 1 1 1 75% 1 1 50% 1 1 73% 1 1 76 % 1 1 70% 1 1 5 1 1 96% 1 1 100% 1 1 94% 1 1 89% 1 1 95% 1 1 10 1 1 100% 1 1 93% 1 1 100% 1 1 100 % 1 1 99% 4. Conclusion 1 In this work, we evaluated the material discrimination ability of muon scattering tomography with short measuring times in experiments at the TUMUTY facility. In these experiments, flour, aluminum, steel, and lead were chosen to represent materials with atomic numbers ranging from low to high. To reduce the statistical fluctuations within short measurements, a strategy based on cluster analysis was utilized. The percentage of the reconstructed non-zero scattering density could be used to discriminate materials from the background with an error rate of 2.31% for 1-min measurement. The mean value and standard deviation of the scattering densities were utilized to perform material discrimination. With SVM classifiers, the overall classification precisions for these four materials from 1-, 5-, and 10-min measurements are 70, 95, and 99%, respectively. However, the detection efficiency was ignored because we wanted to investigate the inherent discrimination ability of muon scattering tomography. If the efficiency is considered, poor results will be achieved due to the poor efficiency of the TUMUTY. However, this problem can be solved by using a facility with higher detection efficiency. Another limitation of this method is that the results obtained with different measuring times are not consistent. Overcoming this limitation requires further research. Reference K.A. Olive , K. Agashe , C. Amsler , et al ., Review of Particle Physics . Chinese Phys. C , 2014 , 38 : 090001 . Doi: 10.1103/PhysRevD.98.030001 http://doi.org/10.1103/PhysRevD.98.030001 . Review of Particle Physics S. Procureur , Muon imaging: Principles, technologies and applications . Nucl. Instruments Methods A 878 , 169 - 179 ( 2017 ). Doi: 10.1016/j.nima.2017.08.004 http://doi.org/10.1016/j.nima.2017.08.004 . Muon imaging: Principles, technologies and applications V. Antonuccio , M. Bandieramonte , D.L. Bonanno et al ., The Muon Portal Project: Design and construction of a scanning portal based on muon tomography . Nucl. Instruments Methods A , 845 , 322 - 325 ( 2017 ). Doi: 10.1016/j.nima.2016.05.006 http://doi.org/10.1016/j.nima.2016.05.006 . The Muon Portal Project: Design and construction of a scanning portal based on muon tomography G. Blanpied , S. Kumar , D. Dorroh et al ., Material discrimination using scattering and stopping of cosmic ray muons and electrons: Differentiating heavier from lighter metals as well as low-atomic weight materials . Nucl. Instrum. Methods A , 784 , 352 - 358 ( 2015 ). Doi: 10.1016/j.nima.2014.11.027 http://doi.org/10.1016/j.nima.2014.11.027 . Material discrimination using scattering and stopping of cosmic ray muons and electrons: Differentiating heavier from lighter metals as well as low-atomic weight materials P. Baesso , D. Cussans , C. Thomay , et al ., Toward a RPC-based muon tomography system for cargo containers . J. Instrum ., 9 , C10041 ( 2014 ). Doi: 10.1088/1748-0221/9/10/C10041 http://doi.org/10.1088/1748-0221/9/10/C10041 . Toward a RPC-based muon tomography system for cargo containers L. Frazão , J. Velthuis , C. Thomay , et al ., Discrimination of high-z materials in concrete-filled containers using muon scattering tomography . J. Instrum ., 11 , P07020 ( 2016 ). Doi: 10.1088/1748-0221/11/07/P07020 http://doi.org/10.1088/1748-0221/11/07/P07020 . Discrimination of high-z materials in concrete-filled containers using muon scattering tomography D. Poulson , J.M. Durham , E. Guardincerri , et al ., Cosmic ray muon computed tomography of spent nuclear fuel in dry storage casks . Nucl. Instrum. Meth. A 842 , 48 - 53 , ( 2017 ). Doi: 10.1016/j.nima.2016.10.040 http://doi.org/10.1016/j.nima.2016.10.040 . Cosmic ray muon computed tomography of spent nuclear fuel in dry storage casks J. M. Durham , D. Poulson , J. Bacon et al ., Verification of spent nuclear fuel in sealed dry storage casks via measurements of cosmic-ray muon scattering . Phys. Rev. Appl ., 9 , 44013 , ( 2018 ). Doi: 10.1103/PhysRevApplied.9.044013 http://doi.org/10.1103/PhysRevApplied.9.044013 Verification of spent nuclear fuel in sealed dry storage casks via measurements of cosmic-ray muon scattering Y. Zheng , Z. Zeng , M. Zeng , et al ., Discrimination of drugs and explosives in cargo inspections by applying machine learning in muon tomography . High Power Laser Part. Beams , 30 , 086002 ( 2018 ). DOI: 10.11884/HPLPB201830.180062 http://doi.org/10.11884/HPLPB201830.180062 Discrimination of drugs and explosives in cargo inspections by applying machine learning in muon tomography X. Wang , M. Zeng , Z. Zeng , et al ., The cosmic ray muon tomography facility based on large scale MRPC detectors . Nucl. Instrum. Meth. A 784 , 390 - 393 , ( 2015 ). Doi: 10.1016/j.nima.2015.01.024 http://doi.org/10.1016/j.nima.2015.01.024 The cosmic ray muon tomography facility based on large scale MRPC detectors L.J. Schultz , Dissertation , Portland State University , 2003 . Dissertation B. Yu , Z. Zhao , X. Wang , et al ., An MAP algorithm with edge-preserving prior for muon tomography . IEEE Nuclear Science Symposium & Medical Imaging Conference . IEEE , 2014 . Doi: 10.1109/NSSMIC.2014.7431083 http://doi.org/10.1109/NSSMIC.2014.7431083 . An MAP algorithm with edge-preserving prior for muon tomography 10.1007/s41365-019-0649-4 2019-02-22 2019-02-22 2019-04-11 2019-07-11 2019-08-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. 2019 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-019-0649-4 http://dx.doi.org/10.1007/s41365-019-0649-4 . 2019 Nuclear Science and Techniques 08 30