5.3.19.7. Energy distribution defined by analytical functions¶
5.3.19.7.1. e-type = 5, (15)¶
Differential spectrum dφ/dE(i) is given by f(x). For 15 case, energy is given by wave length (A).
f(x) : Any analytical function of x, Fortran style. x denotes energy [MeV/n]. One can use intrinsic functions and constants C, e.g., f(x) = exp(-c1*x**2).
nm : number of energy group. If it is given by positive number, linear interpolation is assumed in a bin. If negative, logarithmic interpolation is assumed in a bin. Integrated number of source particles generated in each cell is proportional to \(f(x)\).
eg1 : minimum cut off for energy distribution [MeV/n].
eg2 : maximum cut off for energy distribution [MeV/n].
5.3.19.7.2. e-type = 6, (16)¶
The same energy distribution as in the case of e-type = 5, (15) can be specified. Unlike e-type = 5, (15), the number of source particle generated in each bin is the same for all energy bin, but integrated values of the weight of source particles are adjusted to be proportional to \(f(x)\). The number of source particles generated in each bin can also be changed by specifying p(i). For 16 case, energy is given by wave length (A).
f(x) : Any analytical function of x, Fortran style. x denotes energy (MeV/n). One can use intrinsic functions and constants C, e.g., f(x) = exp(-c1*x**2).
nm : number of energy group. If it is given by positive number, linear interpolation is assumed in a bin. If negative, logarithmic interpolation is assumed in a bin. In the default (p-type=0), equal numbers of particles are generated in each cell. The integrated number of source particles generated in each bin is proportional to p(i).
eg1 : minimum cut off for energy distribution [MeV/n].
eg2 : maximum cut off for energy distribution [MeV/n].
p-type = 0, 1 : (D=0) generation option.
For 0, p(i)=1 for all i is assumed without the following data.
For 1, p(i) must be given from the next line by the format (p(i),i=1,nm).