.. _[t-deposit2]:

[ 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``.

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

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

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

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

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

   * - value
     - explanation
   * - **r1  r2**
     - Specify two cell numbers.

.. include:: ./commontally/part.rst
.. rst-class:: no-caption-number

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

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

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

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

   * - value
     - explanation
   * - (optional)
     - Material id for LET, :math:`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 :math:`dE/dx` for water with 1 g/cm^3 for electrons and positrons.
       Please see ``ParticleTherapy`` in the recommendation setting for more details.


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

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

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

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

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

   * - 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.

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

.. list-table:: **e1-type**
   :header-rows: 0

   * - value
     - explanation
   * - **1, 2, 3, 4, 5**
     - Energy mesh.
       For region **r1**, See :numref:`sec-mesh-type` for mesh subsection format.

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

.. list-table:: **e2-type**
   :header-rows: 0

   * - value
     - explanation
   * - **1, 2, 3, 4, 5**
     - Energy mesh.
       For region **r2**, See :numref:`sec-mesh-type` for mesh subsection format.

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

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

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

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

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

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

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

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

   * - 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.

.. include:: ./commontally/file.rst
.. include:: ./commontally/resfile.rst
.. 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/volume.rst
.. include:: ./commontally/iechrl.rst
.. include:: ./commontally/epsout.rst
.. include:: ./commontally/maxangel.rst
.. include:: ./commontally/ctmin.rst
.. include:: ./commontally/ctmax.rst
.. include:: ./commontally/chmin.rst
.. include:: ./commontally/chmax.rst
.. rst-class:: no-caption-number

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

   * - value
     - explanation
   * - (optional, **D=0.0**)
     - Parameter for representing detector resolution in tally region 1.
       When **dresol1** = :math:`\sigma_r` and **dfano1** = :math:`F`, the deposition energy :math:`E` of each event fluctuates following a Gaussian distribution with standard deviation :math:`\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``.

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

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

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

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

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

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

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

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

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

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

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.