3.13. Card 10: PHITS and ANGEL output controls

Table 3.13.1 PHITSOUT PHITS: other means

Value	Description
(D=0, DP=1)	Toggles writing of PHITS source cards to the *.pht output file(s). [1] Note: This is named IPHTOUT in PHITS.
= 0 or 1	No output for [source] section based on RI source
= 10 or 11	RI source with proj = all
= 20 or 21	RI source with proj = photon
= 30 or 31	RI source with proj = electron
= 40 or 41	RI source with proj = positron
= 50 or 51	RI source with proj = alpha

Table 3.13.2 RADLEV PHITS: same name

Value	Description
(D=0.0, requires abs(PHITSOUT)\(>\)2)	Minimum activity of nuclides included in RI-source output; this is an additional activity filter to ACMIN. [2]
\(>\) 0.0	Threshold radioactivity in [Bq].
\(<\) 0.0	-1 \(\times\) Threshold radioactivity in [Bq/cm\(^3\)].

Table 3.13.3 IANGBPWR PHITS: same name

Value	Description
(D=0, DP=1)	Toggles whether a visualization of the beam power schedule is included in the plots in the *.ang output file.
= 0	Not included
= 1	Included (does nothing if first output time ITOUT is greater than the end of the final time step ITSTEP where beam power is \(>\) 0)

Table 3.13.4 ANGELOUT PHITS: not settable

Value	Description
(D=0, DP=1)	Toggles writing data for all plots to the *.ang output file.
= 0	Some or no plots will be written.
= 1	(and if all ANGELOUT_* variables = 0) All plots will be written

The following variables beginning with * are actually entered into the input deck as ANGELOUT_ followed by that variable name (without the asterisk (*)).

Table 3.13.5 *REGION PHITS: not settable

Value	Description
	Selection of which regions will have data written to the *.ang output file.
= \(n\)	Listed region numbers (IREGS) or names (HRGCMM) are included.
= ALL	All individual regions will be written/plotted.
= SUM	The sum of the data for all regions will be written/plotted.

Table 3.13.6 *NUCLIDES PHITS: not settable

Value	Description
(D=0, maximum=13)	Number of nuclides to have data, summed over all regions (consistent with the behavior of ANGELOUT_REGION= SUM), written to the *.ang output file, followed by the list of nuclides.
= 0	No individual nuclides are plotted (default behavior)
= \(n\)	\(n\) nuclides will have their information written the *.ang file. Immediately after this line, the next \(n\) lines must each contain a nuclide name formatted as TGNZAS or HTGNZAS shown in Card 9b.

The below variables can be set to 0 (don’t draw) or 1 (do draw) for each *.ang plot.

Table 3.13.7 *ACT PHITS: not settable

Value	Description
	Toggles writing activity plot data to the *.ang output file.
= 0	Not written
= 1	Written

Table 3.13.8 *HEAT PHITS: not settable

Value	Description
	Toggles writing total decay heat plot data to the *.ang output file.
= 0	Not written
= 1	Written

Table 3.13.9 *HEAT_ALPHA PHITS: not settable

Value	Description
	Toggles writing \(\alpha\) decay heat plot data to the *.ang output file.
= 0	Not written
= 1	Written

Table 3.13.10 *HEAT_BETA PHITS: not settable

Value	Description
	Toggles writing \(\beta\) decay heat plot data to the *.ang output file.
= 0	Not written
= 1	Written

Table 3.13.11 *HEAT_GAMMA PHITS: not settable

Value	Description
	Toggles writing \(\gamma\) decay heat plot data to the *.ang output file.
= 0	Not written
= 1	Written

Table 3.13.12 *DOSE PHITS: not settable

Value	Description
	Toggles writing dose rate plot data to the *.ang output file.
= 0	Not written
= 1	Written

[1]

While the RI source functionality provides all emissions at exact energies (necessary for some calculations like spectroscopy and detector response), it is still more computationally efficient to use the pure \(\gamma\)-ray source from PHITSOUT = \(\pm1\) if the resolution of DCHAIN’s available photon spectra is sufficient and the photon emissions are the only significant contributor to the secondary calculation’s desired tally, as is usually the case for simple room dose assessments. As electrons, positrons, and alphas are slower to transport and are more likely to be stopped by geometry materials before reaching a tally region of interest, simulations utilizing them will generally require more histories to achieve comparable statistics.

[2]

Depending on the complexity of your problem, it may be necessary or desirable to increase either RADLEV or ACMIN to eliminate writing of excess nuclides whose emission probability would likely be too low to be simulated by PHITS. When using the RI-source functionality in PHITS, it must read and store all emission energies and probabilities for each nuclide, and having a very excessive number of nuclides present can sometimes cause this initialization processes to exceed its limits and terminate.