WO2015128905A1 - 波形弁別装置、波形弁別方法及び波形弁別プログラム - Google Patents

波形弁別装置、波形弁別方法及び波形弁別プログラム Download PDF

Info

Publication number
WO2015128905A1
WO2015128905A1 PCT/JP2014/001106 JP2014001106W WO2015128905A1 WO 2015128905 A1 WO2015128905 A1 WO 2015128905A1 JP 2014001106 W JP2014001106 W JP 2014001106W WO 2015128905 A1 WO2015128905 A1 WO 2015128905A1
Authority
WO
WIPO (PCT)
Prior art keywords
waveform
discrimination
coordinate
coordinate axis
circuit
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/JP2014/001106
Other languages
English (en)
French (fr)
Inventor
青木　徹
昭史 小池
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shizuoka University NUC
Anseen Inc
Original Assignee
Shizuoka University NUC
Anseen Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shizuoka University NUC, Anseen Inc filed Critical Shizuoka University NUC
Priority to PCT/JP2014/001106 priority Critical patent/WO2015128905A1/ja
Priority to JP2016504864A priority patent/JP6414830B2/ja
Priority to US15/121,645 priority patent/US20160356897A1/en
Publication of WO2015128905A1 publication Critical patent/WO2015128905A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

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/20—Measuring radiation intensity with scintillation detectors
    • G01T1/208—Circuits specially adapted for scintillation detectors, e.g. for the photo-multiplier section
    • 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/20—Measuring radiation intensity with scintillation detectors
    • G01T1/202—Measuring radiation intensity with scintillation detectors the detector being a crystal

