SOURCE

Summary

Invokes the use of a user-defined source routine SOURCE to sample the primary particles. ​​

​

See also BEAM, BEAMPOS, POLARIZAti, SPECSOUR, USRICALL, USRGCALL
This option allows to input up to 18 double precision parameters and one character string chosen by the user, to be passed to the user routine SOURCE. To pass more than 6 parameters, two successive SOURCE cards are required.

First card:

WHAT(1)

…WHAT(6): user parameters to be passed to the SOURCE routine as a double precision subarray WHASOU(1)…WHASOU(6) (via COMMON SOURCM).

SDUM

= any user dependent character string (not containing “ & ”), to be passed to the SOURCE routine as a character variable SDUSOU (via COMMON CHEPSR).

Continuation card (if present):​

WHAT(1)

…WHAT(6): user parameters to be passed to the SOURCE routine as a double precision subarray WHASOU(7)…WHASOU(12).

SDUM

= “ & ” in any position in column 71 to 78 (or in the last field if free format is used)

Default

(option SOURCE not given): subroutine SOURCE is not called

Notes

  1. In many simple cases, the primary particle​ properties can be defined by just two input cards: BEAM and BEAMPOS. The two options define the type of the particle, its energy/momentum (monoenergetic or simply distributed) and its starting position and direction (also sampled from a few simple distributions centred around the z-axis). A third option, POLARIZAti, can be used to complete the description of the primaries in case they are totally or partially polarised​.

    However, there are more complex situations where the type of primary particles and their phase space​ coordinates must be sampled from different types of distributions, read from files, or generated by a rule decided by the user. To handle these cases, it is possible to call a user routine SOURCE which can override totally or in part the definitions given in input. The call is activated by option SOURCE. A default version of the routine, which leaves any other input definition unchanged, is present in the FLUKA library. Instructions on how to write, compile and link SOURCE are given in SOURCE: user-written source.

  2. Even when overridden by SOURCE, the momentum or kinetic energy defined by WHAT(1) of option BEAM is meaningful, since it is taken as maximum energy for several scoring facilities and for cross section tabulations. Therefore, it is recommended to input in any case a BEAM card with the maximum energy expected in the current run.

  3. The user has the possibility to write a flexible SOURCE routine which can handle different cases without being recompiled and linked. The 12 WHASOU optional double precision parameters and the SDUSOU character string can be combined to provide a multitude of possible options.

  4. Initialisations (for instance reading data from files, or spectrum normalisation needed for sampling) can be done in SOURCE itself, as explained in SOURCE: user-written source, or in two other user routines USRINI​ and USRGLO​. USRINI is called every time a card USRICALL is found in input (see USRINI: USeR INItialisation), while USRGLO is called before any other initialisation made by FLUKA. Note that more than one primary particle can be sampled in a same call to SOURCE and loaded into stack for later transport. A further user routine, USREIN​ (see USREIN: USeR Event INitialisation (called before each event)), which is called just after SOURCE and before the sampled primary (or primaries) are transported, allows the user to do an initialisation at the beginning of each event. An event​ is defined as all the histories of primaries sampled in a single call to SOURCE and of their descendants. Routine USREOU​ is called instead at the end of each event (see  USREOU: USeR Event OUtput (called at the end of each event)).

  5. In old versions of FLUKA, the call to SOURCE was requested by means of a flag in card START. This feature has been discontinued.

Example 1:

* A user-written SOURCE routine is called without passing any parameter.
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
SOURCE

Example 2:

* Here the user passes to the SOURCE routine 7 numerical values and one
* character string. These can be used as free parameters inside the routine,
* or as flags/switches to select between different options
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
SOURCE         12.58       1.0      -14.  0.987651      100.      365.FLAG18
SOURCE         999.2                                                  &