Help for FERGIE

FERGIE help

(VERSION = 4.24 - 21/06/00)

A program to Facilitate calculation of Energy Resolution Governing Inverse-kinematics Experiments. The program is written in C++ and is currently available on the Sun phx1np/jsw/fergie at the University of Surrey and on WindowsNT at LNS Catania. A description of the physics background may be found in Nucl. Instr. Methods A396 (1997) 147.

Input description

Card 1: TITLE (A60)

Any alphanumeric title. This is printed at the top of the output. Card 2: KTRL (15I1) An array of 15 control integers. The usage is:

KTRL(1) - controls how many columns of output you get.

If 0 (default) you get energy straggling, multiple scattering and energy-loss difference, i.e. the lot, plus their linear and quadratic sums. This won't conveniently fit on an 80-column page (screen or paper). Print in LANDSCAPE orientation.

If 1, you get energy straggling and multiple scattering (only).

If 2, you get multiple scattering and energy-loss difference (only).

If 3, you get energy straggling and energy-loss diffrence (only).

Note: if you get warning messages that a hidden target contribution was bigger than those output, you should try KTRL(1) = 0 or else another new value to check you are not mislead by the numbers currently output.

KTRL(2) - selects multiple scattering formula If 0 (default) use D.K. Scott's formula for multiple scattering.

If 1, use formulae from Marion & Young.

If 2, use formulae from Review of Particle Properties (1994).

KTRL(3) - selects kinematics If 0 use normal first kinematic solution.

If 1 use second kinematic soln. Th(c.m.) is output as 180 - th(c.m.)

If 2, first kinematic soln. with th(c.m.) output as 180 - th(c.m.)

If 3, second kinematic soln., th(c.m.) output as is.

KTRL(4) - controls optional additional output If > 0, gives additional output to a file 'auxilout', including kinematic factors and multiple scattering half angles.

If > 1, gives debug output to stderr. The higher the value of KTRL(4), the more output generated.

KTRL(5) - target options If > 0, read Target Angle (degrees). 0 degrees is perpendicular to the nominal zero degrees of the beam). If the detected particle angle is 90 degrees to the target angle, then the target will appear to be infinitely thick and the energy loss and multiple scattering routines will fail.

If > 1, also read and calculate contribution from a given target non-uniformity.

KTRL(6) - mass excess options If 0, nuclear masses are assumed to be entered exactly in AMU and no additional input is read.

If 1, read mass excesses in keV from input file following Card 7.

If 2, read mass excesses from mass table.

In cases 1 and 2, the baseline nuclear masses are read first (from Card 7). These baseline masses in AMU should include any nuclear excitation energy.

This addition of excitation energy in AMU also has to be done for KTRL(6)=2: then Card 7 should be the integer mass + any excitation energy in MeV/AMU). For KTRL(6)=1, any excitation energy can be more easily added in keV to the mass excess on Card 7a.

KTRL(7) - Detector energy resolution If 1, use a fixed EB resolution (entered in keV), rather than the default inverse fractional energy resolution. Entered as a FWHM. KTRL(8) - Beam energy resolution If 1, read Percent Beam Energy Resolution, to be converted to a spread in excitation energy and added in quadrature with the "energy resolution" contribution. As always, a FWHM. KTRL(9) - Doppler broadening If > 0, read Excitation Energy of detected particle, and calculate the Doppler broadening from the gamma emission in flight and added in quadrature with the "energy resolution" contribution.

If 1, this is done assuming "worst case" distribution of 2 * beta * Egamma.

If 2, there is an empirical reduction by 60%.

See Bohlen et al., Z. Phys. A285 (1978) 379, for proper treatment.

KTRL(10) - Selection of energy straggling formula If 0, use Energy Straggling formulae from D.K. Scott's lecture notes.
Otherwise use formulae from Tschalar a l'INTENSITY. Card 3: TGT_THICK (3F) The target thickness in (mg/cm2), followed [if KTRL(5) or KTRL(11) are set] by the target angle, in degrees, perpendicular to the beam. If KTRL(11) is set, also read per cent target non-uniformity. Card 4: TGT_MATL (A5) Target material as a 5-character string. Any isotopic element such as "12C" or "28SI" is valid. Otherwise, a selection of compounds, but not mixtures, known to the STOPRANGE program can be used, such as

"MYLAR", "KAPTO", "SCINT", "POLYP", "METHA", etc.

Card 5: RESLN_ETH (3F) Detector energy resln, overall angular resln, and, if KTRL(8) is set, the percent beam energy spread (%dW/W). The detector energy resolution is in inverse fractional form if KTRL(7) is not set (e.g. an input of "2000" will give a resolution of 1 in 2000). If KTRL(7) is set, the energy resolution is taken as absolute, and entered in keV. All these to be entered as FWHM. Card 6: THDET (3F) The intitial, final and step size for the detected particle lab angle (in degrees). This governs the loop over which the kinematics and energy resolution calculations are run. If there is no kinematic solution for a given angle in the requested range, a message to that effect is printed and no further calculations are performed for that angle. Card 7: A (4F) Mass of the incident, target, detected and residual particles in the reaction. If KTRL(6) = 0, these should be input exactly in atomic mass units (mass of 12C = 12.0000).

If KTRL(6) = 1, they should be input as integers, and the following additional Card 7a will be read.

If KTRL(6) = 2, the masses should be entered as integers, and the mass excesses will be read from a mass table.

Card 7a (Only read if KTRL(6) = 1): AEX (4F) Mass excesses in keV, to be added to the masses read in Card 7. Card 8: Z (4I) Atomic numbers of the incident, target, detected and residual particles. Card 9: QK (3F) The Q-value, the incident beam energy, and (if KTRL(9) > 0) the excitation energy of the detected particle, all in MeV.

Gallery of Fergie Portraits