Definitions

  • the present invention relates to a waveform discrimination apparatus for discriminating two pulse waveforms having different waveforms from each other, and in particular, two pulse waveforms having different rising characteristics and falling characteristics such as, for example, gamma rays and neutron rays incident on a scintillator. And a waveform discrimination method and a waveform discrimination program for discriminating an electrical signal caused by a physical quantity of.
  • a pulse having an arbitrary intensity having a first waveform similar to each other is generated at an arbitrary timing to form a group
  • a pulse of an arbitrary intensity, which is different from the first waveform but has a second waveform similar to each other but which is different from the first waveform may be generated at any timing to form a group of sets. That is, an optical signal of the first waveform resulting from the input of a gamma ray as an output of the scintillator is generated at an arbitrary timing and at an arbitrary timing to form a set of a group of optical signals.
  • the input photodetector sequentially outputs a first electric signal corresponding to the first waveform, and an optical signal of a second waveform resulting from the input of a neutron beam is generated at an arbitrary timing and at an arbitrary timing.
  • the first waveform in a situation (environment) in which a group of optical signals are composed and a photodetector receiving this group of optical signals sequentially outputs a second electrical signal corresponding to the second waveform In some cases, it may be desirable to distinguish a second waveform different from the first waveform.
  • the output of the waveform shaping amplifier is branched into a peak value observation output and a rise time observation output, and the peak value observation output is directly input to one input terminal of two-input MCA.
  • the rise time observation output is input to the pulse waveform analyzer.
  • the pulse waveform analyzer outputs two signals at timings of 10% and 90% of the rise time, inputs the two signals to the time / amplitude converter, and the time / amplitude converter generates a time difference between the two signals. Is converted to the amplitude of the pulse, and the output of the time / amplitude converter is input to the other input terminal of the two-input MCA, using a very complicated, large and expensive system configuration.
  • Non-Patent Document 1 As shown in FIG. 3 of Non-Patent Document 1, plotting rise time (Rise Time) of gamma ray and neutron beam and its pulse height (Pulse Height) on a coordinate plane. The gamma ray and the neutron ray are separated and displayed. However, in the invention described in Non-Patent Document 1, as shown in FIG. 3 of Non-Patent Document 1, both neutrons and gamma rays simultaneously emitted from calcium 252 ( 252 Cf) as a radiation source are used. The waveforms of (1) overlap in both the peak value axis and the time axis.
  • Non-Patent Document 1 The reason is that the delay due to the transient characteristic (slew rate) of the charge sensitive pre-amplifier shown in FIG. 1 of Non-Patent Document 1 causes the high-speed input signal with a large peak value to have a large rise time. It is for. Therefore, in the invention described in Non-Patent Document 1, when the energy of the incident gamma ray is higher than the energy of the neutron, the discrimination becomes impossible and a counting error occurs.
  • Non-Patent Document 1 in the two-dimensional distribution diagram described in FIG. Within the range of ROI) A, a rectangular counting region of interest B for extracting and suppressing gamma rays which are not to be measured is set.
  • manual adjustment is performed while setting the area of counting interest area B in the coordinate plane that is within the area of counting interest area A while observing the data plot on the two-dimensional distribution map
  • human attention and devise are needed.
  • the plot trajectory of the gamma ray is non-linear, and it is difficult to accurately carry out the separation from the plot by neutrons.
  • Non-Patent Document 1 does not consider the correspondence to the case where the input signal by gamma rays and neutron rays is incident within a short time interval, and gamma rays within a time interval shorter than the time constant of the pulse waveform analyzer. Because the pile-up occurs when there is an input signal by neutron beam, accurate energy can not be measured, and the accuracy decreases.
  • the present inventors have, as an example, a pulse of a physical quantity with an arbitrary intensity having a first waveform, a current waveform by light emission when a gamma ray enters the scintillator, and a second waveform. Attention was focused on the current waveform due to light emission when neutrons entered the scintillator as a pulse of a physical quantity of arbitrary intensity.
  • a waveform detector for inputting a waveform of a measured pulse and converting a physical quantity of the measured pulse into an electrical signal, and (b) a transient response waveform of the electrical signal on a time axis
  • An analog amplifier that expands and amplifies along
  • an AD converter that samples the amplified electrical signal in the rising and falling periods of the electrical signal and converts it into digital data
  • digital data The feature quantity of the rising period is calculated as a point on the first coordinate axis
  • the feature quantity of the falling period is calculated as a point on the second coordinate axis, and further, the point on the first coordinate axis and the second coordinate axis
  • the gist of the present invention is a waveform discrimination apparatus including a signal processing circuit that plots coordinate points on a discriminant plane defined by a first coordinate axis and a second coordinate axis with a set of points as coordinate points.
  • a second aspect of the present invention comprising the steps of (a) inputting the waveform of the pulse to be measured and converting the physical quantity of the pulse to be measured into an electrical signal; The step of expanding and amplifying, (c) sampling the amplified electric signal in the rising period and falling period of the electric signal and converting it into digital data, and (d) using the digital data, the rising period Calculating the feature quantity of the second coordinate axis as a point on the first coordinate axis, and calculating the feature quantity of the falling period as a point on the second coordinate axis; (e) the points on the first coordinate axis and the points on the second coordinate axis The step of plotting the coordinate point on the discriminant plane defined by the first coordinate axis and the second coordinate axis with the set as the coordinate point, and (f) the plot position of the coordinate point, the measured pulse has a first waveform or a first waveform Is a second waveform different from And summarized in that
  • a computer software program for realizing the waveform discrimination method described in the second aspect of the present invention is stored in a computer readable recording medium, and the recording medium is read by a computer system to realize the waveform discrimination of the present invention.
  • the method can be implemented. That is, according to the third aspect of the present invention, (a) a command for causing the waveform of the pulse to be measured to be input to the waveform detector and converting the physical quantity of the pulse to be measured into an electrical signal; A command to expand and amplify the transient response waveform along the time axis, and (c) a command to sample the amplified electric signal in the rising and falling periods of the electric signal and convert it into digital data (D)
  • the difference value calculation circuit, the attenuation amount calculation circuit, and the difference value integration circuit of the signal processing circuit cooperate with each other to calculate the feature value of the rising period as a point on the first coordinate axis using digital data
  • E a two-dimensional coordinate plotting circuit of the signal processing circuit, a set of points
  • the "recording medium” means, for example, a medium capable of recording programs such as an external memory device of a computer, a semiconductor memory, a magnetic disk, an optical disk, an optical magnetic disk, a magnetic tape and the like. Specifically, flexible disks, CD-ROMs, MO disks, cassette tapes, open reel tapes and the like are included in the "recording medium”.
  • the waveform discrimination apparatus according to the first aspect is capable of reducing the size of the apparatus, and in the case of design for miniaturization, the processor embedded in an apparatus such as a microcontroller unit (MCU), the third aspect
  • MCU microcontroller unit
  • the MCU was initially programmed only in assembly language due to the small amount of memory, but when the amount of memory and processing power of the CPU improved, C language was used from the viewpoint of development efficiency.
  • a language processing system such as a BASIC language interpreter is written in ROM in advance, and it is possible to realize a recording medium or the like storing the waveform discrimination program according to the third aspect.
  • FIG. 1 It is a typical block diagram explaining an outline of an important section of a waveform distinction apparatus concerning a 1st embodiment of the present invention. It is a figure which illustrates more concretely an example of the analog amplifier used for the waveform discrimination device concerning a 1st embodiment. It is a three-dimensional drawing which makes the side wall of a housing
  • FIG. 4A when the output of the photomultiplier tube as a photodetector is terminated at 50 ⁇ , when the gamma ray enters the scintillator, the light of the first waveform emitted from the scintillator is transmitted to the photomultiplier tube.
  • FIG. 4B is a diagram showing a pulse waveform obtained by observing the first electric signal shown in FIG. 4A with an oscilloscope when the output of the photomultiplier tube is terminated at 50 k ⁇ .
  • FIG. 6 (a) is an enlarged view of a time axis of the V portion of FIG. 4 (a) as in FIG. 5, and FIG. 6 (b) is a pulse of the second electric signal for comparison. It is a figure which shows a waveform by making the time-axis into common with Fig.6 (a). 7 (a) shows the output waveform of the analog amplifier of the waveform discrimination apparatus according to the first embodiment, and FIG.
  • FIG. 7 (b) shows that the AD converter can digitize the waveform of FIG. 7 (a).
  • the correct wave is generated even when pileup occurs in which the rising waveform of the next pulse signal is superimposed on the falling waveform of the electric signal as the output of the light detector. It is a figure explaining that high value can be acquired, and it is a mimetic diagram notionally explaining the state of processing in an AD converter. It is a conceptual flowchart explaining the flow of the processing using the window part for discrimination in the waveform discrimination method concerning a 1st embodiment, and the straight line showing the linear equation for discrimination. It is a conceptual flowchart explaining the determination method of the window part for discrimination used for the waveform discrimination method which concerns on 1st Embodiment, and the straight line which represents the linear equation for discrimination.
  • the waveform discrimination apparatus receives a waveform of a pulse to be measured, and converts the physical quantity of the pulse to be measured into an electrical signal;
  • An analog amplifier 13 connected to the detector 12 for amplifying the transient response waveform of the electric signal along the time axis for amplification and connected to the analog amplifier 13 and amplified during the rise and fall periods of the electric signal It is connected to AD converter 14 which samples electric signal and converts it into digital data, and AD converter 14, and using digital data, the feature quantity of the rising period is calculated as a point on the first coordinate axis, falling period
  • a signal processing circuit 15 for plotting coordinate points, a display
  • FIG. 1 illustrates the data storage device 18 and the program storage device 19 as if they were individual hardware resources for the sake of convenience, data storage devices may be used as actual physical hardware resources. It does not deny that 18 and program storage devices 19 are each configured as a set of a plurality of storage devices having different functions and sizes.
  • the data storage device 18 can be any combination appropriately selected from a group including a plurality of registers, a plurality of cache memories, a main storage device, and an auxiliary storage device.
  • the cache memory may be a combination of a primary cache memory and a secondary cache memory, and may further have a hierarchy including a tertiary cache memory.
  • the waveform detector 12 shown in FIG. 1 has at least one pulse group included in a first group of pulses forming a group by generating pulses of any intensity having a first waveform as a measured pulse at an arbitrary timing.
  • a pulse is input, and a first electrical signal corresponding to a first waveform is sequentially output, and as another measured pulse, a pulse of any intensity having a second waveform different from the first waveform is generated at an arbitrary timing.
  • at least one pulse included in the second pulse group forming the group of groups is input, and the second electric signal corresponding to the second waveform is sequentially output.
  • the analog amplifier 13 receives at least one of the first and second electrical signals as a discrimination target signal from the waveform detector 12 and the discrimination target signal so as to expand the transient response waveform of the discrimination target signal along the time axis.
  • Amplify the waveform of The analog amplifier 13 has a waveform showing a transient response characteristic so that the fall time will be about 2 ⁇ sec or more even if the first waveform as the pulse to be measured is a nanosecond level half width pulse. It is preferable to expand along the time axis.
  • the sampling interval for the AD converter 14 to acquire digital data can be lengthened, so that a very inexpensive and simple AD converter 14 can be adopted.
  • the AD converter 14 samples the amplified discrimination target signal in the rising period and the falling period of the discrimination target signal, generates discrete data at regular intervals, and converts the discrete data into digital data. Do.
  • the data acquisition circuit 162 (see FIG. 13) of the signal processing circuit 15 sequentially reads discrete digital data sampled and generated by the AD converter 14 from the AD converter 14 and temporarily stores it in the data storage unit 18.
  • the signal processing circuit 15 reads the digital data stored in the data storage device 18, calculates the feature amount Df of the falling period as a point on the second coordinate axis, and sets the feature amount Uf of the rising period to the first coordinate axis Calculate as a point on the orthogonal first coordinate axis. Then, with the set of points on the first coordinate axis and the set of points on the second coordinate axis (Uf, Df) as coordinate points, as shown in FIG.
  • the signal processing circuit 15 AD converts the coordinate points (Uf, Df) As the processing of sampling of the transient response waveform by the unit 14 progresses, it automatically plots in real time on the discrimination plane defined by the first coordinate axis and the second coordinate axis. Since the coordinate points (Uf, Df) are plotted in real time on the discrimination plane shown in FIG. 11, the data storage device 18 functions as a register for temporarily storing digital data output from the AD converter 14 .
  • the second coordinate axis is in the X-axis direction, but the value of the feature amount Df in the falling period is plotted as a value that increases toward the right in the X-axis direction.
  • the first coordinate axis is displayed as the Y axis, but the value of the feature amount Uf in the rising period becomes smaller as it goes upward in the Y axis direction, and becomes lower in the Y axis direction. Since it is plotted as a value that increases as it progresses, the discrimination plane defined by the first coordinate axis and the second coordinate axis is defined as Cartesian coordinates of the third quadrant.
  • the first coordinate axis is displayed as the Y axis and the second coordinate axis is the X axis, but this is merely an example, and either X axis or Y axis is selected as the first coordinate axis, and the remaining one is left. It may be determined as appropriate which of the two coordinate axes is to be selected.
  • the logical hardware resource configuration of the signal processing circuit 15 is expressed as shown in FIG.
  • a microprocessor (MPU) or the like mounted as a microchip can be used.
  • a digital signal processor (DSP) specialized in signal processing with an enhanced arithmetic operation function a microcontroller (microcomputer) mounted with a memory and peripheral circuits and intended for embedded device control, etc. are used. May be Alternatively, the main CPU of the current general-purpose computer may be used as the signal processing circuit 15.
  • the waveform discrimination apparatus determines whether the discrimination target signal uses the first waveform as the generation source or the second waveform as the generation source from the distribution position of the coordinate point on the discrimination plane shown in FIG.
  • the configuration shown in FIGS. 2 and 3 relates to a specific application example of the waveform discrimination apparatus according to the first embodiment, in which the first waveform as the pulse to be measured is from the radiation-light conversion element 11 specific to gamma rays.
  • the light emission waveform and the second waveform as another measured pulse are light emission waveforms from the radiation-light conversion element 11 inherent to the neutron beam.
  • the first waveform is a light emission waveform from the radiation-to-light converter 11 caused by the gamma ray being input to the radiation-to-light converter 11
  • the second waveform is a neutron beam is a radiation-to-light
  • the light emission waveform from the radiation-light conversion element 11 resulting from the input to the conversion element 11 will be exemplarily described. That is, as shown in FIGS. 2 and 3, the waveform discrimination apparatus according to the first embodiment is connected to a radiation-light conversion element 11 for converting neutron rays and gamma rays into light, and to the radiation-light conversion element 11. And a photodetector 12a for converting the light emitted from the radiation-light conversion element 11 into an electric signal.
  • CsLiYCl, LiCaAlF 6 , LiF / ZnS, LiBaF 3 , Li 6 Gd as shown in Table 1 are radiation-light conversion elements 11 for converting neutron beams and gamma rays into light with different physical quantities of transient response characteristics.
  • BO 3 It is possible to use 3 mag.
  • elements serving as light emission centers are CsLiYCl, LiCaAlF 6 , LiF / ZnS, LiBaF 3 , Li 6
  • a scintillator material such as Gd (BO 3 ) 3
  • the transient response characteristics of light emission from the radiation-light conversion device 11 are gamma rays (first waveform) and neutrons. The rising characteristics and the falling characteristics of the line (second waveform) are different.
  • light emission with gamma rays is light emission with a very short light emission time of about several nanoseconds and broad light emission following a sharp peak at the tip.
  • the light emission of the neutron beam is characterized by light emission having a relatively long light emission time of about several hundreds nanoseconds or less, as shown in FIG. It is said that light emission at gamma rays (first waveform) is Cerenkov light emission.
  • scintillators such as CsLiYCl, LiCaAlF 6 , LiF / ZnS, LiBaF 3 , and Li 6 Gd (BO 3 ) 3 emit light with a wavelength of about 190 to 450 nm.
  • Photodetector 12a that converts light emitted by the sensor into an electrical signal, such as a photomultiplier tube (PMT) capable of converting light with a wavelength of about 190 to 450 nm, an electrical signal, a semiconductor photodiode, a photodiode array, Geiger mode parallel readout An APD pixel array or the like can be used.
  • PMT photomultiplier tube
  • the photodetector 12a receives a pulse of a first waveform and outputs a first electric signal corresponding to the first waveform, and receives a pulse of a second waveform and outputs a second electric signal corresponding to the second waveform.
  • the characteristic of outputting a signal is required, and the device performance that maintains linearity between the input and the output of the photodetector 12a is ideal.
  • the radiation-light conversion element 11 is attached to the window of the light detector 12 a, and the light detector 12 a is configured as an upper protruding portion of the housing 21.
  • the inside of the housing 21 is mounted on the circuit board 23 and the circuit board 23 connected to the output of the light detector 12a via the cables 31a and 31b, and through the embedded wiring or the surface wiring in the circuit board 23 ,
  • the analog amplifier 13 and the high voltage power supply 22 electrically connected to the cables 31a and 31b, the circuit board 24 connected to the circuit board 23 via the cables 32a, 32b and 32c, and the circuit board 24
  • the AD converter 14 and the signal processing circuit 15 electrically connected to the cables 32 a, 32 b and 32 c via embedded wiring or surface wiring in the circuit board 24 are incorporated.
  • the case 21 exemplarily has a rectangular parallelepiped shape, and the display device 16 is attached to the top of one side of the case 21.
  • the shape of the case 21 is not limited to the shape of a rectangular parallelepiped, and may be another shape such as a cylindrical shape. In the case of the cylindrical case 21, a part of the circumferential surface of the case 21 is flat Alternatively, the display device 16 may be attached in a structure in which the display device 16 is embedded or in a topology in which a portion protrudes from the circumferential surface of the housing 21.
  • Adjustment knobs 34 a, 34 b, 34 c and 34 d for setting the conditions of the signal processing circuit 15 are provided on the bottom surface of the housing 21.
  • a hole is provided on the bottom surface of the housing 21, and the communication cable 33 connected to the signal processing circuit 15 through the embedded wiring or the surface wiring in the circuit board 24 is conducted to the outside of the housing 21 from this hole. It is
  • the system configuration including the data storage device 18 and the program storage device 19 in FIG. 1 is merely an example, and the data storage device 18 includes the AD converter 14 and signals shown in FIG. It may be present as an internal structure of the processing circuit 15 or the like, and functions of a part of storage devices are distributed as a register or the like as an internal memory of the signal processing circuit 15 and the remaining functions are mounted on the circuit board 24 As a physical structure, the data storage devices 18 may be distributed so as to be executed in the external memory. Alternatively, a physical configuration may be employed in which only the external memory mounted on the circuit board 24 is provided as the data storage device 18, and the data storage disposed outside the housing 21 connected via the communication cable 33 A device 18 may be included.
  • the program storage device 19 of FIG. 1 may also exist as an internal structure of the signal processing circuit 15, the control circuit 17 and the like, and a storage device as an internal memory of the signal processing circuit 15 and the control circuit 17 It may be configured to include both of the storage devices as the memory, or may be a storage device existing as a configuration of only the external memory.
  • the control circuit 17 may be present as an internal structure of the housing 21 as an actual physical structure by dispersing at least a part of the functions of the control circuit 17 of FIG. In some cases, the control circuit 17 may be configured as the internal structure of the signal processing circuit 15, or conversely, the signal processing circuit 15 and the AD converter 14 may be integrated in the internal structure of the control circuit 17 in the form of functional blocks. Such various physical structures may be realized.
  • the waveform discrimination apparatus has a simple configuration, so that the apparatus size can be reduced, and an apparatus such as a microcontroller unit (MCU) can be incorporated.
  • MCU microcontroller unit
  • the MCU incorporates a computer system including the analog amplifier 13, the AD converter 14, the signal processing circuit 15, the control circuit 17 and the like shown in FIG. 1 into one integrated circuit.
  • the MCU is a type of microprocessor that emphasizes self-sufficiency and low cost, it is possible to function as a computer with only one semiconductor chip. If it is configured with an MCU, the number of peripheral parts can be reduced compared to a general-purpose CPU, so it is easy to make the waveform discrimination apparatus according to the first embodiment compact.
  • the input terminal I of the analog amplifier 13 is at the output side of the light detector 12a, and the signal output terminal and the reference potential point terminal of the light detector 12a are respectively the input terminal I of the analog amplifier 13 and the ground It is connected in such a configuration as to be connected between.
  • the input terminal I of the analog amplifier 13 is connected to the non-inverted input terminal of the first operational amplifier U1 constituting the input stage of the analog amplifier 13, and the input resistor R1 is connected between the input terminal I of the analog amplifier 13 and the ground.
  • the inverting input terminal of the first operational amplifier U1 is further connected to the output terminal of the first operational amplifier U1 via a feedback resistor R5.
  • the output terminal of the first operational amplifier U1 is further connected to the inverting input terminal of the second operational amplifier U2 via the transmission resistor R2, and the inverting input terminal of the second operational amplifier U2 is the output of the second operational amplifier U2 via the feedback resistor R3. It is connected to the terminal.
  • the output terminal of the second operational amplifier U2 is further connected to the noninverting input terminal of the third operational amplifier U3 constituting the output stage of the analog amplifier 13 via the transmission resistor R4, and the inverting input terminal of the third operational amplifier U2 is the third operational amplifier
  • the output terminal of the third operational amplifier U3 constitutes the output terminal O of the analog amplifier 13, which is directly connected to the output terminal of U3.
  • the attenuation time constant ⁇ R1 due to the capacitance Cp between the output terminals of the photodetector 12a and the input resistance R1 of the analog amplifier 13.
  • the transient response waveform can be expanded along the time axis by setting the value of Cp as a large value, and the high frequency components of the second waveform (neutron beam) and the first waveform (gamma ray) that are originally intended to be discriminated can be shifted to a low frequency band .
  • the transient response waveform is expanded in the time axis direction by about 1000 times as shown in FIG. 4 (b) as compared to the case where the output is terminated at 50 ⁇ . I understand.
  • FIG. 5 is a diagram showing the pulse waveform shown in the V portion of the first electric signal of FIG. 4A with an enlarged time axis, but for comparison, the pulse waveform of the second electric signal is also shown in FIG.
  • the axes are shown in common.
  • FIG. 5 shows the case where the Ce-added LiCaAlF 6 crystal is used as the radiation-light conversion element (scintillator) 11, and the light emission from the radiation-light conversion element 11 is detected by the photodetector (photoelectron multiplier tube) 12a.
  • Fig. 6 shows pulse waveforms of first and second electrical signals.
  • 6 (a) is a diagram showing the pulse waveform shown in the V portion of the first electric signal of FIG.
  • FIG. 6 (a) shows the pulse waveform of the second electrical signal with the time axis in common with FIG. 6 (a) for comparison with FIG. 6 (a). It can be seen that the signal does not have a steeple portion exhibiting a steep rise / fall characteristic as shown in FIG. 6 (a).
  • the amplifier 13 expands the transient response waveform of the first electric signal of FIG. 4A along the time axis such that the fall time is about 2 ⁇ sec or more as shown in FIG. 7A. That is, the analog amplifier 13 converts at least one of the first and second electrical signals into a signal that can be digitized by the general-purpose AD converter 14 as shown in FIG. 7B.
  • the signal processing circuit 15 connected to the AD converter 14 is configured as shown in FIGS.
  • the signal processing circuit 15 of the waveform discrimination apparatus is used for calibration prior to executing the processing shown in FIGS. 9 and 10.
  • the window boundary condition determination circuit 151 inputs the waveform of the above into the waveform detector 12 to determine the necessary window boundary conditions for discriminating the waveform, and similarly, prior to executing the processing shown in FIG. 9 and FIG.
  • a linear equation determination circuit 152 which inputs a waveform for calibration to the waveform detector 12 and determines a necessary linear equation for discrimination of the waveform, and continuous two of the discrete data of the fixed interval generated by the AD converter 14
  • a difference value integration circuit 153 for calculating the difference between the values and an attenuation amount calculation circuit for calculating the attenuation from the peak value of the waveform realized in the rising period of the discrimination target signal in the falling period of the discrimination target signal 54
  • a difference value integration circuit 155 which integrates the difference values output from the difference value integration circuit 153, and a two-dimensional coordinate which plots coordinate values obtained according to the processing of the flowcharts shown in FIGS.
  • processing circuit 15 determines the peak value of the discrimination target signal, the program counter 161 storing the address on the program storage device 19 currently being executed, the data acquisition circuit 162 acquiring data from the AD converter 14, and And a peak value determination circuit 163 as a functional block.
  • a waveform discrimination determination circuit 158, a waveform point accumulation circuit 159, an accumulated number display instruction circuit 160, a program counter 161, a data acquisition circuit 162, and a peak value determination circuit 163 are connected to one another via a data bus 164.
  • Discrimination determination circuit 158, waveform point accumulation circuit 159, cumulative number display instruction circuit 160, program counter 161, data acquisition circuit 162, peak value determination circuit 163 formally represent hardware resources focusing on logical functions. It does not necessarily mean functional blocks that exist independently as physical areas on a semiconductor chip, but also do not deny the configuration that actually exists.
  • the waveform discrimination apparatus includes an input device that receives an input of data, an instruction, and the like from an operator, an output device that outputs a classification result, and the like. May be further provided.
  • the input device is configured of a keyboard, a mouse, a light pen, a flexible disk device, and the like.
  • the operator of the waveform discrimination from the input device can designate input / output data, and can set individual numerical values, tolerance values, and degrees of error required for the waveform discrimination.
  • analysis parameters such as the form of output data from the input device, and it is also possible to input instructions such as execution or cancellation of calculation.
  • the output device and the display device 16 may be configured by a printer device and a display device, respectively.
  • the waveform discrimination apparatus can be realized by simple and inexpensive hardware resources as shown in FIG. 1 to FIG. 3 and FIG.
  • the whole structure of the waveform discrimination apparatus can be miniaturized by integrating on the circuit board of the above-mentioned structure, so that it is possible to achieve the effect that the portability of the waveform discrimination apparatus is easy.
  • the counting error due to the transient characteristic (slew rate) of the charge sensitive preamplifier used in this system configuration is There was a problem.
  • the waveform discrimination apparatus according to the first embodiment has a remarkable effect that the problem of counting error caused by the charge sensitive preamplifier can be avoided.
  • step S101 of FIG. 9 the window boundary condition determination circuit 151 of the signal processing circuit 15 shown in FIG. 13 determines the discrimination window boundary condition, and the address of the instruction to be read out next from the program storage device 19 by the program counter 161. After counting, the process of the signal processing circuit 15 is advanced to step S102.
  • step S102 of FIG. 9 the linear equation determination circuit 152 of the signal processing circuit 15 shown in FIG. 13 determines the discrimination linear equation, and the program counter 161 advances the processing of the signal processing circuit 15 to step S103.
  • the pulse included in the first pulse group or the pulse included in the second pulse group is input to the waveform detector 12 at an arbitrary timing, Since at least one of the first and second electrical signals is outputted in real time from the waveform detector 12 as a discrimination target signal at an arbitrary timing, the operation progress determination circuit 157 in step S104 operates in the data storage device 18 In the process of reading out the sample value U j stored in, the operation progress determination circuit 157 may directly capture the output of the AD converter 14 without passing through the data storage device 18.
  • step S104 If it is determined in step S104 that the sample value U j is larger than the lower limit identification value LLD (U) of the feature amount in the rising period, the process proceeds to step S105 and is stored in the data storage device 18 shown in FIG.
  • a microprocessor (MPU) register or the like can be used as the data storage device 18.
  • the difference value integration circuit 153 further reads the sample value U j + 1 stored in the data storage device 18, and the process proceeds to step S106.
  • the process of the difference value integration circuit 153 reading out the sample value U j + 1 stored in the data storage device 18 in step S105 is data
  • the sample value U j + 1 is directly taken from the AD converter 14 into the difference value integration circuit 153 in accordance with the timing at which the waveform detector 12 measures the first waveform or the second waveform without passing through the storage device 18 It can be done.
  • step S104 If it is determined in step S104 that the sample value U j is not larger than the lower limit identification value LLD (U) of the feature amount in the rising period, the process proceeds to step S108.
  • step S108 replaces the next sample value U j + 1 stored in the data storage device 18 to the new sample values U j, the new sample value U j arithmetic progression determination circuit 157 fetches the signal processing circuit 15 The process returns to step S104.
  • the calculation progress determination circuit 157 determines whether the difference value ⁇ U j + 1, j is larger than the lower limit identification value LLD (U) of the feature amount of the rising period.
  • step S106 If it is determined in step S106 that the difference value ⁇ U j + 1, j of the discrimination target signal is larger than the lower limit identification value LLD (U) of the feature amount in the rising period, the process proceeds to step S111 and the sample value U j + 1 And the difference value ⁇ U j + 1, j is stored in the data storage unit 18.
  • step S111 the difference value integration circuit 153 further reads out the sample value U j + 2 stored in the data storage device 18, and the process proceeds to step S112. Since the operation of the signal processing circuit 15 proceeds in real time simultaneously with the measurement, the process in which the difference value integration circuit 153 reads out the sample value U j + 2 stored in the data storage device 18 in step S111 is the data storage device 18.
  • the sample value U j + 1 can be taken directly from the AD converter 14 to the difference value integration circuit 153 in accordance with the timing at which the first waveform or the second waveform is measured without the intervention.
  • step S106 If it is determined in step S106 that the difference value ⁇ U j + 1, j is not larger than the lower limit identification value LLD (U) of the feature amount in the rising period, the process proceeds to step S107.
  • step S107 the next sample value Uj + 2 stored in the data storage device 18 is replaced with a new sample value Uj + 1 , and the process proceeds to step S108.
  • step S108 the sample value U j + 1 stored in the data storage unit 18 is replaced with the sample value U j , and this new sample value U j is taken into the operation progress determination circuit 157, and the process returns to step S104.
  • the calculation progress determination circuit 157 reads the difference value ⁇ U j + 1, j stored in the data storage device 18, and the difference value ⁇ U j + 2, j + 1 output from the difference value integration circuit 153 is It is determined whether the difference value ⁇ U j + 1, j is larger than the difference value ⁇ U j + 2, j + 1 is a positive value.
  • step S112 satisfied the difference value ⁇ U j + 2, j + 1 is the difference value ⁇ U j + 1, j is greater than or difference values ⁇ U j + 2, j + 1 is the one of a positive value conditions In this case, the difference value ⁇ U j +2, j + 1 is output to the difference value integration circuit 153 of the signal processing circuit 15, and the process proceeds to step S113.
  • step S112 also the difference value ⁇ U j + 2, j + 1 is the difference value ⁇ U j + 1, j is greater than or difference values ⁇ U j + 2, j + 1 is the one of a positive value conditions If not satisfied, the difference value ⁇ U j +2, j + 1 and the difference value ⁇ U j +1, j are output in parallel or sequentially to the difference value integration circuit 153, and the process proceeds to step S121.
  • step S113 the difference value integration circuit 153 of the signal processing circuit 15 reads the feature amount Us and the difference value ⁇ U j + 1, j from the data storage device 18, and Us + ⁇ U j + 1, j + ⁇ U j + 2, j + 1 Is calculated, and the calculation result is used as a new Us, and the process proceeds to step S114.
  • step S114 the address of the instruction the program counter 161 reads out next from the program storage device 19 is returned from j + 2 to j + 1, and the address of the next sample value U j + 1 stored in the data storage device 18 is a new sample value.
  • the operation progress determination circuit 157 replaces the sample value U j + 1 with the address of U j , reads the new sample value U j + 1 from the data storage device 18, and returns to step S106.
  • step S121 the difference value integration circuit 153 reads the feature amount Us from the data storage device 18, calculates the value of Us + ⁇ U j + 1, j + ⁇ U j + 2, j + 1 as the value Uf of the first coordinate axis, and the step Go to S122.
  • step S122 the peak value determining circuit 163 of the signal processing circuit 15 reads the specimen stored in the data storage device 18 value U j + 1 and sample values Uj + 2, sample values U j + 1 and the sample value U j Compare the size of +2 .
  • the processing of reading the sample value D j stored in the data storage unit 18 in step S202 is The sample value D j can be taken directly from the AD converter 14 to the attenuation amount calculation circuit 154 in accordance with the timing at which the first waveform or the second waveform is measured without passing through the data storage device 18.
  • the calculation progress determination circuit 157 determines whether the amount of attenuation D dj is larger than the lower limit identification value LLD (D) of the feature amount of the fall period.
  • step S202 If it is determined in step S202 that the amount of attenuation D dj is larger than the lower limit identification value LLD (D) of the feature amount in the falling period, the process proceeds to step S203, and the amount of attenuation D dj is stored in the data storage device 18.
  • LLD lower limit identification value
  • the processing for the attenuation amount calculation circuit 154 to read out the sample value D j + 1 stored in the data storage device 18 in step S203 does not involve the data storage device 18 and is performed at the timing when the first waveform or the second waveform is measured.
  • the sample value D j + 1 can be taken directly from the AD converter 14 to the attenuation calculation circuit 154.
  • step S202 If it is determined in step S202 that the attenuation amount D dj is not larger than the lower limit identification value LLD (D) of the feature amount in the falling period, the process proceeds to step S206.
  • step S206 replaces the next sample value D j + 1 stored in the data storage device 18 to the new sample values D j, the new sample value D j uptake attenuation amount calculation circuit 154, the signal processing circuit 15 The process returns to step S202.
  • the calculation progress determination circuit 157 determines whether the attenuation difference value ⁇ D j + 1, j is larger than the lower limit identification value LLD (D) of the feature amount in the fall period, or the attenuation amount D dj + 1. Is determined to be larger than the amount of attenuation D dj .
  • step S204 it is determined that the difference value ⁇ D j + 1, j of the attenuation amount is larger than the lower limit identification value LLD (D) of the feature amount of the fall period, or the attenuation amount D dj + 1 is larger than the attenuation amount D dj If it is determined, the process proceeds to step S 211, and the attenuation amount D dj + 1 and the difference value ⁇ D j + 1, j of the attenuation amount are stored in the data storage device 18.
  • LLD lower limit identification value
  • the process of the attenuation amount calculation circuit 154 reading out the sample value D j +2 stored in the data storage device 18 in step S211 is the attenuation amount calculation circuit from the AD converter 14 at the timing when the first waveform or the second waveform is measured.
  • the sample values D j + 2 can also be taken directly at 154 without going through the data store 18.
  • step S204 it is determined that the difference value ⁇ D j + 1, j of the attenuation amount is not larger than the lower limit identification value LLD (D) of the feature amount of the fall period, or the attenuation amount D dj + 1 is larger than the attenuation amount D dj If it is determined that there is not, the process proceeds to step S205. In step S205, the next sample value D j + 2 stored in the data storage device 18 is replaced with a new sample value D j + 1 , and the process proceeds to step S206.
  • LLD lower limit identification value
  • step S206 the sample value D j + 1 stored in the data storage device 18 is replaced with the sample value D j , the new sample value D j is taken into the attenuation amount calculation circuit 154, and the process returns to step S202.
  • step S212 whether the difference value ⁇ D j + 2, j + 1 of attenuation amounts is larger than the difference value ⁇ D j + 1, j of attenuation amounts or whether the attenuation amount D dj + 2 is larger than the attenuation amount D dj + 1 If one of the conditions is satisfied, the difference ⁇ D j + 2, j + 1 of the attenuation amount is output to the difference value integration circuit 153, and the process proceeds to step S213.
  • step S212 whether the difference value ⁇ D j + 2, j + 1 of attenuation amount is larger than the difference value ⁇ D j + 1, j of attenuation amount or the attenuation amount D dj + 2 is larger than the attenuation amount D dj + 1 If neither of the above conditions is satisfied, the attenuation difference values ⁇ D j + 2, j + 1 and ⁇ D j + 1, j are output in parallel or sequentially to the difference value integration circuit 153, and the process proceeds to step S221. .
  • step S213 the difference value integration circuit 153 reads out the difference value ⁇ D j + 1, j of the feature amount Ds and the attenuation amount from the data storage device 18 , and the value of Ds + ⁇ D j + 1, j + ⁇ D j + 2, j + 1 And calculate the calculation result as a new Ds, and the process proceeds to step S214.
  • step S214 the address of the instruction the program counter 161 reads out next from the program storage device 19 is returned from j + 2 to j + 1, and the address of the next attenuation amount D dj + 1 stored in the data storage device 18 is used as a new attenuation amount.
  • the new sample value D j + 1 is read out from the data storage unit 18 by the attenuation amount calculation circuit 154, and the process returns to step S 204.
  • step S221 the difference value integration circuit 153 reads the feature amount Ds from the data storage device 18, calculates the value of Ds + ⁇ D j + 1, j + ⁇ D j + 2, j + 1 as the value Df of the second coordinate axis, and The value Df of the two coordinate axes is stored in the data storage device 18, and the process proceeds to step S222.
  • FIG. 12 shows a point indicating coordinates (Uf, Df) formed by a combination of the value Df of the second coordinate axis and the value Uf of the first coordinate axis by the two-dimensional coordinate plotting circuit 156 of the signal processing circuit 15. Plot in the discriminant plane defined by the first coordinate axis and the second coordinate axis.
  • step S103 When the coordinates (Uf, Df) are plotted on the discrimination plane, the process returns to step S103, and the feature amount Us of the rising period is reset.
  • the program counter 161 causes the signal processing circuit 15 to synchronize with the clock signal to execute the processing shown in the flowcharts of FIGS. 9 and 10 from time to time.
  • the pulse included in the next first pulse group or the pulse included in the second pulse group is input to the waveform detector 12, and at least one of the first and second electrical signals from the waveform detector 12 is a discrimination target
  • a new point indicating the coordinates (Uf, Df) consisting of a set of the value Uf of the first coordinate axis and the value Df of the second coordinate axis by the two-dimensional coordinate plotting circuit 156 is shown in FIG. Plot in the discriminant plane defined by the first and second coordinate axes as shown.
  • the waveform discrimination method before the signal intensity falls to the baseline during the falling period of the pulse waveform, there is an input to the radiation-light conversion element 11 of gamma rays or neutrons, as shown in FIG.
  • a correct peak value can be obtained even when pile-up occurs in which the rising waveform of the next pulse signal is superimposed on the falling waveform of the electric signal output from the light detector 12a. That is, if pile-up occurs and it is determined in step S212 that the attenuation amount D dj +2 is smaller than the attenuation amount D dj +1 , the process returns to step S103 via steps S221 and S222.
  • FIG. 14 shows the case where pileup occurs in two places for the falling waveform, it can be determined that pileup has occurred in step S212 each time pileup occurs.
  • Us the rising period
  • one of the features of the waveform discrimination apparatus is that the analog amplifier 13 shown in FIG. 1 outputs the transient response waveform of the discrimination target signal output from the waveform detector 12 in time. It is to amplify to expand along an axis. By extending the fall time of the discrimination target signal along the time axis, it is possible to lengthen the sampling interval for the AD converter 14 to acquire digital data, so waveform discrimination according to the first embodiment is realized. According to the apparatus, it is possible to adopt a very inexpensive and simple AD converter 14. However, if the fall time is extended too much along the time axis, the probability of pile-up as shown in FIG.
  • the value of the input resistor R1 connected between the input terminal I of the analog amplifier 13 shown in FIG. 2 and the ground is appropriately selected according to the characteristics of the physical quantity of the pulse to be measured within the range of about 5 k ⁇ to 1 M ⁇ . Adjust to the value.
  • the value of the input resistance R1 is adjusted by further providing an input resistance adjusting knob similar to the adjusting knobs 34a, 34b, 34c and 34d provided on the bottom of the housing 21 shown in FIG.
  • the versatility may be enhanced by variably adjusting it while observing the characteristics of the physical quantity of the pulse to be measured.
  • the configuration of the waveform discrimination apparatus for realizing the waveform discrimination method according to the first embodiment is based on simple and inexpensive hardware resources as shown in FIGS. 1 to 3 and 13. As a result, the cost required for measurement can also be reduced. Further, since the waveform discrimination apparatus used for measurement is integrated on a small circuit board to be easily portable, the remarkable effect of improving the workability can be achieved.
  • the waveform point accumulation circuit 159 of the signal processing circuit 15 continuously repeats the feedback loop from step S222 of FIG. 10 to step S103 of FIG. 9 as long as the power for driving the signal processing circuit 15 is on.
  • new points are sequentially accumulated on the discrimination plane according to the repetition of the loop along the series of flows shown in FIGS. 9 and 10, a large number of coordinate points in the discrimination plane form the first pulse group. It is localized and plotted depending on whether it is the waveform of the included pulse or the waveform of the pulse included in the second pulse group. As shown in FIG.
  • the first pulse is generated by classifying and analyzing the localized area according to the flowchart shown in FIG. Whether it is the waveform of the pulse included in the group or the waveform of the pulse included in the second pulse group can be discriminated.
  • the waveform discrimination determination circuit 158 of the signal processing circuit 15 shown in FIG. 13 determines whether or not the position of the coordinate point is located inside the discrimination window. 12, the lower limit identification value LLD (D) of the feature amount Df in the falling period and the upper limit identification value ULD of the feature amount Df in the falling period along the X axis which is the second coordinate axis in FIG. D) is determined, and the lower limit identification value LLD (U) of the feature amount Uf in the rising period and the upper limit identification value ULD (U) of the feature amount Uf in the rising period are determined along the Y axis that is the first coordinate axis.
  • the lower limit identification value LLD (D), the upper limit identification value ULD (D), the lower limit identification value LLD (U) and the upper limit identification value ULD (U) are respectively determined in advance according to the procedure shown in FIG. It may be stored in the storage unit 18 and read from the data storage unit 18 when the position of the window portion is determined. That is, in FIG.
  • step S301 it is determined whether or not the distribution of coordinate points (Uf, Df) defined by a set of the value Uf of the first coordinate axis and the value Df of the second coordinate axis is located inside the window for discrimination.
  • the waveform discrimination determination circuit 158 makes this determination. If it is determined in step S301 that the distribution of coordinate points (Uf, Df) is not located inside the window for discrimination, the waveform discrimination determination circuit 158 outputs the waveform detector 12 in step S304. It is determined that the discrimination target signal is a signal whose source is the first waveform. On the other hand, if it is determined in step S301 that the distribution of coordinate points (Uf, Df) is located inside the window for discrimination, the process proceeds to step S302.
  • step S302 the distribution of coordinate points (Uf, Df) defined by the combination of the value Uf of the first coordinate axis and the value Df of the second coordinate axis is closer to the second coordinate axis than the straight line representing the discrimination linear equation.
  • the waveform discrimination determination circuit 158 determines whether or not it exists in the area.
  • the discrimination linear equation is expressed by a linear function of inclination a and an intercept b of the first coordinate axis as shown in FIG.
  • the values of the slope a and the intercept b of these discrimination linear equations are determined in advance according to the procedure shown in FIG. 16 and stored in the data storage unit 18 to determine the position of the window portion. It may be read from the storage device 18. That is, in FIG. 12, a discrimination linear equation is defined on the discrimination plane using the values of the inclination a and the intercept b stored in the data storage device 18. If it is determined in step S301 that the distribution of the coordinate points (Uf, Df) is not located closer to the second coordinate axis than the straight line representing the discrimination linear equation, the waveform discrimination determination circuit 158 detects the waveform in step S304. It is determined that the discrimination target signal output from the unit 12 is a signal whose source is the first waveform.
  • step S301 if it is determined in step S301 that the distribution of coordinate points (Uf, Df) is positioned closer to the second coordinate axis than the straight line representing the discrimination linear equation, the process proceeds to step S303, and the waveform discrimination determination circuit 158 It is determined that the discrimination target signal output from the waveform detector 12 is a signal having the second waveform as a generation source.
  • the waveform point accumulation circuit 159 can count the cumulative number of coordinates corresponding to the first waveform and the cumulative number of coordinates corresponding to the first waveform.
  • the accumulated number display command circuit 160 of the signal processing circuit 15 displays the accumulated number of coordinates corresponding to the first waveform and the second waveform respectively accumulated and counted by the waveform point accumulation circuit 159 as shown in FIGS. 1 and 3
  • the display 16 can be displayed on the display 16 by transmitting a display command and data necessary for the display.
  • step S401 in FIG. 16 the pulse included in the second pulse group for calibration whose waveform is known is input to the waveform detector 12.
  • the second electric signal output from the waveform detector 12 is sequentially output from the waveform detector 12 as a discrimination target signal, and the analog amplifier 13 expands the transient response waveform of the discrimination target signal along the time axis to perform AD conversion
  • the unit 14 samples the amplified discrimination target signal and converts it into digital data.
  • step S401 a plurality of digital data whose source is the second waveform for calibration is sequentially input in real time to the window boundary condition determination circuit 151 of the signal processing circuit 15 shown in FIG. When the second calibration waveform is measured, the process proceeds to step S402.
  • step S402 the window boundary condition determination circuit 151 sets the peak value of the rising period of the plurality of second electrical signals output from the waveform detector 12 to correspond to the plurality of calibration second waveforms, respectively.
  • a search is made by statistical processing using digital data sequentially converted by the AD converter 14, and the process proceeds to step S403.
  • step S403 the window boundary condition determination circuit 151 determines the lower limit identification value LLD (D) of the feature amount Df of the fall period and the feature amount of the fall period using the peak value of the rise period searched in step S402.
  • LLD lower limit identification value
  • the upper limit identification value ULD (D) of Df, the lower limit identification value LLD (U) of the feature amount Uf of the rising period, and the upper limit identification value ULD (U) of the feature amount Uf of the rising period are determined.
  • the values of LLD (D), ULD (D), LLD (U) and ULD (U) determined in step S403 are stored in the data storage device 18 in step S404.
  • step S411 the processing of the signal processing circuit 15 proceeds to step S411 by the program counter 161.
  • step S411 the pulse included in the first pulse group whose waveform is already known for calibration is input to the waveform detector 12, and a plurality of calibration first waveforms are measured.
  • step S412 the first electrical signal output from the waveform detector 12 is sequentially output as a discrimination target signal from the waveform detector 12, and the analog amplifier 13 expands the transient response waveform of the discrimination target signal along the time axis.
  • the AD converter 14 samples the amplified discrimination target signal and converts it into digital data. Then, according to the flowcharts shown in FIGS. 9 and 10, in step S412, the rising feature amount Uf and the falling feature amount Df are calculated.
  • step S415 the linear equation determination circuit 152 determines the intercept b of the discrimination linear equation.
  • the series of waveform discrimination operations shown in FIG. 9, FIG. 10, FIG. 15 and FIG. 16 are the waveform discrimination shown in FIG. 1 by a program for executing an algorithm equivalent to FIG. 9, FIG. 10, FIG. It can control and execute the device.
  • the waveform discrimination program may be stored in the program storage device 19 shown in FIG. Further, the waveform discrimination program is stored in a computer readable recording medium, and the program storage device 19 reads the recording medium to execute a series of waveform discrimination operations according to the first embodiment. Can.
  • the "computer readable recording medium” refers to, for example, a medium capable of recording various programs such as an external memory device of a microprocessor, a semiconductor memory, a magnetic disk, an optical disk, an optical magnetic disk, and a magnetic tape. It does not matter if it is. Specifically, a flexible disk, a CD-ROM, an MO disk, a cassette tape, an open reel tape and the like are included in the “computer readable recording medium”.
  • the waveform discrimination program is: (a) an instruction to input the waveform of the measured pulse to the waveform detector 12 and convert the physical quantity of the measured pulse into an electrical signal; (b) a command to cause the analog amplifier 13 to expand and amplify the transient response waveform of the electric signal along the time axis; (c) a command to cause the AD converter 14 to sample the amplified electric signal and convert it into digital data in the rising period and the falling period of the electric signal; (d) The difference value calculation circuit 153, the attenuation amount calculation circuit 154, and the difference value integration circuit 155 of the signal processing circuit 15 cooperate with each other to use the digital data to set the feature amount Uf of the rising period as a point on the first coordinate axis An instruction to calculate and calculate the feature amount Df of the falling period as a point on the second coordinate axis; (e) In the two-dimensional coordinate plotting circuit 156 of the signal processing circuit 15, a set of points on the first coordinate axis and points on the
  • the control circuit 17 and the signal processing circuit 15 of the waveform discrimination apparatus can be configured to internally or externally connect, for example, a flexible disk drive (flexible disk drive) and an optical disk drive (optical disk drive).
  • the flexible disk drive is inserted into the flexible disk drive and the CD-ROM is inserted into the optical disk drive from the insertion slot, and the waveform discrimination program stored in these recording media is performed by performing a predetermined read operation.
  • It can be installed in the program storage device 19 that constitutes the waveform discrimination device.
  • this waveform discrimination program can be stored in the program storage device 19 via an information processing network such as the Internet.
  • the present invention has been described according to the first embodiment, but it should not be understood that the statements and drawings that form a part of this disclosure limit the present invention. Various alternative embodiments, examples and operation techniques will be apparent to those skilled in the art from this disclosure.
  • the first waveform is the emission waveform from the radiation-light conversion element 11 specific to gamma rays
  • the second waveform is the radiation waveform from the radiation-light conversion element 11 specific to neutron rays.
  • the present invention is not limited to the description of the first embodiment.
  • the waveform detector 12 receives a sound wave having a first waveform and outputs a first electric signal, and receives a sound wave having a second waveform and outputs a second electric signal. It does not matter. Thus, it is a matter of course that the present invention includes various embodiments and the like which are not described herein. Accordingly, the technical scope of the present invention is defined only by the invention-specifying matters according to the scope of claims appropriate from the above description.
  • the present invention relates to a waveform discrimination apparatus, a waveform discrimination method, and a waveform discrimination program for discriminating two pulse waveforms having different waveforms from each other, for example, gamma rays and neutrons generated from radioactive materials which do not exist in nature used for nuclear power generation etc. Can be used to separate correctly.
  • a waveform discrimination apparatus for example, gamma rays and neutrons generated from radioactive materials which do not exist in nature used for nuclear power generation etc.
  • a waveform discrimination program for discriminating two pulse waveforms having different waveforms from each other, for example, gamma rays and neutrons generated from radioactive materials which do not exist in nature used for nuclear power generation etc. Can be used to separate correctly.
  • there is an industrial application value of clearly separating echoes from foreign substances which can not be found by measurement of propagation time in ultrasonic flaw detection, by providing means for discriminating two echo pulse waveforms different in waveform from each other.

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)
  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Measurement Of Radiation (AREA)

