Skip to content
/ PACIAE Public

A parton and hadron cascade model PACIAE (Parton And-hadron China Institute of Atomic Energy).

License

GPL-2.0 license
4 stars 1 fork Branches Tags Activity
Star
Notifications You must be signed in to change notification settings

ArcsaberHep/PACIAE

BranchesTags

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

68 Commits
.gitignore
.gitignore
 
 
LICENSE
LICENSE
 
 
PACIAE.sh
PACIAE.sh
 
 
README.md
README.md
 
 
analy.f
analy.f
 
 
coales_30.f
coales_30.f
 
 
hadcas_30.f
hadcas_30.f
 
 
main_30.f
main_30.f
 
 
p_30.f
p_30.f
 
 
parcas_30.f
parcas_30.f
 
 
parini_30.f
parini_30.f
 
 
rms_analysis.f90
rms_analysis.f90
 
 
sfm_30.f
sfm_30.f
 
 
usu.dat
usu.dat
 
 

Repository files navigation

The parton and hadron cascade model PACIAE

PACIAE model (Parton And-hadron China Institute of Atomic Energy) is a multipurpose Monte Carlo event generator developed to describe a wide range of collisions, including lepton-lepton, lepton-hadron, lepton-nucleus, hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions. It is built based on PYTHIA-6.428 and incorporates parton and hadron rescattering stages to take care of the nuclear medium effects.

Installation

Just download the file and decompress it. Then you can get the PACIAE source code directly.

Usage

We encourage users to run the program on LINUX. There are two ways to run the program.

  1. Use the normal compiler to complie the soruce code and run the PACIAE program. Take GFortran on LINUX as an example:

    • Compile and link the programs by the command: 
          gfortran -O -C *.f -o paciae.x
    • Modify the input file of usu.dat according to your wish.
    • Run the program by the command: 
          time ./paciae.x
      One could use the following command to run the program in the background and record the log information. 
          nohup time ./paciae.x > paciae.log &

  2. Use the PACIAE.sh shell-script to compile and link the code, generate usu.dat file and run program automatically.

    • Modify the PACIAE.sh file as needed.
    • Give executable permission to the PACIAE.sh file by command: 
          chmod +x PACIAE.sh
    • Run the PACIAE.sh script by the command: 
          time ./PACIAE.sh
      The more recommended command is: 
         time ./PACIAE.sh | tee $(date "+%Y%m%d%H%M%S").log
      which records the screen information to a log file with date and time. Of course, run them in the background: 
         nohup time ./PACIAE.sh &

         nohup time ./PACIAE.sh | tee $(date "+%Y%m%d%H%M%S").log &
    • "rms_analysis.f90" is a stand-alone program used in conjunction with the PACIAE.sh script to average the internal analysis files "rms.out" generated by each PACIAE simulation.

    It is also worth mentioning that tasks can be submitted to computer clusters and supercomputers using the PACIAE.sh script. (It is currently only available in SLRUM and LSF scheduling systems) More detailed information and usage please read the PACIAE.sh file.

Maintainers

@ArcsaberHep @ArcsaberkxL

Contributing

Feel free to dive in! Any bug reports, comments and suggestions are welcome. Please do not hesitate to contact us.

Contributors

License

GPL v2.0

Released programs

The code now are hosted on https://github.com/ArcsaberHep/PACIAE and https://gitee.com/arcsaberhep/PACIAE.

The released code are available on https://github.com/ArcsaberHep/PACIAE/releases and https://gitee.com/arcsaberhep/PACIAE/releases.

Released papers





Relevant papers

Update Notes:

...(waiting for updating)

01-03/2024: In version PACIAE 3.0.01

  • In "coales_30.f", a new coalescence scheme has been introduced by Yu-Liang with some fixes and improvements from Ben-Hao and An-Ke, which gives faster and more direct coalescence. A new parameter "bmrat" (ratio of baryon to meson) was introduced.
  • In "hadcas_30.f", the heavy flavored D meson channels (D pi, D rho and DN) have been added. The AQM cross sections were adopted. Four subroutine "DMeson_coll, Dpn_coll, DMeson1, Dmeson2" and one function "sigma_AQM" were added. The "BLOCK DATA PYCIDATA" and related comments were fixed and updated.
  • "rms0.out" output updated.
  • Some comments updated in program.
  • "usu.dat" updated.
  • "PACIAE.sh" updated.
  • Bug fixed.

