.. _sec-libout:

[ Libout ] section
==================================================

The [libout] section is used to specify parameters for the cross-section and kerma factor output mode ( **inucr = 100** ) when the cross-section and kerma factor output mode ( **icntl = 1** ) is selected.

Specify the parameters in the following format.

.. code-block:: text

  [ Libout ]
    part = neutron
    eneinc = 100
    e-type = 3
    ne = 50
    emin = 1
    emax = 100
    epsout = 1
    a-type = -3
    na = 9
    12.5 17.5
    27.5 32.5
    42.5 47.5
    57.5 62.5
    72.5 77.5

The parameters that can be specified are as follows.

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

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

  * - Value
    - Explanation
  * - Number (optional) 208Pb
    - Specify the number of nuclides for which cross sections are output. If you want to output all the nuclides written in material number 1, you do not need to set this parameter.

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

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

  * - Value
    - Explanation
  * - proton, deuteron, triton, 3He, alpha, photon
    - Specify the emitted particle for outputting the differential cross section; when part=photon is specified, the differential cross section is output for the same gamma ray as when gamout=1.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - value (D=15)
    - Specify the incident particle energy [MeV/u] at which to output the cross section.

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

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

  * - Value
    - Explanation
  * - 0, 1 (D=1)
    - Output input echoes of the [libout] section and basic information on libraries.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - 0, 1, 2 (D=2)
    - When a non-zero number is specified, 1: differential cross section [mb/MeV] for the energy of the emitted particle, 2: double differential cross section [mb/sr/MeV] for the energy and angle of the emitted particle are output. e-type must be specified for 1, and e-type and a-type must be specified for 2. maxcas should be set to a large number, as Monte Carlo method is used for this calculation.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - 1, 2, 3, 4, 5
    - Energy mesh for emitted particles. See :numref:`sec-mesh-type` for mesh subsection format.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - 1, 2, 3, 4, -1, -2, -3, -4
    - Angle mesh of emitted particles (1, 2:cos, -1, -2:degree). See :numref:`sec-mesh-type` for mesh subsection format. when 3, -3 is specified, PHITS reads the angle segment points in the same format as for 1, -1 and output results only for odd-numbered angle bins; when 4, -4 is specified, PHITS reads angle segment points :math:`a(i)` in the same format as for 1, -1 and output differential cross sections for angle bins :math:`a(i) \pm \mathrm{delcos}` .

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

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

  * - Value
    - Explanation
  * - value (D=0.015)
    - Specify the angle bin width (cos) when a-type=4 or -4 is specified.

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

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

  * - Value
    - Explanation
  * - Number (D=0)
    - If a non-zero number is specified, the cross section for each MT is output; the MT is specified by number in the next line.

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

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

  * - Value
    - Explanation
  * - 0, 1 (D=0)
    - When 1 is specified, the differential cross section [mb/MeV] for the energy of the emitted gamma rays is output; when part=photon, it is not necessary to specify this parameter.

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

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

  * - Value
    - Explanation
  * - 1, 2, 3, 4, 5
    - Energy mesh for emitted gamma radiation. See :numref:`sec-mesh-type` for mesh subsection format.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - 0, 1, 2, -1, -2 (D=2)
    - If a non-zero number is specified, the cross section [b] of the integral value is output. 1: total only, 2: total, absorption, elastic, non-elastic When positive, PHITS outputs the cross section at the energy point in the library closest to the energy point specified in ei-type. When negative, cross sections are output at all energy points given in the library.


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

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

  * - Value
    - Explanation
  * - 1, 2, 3, 4, 5
    - Specify the incident energy point [MeV/u] at which to output the cross section of the integral value. See :numref:`sec-mesh-type` for mesh subsection format.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - filename
    - When the cross section to be output is one of integral cross section, differential cross section, or gamma differential cross section, the cross section is output to a file named filename. When the data of the type more than one are output, the output file name is given as  filename_sig.xxx, filename_ddx.xxx, or filename_gam.xxx.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - 0, 1 (D=0)
    - When 1 is specified, parameters for ANGEL drawing are added to the output file.

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

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

  * - Value
    - Explanation
  * - ANGEL parameter
    - Specify the ANGEL parameter in output of the integral cross section.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - ANGEL parameter
    - Specify the ANGEL parameter in output of the differential cross section.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - ANGEL parameter
    - Specify the ANGEL parameter in output of the differential cross section of gamma emissions.
    
.. rst-class:: no-caption-number

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

  * - Value
    - Explanation
  * - 0, 1 (D=0)
    - When 1 is specified, the results of the double-differential cross section for each angle are changed and output together in one table.

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

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

  * - Value
    - Explanation
  * - 0, 1 (D=0)
    - When 1 is specified, eps file is automatically generated.

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

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

  * - Value
    - Explanation
  * - filename
    - When specified, the differential cross section table is read from the filename file, and the read data are output side by side when outputting differential cross sections. This is used, for example, to create a figure comparing multiple cross sections. Supported files are PHITS tally output files with axis = eng and ddx files created with this function. Both the energy mesh and the angle mesh in the file must match those in another file to be compared.