Advanced Computing Platform for Theoretical Physics

Commit 5f40f970 authored by ralfulrich's avatar ralfulrich
Browse files

Merge branch 'master' of gitlab.ikp.kit.edu:AirShowerPhysics/crmc

parents e1a48ce5 12621678
-------------------------------------------------------------------------------
CRMC version 2.0.0 (08.2021)
-------------------------------------------------------------------------------
Add compatibility with HEPMC3 and direct link to RIVET.
New interface to link CRMC with GEANT4.
Improve precision of energy rescaling to avoid energy violation in lab frame at low energy from nuclear fragments for all models when necessary
Correct Hijing and DPMJET interface.
Improve build system.
-------------------------------------------------------------------------------
CRMC version 1.8.0 (01.2021)
-------------------------------------------------------------------------------
Update DPMJETIII versions (3.06, 3-17.1 and 3-19.1) and limit mass range to light nuclei to avoid problems with old DPMJET 3.0-6.
Update Sibyll to version 2.3d
Add a test on pz conservation and extend array size for cluster in EPOS to avoid stops.
-------------------------------------------------------------------------------
CRMC version 1.7.1 (10.2019)
-------------------------------------------------------------------------------
Extend array in EPOS for safe central heavy ion collisions (some stop otherwise)
Add possibility of sexaquark production in EPOS.
-------------------------------------------------------------------------------
CRMC version 1.7.0 (03.2018)
-------------------------------------------------------------------------------
added small tool to produce cross section tables
update to Sibyll2.3c
update to DPMJETIII-17.1
fix phi particle width in EPOS
-------------------------------------------------------------------------------
CRMC version 1.6.0 (04.2016)
-------------------------------------------------------------------------------
Update Sibyll to version 2.3 including charm production now.
keep original mass of QII to have correct energy conservation (since particles are put on-shell in double precision in crmc now, previous problems with mass in QII are probably solved now)
use the same cross-section in hepmc file and in screen output (sigine in hp and sigineaa in pA/AB)
updated event type definition for sibyll
correct bug in DPMJET interface (impact parameter problem)
-------------------------------------------------------------------------------
CRMC version 1.5.7 (04.2016)
-------------------------------------------------------------------------------
Same as v1.6.0 but with Sibyll 2.1
-------------------------------------------------------------------------------
CRMC version 1.5.6 (10.2015)
-------------------------------------------------------------------------------
Made ROOT-output work again
New option ("-S") to give centre-of-mass energy as input
When compiling, cmake automatically tries to use c++11
For DPMJET, some produced particles could make the program stop during decay
-------------------------------------------------------------------------------
CRMC version 1.5.5 (09.2015)
-------------------------------------------------------------------------------
Restored compatibility with f77 compiler (unsupported trim function)
ROOT dependency is now optional
New option ("-x") to only display cross section for model with set energy and exit
When "ibreit" option is used (for decays with Breit-Wigner masses) the energy rescaling is done after smearing of the energy
Energy balance corrected for QGSJetII (<1 GeV after using original particle masses)
For QGSJetII input energy given as E_lab not E_kin
Functionality for running HIJING model
-------------------------------------------------------------------------------
CRMC version 1.5.3 (07.2014)
-------------------------------------------------------------------------------
HepMC signal vertex propely set
When running DPMJET, use EPOS for decays that have a known Breit-Wigner width to obtain mother/daughter information
Bug correction, DPMJET used impact parameter of first event, always
-------------------------------------------------------------------------------
CRMC version 1.4 (10.2013)
-------------------------------------------------------------------------------
PYTHIA (and models depending on it) now use decay definition as set in crmc
Mother/daughter information can be adjusted for EPOS. istmax=0 only final state, istmax=1 all particles saved
CRMC v1.8.0 Last modifications 2020/12/01
CRMC v2.0.0 Last modifications 2021/08/23
**********************************************************
The Program "crmc"
......@@ -16,7 +16,7 @@ Version with extended arrays for heavy ion collisions.
Authors
**********************************************************
C++ interface : Colin Baus (colin.baus@kit.edu)
C++ interface : Colin Baus
Ralf Ulrich (ralf.ulrich@kit.edu)
Fortran interface : Tanguy Pierog (tanguy.pierog@kit.edu)
......@@ -41,37 +41,52 @@ Reference : C. Baus, T. Pierog and R. Ulrich.
* Post LHC :
EPOS LHC (-m0) : T. Pierog et al.
EPOS LHC (-m 0) : T. Pierog et al.
Phys.Rev. C92 (2015) 034906
QGSJETII-04 (-m7) : S.Ostapchenko,
QGSJETII-04 (-m 7) : S.Ostapchenko,
Phys.Rev. D83 (2011) 014018
SIBYLL2.3d (-m6) : F. Riehn, R. Engel, A. Fedynitch, T.K. Gaisser, and T. Stanev,
hep-ph: 1912.03300
SIBYLL2.3d (-m 6) : F. Riehn, R. Engel, A. Fedynitch, T.K. Gaisser, and T. Stanev,
hep-ph: 1912.03300
and
E.-J. Ahn, R. Engel, T.K. Gaisser, P. Lipari, and T. Stanev,
Phys.Rev. D80 (2009) 094003
E.-J. Ahn, R. Engel, T.K. Gaisser, P. Lipari, and T. Stanev,
Phys.Rev. D80 (2009) 094003
DPMJETIII 2017-1 : Fedynitch at al.
DPMJETIII 2017-1 (-m 12) : Fedynitch at al.
DPMJETIII 2019-1 : Fedynitch at al.
DPMJETIII 2019-1 (-m 12) : Fedynitch at al.
* Pre LHC :
DPMJET 3.06 : 10.1103/PhysRevC.77.014904 (not linked any more)
DPMJET 3.06 (-m 12) : 10.1103/PhysRevC.77.014904 (not linked any more)
EPOS 1.99 (-m1) : K. Werner, F.M. Liu and T. Pierog,
Phys.Rev. C74 (2006) 044902
EPOS 1.99 (-m 1) : K. Werner, F.M. Liu and T. Pierog,
Phys.Rev. C74 (2006) 044902
and
T. Pierog and K. Werner,
Nucl.Phys.Proc.Suppl. 196 (2009) 102-105
T. Pierog and K. Werner,
Nucl.Phys.Proc.Suppl. 196 (2009) 102-105
QGSJET01 (-m2) : N.N. Kalmykov, S. Ostapchenko, and A.I. Pavlov,
Nucl.Phys. B (Proc. Suppl.) 52B (1997) 17
QGSJET01 (-m 2) : N.N. Kalmykov, S. Ostapchenko, and A.I. Pavlov,
Nucl.Phys. B (Proc. Suppl.) 52B (1997) 17
QGSJETII-03 (-m11) : S. Ostapchenko,
Nucl.Phys.Proc.Suppl. 151 (2006) 143
QGSJETII-03 (-m 11) : S. Ostapchenko,
Nucl.Phys.Proc.Suppl. 151 (2006) 143
* Other models :
PHOJET 1.12 (-m 8) : R. Engel, J. Ranft and S. Roesler,
Phys. Rev. D52 (1995) 1459
and
F. W. Bopp, R. Engel and J. Ranft
arXiv:hep-ph/9803437.
PYTHIA 6.4.28 (-m 4) : T. Sjostrand, S. Mrenna and P. Skands,
JHEP 0605 (2006) 026
HIJING 1.38 (-m 5) : obsolete (please use a more modern version like HIJING++)
GHEISHA (-m 3) : obsolete
**********************************************************
Installation
......@@ -79,7 +94,9 @@ QGSJETII-03 (-m11) : S. Ostapchenko,
To install please get the latest version from the authors and install
the following pre-requisites: BOOST, HepMC3, CMAKE (*). Consider also
the extensions provided by Rivet (>v3) and fastjet.
the extensions provided by Rivet (>v3) and fastjet. We recommend to use
a LTS Ubuntu version to avoid compatibility problems (latest tests done
using 20.4).
In order to compile with Fastjet, set the Fastjet_DIR environment
variable to the location of your Fastjet installation.
......@@ -176,16 +193,16 @@ The environment variable $CRMC_OUT can be set to define the path the path of the
otherwise $PWD is used as default path.
**Example to generate 100 7 TeV pp collisions with EPOS LHC:
bin/crmc -o hepmc -S7000 -n100 -m0
bin/crmc -o hepmc -S 7000 -n 100 -m 0
**Example to generate 100 1.38 ATeV PbPb collisions with EPOS 1.99:
bin/crmc -o hepmc -p1380 -P-1380 -n100 -i208 -I208 -m1
bin/crmc -o hepmc -p 1380 -P -1380 -n 100 -i 208 -I 208 -m 1
**Example to generate 100 5.02 ATeV pPb collisions with QGSJetII-04:
bin/crmc -o hepmc -p4000 -P-1577 -n100 -m7 -i2212 -I822080
bin/crmc -o hepmc -p 4000 -P -1577 -n 100 -m 7 -i 2212 -I 822080
**Example to test Sibyll2.3
bin/crmc -T -m6
**Example to test Sibyll2.3d
bin/crmc -T -m 6
** On Rivet output
if you select Rivet as output CRMC will produce a yoda file with analysis output for you. Of course,
......@@ -203,7 +220,7 @@ In this file the 3 first commented options are valid for EPOS 1.99 and EPOS LHC
- By defaults EPOS models use a simplified treatment of QGP in events where the energy density is high enough (including in pp). If you
uncomment the line "!switch fusion off" this will be disable and running time will be shorter (but data description will be worth since
the physics model will be incomplete). In that mode EPOS is comparable to PYTHIA model (no final state interaction).
the physics model will be incomplete). In that mode EPOS is comparable to PYTHIA model without color reconnection (no final state interaction).
- By default only the final particles are recorded in the output file like for other cosmic ray models.
Uncommenting "!set istmax 1" allows the user to have the full chain of mother/daughter from the beam to the final particles with EPOS models.
......@@ -230,15 +247,14 @@ cross-section will be always correct but it takes several minutes to compute it
Note on Simulating Heavy Ion (HI) Events
******I****************************************************
QGSJET01 and SIBYLL2.3 are limited to hA collisions with A<64 (A<19) and h being a pion, a kaon or a nucleon.
QGSJET01 and SIBYLL2.3d are limited to hA collisions with A<64 (A<19) and h being a pion, a kaon or a nucleon.
Only EPOS (1.99 and LHC) and QGSJETII (03 and 04) models can run (h)AB collisions with A and B up to 208.
QGSJETII was never designed for HI collisions, so results should be interpreted with care (no final
state interactions).
QGSJETII was never designed for HI collisions, so results should be interpreted with care (no final state interactions).
EPOS 1.99 and LHC include a simplified treatment of final state interaction but doesn't include full hydro
simulations as already possible in EPOS 2. It gives a good overall description of HI data but it should not be
used to interpret QGP related observables (flow, jet quenching, etc ...) since the model is oversimplified in
simulations as already possible in EPOS 4. It gives a good overall description of HI data but it should not be
used to interpret QGP related observables (correlation, flow, jet quenching, etc ...) since the model is oversimplified in
this distribution.
* For detailed simulations of HI collisions within EPOS 2, please contact K. Werner (werner@subatech.in2p3.fr).*
* For detailed simulations of HI collisions within EPOS 4, please contact K. Werner (werner@subatech.in2p3.fr).*
Concerning the cross-section calculation in EPOS, to avoid unnecessary long time calculation needed
for HI interactions, the default type of cross-section calculation is fast and good for hA but not reliable for
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment