• Halflife fitting: halflife.c (Dr. Carl Wheldon)

    Fits a convolution of a Gaussian and an exponential decay (derived here gauss_exp_conv.pdf) to properly take account of finite time resolution in half-life measurements. This program uses non-linear curve fitting via a modified version of CURFIT from Bevington, "Data Reduction and Error Analysis for the Physical Sciences" (p237) as extracted from RadWare developed by David Radford at ORNL. Constant background is taken into account through χ2 minimisation. Input is from the command line. A range of input formats, including 1, 2 and 3 column ASCII data and Ortec .Spe files can be processed. Data compression and ASCII output are implemented. Example fits can be found here.

  • Kinematics and Q-values: ckin.c (Dr. Carl Wheldon)

    Calculates two-body kinematics and Q-values, including decay Q-values. Input is from the command line. This program uses the 2012 mass-table file from AMDC Atomic Mass Data Center. There is now the option (see top of source code) to compile with the Cern Libraries or libwclbes.a (compiled with gfortran, linux, x86_64) library file to enable the calculation of barrier penetration probabilities. To compile libwclbes.a for your system use the tarball wclbes.tar.gz (readme file with instructions for gfortran and f90/f77 inside). The program accepts reactions entered with 2 different formats and an arbitrary number of ejectiles, as well as accepting common notation, such as d (2H), a (4He) etc. Ckin outputs both relativistic and non-relativistic solutions. For two-body kinematics and other related derivations see kinematics.pdf.