Energy response correction for an electronic personal dosimeter using the channel ratio method

Energy response correction for an electronic personal dosimeter using the channel ratio method 1 2 corresp first-author Corresponding author, zhangyaohp@163.com Corresponding author, zhangyaohp@163.com zhangyaohp@163.com 1 1 2 1 3 1 2 1 3 National Institute of Measurement and Testing Technology, Chengdu 610021, China National Institute of Measurement and Testing Technology, Chengdu 610021, China Sichuan Radiation Technology Co., Ltd, Chengdu 610021, China Sichuan Radiation Technology Co., Ltd, Chengdu 610021, China Nimtt Measurement and Testing Technology, Ltd, Chengdu 610021, China Nimtt Measurement and Testing Technology, Ltd, Chengdu 610021, China An electronic personal dosimeter (EPD) mainly uses a Si-PIN photodiode as X- and gamma-ray detectors. The photon energy response of this instrument is inconsistent in the case of no correction, which seriously affects the accurate monitoring of personal dose equivalent H p (10) parameters for radiation workers. For this reason, in this paper we propose a method of combining composite screen detection technology, multichannel measurement technology, and the channel ratio method (CRM) to achieve accurate measurement of the personal dose equivalent parameters. According to China National Standard GB/T 13161-2003 and National Verification Regulation JJG 1009-2006, the instrument was tested in the energy range between 48 keV and 1.25 MeV. The experimental results showed that the difference of energy response to 137 Cs corrected by the new method was almost constant within ±6.0%, which fulfilled the ±30% requirement of GB/T 13161-2003 and JJG 1009-2006. Meanwhile, the method proposed obtained energy information regarding the radiation field. Electronic personal dosimeter (EPD) Si-PIN detector Photon energy response Channel ratio method (CRM). 1 Introduction 1 The electronic personal dosimeter (EPD), characterized by miniaturization, multi-functionality, wide measuring range, real-time measurement, and low power consumption, has been widely used in personal dose equivalent H p (10) monitoring in nuclear power environments, nuclear reactors, nuclear medicine, and other workplaces in which humans are exposed to radiation [ 1 1 - 3 3 1-3 ]. Since the Si-PIN semiconductor detector has advantages in sensitivity, volume, power consumption, price, etc., it is one of the key components of the EPD system [ 4 4 - 6 6 4-6 ]. At the same time, the Si-PIN detector has difficulty forming the equilibrium condition of charged particles because the atomic coefficient of the detector is higher than that of the tissue material and the sensitive layer is thin. For these reasons, the photon energy response of a Si-PIN detector is not consistent in the energy range from 48 keV to 1.25 Me V, especially in the low-energy region (less than 100 keV), where the photon energy response to 137 Cs is high. Because the energy response of a Si-PIN detector is not consistent, using the traditional single-channel counting method seriously affects the measurement accuracy of the instrument, so it is necessary to solve this problem by other methods. The method of correcting the energy response of a Si-PIN detector is one of the core technologies of an EPD. EPDs based on a Si-PIN detector mainly include three kinds of methods for energy response correction, namely the compensation material method, the spectrometer–dosimeter method, and the channel ratio method (CRM). The related parameters of an energy-compensated material can be calculated by numerical calculation or Monte Carlo simulation [ 7 7 - 10 10 7-10 ]. The core idea of this method is to design a kind of composite material that can achieve the purpose of off-setting the high response of the Si-PIN detector in the low-energy region [ 11 11 - 14 14 11-14 ]. The difficulty of this method resides in its inability to process bulk materials to meet performance requirements. At the same time, the measurement error of this method is correspondingly large, and the performance requirements of a Si-PIN detector are relatively high. The spectrometer–dosimeter method measures the pulse-height spectra of the Si-PIN detector, and the linear relationship between the pulse amplitude and the absorbed energy at the Si-PIN detector is obtained by experiments [ 15 15 , 16 16 ]. Because the influence of photons of different energies on the personal dose equivalent ( H p (10)) is different, the energy response correction of the Si-PIN detector can be realized by using the counts of each channel multiplied by weighting factors [ 17 17 ]. However, this method does not meet the basic requirements of the EPD because of its power consumption, size, cost, and other factors. Historically, some authors have put forward a simplified multi-channel spectrometric system by using a few comparators, and the energy response correction of the detector was done in the same way. This method can reduce the power consumption and the cost of the instrument, but it is difficult to adjust the threshold parameters, and the measurement accuracy is limited. Another method was called the dose conversion unit (DCU), which considered the pulse amplitude as well as the pulse number simultaneously. The DCU used a microprocessor to control a digital-to-analog converter (DAC), and the output of the DAC was fed into a single-channel analyzer (SCA) as low-level discriminator (LLD) values based on the dose conversion function. The final output of the DCU system was the count rates of pulses, which were directly proportional to the exposure dose rate in air [ 18 18 , 19 19 ]. From the relevant literature [ 17 17 - 19 19 17-19 ], it was found that this method had some limitations in the measuring range; in addition, the system structure and the calculation methods were relatively complex. The basic idea of the CRM is that there is a relationship between the channel ratio value and the beam energy. By changing the thickness of the additional filter material of the two channels, the beam energy can be identified using the CRM. By utilizing the conversion factor under the specific energy conditions, the personal dose equivalent parameters can be obtained. So far, the CRM has not been successfully applied in the EPD, but the method has potential application value in neutron-gamma resolution and beta measurements [ 20 20 - 22 22 20-22 ]. In summary, the CRM is more advantageous in terms of instrument structure than the spectrometer–dosimeter method, and its measurement accuracy is higher than the compensation material method. Therefore, in this paper we focus on the principle and method of using the CRM to correct the photon energy response of an EPD. 2 Mathematical relationship between channel ratio and beam energy 1 Figure 1 Figure 1 is a schematic of the two-channel ratio method. CM 1 and CM 2 are composite filter materials, the thicknesses are T 1 (cm) and T 2 (cm), the linear attenuation coefficients are μ 1 ( E ) (cm -1 ) and μ 2 ( E ) (cm -1 , and the sensitive area is S (cm 2 . D 1 and D 2 are the Si-PIN detectors, the thickness of the sensitive layer is T 3 (cm), the linear attenuation coefficient is μ 3 ( E ) (cm -1 , and the sensitive area is S (cm 2 . Channel 1 (CH 1 ) consists of a composite material CM 1 and a Si-PIN detector D 1 ; channel 2 (CH 2 ) consists of a composite material CM 2 and a Si-PIN detector D 2 . The back end of each channel also includes a charge-sensitive amplifier, a main amplifier, and a discriminator circuit (shown in Fig. 2 Fig. 2 ). The multichannel design method is based on the extension of the measuring range, and the correction principle of the small detector is the same as that of the large detector. Assuming that the particle fluence at the detector position is D S (cm -2 ), the particle energy is E (keV), and the particle numbers interacting with the Si-PIN detector are DM 1 and DM 2 . Assuming that particles interacting with the Si-PIN detector are recorded by the subsequent circuit, the following relationship exists: <math display="block" id="M1"><mrow><mrow><mo>{</mo><mtable columnalign="left"><mtr><mtd><mi>D</mi><msub><mi>M</mi><mtext>1</mtext></msub><mo>=</mo><msub><mi>D</mi><mi>s</mi></msub><msup><mi>e</mi><mrow><mo>−</mo><msub><mi>μ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><msub><mi>T</mi><mn>1</mn></msub></mrow></msup><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><msub><mi>μ</mi><mn>3</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><msub><mi>T</mi><mn>3</mn></msub></mrow></msup><mo stretchy="false">)</mo><mi>S</mi><mo>;</mo></mtd></mtr><mtr><mtd><mi>D</mi><msub><mi>M</mi><mn>2</mn></msub><mo>=</mo><msub><mi>D</mi><mi>s</mi></msub><msup><mi>e</mi><mrow><mo>−</mo><msub><mi>μ</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><msub><mi>T</mi><mn>2</mn></msub></mrow></msup><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><msub><mi>μ</mi><mn>3</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><msub><mi>T</mi><mn>3</mn></msub></mrow></msup><mo stretchy="false">)</mo><mi>S</mi><mo>.</mo></mtd></mtr></mtable></mrow></mrow></math> When μ 1 ( E ) is equal to μ 2 ( E ), Eq. (1) is further simplified as follows: <math display="block" id="M2"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>μ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mo>=</mo><mi>ln</mi><mo stretchy="false">(</mo><mi>D</mi><msub><mi>M</mi><mn>1</mn></msub><mo>/</mo><mi>D</mi><msub><mi>M</mi><mn>2</mn></msub><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><msub><mi>T</mi><mn>2</mn></msub><mo>−</mo><msub><mi>T</mi><mn>1</mn></msub><mo stretchy="false">)</mo><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> The linear attenuation coefficient μ 1 ( E ) can be calculated using Eq. (2). If the corresponding relation between μ 1 ( E ) and E is successfully established in a table, the energy information of the radiation field is obtained using the look-up table (LUT) method. Meanwhile, the channel ratio value and the beam energy should have a one-to-one relationship in the measurement range. μ 1 ( E ) should satisfy the monotone function characteristic in the energy range from 48 keV to 1.25 MeV. When the physical form of the filter material is considered, the photon mass attenuation coefficient [ μ 1 ( E )/ ρ ] is often used to describe the reduction of the photon fluence, where ρ is the density of the composite filter material (g/cm 3 ) [ 23 23 , 24 24 ]. After utilizing the mass attenuation coefficient, the basic feature remains unchanged. The calculation model under the monotone energy radiation field is now discussed. In the continuous spectrum of the radiation field (such as the X-ray diagnostic radiation source), Eq. (2) will be complex, and is transformed to http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10955375&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10955371&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=10955375&type= where ϕ ( E 1 ), ϕ ( E 2 ),..., ϕ ( E n ) are the particle fluences (cm -2 ), E 1 , E 2 ,..., E n are the energy of the X-ray particles (keV), and n indicates the particles of different energy. Although DM 1 , DM 2 , and the thickness of the composite material are known, the value μ 1 ( E ) is difficult to obtain using Eq. (3). Equations (1)-(3) are established on the premise that the detection system meets the narrow-beam geometry. Under broad-beam geometry (BBG) conditions, the influence of the scattered photons must be considered, and the μ 1 ( E ) value of the filter material is not accurately computed using the channel ratio value. Owing to the EPD having performance requirements for angle response and sensitivity, the detector unit is in the BBG condition, assuming that the scattered photons are recorded by the subsequent circuit under the BBG condition. When using the buildup factor to correct the scattered photons, Eq. (2) will be transformed into Eq. (4), and Eq. (3) will be even more complex: <math display="block" id="M3"><mrow><mtable columnalign="left"><mtr columnalign="left"><mtd columnalign="left"><mrow><msub><mi>μ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mo>=</mo><mi>ln</mi><mfrac><mrow><mi>D</mi><msub><mi>M</mi><mn>1</mn></msub><msub><mi>Δ</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><msub><mi>T</mi><mn>2</mn></msub><mn>,</mn><msub><mi>μ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow><mrow><mi>D</mi><msub><mi>M</mi><mn>2</mn></msub><msub><mi>Δ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mn>,</mn><msub><mi>T</mi><mn>1</mn></msub><mn>,</mn><msub><mi>μ</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow></mfrac><mo>/</mo><mo stretchy="false">(</mo><msub><mi>T</mi><mn>2</mn></msub><mo>−</mo><msub><mi>T</mi><mn>1</mn></msub><mo stretchy="false">)</mo><mn>.</mn></mrow></mtd></mtr></mtable></mrow></math> where Δ 1 ( E , T 1 , μ 1 ( E )) and Δ 2 ( E , T 2 , μ 1 ( E ))) are the buildup factors [ Δ 2 ( E , T 2 , μ 1 ( E )) / Δ 1 ( E , T 1 , μ 1 ( E )) > 1] for the scattered photons, which are related to the beam energy E , material properties μ 1 ( E ), and the thickness of the material. From Eq. (4), it is not possible to calculate μ 1 ( E ) directly under the BBG condition. Because Eq. (3) cannot be solved, μ 1 ( E ) under the continuous spectrum condition and the further derivation of Eq. (3) under the BBG condition is ignored. In summary, using the channel ratio value to solve for μ 1 ( E ), and then using μ 1 ( E ) to solve for the beam energy E , the measurement error of E under the BBG condition will increase. In practical application, through the narrow-spectrum series produced by an X-ray machine and a mono-energetic isotope radiation field, the relation between channel ratio value and beam energy can be obtained using the calibration method. The calibration method allows one to by-pass the process of the Eqs. (2)-(4) to solve for μ 1 ( E ) using the CRM, and the beam energy can be directly obtained using the LUT method. 3 Energy response correction using the channel ratio method 1 In practical application, it is essential to consider the consistency of aspects of the instrument, such as the filter material, the Si-PIN detector, and the circuit parameters. At the same time, the depletion layer of the Si-PIN detector and the surrounding material are not tissue equivalent. Calibration of the EPD under the standard radiation field can further simplify the mathematical calculation model. It is assumed that the counts per second (CPS) detected by CH 1 and CH 2 are CT 1 and CT 2 , and the CPS is linear with the dose rate. Assuming the dose rate at the position of the detector is D s ( μ Sv/h), the following relationship exists: <math display="block" id="M4"><mrow><mrow><mo>{</mo><mtable columnalign="left"><mtr><mtd><mi>C</mi><msub><mi>T</mi><mn>1</mn></msub><mo>=</mo><msub><mi>a</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><msub><mi>D</mi><mi>s</mi></msub><mo>;</mo></mtd></mtr><mtr><mtd><mi>C</mi><msub><mi>T</mi><mn>2</mn></msub><mo>=</mo><msub><mi>a</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><msub><mi>D</mi><mi>s</mi></msub><mn>.</mn></mtd></mtr></mtable></mrow></mrow></math> The sensitivity coefficients a 1 ( E ) [CPS/( μ Sv/h)] and a 2 ( E ) [CPS/( μ Sv/h)] of the two-channel detectors are calibrated under the condition of energy E in the standard radiation field. In practical application, the sensitivity of the detector, the measuring range, and other technical parameters must also be considered. The composite material described in this paper is mainly composed of Al, Cu, and scintillation materials. A simplified model of Eq. (6) can be used to overcome the uncertain properties of materials such as size and density. There is no exact analytical formula for the function y in Eq. (6), which can be solved using the calibration method in a radiation protection standard laboratory. <math display="block" id="M5"><mrow><mi>E</mi><mo>=</mo><mi>y</mi><mo stretchy="false">(</mo><mi>C</mi><msub><mi>T</mi><mn>2</mn></msub><mo>/</mo><mi>C</mi><msub><mi>T</mi><mn>1</mn></msub><mo stretchy="false">)</mo><mn>.</mn></mrow></math> Equation (6) can be used to calculate the energy information of the radiation field using the LUT method. Assuming b 1 ( E ) = 1/ a 1 ( E ) and b 2 ( E ) = 1/ a 2 ( E ), the accurate measurement of the personal dose rate equivalent can be accomplished using <math display="block" id="M6"><mrow><mrow><mo>{</mo><mtable columnalign="left"><mtr><mtd><msub><mi>D</mi><mi>s</mi></msub><mo>=</mo><mi>C</mi><msub><mi>T</mi><mn>1</mn></msub><msub><mi>b</mi><mn>1</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mo>;</mo></mtd></mtr><mtr><mtd><msub><mi>D</mi><mi>s</mi></msub><mo>=</mo><mi>C</mi><msub><mi>T</mi><mn>2</mn></msub><msub><mi>b</mi><mn>2</mn></msub><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mn>.</mn></mtd></mtr></mtable></mrow></mrow></math> The conversion factors b 1 ( E ) and b 2 ( E ) can be calculated using the LUT method. After using Eq. (7) in the energy range between 48 keV and 1.25 MeV, the photon energy response of the EPD is corrected. 4 Characteristics of reference radiation field 1 A medium-energy reference radiation field is based on the narrow-spectrum series produced by an X-ray machine. A high-energy reference radiation field is based on 137 Cs and 60 Co isotope radiation sources. The test conditions are summarized in Table 1 Table 1 , and the relevant standard values were obtained in a 1-liter spherical ionization chamber (type 32002) at PTW Freiburg GmbH. 1 3   1 2 Tube voltage(kV) 1 2 Tube current(mA) 4 1 Additional filtration 1 3 Half-value layer (HVL) for mm Cu 1 2 Distance (m) 1 2 Mean energy (keV) 1 2 Standard value ( μ Sv/h) 4 1 (mm) 1 1   1 1   1 1 Al 1 1 Cu 1 1 Pb 1 1 Sn 1 1   1 1   1 1   1 1 Narrow spectrum 1 1 60 1 1 1.5 1 1 4.0 1 1 0.61 1 1 / 1 1 / 1 1 0.24 1 1 3.0 1 1 47.9 1 1 1549 1 1   1 1 80 1 1 2.3 1 1 4.0 1 1 1.95 1 1 / 1 1 / 1 1 0.58 1 1 3.0 1 1 65.0 1 1 1498 1 1   1 1 100 1 1 4.3 1 1 4.0 1 1 4.70 1 1 / 1 1 / 1 1 1.10 1 1 3.0 1 1 83.1 1 1 1496 1 1   1 1 120 1 1 3.6 1 1 4.0 1 1 4.17 1 1 / 1 1 1.03 1 1 1.68 1 1 3.0 1 1 100.0 1 1 1515 1 1   1 1 150 1 1 0.6 1 1 4.0 1 1 0.65 1 1 / 1 1 2.11 1 1 2.33 1 1 3.0 1 1 117.7 1 1 1563 1 1   1 1 200 1 1 1.1 1 1 4.0 1 1 1.98 1 1 1.14 1 1 1.47 1 1 3.95 1 1 3.0 1 1 164.0 1 1 1436 1 1   1 1 250 1 1 0.8 1 1 4.0 1 1 / 1 1 1.98 1 1 2.95 1 1 5.15 1 1 3.0 1 1 207.5 1 1 1478 1 1 Isotope source 1 1 137 Cs 1 1   1 1   1 1   1 1   1 1   1 1 / 1 1 1.0 1 1 661.6 1 1 1148 1 1   1 1 60 Co 1 1   1 1   1 1   1 1   1 1   1 1 / 1 1 1.0 1 1 1250.0 1 1 7334 The narrow-spectrum series used in the research described in this paper mainly include N60, N80, N100, N120, N150, N200, and N250 [ 25 25 , 26 26 ]. The energy spectrum is shown in Fig. 3 Fig. 3 (normalized to the maximum counts). It can be seen from Fig. 3 Fig. 3 that the narrow-spectrum series have the characteristics of a continuous spectrum, but the attenuation law is similar to that of mono-energetic rays. The testing setup in the standard laboratory is shown in Fig. 4 Fig. 4 . The prototype model of the EPD tested was Nt2000. The EPD was located in the front of the water phantom and the Si-PIN detector was coincident with the center of the radiation field. 5 RESULTS AND DISCUSSION 1 The counts per second (CPS) detected by the two-channel detectors are shown in Table 2. Table 2. There were 10 measurements made in each energy condition, and the CPS values in Table 2 Table 2 are the arithmetic means of the 10 measured values in each energy condition. It can be seen from the results presented in the table that the photon energy responses of the two-channel detectors are inconsistent in the energy range between 48 keV and 1.25 MeV. 1 1 Mean energy (keV) 1 1 47.9 1 1 65.0 1 1 83.1 1 1 100.0 1 1 117.7 1 1 164.0 1 1 207.5 1 1 661.6 1 1 1250.0 1 1 Standard value ( μ Sv/h) 1 1 1549 1 1 1498 1 1 1496 1 1 1515 1 1 1563 1 1 1436 1 1 1478 1 1 1148 1 1 7334 1 1 CT 1 (CPS) 1 1 2895 1 1 2703 1 1 2251 1 1 2153 1 1 2081 1 1 1849 1 1 1735 1 1 732 1 1 3305 1 1 CT 2 (CPS) 1 1 809 1 1 1032 1 1 1001 1 1 1084 1 1 1187 1 1 1278 1 1 1355 1 1 680 1 1 3808 The EPD uses a single-channel detector (CH 1 ) to obtain the counts-dose rate conversion factor b 1 ( E ) [( μ Sv/h)/ CPS] (as shown in Fig. 5 Fig. 5 ) under the different beam energy conditions in the standard laboratory. The standard deviation of the counts-dose rate conversion factors of all energy conditions are between 0.016 and 0.039, and the standard deviation at the kink point of 207.5 keV is 0.023. If the photon energy response of CH 1 is not corrected, the difference response from low energy to high energy can be up to approximately 5 times. This situation will seriously affect the accurate measurement of the personal dose equivalent for radiation field staff. The plot shown in Fig. 6 Fig. 6 was calculated in accordance with Eq. (2). Because the filter material, detector, and back-end circuit parameters may be inconsistent, the calculation results in Fig. 6 Fig. 6 did not consider the density and thickness of the composite material. The standard deviation of the natural logarithm of the CT 1 / CT 2 of all energy conditions is between 0.022 and 0.046, and the standard deviation at the kink point of 207.5 keV is 0.024. Under the 137 Cs and 60 Co conditions, the CT 1 / CT 2 values are 1.077 and 0.868, respectively. As a result, the natural logarithm of the CT 1 / CT 2 is negative in the 60 Co condition. The reason for the above phenomenon is that the influence of the scattered photons in matter is neglected in Eq. (2). Because the purpose of this article is to determine the relationship between the channel ratio value and the beam energy, as long as they have a one-to-one relationship, the above problems have no effect on practical application. It can be seen from the curve of Fig. 6 Fig. 6 that the natural logarithm of the CT 1 / CT 2 satisfies the monotone function characteristic in the energy range from 48 keV to 1.25 MeV, which further proves the feasibility of using the channel ratio value to identify the beam energy. Figure 7 Figure 7 plots the discrete data points of the function y in Eq. (6). The standard deviation of the CT 2 / CT 1 of all energy conditions is between 0.013 and 0.019, and the standard deviation at the kink point of 207.5 keV is 0.019. It can be seen from Fig. 7 Fig. 7 that, although there is no linear or exponential relationship between the beam energy and the channel ratio value, there is a one-to-one relationship between them, which meets the requirement of Eq. (2). The relationship between the beam energy and the channel ratio value can be realized by polynomial fitting or by interpolation in the energy range between 48 keV and 1.25 MeV. Using the relation curve of Fig. 7 Fig. 7 , the beam energy can be directly identified by the channel ratio value. Figure 8 Figure 8 shows the relationship between the channel ratio value ( CT 2 / CT 1 ) and the counts-dose rate conversion factor of CH 1 . Owing to the change tendencies of the detector’s conversion factor over the entire energy range being different, the conversion factors in the low- and high-energy regions can be fitted. It can be seen from Fig. 8 Fig. 8 that the counts-dose rate conversion factor [ b 1 ( E )] of CH 1 can also be directly obtained from the channel ratio value ( CT 2 / CT 1 ). It can be seen from Figs. 6 Figs. 6 , 7 7 , and 8 8 that the standard deviation exhibits no obvious differences between the kink point and the other points under the experimental conditions. Table 3 Table 3 presents the measurement results after the energy response correction of the EPD. Ten measurements were made for each condition. The measured values in Table 3 Table 3 are the arithmetic mean of the 10 values in each energy condition. The personal dose rate equivalent H p (10) parameters can be calculated using the CRM, and the measurement accuracy compared to the traditional single channel method is significantly improved. 1 1 Energy(keV) 1 1 47.9 1 1 65.0 1 1 83.1 1 1 100.0 1 1 117.7 1 1 164.0 1 1 207.5 1 1 661.6 1 1 1250.0 1 1 Standard value ( μ Sv/h) 1 1 1549 1 1 1498 1 1 1496 1 1 1515 1 1 1563 1 1 1436 1 1 1478 1 1 1148 1 1 7334 1 1 Measured value ( μ Sv/h) 1 1 1510 1 1 1478 1 1 1486 1 1 1530 1 1 1503 1 1 1465 1 1 1501 1 1 1140 1 1 6850 1 1 Relative error (%) 1 1 -2.52 1 1 -1.34 1 1 -0.67 1 1 0.99 1 1 -3.84 1 1 2.02 1 1 1.56 1 1 -0.70 1 1 -6.60 1 1 Relative response to 137 Cs 1 1 0.98 1 1 0.99 1 1 1.00 1 1 1.02 1 1 0.97 1 1 1.03 1 1 1.02 1 1 1.00 1 1 0.94 1 1 Response difference to 137 Cs (%) 1 1 -1.83 1 1 -0.64 1 1 0.03 1 1 1.70 1 1 -0.32 1 1 2.74 1 1 2.27 1 1 0.00 1 1 -5.94 It can be seen from Fig. 9 Fig. 9 that the photon energy response to 137 Cs using the CRM is similar to a horizontal line. The standard deviation of the corrected energy response to 137 Cs is between 0.016 and 0.040.The relative error of the personal dose rate equivalent is within ±6.6%. The difference of the energy response to 137 Cs is within ±6.0%, which fulfills the ±30% requirement of GB/T 13161-2003 and JJG 1009-2006. 6 Conclusion 1 In this paper, we described the use of composite screen detection technology, multi-channel acquisition technology, and the CRM, by identifying the beam energy information, to achieve photon energy response correction for an EPD in the energy range from 48 keV to 1.25 MeV. At the same time, the principle and method of the photon energy response correction were derived, and the relevant features of the standard radiation field were described. The proposed method has the advantages of simple structure, easy implementation, low cost, low power consumption, and high measurement accuracy. The CRM was used to correct the photon energy response of the EPD, and the difference of the energy response to 137 Cs was almost constant within ±6.0%, which fulfilled the energy response requirement of GB/T 13161-2003 and JJG 1009-2006. From the related test results, it was shown that the proposed method can realize the accurate measurement of the personal dose equivalent of radiation workers, and the method has practical application value in a personal dose monitoring system. Although the channel ratio value and the beam energy have a one-to-one relationship, the relationship between them is not linear. The CRM is thus more suitable for personal dose equivalent measurement in mono-energetic conditions. When the EPD is used in a mixed radiation field, the measurement error of the personal dose equivalent will increase. In the future, it will be necessary to research the measurement method in a standard mixed radiation field. References L.A. Benevides , R.K. Piper , A. Romanyukha , A performance comparison of Thermo Fisher EPD-MK2 and TLD (LiF:Mg,Cu,P) as part of accreditation proficiency testing . Radiation Measurements . 71 , 183 - 186 ( 2014 ). doi: 10.1016/j.radmeas.2014.03.019 http://doi.org/10.1016/j.radmeas.2014.03.019 A performance comparison of Thermo Fisher EPD-MK2 and TLD (LiF:Mg,Cu,P) as part of accreditation proficiency testing J.G. Alves , Developments in standards and other guidance for individual monitoring . Radiation Measurements . 43 , 558 - 564 ( 2008 ). doi: 10.1016/j.radmeas.2007.12.044 http://doi.org/10.1016/j.radmeas.2007.12.044 Developments in standards and other guidance for individual monitoring F. Wang , M.Y. Wang , F.S. Tian et al ., Study on two-dimensional distribution of X-ray image based on improved Elman algorithm . Radiation Measurements . 77 , 1 - 4 ( 2015 ). doi: 10.1016/j.radmeas.2015.03.012 http://doi.org/10.1016/j.radmeas.2015.03.012 Study on two-dimensional distribution of X-ray image based on improved Elman algorithm J.F. Qian , T.Y. Shen , Nuclear Radiation Dosimetry ( Beijing , National Defense Industry Press , 2009 ), p, 199 - 278 (in Chinese) Nuclear Radiation Dosimetry Z. Li , X.G. Tuo , R. Shi et al ., Analytic fitting and simulation methods for characteristic X-ray peaks from Si-PIN detector . Nucl. Sci. Tech . 24 , 060206 ( 2013 ). doi: 10.13538/j.1001-8042/nst.2013.06.007 http://doi.org/10.13538/j.1001-8042/nst.2013.06.007 Analytic fitting and simulation methods for characteristic X-ray peaks from Si-PIN detector T. Ishikura , K. Aoyama , T. Fujimoto et al ., Development of a small semiconductor detector using a printed circuit silicon-board for a personal dosemeter . Radiation Measurements . 43 , 546 - 549 ( 2008 ). doi: 10.1016/j.radmeas.2008.01.009 http://doi.org/10.1016/j.radmeas.2008.01.009 Development of a small semiconductor detector using a printed circuit silicon-board for a personal dosemeter H. Sato , O. Maida , T. Ito , Structure optimization of diamond personal dosimeters based on Monte Carlo simulations . Diamond & Related Materials . 20 , 140 - 144 ( 2011 ). doi: 10.1016/j.diamond.2010.11.015 http://doi.org/10.1016/j.diamond.2010.11.015 Structure optimization of diamond personal dosimeters based on Monte Carlo simulations J. Barthe , Electronic dosimeters based on solid state detectors . Nucl. Instrum. Methods Phys. Res., Sect. B . 184 , 158 - 189 ( 2001 ). doi: 10.1016/S0168-583X(01)00711-X http://doi.org/10.1016/S0168-583X(01)00711-X Electronic dosimeters based on solid state detectors Z.G. Jiang , Y.G. Yuan , H.Y. Wang et al ., Development of a spherical tissue equivalent proportional counter for neutron monitoring . Nucl. Sci. Tech . 26 , 060403 ( 2015 ). doi: 10.13538/j.1001-8042/nst.26.060403 http://doi.org/10.13538/j.1001-8042/nst.26.060403 Development of a spherical tissue equivalent proportional counter for neutron monitoring P. Sun , Y.C. Fu , J. Hu et al ., Electron irradiation effects of radiochromic PCDA vesicle gel dosimeters . Nucl. Sci. Tech . 26 , 050301 ( 2015 ). doi: 10.13538/j.1001-8042/nst.26.050301 http://doi.org/10.13538/j.1001-8042/nst.26.050301 Electron irradiation effects of radiochromic PCDA vesicle gel dosimeters C.R. Chen , S.H. Jiang , H.H. Tseng et al ., Energy response and filter compensation of Si PIN photodiode for personal dosimetry application . IEEE T NUCL SCI . 40 , 857 - 862 ( 1993 ). doi: 10.1109/23.256674 http://doi.org/10.1109/23.256674 Energy response and filter compensation of Si PIN photodiode for personal dosimetry application H.H. Tseng , H.F. Chang , Energy compensation characteristics of Cu and Sn on low cost Si photodiodes for γ(x) dosimetry applications . IEEE T NUCL SCI . 39 , 1523 - 1527 ( 1992 ). doi: 10.1109/23.173237 http://doi.org/10.1109/23.173237 Energy compensation characteristics of Cu and Sn on low cost Si photodiodes for γ(x) dosimetry applications F. Wang , M.Y. Wang , Y.F. Liu et al ., Obtaining low energy dose with CMOS sensors . Nucl. Sci. Tech . 25 , 060401 ( 2014 ). doi: 10.13538/j.1001-8042/nst.25.060401 http://doi.org/10.13538/j.1001-8042/nst.25.060401 Obtaining low energy dose with CMOS sensors C.M. Deng , Z.S. Liu , Z.Y. Zhang et al ., The Research of Semiconductor Detector Test and Circuit for the SDM2000 Personal Dosimeter. China Nuclear Science and Technology Report. 2 , 122 - 132 ( 2004 ). (in Chinese) The Research of Semiconductor Detector Test and Circuit for the SDM2000 Personal Dosimeter. China Nuclear Science and Technology Report. 2 H.Z. Dizaji , T. Kakavand , F.A. Davani , Spectrometry and dosimetry of fast neutrons using pin diode detectors . Nucl. Instrum. Methods Phys. Res., Sect. A . 741 , 84 - 87 ( 2014 ). doi: 10.1016/j.nima.2013.12.018 http://doi.org/10.1016/j.nima.2013.12.018 Spectrometry and dosimetry of fast neutrons using pin diode detectors J.H. Wu , Y.J. Xu , S.L. Liu et al ., Energy spectrum measurement and dose rate estimation of natural neutrons in Tibet region . Nucl. Sci. Tech . 26 , 060202 ( 2015 ). doi: 10.13538/j.1001-8042/nst.26.060202 http://doi.org/10.13538/j.1001-8042/nst.26.060202 Energy spectrum measurement and dose rate estimation of natural neutrons in Tibet region M. Slapa , M. Traczyk , A New Concept of Dosimeter with Silicon Photodiodes . IEEE T NUCL SCI . 43 , 1855 - 1859 ( 1996 ). doi: 10.1109/NSSMIC.1995.510394 http://doi.org/10.1109/NSSMIC.1995.510394 A New Concept of Dosimeter with Silicon Photodiodes W. Lee , G. Cho , H.D. Kim , A Radiation Monitoring System With Capability of Gamma Imaging and Estimation of Exposure Dose Rate . IEEE T NUCL SCI . 49 , 1547 - 1551 ( 2002 ). doi: 10.1109/TNS.2002.1039698 http://doi.org/10.1109/TNS.2002.1039698 A Radiation Monitoring System With Capability of Gamma Imaging and Estimation of Exposure Dose Rate B.J. Lee , W. Lee , G. Cho et al ., Solid-state personal dosimeter using dose conversion algorithm . Nucl. Instrum. Methods Phys. Res., Sect. A . 505 , 403 - 406 ( 2003 ). doi: 10.1016/S0168-9002(03)01108-2 http://doi.org/10.1016/S0168-9002(03)01108-2 Solid-state personal dosimeter using dose conversion algorithm P. Huang , F. Fang , Y.D. Zhang et al ., The Detection Technology of PIN Silicon Photodiode Applied in the n / γ Mixed Radiation Field . Nuclear Electronics & Detection Technology . 31 , 1341 - 1343 ( 2011 ). (in Chinese) doi: 10.3969/j.issn.0258-0934.2011.12.010 http://doi.org/10.3969/j.issn.0258-0934.2011.12.010 The Detection Technology of PIN Silicon Photodiode Applied in the n / γ Mixed Radiation Field M. Jung , C. Teissier , P. Siffert , A proposal for fast neutron personal real time dosimetry in mixed n, γ fields . Nucl. Instrum. Methods Phys. Res., Sect. A . 431 , 264 - 272 ( 1999 ). doi: 10.1016/S0168-9002(99)00241-7 http://doi.org/10.1016/S0168-9002(99)00241-7 A proposal for fast neutron personal real time dosimetry in mixed n, γ fields M.C Chen , L.H. Tian , Mass thickness measurements for dual-component samples utilizing equivalent energy of X-rays . Nucl. Sci. Tech . 25 , 040202 ( 2014 ). doi: 10.13538/j.1001-8042/nst.25.040202 http://doi.org/10.13538/j.1001-8042/nst.25.040202 Mass thickness measurements for dual-component samples utilizing equivalent energy of X-rays H. Michael , An Investigation into Silicon PIN Diode Detectors for Dosimetry Applications. Master Thesis , University of Surrey ( 2010 ). An Investigation into Silicon PIN Diode Detectors for Dosimetry Applications. Master Thesis M.J. Berger , XCOM: Photon cross sections database , v 3.1 , National Institute of Standards and Technology (NIST), U.S.A . ( April 2010 ) http://www.nist.gov/physlab/data/xcom/index.cfm http://www.nist.gov/physlab/data/xcom/index.cfm [accessed August 2015 ]. XCOM: Photon cross sections database J.X. Chen , Y.F. Li , S.L. Liu et al ., JJG 1009-2006, Direct reading individual personal dose equivalent (rate) monitor for X and gamma radiation . Beijing, General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China (2006) . (in Chinese) Beijing, General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China (2006) Y. Liu , GB/T 13161-2003, Direct reading personal dose equivalent (rate) monitor—X, gamma radiation. Beijing, General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China ( 2003 ). (in Chinese) 10.1007/s41365-017-0256-1 This work was supported by the National Key Scientific Instruments To Develop Dedicated (2013YQ090811) 2017-02-17 2017-05-16 2017-05-18 2017-06-07 2017-07-01 2021-04-25 © 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-017-0256-1 http://dx.doi.org/10.1007/s41365-017-0256-1 . 2017 Nuclear Science and Techniques 07 28