Back to home page

sPhenix code displayed by LXR

 
 

    


Warning, /JETSCAPE/external_packages/trento/README.rst is written in an unsupported language. File is not indexed.

0001 ====================
0002 T\ :sub:`R`\ ENTo 3D
0003 ====================
0004 
0005 *Three-dimensional Extended Reduced Thickness Event-by-event Nuclear Topology*
0006 
0007 A three-dimensional initial condition model for high-energy nuclear collisions. 
0008 For a documentation of the original boost-invariant T\ :sub:`R`\ ENTo model, please read the docs at `qcd.phy.duke.edu/trento <http://qcd.phy.duke.edu/trento>`_. For model details please refer to `PRC 92 011901 <https://doi.org/10.1103/PhysRevC.92.011901>`_, `PRC 96 044912 <https://doi.org/10.1103/PhysRevC.96.044912>`_.
0009 
0010 |
0011 
0012    |image1| |image2|
0013 
0014    **Figure 1**: *An event from* T\ :sub:`R`\ ENTo 3D. *Left: the event projected on to* :math:`x-y` *plane at midrapidty; right: same event projected on to* *y*\ -η *plane.*
0015 
0016 |
0017 |
0018 
0019 **Installation**
0020 
0021 First, install the dependencies: a C++11 compiler, Boost, GSL and HDF5. Then download the source code and compile with CMake:
0022 
0023 .. code-block:: shell
0024 
0025    mkdir build && cd build
0026    cmake ..
0027    make install
0028 
0029 |
0030 |
0031 
0032 **Usage**: 
0033 
0034 T\ :sub:`R`\ ENTo 3D keeps all the options in the original boost-invariant model with additional options listed in Table 1. *Note that to work in three-dimensional mode*, :code:`--eta-max` *must be set to a positive value.*
0035 
0036 Examples:
0037 
0038 * Genereate 10 Pb Pb events at *s*\ :sup:`1/2` = 5020 GeV in boost-invariant mode and output to a :code:`.hdf5` file
0039 
0040 .. code-block:: shell
0041 
0042    trento3d Pb Pb 10 -e 5020 -o PbPb.hdf5
0043 
0044 * Same as above but generate three-dimensional initial condition with -10<η<10, *d*\ η=0.2.
0045 
0046 .. code-block:: shell
0047 
0048    trento3d Pb Pb 10 -e 5020 --eta-max=10.0 --eta-step=0.2 -o PbPb.hdf5
0049 
0050 * Use the absolute-skewness parametrization (see Table 2) instead of the relative-skewness parametrization, with skew coefficient γ\ :sub:`0`\ = 1.0.
0051 
0052 .. code-block:: shell
0053 
0054    trento3d Pb Pb 10 -e 5020 -r 2 -t 1.0 --eta-max=10.0 --eta-step=0.2 -o PbPb.hdf5
0055 
0056 .. csv-table:: **Table 1**: Additional program options
0057    :header: "Options", "Default", "Description"
0058    :widths: 10, 10, 35
0059    :align: center
0060 
0061    "-m, --mean-coeff", 1.0 (float>0), "rapidity mean coefficient *μ*\ :sub:`0` "
0062    "-s, --std-coeff", 3.0 (float>0), "rapidity std coefficient *σ*\ :sub:`0`"
0063    "-t, --skew-coeff", 0.0 (float>0), "rapidity skew coefficient *γ*\ :sub:`0`"
0064    "-r, --skew-type", 1 (int), "
0065                                         1 = relative skewness
0066 
0067                                         2 = absolute skewness
0068                         
0069                                         else = no skewness"
0070    "-j, --jacobian", 0.8 (float>0), "<\ *p*\ :sub:`t`\ /\ *m*\ :sub:`t`\ > used in Jacobian *dy/d*\ η"
0071    "-e, --beam-energy", 2760 (float>0), "collision beam energy  *s*\ :sup:`1/2` [GeV], initializes cross section"
0072    "--xy-max",  10.0 (float) , "transverse x [fm] and y [fm] maximum (x,y grid from -max to +max)"
0073    "--xy-step",  0.2 (float), "transverse x [fm] and y [fm] step size"
0074    "--eta-max",  0.0 (float) , "space-time rapidity maximum (η grid from -max to +max)"
0075    "--eta-step",  0.5 (float), "space-time rapidity step size"
0076 
0077 |
0078 |
0079 
0080 **Longitudinal extension**: 
0081 
0082 T\ :sub:`R`\ ENTo 3D reproduces T\ :sub:`R`\ ENTo at midrapidity (η=0) exactly. At finite space-time rapidity, the entropy production is the product of its midrapidity value and a longitudinal profile function that varies at each transvese location. The profile is characterized by its first η-cumulants: mean, stadard deviationa and skewness. They are parametrized in terms of nuclear thickness function:
0083 
0084 .. csv-table:: **Table 2**: cumulant parametrization
0085    :header: "Cumulants", "Parametrization"
0086    :widths: 15, 30
0087    :align: center
0088 
0089    "mean", "\ *μ*\ :sub:`0`\ /2 log [(\ *T*\ :sub:`A` *e*\ :sup:`Y`\ + \ *T*\ :sub:`B` *e*\ :sup:`-Y` ) / (\ *T*\ :sub:`A` *e*\ :sup:`-Y`\ + \ *T*\ :sub:`B` *e*\ :sup:`Y`\ )], *Y* is the beam rapidity"
0090    "standard deviation", *σ*\ :sub:`0`
0091    "skewness",  "Relative skewness, *γ*\ :sub:`0` (\ *T*\ :sub:`A`\ - \ *T*\ :sub:`B`\  )/(\ *T*\ :sub:`A`\ + \ *T*\ :sub:`B`\ )
0092    
0093    Absolute skewness, *γ*\ :sub:`0` (\ *T*\ :sub:`A`\ - \ *T*\ :sub:`B`\  )"
0094 
0095 .. |image1| image:: doc/_static/event.png
0096    :width: 30%
0097 .. |image2| image:: doc/_static/event-eta.png
0098    :width: 30%
0099