.. _[t-4dtrack]: [ T-4Dtrack ] section ================================================== This tally outputs particle track information (position and time). By loading the output file into the **Particle tracks** tab of PHIG-3D, the trajectories can be visualized on the geometry. .. rst-class:: no-caption-number .. list-table:: **file** :header-rows: 0 * - value - explanation * - ``file name`` - Output file name. .. rst-class:: no-caption-number .. list-table:: **reg** :header-rows: 0 * - value - explanation * - **all** (default) - Specify the region number for writing out tracks. Refer to the common tally parameter description for the format. .. rst-class:: no-caption-number .. list-table:: **part** :header-rows: 0 * - value - explanation * - **all** (default) - Specify the particles for which to write out tracks. If not specified, tracks for all particles will be written out. If **part** includes **nucleus**, unstable nuclei will also be included in the output. However, only nuclei with half-lives between **HalfLifeMin** and **HalfLifeMax** will be considered. Refer to :numref:`sec-tallypara-part` for the format. .. rst-class:: no-caption-number .. list-table:: **stepminEM** :header-rows: 0 * - value - explanation * - 0.0 (default) - Minimum step distance (cm) for electrons, positrons, and photons. Adjusting this value can reduce the data size significantly without compromising the quality of the visualizing using PHIG-3D. See the following 1 for details. .. rst-class:: no-caption-number .. list-table:: **SkipRate** :header-rows: 0 * - value - explanation * - 0.0 (default) - Specify the ratio for skipping the output of particle track endpoints. This setting applies to all particles. Set a value between 0 and 1. For example, if **SkipRate=0.9** is specified, the output of endpoints will be skipped with a 90% probability. The size of the output file will be approximately (1 - **SkipRate**) times the original size. For exceptions to skipping, see 2 below. .. rst-class:: no-caption-number .. list-table:: **SkipRate-of(kf-code)** :header-rows: 0 * - value - explanation * - 0.0 (default) - Apply **SkipRate** only to particles specified by kf-code. For example, if **SkipRate-of(2112)=0.9** is specified, the output of endpoints for neutrons will be skipped with a 90% probability. If both **SkipRate** and **SkipRate-of** are set simultaneously, the value of **SkipRate-of** takes precedence. 1. A particle track is essentially composed of continuous steps. The mean free paths of electrons, positrons, and photons are generally very small, so writing out every step would result in a very large data size. Therefore, to save data size, steps whose distance from the starting point of a step is less than **stepminEM** are not written out. 2. The output of track points is not skipped if any of the following three conditions are met: 1. The particle reaches the energy cutoff. 2. The particle reaches the boundary of a region. 3. An interaction occurs, and the number of final states is 2 or more. .. rst-class:: no-caption-number .. list-table:: **HistoryMax** :header-rows: 0 * - value - explanation * - None (default) - The maximum number of events to write out. If omitted, all events are written out. .. rst-class:: no-caption-number .. list-table:: **HistoryFile** :header-rows: 0 * - value - explanation * - ``file name`` - Specifies the file name containing the event numbers to be written out. If not specified, all events are written out. File structure: Write the **NOCAS** (event number) and **NOBCH** (batch number) of the desired event on each line (Example 1). The numbers can be read in real number format, so dump data files specifying **NOCAS** and **NOBCH** can also be loaded. .. rst-class:: no-caption-number .. list-table:: **format** :header-rows: 0 * - value - explanation * - **t4d** (default) - Specify the output format. Output in **t4d** format. For details, refer to ``/phits/utility/t-4dtrack/T4D-format/readme-en.pptx``. * - **pict** - Output in **pict** format used by EGS5. .. rst-class:: no-caption-number .. list-table:: **verbose** :header-rows: 0 * - value - explanation * - 0 - Specify the output level. For details, refer to ``/phits/utility/t-4Dtrack/format/readme-en.pptx``. Output the minimum necessary information for track display in PHIG-3D. * - 1 (default) - Also output track cell information and other details. For examples of how to use cell information, see ``/phits/utility/t-4Dtrack/OriginTrackFinder/readme-en.pptx``. * - 2 - Output counter information and more. When outputting counter values, it is necessary to use **ctmin** or **ctmax** in a separate tally. .. rst-class:: no-caption-number .. list-table:: **precision** :header-rows: 0 * - value - explanation * - 6 (default) - Number of digits for outputting real numbers. An example of the input for the file specified by **HistoryFile**: .. code-block:: text :caption: Example of **HistoryFile** input # NOCAS NOBCH 1 1 2 1 3 1 1 2 2 2 3 2 .. rst-class:: no-caption-number .. list-table:: **PinnedOption** :header-rows: 0 * - value - explanation * - 0 (default) - Choose whether to keep (pin) a specified particle in the visualization after it stops due to falling below the energy cutoff, instead of immediately removing it. The duration for which it remains can be specified. Do not pin. * - 1 - Pin generated nuclei. The target nuclei can be specified by their half-life. * - 2 - Pin particles specified in **PinnedParticles**. .. rst-class:: no-caption-number .. list-table:: **PinnedTimeOption** :header-rows: 0 * - value - explanation * - 0 (default) - Specify the type of duration for pinning. Fixed duration. Specified by **PinnedTimeFix**. * - 1 - **aHalf** x (the half-life of an unstable nucleus) * - 2 - max(0, **aHalf** x log10(half-life/s)+**bHalf**) (ns) .. rst-class:: no-caption-number .. list-table:: **PinnedTimeFix** :header-rows: 0 * - value - explanation * - 10.0 (default) - Fixed pinning duration (ns). Used with **PinnedTimeOption = 0**. .. rst-class:: no-caption-number .. list-table:: **aHalf** :header-rows: 0 * - value - explanation * - 1.0 (default) - Scale parameter for half-life. Used with **PinnedTimeOption = 1,2**. .. rst-class:: no-caption-number .. list-table:: **bHalf** :header-rows: 0 * - value - explanation * - 0.0 (default) - Scale parameter for the logarithmic term of the half-life. Used with **PinnedTimeOption = 2**. .. rst-class:: no-caption-number .. list-table:: **HalfLifeMin** :header-rows: 0 * - value - explanation * - 1.0e-22 s (default) - Minimum half-life for targeted unstable nuclei. Used with **part = nucleus**, **PinnedTimeOption = 1,2**. Enter a number followed by at least one space, then enter the unit of time. Available units are **s, min, h, d, y**, representing seconds, minutes, hours, days, and years, respectively. .. rst-class:: no-caption-number .. list-table:: **HalfLifeMax** :header-rows: 0 * - value - explanation * - 3.0e+32 s (default) - Maximum half-life for targeted unstable nuclei. Used with **part = nucleus**, **PinnedTimeOption = 1,2**. The method of entering the value is the same as for **HalfLifeMin**. .. rst-class:: no-caption-number .. list-table:: **PinnedParticles** :header-rows: 0 * - value - explanation * - None (default) - Enter the names or kf codes of the particles to be pinned. If entering multiple, separate them with spaces.