Abstract

　小型の回路基板上に集積化でき、簡単、安価でポータブル化が容易な波形弁別装置を提供する。被測定パルスの波形を入力し、電気信号に変換する波形検出器（１２）と、電気信号の過渡応答波形を時間軸に沿って拡大するアナログ増幅器（１３）と、電気信号の立ち上がり及び立ち下がり期間において、電気信号をデジタルデータに変換するＡＤ変換器（１４）と、デジタルデータを用いて、立ち上がり期間の特徴量を第１座標軸上の点として計算し、立ち下がり期間の特徴量を第２座標軸上の点として計算し、第１座標軸上の点及び第２座標軸上の点の組を座標点として、弁別平面に、座標点をプロットする信号処理回路（１５）とを備え、座標点のプロット位置から、被測定パルスが第１波形、或いは第１波形とは異なる第２波形であるのかを弁別する波形弁別装置。

Description

波形弁別装置、波形弁別方法及び波形弁別プログラム
　本発明は互いに波形の異なる２つのパルス波形を弁別する波形弁別装置に係り、特に、例えばガンマ線と中性子線がシンチレータに入射した場合のような、立ち上がり特性と立ち下がり特性がそれぞれ異なる２つのパルス波形の物理量を原因とする電気信号を弁別する波形弁別装置、波形弁別方法及び波形弁別プログラムに関する。
　ガンマ線と中性子線がシンチレータに入射し、シンチレータ内でそれぞれ異なる波形の光を発生する場合のように、互いに相似形をなす第１波形を有する任意の強度のパルスが任意のタイミングで発生して一群の集合を構成し、第１波形とは異なるが、互いに相似形をなす第２波形を有する任意の強度のパルスが任意のタイミングで発生して一群の集合を構成する場合がある。
　即ち、シンチレータの出力として、ガンマ線が入力されたことに起因する第１波形の光信号が任意の強度で任意のタイミングで発生して一群の光信号の集合を構成し、この一群の光信号を入力した光検出器が、第１波形に対応した第１の電気信号を逐次出力し、中性子線が入力されたことに起因する第２波形の光信号が任意の強度で任意のタイミングで発生して一群の光信号の集合を構成し、この一群の光信号を入力した光検出器が、第２波形に対応した第２の電気信号を逐次出力するような状況（環境）において、第１波形と、第１波形とは異なる第２波形を弁別したい場合がある。
　Ｃｅを添加したＬｉＣａＡｌＦ6結晶からなるシンチレータにガンマ線と中性子線が同時に入射した場合、シンチレータからは、ガンマ線に固有な波形の発光と中性子線に固有な波形の発光が発生することが知られている。このため、シンチレータからの２種類の発光を光電子増倍管で電気信号に変換し、光電子増倍管の出力を電荷感応型前置増幅器と波形整形増幅器（シェイピングアンプ）で増幅して２入力マルチチャネルアナライザ（ＭＣＡ）で分析するシステムが提案されている（非特許文献１参照。）。
　非特許文献１に記載の発明においては、波形整形増幅器の出力を波高値観測用出力と立ち上がり時間観測用出力に分岐し、波高値観測用出力を直接２入力ＭＣＡの一方の入力端子に入力し、立ち上がり時間観測用出力をパルス波形分析装置に入力している。そしてパルス波形分析装置から立ち上がり時間の１０％と９０％のタイミングの２つの信号を出力し、この２つの信号を時間／振幅変換器に入力し、時間／振幅変換器によって、２つの信号の時間差をパルスの振幅に変換し、時間／振幅変換器の出力を２入力ＭＣＡの他方の入力端子に入力するという、極めて複雑かつ大型で高価なシステム構成を用いている。
　非特許文献１に記載されている発明では、非特許文献１の図３に示されるようにガンマ線と中性子線の立ち上がり時間（Rise Time）とその波高値（Pulse Height）を座標面にプロットすることによって、ガンマ線と中性子線を分離して表示している。しかしながら、非特許文献１に記載されている発明では、非特許文献１の図３に示されるように、放射線の発生源としてのカルフォルニウム２５２（252Ｃｆ）から同時に放射される中性子とガンマ線により両者の波形が波高値軸でも時間軸でも重なっている。
　この理由は、非特許文献１の図１に示される電荷感応型前置増幅器の過渡特性（スルーレート）に起因する遅延によって高速な入力信号波高値が大きい信号は立上時間も大きくなってしまうためである。このため、非特許文献１に記載されている発明では、入射するガンマ線のエネルギーが中性子のエネルギー以上の場合、弁別が不可能となり、計数誤差を生じる。
　又、非特許文献１に記載されている発明が用いている測定システムでは、非特許文献１の図５に記載された２次元分布図において、本来抽出したい中性子を抽出する矩形の計数関心領域（ＲＯＩ)Ａの範囲内に、非測定対象であるガンマ線を抽出して抑制するための矩形の計数関心領域Ｂが設定されている。非特許文献１に記載されている発明では、計数関心領域Ａの領域内となる座標面において、計数関心領域Ｂの領域設定を、２次元分布図上でのデータのプロットを見ながら行う手動調整で行っており、誤差を少なくするためには人間の注意力と工夫が必要になる。しかし、非特許文献１の図５に示されるように、ガンマ線のプロット軌跡は非線形であり、中性子によるプロットとの分離を正確に実施することは困難である。
　即ち、非特許文献１に記載の発明のようなデスクトップサイズのＭＣＡを用いた大型で高価かつ複雑なシステム構成を用いたとしても、従来の技術によっては、ランダムにそれぞれ独立なイベントとして入射するガンマ線と中性子線とを、それぞれの線量がリアルタイムで計数できるように、それぞれの入力波形を弁別することはできない。又、ガンマ線と中性子線の線量をリアルタイムで正確に測定することはできない。
　又、非特許文献１では、ガンマ線と中性子線による入力信号が、短い時間間隔内で入射する場合への対応は考慮されておらず、パルス波形分析装置の時定数よりも短い時間間隔内でガンマ線と中性子線による入力信号があった場合にパイルアップを起こすことから、正確なエネルギーが計測できなくなり、精度が下がる。
山崎淳　他１１名著、「Ｃｅ：ＬｉＣａＡｌＦ６シンチレータを用いたパルス波形の弁別を基礎とした中性子／ガンマ線弁別(Neutron-gamma discrimination based on pulse shape discrimination in a Ce:LiCaAlF6 scintillator)」、 ニュークリア・インスツルメンツ・イン・フィジックス・リサーチ、Ａ６５２、２０１１年、ｐ．４３５－４３８
　上記事情を鑑み、本発明は、小型の回路基板上に集積化できるような簡単かつ安価な構成によって、ポータブル化が容易な波形弁別装置、波形弁別方法及び波形弁別プログラムを提供することを目的とする。
　上記目的を達成するために、本発明者らは、一例として、第１波形を有する任意の強度の物理量のパルスとして、ガンマ線がシンチレータに入射した際の発光による電流波形を、第２波形を有する任意の強度の物理量のパルスとして中性子がシンチレータに入射した際の発光による電流波形に着目した。
　即ち、例示的な検討例ではあるが、ガンマ線と中性子とが、同時又は一方がシンチレータに入射した際の発光による電流波形を比較すると、ガンマ線については、立ち上がり部の信号強度と立ち下がり部の減衰強度が線形比例すること、中性子については、非線形関係となる物理現象を応用して、第１波形と第２波形を有する任意の強度の物理量を正確に分離して計数する波形弁別装置、波形弁別方法及び波形弁別プログラムを見いだした。
　本発明の第１の態様は、(a)被測定パルスの波形を入力し、被測定パルスの物理量を電気信号に変換する波形検出器と、(b)電気信号の過渡応答波形を時間軸に沿って拡大して増幅するアナログ増幅器と、(c)電気信号の立ち上がり期間及び立ち下がり期間において、増幅された電気信号を標本化してデジタルデータに変換するＡＤ変換器と、(d)デジタルデータを用いて、立ち上がり期間の特徴量を第１座標軸上の点として計算し、立ち下がり期間の特徴量を第２座標軸上の点として計算し、更に、第１座標軸上の点及び第２座標軸上の点の組を座標点として、第１座標軸と第２座標軸が定義する弁別平面に、座標点をプロットする信号処理回路とを備える波形弁別装置であることを要旨とする。第１の態様に係る波形弁別装置は、座標点のプロット位置から、被測定パルスが第１波形、或いは第１波形とは異なる第２波形であるのかを弁別する。
　本発明の第２の態様は、(a)被測定パルスの波形を入力し、被測定パルスの物理量を電気信号に変換するステップと、(b)電気信号の過渡応答波形を時間軸に沿って拡大して増幅するステップと、(c)電気信号の立ち上がり期間及び立ち下がり期間において、増幅された電気信号を標本化してデジタルデータに変換するステップと、(d)デジタルデータを用いて、立ち上がり期間の特徴量を第１座標軸上の点として計算し、立ち下がり期間の特徴量を第２座標軸上の点として計算するステップと、(e)第１座標軸上の点及び第２座標軸上の点の組を座標点として、第１座標軸と第２座標軸が定義する弁別平面に、座標点をプロットするステップと、(f)座標点のプロット位置から、被測定パルスが第１波形、或いは第１波形とは異なる第２波形であるのか弁別するステップとを含む波形弁別方法であることを要旨とする。
　本発明の第２の態様で述べた波形弁別方法を実現するためのコンピュータソフトウェアプログラムは、コンピュータ読取り可能な記録媒体に保存し、この記録媒体をコンピュータシステムによって読み込ませることにより、本発明の波形弁別方法を実行することができる。
　即ち、本発明の第３の態様は、(a)波形検出器に被測定パルスの波形を入力させ、被測定パルスの物理量を電気信号に変換させる命令と、(b)アナログ増幅器に電気信号の過渡応答波形を時間軸に沿って拡大して増幅させる命令と、(c)ＡＤ変換器に電気信号の立ち上がり期間及び立ち下がり期間において、増幅された電気信号を標本化してデジタルデータに変換させる命令と、(d)信号処理回路の差分値計算回路、減衰量計算回路及び差分値積算回路を互いに連携させ、デジタルデータを用いて、立ち上がり期間の特徴量を第１座標軸上の点として計算させ、立ち下がり期間の特徴量を第２座標軸上の点として計算させる命令と、(e)信号処理回路の２次元座標プロット回路に、第１座標軸上の点及び第２座標軸上の点の組を座標点として、第１座標軸と第２座標軸が定義する弁別平面に、座標点をプロットさせる命令と、(f)信号処理回路の波形弁別判定回路に、座標点のプロット位置から、被測定パルスが第１波形、或いは第１波形とは異なる第２波形であるのか弁別させる命令とを含む一連の命令を制御回路に実行させる波形弁別プログラムであることを要旨とする。
　ここで、「記録媒体」とは、例えばコンピュータの外部メモリ装置、半導体メモリ、磁気ディスク、光ディスク、光磁気ディスク、磁気テープなどのプログラムを記録することができるような媒体などを意味する。具体的には、フレキシブルディスク、ＣＤ－ＲＯＭ，ＭＯディスク、カセットテープ、オープンリールテープなどが「記録媒体」に含まれる。第１の態様に係る波形弁別装置は、装置サイズを小型化することが可能であり、小型化の設計にさいしては、マイクロコントローラユニット（MCU）等の機器組み込み用プロセッサで、第３の態様に係る波形弁別プログラムを格納する記録媒体等の構成を実現することができる。ＭＣＵは当初、搭載メモリの少なさからアセンブリ言語でのみプログラムが組まれていたが、メモリ量やＣＰＵの処理能力が向上すると、開発効率の観点からＣ言語が使われるようになった。ＢＡＳＩＣ言語インタプリタなどの言語処理系が予めＲＯＭに書き込まれた半完成製品も存在し、第３の態様に係る波形弁別プログラムを格納する記録媒体等を実現することが可能である。
　本発明によれば、小型の回路基板上に集積化できるような簡単かつ安価な構成によって、ポータブル化が容易な波形弁別装置、波形弁別方法及び波形弁別プログラムを提供することができる。
本発明の第１の実施の形態に係る波形弁別装置の主要部の概略を説明する模式的なブロック図である。 第１の実施の形態に係る波形弁別装置に用いられるアナログ増幅器の一例を、より具体的に説明する図である。 第１の実施の形態に係る波形弁別装置の物理的な実装構造の主要部の概略を、筐体の側壁を透明にして模式的に説明する立体図である。 図４(a)は、光検出器としての光電子増倍管の出力を５０Ωで終端した場合において、ガンマ線がシンチレータに入射した際に、シンチレータから出射した第１波形の光が光電子増倍管に入射したことによる、光電子増倍管の出力としての第１の電気信号をオシロスコープで観測したパルス波形を示す図である。図４(ｂ)は、光電子増倍管の出力を５０ｋΩで終端した場合において、図４(a)に示した第１の電気信号をオシロスコープで観測したパルス波形を示す図である。 図４(a)の第１の電気信号のV部で示したパルス波形を時間軸を拡大して示す図であるが、比較のために第２の電気信号のパルス波形も第１の電気信号と時間軸を共通にして示している。 図６（ａ）は、図５と同様に、図４(a)のV部を時間軸を拡大して示す図で、図６（ｂ）は、比較のために第２の電気信号のパルス波形を、図６（ａ）と時間軸を共通にして示す図である。 図７（ａ）は、第１の実施の形態に係る波形弁別装置のアナログ増幅器の出力波形を示し、図７（ｂ）は、ＡＤ変換器が図７（ａ）の波形をデジタイズが可能な信号に変換している状態を概念的な理解を補助する模式図である。 ＡＤ変換器が出力したデジタルデータに対して、図９及び図１０の処理を実行し、立ち上がり期間の特徴量及び立ち下がり期間の特徴量を計算する場合の、演算の内容や前提とする定義を説明するための模式図である。 本発明の第１の実施の形態に係る波形弁別方法の基本となる２次元分布の作成方法の一例を説明する概念的なフローチャートである。 図９に続く第１の実施の形態に係る２次元分布の作成方法を説明する概念的なフローチャートである。 立ち上がり期間の特徴量を第１座標軸上の点、立ち下がり期間の特徴量第２座標軸上の点として、第１座標軸と第２座標軸が定義する弁別平面にプロットすると、ガンマ線（第１波形）と中性子線（第２波形）の分布領域が局在し、分類整理できることを説明する図である。 第１の実施の形態に係る波形弁別方法に用いる弁別用の窓部と、弁別用線形方程式を表す直線を説明する図である。 第１の実施の形態に係る波形弁別装置を構成する信号処理回路の内部構造を、論理的なハードウェア資源の組み合わせとして概念的に説明する図である。 第１の実施の形態に係る波形弁別方法によれば、光検出器の出力としての電気信号の立ち下がり波形に対し、次のパルス信号の立ち上がり波形が重畳するパイルアップが発生した場合でも正しい波高値を取得ができることを説明する図で、ＡＤ変換器における処理の状態を概念的に説明する模式図である。 第１の実施の形態に係る波形弁別方法における弁別用の窓部と、弁別用線形方程式を表す直線を用いた処理の流れを説明する概念的なフローチャートである。 第１の実施の形態に係る波形弁別方法に用いられる弁別用の窓部と、弁別用線形方程式を表す直線の決定方法を説明する概念的なフローチャートである。
　次に、図面を参照して、本発明の第１の実施の形態を説明する。以下の図面の記載において、同一又は類似の部分には同一又は類似の符号を付している。ただし、図面は模式的なものであり、具体的な厚みや寸法は以下の説明を参酌して判断すべきものである。又、図面相互間においても互いの寸法の関係や比率が異なる部分が含まれていることは勿論である。
　又、以下に示す第１の実施の形態は、本発明の技術的思想を具体化するための装置や方法を例示するものであって、本発明の技術的思想は、構成部品の材料、形状、構造、配置等を下記のものに特定するものでない。本発明の技術的思想は、特許請求の範囲に記載された請求項が規定する技術的範囲内において、種々の変更を加えることができる。
（波形弁別装置の構成）
　本発明の第１の実施の形態に係る波形弁別装置は、図１に示すように、被測定パルスの波形を入力し、被測定パルスの物理量を電気信号に変換する波形検出器１２と、波形検出器１２に接続され、電気信号の過渡応答波形を時間軸に沿って拡大して増幅するアナログ増幅器１３と、アナログ増幅器１３に接続され、電気信号の立ち上がり期間及び立ち下がり期間において、増幅された電気信号を標本化してデジタルデータに変換するＡＤ変換器１４と、ＡＤ変換器１４に接続され、デジタルデータを用いて、立ち上がり期間の特徴量を第１座標軸上の点として計算し、立ち下がり期間の特徴量を第２座標軸上の点として計算し、更に、第１座標軸上の点及び第２座標軸上の点の組を座標点として、第１座標軸と第２座標軸が定義する弁別平面に、座標点をプロットする信号処理回路１５と、信号処理回路１５に接続された表示装置１６及びデータ記憶装置１８と、波形検出器１２、アナログ増幅器１３、ＡＤ変換器１４、信号処理回路１５及び表示装置１６に接続された制御回路１７と、制御回路１７に接続されたプログラム記憶装置１９とを備える。
　なお、図１にはデータ記憶装置１８及びプログラム記憶装置１９が、便宜上、それぞれ単独のハードウェア資源であるかのように図示されているが、現実の物理的ハードウェア資源としては、データ記憶装置１８及びプログラム記憶装置１９が、それぞれ、機能及び規模を異にする複数の記憶装置の集合として構成されている構成を否定するものではない。例えば、データ記憶装置１８は、複数のレジスタ、複数のキャッシュメモリ、主記憶装置、補助記憶装置を含む一群の内から適宜選択された任意の組み合わせとすることも可能である。又、キャッシュメモリは１次キャッシュメモリと２次キャッシュメモリの組み合わせとしてもよく、更に３次キャッシュメモリを備えるヒエラルキーを有しても構わない。
　図１に示す波形検出器１２は、被測定パルスとしての第１波形を有する任意の強度のパルスが任意のタイミングで発生して一群の集合を構成した第１のパルス群に含まれる少なくとも１つのパルスを入力し、第１波形に対応した第１の電気信号を逐次出力し、他の被測定パルスとして、第１波形とは異なる第２波形を有する任意の強度のパルスが任意のタイミングで発生して一群の集合を構成した第２のパルス群に含まれる少なくとも１つのパルスを入力し、第２波形に対応した第２の電気信号を逐次出力する。
　アナログ増幅器１３は、第１及び第２の電気信号の少なくとも一方を、波形検出器１２から弁別対象信号として入力し、弁別対象信号の過渡応答波形を時間軸に沿って拡大するように弁別対象信号の波形を増幅する。アナログ増幅器１３は、被測定パルスとしての第１波形がナノ秒レベルの半値幅のパルスであっても、立ち下がり時間を２μ秒程度以上の長さにとなるように過渡応答特性を示す波形を時間軸に沿って拡大することが好ましい。立ち下がり時間がマイクロ秒レベルになると、ＡＤ変換器１４がデジタルデータを取得するためのサンプリング間隔を長くすることができるので、非常に安価かつ簡単なＡＤ変換器１４が採用可能となる。そして、ＡＤ変換器１４は、弁別対象信号の立ち上がり期間及び立ち下がり期間において、増幅された弁別対象信号を標本化し、それぞれ一定間隔の離散的データを生成し、この離散的データをデジタルデータに変換する。
　ＡＤ変換器１４が標本化して生成した離散的なデジタルデータを、信号処理回路１５のデータ取込回路１６２（図１３参照。）が、ＡＤ変換器１４から順次読み込み、データ記憶装置１８に一時的に格納する。更に、信号処理回路１５は、データ記憶装置１８に格納されたデジタルデータを読み出し、立ち下がり期間の特徴量Ｄfを第２座標軸上の点として計算し、立ち上がり期間の特徴量Ｕfを第１座標軸に直交する第１座標軸上の点として計算する。そして、第１座標軸上の点及び第２座標軸上の点の組（Ｕf，Ｄf）を座標点として、図１１に示すように、信号処理回路１５は座標点（Ｕf，Ｄf）を、ＡＤ変換器１４による過渡応答波形の標本化の処理の進行とともに、第１座標軸と第２座標軸が定義する弁別平面にリアルタイムで自動的にプロットする。座標点（Ｕf，Ｄf）を図１１に示す弁別平面にリアルタイムでプロットしているので、データ記憶装置１８は、ＡＤ変換器１４が出力するデジタルデータを一時的に記憶するレジスタとして機能している。
　図１１及びこの図１１に続く図１２では、第２座標軸をＸ軸方向としているが、立ち下がり期間の特徴量Ｄfの値は、Ｘ軸方向の右側に進むに従い大きくなる値としてプロットされている。一方、図１１及び図１２では、第１座標軸をＹ軸として表示しているが、立ち上がり期間の特徴量Ｕfの値は、Ｙ軸方向の上側に進むに従い小さくなり、Ｙ軸方向の下側に進むに従い大きくなる値としてプロットされているので、第１座標軸と第２座標軸が定義する弁別平面は第３象限のデカルト座標として定義されている。なお、図１１及び図１２では、第１座標軸をＹ軸、第２座標軸をＸ軸として表示しているが、例示に過ぎず、Ｘ軸とＹ軸のいずれを第１座標軸として選択し、残るいずれを第２座標軸として選択するかは適宜決めてよい。
　信号処理回路１５の論理的なハードウェア資源の構成は図１３のように表現される。信号処理回路１５は、マイクロチップとして実装されたマイクロプロセッサ（ＭＰＵ）等が使用可能である。又、信号処理回路１５として、算術演算機能を強化し信号処理に特化したデジタルシグナルプロセッサ（ＤＳＰ）や、メモリや周辺回路を搭載し組込み機器制御を目的としたマイクロコントローラ（マイコン）等を用いてもよい。或いは、現在の汎用コンピュータのメインＣＰＵを信号処理回路１５に用いてもよい。
　第１の実施の形態に係る波形弁別装置は、図１１に示す弁別平面上の座標点の分布位置から、弁別対象信号が第１波形を発生源とするのか、或いは第２波形を発生源とするのかをリアルタイムで自動的に弁別することができる。特に、図１１に示す弁別平面に対し、図１２に示すように、弁別用の窓部及び弁別用線形方程式を定義することにより、弁別対象信号が第１波形を発生源とするのか、或いは第２波形を発生源とするのかを、コンピュータソフトウェアの処理にしたがって、リアルタイムで自動的に弁別することができる。
　図２及び図３に示す構成は、第１の実施の形態に係る波形弁別装置の具体的応用例に係り、被測定パルスとしての第１波形がガンマ線に固有な放射線－光変換素子１１からの発光波形、他の被測定パルスとしての第２波形が中性子線に固有な放射線－光変換素子１１からの発光波形の場合である。即ち、以後の説明では、第１波形が、ガンマ線が放射線－光変換素子１１に入力されたことに起因する放射線－光変換素子１１からの発光波形、第２波形が、中性子線が放射線－光変換素子１１に入力されたことに起因する放射線－光変換素子１１からの発光波形として例示的に説明する。
　即ち、第１の実施の形態に係る波形弁別装置は、図２及び図３に示すように、中性子線及びガンマ線を光に変換する放射線－光変換素子１１と、放射線－光変換素子１１に接続された放射線－光変換素子１１の発光する光を電気信号に変換する光検出器１２ａとを備えている。中性子線及びガンマ線をそれぞれ過渡応答特性の波形が異なる物理量である光に変換する放射線－光変換素子１１としては、表１に示すようなＣｓＬｉＹＣｌ，ＬｉＣａＡｌＦ6，ＬｉＦ／ＺｎＳ，ＬｉＢａＦ3，Ｌｉ6Ｇｄ（ＢＯ3）3等を使用することが可能である。　
Figure JPOXMLDOC01-appb-T000001
　表１に示すように、より強い発光のためには発光中心となる元素、例えばＹ，Ｃｅ，Ｐｒ，Ｓｍ，Ｅｕ，ＴｂやＭｎなどをＣｓＬｉＹＣｌ，ＬｉＣａＡｌＦ6，ＬｉＦ／ＺｎＳ，ＬｉＢａＦ3，Ｌｉ6Ｇｄ（ＢＯ3）3等シンチレータ材料に添加するのが好ましい。例えば、Ｃｅを添加したＬｉＣａＡｌＦ6を放射線－光変換素子１１とする場合は、図５に示すように、放射線－光変換素子１１からの発光の過渡応答特性は、ガンマ線（第１波形）と中性子線（第２波形）とで、立ち上がり特性及び立ち下がり特性がそれぞれ異なる。
　図５において、ガンマ線（第１波形）での発光は、数ナノ秒程度の非常に発光時間の短い発光と、先端の鋭いピークに続くブロードな発光をしている。一方、中性子線（第２波形）の発光は、図５に示すように、数百ナノ秒以下程度の比較的発光時間の長い発光に特徴がある。ガンマ線（第１波形）での発光はチェレンコフ発光であると言われている。
　表１に示すように、ＣｓＬｉＹＣｌ，ＬｉＣａＡｌＦ6，ＬｉＦ／ＺｎＳ，ＬｉＢａＦ3，Ｌｉ6Ｇｄ（ＢＯ3）3等のシンチレータは、波長１９０～４５０ｎｍ程度の光を発光するので、放射線－光変換素子１１の発光する光を電気信号に変換する光検出器１２ａとしては、波長１９０～４５０ｎｍ程度の光を電気信号に変換できる光電子増倍管（ＰＭＴ）、半導体フォトダイオード、フォトダイオードアレイ、ガイガーモード並列読み出しＡＰＤピクセルアレイ等が使用可能である。光検出器１２ａには、第１波形のパルスを入力して第１波形に対応した第１の電気信号が出力され、第２波形のパルスを入力して第２波形に対応した第２の電気信号を出力する特性が要求され、光検出器１２ａの入力と出力に線形性が保たれるデバイス性能が理想である。
　図３に示すように、放射線－光変換素子１１は、光検出器１２ａの窓部に貼り付けられ、光検出器１２ａは、筐体２１の上部突出部として構成されている。筐体２１の内部には、光検出器１２ａの出力にケーブル３１ａ，３１ｂを介して接続された回路基板２３と、回路基板２３に搭載され、回路基板２３中の埋め込配線又は表面配線を介して、ケーブル３１ａ，３１ｂに電気的に接続されるアナログ増幅器１３及び高圧電源２２と、回路基板２３にケーブル３２ａ，３２ｂ、３２ｃを介して接続された回路基板２４と、回路基板２４に搭載され、回路基板２４中の埋め込配線又は表面配線を介して、ケーブル３２ａ，３２ｂ、３２ｃに電気的に接続されるＡＤ変換器１４及び信号処理回路１５が内蔵されている。
　図３においては、筐体２１は例示的に直方体の形状をなし、筐体２１の１側面の上部には表示装置１６が取り付けられている。筐体２１の形状は、直方体の形状に限定されるものではなく、円筒形等の他の形状でもよく、円筒形の筐体２１であれば、筐体２１の円周面の一部を平坦にして、表示装置１６を埋め込む構造や、一部が筐体２１の円周面から突出するようなトポロジーで取り付けても構わない。
　筐体２１の底面には信号処理回路１５の条件を設定する調整つまみ３４ａ，３４ｂ，３４ｃ，３４ｄが設けられている。筐体２１の底面には孔が設けられ、この孔から、信号処理回路１５に回路基板２４中の埋め込配線又は表面配線を介して接続される通信用ケーブル３３が筐体２１の外部に導かれている。
　図３には図示を省略しているが、図１のデータ記憶装置１８及びプログラム記憶装置１９を含むシステム構成は例示に過ぎず、データ記憶装置１８が図３に示したＡＤ変換器１４や信号処理回路１５等の内部構造として存在してもよく、信号処理回路１５の内部メモリとして一部の記憶装置の機能をレジスタ等の形態で分散させ、残余の機能を回路基板２４上に実装される外部メモリにおいて実行させるように、物理的な構造としては、データ記憶装置１８を分散配置してもよい。或いは、回路基板２４上に実装される外部メモリのみをデータ記憶装置１８として備えるような物理的な構成でもよく、通信用ケーブル３３を介して接続された筐体２１の外部に配置されたデータ記憶装置１８を含んでもよい。
　同様に、図１のプログラム記憶装置１９についても、信号処理回路１５や制御回路１７等の内部構造として存在してもよく、信号処理回路１５や制御回路１７の内部メモリとしての記憶装置と、外部メモリとしての記憶装置の両方を含むように構成してもよく、外部メモリのみの構成として存在する記憶装置でもよい。
　更に、図１の制御回路１７についても、その少なくとも一部の機能を分散させることにより、現実の物理的な構造として制御回路１７を筐体２１の内部構造として存在させてもよい。場合によっては、制御回路１７を信号処理回路１５の内部構造として構成することや、逆に、制御回路１７の内部構造に信号処理回路１５やＡＤ変換器１４を機能ブロック等の形で集積化するような種々の物理的な構造を実現してもよい。
　図３に例示したように、第１の実施の形態に係る波形弁別装置は、簡単な構成であるので、装置サイズを小型化することが可能であり、マイクロコントローラユニット（MCU）等の機器組み込み用プロセッサで実現することができる。ＭＣＵは、図１に示したアナログ増幅器１３、ＡＤ変換器１４、信号処理回路１５及び制御回路１７等を含むコンピュータシステムを１つの集積回路に組み込んだものである。ＭＣＵは、自己充足性と低価格性を重視したタイプのマイクロプロセッサと言え、半導体の１チップのみでコンピュータとして機能させるということが可能である。ＭＣＵで構成すれば、汎用ＣＰＵと比較した場合に周辺部品が少なくて済むため、第１の実施の形態に係る波形弁別装置をコンパクトに組み立てるのが容易である。
　図２に示すように、光検出器１２ａの出力側にはアナログ増幅器１３の入力端子Ｉが、光検出器１２ａの信号出力端子と基準電位点端子が、各々アナログ増幅器１３の入力端子Ｉと接地間に接続されるような構成で、接続されている。アナログ増幅器１３の入力端子Ｉはアナログ増幅器１３の入力段を構成する第１オペアンプＵ１の非反転入力端子に接続され、アナログ増幅器１３の入力端子Ｉと接地間には入力抵抗Ｒ１が接続されている。入力端子Ｉと光検出器１２ａの出力電流の基準電位とする回路の接地間に接続される入力抵抗Ｒ１の値は、５ｋΩ以上が好ましく、例えば、Ｒ１＝５０ｋΩ～１ＭΩ程度の大きな値とすること好ましい。通常、アナログ増幅器として用いるオペアンプの入力抵抗（Ｉｍｐ）は、１０Ｍオーム乃至１ＴΩであるので、入力抵抗Ｒ１の最大値は、アナログ増幅器１３の入力段を構成する第１オペアンプＵ１の入力抵抗を考慮して決めればよい。第１オペアンプＵ１の反転入力端子がバイアス補償抵抗Ｒ６を介して接地されている。第１オペアンプＵ１の反転入力端子は更に帰還抵抗Ｒ５を介して第１オペアンプＵ１の出力端子に接続されている。第１オペアンプＵ１の出力端子は更に、伝送抵抗Ｒ２を介して第２オペアンプＵ２の反転入力端子に接続され、第２オペアンプＵ２の反転入力端子は、帰還抵抗Ｒ３を介して第２オペアンプＵ２の出力端子に接続されている。
　第２オペアンプＵ２の出力端子は更に、伝送抵抗Ｒ４を介してアナログ増幅器１３の出力段を構成する第３オペアンプＵ３の非反転入力端子に接続され、第３ペアンプＵ２の反転入力端子は第３オペアンプＵ３の出力端子に直結され、第３オペアンプＵ３の出力端子がアナログ増幅器１３の出力端子Oを構成している。
　図２に示すようなアナログ増幅器１３の回路を構成することで、第１の実施の形態に係る波形弁別装置では、放射線－光変換素子１１で変換された光を光検出器１２ａで第１又は第２の電気信号に逐次変換し、アナログ増幅器１３にて、第１又は第２の電気信号を、対応する第１又は第２の電圧信号への変換と増幅を行うことができる。
　図２においては、入力抵抗Ｒ１の値を、例えば、５０ｋΩ程度の大きな値に設定することにより、光検出器１２ａの出力端子間容量Ｃｐとアナログ増幅器１３の入力抵抗Ｒ１による減衰時定数τ＝Ｒ１・Ｃｐの値を大きな値として、過渡応答波形を時間軸に沿って拡大し、本来弁別したい第２波形（中性子線）と第１波形（ガンマ線）の高周波成分を低周波数帯シフトさせることができる。
　図４(a)は、光検出器１２ａの出力を５０Ωで終端した場合において、第１波形（ガンマ線）の光が光検出器１２ａに入射したことにより光検出器１２ａから出力される第１の電気信号のパルス波形である。光検出器１２ａの出力を５０ｋΩで終端した場合は、５０Ωで終端した場合に比して、図４（ｂ）に示すように１０００倍程度に過渡応答波形が時間軸方向に拡大されることが分かる。
　図５は、図４(a)の第１の電気信号のV部で示したパルス波形を時間軸を拡大して示す図であるが、比較のために第２の電気信号のパルス波形も時間軸を共通にして示している。図５は、Ｃｅを添加したＬｉＣａＡｌＦ6結晶を放射線－光変換素子（シンチレータ）１１として用い、放射線－光変換素子１１からの発光を光検出器（光電子増倍管）１２ａで検出した場合の第１及び第２の電気信号のパルス波形を示す。
　図６（ａ）は、図５と同様に、図４(a)の第１の電気信号のV部で示したパルス波形を時間軸を拡大して示す図である。図６（ａ）に示すように第１波形（ガンマ線）は、半値幅４ｎｓ程度の急峻な立ち上がり／立ち下がり特性を示す尖塔部分と、この尖塔部分が立ち上がったのち再び緩やかな立ち上がり／立ち下がり特性を示す丘状部分の２つの部分を有することが分かる。図６（ｂ）は、図６（ａ）との比較のために、第２の電気信号のパルス波形を、図６（ａ）と時間軸を共通にして示しているが、第２の電気信号には、図６（ａ）に示すような急峻な立ち上がり／立ち下がり特性を示す尖塔部分が存在しないことが分かる。
　第１の実施の形態に係る波形弁別装置では、アナログ増幅器１３の入力抵抗Ｒ１の値を５０ｋΩ程度の大きな値に設定して減衰時定数τ＝Ｒ１・Ｃpの値を大きな値としているので、アナログ増幅器１３は、図４(a)の第１の電気信号の過渡応答波形を、図７（ａ）に示すように立ち下がり時間が２μ秒程度以上となるように時間軸に沿って拡大する。即ち、アナログ増幅器１３は、汎用のＡＤ変換器１４が図７（ｂ）に示すようなデジタイズが可能な信号に、第１及び第２の電気信号の少なくとも一方を変換する。
　ＡＤ変換器１４は、第１及び第２の電気信号の波高値だけでなく減衰時間の違いを取得することで、ＡＤ変換器１４に接続された信号処理回路１５が、図９及び図１０に示すフローチャートにしたがって、図１２に示すような弁別平面の上に座標点を逐次生成し、２次元分布が作成される。第１座標軸と第２座標軸が定義する２軸間の相関の違いから、図１２に示すように２次元分布の領域が決定され、第１波形（ガンマ線）と第２波形（中性子線）の立ち上がり特性及び立ち下がり特性の違いが判断できるので、第１波形（ガンマ線）と第２波形（中性子線）とを弁別することが可能となる。
　図１３に論理的なハードウェア資源の構成を示すように、第１の実施の形態に係る波形弁別装置の信号処理回路１５は、図９及び図１０に示す処理を実行するに先立ち、校正用の波形を波形検出器１２に入力させて、波形の弁別に必要な窓部境界条件を決定する窓部境界条件決定回路１５１、同様に、図９及び図１０に示す処理を実行するに先立ち、校正用の波形を波形検出器１２に入力させて、波形の弁別に必要な線形方程式を決定する線形方程式決定回路１５２と、ＡＤ変換器１４が生成した一定間隔の離散的データの互いに連続する２値間の差分を計算する差分値積算回路１５３と、弁別対象信号の立ち下がり期間において、弁別対象信号の立ち上がり期間において実現された波形のピーク値からの減衰量を計算する減衰量計算回路１５４と、差分値積算回路１５３が出力した差分値を積算する差分値積算回路１５５と、図９及び図１０に示すフローチャートの処理にしたがって得られた座標値を２次元空間にプロットする２次元座標プロット回路１５６と、図９及び図１０に示すフローチャートの処理に伴う演算の進行を判定し、分岐の方向を決定する演算進行判定回路１５７と、２次元空間にプロットされた座標点の分布から波形を弁別して判定する波形弁別判定回路１５８と、波形弁別判定回路１５８が判定した波形点を累積する波形点累積回路１５９と、波形点累積回路１５９が計数した累積数を表示するように命令する累積数表示命令回路１６０と、図１のプログラム記憶装置１９に格納された実行すべき命令を信号処理回路１５がどこまで実行したか、或いは信号処理回路１５が、現在実行しているプログラム記憶装置１９上のアドレスを記憶するプログラムカウンタ１６１と、ＡＤ変換器１４からデータを取込むデータ取込回路１６２と、弁別対象信号のピーク値を決定するピーク値決定回路１６３とを機能ブロックとして備える。
　図１３に示すように、窓部境界条件決定回路１５１、線形方程式決定回路１５２、差分値積算回路１５３、減衰量計算回路１５４、差分値積算回路１５５、２次元座標プロット回路１５６、演算進行判定回路１５７、波形弁別判定回路１５８、波形点累積回路１５９、累積数表示命令回路１６０、プログラムカウンタ１６１、データ取込回路１６２、ピーク値決定回路１６３は、データバス１６４を介して互いに接続されている。図１３に示す窓部境界条件決定回路１５１、線形方程式決定回路１５２、差分値積算回路１５３、減衰量計算回路１５４、差分値積算回路１５５、２次元座標プロット回路１５６、演算進行判定回路１５７、波形弁別判定回路１５８、波形点累積回路１５９、累積数表示命令回路１６０、プログラムカウンタ１６１、データ取込回路１６２、ピーク値決定回路１６３は、論理的な機能に着目したハードウェア資源を形式的に表現しているのであって、必ずしも、半導体チップ上に物理的な領域としてそれぞれ独立して存在する機能ブロックを意味するものではないが、現実に存在する構成を否定するものでもない。
　図１及び図１３では図示を省略しているが、第１の実施の形態に係る波形弁別装置は、操作者からのデータや命令などの入力を受け付ける入力装置、分別結果を出力する出力装置等を更に備えるようにしてもよい。入力装置はキーボード、マウス、ライトペン又はフレキシブルディスク装置などで構成される。入力装置より波形弁別の実行者は、入出力データを指定したり、波形弁別に必要な個別の数値や許容誤差の値及び誤差の程度を設定したりすることができる。更に、入力装置より出力データの形態等の解析パラメータを設定することも可能で、又、演算の実行や中止等の指示の入力も可能である。又出力装置及び表示装置１６は、それぞれプリンタ装置及びディスプレイ装置等により構成してもよい。
　第１の実施の形態に係る波形弁別装置によれば、図１～図３及び図１３に示したような簡単かつ安価なハードウェア資源で実現可能であるので、波形弁別装置の主要部を小型の回路基板上に集積化して、波形弁別装置の全体の構造を小型化できるので、波形弁別装置のポータブル化が容易であるという効果を奏することができる。
　特に、放射線計測応用の分野においては、非特許文献１に記載の発明のシステム構成では、このシステム構成に用いられていた電荷感応型前置増幅器の過渡特性（スルーレート）に起因する計数誤差の問題があった。第１の実施の形態に係る波形弁別装置によれば、電荷感応型前置増幅器に起因した計数誤差の問題を回避することができるという顕著な効果を奏することができる。
（２次元分布の作成）
　図９～図１３を用いて、本発明の第１の実施の形態に係る波形弁別方法の基本となる２次元分布の作成方法を説明する。なお、以下に述べる２次元分布の作成方法は、一例であり、特許請求の範囲に記載した趣旨の範囲内であれば、この変形例を含めて、これ以外の種々の２次元分布の作成方法が実現可能であることは勿論である。
　図９のステップＳ１０１において、図１３に示した信号処理回路１５の窓部境界条件決定回路１５１が弁別窓部境界条件を決定し、プログラムカウンタ１６１がプログラム記憶装置１９から次に読み出す命令のアドレスをカウントして、ステップＳ１０２に信号処理回路１５の処理を進める。図９のステップＳ１０２において、図１３に示した信号処理回路１５の線形方程式決定回路１５２が弁別用線形方程式を決定し、プログラムカウンタ１６１がステップＳ１０３に信号処理回路１５の処理を進める。
　ステップＳ１０３において、図１３に示した信号処理回路１５の演算進行判定回路１５７が立ち上がり期間の特徴量Ｕｓをリセットし、特徴量Ｕｓ＝０の値を図１３に示したデータ記憶装置１８に格納する。その後、プログラムカウンタ１６１によってステップＳ１０４に信号処理回路１５の処理が進む。ステップＳ１０４において、演算進行判定回路１５７が弁別対象信号の標本値Ｕjをデータ記憶装置１８から読み出し、標本値Ｕjが立ち上がり期間の特徴量の下限識別値ＬＬＤ（Ｕ）より大きいか否かを判定する。図８の例ではｊ＝ｍ－１として、標本値Ｕmが立ち上がり期間の特徴量の下限識別値ＬＬＤ（Ｕ）より大きいか否かを判定する。
　図９及び図１０のフローチャートに示した信号処理回路１５の動作は、第１のパルス群に含まれるパルス又は第２のパルス群に含まれるパルスが任意のタイミングで波形検出器１２に入力され、第１及び第２の電気信号の少なくとも一方が、波形検出器１２から弁別対象信号として任意のタイミングで出力されることによりリアルタイムで進行するので、ステップＳ１０４における演算進行判定回路１５７がデータ記憶装置１８に格納された標本値Ｕjを読み出す処理は、データ記憶装置１８を介さずに、ＡＤ変換器１４の出力を直接、演算進行判定回路１５７が取り込むようにしてもよい。
　ステップＳ１０４で、標本値Ｕjが立ち上がり期間の特徴量の下限識別値ＬＬＤ（Ｕ）より大きいと判断された場合は、ステップＳ１０５に進み、図１３に示したデータ記憶装置１８に格納する。データ記憶装置１８としては、マイクロプロセッサ（ＭＰＵ）のレジスタ等が使用可能である。ステップＳ１０５では更に、差分値積算回路１５３が、データ記憶装置１８に格納された標本値Ｕj+1を読み出してステップＳ１０６に進む。上述したとおり、信号処理回路１５の動作は、測定と同時にリアルタイムで進行するので、ステップＳ１０５におけるデータ記憶装置１８に格納された標本値Ｕj+1を差分値積算回路１５３が読み出す処理は、データ記憶装置１８を介さずに、波形検出器１２が第１波形又は第２波形を測定したタイミングに対応して、ＡＤ変換器１４から差分値積算回路１５３に標本値Ｕj+1が直接取り込まれるようにできる。
　ステップＳ１０４で、標本値Ｕjが立ち上がり期間の特徴量の下限識別値ＬＬＤ（Ｕ）より大きくないと判断された場合は、ステップＳ１０８に進む。ステップＳ１０８では、データ記憶装置１８に格納された次の標本値Ｕj+1を新たな標本値Ｕjに置き換え、この新たな標本値Ｕjを演算進行判定回路１５７が取り込み、信号処理回路１５の処理はステップＳ１０４に戻る。
　ステップＳ１０６では差分値積算回路１５３が、データ記憶装置１８に格納された標本値Ｕjを読み出し、差分値ΔＵj+1,j＝Ｕj+1－Ｕjを計算し、計算結果を演算進行判定回路１５７に出力する。図８のｊ＝ｍの例では、差分値ΔＵm+1,m＝Ｕm+1－Ｕmを計算し、計算結果を演算進行判定回路１５７に出力する。ステップＳ１０６において、演算進行判定回路１５７は、差分値ΔＵj+1,jが立ち上がり期間の特徴量の下限識別値ＬＬＤ（Ｕ）より大きいか否かを判定する。
　ステップＳ１０６で、弁別対象信号の差分値ΔＵj+1,jが立ち上がり期間の特徴量の下限識別値ＬＬＤ（Ｕ）より大きいと判断された場合は、ステップＳ１１１に進み、標本値Ｕj+1及び差分値ΔＵj+1,jをデータ記憶装置１８に格納する。ステップＳ１１１では更に、差分値積算回路１５３が、データ記憶装置１８に格納された標本値Ｕj+2を読み出してステップＳ１１２に進む。信号処理回路１５の動作は、測定と同時にリアルタイムで進行するので、ステップＳ１１１におけるデータ記憶装置１８に格納された標本値Ｕj+2を差分値積算回路１５３が読み出す処理は、データ記憶装置１８を介さずに、第１波形又は第２波形が測定されたタイミングに対応して、ＡＤ変換器１４から差分値積算回路１５３に標本値Ｕj+1が直接取り込まれるようにできる。
　ステップＳ１０６で、差分値ΔＵj+1,jが立ち上がり期間の特徴量の下限識別値ＬＬＤ（Ｕ）より大きくないと判断された場合は、ステップＳ１０７に進む。ステップＳ１０７では、データ記憶装置１８に格納された次の標本値Uj+2を新たな標本値Ｕj+1に置き換え、ステップＳ１０８に進む。ステップＳ１０８では、データ記憶装置１８に格納された標本値Ｕj+1を標本値Ｕjに置き換え、この新たな標本値Ｕjを演算進行判定回路１５７に取り込ませて、ステップＳ１０４に戻る。
　ステップＳ１１２で、差分値積算回路１５３が、データ記憶装置１８に格納された標本値Ｕj+1を読み出し、差分値ΔＵj+2,j+1＝Ｕj+2－Ｕj+1を計算し、計算結果を演算進行判定回路１５７に出力する。ステップＳ１１２において、演算進行判定回路１５７は、データ記憶装置１８に格納された差分値ΔＵj+1,jを読み出し、差分値積算回路１５３が出力した差分値ΔＵj+2,j+1が、差分値ΔＵj+1,jより大きいか、又は差分値ΔＵj+2,j+1が正の値であるか否かを判定する。ステップＳ１１２で、差分値ΔＵj+2,j+1が差分値ΔＵj+1,jより大きい又は差分値ΔＵj+2,j+1が正の値であるのいずれかの条件を満足した場合は、差分値ΔＵj+2,j+1を信号処理回路１５の差分値積算回路１５３に出力して、ステップＳ１１３に進む。一方、ステップＳ１１２で、差分値ΔＵj+2,j+1が差分値ΔＵj+1,jより大きい又は差分値ΔＵj+2,j+1が正の値であるのいずれの条件をも満足しない場合は、差分値ΔＵj+2,j+1及び差分値ΔＵj+1,jを、差分値積算回路１５３に並列又は順次出力して、ステップＳ１２１に進む。
　ステップＳ１１３では、信号処理回路１５の差分値積算回路１５３が、特徴量Ｕｓ及び差分値ΔＵj+1,jをデータ記憶装置１８から読み出し、Ｕｓ＋ΔＵj+1,j＋ΔＵj+2,j+1の値を計算し、計算結果を新たなＵｓとしてステップＳ１１４に進む。
　ステップＳ１１４では、差分値積算回路１５３が、新たな特徴量Ｕｓの値（＝Ｕｓ＋ΔＵj+1,j＋ΔＵj+2 ,j+1）及び標本値Ｕj+2をデータ記憶装置１８に格納する。ステップＳ１１４では、プログラムカウンタ１６１がプログラム記憶装置１９から次に読み出す命令のアドレスをｊ＋２からｊ＋１に戻し、更にデータ記憶装置１８に格納された次の標本値Ｕj+1のアドレスを新たな標本値Ｕjのアドレスに置き換え、新たな標本値Ｕj+1を演算進行判定回路１５７がデータ記憶装置１８から読み出し、ステップＳ１０６に戻る。
　ステップＳ１２１では、差分値積算回路１５３が、特徴量Ｕｓをデータ記憶装置１８から読み出し、Ｕｓ＋ΔＵj+1,j＋ΔＵj+2,j+1の値を第１座標軸の値Ｕｆとして計算し、ステップＳ１２２に進む。ステップＳ１２２では、信号処理回路１５のピーク値決定回路１６３が、データ記憶装置１８に格納された標本値Ｕj+1及び標本値Uj+2を読み出し、標本値Ｕj+1と標本値Ｕj+2の大きさを比較する。ピーク値決定回路１６３がＵj+2＞Ｕj+1と判断した場合は、標本値Ｕj+2の値をピーク値Ｕpと判断し、ピーク値Ｕp＝Ｕj+2の値と、差分値積算回路１５３が決定した第１座標軸の値Ｕｆをデータ記憶装置１８に格納し、ステップＳ２０１に進む。ピーク値決定回路１６３がＵj+2＜Ｕj+1と判断した場合は、標本値Ｕj+1の値をピーク値Ｕpと判断し、ピーク値Ｕp＝Ｕj+1の値をデータ記憶装置１８に格納する。更に、差分値積算回路１５３は、ステップＳ１２１で決定した第１座標軸の値ＵｆをＵｆ＝Ｕｓ＋ΔＵj+1,jと補正して、補正したＵｆの値をデータ記憶装置１８に格納し、ステップＳ２０１に進む。
　図１０のステップＳ２０１において、演算進行判定回路１５７が立ち下がり期間の特徴量Ｄｓをリセットし、特徴量Ｄｓ＝０の値をデータ記憶装置１８に格納する。その後、プログラムカウンタ１６１がプログラム記憶装置１９から次に読み出す命令のアドレスをカウントして、ステップＳ２０２に進む。ステップＳ２０２において、減衰量計算回路１５４が標本値Ｄjとピーク値Ｕpをデータ記憶装置１８から読み出し、減衰量Ｄdj＝Ｕp－Ｄjを計算し、減衰量Ｄdjを演算進行判定回路１５７に出力する。図８には、例としてｊ＝ｎの場合についての定義が示され、減衰量Ｄdn＝Ｕp－Ｄnが、ピーク値Ｕp＝Ｕmaxに対して計算されることが示されている。
　なお、信号処理回路１５の動作は、波形検出器１２による測定と同時にリアルタイムで進行するので、ステップＳ２０２におけるデータ記憶装置１８に格納された標本値Ｄjを減衰量計算回路１５４が読み出す処理は、データ記憶装置１８を介さずに、第１波形又は第２波形が測定されたタイミングに対応して、ＡＤ変換器１４から減衰量計算回路１５４に標本値Ｄjが直接取り込まれるようにできる。
　演算進行判定回路１５７は、ステップＳ２０２において、減衰量Ｄdjが立ち下がり期間の特徴量の下限識別値ＬＬＤ（Ｄ）より大きいか否かを判定する。ステップＳ２０２で、減衰量Ｄdjが立ち下がり期間の特徴量の下限識別値ＬＬＤ（Ｄ）より大きいと判断された場合は、ステップＳ２０３に進み、減衰量Ｄdjをデータ記憶装置１８に格納する。
　ステップＳ２０３では更に、減衰量計算回路１５４が標本値Ｄj+1とピーク値Ｕpをデータ記憶装置１８から読み出し、減衰量Ｄdj+1＝Ｕp－Ｄj+1を計算し、ステップＳ２０４に進む。ステップＳ２０３におけるデータ記憶装置１８に格納された標本値Ｄj+1を減衰量計算回路１５４が読み出す処理は、データ記憶装置１８を介さずに、第１波形又は第２波形が測定されたタイミングに対応して、ＡＤ変換器１４から減衰量計算回路１５４に標本値Ｄj+1が直接取り込まれるようにできる。
　ステップＳ２０２で、減衰量Ｄdjが立ち下がり期間の特徴量の下限識別値ＬＬＤ（Ｄ）より大きくないと判断された場合は、ステップＳ２０６に進む。ステップＳ２０６では、データ記憶装置１８に格納された次の標本値Ｄj+1を新たな標本値Ｄjに置き換え、この新たな標本値Ｄjを減衰量計算回路１５４が取り込み、信号処理回路１５の処理はステップＳ２０２に戻る。
　ステップＳ２０４では差分値積算回路１５３が、データ記憶装置１８に格納された減衰量Ｄdjを読み出し、減衰量の差分値ΔＤj+1,j＝Ｄdj+1－Ｄdjを計算し、計算結果を演算進行判定回路１５７に出力する。ステップＳ２０４において、演算進行判定回路１５７は、減衰量の差分値ΔＤj+1,jが立ち下がり期間の特徴量の下限識別値ＬＬＤ（Ｄ）より大きいか否か、又は減衰量Ｄdj+1が減衰量Ｄdjより大きいか否か、を判定する。
　ステップＳ２０４で、減衰量の差分値ΔＤj+1,jが立ち下がり期間の特徴量の下限識別値ＬＬＤ（Ｄ）より大きいと判断、又は減衰量Ｄdj+1が減衰量Ｄdjより大きいと判断された場合は、ステップＳ２１１に進み、減衰量Ｄdj+1及び減衰量の差分値ΔＤj+1,jをデータ記憶装置１８に格納する。ステップＳ２１１では更に、減衰量計算回路１５４が標本値Ｄj+2とピーク値Ｄpをデータ記憶装置１８から読み出し、減衰量Ｄdj+2＝Ｕp－Ｄj+2を計算し、ステップＳ２１２に進む。ステップＳ２１１におけるデータ記憶装置１８に格納された標本値Ｄj+2を減衰量計算回路１５４が読み出す処理は、第１波形又は第２波形が測定されたタイミングでＡＤ変換器１４から減衰量計算回路１５４に標本値Ｄj+2が、データ記憶装置１８を介さずに直接取り込まれるようにすることもできる。
　ステップＳ２０４で、減衰量の差分値ΔＤj+1,jが立ち下がり期間の特徴量の下限識別値ＬＬＤ（Ｄ）より大きくないと判断、又は減衰量Ｄdj+1が減衰量Ｄdjより大きくないと判断された場合は、ステップＳ２０５に進む。ステップＳ２０５では、データ記憶装置１８に格納された次の標本値Ｄj+２を新たな標本値Ｄj+1に置き換え、ステップＳ２０６に進む。ステップＳ２０６では、データ記憶装置１８に格納された標本値Ｄj+1を標本値Ｄjに置き換え、この新たな標本値Ｄjを減衰量計算回路１５４に取り込ませて、ステップＳ２０２に戻る。
　ステップＳ２１２では、差分値積算回路１５３が、減衰量の差分値ΔＤj+2,j＝Ｄdj+2－Ｄdj+1を計算し、減衰量の差分値ΔＤj+2,j+1が減衰量の差分値ΔＤj+1,jより大きいか否か、又は減衰量Ｄdj+2が減衰量Ｄdj+1より大きいか否かを判定する。ステップＳ２１２で、減衰量の差分値ΔＤj+2,j+1が減衰量の差分値ΔＤj+1,jより大きい又は減衰量Ｄdj+2が減衰量Ｄdj+1より大きいか、のいずれかの条件を満足した場合は、減衰量の差分値ΔＤj+2,j+1を差分値積算回路１５３に出力して、ステップＳ２１３に進む。
　一方、ステップＳ２１２で、減衰量の差分値ΔＤj+2,j+1が減衰量の差分値ΔＤj+1,jより大きい又は減衰量Ｄdj+2が減衰量Ｄdj+1より大きいか、のいずれの条件をも満足しない場合は、減衰量の差分値ΔＤj+2,j+1及びΔＤj+1,jを差分値積算回路１５３に並列又は順次出力して、ステップＳ２２１に進む。
　ステップＳ２１３では、差分値積算回路１５３が、特徴量Ｄｓ及び減衰量の差分値ΔＤj+1,jをデータ記憶装置１８から読み出し、Ｄｓ＋ΔＤj+1,j＋ΔＤj+2,j+1の値を計算し、計算結果を新たなＤｓとしてステップＳ２１４に進む。ステップＳ２１４では、差分値積算回路１５３が、新たな特徴量Ｄｓの値（＝Ｄｓ＋ΔＤj+1,j＋ΔＤj+2,j+1）及び減衰量Ｄdj+2をデータ記憶装置１８に格納する。ステップＳ２１４では、プログラムカウンタ１６１がプログラム記憶装置１９から次に読み出す命令のアドレスをｊ＋２からｊ＋１に戻し、更にデータ記憶装置１８に格納された次の減衰量Ｄdj+1のアドレスを新たな減衰量Ｄdjのアドレスに置き換え、新たな標本値Ｄj+1を減衰量計算回路１５４がデータ記憶装置１８から読み出し、ステップＳ２０４に戻る。
　ステップＳ２２１では、差分値積算回路１５３が、特徴量Ｄｓをデータ記憶装置１８から読み出し、Ｄｓ＋ΔＤj+1,j＋ΔＤj+2,j+1の値を第２座標軸の値Ｄｆとして計算し、第２座標軸の値Ｄｆをデータ記憶装置１８に格納して、ステップＳ２２２に進む。
　ステップＳ２２２では、信号処理回路１５の２次元座標プロット回路１５６が、第２座標軸の値Ｄｆと第１座標軸の値Ｕｆとの組からなる座標（Ｕf，Ｄf）を示す点を、図１２に示したような、第１座標軸と第２座標軸が定義する弁別平面にプロットする。座標（Ｕf，Ｄf）が弁別平面にプロットされたら、ステップＳ１０３に戻り、立ち上がり期間の特徴量Ｕｓをリセットする。ステップＳ１０３に戻ると、プログラムカウンタ１６１が信号処理回路１５が、クロック信号に同期させて、図９及び図１０に示すフローチャートに示された処理を時々刻々を実行させる。次の第１のパルス群に含まれるパルス又は第２のパルス群に含まれるパルスが波形検出器１２に入力され、第１及び第２の電気信号の少なくとも一方が、波形検出器１２から弁別対象信号として出力されることにより、２次元座標プロット回路１５６が、第１座標軸の値Ｕｆと第２座標軸の値Ｄｆとの組からなる座標（Ｕf，Ｄf）を示す新たな点を、図１２に示したような、第１座標軸と第２座標軸が定義する弁別平面にプロットする。
　第１の実施の形態に係る波形弁別方法によれば、パルス波形の立ち下がり期間において信号強度がベースラインに落ちる前に、ガンマ線や中性子線の放射線－光変換素子１１に対する入力があり、図１４に示すように、光検出器１２ａから出力される電気信号の立ち下がり波形に対し、次のパルス信号の立ち上がり波形が重畳するパイルアップが発生した場合も正しい波高値を取ることができる。
　即ち、パイルアップが発生し、ステップＳ２１２で減衰量Ｄdj+2が減衰量Ｄdj+1より小さいと判断された場合は、ステップＳ２２１及びステップＳ２２２を経由してステップＳ１０３に戻る。図１４に示すようなパイルアップが発生した場合は、ステップＳ１０３において、演算進行判定回路１５７が立ち上がり期間の特徴量Ｕｓをリセットし（Ｕｓ＝０）、ステップＳ１０４に信号処理回路１５の処理が進むことによりパイルアップ箇所の波形が測定できる。
　図１４では、立ち下がり波形に対し、２箇所でパイルアップが発生した場合が示されているが、パイルアップが発生すると、その都度、ステップＳ２１２で、パイルアップが発生したと判断できるので、ステップＳ２２１及びステップＳ２２２を経由してステップＳ１０３に戻り、ステップＳ１０３において、演算進行判定回路１５７が立ち上がり期間の特徴量Ｕｓをリセットし（Ｕｓ＝０）、ステップＳ１０４以降の一連のステップが実行されることにより、連続したパイルアップが発生した場合であっても正しい波形を計測することができる。
　既に述べたとおり、第１の実施の形態に係る波形弁別装置の特徴の一つは、図１に示したアナログ増幅器１３が、波形検出器１２から出力される弁別対象信号の過渡応答波形を時間軸に沿って拡大するように増幅することにある。弁別対象信号の立ち下がり時間を時間軸に沿って拡大することにより、ＡＤ変換器１４がデジタルデータを取得するためのサンプリング間隔を長くすることができるので、第１の実施の形態に係る波形弁別装置によれば、非常に安価かつ簡単なＡＤ変換器１４が採用可能となる。しかしながら、立ち下がり時間を時間軸に沿ってあまり拡大しすぎると、被測定パルスの物理量の特性によっては、図１４に示すようなパイルアップの確率が高くなり、パイルアップの間隔が短くなりすぎて、波形弁別に必要な標本値がとれなくなり精度が低下する恐れもある。したがって、図２に示したアナログ増幅器１３の入力端子Ｉと接地間に接続する入力抵抗Ｒ１の値は、５ｋΩ～１ＭΩ程度範囲内で、被測定パルスの物理量の特性に合わせて適宜選択し、最適値に調整すればよい。この入力抵抗Ｒ１の値の調整は、図３に示した筐体２１の底面に設けられた調整つまみ３４ａ，３４ｂ，３４ｃ，３４ｄと同様な入力抵抗調整つまみを更に設け、入力抵抗Ｒ１の値を可変にして、被測定パルスの物理量の特性をみながら調整するようにして、汎用性を高めるようにしてもよい。
　そして、第１の実施の形態に係る波形弁別方法を実現するための波形弁別装置の構成は、図１～図３及び図１３に示したような簡単かつ安価なハードウェア資源を基礎としているので、結果として、測定に必要な費用も安価にできる。又、測定に用いる波形弁別装置を小型の回路基板上に集積化してポータブル化が容易であるので、作業性が向上するという顕著な効果を奏することができる。
（パルスの波形の弁別）
　信号処理回路１５の波形点累積回路１５９は、図１０のステップＳ２２２から図９のステップＳ１０３に戻る帰還ループを、信号処理回路１５を駆動する電源が入っている限り、継続的に繰り返す。図９及び図１０に示す一連のフローに沿ったループの繰り返しに応じて、新たな点が弁別平面上に逐次累積されるので、弁別平面には多数の座標点が、第１のパルス群に含まれるパルスの波形であるのか、又は第２のパルス群に含まれるパルスの波形であるのかに依拠して、局在してプロットされる。図１２に示すように、弁別平面上に複数の座標点が局在した領域に分布するので、図１５に示すフローチャートにしたがって、局在した領域を分類し、解析することにより、第１のパルス群に含まれるパルスの波形であるのか、第２のパルス群に含まれるパルスの波形であるのかが、弁別できる。　
　先ず、図１５のステップＳ３０１において、図１３に示した信号処理回路１５の波形弁別判定回路１５８が座標点の位置が弁別用の窓部の内部に位置するか否かを判定する。弁別用の窓部は図１２においては、第２座標軸であるＸ軸に沿って立ち下がり期間の特徴量Ｄｆの下限識別値ＬＬＤ（Ｄ）と立ち下がり期間の特徴量Ｄｆの上限識別値ＵＬＤ（Ｄ）が定められ、第１座標軸であるＹ軸に沿って立ち上がり期間の特徴量Ｕｆの下限識別値ＬＬＤ（Ｕ）と立ち上がり期間の特徴量Ｕｆの上限識別値ＵＬＤ（Ｕ）が定められる。これらの下限識別値ＬＬＤ（Ｄ）、上限識別値ＵＬＤ（Ｄ）、下限識別値ＬＬＤ（Ｕ）及び上限識別値ＵＬＤ（Ｕ）は、それぞれ図１６に示す手順で予め決定しておき、それぞれデータ記憶装置１８に格納しておき、窓部の位置を決定する際にデータ記憶装置１８から読み出せばよい。即ち、図１２においては、データ記憶装置１８に格納されたデータをそれぞれ用いて、第２座標軸に垂直で、第２座標軸に対する切片がＬＬＤ（Ｄ）とＵＬＤ（Ｄ）の２本の平行な直線（縦線）と、第１座標軸に垂直で、第１座標軸に対する切片が値ＬＬＤ（Ｕ）とＵＬＤ（Ｕ）である２本の平行な直線（水平線）で囲まれる矩形領域で弁別用の窓部が定義される。
　ステップＳ３０１においては、第１座標軸の値Ｕｆと第２座標軸の値Ｄｆとの組で定義される座標点（Ｕf，Ｄf）の分布が、弁別用の窓部の内部に位置するか否かを波形弁別判定回路１５８が判断する。ステップＳ３０１で座標点（Ｕf，Ｄf）の分布が、弁別用の窓部の内部に位置しないと判断された場合は、波形弁別判定回路１５８は、ステップＳ３０４において、波形検出器１２から出力された弁別対象信号が第１波形を発生源とする信号であると判断する。一方、ステップＳ３０１において、座標点（Ｕf，Ｄf）の分布が、弁別用の窓部の内部に位置すると判断された場合は、ステップＳ３０２に進む。
　ステップＳ３０２においては、第１座標軸の値Ｕｆと第２座標軸の値Ｄｆとの組で定義される座標点（Ｕf，Ｄf）の分布が、弁別用線形方程式を表す直線よりも第２座標軸側の領域に存在するか否かを波形弁別判定回路１５８が判断する。弁別用線形方程式は図１２に示すように傾きａ、第１座標軸の切片ｂの１次関数で表現される。
　これらの弁別用線形方程式の傾きａ及び切片ｂの値は、それぞれ、図１６に示す手順で予め決定しておき、データ記憶装置１８に格納しておき、窓部の位置を決定する際にデータ記憶装置１８から読み出せばよい。即ち、図１２においては、データ記憶装置１８に格納された傾きａ及び切片ｂの値を用いて、弁別用線形方程式が弁別平面上に定義される。
　ステップＳ３０１で座標点（Ｕf，Ｄf）の分布が、弁別用線形方程式を表す直線よりも第２座標軸側に位置しないと判断された場合は、波形弁別判定回路１５８は、ステップＳ３０４において、波形検出器１２から出力された弁別対象信号が第１波形を発生源とする信号であると判断する。一方、ステップＳ３０１において、座標点（Ｕf，Ｄf）の分布が、弁別用線形方程式を表す直線よりも第２座標軸側に位置すると判断された場合は、ステップＳ３０３に進み、波形弁別判定回路１５８は、波形検出器１２から出力された弁別対象信号が第２波形を発生源とする信号であると判断する。
　このように、図１５に示すフローチャートにしたがって、第１のパルス群に含まれるパルスの波形であるのか、第２のパルス群に含まれるパルスの波形であるのかを座標点（Ｕf，Ｄf）の分布位置によって弁別することにより、波形点累積回路１５９は、第１波形に対応した座標の累積数及び第１波形に対応した座標の累積数をそれぞれ計数することができる。
　波形点累積回路１５９が累積し、計数した第１波形及び第２波形にそれぞれ対応した座標の累積数は、信号処理回路１５の累積数表示命令回路１６０が、図１及び図３に示す表示装置１６に対し、表示命令及び表示に必要なデータを送信することによって、表示装置１６に表示させることができる。
（弁別用窓部及び弁別用線形方程式の決定）
　本発明者らは、一例としてガンマ線と中性子線の波形を互いに分離する第１の実施の形態に係る波形弁別方法の応用において、図１１に示すように、ガンマ線の立ち上がり特徴量Ｕfとガンマ線の立ち下がり特徴量Ｄfとが線形比例することを見いだした。一方、図１１に示すように、中性子線の立ち上がり特徴量Ｕfとガンマ線の立ち下がり特徴量Ｄfとの間にも、ガンマ線の場合に比すと弱いが、その座標点の分布領域のトポロジーには同様な線形比例の関係が認められることを見いだした。ガンマ線の立ち上がり特徴量Ｕfと立ち下がり特徴量Ｄfとの強い線形比例関係を、予め弁別用線形方程式Ｕf＝ａＤf＋ｂとして求めておくことにより、第１波形と第２波形とが正確に弁別できる。
　先ず、図１６のステップＳ４０１において、波形が既知の校正用となる第２のパルス群に含まれるパルスを波形検出器１２に入力させる。波形検出器１２から出力される第２の電気信号は逐次、波形検出器１２から弁別対象信号として出力され、アナログ増幅器１３が弁別対象信号の過渡応答波形を時間軸に沿って拡大し、ＡＤ変換器１４が増幅された弁別対象信号を標本化し、デジタルデータに変換する。ステップＳ４０１では、図１３に示した信号処理回路１５の窓部境界条件決定回路１５１に、このようにして校正用第２波形を発生源とするデジタルデータがリアルタイムで逐次、複数個入力され、複数の校正用第２波形が計測されると、ステップＳ４０２に進む。
　ステップＳ４０２において、窓部境界条件決定回路１５１は、複数の校正用第２波形にそれぞれ対応して波形検出器１２から出力された複数の第２の電気信号について、その立ち上がり期間のピーク値を、ＡＤ変換器１４が逐次変換したデジタルデータを用いて、統計処理により探索し、ステップＳ４０３に進む。
　ステップＳ４０２において探索された立ち上がり期間のピーク値を用いて、ステップＳ４０３では、窓部境界条件決定回路１５１が、立ち下がり期間の特徴量Ｄｆの下限識別値ＬＬＤ（Ｄ）、立ち下がり期間の特徴量Ｄｆの上限識別値ＵＬＤ（Ｄ）、立ち上がり期間の特徴量Ｕｆの下限識別値ＬＬＤ（Ｕ）及び立ち上がり期間の特徴量Ｕｆの上限識別値ＵＬＤ（Ｕ）を決定する。ステップＳ４０３において決定されたＬＬＤ（Ｄ），ＵＬＤ（Ｄ），ＬＬＤ（Ｕ）及びＵＬＤ（Ｕ）の値は、ステップＳ４０４において、データ記憶装置１８に格納される。
　その後、プログラムカウンタ１６１によって信号処理回路１５の処理がステップＳ４１１に進む。ステップＳ４１１において、波形が既知の校正用となる第１のパルス群に含まれるパルスを波形検出器１２に入力させ、複数の校正用第１波形が計測される。ステップＳ４１２において、波形検出器１２から出力される第１の電気信号は逐次、波形検出器１２から弁別対象信号として出力され、アナログ増幅器１３が弁別対象信号の過渡応答波形を時間軸に沿って拡大し、ＡＤ変換器１４が増幅された弁別対象信号を標本化し、デジタルデータに変換する。そして図９及び図１０に示したフローチャートにしたがって、ステップＳ４１２において、立ち上がり特徴量Ｕf、立ち下がり特徴量Ｄfがそれぞれ算出される。　
　更に、図９及び図１０に示したフローチャートにしたがって、ステップＳ４１３において、座標点（立ち上がり特徴量Ｕf、立ち下がり特徴量Ｄf）がそれぞれ算出され、図１１に示したのと同様に、複数個の座標点が弁別平面にプロットされる。その後、プログラムカウンタ１６１によって信号処理回路１５の処理がステップＳ４１４に進む。ステップＳ４１４において、信号処理回路１５の線形方程式決定回路１５２が、弁別平面上にプロットされた座標点の分布から弁別用線形方程式Ｕ＝ａＤ＋ｂの平均傾きａを算出する。
　その後、プログラムカウンタ１６１によって信号処理回路１５の処理がステップＳ４１５に進む。ステップＳ４１５において、線形方程式決定回路１５２が、弁別用線形方程式の切片ｂを決定する。線形方程式決定回路１５２は、ステップＳ４１６において、弁別用線形方程式Ｕ＝ａＤ＋ｂの平均傾きａ及び切片ｂの値をデータ記憶装置１８に格納する。
（波形弁別プログラム）
　図９、図１０，図１５及び図１６に示した一連の波形弁別の操作は、図９、図１０，図１５及び図１６と等価なアルゴリズムを実行させるプログラムにより、図１に示した波形弁別装置を制御して実行できる。この波形弁別プログラムは、図１に示したプログラム記憶装置１９に記憶させればよい。又、この波形弁別プログラムは、コンピュータ読取り可能な記録媒体に保存し、この記録媒体をプログラム記憶装置１９に読み込ませることにより、第１の実施の形態に係る一連の波形弁別の操作を実行することができる。
　ここで、「コンピュータ読取り可能な記録媒体」とは、例えばマイクロプロセッサの外部メモリ装置、半導体メモリ、磁気ディスク、光ディスク、光磁気ディスク、磁気テープなどの種々のプログラムを記録することができるような媒体であれば構わない。具体的には、フレキシブルディスク、ＣＤ－ＲＯＭ，ＭＯディスク、カセットテープ、オープンリールテープなどが「コンピュータ読取り可能な記録媒体」に含まれる。
　即ち、第１の実施の形態に係る波形弁別プログラムは：
　(a)波形検出器１２に被測定パルスの波形を入力させ、被測定パルスの物理量を電気信号に変換させる命令；
　(b)アナログ増幅器１３に電気信号の過渡応答波形を時間軸に沿って拡大して増幅させる命令；
　(c)ＡＤ変換器１４に電気信号の立ち上がり期間及び立ち下がり期間において、増幅された電気信号を標本化してデジタルデータに変換させる命令；
　(d)信号処理回路１５の差分値計算回路１５３、減衰量計算回路１５４及び差分値積算回路１５５を互いに連携させ、デジタルデータを用いて、立ち上がり期間の特徴量Ｕfを第１座標軸上の点として計算させ、立ち下がり期間の特徴量Ｄfを第２座標軸上の点として計算させる命令；
　(e)信号処理回路１５の２次元座標プロット回路１５６に、第１座標軸上の点及び第２座標軸上の点の組を座標点として、第１座標軸と第２座標軸が定義する弁別平面に、座標点をプロットさせる命令；
　(f)信号処理回路１５の波形弁別判定回路１５８に、座標点のプロット位置から、被測定パルスが第１波形、或いは第１波形とは異なる第２波形であるのか弁別させる命令等を含む一連の命令を、図１に示した制御回路１７に実行させる波形弁別プログラムである。
　第１の実施の形態に係る波形弁別装置の制御回路１７や信号処理回路１５は、例えばフレキシブルディスク装置（フレキシブルディスクドライブ）及び光ディスク装置（光ディスクドライブ）を内蔵若しくは外部接続するように構成できる。フレキシブルディスクドライブに対してはフレキシブルディスクを、又光ディスクドライブに対してはＣＤ－ＲＯＭをその挿入口から挿入し、所定の読み出し操作を行うことにより、これらの記録媒体に格納された波形弁別プログラムを波形弁別装置を構成するプログラム記憶装置１９にインストールすることができる。更に、インターネット等の情報処理ネットワークを介して、この波形弁別プログラムをプログラム記憶装置１９に格納することが可能である。
（その他の実施の形態）
　上記のように、本発明は第１の実施の形態によって記載したが、この開示の一部をなす論述及び図面は本発明を限定するものであると理解すべきではない。この開示から当業者には様々な代替実施の形態、実施例及び運用技術が明らかとなろう。
　既に述べた第１の実施の形態の説明においては、第１波形がガンマ線に固有な放射線－光変換素子１１からの発光波形、第２波形が中性子線に固有な放射線－光変換素子１１からの発光波形の場合において、波形検出器１２が第１波形の光パルスを入力して第１の電気信号を出力し、第２波形の光パルスを入力して第２の電気信号を出力する光検出器である場合について例示的に説明したが、第１の実施の形態の説明に限定されるものではない。例えば、波形検出器１２が第１波形を有する音波を入力して第１の電気信号を出力し、第２波形を有する音波を入力して第２の電気信号を出力する音響電気変換素子であっても構わない。
　このように、本発明はここでは記載していない様々な実施の形態等を含むことは勿論である。したがって、本発明の技術的範囲は上記の説明から妥当な特許請求の範囲に係る発明特定事項によってのみ定められるものである。
　本発明は、互いに波形の異なる２つのパルス波形を弁別する波形弁別装置、波形弁別方法及び波形弁別プログラムであって、例えば、原子力発電等に用いられる自然界には存在しない放射性物質から生じるガンマ線と中性子を正確に分離するために利用することができる。又、例えば超音波探傷において伝播時間による計測では見いだせない異物質によるエコーを互いに波形の異なる２つのエコーパルス波形を弁別する手段を提供して明確に分離する産業上の利用価値がある。
　１１…光変換素子
　１２…波形検出器
　１２ａ…光検出器
　１３…アナログ増幅器
　１４…ＡＤ変換器
　１５…信号処理回路
　１６…表示装置
　１７…制御回路
　１８…データ記憶装置
　１９…プログラム記憶装置
　２１…筐体
　２２…高圧電源
　２３…回路基板
　２４…回路基板
　３１ａ，３１ｂ，３２ａ，３２ｂ，３２ｃ…ケーブル
　３３…通信用ケーブル
１５１…窓部境界条件決定回路
１５２…線形方程式決定回路
１５３…差分値積算回路
１５４…減衰量計算回路
１５５…差分値積算回路
１５６…２次元座標プロット回路
１５７…演算進行判定回路
１５８…波形弁別判定回路
１５９…波形点累積回路
１６０…累積数表示命令回路
１６１…プログラムカウンタ
１６２…データ取込回路
１６３…ピーク値決定回路
１６４…データバス

