7.14. [ T-Interact ] section (formerly named [ T-Star ])

This tally gives the number of interactions occurring in the specified regions. It was formerly named [t-star] because it can calculate the star density, but we changed the name from version 3.04 in order to explicitly state that it is the tally for counting the number of interactions. Note that the former name, [t-star], is still accepted even after version 3.04. To calculate the number of atomic interactions using this tally, use the EGS5 mode, negs=1,2. To calculate the number of nuclear interactions using nuclear data library, use the Event Generator mode, e-mode>=1.

Table 7.14.1 mesh

value	explanation
reg, r-z, xyz	Mesh type. A mesh type subsection is required below this option.

Table 7.14.2 MorP

value	explanation
mean, prob (omissible)	Output the mean number (mean) or probability density (prob) of interactions. Only when mesh=reg, select prob.

Table 7.14.3 part

value	explanation
all (default), particle name	Tally particle.

Table 7.14.4 material

value	explanation
(optional)	Specify materials for scoring.
all	all is the default, same as no definition.
number of materials	To set number of materials, define the material numbers in the next line. The number of materials can be set as a negative value, in which case the specified materials are not included for scoring.
(next line)	2 5 8 Material numbers. If the list does not fit on one line, it can be split across multiple lines without using the line-connection symbol.

Table 7.14.5 mother

value	explanation
(optional)	Specify mother nuclei.
all	all is the default, same as no definition.
number of mother nuclei	When setting the number of mother nuclei, define the mothers in the next line. The number of mother can be set as a negative value, in which case the specified mothers are not included for scoring. When output=atomic, mother cannot be specified.
(next line)	208Pb Pb Nucleus if specified using mass. Without mass, all nuclei are isotopes of Pb. To specify multiple mother groups, use multiple [t-interact] tallies.

Table 7.14.6 e-type

value	explanation
1, 2, 3, 4, 5	Energy mesh. See Section 6.6.1 for mesh subsection format.

Table 7.14.7 t-type

value	explanation
1, 2, 3, 4, 5 (optional)	Time mesh. See Section 6.6.1 for mesh subsection format.

Table 7.14.8 unit

value	explanation
1, 2	1: [1/cm^3/source] 2: [1/cm^3/MeV/source]
11, 12	11: [1/cm^3/nsec/source] 12: [1/cm^3/nsec/MeV/source]

Table 7.14.9 axis

value	explanation
eng, reg, x, y, z, r	x axis value of output data.
xy, yz, xz, rz	2-dimensional.
t	Time axis.
act	Number of interaction, for MorP=prob.

Table 7.14.10 maxact

value	explanation
(optional, D=100)	Maximum number of interactions per history for axis=act. If the number of interactions exceeds this value, the history is counted as ‘overflow,’ whose probability is output just above the ‘sum over’ line.

Table 7.14.11 samepage

value	explanation
(optional, D=part)	The type of data to be displayed on the same page of the image output file. Parameters that can be defined in axis can be specified.

Table 7.14.12 file

value	explanation
file name	Define output file names. This is required by each setting of axis.

Table 7.14.13 resfile

value	explanation
(optional, D=file)	Define a file name of the past tally in the restart calculation. Even if several axis parameters were defined, specify only one resfile.

Table 7.14.14 factor

value	explanation
(optional, D=1.0)	Normalization factor.

Table 7.14.15 title

value	explanation
(optional)	Title.

Table 7.14.16 angel

value	explanation
(optional)	ANGEL parameters.

Table 7.14.17 sangel

value	explanation
(optional)	Special format for ANGEL parameters.

Table 7.14.18 2d-type

value	explanation
1, 2, 3, 4, 5, 6, 7	Options for 2-dimensional plot.
(optional, D=3)

Table 7.14.19 x-txt

value	explanation
(optional)	\(x\) axis title.

Table 7.14.20 y-txt

value	explanation
(optional)	\(y\) axis title.

Table 7.14.21 z-txt

value	explanation
(optional)	\(z\) axis title.

Table 7.14.22 output

value (corresponding reaction ID number in parentheses)	explanation
	Reaction type. When specified by a string keyword, only one can be set. When specified by reaction ID numbers, multiple numbers separated by spaces can be set (see below).
