5.3.19.3. Continuous energy distribution with double (energy & angular) differential type¶
5.3.19.3.1. e-type = 41, (51)¶
Any energy distribution can be specified by giving data set of energy bins e(i) and differential probabilities of the particle generation \(\left[\frac{d^2\varphi}{dE/d\Omega}\right](i,k)\) by hand. The integrated number of the particle generation in the bin is proportional to \(\left[\frac{d^2\varphi}{dE/d\Omega}\right](i)*\{e(i+1)-e(i)\}*\{a(k+1)-a(k)\}\) , and the specified energy distribution is statistically described. For 51 case, energy is given by wave length (A).
ne : Number of energy group. If it is given by positive number, source particles are generated so that the energy differential fluxes in the unit of [1/MeV/sr] become constant in each bin. On the other hand, if ne is negative, the fluxes in the unit of [1/Lethargy/sr] become constant in each bin. Data must be given from the next line by the format \((e(i), (\left[\frac{d^2\varphi}{dE/d\Omega}\right](i,k),k=1,|na|), i=1,\lvert ne \rvert),\; e(\lvert ne \rvert + 1)\) . The integrated number of the particle generation in each energy bin is proportional to \(\left[\frac{d^2\varphi}{dE/d\Omega}\right](i)*\{e(i+1)-e(i)\}*\{a(k+1)-a(k)\}\) .
5.3.19.3.2. e-type = 42, (52)¶
The same energy distribution as in the case of e-type = 41, (51) can be specified. Unlike e-type = 41, (51), the distribution is described by giving data set of energy bins e(i) and weights of the source particle w(i) by hand. The number of source particles 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 \(w(i,k)*\{e(i+1)-e(i)\}*\{a(k+1)-a(k)\}\) . The number of source particles generated in each bin can also be changed by specifying p(i). For 52 case, energy is given by wave length (A).
ne : Number of energy group. If it is given by positive number, source particles are generated so that the energy differential fluxes in the unit of [1/MeV/sr] become constant in each bin. On the other hand, if ne is negative, the fluxes in the unit of [1/Lethargy/sr] become constant in each bin. Data must be given from the next line by the format \((e(i),( w(i,k),k=1,|na|),i=1,|ne|), \lvert ne \rvert),\; e(\lvert ne \rvert + 1)\) . In default (p-type=0), equal number of particle is generated in each cell. The integrated number of source particles generated in each bin is proportional to p(i).
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,ne).It is not possible to specify different values for each angular bin.