Claims (26)

  1. 　被測定パルスの波形を入力し、前記被測定パルスの物理量を電気信号に変換する波形検出器と、
    　前記電気信号の過渡応答波形を時間軸に沿って拡大して増幅するアナログ増幅器と、
    　前記電気信号の立ち上がり期間及び立ち下がり期間において、増幅された前記電気信号を標本化してデジタルデータに変換するＡＤ変換器と、
    　前記デジタルデータを用いて、前記立ち上がり期間の特徴量を第１座標軸上の点として計算し、前記立ち下がり期間の特徴量を第２座標軸上の点として計算し、更に、前記第１座標軸上の点及び前記第２座標軸上の点の組を座標点として、前記第１座標軸と前記第２座標軸が定義する弁別平面に前記座標点をプロットする信号処理回路と、
    　を備え、前記座標点のプロット位置から、前記被測定パルスが第１波形、或いは前記第１波形とは異なる第２波形であるのか弁別することを特徴とする波形弁別装置。
  2. 　中性子線とガンマ線とで発光特性が異なり、前記中性子線及び前記ガンマ線を光に変換する放射線－光変換素子を更に含み、
    　前記波形検出器は、前記光を電気信号に変換する光検出器であることを特徴とする請求項１に記載の波形弁別装置。
  3. 　前記放射線－光変換素子は、ＣｓＬｉＹＣｌ，ＬｉＣａＡｌＦ6，ＬｉＦ／ＺｎＳ，ＬｉＢａＦ3，Ｌｉ6Ｇｄ（ＢＯ3）3のいずれかを材料とするシンチレータであることを特徴とする請求項２に記載の波形弁別装置。
  4. 　前記波形検出器は、波長１９０～４５０ｎｍの光を電気信号に変換する光検出器でることを特徴とする請求項２又は３に記載の波形弁別装置。
  5. 　前記光検出器は、光電子増倍管、半導体フォトダイオード、フォトダイオードアレイ、ガイガーモード並列読み出しＡＰＤピクセルアレイのいずれかであることを特徴とする請求項２～４のいずれか１項に記載の波形弁別装置。
  6. 　前記光検出器は光電子増倍管であり、該光電子増倍管の信号出力端子と基準電位点端子が、各々前記アナログ増幅器の入力端子と接地間に接続され、更にアナログ増幅器の入力端子と接地間に５ｋΩ以上の入力抵抗が接続されていることを特徴とする請求項５に記載の波形弁別装置。
  7. 　前記アナログ増幅器は、前記電気信号の立ち下がり時間が２μ秒以上になるように、前記電気信号の過渡応答波形を時間軸に沿って拡大することを特徴とする請求項１～６のいずれか１項に記載の波形弁別装置。
  8. 　前記信号処理回路は、予め波形が既知の校正用物理量を前記波形検出器に入力することによって、前記弁別平面上に定めた弁別用の窓部の内部に前記座標点のプロット位置が存在するか否かを判定する波形弁別判定回路を含むことを特徴とする請求項７に記載の波形弁別装置。
  9. 　前記波形弁別判定回路は、立ち下がり期間の特徴量の下限識別値、立ち下がり期間の特徴量の上限識別値、立ち上がり期間の特徴量の下限識別値及び立ち上がり期間の特徴量の上限識別値で囲まれた矩形の領域を前記弁別用の窓部とすることを特徴とする請求項８に記載の波形弁別装置。
  10. 　前記波形弁別判定回路は、前記弁別用の窓部の内部に前記座標点のプロット位置が存在しない判定された場合、前記被測定パルスが前記第１波形であると判別し、
    　前記弁別用の窓部の内部に前記座標点のプロット位置が存在すると判定された場合、前記被測定パルスが弁別用線形方程式を表す直線よりも第２座標軸側の領域に存在するか否かを判定することを特徴とする請求項８又は９に記載の波形弁別装置。
  11. 　前記信号処理回路は、前記立ち上がり期間の互いに連続する２つの前記デジタルデータの差分値を計算する差分値計算回路を更に含むことを特徴とする請求項９～１０のいずれか１項に記載の波形弁別装置。
  12. 　前記信号処理回路は、前記差分値を積算して前記立ち上がり期間の特徴量を決定計算する差分値積算回路を更に含むことを特徴とする請求項１１に記載の波形弁別装置。
  13. 　前記信号処理回路は、前記立ち上がり期間のピーク値と前記立ち下がり期間の前記デジタルデータの差分により、前記立ち下がり期間の減衰量を計算する減衰量計算回路を更に含むこことを特徴とする請求項１２に記載の波形弁別装置。
  14. 　前記差分値計算回路は、前記立ち下がり期間の互いに連続する２つの前記減衰量の差分値を計算することを特徴とする請求項１３に記載の波形弁別装置。
  15. 　前記差分値積算回路は、前記減衰量の差分値を積算して前記立ち下がり期間の特徴量を決定することを特徴とする請求項１４に記載の波形弁別装置。
  16. 　被測定パルスの波形を入力し、前記被測定パルスの物理量を電気信号に変換するステップと、
    　前記電気信号の過渡応答波形を時間軸に沿って拡大して増幅するステップと、
    　前記電気信号の立ち上がり期間及び立ち下がり期間において、増幅された前記電気信号を標本化してデジタルデータに変換するステップと、
    　前記デジタルデータを用いて、前記立ち上がり期間の特徴量を第１座標軸上の点として計算し、前記立ち下がり期間の特徴量を第２座標軸上の点として計算するステップと、
    　前記第１座標軸上の点及び前記第２座標軸上の点の組を座標点として、前記第１座標軸と前記第２座標軸が定義する弁別平面に、前記座標点をプロットするステップと、
    　前記座標点のプロット位置から、前記被測定パルスが第１波形、或いは前記第１波形とは異なる第２波形であるのか弁別するステップと、
    　を含むことを特徴とする波形弁別方法。
  17. 　前記弁別するステップでは、予め波形が既知の校正用物理量を入力して測定することによって、前記弁別平面上に定めた弁別用の窓部の内部に前記座標点のプロット位置が存在するか否かを判定することにより、前記第１波形と前記第２波形が弁別されることを特徴とする請求項１６に記載の波形弁別方法。
  18. 　前記弁別用の窓部は立ち下がり期間の特徴量の下限識別値、立ち下がり期間の特徴量の上限識別値、立ち上がり期間の特徴量の下限識別値及び立ち上がり期間の特徴量の上限識別値で囲まれた矩形の領域であることを特徴とする請求項１７に記載の波形弁別方法。
  19. 　前記弁別用の窓部の内部に前記座標点のプロット位置が存在しない判定された場合、前記被測定パルスが前記第１波形であると判別し、
    　前記弁別用の窓部の内部に前記座標点のプロット位置が存在すると判定された場合、前記被測定パルスが弁別用線形方程式を表す直線よりも第２座標軸側の領域に存在するか否かを判定することを特徴とする請求項１７又は１８に記載の波形弁別方法。
  20. 　前記立ち上がり期間の特徴量を計算するステップは、前記立ち上がり期間の互いに連続する２つの前記デジタルデータの差分値を計算する段階を含むことを特徴とする請求項１６～１９のいずれか１項に記載の波形弁別方法。
  21. 　前記立ち上がり期間の特徴量を計算するステップは、前記差分値を積算して前記立ち上がり期間の特徴量を決定する段階を更に含むことを特徴とする請求項２０に記載の波形弁別方法。
  22. 　互いに連続する２つの前記差分値を比較することにより、前記立ち上がり期間のピーク値を決定することを特徴とする請求項２１に記載の波形弁別方法。
  23. 　前記立ち下がり期間の特徴量を計算するステップは、前記立ち上がり期間のピーク値と前記立ち下がり期間の前記デジタルデータとの差分により、前記立ち下がり期間の減衰量を計算する段階を含むことを特徴とする請求項２２に記載の波形弁別方法。
  24. 　前記立ち下がり期間の特徴量を計算するステップは、前記立ち下がり期間の互いに連続する２つの前記減衰量の差分値を計算する段階を含むことを特徴とする請求項２３に記載の波形弁別方法。
  25. 　前記立ち下がり期間の特徴量を計算するステップは、前記減衰量の差分値を積算して前記立ち下がり期間の特徴量を決定する段階を更に含むことを特徴とする請求項２４に記載の波形弁別方法。
  26. 　波形検出器に被測定パルスの波形を入力させ、前記被測定パルスの物理量を電気信号に変換させる命令と、
    　アナログ増幅器に前記電気信号の過渡応答波形を時間軸に沿って拡大して増幅させる命令と、
    　ＡＤ変換器に前記電気信号の立ち上がり期間及び立ち下がり期間において、増幅された前記電気信号を標本化してデジタルデータに変換させる命令と、
    　信号処理回路の差分値計算回路、減衰量計算回路及び差分値積算回路を互いに連携させ、前記デジタルデータを用いて、前記立ち上がり期間の特徴量を第１座標軸上の点として計算させ、前記立ち下がり期間の特徴量を第２座標軸上の点として計算させる命令と、
    　前記信号処理回路の２次元座標プロット回路に、前記第１座標軸上の点及び前記第２座標軸上の点の組を座標点として、前記第１座標軸と前記第２座標軸が定義する弁別平面に、前記座標点をプロットさせる命令と、
    　前記信号処理回路の波形弁別判定回路に、前記座標点のプロット位置から、前記被測定パルスが第１波形、或いは前記第１波形とは異なる第２波形であるのか弁別させる命令と、
    　を含む一連の命令を制御回路に実行させることを特徴とする波形弁別プログラム。
PCT/JP2014/001106 2014-02-28 2014-02-28 波形弁別装置、波形弁別方法及び波形弁別プログラム Ceased WO2015128905A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
PCT/JP2014/001106 WO2015128905A1 (ja) 2014-02-28 2014-02-28 波形弁別装置、波形弁別方法及び波形弁別プログラム
JP2016504864A JP6414830B2 (ja) 2014-02-28 2014-02-28 波形弁別装置、波形弁別方法及び波形弁別プログラム
US15/121,645 US20160356897A1 (en) 2014-02-28 2014-02-28 Waveform discrimination device, waveform discrimination method, and waveform discrimination program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2014/001106 WO2015128905A1 (ja) 2014-02-28 2014-02-28 波形弁別装置、波形弁別方法及び波形弁別プログラム

Publications (1)

Publication Number Publication Date
WO2015128905A1 true WO2015128905A1 (ja) 2015-09-03

Family

ID=54008279

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2014/001106 Ceased WO2015128905A1 (ja) 2014-02-28 2014-02-28 波形弁別装置、波形弁別方法及び波形弁別プログラム

Country Status (3)

Country Link
US (1) US20160356897A1 (ja)
JP (1) JP6414830B2 (ja)
WO (1) WO2015128905A1 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018235810A1 (ja) * 2017-06-23 2018-12-27 浜松ホトニクス株式会社 光検出器及び光検出装置
JP2022036679A (ja) * 2020-08-24 2022-03-08 株式会社東芝 放射線検出システムおよび放射線検出方法
JP2023013476A (ja) * 2021-07-16 2023-01-26 株式会社日立ハイテク 信号検出装置、測定装置、および質量分析装置
CN116819599A (zh) * 2022-12-26 2023-09-29 成都理工大学工程技术学院 一种中子-伽马射线甄别方法、系统、设备及介质

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6628701B2 (ja) * 2016-08-05 2020-01-15 三菱電機株式会社 放射線測定装置
JP2019052854A (ja) * 2017-09-12 2019-04-04 株式会社東芝 放射線検出装置および放射線検出方法
CN113835114B (zh) * 2021-08-25 2024-04-26 吉林大学 一种紧凑式高能伽玛射线反符合叠层探测器
CN117749141B (zh) * 2024-02-20 2024-05-07 成都工业学院 堆积脉冲信号识别与成形方法、计算机程序产品及终端

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57198878A (en) * 1981-05-30 1982-12-06 Toshiba Corp Wide-ranged measuring system for neutron
JPS6249282A (ja) * 1985-08-29 1987-03-03 Nippon Atom Ind Group Co Ltd 放射能濃度モニタ
JP2000258536A (ja) * 1999-03-05 2000-09-22 Commiss Energ Atom ガンマ線検出器からの検出結果をリアルタイムで選別するための方法および装置ならびに検出素子の一様化のための補正を行うための方法および装置
JP2008026195A (ja) * 2006-07-24 2008-02-07 Hitachi Ltd 中性子強度分布測定装置
JP2010133832A (ja) * 2008-12-04 2010-06-17 Hitachi Ltd 放射線モニタリング装置及び放射線モニタリングシステム
WO2011145651A1 (ja) * 2010-05-18 2011-11-24 株式会社トクヤマ 中性子線検出器、中性子線検出用シンチレーターおよび中性子線とγ線とを弁別する方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US2994781A (en) * 1959-08-27 1961-08-01 Atomic Energy Authority Uk Nuclear particle discriminators
US5347129A (en) * 1993-04-14 1994-09-13 University Of Missouri-Columbia System for determining the type of nuclear radiation from detector output pulse shape
US6953937B2 (en) * 2003-06-26 2005-10-11 Battelle Energy Alliance, Llc Method and apparatus for the detection of neutrons and gamma rays

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS57198878A (en) * 1981-05-30 1982-12-06 Toshiba Corp Wide-ranged measuring system for neutron
JPS6249282A (ja) * 1985-08-29 1987-03-03 Nippon Atom Ind Group Co Ltd 放射能濃度モニタ
JP2000258536A (ja) * 1999-03-05 2000-09-22 Commiss Energ Atom ガンマ線検出器からの検出結果をリアルタイムで選別するための方法および装置ならびに検出素子の一様化のための補正を行うための方法および装置
JP2008026195A (ja) * 2006-07-24 2008-02-07 Hitachi Ltd 中性子強度分布測定装置
JP2010133832A (ja) * 2008-12-04 2010-06-17 Hitachi Ltd 放射線モニタリング装置及び放射線モニタリングシステム
WO2011145651A1 (ja) * 2010-05-18 2011-11-24 株式会社トクヤマ 中性子線検出器、中性子線検出用シンチレーターおよび中性子線とγ線とを弁別する方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018235810A1 (ja) * 2017-06-23 2018-12-27 浜松ホトニクス株式会社 光検出器及び光検出装置
CN110785679A (zh) * 2017-06-23 2020-02-11 浜松光子学株式会社 光检测器和光检测装置
US11194061B2 (en) 2017-06-23 2021-12-07 Hamamatsu Photonics K.K. Optical detector and optical detection device
CN110785679B (zh) * 2017-06-23 2023-11-07 浜松光子学株式会社 光检测器和光检测装置
JP2022036679A (ja) * 2020-08-24 2022-03-08 株式会社東芝 放射線検出システムおよび放射線検出方法
JP7438893B2 (ja) 2020-08-24 2024-02-27 株式会社東芝 放射線検出システムおよび放射線検出方法
JP2023013476A (ja) * 2021-07-16 2023-01-26 株式会社日立ハイテク 信号検出装置、測定装置、および質量分析装置
JP7543221B2 (ja) 2021-07-16 2024-09-02 株式会社日立ハイテク 信号検出装置、測定装置、および質量分析装置
CN116819599A (zh) * 2022-12-26 2023-09-29 成都理工大学工程技术学院 一种中子-伽马射线甄别方法、系统、设备及介质
CN116819599B (zh) * 2022-12-26 2024-01-23 成都理工大学工程技术学院 一种中子-伽马射线甄别方法、系统、设备及介质

Also Published As

Publication number Publication date
JP6414830B2 (ja) 2018-10-31
US20160356897A1 (en) 2016-12-08
JPWO2015128905A1 (ja) 2017-03-30

Similar Documents

Publication Publication Date Title
JP6414830B2 (ja) 波形弁別装置、波形弁別方法及び波形弁別プログラム
JP5027124B2 (ja) パルス波形分析による単一トランスデューサ内同時発生放射線の検出のための方法及び装置
JP4160275B2 (ja) エネルギー測定方法及び測定装置
US6362472B1 (en) Method for calibrating a radiation detection system
JP4486623B2 (ja) コンプトン撮像カメラ
US20070210247A1 (en) Stabilization of a Scintillation Detector
JP2008111677A (ja) ＺｎＳ蛍光体を用いたシンチレーション中性子検出器の中性子／ガンマ線弁別方法
KR101864716B1 (ko) 알파·베타·감마 방사선 검출기능을 갖는 스마트 슬림형 플라스틱 섬광계수기
US2686266A (en) Improvement in radiation detectors
JP2005208057A (ja) ガンマ線検出器及びガンマ線撮像装置
CN110837098A (zh) 甄别脉冲信号的方法、fpga、装置和存储介质
Unzueta et al. Position sensitive alpha detector for an associated particle imaging system
CN116299640A (zh) 一种便携式宇生缪子探测系统
KR101330117B1 (ko) 다채널 광전소자를 이용하는 양전자방출 단층촬영장치
KR102280128B1 (ko) 방사성 핵종 판별 장치 및 방법
TW201809722A (zh) 計數方法以及放射線檢測裝置
KR102644122B1 (ko) 핵종판별을 위한 방사선 스펙트럼의 온도보상방법
JP2001194460A (ja) 放射線モニタ
JP6938239B2 (ja) 光検出器及び光検出装置
JP6818579B2 (ja) 土壌放射能汚染検査装置
Flaska et al. Optimization of an offline pulse-shape discrimination technique for the liquid scintillator BC-501A
JP7805472B2 (ja) 放射能分析装置及び放射能分析方法
USRE25634E (en) Nuclear borehole logging utilizing a
JP7438893B2 (ja) 放射線検出システムおよび放射線検出方法
RU215240U1 (ru) Чувствительный элемент блока механической калибровки

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2016504864

Country of ref document: JP

Kind code of ref document: A

WWE Wipo information: entry into national phase

Ref document number: 15121645

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 14883852

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 14883852

Country of ref document: EP

Kind code of ref document: A1