WO2012163239A1 - 核辐射脉冲幅度的数字化方法及系统 - Google Patents
核辐射脉冲幅度的数字化方法及系统 Download PDFInfo
- Publication number
- WO2012163239A1 WO2012163239A1 PCT/CN2012/075885 CN2012075885W WO2012163239A1 WO 2012163239 A1 WO2012163239 A1 WO 2012163239A1 CN 2012075885 W CN2012075885 W CN 2012075885W WO 2012163239 A1 WO2012163239 A1 WO 2012163239A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- signal
- time
- circuit
- measured
- threshold
- Prior art date
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/16—Measuring radiation intensity
- G01T1/17—Circuit arrangements not adapted to a particular type of detector
Definitions
- the present invention relates to the field of nuclear radiation pulse signal processing technology, and more particularly to a digital method and system for nuclear radiation pulse amplitude.
- a nuclear case energy detector is an instrument that converts an incident pulse amplitude of a measured nuclear radiation into a measurable electrical signal, and the electrical signal is a series of separated pulse signals, wherein the amplitude value of the pulse signal represents The energy of the measured ray.
- the traditional method of measuring the amplitude of a nuclear radiation pulse signal is to use an analog-to-digital converter (ADC) to perform AD conversion on the input analog signal.
- ADC analog-to-digital converter
- This method is not suitable for large-channel-counting nuclear detector systems that require multi-channel parallel readout (such as nuclear imaging detector systems) because each channel uses a high-speed ADC that enables signal readout and digital electronics systems.
- the large scale is unacceptable, so a new type of nuclear radiation pulse amplitude digitization method is needed to meet the needs of large channel number parallel readout and digitization, while achieving high integration.
- the present invention converts the voltage amplitude value of a nuclear radiation pulse signal into a threshold time amount by using a time-to-digital conversion technique without significantly increasing the front-end readout circuit of the nuclear event energy detector, thereby realizing energy.
- the digitization of the signal Since the time measurement technique is relatively simple and mature, the present invention can not only achieve high-precision measurement for a digital system of large-channel digital nuclear radiation pulse signals, but most importantly, a large-channel number system can be realized because no additional front-end readout circuit is added. High integration.
- the present invention provides a possibility for digital integrated design of multi-channel core signals based on a single FPGA or ASIC.
- the present invention provides a nuclear pulse amplitude digitization method that utilizes a nuclear case energy detector, a constant ratio timing circuit, and a shaping circuit, and includes the following steps: Performing a constant ratio timing processing on the measured pulse signal outputted by the nuclear case energy detector by using a constant ratio timing circuit, and outputting a constant ratio timing signal, wherein the constant ratio timing signal is used to mark the arrival time Td of the nuclear case;
- the value of the measured shaped signal and the dynamic threshold signal are compared at the same time.
- the time is taken as the threshold time of the nuclear case.
- the dynamic threshold signal satisfies the following formula:
- Votl is the amplitude of the shaped signal
- Tot2 is the time of the over-threshold point obtained by comparing the other shaped signal with the dynamic threshold signal
- Vot2 is the amplitude of the other shaped signal.
- the invention also provides a nuclear pulse amplitude digitization system, comprising a nuclear case energy detector, a shaping circuit, a constant ratio timing unit, a comparator circuit and a time-to-digital conversion TDC unit, wherein the output of the nuclear event energy detector is respectively formed
- the circuit is coupled to an input of a constant ratio timing circuit, the constant ratio timing unit comprising a constant ratio timing circuit, the time digital conversion TDC unit comprising a TDC circuit;
- the forming circuit is configured to perform forming processing on the measured pulse signal output by the nuclear case energy detector, and output the measured forming signal;
- the constant ratio timing circuit is configured to perform a constant ratio timing processing on the measured pulse signal outputted by the nuclear case energy detector, and output a constant ratio timing signal, where the constant ratio timing signal is used to mark the arrival of the nuclear case Time Td ;
- the comparator circuit is configured to compare the value of the measured shaped signal and the dynamic threshold signal at the same time from the time of To, and output the pulse signal when the latter is changed from less than the former to not less than the former.
- the pulse signal outputted by the comparator circuit is used to mark the over-threshold time Tot of the nuclear case, and the Td and To are different by a preset duration;
- the TDC circuit records the Td, records the Tot, and converts the length of time between Td and Tot into a digital quantity;
- the dynamic threshold signal satisfies the following formula:
- Votl is the amplitude of the shaped signal
- Tot2 is the time of the over-threshold point obtained by comparing the other shaped signal with the dynamic threshold signal
- Vot2 is the amplitude of the other shaped signal.
- the invention also provides a digitization method for a nuclear radiation pulse amplitude, which linearly converts a nuclear radiation pulse amplitude into a time length, and realizes digitization of a nuclear radiation pulse signal by measuring the length of time, wherein The method includes the following steps:
- the above dynamic threshold signal needs to satisfy the following formula:
- Totl is an over-threshold point when any measured shaped signal is compared with the dynamic threshold signal
- Tot2 is an over-threshold time obtained by comparing another measured shaped signal with the dynamic threshold signal
- the invention also provides a digitizing system for a nuclear radiation pulse amplitude, comprising a nuclear case energy detector, a shaping circuit, a constant ratio timing unit, a comparator circuit and a time-to-digital conversion TDC unit, wherein the output of the nuclear case energy detector is respectively Forming circuit and constant ratio timing circuit
- the constant ratio timing unit includes a constant ratio timing circuit
- the time digital conversion TDC unit includes a TDC circuit
- the forming circuit is configured to shape a measured current signal output by the nuclear case energy detector, and output a measured forming signal
- the constant ratio timing circuit is configured to perform a constant ratio timing processing on the measured current signal output by the nuclear case energy detector, and output a constant ratio timing signal, where the constant ratio timing signal is used to mark the arrival time of the nuclear case Td;
- the comparator circuit is configured to compare values of the measured shaped signal and a dynamic threshold signal at the same time from a certain time To, when the dynamic threshold signal is less than the measured time at a certain time Tot
- the pulse signal is output, and the time Tot is the time exceeding the threshold point.
- the pulse signal is used to mark the threshold time Tot, and the Td and To are different by the preset duration. ⁇ ;
- the TDC circuit for recording the times Td and Tot and converting the length of time between the times Td and Tot into digital quantities
- the dynamic threshold signal satisfies the following formula:
- Totl is an over-threshold point when any measured shaped signal is compared with the dynamic threshold signal
- Tot2 is an over-threshold time obtained by comparing another measured shaped signal with the dynamic threshold signal
- FIG. 1 is a current waveform of 511 keV gamma ray according to an embodiment of the present invention
- FIG. 2 is a voltage waveform obtained by integrating current waveforms of different energy nucleus examples according to an embodiment of the present invention
- FIG. 3 is a schematic structural diagram of a constant ratio timing circuit according to an embodiment of the present invention.
- FIG. 4 is a flowchart of a method for digitizing a nuclear radiation pulse amplitude according to an embodiment of the present invention
- FIG. 5 is a schematic structural diagram of a nuclear radiation pulse amplitude digitization system according to an embodiment of the present invention
- FIG. 5 is a schematic structural diagram of a nuclear radiation pulse amplitude digitization system according to an embodiment of the present invention
- FIG. 6 is a schematic structural diagram of an active RC integration circuit according to an embodiment of the present invention
- FIG. 7 is a waveform diagram of an output signal of an active RC integration circuit according to an embodiment of the present invention
- FIG. 9 is a schematic diagram of a digital threshold voltage value according to an embodiment of the present invention
- FIG. 9b is a diagram showing another generation manner of a digital threshold voltage value according to an embodiment of the present invention
- Figure 9d is a schematic diagram of a circuit principle for generating a dynamic threshold using a dynamic threshold generation circuit in accordance with an embodiment of the present invention.
- Figure 9e is a specific circuit diagram of the circuit schematic design of the embodiment shown in Figure 9d.
- Figure 9f is a diagram showing the correspondence between the over-threshold time measured in accordance with the two dynamic threshold curves generated by the gain of the circuit of Figure 9e controlling the gain of the last-stage amplifying circuit and the amplitude value of the output signal of the nuclear case energy detector.
- FIG. 10 is still another schematic structural diagram of a nuclear radiation pulse amplitude digitization system according to an embodiment of the present invention.
- FIG. 11 is a schematic structural diagram of a dynamic threshold voltage generating unit based on an FPGA chip according to an embodiment of the present invention.
- FIG. 12 is a structural diagram of a constant ratio timing discriminating circuit according to an embodiment of the present invention.
- FIG. 13 is still another structural diagram of a constant ratio timing discriminating circuit according to an embodiment of the present invention
- FIG. 14 is a timing diagram of an embodiment of the present invention.
- the nuclear case energy detector of the present invention is mainly used for detecting the energy of gamma rays, and the nuclear case energy detector is used for converting a nuclear radiation pulse into a pulse signal to be measured, which is generally composed of a scintillation crystal such as LYSO (yttrium silicate). a scintillation crystal, a sodium iodide crystal, and a photoelectric converter (such as a photomultiplier tube, ⁇ ) are combined, and the measured pulse signal is output by a photoelectric converter, and the measured pulse signal is usually a current signal, so in the present invention It is also referred to as the measured current signal.
- LYSO yttrium silicate
- timing current signal is a current signal that is output by a nuclear event energy detector and is output separately for timing use while outputting the measured current signal.
- the timing current signal can be directly used for the constant ratio timing processing mentioned later to mark the arrival time of the nuclear case.
- Fig. 1 shows an example of a waveform of a measured current signal outputted by ⁇ corresponding to a detected output of 511 keV gamma rays.
- the rise time of the leading edge of the current waveform is less than Ins, and the fall time of the trailing edge is determined by the decay time of the LYSO. In one embodiment of the invention, it can be set to 42 ns, and the area of the current waveform represents the energy value detected by the secondary core case.
- the shaping circuit is configured to integrate the measured current signal (that is, the pulse current signal) outputted by the nuclear energy detector to obtain a voltage signal (or referred to as a shaped signal, the amplitude of which is related to the energy of the incident gamma ray) In proportion to).
- a voltage signal or referred to as a shaped signal, the amplitude of which is related to the energy of the incident gamma ray
- the voltage signal can be further low-pass filtered or amplified, or reshaped by other forming circuits commonly used in the art.
- the horizontal axis represents time, the unit is nanosecond, and the vertical axis represents voltage in volts.
- the measured current signal of the nuclear case of different energies The waveform of the integrated voltage signal differs only in the peak value (ie, the amplitude value) of the voltage signal waveform, and their peak time Tf, rise time, and fall time are the same.
- the rise time of the above shaped signal mainly depends on the luminescence attenuation characteristic of the nuclear case energy detector, and the fall time depends on the time constant of the shaping circuit. That is, after selecting the nuclear case energy detector and the shaping circuit, the rise time, peak time, and fall time of the shaped signal are fixed values.
- the “constant ratio timing” referred to in the present invention refers to a timing method in which a height is equivalent to a point at which the signal amplitude is a fixed ratio on a rising edge of the signal as a timing instant.
- the “constant ratio timing processing” in the present invention means that a constant ratio timing signal is generated and output at a constant proportional point of the rising edge of the input signal by a constant ratio timing circuit.
- the pulse current signal outputted by the nuclear energy detector can be divided into two paths, one is sent to the constant ratio timing circuit for processing to obtain the constant ratio timing signal, and the other is used as the measured current signal.
- the signal received by the constant ratio timing circuit may also be a shaped signal after the formed current signal is shaped. Since the principle of the constant ratio timing is to generate a timing signal at a constant proportional point of the input signal amplitude, the timing signal generated by the measured current signal after the forming by the constant ratio timing processing does not vary with the magnitude of the measured current signal. As can be seen from FIG. 2, for the shaped signals of different amplitudes, the time corresponding to the point at which the height of the rising edge corresponds to a fixed ratio of the amplitude of the pulse signal is taken as the timing time T1, and the time difference between the peak time Tf of the shaped signal and the peak time Tf of the shaped signal is A fixed value, ⁇ ⁇ ⁇ 1.
- the measured current signal is directly used to generate a constant ratio timing signal, but those skilled in the art should understand that the timing current signal output by the nuclear case energy detector can also be directly used, or The shaped signal of the measured current signal is used to generate a constant ratio timing signal.
- FIG. 3 is a schematic structural diagram of a constant ratio timing circuit according to an embodiment of the present invention.
- the constant ratio timing circuit includes an amplifier, a delay circuit, an attenuation circuit, and a comparator.
- the measured current signal output by the nuclear energy detector is divided into two paths through the amplifier, and one path is delayed by the delay circuit 31.
- the attenuation circuit 32 attenuates, and then the two signals are compared by the comparison circuit B1, and the timing of the rising edge of the constant ratio timing signal outputted by the comparison circuit B1 is at the current waveform.
- the value of the upper value is a fixed ratio to the maximum amplitude value of the current waveform, and the magnitude of the ratio is determined by the delay amount of the delay circuit 31, the attenuation coefficient of the attenuation circuit 32, and the rise time of the current waveform, which can be changed by changing the delay amount.
- the ratio is adjusted to achieve the purpose of adjusting the timing of the rising edge of the constant ratio timing signal.
- an embodiment of the present invention provides a method for digitizing a pulse amplitude of a nuclear radiation to realize digitization of the amplitude of the pulse signal to be measured outputted by the energy detector of the nuclear case.
- the method includes at least the following steps:
- the above dynamic threshold signal needs to satisfy the following formula:
- Vot2 Votl, ⁇ work which, over the threshold point in time (i.e., two waveform crossing point in time) Totl of Renyi Cheng-shaped signal and the dynamic threshold signal comparison obtained, Tot2 another shaped signal and the dynamic threshold signal
- the obtained threshold crossing time (that is, the time point at which the two waveforms intersect) is compared.
- the length of time between the above Td and Tot is referred to as a threshold time.
- the invention also discloses a nuclear radiation pulse amplitude digitization system capable of implementing the above method.
- Figure 5 illustrates the structure of one embodiment of the system, including:
- the input terminals of 3 are connected, the constant ratio timing unit includes the above-described constant ratio timing circuit, and the time digital conversion TDC unit 5 includes a TDC circuit.
- the forming circuit 2 is configured to form a processed current signal outputted by the nuclear case energy detector 1 and output a measured forming signal;
- the constant ratio timing circuit is configured to perform a constant ratio timing processing on the measured current signal outputted by the nuclear case energy detector 1, and output a constant ratio timing signal, and the arrival time of the constant ratio timing signal is Td;
- the comparator circuit 4 is configured to compare the value of the measured shaped signal and the dynamic threshold signal at the same time from a certain time To, when the latter changes from less than the former to not less than the former when at a time Tot And outputting a pulse signal, wherein the time To is different from the arrival time of the constant ratio timing signal by Td by a preset duration ⁇ , and the time Tot is a threshold crossing time;
- Vot ⁇ Vot2 ⁇ Votl the amplitude Vot of the measured shaped signal is linearly proportional to (Tot-Td+ ⁇ ) or (Tot-Td- ⁇ ). Since both Td and Tot can be measured, and ⁇ is a predetermined known value, the length of time between Td and Tot is digitized, and the signal amplitude can be digitized.
- the forming circuit may be an active integrating circuit (e.g., an active RC integrating circuit).
- Fig. 6 shows a structure of an active RC integration circuit.
- the active RC integration circuit comprises a parallel resistor R and a capacitor C. Both ends of the resistor R and the capacitor C are connected to one input terminal Vm and an output terminal Vout of an operational amplifier, and the other input terminal of the operational amplifier Ground.
- the active RC integration circuit integrates the current signal output by the nuclear instance energy detector 1, and the output voltage waveform can be seen in FIG.
- the resistance ⁇ 5 ⁇ ⁇
- the capacitance C 18 pf in the active RC integration circuit
- the peak time of the output voltage waveform is 60 ns
- the time to return to the baseline after the peak is greater than 300 ns.
- the difference in voltage waveform after integration of different energy nuclear cases is only the peak value of the voltage waveform, that is, the amplitude of the voltage waveform, and their peak time and fall time are the same.
- the method further includes the following step S5: generating a dynamic threshold signal.
- the present invention there are two methods for generating the above dynamic threshold signal, one of which is generated based on a series of digital threshold voltage values stored in advance; the other is generated in real time by a dynamic threshold generating circuit.
- the dynamic threshold signal is generated before the time of To, and the comparator circuit is input as one of the inputs at the time of To. Since the value of the threshold corresponding to the shaped signal of any amplitude can be calculated in advance according to the above formula 1, the values can be stored in a threshold sequence (ie, the series of digital threshold voltage values described above).
- FIG. 9 is a schematic diagram of a digital threshold voltage value according to an embodiment of the present invention
- FIG. 9b is a diagram of another generation manner of a digital threshold voltage value according to an embodiment of the present invention.
- the series of digital threshold voltage values are generated as follows: 1. Determine the time interval ⁇ and the number of the above-mentioned series of digital threshold voltage values N, N are positive integers not less than 1.
- N is equal to five.
- N shaped signals Si Sw are selected such that the amplitude of S! ⁇ is A, the amplitude of Si is A/N, and the amplitude values of S 2 to S N _i are sequentially increased by A/N.
- the amplitudes of S Ss are 0.5V, 1.0V, 1.5V, 2.0V, and 2.5V, respectively.
- the rise time, peak time and fall time of the shaped signal are fixed values. Therefore, N shaped signals can be easily generated using a signal generator or software simulation, the N shaped signals having the same shape as the measured shaped signals.
- the peak time at which the shaped signal is selected is the time origin.
- the value of the voltage corresponding to the shaping signal S M at the time ⁇ + ⁇ is stored as a digital threshold voltage value, the ⁇ is a positive integer, and 1 ⁇ ⁇ , the series of digital threshold voltage values are obtained.
- the series of digital threshold voltage values is, for example, five, which may constitute a discrete dynamic threshold curve.
- the above series of digital threshold voltage values are generated as follows: 1. Determine the maximum time measurement value Tmax.
- the maximum time measurement Tmax refers to the desired threshold time corresponding to the maximum measurable signal amplitude value of the system.
- the maximum time measurements corresponding to the three dynamic threshold curves D1-D3 of the embodiment shown in Figure 9b are 156 ns, 312.5 ns, and 625 ns, respectively.
- the peak time of the measured shaped signal is selected as the time origin.
- N is a positive integer not less than 1.
- N is equal to 5.
- those skilled in the art can flexibly design the values of the time intervals ⁇ and N according to the precision and other needs, and will not be described here. 5. Selecting N shaped signals
- the amplitude of the S N is A
- the amplitude of the S1 is A/N
- the amplitude value of the S 2 to is incremented by A/N.
- the voltage value, the above N digital threshold voltage values may constitute a discrete dynamic threshold curve.
- ⁇ can be selected to be smaller.
- a 12-bit DAC is used, the conversion clock is 160 MHz, and the clock period is 6.25 ns. Still taking the above 156 ns, 312.5 ns, and 625 ns as an example, they are divided into 25, 50, and 100 data calculation points by equal intervals of 6.25 ns (for convenience of viewing, FIG. 9b shows only five computer points, corresponding thereto, It is also necessary to select 25, 50 and 100 shaped signals of different amplitudes.
- the threshold time has a good linear relationship with the amplitude value of the output signal of the shaping circuit.
- An arbitrary waveform signal generator can be used to generate the same signal waveform as the above-mentioned nuclear case energy detector output waveform, and the amplitude of the waveform can be controlled to obtain no The same measured value.
- Figure 9c shows the correspondence between the over-threshold time measured according to the three dynamic threshold curves D1-D3 shown in Figure 9b and the amplitude value of the nuclear event energy detector output signal. It can be observed that, using the technical solution provided by the present invention, the threshold time and the amplitude value of the output signal of the forming circuit do have a good linear relationship.
- FIG. 11 is a schematic diagram of a dynamic threshold generating circuit 6 for generating a dynamic threshold from a series of digital threshold voltage values stored in advance based on an FPGA (Field Programmable Gate Array) chip.
- the dynamic threshold generating circuit 6 includes a static memory (RAM), a digital-to-analog conversion circuit (DAC), and a low-pass filter.
- the static memory is used to store a pre-computed series of digital threshold voltage values and is coupled to a write read logic module.
- the write read logic module accepts timing signals from a CFD (constant time timing circuit) for controlling the writing and reading of the static memory.
- the digital to analog conversion circuit is coupled to the static memory and receives a series of digital threshold voltage values read from the static memory for converting the series of digital threshold voltage values to analog signals for output to the low pass filter.
- the low pass filter is coupled to the digital to analog conversion circuit for low pass filtering the analog signal to produce a dynamic threshold voltage signal.
- a series of digital threshold voltage values calculated by the PC are stored in the static memory of the FPAG, and the write read logic module sets the above one according to a preset time interval (for example, 6.25 ns).
- the series digital threshold voltage values are read out one by one, and after a digital analog conversion circuit and a low-pass filter circuit, a dynamic threshold voltage signal of continuous time and continuous amplitude is generated.
- a higher resolution DAC chip can be selected, and a higher conversion clock rate can be selected.
- the dynamic threshold signal is generated by the dynamic threshold generation circuit in real time.
- the hardware circuit can be used to simulate the signal acquisition process of the detector detecting particles, and generate a voltage signal or a current signal having the same shape as the measured signal output by the nuclear instance energy detector. , and after the nuclear case
- the shaping circuit of the detector output signal is formed by the same forming circuit, multiplied by the linear ramp voltage signal generated by the other circuit at the peak of the signal, and then amplified by the amplifier circuit to be a continuous dynamic threshold signal.
- the embodiment of the present specification is a digitization of a 511 keV gamma ray energy signal
- the nuclear case energy detector is a LYSO crystal based scintillation detector
- the output current signal of the nuclear event energy detector is an exponential decay signal (see FIG. 1).
- the decay time constant is equal to the luminescence decay time constant of the LYSO crystal.
- the analytical expression of the output waveform after the signal passes through the active RC integration circuit (Fig. 6) is as shown in Equation 2:
- the time t after the peak time (that is, t in the following equations is the time origin of the peak time), and the signal waveform is determined by (Formula 3).
- the generation of the dynamic threshold is how to implement the signal output of (Equation 4) with an analog circuit. It is essentially two signals Multiplication: One is a linear ramp voltage signal and the other is a double exponential signal in parentheses. Comparing Equation 3 and Equation 4, it can be seen that the signal is the part of the detector output signal after the peak time.
- Equation 4 can be determined by the gain of the last stage of the amplification circuit of the dynamic threshold generation circuit. By selecting different gains, a dynamic threshold curve with different rising slopes can be obtained.
- Figure 9d is a schematic diagram of a circuit principle for generating a dynamic threshold using a dynamic threshold generation circuit in accordance with an embodiment of the present invention.
- the dynamic threshold generation circuit includes two sections, a first portion of the circuit for generating a shaped signal of the same shape as the shaped signal being measured, and a second portion for generating a linear ramp voltage signal.
- the dynamic threshold generating circuit further includes an analog multiplier for multiplying the shaped signal by the linear ramp voltage signal, and the amplifying circuit for amplifying the multiplied signal to obtain a desired Dynamic threshold signal.
- the shaped signal is multiplied by the linear ramp voltage signal from its peaking time, and the linear ramp voltage signal is zero when the shaped signal is at its peaking moment.
- the first part of the schematic includes a differential circuit, a voltage-current conversion circuit, and an integration circuit that are sequentially connected
- the second portion includes a voltage-current conversion circuit and an integration circuit that are sequentially connected. Both the first part and the second part receive an input of a square wave voltage signal.
- the square wave voltage signal is for example from an FPGA.
- FIG. 9e is a specific circuit of the circuit schematic design of the embodiment shown in FIG. 9d.
- the upper part corresponds to the first part of FIG. 9d, and the square wave voltage signal output by the FPGA is generated by the differential circuit and the nuclear case energy detector (in the In the embodiment, the LYSO scintillation crystal) outputs a voltage signal of the same shape as the current signal.
- the nuclear case energy detector in the In the embodiment, the LYSO scintillation crystal
- the voltage-current conversion circuit converts into a current signal having the same shape as the current signal output from the detector, and then passes through the same shaping circuit as the forming circuit that generates the measured signal to be formed, and produces the same shape as the shaped signal to be measured. Forming signal.
- the forming circuit is the integrating circuit.
- the lower half of Figure 9e corresponds to the second part of Figure 9d, and the square wave output by the FPGA
- the voltage signal is converted to a current signal by a voltage-current conversion circuit, and the current signal is passed through an integration circuit of the second portion to generate the linear ramp voltage signal.
- the two parts of the signal are multiplied by an analog multiplier circuit that outputs a dynamic threshold signal that satisfies the requirements.
- an analog signal delay circuit is provided prior to the generation of the linear ramp voltage signal to assist in adjusting the timing at which the two signals begin to multiply.
- an arbitrary waveform signal generator is used to generate a signal waveform identical to that of the above-described nuclear case energy detector output, and the amplitude of the waveform is controlled to obtain different measured values.
- Figure 9f shows the correspondence between the measured over-threshold time and the amplitude value of the nuclear event energy detector output signal in accordance with Figure 9e, which controls the two dynamic threshold curves produced by the gain of the last stage of amplification. It can be observed that, using the analog circuit real-time generating technical solution provided by the present invention, the threshold time and the amplitude value of the output signal of the measured forming circuit do have a good linear relationship.
- a pre-screening circuit can be used in conjunction with a constant ratio timing circuit to perform a constant ratio timing output of an active instance.
- the pre-screening circuit is a simple leading edge timing circuit that compares the input signal with a fixed leading edge timing threshold (ie, always galvanic), and only when the input signal has a value that exceeds the leading edge timing threshold. Screen the signal output.
- the purpose of introducing a pre-screening circuit is to prevent noise from causing the output of a constant-ratio timing signal.
- the leading-edge timing threshold of the pre-screening circuit is given by a digital-to-analog conversion (DAC) chip output. It should be noted that the output of the pre-screening signal precedes the output of the constant ratio timing signal.
- the output of the constant ratio timing circuit can only be output through the AND gate if the pre-screening circuit has an output.
- the outputs of the constant ratio timing circuit and the pre-screening circuit are both connected to the input end of the FPGA chip (see FIG. 13), and the necessary logic selection is implemented inside the FPGA (the logic is selected as described above). Similar to the function of the gate), the output of the constant ratio timing circuit is considered to be the arrival of a nuclear case signal only when the pre-screening signal is received, and the arrival time of the constant-ratio timing signal is recorded.
- the TDC unit 5 may further include an indication unit, and the indication unit may instruct the TDC circuit to operate when receiving the pre-screening signal.
- the basic principle of the currently used time measurement technology is to use the counter to calculate from the beginning.
- the number of clocks between the end time (in the present invention, the start time is the arrival time Td of the constant ratio timing signal, and the end time is the over-threshold time Tot), and increasing the clock frequency is a method for improving the measurement accuracy.
- the four-phase clock can be used to reduce the time measurement error to one-fourth of a clock cycle when using the FPGA to implement the TDC circuit.
- the constant ratio timing signal in Fig. 14 arrives (the timing signal in Fig. 14 represents the constant ratio timing signal and the pulse signal output from the comparator circuit described below), the calculated value of the 9-bit counter is locked together with the state of the four phases.
- the computer can finally digitize the amplitude according to subtracting ⁇ from the threshold time.
- the time measurement accuracy obtained by using the four-phase clock is one quarter of the 160 MHz clock cycle, that is, 1.5625 ns.
- the main factors affecting the final measurement accuracy are constant ratio timing accuracy, TDC measurement accuracy, accuracy of the dynamic threshold generation circuit, etc., and also depend on the correspondence between the TDC measurement range and the measured signal amplitude range, and the core is The characteristics of the case energy detector and these factors are considered together to realize the high-precision nuclear radiation pulse amplitude digital measurement.
- the present invention can be used on the basis of not significantly increasing the front-end readout circuit of the nuclear case energy detector.
- the time digital conversion technology converts the voltage amplitude value of the nuclear radiation pulse signal into a threshold time amount, thereby realizing the digitization of the energy signal.
- FPGA chip technology or semiconductor integrated circuit technology is also utilized.
- the most prominent advantage of using the above technology is that parallel readout of a large number of channels can be realized, and high integration of multi-channel core case energy signal measurement can be realized. degree. Therefore, it is especially suitable for nuclear imaging detectors (nuclear imaging detectors require simultaneous reading of nuclear case energy information detected by tens or even hundreds or tens of thousands of channels), which can satisfy the laboratory research of nuclear imaging detectors. The need for practical applications.
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Radiation (AREA)
Abstract
Description
Claims
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP12793792.8A EP2592437B1 (en) | 2011-06-02 | 2012-05-22 | Method and system for digitalizing nuclear radiation pulse amplitude |
US13/807,765 US9182500B2 (en) | 2011-06-02 | 2012-05-22 | Method and system for amplitude digitization of nuclear radiation pulses |
JP2013517000A JP5698352B2 (ja) | 2011-06-02 | 2012-05-22 | 核放射パルス幅のディジタル化方法及びシステム |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110147423.2 | 2011-06-02 | ||
CN2011101474232A CN102338880B (zh) | 2011-06-02 | 2011-06-02 | 核脉冲幅度数字化方法及系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2012163239A1 true WO2012163239A1 (zh) | 2012-12-06 |
Family
ID=45514707
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/CN2012/075885 WO2012163239A1 (zh) | 2011-06-02 | 2012-05-22 | 核辐射脉冲幅度的数字化方法及系统 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9182500B2 (zh) |
EP (1) | EP2592437B1 (zh) |
JP (1) | JP5698352B2 (zh) |
CN (1) | CN102338880B (zh) |
WO (1) | WO2012163239A1 (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015021866A (ja) * | 2013-07-19 | 2015-02-02 | 独立行政法人理化学研究所 | 放射線検出器のための信号データ処理方法、信号データ処理装置、および放射線検出システム |
CN109959962A (zh) * | 2017-12-14 | 2019-07-02 | 中国核动力研究设计院 | 基于脉冲型中子探测器信号特性的核信号发生器 |
CN113359181A (zh) * | 2021-07-01 | 2021-09-07 | 成都理工大学 | 一种新型平头尖顶脉冲成形系统及方法 |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102338880B (zh) * | 2011-06-02 | 2012-12-26 | 中国科学技术大学 | 核脉冲幅度数字化方法及系统 |
US8981756B2 (en) * | 2012-04-11 | 2015-03-17 | Rockwell Automation Technologies, Inc. | Adaptive threshold voltage for frequency input modules |
CN102830420B (zh) * | 2012-08-31 | 2014-11-26 | 中国原子能科学研究院 | 核辐射探测器工作电压标定方法及系统 |
CN102841366B (zh) * | 2012-08-31 | 2014-11-19 | 中国原子能科学研究院 | 脉冲幅度甄别器甄别阈测定方法及系统 |
CN102841368B (zh) * | 2012-08-31 | 2014-11-19 | 中国原子能科学研究院 | 气体核辐射探测器收集的电荷数与外加电压的关系曲线测量方法及系统 |
CN102916683B (zh) * | 2012-10-18 | 2016-09-14 | 成都理工大学 | 一种可调参数核脉冲模拟方法 |
CN103675891B (zh) * | 2013-12-11 | 2016-03-02 | 成都理工大学 | 基于双线性变换法的数字核脉冲高斯成形方法 |
CN104242877B (zh) * | 2014-07-16 | 2017-04-05 | 成都理工大学 | 基于模拟cr‑rc电路的核脉冲信号数字高斯成形方法 |
CN104155680B (zh) * | 2014-08-28 | 2016-05-04 | 东华理工大学 | 基于ads1605的多道能谱分析仪数据采集装置 |
CN104777508B (zh) * | 2015-03-30 | 2017-06-13 | 东南大学 | 一种基于模型库的数字脉冲重叠峰分离算法 |
CN106160706A (zh) * | 2015-04-24 | 2016-11-23 | 南京理工大学 | 一种双通道时刻鉴别电路 |
CN105141294B (zh) * | 2015-08-17 | 2018-03-30 | 西北核技术研究所 | 双预置甄别门控型恒比定时甄别器及数字恒比定时甄别方法 |
FR3058230B1 (fr) * | 2016-10-27 | 2019-03-15 | Detection Technology Sas | Dispositif de spectrometrie |
JP7057630B2 (ja) * | 2017-06-23 | 2022-04-20 | 浜松ホトニクス株式会社 | 放射線検出装置 |
RU2705933C1 (ru) * | 2018-12-26 | 2019-11-12 | Федеральное государственное автономное образовательное учреждение высшего образования "Национальный исследовательский Томский политехнический университет" | Двухканальный сцинтилляционный счетчик ионизирующего излучения |
CN110133708B (zh) * | 2019-04-24 | 2021-03-09 | 中国科学技术大学 | 一种核脉冲信号的测量装置及测量方法 |
CN110531404B (zh) * | 2019-06-03 | 2021-07-09 | 中国科学技术大学 | 核脉冲电荷时间转换方法与系统 |
CN110347965B (zh) * | 2019-07-17 | 2023-02-03 | 四川新先达测控技术有限公司 | 核脉冲信号处理方法及装置 |
CN110568468B (zh) * | 2019-08-27 | 2023-01-10 | 福州智元仪器设备有限公司 | 一种辐射脉冲计数突变算法 |
CN110954938B (zh) * | 2019-12-11 | 2021-08-31 | 湖北锐世数字医学影像科技有限公司 | 比较器的阈值校正方法、装置和采样电路 |
CN111308223B (zh) * | 2019-12-12 | 2022-05-13 | 江苏科技大学 | 基于传输响应的脉冲信号幅度测量方法 |
CN113126140A (zh) * | 2019-12-30 | 2021-07-16 | 中核控制系统工程有限公司 | 一种高速数字n/γ波形实时甄别的系统及方法 |
RU198513U1 (ru) * | 2019-12-30 | 2020-07-14 | Федеральное государственное автономное образовательное учреждение высшего образования "Новосибирский национальный исследовательский государственный университет" (Новосибирский государственный университет, НГУ) | Стенд для проведения испытаний сцинтилляционных датчиков |
CN111487663B (zh) * | 2020-03-16 | 2022-06-03 | 中国石油天然气集团有限公司 | 一种核测井仪器能量电流谱获取方法 |
JP7386577B2 (ja) | 2020-07-03 | 2023-11-27 | 国立大学法人 東京大学 | 検出回路 |
CN112255662A (zh) * | 2020-10-14 | 2021-01-22 | 南昌华亮光电有限责任公司 | 基于智能时间采样的伽马相机及探测方法 |
CN112327347B (zh) * | 2020-10-29 | 2022-11-18 | 中广核久源(成都)科技有限公司 | 曲率可调的数字化核脉冲成形系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006300551A (ja) * | 2005-04-15 | 2006-11-02 | Sumitomo Heavy Ind Ltd | 放射線検出回路、放射線検出器および放射線検査装置 |
CN101034120A (zh) * | 2007-03-09 | 2007-09-12 | 中国科学院上海光学精密机械研究所 | 脉冲形状的测量装置及测量方法 |
US20100084559A1 (en) * | 2008-09-11 | 2010-04-08 | Siemens Medical Solutions Usa, Inc. | Use of crystal location in nuclear imaging apparatus to minimize timing degradation in a photodetector array |
CN101852861A (zh) * | 2010-04-09 | 2010-10-06 | 罗康林 | 基于峰值侦测的多道脉冲幅度分析仪 |
CN102338880A (zh) * | 2011-06-02 | 2012-02-01 | 中国科学技术大学 | 核脉冲幅度数字化方法及系统 |
Family Cites Families (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4217497A (en) * | 1978-03-14 | 1980-08-12 | The United States Of America As Represented By The Department Of Health, Education And Welfare | Portable instrument for measuring neutron energy spectra and neutron dose in a mixed n-γ field |
JPS6085382A (ja) * | 1983-10-17 | 1985-05-14 | Hitachi Ltd | 広レンジ信号処理方法 |
DE4223773C2 (de) * | 1992-07-18 | 1997-07-10 | Berthold Lab Prof Dr | Verfahren zur Unterscheidung und gleichzeitigen oder getrennten Messung von Einzel- und Mehrelektronenereignissen in einem optoelektronischen Detektor |
JPH10197643A (ja) * | 1997-01-06 | 1998-07-31 | Toshiba Corp | 放射線計測装置及び放射線計測システム |
JP3103047B2 (ja) * | 1997-09-01 | 2000-10-23 | シュランベルジェ、ホールディング、リミテッド | 核分光システムにおけるパルス波形調整及び弁別方法 |
US6252232B1 (en) * | 1998-06-29 | 2001-06-26 | General Electric Company | Digital integrator |
US6160259A (en) * | 1998-10-23 | 2000-12-12 | Adac Laboratories | Channel-specific control of pulse integration in a gamma camera system in response to pulse pile-up |
US6232604B1 (en) * | 1998-11-06 | 2001-05-15 | General Electric Company | Analog time adjustment for coincidence detection electronics |
JP4064009B2 (ja) * | 1999-07-30 | 2008-03-19 | 株式会社東芝 | 線種弁別型放射線検出装置 |
US6459086B1 (en) * | 2000-11-24 | 2002-10-01 | Koninklijke Philips Electronics, N.V. | Digital peak detector for radiation detection systems |
JP4160275B2 (ja) * | 2001-05-28 | 2008-10-01 | 浜松ホトニクス株式会社 | エネルギー測定方法及び測定装置 |
US6781426B2 (en) * | 2001-07-26 | 2004-08-24 | The Regents Of The University Of California | Nanosecond monolithic CMOS readout cell |
US7818047B2 (en) * | 2001-11-09 | 2010-10-19 | Nova R&D, Inc. | X-ray and gamma ray detector readout system |
US6664543B2 (en) * | 2002-01-28 | 2003-12-16 | Cti Pet Systems, Inc. | Continuous sampling and digital integration for PET scintillation |
US6903344B2 (en) * | 2003-03-25 | 2005-06-07 | Cti Pet Systems, Inc. | Baseline correction in PET utilizing continuous sampling ADCs to compensate for DC and count rate errors |
US6822506B2 (en) * | 2003-04-17 | 2004-11-23 | Concorde Microsystems, Inc. | Amplitude and rise-time sensitive timing-shaping filters with built-in pulse-tail cancellation for high count-rate operation |
US7126126B2 (en) * | 2003-10-16 | 2006-10-24 | Brookhaven Science Associates, Llc | Compact conscious animal positron emission tomography scanner |
US7045802B2 (en) * | 2003-11-26 | 2006-05-16 | General Electric Company | Method and apparatus for coincidence imaging digital triggering |
JP4619674B2 (ja) * | 2004-03-24 | 2011-01-26 | 株式会社東芝 | 磁気共鳴イメージング装置 |
TWI313758B (en) * | 2004-04-02 | 2009-08-21 | Nat Health Research Institutes | A device and method for characterizing energy,peak time and decay time constant of gamma ray |
JP4611106B2 (ja) * | 2005-03-11 | 2011-01-12 | 住友重機械工業株式会社 | 放射線検出回路及び放射線検査装置 |
EP2076790B1 (en) * | 2006-10-04 | 2014-02-12 | CERN - European Organization For Nuclear Research | Readout circuit for use in a combined pet-ct apparatus |
US7800439B2 (en) * | 2006-10-27 | 2010-09-21 | Ut-Battelle, Llc | High speed preamplifier circuit, detection electronics, and radiation detection systems therefrom |
US7605375B2 (en) * | 2007-04-26 | 2009-10-20 | Oy Ajat Ltd. | Multi-functional radiation/photon identifying and processing application specific integrated circuit and device |
EP2183578B1 (en) * | 2007-08-03 | 2018-08-29 | Judith B. Mott Revocable Trust Dated June 6, 2006, As Amended And Restated | Pileup rejection in an energy-dispersive radiation spectrometry system |
US7683334B2 (en) * | 2007-08-07 | 2010-03-23 | The State Of Oregon Acting By And Through The State Board Of Higher Education On Behalf Of Oregon State University | Simultaneous beta and gamma spectroscopy |
US8340377B2 (en) * | 2007-09-17 | 2012-12-25 | Siemens Medical Solutions Usa, Inc. | Method for energy calculation and pileup determination for continuously sampled nuclear pulse processing |
JP2009079969A (ja) * | 2007-09-26 | 2009-04-16 | Toshiba Corp | 放射線スペクトル計測システム |
KR101583272B1 (ko) * | 2007-11-02 | 2016-01-07 | 유니버시티 오브 워싱톤 스루 이츠 센터 포 커머셜리제이션 | 양전자 방출 단층촬영을 위한 데이터 취득 |
US8258480B2 (en) * | 2008-03-03 | 2012-09-04 | The Board Of Trustees Of The Leland Stanford Junior University | High energy photon detection using pulse width modulation |
US8492705B2 (en) * | 2009-10-23 | 2013-07-23 | Saint-Gobain Ceramics & Plastics, Inc. | Radiation detector and method of using a radiation detector |
US9223031B2 (en) * | 2012-01-04 | 2015-12-29 | Kabushiki Kaisha Toshiba | Method for timing-pick-off of undersampled pulses from radiation detectors using a photosensor-based radiation detection system applying a prototype waveform and a weighting function |
-
2011
- 2011-06-02 CN CN2011101474232A patent/CN102338880B/zh not_active Expired - Fee Related
-
2012
- 2012-05-22 WO PCT/CN2012/075885 patent/WO2012163239A1/zh active Application Filing
- 2012-05-22 EP EP12793792.8A patent/EP2592437B1/en active Active
- 2012-05-22 JP JP2013517000A patent/JP5698352B2/ja not_active Expired - Fee Related
- 2012-05-22 US US13/807,765 patent/US9182500B2/en not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006300551A (ja) * | 2005-04-15 | 2006-11-02 | Sumitomo Heavy Ind Ltd | 放射線検出回路、放射線検出器および放射線検査装置 |
CN101034120A (zh) * | 2007-03-09 | 2007-09-12 | 中国科学院上海光学精密机械研究所 | 脉冲形状的测量装置及测量方法 |
US20100084559A1 (en) * | 2008-09-11 | 2010-04-08 | Siemens Medical Solutions Usa, Inc. | Use of crystal location in nuclear imaging apparatus to minimize timing degradation in a photodetector array |
CN101852861A (zh) * | 2010-04-09 | 2010-10-06 | 罗康林 | 基于峰值侦测的多道脉冲幅度分析仪 |
CN102338880A (zh) * | 2011-06-02 | 2012-02-01 | 中国科学技术大学 | 核脉冲幅度数字化方法及系统 |
Non-Patent Citations (1)
Title |
---|
SUN, JIAN ET AL.: "Development of Digital Pulse Processing Spectrometer for Online Analysis of Nuclear Physics Detectors", PROCEEDINGS OF THE 15TH NATIONAL ACADEMIC ANNUAL CONFERENCE ON NUCLEAR ELECTRONICS & DETECTION TECHNOLOGY, 2010, pages 128 - 132, XP008171067 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015021866A (ja) * | 2013-07-19 | 2015-02-02 | 独立行政法人理化学研究所 | 放射線検出器のための信号データ処理方法、信号データ処理装置、および放射線検出システム |
CN109959962A (zh) * | 2017-12-14 | 2019-07-02 | 中国核动力研究设计院 | 基于脉冲型中子探测器信号特性的核信号发生器 |
CN109959962B (zh) * | 2017-12-14 | 2022-07-26 | 中国核动力研究设计院 | 基于脉冲型中子探测器信号特性的核信号发生器 |
CN113359181A (zh) * | 2021-07-01 | 2021-09-07 | 成都理工大学 | 一种新型平头尖顶脉冲成形系统及方法 |
CN113359181B (zh) * | 2021-07-01 | 2022-11-04 | 成都理工大学 | 一种新型平头尖顶脉冲成形系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
US9182500B2 (en) | 2015-11-10 |
JP5698352B2 (ja) | 2015-04-08 |
EP2592437A1 (en) | 2013-05-15 |
EP2592437A4 (en) | 2015-02-11 |
EP2592437B1 (en) | 2019-07-10 |
JP2013534629A (ja) | 2013-09-05 |
CN102338880A (zh) | 2012-02-01 |
US20130161528A1 (en) | 2013-06-27 |
CN102338880B (zh) | 2012-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2012163239A1 (zh) | 核辐射脉冲幅度的数字化方法及系统 | |
US10817588B2 (en) | Ratio-reference measurement of the arrival time and/or amplitude of a digitized electronic pulse | |
JP6096408B2 (ja) | ガンマ線検出システム及びガンマ線検出方法 | |
Deng et al. | Development of an eight-channel time-based readout ASIC for PET applications | |
JP2013534629A5 (zh) | ||
Feinstein et al. | The analogue ring sampler: A front-end chip for ANTARES | |
Ugur et al. | Field programmable gate array based data digitisation with commercial elements | |
Di Lorenzo et al. | PARISROC, an autonomous front-end ASIC for triggerless acquisition in next generation neutrino experiments | |
Kong et al. | An FPGA-based fast linear discharge readout scheme enabling simultaneous time and energy measurements for TOF-PET detectors | |
US9696432B2 (en) | Device for processing a signal delivered by a radiation detector | |
US11635531B2 (en) | Apparatus for measuring photon information and photon measurement device | |
Fang et al. | IMOTEPAD: A mixed-signal 64-channel front-end ASIC for small-animal PET imaging | |
US10473521B2 (en) | Analog photon counting | |
Breton et al. | Picosecond time measurement using ultra fast analog memories | |
Ferry et al. | SCOTT: A time and amplitude digitizer ASIC for PMT signal processing | |
Yonggang et al. | A novel nuclear pulse digitizing scheme using time over dynamic threshold | |
Sun et al. | Energy and timing measurement with time-based detector readout for PET applications: Principle and validation with discrete circuit components | |
Calò et al. | A CMOS front-end for timing and charge readout of silicon photomultipliers | |
CN112946726B (zh) | 基于fpga-adc模块输出信号的非线性修正方法、系统 | |
CN112649837B (zh) | 粒子识别方法、电子学读出系统和中子探测器 | |
Rivetti et al. | Experimental results from a pixel front-end for the NA62 experiment with on pixel constant fraction discriminator and 100 ps time to digital converter | |
Ronconi et al. | Timestamp and Amplitude Measurement Solution for Radiation Detectors | |
JP6969734B2 (ja) | 電荷検出回路及びそれを含む放射線検出装置 | |
Lai | The Panda Strip Asic: Pasta | |
Zhang et al. | The development of a readout ASIC for TPC based X-ray polarimeter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
ENP | Entry into the national phase |
Ref document number: 2013517000 Country of ref document: JP Kind code of ref document: A |
|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 12793792 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2012793792 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 13807765 Country of ref document: US |
|
NENP | Non-entry into the national phase |
Ref country code: DE |