all	All reactions.
decay (1000)	Decay reaction.
elastic (2002)	Elastic scattering.
nuclear	Nuclear reaction excluding elastic scattering. This corresponds to nonela in [t-product].
fission (2018)	Fission.
absorption	Absorption.
heavyion	Heavy ion reaction.
transmut	Reactions inducing transmutation of target nuclei.
atomic	Atomic interactions excluding electron multiple scattering and Rayleigh scattering.
deltaray (3001)	Delta-ray production.
atmflu (3002)	Atomic fluorescence x-ray emission.
auger (3003)	Auger electron emission.
brems (3004)	Bremsstrahlung.
photoelec (3005)	Photoelectric effect.
compton (3006)	Compton scattering.
pairprod (3007)	Electron-positron pair production.
annih (3008)	Positron annihilation.
knockelec	Production of knock-on electrons by electrons and positrons.
ts_hit (4000)	All events in track-structure mode.
ts_ioniz (4002)	Ionization in track-structure mode.
ts_bio (4002 4003 4004 4009)	Sum of the events that can cause DNA damage, ionization, electronic-excitation, electron capture, and dissociative electron attachment, in track-structure mode.
ts_w (4002 4010)	Sum of the events that generate electrons, ionization and electron stripping, in track-structure mode. The w-value can be deduced from the quotient of deposition energy and this value.
ets_elas	Elastic scattering in PHITS-etsmode [1].
ets_ioniz	Ionization in PHITS-etsmode.
ets_e-exc	Electronic excitation in PHITS-etsmode.
ets_dea	Dissociative electron attachment in PHITS-etsmode.
ets_v-exc	Vibration excitation in PHITS-etsmode.
ets_p-exc	Photon excitation in PHITS-etsmode.
ets_r-exc	Rotation excitation in PHITS-etsmode.
ets_plasmon	Plasmon excitation in PHITS-etsmode.
ets_ioniz_e-exc	Sum of the ionization and electronic excitation in track-structure mode.
ets_hit	All reactions in PHITS-etsmode.
kurbuc_elas	Elastic scattering in PHITS-KURBUC [2].
kurbuc_ioniz	Ionization in PHITS-KURBUC.
kurbuc_e-exc	Electronic excitation in PHITS-KURBUC.
kurbuc_e-cap	Electron capture in PHITS-KURBUC.
kurbuc_ioniz_e-exc_e-cap	Ionization, electronic excitation, and electron capture in PHITS-KURBUC.
kurbuc_e-stp	Electron stripping in PHITS-KURBUC.
kurbuc_hit	All reactions in PHITS-KURBUC.
its_elas	Elastic scattering in ITSART [3].
its_ioniz	Ionization in ITSART.
its_e-exc	Electronic excitation in ITSART.
its_ioniz_e-exc	Ionization and electronic excitation in ITSART.
its_hit	All reactions in ITSART.

[1]

PHITS-etsmode is track-structure mode in PHITS.

[2]

PHITS-KURBUC is proton and carbon track-structure mode in PHITS.

[3]

ITSART is generalized track-structure mode in PHITS.

From version 3.37, output can be defined by using reaction ID numbers. When reaction ID number 3000 is specified, the entire atomic-interaction category is tallied; it therefore covers a broader range than the atomic keyword. Specifically, it also includes events in which only electron multiple scattering or Rayleigh scattering (which the atomic keyword excludes) occurred; thus, 3000 is not equivalent to the atomic keyword. Likewise, reaction ID number 3009 covers Rayleigh scattering and electron multiple scattering, including events in which only these occurred. Reaction ID number 2000 tallies the entire nuclear-reaction category (including elastic scattering) but does not include optical photon interactions. Reactions not listed in this table can also be used by specifying their reaction ID numbers. Note, however, that reaction ID numbers of the form 2000+MT are evaluated from product-based MT estimation, so only MT = 2, 11, 16, 17, 22, 23, 24, 25, 28, 29, 30, 32, 33, 34, 35, 36, 37, 41, 42, 44, 45, 102–109, and 111–117 (i.e. 2002, 2011, 2016, …, 2117) actually match. Even within this range, other MT values (such as MT=4, 5, and 51–91) cannot be selected. In particular, 2003 (MT=3) cannot be selected by the product-based judgment and therefore cannot be used in place of the nuclear keyword. Unlike string keywords, reaction ID numbers can be specified as multiple numbers separated by spaces; in that case, the sum of the matched reactions (a logical OR) is output. For example, 3005 3006 tallies the sum of the photoelectric effect and Compton scattering reactions. To tally multiple reactions individually, define multiple [t-interact] sections and specify a reaction ID number in output for each section.

Table 7.14.23 gshow

value	explanation
0 (default), 1, 2, 3, 4, 5	When mesh=xyz and axis=xy,yz,xz, region border (1), material name (2), region name (3), and LAT number (4) are plotted using this option. gshow=5 outputs only material colors in pixel style when icntl=8.

Table 7.14.24 rshow

value	explanation
0 (default), 1, 2, 3	When mesh=xyz and axis=xy,yz,xz, region border (1), material name (2), and region name (3) are plotted using this option. A xyz mesh section must be added below this option.

Table 7.14.25 volume

value	explanation
(optional)	This option defines the volume for each region for reg mesh. Volume definitions are required below this option. Values defined in [volume] are used in the case of no definition. If special description such as using (  ) is used to specify a region in reg=, its internally defined region number is output in input echo when this volume subsection is not defined.
reg vol	Volume definition. For details see the [volume] section.

