Computer Programs
NEA-0914 REFIT-2009.
last modified: 28-JUN-2022 | catalog | categories | new | search |

NEA-0914 REFIT-2009.

REFIT-2009, Multilevel Resonance Parameter Least Square Fit of N Transmission, Capture, Fission & Self Indication Data

top ]
1. NAME OR DESIGNATION OF PROGRAM

REFIT-2009.

top ]
2. COMPUTERS

To submit a request, click below on the link of the version you wish to order. Rules for end-users are available here.

Program name Package id Status Status date
REFIT-2009 NEA-0914/09 Tested 28-JUN-2022

Machines used:

Package ID Orig. computer Test computer
NEA-0914/09 PC Windows Gitlab,Linux-based PC,PC Windows
top ]
3. DESCRIPTION OF PROGRAM OR FUNCTION

REFIT carries out fits to measured neutron cross-section data, by adjusting the nuclear parameters used in the multi-level R-matrix formalism as well as experimental parameters. The adjustment of the parameters is continued until the calculated transmission and/or reaction yields from neutron Time Of Flight (TOF) measurements agree with the observed data within the limits of the measured uncertainties. In its present form it can perform simultaneous fits on up to 30 sets of measured data. The data sets can be of different types, refer to different target sample thickness and refer to different sample isotopic composition. The types of measurement include transmission and the reaction types: capture, fission, scattering and self-indication. The output gives details of the analysis and all the fitted parameters with their uncertainties. The code also outputs a new list of resonance parameters in ENDF6 format suitable for putting in an evaluation file.

top ]
4. METHODS

The nuclear cross section for each nucleus is calculated from initial parameters using an R-matrix multi-level formalism. The cross section is then Doppler broadened using either an ideal gas model or Moxon's simple Phonon Model, or Meister's Einstein Crystal Model, or Naberejnev's Harmonic Crystal Model. Some of the parameters associated with the Doppler broadening can be adjusted in the fit e.g. in the gas model the effective temperature. The transmission or reaction yield for a given sample is then calculated from the Doppler broadened cross sections, given abundances of the nuclei and the thickness of the sample in the neutron beam. In the case of reaction data the calculated reaction yield includes the effects of a multiple collision of the neutrons in the sample and the detection of scattered neutrons. The theoretical curves are then folded with the resolution function which can be calculated from parameters that include: the initial pulse width, the decay of neutrons leaving the source and moderator, etc or from tabulated values. The nuclear parameters are the resonance energies, neutron widths, radiation widths and effective nuclear interaction radii. Provision has also been made to fit experimental parameters such as zero delay, flight path length, background, normalisation, resolution parameters etc. as well as the nuclear parameters.

top ]
5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

The program in its present form is set up for a maximum of 30 different samples, 40 different isotopes with a total of only 8000 resonance parameters, 185 variables and a total of 10000 data points.

top ]
6. TYPICAL RUNNING TIME

The running time is variable and case dependent. The test cases run with the following elapsed time on the PC:
Copper test2 takes a few seconds.
Resium theoretical tests 1 and 2 take a few seconds and stop when measurement data are required.
U238 test takes 2.5 hours.

top ]
8. RELATED OR AUXILIARY PROGRAMS

A modified version of the Harwell computer library minimising routine VA02A is used to obtain the best fit to the data by varying the required parameters, and outputs the parameters at the minimum value of chi-squared. The modified library routine SV01AS outputs the uncertainties (one standard deviation) together with the correlation coefficients matrix. The source of these routines is present with the main source code.

 

Also contained in the package is the following auxiliary program:
RESCON-2009 converts ENDF format into REFIT format.

top ]
9. STATUS
Package ID Status date Status
NEA-0914/09 28-JUN-2022 Tested at NEADB
top ]
10. REFERENCES
NEA-0914/09, included references:
- M.C. Moxon, T.C. Ware, C.J. Dean:
REFIT-2009: A Least-Square Program for Resonance Analysis of Neutron
Transmission, Capture, Fission and Scattering Data
Users' Guide for REFIT-2009-2009-10 (April 2010)
top ]
11. HARDWARE REQUIREMENTS

WINDOWS (done in 2011):

  • COMPUTER : Dell INTEL Duo E6550,  2.33 GHz (for Windows)

  • OPERATING SYSTEM : DOS, under Windows XP Pro SP 3

  • COMPILER: lahey/ Fujitsu F 95  v7.1  (Windows)

 

LINUX (done in 2022):

  • MACHINE: CPS GitLab Cloud Service  

  • OPERATING SYSTEM: Dockerized Ubuntu 16.04, Ubuntu 18.04, and Ubuntu 20.04 respectively.

  • COMPILER: gfortran v5, gfortran v7, and gfortran v9 for each respective Ubuntu version tested.

 

OTHER: GitLab repository located at http://git.oecd-nea.org/databank/cps/refit

 

top ]
12. PROGRAMMING LANGUAGE(S) USED
Package ID Computer language
NEA-0914/09 FORTRAN-90, FORTRAN-95
top ]
15. NAME AND ESTABLISHMENT OF AUTHORS

REFIT-2007 (NEA-0914/07) and REFIT-2009 (NEA-0914/08/09) were developed by:

 

M. C. Moxon (contractor), T. C. Ware (contractor) and C. J. Dean

Serco Technical and Assurance Services

Building A32, Winfrith Technology Centre, Winfrith,

Dorchester, Dorset DT2 8DH

United Kingdom

 

REFIT was originally developed by Mick Moxon of the Nuclear Physics Division of the United Kingdom Atomic Energy Authority (UKAEA). The division was transferred to AEA Technology (AEAT) in 1990 before AEAT separated from UKAEA in 1996. The division was part of AEAT sold to Serco in 2001, who then sold it to AMEC in 2012, then to Jacobs in 2020. Mick Moxon retired from AEAT in 1993 and developed REFIT-2007 and REFIT-2009 under contract to Serco (2006-2010).

 

The authors can now be contacted as follows:

 

T. C. Ware

The ANSWERS Software Service

Jacobs

Kings Point House, Queen Mother Square, Poundbury,

Dorchester, Dorset DT1 3BW

United Kingdom

E-mail: tim.ware(at)jacobs.com

 

Sadly, Christopher Dean passed away in 2012 and Mick Moxon passed away in 2020.

top ]
16. MATERIAL AVAILABLE
NEA-0914/09
REFIT_V2009.inf     Information file
REFIT-2009-10.exe   REFIT executable (provided by the authors)
RESCON-2009-10.exe  RESCON executable (provided by the authors)
REFIT_NEA_DB.exe    REFIT executable as obtained at the NEA Data Bank in 2011
\data\              Resonance and measurement data for test cases
\doc\REFIT-2009_Manual.pdf Users' Guide for REFIT-2009-2009-10 (April 2010)
\rescon\RESCON-2009-10.FOR RESCON source (auxiliary code)
\source\            Source code directory
DOPINP-2009-10.FOR  Part of REFIT source used to model Doppler Broadening
                    at low neutron energies
REFIT-2009-10.FOR   REFIT source
*.REF               Common files used by REFIT (74 files)
\test\              test cases
*.cmm               Command files
*.cnt               Control files
*.out               Formatted output files
*.oux               Summary output files
*.plt               Formatted output file suitable for plotting with
                    common data graphing software packages
*.fit0; *.fitx      Unformatted output files
top ]
17. CATEGORIES
  • M. Data Management
  • O. Experimental Data Processing

Keywords: Doppler broadening, R matrix, capture, least square fit, multilevel analysis, neutron cross sections, neutron reactions, neutron spectra, resolved region, resonance integrals, transmission.