PACIAE 3.0 is now released!

12/2023: In version PACIAE 3.0

  • The name "-loop" was renamed as "-framework".
  • Some comments updated in program.
  • Bug fixed.

11/2023: In version PACIAE 3.0

  • In subroutines "break, break_sbe" of "parini_30.f, main_30.f", statements have been rewritten to get KF codes from corresponding diquark directly.
  • Some debug statements were removed in "parini_30.f" and "parcas_30.f".
  • The kinematical masses of quarks were used in PACIAE3.0 now.
  • Redundant statements deleted / modified.
  • Bug fixed.

10/2023: In version PACIAE 3.0

  • In subroutine "xevent, updtlp" of "parini_30.f", bugs were fixed for loop-B to avoid potential dead-loops.
  • Redundant statements deleted / modified.
  • Bug fixed.

09/2023: In version PACIAE 3.0

  • In subroutine "xevent" of "parini_30.f", the frame "5MOM" was introduced to treat the collisions between those possible hadrons from diffractive events.
  • The program was re-named to version PACIAE 3.0.
  • "usu.dat" updated.
  • "PACIAE.sh" updated.
  • Redundant statements deleted / modified.
  • Bug fixed.

08/2023: In version PACIAE 2.3

  • In "main_23.f", the primary long-lived particles were defined according to ALICE-PUBLIC-2017-005. A new array "KF_woDecay" was introduced to specify particles not to decay from "usu.dat" directly.
  • The subroutine "break_f" of "coales_23.f" has been rewritten to use probability ratios directly. COMMON/sa38/ has been modified, correspondingly.
  • In COMMON/sa25/, the parameter "mstj1_1" and "mstj1_2" were renamed/replaced as "i_inel_proc" and "i_time_shower", respectively.
  • In "main_23.f" and "analy.f", the parameter "parp22" was renamed/replaced as "yOrEta". "parp22" was removed from "parini_23.f". The parameter "mstj2" was renamed/replaced as "iMode".
  • In subroutine "sysini" of "parini_23.f", the fitting formula of the energy dependent NN total cross section was introduced, according to PRC 97 (2018) 054910, Constantin Loizides et al.
  • An individual analyzing program "rms_analysis.f90" was introduced to average multiple "rms.out" files when "PACIAE.sh" was used.
  • The output format of "rms.out" were optimized.
  • In subroutine "break_f" of "coales_23.f", the probabilities of "u:d:s" were replaced by Quantum Tunneling probabilities.
  • In subroutine "oscar" of "main_23.f", a bug was fixed for the case nosc=0.
  • In subroutine "tabh" of "coales_23.f", the probabilities of mixing-state measons were corrected.
  • PACIAE.sh updated. A "sim" folder (simulation) would be generated to storing the results of simulations.

07/2023: In version PACIAE 2.3

  • In "main_23.f" and "parini_23.f", the parameters "dtt" and "smtj3" were renamed/replaced as "x_ratio" and "decpro" directly.
  • In "parini_23.f", several subroutines about leading-proton reconstruction were introduced. TODO: It still need detailed and complete treatment.
  • In "PAPTDI", "deexcitation_E", and "deexcitation_EP" of "coales_23.f", minor optimizations were made.
  • In "coales_23.f", two new implementations from "deexcitation_EP" were introduced, in which the excited quark-antiquark would be local pT compensation. TODO: It still need detailed and complete treatment.
  • "usu.dat" updated.
  • "PACIAE.sh" updated.
  • Redundant statements deleted / modified.
  • Bug fixed.

