.. _[t-point]:

[ T-Point ] セクション
==================================================

**[t-point]** タリーは、ある点（point detector）やリング状の線分
（ring detector）における fluence を計算することができます。
PHITS は、通常、設定した領域に確率的に入ってくる放射線をタリーするため、
タリー領域が小さくなると、十分な統計精度を得るためには極めて長い計算時間
を要する場合があります。
極端な例として、タリー領域を点や線分に設定した場合、どれだけ計算時間を費
やしても結果を得ることはできません。
**[t-point]** タリーは、 **[t-track]** や **[t-cross]** タリーの手法とは異
なり、線源位置、もしくは、散乱によって粒子が発生した地点で、タリー位置へ
の方向確率と透過確率を計算して評価することができるため、粒子の存在確率が
小さい領域での測定に有効に機能します。
ただし、散乱点での生成粒子の角度分布、エネルギー分布があらかじめ分かって
いないと評価できません。
また、物質中での荷電粒子のエネルギー損失があると透過確率の評価が難しいの
で、対象粒子は、核データの存在する中性子と光子に限られます。
したがって、 **[t-point]** タリーの利用条件は次の通りです。

1. 輸送する粒子はライブラリを使う上限エネルギー(dmax)以下の中性子・光子のみとする。 [#]_
2. 検出可能な粒子は中性子・光子のみとする。
3. Event Generator mode 及び EGS5 mode を使用しない。 [#]_
4. 面定義で全反射や白色反射の面を利用しない。

本タリーの計算原理や具体的な利用方法については、
``\phits\utility\tpoint`` フォルダにある readme やサンプルインプットをご
参照ください。

.. rst-class:: no-caption-number

.. list-table:: **point**
   :header-rows: 0

   * - value
     - explanation
   * - データ数
     - point detector の個数を指定するオプション。
       このオプションの下にサブセクションが必要です。

.. rst-class:: no-caption-number

.. list-table:: **ring**
   :header-rows: 0

   * - value
     - explanation
   * - データ数
     - ring detector の個数を指定するオプション。
       このオプションの下にサブセクションが必要です。

.. rst-class:: no-caption-number

.. list-table:: **part**
   :header-rows: 0

   * - value
     - explanation
   * - 粒子名（省略可）
     - タリーする粒子。
       省略時は **neutron photon** です。

.. include:: ./commontally/e-type.rst
.. include:: ./commontally/t-type.rst
.. rst-class:: no-caption-number

.. list-table:: **unit**
   :header-rows: 0

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

.. rst-class:: no-caption-number

.. list-table:: **axis**
   :header-rows: 0

   * - value
     - explanation
   * - **eng, t**
     - 出力データの x 軸。

.. rst-class:: no-caption-number

.. list-table:: **samepage**
   :header-rows: 0

   * - value
     - explanation
   * - (省略可, **D=part**)
     - 画像出力ファイルで同じページに表示するデータの種類を指定します。
       **axis** で指定できるパラメータに加えて、**point** と **ring** も指定できます。

.. include:: ./commontally/file.rst
.. include:: ./commontally/resfile.rst
.. include:: ./commontally/factor.rst
.. include:: ./commontally/title.rst
.. include:: ./commontally/angel.rst
.. include:: ./commontally/sangel.rst

**unit = 3, 13** の場合の Lethargy はエネルギーに関する自然対数目盛を表し
ています。
各エネルギービンの上限と下限がそれぞれ :math:`E_{\rm high}`,
:math:`E_{\rm low}` のときに、各 Lethargy の幅を
:math:`\ln(E_{\rm high}/E_{\rm low})` で与えます。

通常のタリーで指定する **mesh** パラメータの代わりに、 **[t-point]** タ
リーでは、 **point** もしくは **ring** パラメータを指定します。
まず、point detector の場合、入力書式は次のようになります。

.. code-block:: text
   :caption: point detector の定義例

     [ T-point ]
     point =    1     # number of point detectors
     non     x        y       z        r0
       1     10.0     0.0     50.0     1.0

**point=n** で point detector の個数 :math:`n` を指定します。
ひとつの **[t-point]** タリーセクションで定義できるポイントの上限数は 20
です。
これ以上のポイントを定義するときは、複数のセクションを用意してください。
次の行には、ポイントを指定する座標（x, y, z）、ポイント周辺の特異領域の
半径（r0）のデータ順を定義します。
読み飛ばしのコラムの指標 **non** も使えます。
これらのデータの並びを変えるときは、この行で指定してください。
この定義文の下に **point=** で指定した数だけデータを定義します。
座標は x, y, z 座標 [cm] で指定し、特異領域の半径（r0）も [cm] で指定して
ください。
ここで、特異領域の半径については、
``\phits\utility\tpoint`` フォルダにある readme をご参照ください。

次に、ring detector の場合は、次のような入力書式になります。

.. code-block:: text
   :caption: ring detector の定義例

     [ T-point ]
     ring =    1     # number of ring detectors
     non   axis    ar       rr       r0
       1    z      50.0     10.0     1.0

**ring=n** で ring detector の個数 :math:`n` を指定します。
次の行には、リングの軸を x, y, z の文字で指定する **axis** 、原点からリン
グの位置までの距離（ar）、リングの半径（rr）、ポイント周辺の特異領域の半
径（r0）のデータ順を定義します。
読み飛ばしのコラムの指標 **non** も使えます。
これらのデータの並びを変えるときは、この行で指定してください。
この定義文の下に **ring=** で指定した数だけデータを定義します。
リング軸は x, y, z の文字で指定し、リングまでの距離（ar）とリングの半径
（rr）は [cm] で指定します。
特異領域の半径（r0）も [cm] で指定してください。
ここで、特異領域の半径については、
``\phits\utility\tpoint`` フォルダにある readme をご参照ください。

.. rst-class:: no-caption-number

.. list-table:: **epsout**
   :header-rows: 0

   * - value
     - explanation
   * - **0** (default), **1**, **2**
     - **epsout=1** では、結果を eps ファイルとして出力します。
       ファイル名は拡張子を ``.eps`` に置き換えたものになります。
       **epsout=2** では、2次元タイプの **axis=xy, yz, xz, rz** を除いて、eps ファイルにエラーバーも表示します。

.. include:: ./commontally/maxangel.rst
.. include:: ./commontally/ctmin.rst
.. include:: ./commontally/ctmax.rst
.. include:: ./commontally/chmin.rst
.. include:: ./commontally/chmax.rst
.. include:: ./commontally/stdcut.rst
.. rst-class:: no-caption-number

.. list-table:: **multiplier**
   :header-rows: 0

   * - value
     - explanation
   * - 物質数
     - 物質ごとに multiplier を指定します。
   * - （省略可）
     - このオプションの下に multiplier サブセクションが必要です。
       詳細は multiplier サブセクションを参照してください。

.. rst-class:: no-caption-number

.. list-table:: **dump**
   :header-rows: 0

   * - value
     - explanation
   * - データ数
     - Dump機能により出力する粒子情報の数を定義します。
   * - 次行: データ定義
     - 出力する粒子情報の順番を :numref:`tbl-tally-dump` に示す粒子情報番号で定義します。

.. include:: ./commontally/iextstat.rst
.. include:: ./commontally/prodenmn.rst
.. include:: ./commontally/prodenmx.rst
.. include:: ./commontally/nbproden.rst


.. [#] Version 3.36より荷電粒子ライブラリや光核反応ライブラリも利用可能となりました。

.. [#] EGSモードを使うと光子フルエンスは正しく評価できませんが、光核反応ライブラリと組み合わせて中性子フルエンスを評価することは可能です。EGSモードを使わないと電子・陽電子の輸送を考慮できませんので、電子加速器からの中性子発生を評価する際は、EGSモードと光核反応ライブラリと組み合わせてご利用下さい。