sesh
calculates corrections for time-of-flight experimental data.
The code sesh
presented here is a modernized version of SESH re-written in Modern Fortran.
Original code is by F.H. Froehner published in the report: "SESH - A Fortran IV Code For Calculating
the Self-shielding and Multiple Scattering Effects For Neutron Cross Section Data Interpretation in
the Unresolved Resonance Region", Report GA-8380 (1968). This code has been adapted over the years
from the input of several nuclear data researchers and their contributions are appreciated.
sesh
solves the problem of resonance self-shielding, the experimental effect seen in time-of-flight
experiments where the measured time/energy-averaged transmission data cannot be compared directly to
theoretical cross section due to the inequality:
for non-linear equations.
Requirements:
- Fortran compiler:
gfortran
(gcc 10.2)
It should work with other compilers, but has only recently been tested with gfortran
.
To compile:
Run the following in a bash-like terminal starting in the sesh main project directory:
mkdir build
cd build
../compile-sesh.sh ../
make
This will put the sesh
executable in your sesh/build/
directory.
To run executable:
Open bash-like terminal and run:
sesh
which will prompt you for input and output file names. If you have already defined these names
in a dedicated file like auto-in
, you can simply run sesh < auto-in
. To run the example file
provided, from the sesh/build/
directory run:
sesh < ../auto-in
An input file called ta181.inp
is given as an example for SESH input. The SESH input is more
fully defined in the report by F. Froehner described here.
Input file notes to augment Frohner's report:
- After the energy groups and number of histories card you must include a line with a space in it
- An update was made to include a description of beam radius on sample (at/b), this follows geometry card 3.3