06/2023: In version PACIAE 2.3

  • In "main_23.f", three input parameter "parp82", "i_coord_recover", "i_tune" were introduced. "parp82" corresponds to PARP(82) in PYTHIA. "i_coord_recover" controls recovering the positions of partons after parton rescattering to those before or not. "i_tune" corresponds to MSTP(5) in PYTHIA which gives easy access to different tunes encoded in PYTHIA.
  • In "main_23.f", the 4-position of one of the rest partons would be assigned to the first parton in the string. This treatment would give random 3-coordinates to produced hadrons that surround the first parton after PYTHIA SFM, i.e. more random position distribution for produced hadrons.
  • In "main_23.f", the output format of "rms0.out" was optimized.
  • The subroutine "oscar" of "main_23.f", the OSCAR output was optimized.
  • In "main_23.f", "parini_23.f", and "analy.f", the subroutine "PASTAT" was introduced to output parton-parton level cross-sections to "main.out" file. However, one should note that it is just for the basic check.
  • In "parini_23.f", the leading-proton reconstruction was improved. TODO: It still need detailed and complete treatment.
  • In subroutine "xevent" of "parini_23.f", now when simulate pA and Ap collision, the junction-type NN event would be re-generated.
  • In subroutine "bream" of "parini_23.f", the broken quarks would be massless if breaking diquark is massless from PYTHIA treatment.
  • In subroutines "updtlp" and "updatl" of "parini_23.f", minor optimizations were made.
  • In subroutines "ctlcre_par", "ctlcre_para", "his_p", and "updpli_p" of "parcas_23.f", minor optimizations were made.
  • In subroutines "ctlcre_h" of "hadcas_23.f", one minor optimization was made.
  • "usu.dat" updated.
  • "PACIAE.sh" updated.
  • Redundant statements deleted / modified.
  • Bug fixed.

05/2023: In version PACIAE 2.3

  • In "main_23.f" and "analy.f", the internal online analyzing module was improved and optimized. The 5 distributions was extended to 6, i.e. the transverse momentum spectra dN/dpT. The format of user output file "rms.out" was optimized, too.
  • In ""main_23.f" and "parini_23.f", the diffractive NN event without parton generation after PYTHIA calling would not be thrown away now.
  • In subroutine "scat" of "parini_23.f", a minor bug of "m1 -> mm1" was fixed.
  • In "parcas_23.f", the max simulated volume was introduced, controlled by "PARAM(10)" (para10).
  • In "parcas_23.f", the inelastic processes and time-like radiation was fixed and improved to consider special pure-gluon strings (glue) and junction-type strings properly. TODO: It still need detailed and complete treatment.
  • In "analy.f", the usbroutine "analy_parton" and "stati_parton" were introduced to analyze partons. The subroutines "output_hadron_distribution" and "output_parton_distribution" were introduced to output analyzing information of hadrons and partons, respectively. The analyzing module was improved and optimized to 6 distributions.
  • "usu.dat" updated.
  • "PACIAE.sh" updated.
  • Redundant statements deleted / modified.
  • Bug fixed.

