Abstract:
A method for improving the measurement accuracy of scintillation based radiation detectors for ionizing radiation pulses when the multi-channel analyzer of the radiation detector is exposed to continuously varying temperatures.

Description:
CROSS REFERENCE TO RELATED APPLICATIONS 
       [0001]    This application claims the priority of previously filed U.S. provisional 62/067,769 filed Oct. 23, 2014. 
     
    
     BACKGROUND 
       [0002]    Radiation detection devices are used in a variety of industrial, scientific, military, and government applications. Exemplary scintillator detectors have scintillator crystals made of activated sodium iodide or cesium iodide, or other materials that are effective for detecting gamma rays. 
         [0003]    Generally, the scintillator crystals are enclosed in casings or sleeves that include a window to permit radiation-induced scintillation light to pass out of the crystal package. The light passes to a light-sensing device such as a photomultiplier tube (PMT), and the photomultiplier tube converts the light photons emitted from the crystal into electrical pulses. The electrical pulses are shaped and digitized by associated electronics and may be registered as counts that are transmitted to analyzing equipment. 
         [0004]    Scintillators are used for nuclear and X-ray radiation detection. In response to a pulse of ionizing radiation they create a light flash. That light flash is recorded and analyzed by suitable instrumentation. 
         [0005]    Existing instruments convert the light pulse into a concomitant electronic pulse. Analog instruments will transform the electronic pulse to create a new pulse shape with a pulse-height that is proportional to the integral over the original pulse. Digital instruments will attempt to perform a direct integration of the electronic pulse without applying that transformation. 
         [0006]    The performance of both types of instruments degrades when the pulse shape of the scintillator light changes due to an external effect, such as a change in temperature. 
         [0007]    Multichannel analyzers are affected by a change in pulse shape in three distinct ways. 
         [0008]    First, there is a maximum electronic pulse height that can be processed by the electronics. It corresponds to the maximum radiation pulse energy (E_max). By the same token, there is a minimum electronic pulse height, corresponding to a minimal radiation pulse energy that can be detected (E_min). For practical applications it is mandatory to maintain a certain range of [E_min, E_max] over the entire range of measurement conditions, such as environmental temperature. If the pulse shapes change, E_min and E_max will be affected. For example, the pulses may lengthen even if the amount of light contained in the scintillator pulses does not change. In that case the resulting electronic pulses will have lower amplitudes, and small pulses that previously made it above the trigger threshold, will no longer be tall enough to trigger the data acquisition. 
         [0009]    Traditional analog or digital MCA&#39;s (multichannel analyzers) measure deposited energies using a fixed shaping time (analog MCA) or a fixed integration time (digital MCA). Both devices rely on the pulse shape being constant, except for statistical electronic noise fluctuations around the average pulse shape. If the pulse shape changes due to an external parameter such as temperature, traditional MCAs will not measure energies as accurately as possible. 
         [0010]    There is a second impact a changing pulse shape has on the accuracy of the energy measurement. In practice the shaping time or integration time is chosen such as to use about 90% of the light from the scintillation pulse for the measurement of the energy. Using less will deteriorate the precision of the measurement. If the scintillation pulses lengthen, for example due to a temperature change, they may become too long for the given operating parameters of the MCA. As a result, the MCA will measure energy values that are systematically smaller than the correct values. It will also lose accuracy of measurement (energy resolution) as it samples less than 90% of the scintillation light. 
         [0011]    The third way in which the pulse shape changes affect the performance of the MCA has to do with pile up rejection. An MCA will typically implement a method to recognize a set of two or more pulses that occur too close together in time to be measured independently and precisely. This situation is referred to as a pile up. Clearly, any measure of “too close in time” must be related to how long the pulses are. Hence, if the pulses lengthen, the MCA performance may deteriorate as it either fails to recognize pile up or falsely rejects well-separated pulses as being piled-up. In both cases, the resulting energy histograms will be distorted and reported count rates may consistently be too low. 
         [0012]    Advanced digital MCA&#39;s, as described here, can improve the accuracy of gamma-ray spectroscopy in scintillators where the pulse shape changes with temperature or other environmental parameters. 
         [0013]    This supports using scintillator detectors over a wider range of environmental conditions than would ordinarily be possible. 
       SUMMARY 
       [0014]    The need is addressed by a method that uses one or more sensor data, for example from temperature sensors, to adjust the signal-processing parameters to maintain the MCA performance over a wide temperature range. In particular, the method ensures a guaranteed measurement range [E_min, E_max], a near-constant conversion gain, near-optimal energy resolution, and a near-constant pile-up rejection performance. This supports precision radiation measurements with scintillator-based detectors over a wider temperature range than otherwise possible. 
         [0015]    The need can be met by a method to stabilize the performance of a multichannel analyzer for scintillation detectors in varying temperature environments that includes at least: compensating for the temperature induced change of scintillator light pulse shape by always measuring most of the light within a scintillator pulse by varying the integration time, maintaining the pile up rejection performance when the scintillator light pulse shape changes by varying the pileup rejection parameters, and optionally, compensating for the temperature induced gain shifts of the light-to-electronics converter by adjusting the gain of the light-to-electric converter. 
         [0016]    The need can be met by a method for performance stabilization for scintillator-based radiation detectors including at least: receiving ionizing radiation pulses into a scintillator utilizing activated crystals that generate light pulses; converting the generated light pulses into analog electrical pulses; converting the analog electrical pulses into digital pulses; and processing the digital pulses to provide a measure of the energy E received in the scintillator. 
         [0017]    The need can be further met by the method wherein the step of processing the digital electrical pulses to provide a measure of the total pulse energy (E) received in the scintillator is computed from a sum of y_e[n]=(y[n]−y_dc) data where y[n] is the output of the ADC n clock samples after a trigger and y_dc is a direct current baseline in the absence of a pulse. 
         [0018]    The need can also be met by an apparatus for compensating for pulse shape variation of light generated by scintillator detectors including at least: a scintillator for receiving radiation pulses and generating light pulses; a light-to-electric pulse converter to convert the light pulses to analog electrical pulses; an optional amplifier to amplify the analog electrical pulses; an analog to digital converter to convert the analog electrical pulses to digital pulses; a pulse-processing unit to measure the energy of incoming digital pulses; an optional computational unit for managing the pulse processing unit and gathering information from the pulse processing unit. 
         [0019]    In one embodiment this application anticipates using thallium-doped Nal as the scintillator. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0020]      FIG. 1  illustrates an apparatus that could be used to carry out the method of this disclosure. 
           [0021]      FIG. 2  illustrates scintillator pulses at low and high temperature as converted by the ADC and with the DC-offset (y_dc) subtracted. 
           [0022]      FIG. 3 a    illustrates the impact of pulse lengthening on gain and energy resolution when measuring a Cs-137 radioactive source if the integration time IT is not adjusted with temperature. 
           [0023]      FIG. 3 b    shows the remaining effect after the gain shift has been corrected for. 
           [0024]      FIG. 4 a    shows a normal Nal-scintillator pulse. 
           [0025]      FIG. 4 b    shows two pulses occurring too close together—a pile up. 
           [0026]      FIG. 4 c    shows a slowed-down Nal-pulse at low temperature. 
           [0027]      FIG. 4 d    shows the integrated pulses from a time shortly before the trigger to an integration time of 1.2 micro-seconds. 
       
    
    
     DETAILED DESCRIPTION 
       [0028]    This need can be met with a digital multichannel analyzer (MCA) used in conjunction with a scintillator where the scintillator pulse shape changes in response to an environmental parameter such as temperature, the MCA including at least analog to digital circuitry that converts a momentary electronic pulse amplitude from the scintillator into a digital number, a digital processing unit for performing computations on a sequence of such digital numbers, a memory device for storing operations parameters or a means to compute the parameters when needed, wherein the method used by the MCA is described below. The parameters may include the integration time, a hold-off time, pileup-rejection parameters, operation voltage for the light-to-electric converter and a digital gain factor. 
         [0029]    The device is illustrated in  FIG. 1 , beginning with a scintillator  10 . This disclosure anticipates a number of different scintillator crystals that can be used, especially those for which the pulse shape changes with temperature, including at least activated sodium iodide or cesium iodide, or other materials that are effective for detecting gamma rays. Also anticipated are europium-doped alkaline earth halides, such as Srl2(Eu), Bal2(Eu) or mixtures thereof. In addition plastic scintillator containing  2 , 5 -diphenyl oxazole (PPO) as the primary scintillating materials and diphenyl anthracene (DPA) as the wavelength shifter to match the emission light wavelengths to the photomultiplier tube sensitivity. 
         [0030]    Scintillator  10  may receive a radiation pulse  5 . The radiation pulse may be energetic photons &gt;10 eV, atomic or nuclear particles such as electrons, positrons, protons, neutrons, alpha-particles, or others. When the scintillator crystal receives a radiation pulse it creates a light pulse  15  that passes to a light-to-electrical pulse converter  20 . This may be a photomultiplier, SiPMT, an avalanche photodiode, or a photodiode. Device  20  then creates an electrical pulse  25  that substantially has the same amplitude vs time characteristic (pulse shape) as the original light pulse. Electrical pulse  25  passes to a series of subsystems ( 30 ,  40 ,  50 ) that represent a multichannel analyzer (MCA). Electrical pulse  25  first passes to an optional electronic amplifier  30  that amplifies electrical pulse  25  and substantially preserves the electronic pulse shape in creating an amplified electronic pulse  35 . In some embodiments the amplifier  30  may be omitted. Pulse  35  passes to an analog to digital converter (ADC)  40  that converts the signal to a digital signal  45  that substantially preserves the electronic pulse shape. A pulse-processing unit  50  applies methods to recognize incoming pulses and measures their energy taking their pulse shape into account. An optional computer or microcontroller or field programmable gate array  60  issues commands and exchanges data with the pulse processing or signal processing unit  50  to retrieve pulse energy and pulse type information. 
         [0031]    For each ionizing-radiation pulse, the MCA receives an electronic pulse proportional, or nearly proportional, to the momentary light output from the scintillator. The electronic pulse is amplified if needed and presented to the analog to digital circuitry. This creates a digital image of the electronic pulse and therefore of the momentary scintillation light pulse. 
         [0032]    A digital pulse processing unit monitors the output from the analog to digital circuitry. In the absence of a pulse, it measures the direct current (DC) baseline (y_dc). The DC-baseline is always subtracted from the analog to digital circuitry samples (y[n]) to create the y-samples (y_e[n]). When a pulse is recognized by the trigger logic, the digital shaper is activated. The digital shaper includes a method to compute a sum of (y_e[n]) values. On trigger, the digital shaper performs the computation: 
         [0000]        y _ e[n ]=( y[n]−y _ dc ) 
         [0033]    where n is the sample count since the last trigger. A predetermined integration time number IT of y_e samples are summed to provide a measure of the energy (E) deposited in the scintillator. The digital shaper performs many of these computations in a parallel manner. 
         [0034]    The sum over IT samples of y_e serves as a measurement of the original ionizing pulse energy. A judicious choice of the parameter IT as a function of temperature significantly improves the accuracy and energy resolution of the radiation pulse measurement, compared to prior art. To accomplish this, in one embodiment, it is assumed that a temperature sensor (not shown) measures the scintillation detector temperature and a signal processor uses the temperature to determine the various parameter values that stabilize detector performance. In another embodiment an external sensor and processor measures the temperature and determines the temperature-dependent coefficients and writes the updated coefficients to the signal processor of the scintillator. 
         [0035]    Other operating parameters, such as the hold-off time and the pileup rejection controls may also be chosen as a function of temperature. The hold-off time avoids retriggering of the system on the falling edge of a pulse and it will typically be equal or slightly larger than the integration time (IT). 
         [0036]    In one embodiment the pile up rejection may be governed by a single parameter: PUT=the time at which 40% of the light is collected. That parameter would then be selected as a function of temperature. 
         [0037]    To maintain the active range of measurable radiation pulse energies, [E_min, E_max], the method would adjust the gain of the light-to-electric converter by adjusting one or more parameters. Typically, only one parameter, the operating voltage, needs to be adjusted. 
         [0038]    The MCA delivers event data (sums of y_e samples). A plurality of event data can be graphically represented in a category plot, often called a histogram. It depicts the frequency at which the energies (E) of measured events fall into given ranges, called bins. In a typical histogram the bins have a width (dE) and the condition for an event to fall into the n′th bin is expressed as: 
         [0000]        n*dE←E &lt;( n+I )* dE.    
         [0039]    The content of a bin is a number. At the beginning of a histogram data acquisition the content of all bins is set to zero. During data acquisition, whenever an event&#39;s energy falls into a certain bin, that bin&#39;s content is incremented by 1. Over time a histogram builds up. 
         [0040]    The performance of a radiation detection system is judged by the contents of the histogram after some acquisition time. In the case of measuring radiation from the nuclear isotope Cs-137 there will be a complex histogram. 
         [0041]    The highest energy peak, a bell-shaped curve, will correspond to gamma-rays with energies of approximately 662 keV. Its full width at half-maximum height is an indicator of scintillator performance. The smaller the width, the better the performance. 
         [0042]    Another important measure of the system performance is by how much the peak position in the energy histogram will drift as the temperature changes. The requirement is that any peak will drift as little as possible. 
         [0043]    Turning to  FIG. 2 , a graph shown generally as the numeral  200 , shows a Nal(TI) scintillator pulse at room temperature  210  and at −20 deg C.  220 . The vertical stroke at 1.0 micro-seconds may indicate the fixed integration time. The pulse will be summed or integrated from 0 to this limit. Any pulse fraction beyond this limit is ignored. Hence, information is lost for the slow pulse. This creates an apparent gain shift and reduces the accuracy of the energy measurement. By careful calibration work, the gain can be adjusted as a function of temperature to overcome this. 
         [0044]    Turning now to  FIG. 3 a   , shown generally as the numeral  300 , an energy spectrum of Cs-137 from a Na(TI) detector calibrated at room temperature  310 . The full-energy peak at 0.662 keV has an energy resolution of 6.5% fwhm (full width at half maximum). The solid line  320  shows the energy spectrum obtained from the same detector at −20 deg C. while keeping the integration time constant and accounting for the gain drift of the photomultiplier and the electronics. Due to the lengthening of the pulse, only a part of it is being integrated, leading to an apparent loss of gain and a poorer energy resolution (8.5% fwhm for the full-energy peak). 
         [0045]    In  FIG. 3 b   , we make up for the apparent gain-loss of the lengthened pulses by applying a separate gain factor, so the peaks measured at cold  370  and room temperature  360  overlap. This gain factor can be developed by experimentation to be applied for any temperature. It can then be used to adjust the gain of the light-to electric converter according to an experimentally predetermined function or based on a look-up table of gain versus temperature. But still, the two histograms are quite different, because of the poorer energy resolution at −20 deg C. 
         [0046]    There are many algorithms for radio-isotope identification in use by industrial and military operators. All algorithms are carefully tuned and tested for a particular detector type and size. A change of energy resolution is highly undesirable. Hence, the integration (IT) time needs to be adjusted as a function of temperature. 
         [0047]    As previously discussed, the second impact a changing pulse shape has is on the accuracy of the energy measurement. In practice the shaping time or integration time is chosen such as to use about 90% of the light from the scintillation pulse for the measurement of the energy. Using less will deteriorate the precision of the measurement. If the scintillation pulses lengthen, for example due to a temperature change, they may become too long for the given operating parameters of the MCA. As a result, the MCA will measure energy values that are systematically smaller than the correct values. It will also lose accuracy of measurement (energy resolution) as it samples less than 90% of the scintillation light. The solution is to adjust the integration time as a function of temperature so that a longer integration time is used at lower temperatures. This change in integration time can be studied experimentally, so that a required integration time as a function of temperature is known and can be applied automatically as a function of the temperature. 
         [0048]    In  FIGS. 4 a  through 4 d    the concept of pile up rejection is illustrated.  FIG. 4 a   , shown generally by the numeral  400 , is the shape of a single isolated Nal(TI) light pulse at room temperature. 
         [0049]    In  FIG. 4 b    is an illustration for the case of 2 close-by pulses at room temperature. They occur within 0.40 microseconds of each other. The second pulse  440  piles on top of the tail end of the previous pulse  430 . They occur too close together (less than the integration time of 1.2 microsecond) and therefore cannot be measured accurately. Hence, they must be rejected as “piled up”. 
         [0050]    In  FIG. 4 c   , shown generally by the numeral  450 , a Nal(TI) pulse at −20 deg C. is illustrated. It is much slower than a room-temperature pulse  410  of  FIG. 4 a    and does not return to near-baseline values at 1 microsecond as the room temperature pulse does. 
         [0051]      FIG. 4 d   , shown generally as the numeral  470  illustrates an overall algorithm to determine pile up. Plotted are time integrals over the pulses shown in  FIGS. 4 a , 4 b  and 4 c   . For example the time integral of the single room temperature pulse  410  of  FIG. 4 a    is shown as the integral curve  480 . The time integral of the single cold temperature pulse  460  of  FIG. 4 c    is shown as the integral curve  490 . And curve  500  represents the time integral of the two piled up pulses of  FIG. 4 b   . The time integrals have been normalized to unity at 1.2 microseconds. We use the time (T 50 ) at which the integral crosses the 50% threshold as an indicator of pulse pile up. For a room temperature Nal(TI) pulse the time is at 0.28 microsecond, represented by the vertical line at 0.28 microseconds. For the piled-up pulse of  FIG. 4 b    it is 0.43 microsecond. For the slow pulse shown in  FIG. 4 c    it is 0.46 microsecond. 
         [0052]    A simple algorithm may measure T 50  for each recognized pulse and accept pulses for which T 50 &lt;0.35 microsecond and reject pulses for which T 50 &gt;=0.35 microsecond as “piled-up”, but such an algorithm would incorrectly reject the slow pulses as piled up. Hence it is necessary to adjust the integration time to the expected pulse length to maintain a valid pile-up rejection. In addition, one or more pile up parameters may have to be adjusted in response. For example, the pile up inspection time PUT covers a certain fraction of the light pulse energy. That fraction may become temperature dependent. 
         [0053]      FIGS. 4 a , 4 b , 4 c , 4 d    illustrate the concept of pile up recognition. PUT&lt;IT represents the time by which 40% of the light in the pulse should be accumulated for a regular pulse. Pileup may be recognized if less than 40% of the light is collected during the time PUT compared to the integration time IT. 
         [0054]    Although certain embodiments and their advantages have been described herein in detail, it should be understood that various changes, substitutions and alterations could be made without departing from the coverage as defined by the appended claims. Moreover, the potential applications of the disclosed techniques is not intended to be limited to the particular embodiments of the processes, machines, manufactures, means, methods and steps described herein. As a person of ordinary skill in the art will readily appreciate from this disclosure, other processes, machines, manufactures, means, methods, or steps, presently existing or later to be developed that perform substantially the same function or achieve substantially the same result as the corresponding embodiments described herein may be utilized. Accordingly, the appended claims are intended to include within their scope such processes, machines, manufactures, means, methods or steps.