7.6. [ T-Deposit2 ] section

This tally scores deposit energy distributions in two regions and plots the correlation between two deposit energies. By this, one can simulate, for example, dE-E counters and plot the correlations in a 2-dimensional graph. In the calculation of deposition energies by ionization of charged particles, the same as [t-deposit], any factor can be multiplied using the user-defined subroutines usrdfn1.f and usrdfn2.f.

Table 7.6.1 mesh

value	explanation
reg	Mesh type. A reg mesh type subsection is required below this option.

Table 7.6.2 reg

value	explanation
r1 r2	Specify two cell numbers.

Table 7.6.3 part

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

Table 7.6.4 letmat1

value	explanation
(optional)	Material id for LET, \(dE/dx\), of region r1. If omitted, a real material is assumed.

Table 7.6.5 letmat2

value

explanation

(optional)

Material id for LET, \(dE/dx\), of region r2. If omitted, a real material is assumed. If a material not used in the geometry is selected, its material density must be defined in [material]. To calculate LET in water, define water with 1 g/cm^3 in [material]. When letmat1 < 0 or letmat2 < 0 is set, PHITS automatically calculates \(dE/dx\) for water with 1 g/cm^3 for electrons and positrons. Please see ParticleTherapy in the recommendation setting for more details.

Table 7.6.6 dedxfnc1

value	explanation
(optional, D=0)	For region r1. 0: without. 1: use usrdfn1.f. 2: use usrdfn2.f.

Table 7.6.7 dedxfnc2

value	explanation
(optional, D=0)	For region r2. 0: without. 1: use usrdfn1.f. 2: use usrdfn2.f. As examples, the default program of usrdfn1.f returns the dose equivalent calculated from deposit energy multiplied by the Q(L) relationship defined in the ICRP60, while that of usrdfn2.f estimates biological dose on the basis of the Microdosimetric Kinetic Model.

Table 7.6.8 e1-type

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

Table 7.6.9 e2-type

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

Table 7.6.10 t-type

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

Table 7.6.11 unit

value	explanation
1, 2	1: Number [1/source] 2: Number [1/nsec/source]

Table 7.6.12 axis

value	explanation
eng1, eng2, t	x axis value of output data.
e12, e21, t-e1, t-e2, e1-t, e2-t	2-dimensional.

Table 7.6.13 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.6.14 file

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

Table 7.6.15 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.6.16 factor

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

Table 7.6.17 title

value	explanation
(optional)	Title.

Table 7.6.18 angel

value	explanation
(optional)	ANGEL parameters.

Table 7.6.19 sangel

value	explanation
(optional)	Special format for ANGEL parameters.

Table 7.6.20 2d-type

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

Table 7.6.21 x-txt

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

Table 7.6.22 y-txt

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

Table 7.6.23 z-txt

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

Table 7.6.24 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.6.25 iechrl

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

Table 7.6.26 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.6.27 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.6.28 ctmin(i)

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

Table 7.6.29 ctmax(i)

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

Table 7.6.30 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.6.31 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.6.32 dresol1

value	explanation
(optional, D=0.0)	Parameter for representing detector resolution in tally region 1. When dresol1 = \(\sigma_r\) and dfano1 = \(F\), the deposition energy \(E\) of each event fluctuates following a Gaussian distribution with standard deviation \(\sqrt{\sigma_r^2 + FE}\). If dresol1 < 0, the energy resolutions, energy dependencies and spectral peak shapes, are overwritten by those defined in usrresol.f.

Table 7.6.33 dresol2

value	explanation
(optional, D=0.0)	Parameter for representing detector resolution in tally region 2. See the description of dresol1 for more details.

Table 7.6.34 dfano1

value	explanation
(optional, D=0.0)	Parameter for representing detector resolution in tally region 1. See the description of dresol1 for more details.

Table 7.6.35 dfano2

value	explanation
(optional, D=0.0)	Parameter for representing detector resolution in tally region 2. See the description of dresol1 for more details.

Table 7.6.36 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.

To tally the energy loss for each projectile particle entering the tally region, the counter should be defined using part in [counter], and ctmin and ctmax in this tally section.