04/2023: In version PACIAE 2.3

  • In "main_23.f", a energy-dependent "smtj3" was introduced for low-energy loop-A at $\sqrt{s_{NN}} < 3$ GeV.
  • In "parini_23.f", the loop-A was improved.
  • In "parini_23.f", the hadron-hadron cross-sections inspired by additive quark model (AQM) wee introduced.
  • In all of ".f" code files, many "DO-ENDDO" statements and some statements of variable initialization were optimized to save running time (memory optimization).
  • In "main_23.f", different code blocks were separated by dashed lines "-----" with corresponding comments to explain their functions, i.e. "what we are doing here".
  • In "main_23.f", three subroutines "rest_hadronization", "rest_sfm" and "rest_coal" were introduced to hadronize the rest partons that failed during the normal hadronization process. The re-hadronization is set after normal hadronization now, instead of calling "pa2evnt" after the whole simulation as before. "pa2evnt" was also fixed and improved again.
  • In "main_23.f" and "parini_23.f" (near "PYINIT" calling), the re-generation of the NN binary collision were introduced to deal with case where charge or 4-momentum was not conserved or any errors occurred after calling PYTHIA.
  • In "main_23.f", the Coalescence Mode 2 was introduced. In this mode, the gluon splitting and energetic quark deexcitation would be performed before parton rescattering. It means that there will be no gluons into parton rescattering.
  • In "main_23.f", a subroutine "prt_final_info" was introduced to print final 4-momentum information to "main.out".
  • In subroutines "decmom_sbe" of "main_23.f" and "decmom" of "parini_23.f", minor bugs were fixed.
  • The subroutine "oscar" of "main_23.f" was improved to print event particle history correctly.
  • The spectator nucleons were moved after parton rescattering and hadron rescattering correctly avoiding the interactions between them and other hadrons, i.e. "sbh moving" blocs in "main_23.f".
  • The input parameter adj1(29) could select fragmentation function for both string fragmentation model (SFM) and deexcitation of coalescence model (Coal) now. A series of parameters of fragmentation functions were introduced.
  • In subroutine "scat" of "parini_23.f", gamma66 removing block was moved to do the job correctly.
  • In subroutine "ctlcre_par" of "parcas_23.f", minor bug were fixed.
  • In "coales_23.f", the coalescence model was improved.
  • In subroutine "tabh" of "coales_23.f", the specie of hadron is extended to 200 now.
  • In subroutine "coal" of "coales_23.f", the reconstruction of parton list after normal "coal" calling was fixed and improved. A process of "final coalescence try" for the quarks that failed in normal coalescence was introduced.
  • "usu.dat" updated.
  • "PACIAE.sh" updated.
  • Redundant statements deleted / modified.
  • Bug fixed.

03/2023: In version PACIAE 2.3

  • The "stahad_23.f" file has beed removed.

  • The simulation mode in PACIAE has been classified into three modes, controlled by switch "mstj2":

    • =1, (low-energy) pure hadron simulation loop-A;
    • =2, PYTHIA-like pure hadron simulation loop-B;
    • =3, parton-hadron cascade simulation loop-C.

  • In "main_23.f, parini_23.f", the low-energy simulation loop-A has been improved considering endothermic, exothermic processes and Delta particle decay more carefully. A COMMON BLOCK /delt/ has been added for processing Delta particle in loop-A.

  • Parameter "mstj3" was renamed as "smtj3" and "smtj3/10" was set as the probability of the Delta particle decay. Parameter "dtt" is the ratio of the inelastic cross-section to the total one for hadron scattering now (D=0.85 for high energy, 0.1 for low-energy loop-A).

  • In "main_23.f", two subroutines "share_p_PYJETS" and "share_p_sbe" are introduced to share the 4-momentum for the conservation adjustment. The array "throe_p" of COMMON BLOCK /sa16/ is used to store the lost 4-momentum now.

  • In "main_23.f", the subroutine "pa2ent" has been fixed.

  • In "remo" of "parini_23.f", the junctions were removed correctly now.

  • In "coales_23.f", the quark energy sorting, i.e. the subroutine "eord", was closed.

  • In "coales_23.f", a subroutine "PAPTDI" has been introduced for transverse momentum sampling. The subroutine "ffm" for high-energy deexcitation mechanism has been renamed and improved as "deexcitation" (two implementations "deexcitation_E" and "deexcitation_EP" for energy mode and light-cone variable mode).

  • The subroutine "break_glu" of "coales_23.f" has been fixed to break-up gluons correctly.

  • The subroutine "funcz" of "coales_23.f" was improved.

  • The COMMON BLOCK /sa18/ was re-set as switches and parameters about coalescence.

  • The "analy.f" has been improved for output the partial phase-space statistics. The forth distribution is transverse mass distribution now.

  • In "parini_23.f", "parcas_23.f", "coales_23.f", and "hadcas_23.f", all of the "TAB character" were replaced by corresponding "space character" safely.

  • "usu.dat" updated.

  • "PACIAE.sh" updated.

  • Redundant statements deleted / modified.

  • Bug fixed.

02/2023: In version PACIAE 2.3

  • In "main_23.f, analy.f", the parameter "parp22" was used to select y/eta in partial phase-space statistics.
  • In "main_23.f", The mistake-proofing statements were added.
  • In subroutine "oscar" of "main_23.f", the OSCAR1997A/1999A/2013A output was rewritten. 1997A for final particle information, 1999A for full event history (initial nucleon, initial parton state, final partonic state, initial hadron state and final particle state), and 2013A dummy up to now. Some related statements were added in "main_23.f" and "parini_23.f".
  • In "parini_23.f", some pre-statements for statistical hadronization was added.
  • In "coal_23.f", the separate position/momentum phase-space constraint was introduced.
  • In "p_23.f", the size of COMMON BLOCK /HEPEVT/ was extended to 80000, also with local /PYCBLS/. In subroutine "PYTIME", DATE_AND_TIME was activated. (However, a potential bug is MSTU(5). Still 10000 in p_23.f now. A possible way may be using the compiling option "-fdefault-integer-8" when one uses GFortran.) The modification of PYTHIA6 were indicated at the beginning of the file.
  • In "usu.dat", comments are updated and default values are given in (D=) (not yet tuned). In addition, the KF code list was attached at the end of the file directly, which was printed from PYTHIA6.
  • A file "stahad_23.f" was introduced for the statistical hadronization in the future. Dummy up to now.
  • "PACIAE.sh" updated.
  • Redundant statements deleted.
  • Bug fixed
  • In "main_23.f", "sfm_23.f", "stahad_23.f", "analy.f", and "p_23.f", all of the "TAB character" were replaced by corresponding "space character" safely. Partly done in "parini_23.f", "parcas_23.f", and "hadcas_23.f"

01/2023: In version PACIAE 2.3

  • In "main_23.f" and "parini_23.f", the low-energy simulation was introduced.
  • "PACIAE.sh" updated.
  • Bug fixed.

12/2022: In version PACIAE 2.3

  • In "parcas_23.f", the working COMMON BLOCK /parlist/ was replaced by /PYJETS/.
  • Output display optimized for "rms.out".
  • Redundant statements deleted.
  • The "LICENSE" file was introduced. (GPL v2.0 based on PYTHIA 8)
  • Bug fixed.

11/2022: In version PACIAE 2.3

  • In "parini_23.f", the selections of distributing nucleons into the overlapping region forcedly was introduced, controlled by adj(30).
  • In "parcas_23.f", all of six quarks and gluons were included now. (u, d, s, g old)
  • "eps09.f" has been removed now.
  • Redundant statements deleted.
  • Bug fixed.

10/2022: In version PACIAE 2.3

  • In subroutine "main" of "main_23.f", the real-time clock random seed was introduced.
  • In subroutine "scat" of "parini_23.f", the long-written statement about executing the binary collision by calling PYEVNW / PYEVNT was replaced by a new subroutine "xevent".
  • In "parini_23.f", the CME was introduced in PACIAE 2.3 based on Zhi-Lei's improvement in PACIAE 2.2.1b and PACIAE 2.2.1c. The statements about calculation of the eccentricity were corrected.
  • In "main_23.f", "eps09.f" and "p_23.f", the statements with old syntax were rewritten by more modern-style one.
  • A shell-script file "PACIAE.sh" was introduced for automatic compilation, building and running (pseudo-parallel running possible).
  • In "usu.dat", the gluon-splitting possibility were introduced for "coal_23.f".
  • The "README.md" file was introduce.
  • Bug fixed.

04/2021: In version PACIAE 2.2.1 b and c

  • The CME was introduced by Zhi-Lei She et al.

01/2020: In version PACIAE 2.3

  • The EPS09 nuclear shadowing was introduced by Liang, whose subroutine is called "shanul_eps09". An extra file named "eps09.f" was added.

About

A parton and hadron cascade model PACIAE (Parton And-hadron China Institute of Atomic Energy).

Topics

monte-carlo-simulation high-energy-physics heavy-ion-collisions nuclear-physics

Resources

Readme

License

GPL-2.0 license
Activity

Stars

4 stars

Watchers

1 watching

Forks

1 fork
Report repository

Releases 8

PACIAE 3.0.01 Latest
Aug 1, 2024
+ 7 releases

Packages

No packages published

Languages