Table 7.14.26 iechrl

value	explanation
72 (default)	Maximum number of columns for volume input echo.

Table 7.14.27 resol

value	explanation
1 (default)	This option multiplies the region line resolution by a factor of resol with the gshow or rshow option set to define the line thickness.

Table 7.14.28 width

value	explanation
0.5 (default)	The option defines the line thickness.

Table 7.14.29 volmat

value	explanation
(optional, D=9)	This option corrects the volume value for each mesh when material is defined by the xyz mesh. 0 means no correction. Value of volmat is the number of scans per one xyz mesh side.

Table 7.14.30 epsout

value	explanation
0 (default), 1, 2	When epsout=1, results are plotted into eps files. The eps file is named by replacing the extension with .eps. When epsout=2, error bars are also displayed in the eps file, except for the 2-dimensional type, axis=xy, yz, xz, rz.

Table 7.14.31 bmpout

value	explanation
0 (default), 1	Generate Bitmap figure of 2-dimensional tally output. This file is named by replacing the extension with .bmp. When mesh=xyz and axis=xy, yz, xz, it is available.

Table 7.14.32 vtkout

value

explanation

0 (default), 1

Output the tally results in the xyz mesh in the input format of ParaView. This file is named by replacing the extension with .vtk. This option is available when mesh=xyz and axis=xy, yz, xz. In addition to the tally values, the VTK file contains the region number (region) and material number (material) assigned to each mesh cell. When gshow=0 (default), these numbers are determined from the geometry at the center coordinates of each mesh cell. Each mesh cell is evaluated independently, and this evaluation does not depend on the value of resol. This improves the accuracy near region and material boundaries and the robustness of VTK output for complex CSG and tetrahedral-mesh geometries. Note that this method judges whether the center of each cell lies inside a component, and thus components that are small compared with the mesh cells may be overlooked. To reduce such omissions, increase the number of meshes instead of resol. When gshow is set to 1 or greater, boundary geometry is also written in VTK format together with the region and material numbers. In this case, resol is used for the resolution of the boundary lines.

Table 7.14.33 vtkfmt

value	explanation
0 (default), 1	Format of output file for ParaView. 0: ascii, 1: binary.

Table 7.14.34 maxangel

value	explanation
Number of part (default)	Specify the number of part shown in the eps file. This parameter limits the number of particles shown in a tally eps file, while it does not limit them in a numerical data file.

Table 7.14.35 ctmin(i)

value	explanation
(optional, D=-9999)	Minimum value for the i-th counter.

Table 7.14.36 ctmax(i)

value	explanation
(optional, D=9999)	Maximum value for the i-th counter.

Table 7.14.37 chmin(i)

value	explanation
(optional, D=-9999)	Minimum value for the i-th history-counter. This parameter cannot be specified in the batch variance mode, istdev=1.

Table 7.14.38 chmax(i)

value	explanation
(optional, D=9999)	Maximum value for the i-th history-counter. This parameter cannot be specified in the batch variance mode, istdev=1.

Table 7.14.39 trcl

value	explanation
(optional)	Coordinate transformation number or definition for r-z or xyz mesh.

Table 7.14.40 gslat

value	explanation
2 (default)	Option to draw lattice or tetrahedral geometry boundary lines when gshow or rshow is specified.
0	Not drawing.
1	Drawing. Note that the line of voxel phantoms may not be drawn clearly.
2	Not drawing the lines in the same cell.
3	Not drawing the lines in the same material. In this case, the boundary of two adjacent cells filled with the same material is not also drawn, even if not lattice or tetrahedral geometry. When gshow>=3, cell and lattice numbers are now shown.

Table 7.14.41 stdcut

value	explanation
(optional, D=-1)	Threshold value of STD cut off.

When specifying stdcut, PHITS automatically stop the calculation depending on values of STD, standard deviation. This function is available when stdcut is positive and itall=0,1 is set in [parameters] section. When all relative values of STD of the tally result are larger than 0 and smaller than stdcut at the last of one batch, the calculation is stopped. If stdcut in two or more tally sections is set, all the results of the tally sections have to satisfy the conditions in order for the function to work.

Table 7.14.42 nlatcel

value	explanation
(optional)	By setting the cell number which has the LAT parameter, the number of interactions in cells that have the same cell number but are placed in different lattice coordinates are processed as different cells. It works only for MorP=prob.

Table 7.14.43 nlatmem

value	explanation
(optional, D=1000)	The maximum number of cells dealt with separately in one history.

When coll or reaction events belonging to coll are specified in [counter], the counter value changes before the particle information is tallied. Therefore, to tally the information at the moment the reaction event occurs, specify the counter value after the change.