RSIC COMPUTER CODE PSR-092
1. NAME AND TITLE
FORIST: Neutron Spectrum Unfolding Code System -- Iterative Smoothing Technique.
FORIST is a modification of PSR-17/COOLC. The smoothing function used in COOLC has been
redefined and an iterative procedure for selecting the widths of the smoothing functions has been
added. This iterative procedure optimizes the energy resolution for a given statistical error of the
Nuclear Engineering Program, University of Illinois at Urbana-Champaign, Urbana, Illinois.
McMaster University, Hamilton, Ontario, Canada.
3. CODING LANGUAGE AND COMPUTER
Fortran IV; IBM 360/370 (A), CYBER (B).
4. NATURE OF PROBLEM SOLVED
FORIST solves the matrix unfolding problem for neutron spectrum measurements using NE-213
5. METHOD OF SOLUTION
A constrained least-squares unfolding procedure is used. The response function used in unfolding
includes the resolution of the spectrometer; therefore, the unfolded spectrum has large errors initially.
The unfolded spectrum must be smoothed to reduce these errors. The smoothing functions are
determined iteratively to optimize energy resolution.
6. RESTRICTIONS OR LIMITATIONS
FORIST is currently dimensioned for a 113 by 77 response matrix. It is written for measurements
consisting of two different gain runs. The iterative procedure used in smoothing is applicable for
unfolding neutron spectra; a slight change in the smoothing procedure may be necessary to apply this
iterative procedure to unfolding gamma-ray spectra.
7. TYPICAL RUNNING TIME
The sample problem, which contains two spectra, takes approximately 30 seconds on an IBM
8. COMPUTER HARDWARE REQUIREMENTS
FORIST is operable on the IBM 360/370 or the CYBER computers. Approximately 270 K bytes
of storage are needed. A scratch disc or tape, printer, and punch are required.
9. COMPUTER SOFTWARE REQUIREMENTS
A Fortran IV compiler is required.
R. H. Johnson, D. T. Ingersoll, B. W. Wehring, and J. J. Dorning, "NE-213, Neutron Spectrometry System for Measurements from 1.0 to 20 MeV," Nucl. Instr. and Meth. 145 (1977) 337-346.
R. H. Johnson, "A User's Manual for COOLC and FORIST," PNE-75-107 (December 1975).
R. H. Johnson, J. J. Dorning, and B. W. Wehring, "Integral Test of Cross Sections Using Neutron Leakage Spectra from Spheres of Iron, Niobium, Beryllium, and Polyethylene," Proceedings of the Conference on Nuclear Cross Sections and Technology, Washington, D.C. (March 1975). [To be published by National Bureau of Standards.]
D. T. Ingersoll and B. W. Wehring, "Gamma-Ray Pulse-Height Response of an NE-213 Scintillation Detector," Nucl. Instr. and Meth. 147 (1977) 551-561.
D. T. Ingersoll, Informal Note on FORIST (July 12, 1976).
11. CONTENTS OF CODE PACKAGE
Included are the referenced documents and one (1.2MB) DOS diskette which contains the source
code and sample problem input and output.
12. DATE OF ABSTRACT
March 1984; reviewed May 1984, updated August 1991.
KEYWORDS: OPTIMIZATION; UNFOLDING; SCINTILLATION DETECTOR