.. _[t-heat]:

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

カーマ近似を用いて付与エネルギー（発熱）を出力するタリーです。
ただし、バージョン 3.05 以降 **[t-deposit]** でも同様の計算が可能となったため、それ以降のバージョンでは利用を奨励していません。

なお、イベントジェネレータモードを用いた場合、 **e-mode>=1** では中性子のカーマ近似は用いません。
また、 **electron=1** かつ電子の輸送を含んでいる時は、光子のカーマ近似は用いずに電子の付与エネルギーから計算します。

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

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

   * - value
     - explanation
   * - **reg, r-z, xyz**
     - メッシュ型。
       メッシュ型サブセクションが必要です。

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

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

   * - value
     - explanation
   * - **reg, x, y, z, r**
     - 出力データの :math:`x` 軸。
   * - **xy, yz, xz, rz**
     - 2次元表示。

.. include:: ./commontally/samepage.rst
.. include:: ./commontally/file.rst
.. include:: ./commontally/resfile.rst
.. include:: ./commontally/material.rst
.. rst-class:: no-caption-number

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

   * - value
     - explanation
   * - **heat**
     - **total**: 全付与エネルギー
       **leakage**: 外部ボイドに出ていった粒子の運動エネルギー
       **recoil**: 計算打切エネルギー **emin(15-19)** 設定時の残留核の運動エネルギー
       **ionization**: 荷電粒子のエネルギー損失による付与エネルギー
       **low neutron**: 中性子カーマ係数を用いた際の付与エネルギー
       **photon**: 光子カーマ係数を用いた際の付与エネルギー
       **electron=1** のときは、カットオフ電子による付与エネルギー
       **others**: 残留核の残留エネルギー
       **igamma=1** の場合、励起状態にある残留核から光子が放出されるのでゼロに近づきます。
   * - **simple**
     - **total**: 全付与エネルギー
       **leakage**: 外部ボイドに出ていった粒子の運動エネルギー
       **recoil**: 計算打切エネルギー **emin(15-19)** 設定時の残留核の運動エネルギー
       **ionization**: 荷電粒子のエネルギー損失による付与エネルギー
       **low neutron**: 中性子カーマ係数を用いた際の付与エネルギー
       **photon**: 光子カーマ係数を用いた際の付与エネルギー
       **electron=1** のときは、カットオフ電子による付与エネルギー
       **others**: 残留核の残留エネルギー
       **igamma=1** の場合、励起状態にある残留核から光子が放出されるのでゼロに近づきます。
   * - **all**
     - **simple** の結果に加えて、recoil に対する :math:`d`, :math:`t`, :math:`^3\mathrm{He}`, :math:`\alpha`, residual nuclei の寄与を書き出します。
       ionization に対する :math:`p`, :math:`\pi^+`, :math:`\pi^-`, others の寄与を書き出します。
       ionization に関しては **part** で指定した粒子の寄与も出力しますが、eps ファイルには出力されません。
       **stopped particle**: 物質中で停止した粒子の運動エネルギーの proton, neutron, :math:`\pi^+`, :math:`\pi^-`, others の寄与を書き出します。
       stopped particle に関しては **part** で指定した粒子の寄与も出力しますが、eps ファイルには出力されません。
       **others**: 残留核の励起エネルギーと fission の成分が出力されます。
       **axis** が 2 次元表示のときは、 **all** は **simple** と同じで、total, recoil, ionization, low neutron, electron, others が出力されます。

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

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

   * - value
     - explanation
   * - 粒子名 (省略可)
     - **output=all** の場合に、ここで指定した荷電粒子による付与エネルギーである ionization 成分と、物質中で停止した粒子のエネルギー量である stopped particle 成分を出力します。
       ただし、eps ファイルには出力されません。
       また、カーマ近似を使用した計算において光子と中性子による寄与を調べる場合、 **part=photon neutron** としてもこれらの寄与を区別できません。
       各々の寄与を求める場合は **output=simple** として出力される内容を確認してください。

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

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

   * - value
     - explanation
   * - **0, 1, 2**
     - **0**: [Gy/source]
       **1**: [MeV/cm :math:`^3` /source]
       **2**: [MeV/source]

.. include:: ./commontally/factor.rst
.. include:: ./commontally/title.rst
.. include:: ./commontally/angel.rst
.. include:: ./commontally/sangel.rst
.. include:: ./commontally/2d-type.rst
.. include:: ./commontally/x-txt.rst
.. include:: ./commontally/y-txt.rst
.. include:: ./commontally/z-txt.rst
.. include:: ./commontally/gshow.rst
.. include:: ./commontally/rshow.rst
.. include:: ./commontally/resol.rst
.. include:: ./commontally/width.rst

**[t-heat]** について全般的に言えることですが、最終的に熱に変換するエネルギーは、荷電粒子の ionization の過程を通してです。
しかしながら、輸送コードの中では cutoff エネルギーがあり、最後の過程まで追わずに輸送を終了します。
従って、熱の成分として recoil や others などの成分が出力されます。
これらの成分は、輸送のパラメータにより変化します。

.. include:: ./commontally/volume.rst
.. include:: ./commontally/iechrl.rst
.. include:: ./commontally/volmat.rst
.. include:: ./commontally/epsout.rst
.. include:: ./commontally/bmpout.rst
.. include:: ./commontally/vtkout.rst
.. include:: ./commontally/vtkfmt.rst
.. rst-class:: no-caption-number

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

   * - value
     - explanation
   * - **0** (default), **1**
     - 電子の寄与のオプションです。
       **0**: 光子のカーマ係数を使います。
       **1**: 電子のエネルギー付与から計算します。
       この場合、電子の輸送が必要です。

.. include:: ./commontally/ctmin.rst
.. include:: ./commontally/ctmax.rst
.. include:: ./commontally/chmin.rst
.. include:: ./commontally/chmax.rst
.. include:: ./commontally/trcl.rst
.. include:: ./commontally/gslat.rst

**unit=0** とした場合、 [Gy/source] を単位とする吸収線量が出力されます。
**mesh=reg** の場合は、 **volume** パラメータを用いるか、 **[volume]** セクションにおいて各領域の体積を設定する必要があります。
ただし、吸収線量は不可算量ですので、各 :math:`x` 軸に関する積分値は出力されません。
また、領域内に複数の物質が混在する場合、その領域の吸収線量は全体の平均値とはなりません。
例えば、領域 1 と 2 にある質量 :math:`M_1` と :math:`M_2` の物質にそれぞれ :math:`E_1` と :math:`E_2` の熱量が吸収された場合、その平均吸収線量は :math:`\frac{E_1+E_2}{M_1+M_2}` ですが、PHITS で計算する値は :math:`\frac{E_1}{M_1}\frac{V_1}{V_1+V_2}+\frac{E_2}{M_2}\frac{V_1}{V_1+V_2}` となります。
ここで、 :math:`V_1` と :math:`V_2` はそれぞれ領域 1 と 2 が占める体積を表します。

.. include:: ./commontally/stdcut.rst