100M events. UrQMD (eos 1) + Geant-4 based general-purpose simulation project for minbias (b = 0-15 fm) Xe(54/124) + W(74/184) (fixed target) collisions at T = 2.5 GeV/n. The wire-target is installed at z = -85 cm, wire diameter is 100 um, wire is shifted by 1 cm upwards in Y.
Please find below the request details:
inputfile.txt (371 Bytes) - v.2
convertUrQMD.C (10.2 KB) - v.2
MpdUnigenGenerator.cxx (10.3 KB)
MpdUnigenGenerator.h (2.6 KB)
runMC.C (11.6 KB) - v.2
runReco.C (8.7 KB)
Please use -‘global_polarization’ branch of the MpdRoot and upload files MpdUnigenGenerator.cxx and MpdUnigenGenerator.h into mpdroot/simulation/generators/mpdGen/
Configuration:
-
inputfile - UrQMD input file that simulates minbias (0-15 fm) Xe+W(T=2.5 GeV/n) collisions, sets resonances (rho, Omega, phi, K*, K0, Lambda*, Sigma*) as stable particles to be later decayed by Geant-4. Important №1: The value of ‘rsd’ parameter should be set as a unique integer for each UrQMD run. Important №2: Output data format should be set as f20 (#f20 in inputfile)
-
The output file of the UrQMD should be converted to mcini with convertUrQMD.C:
root -l -b -q convertUrQMD.C’(TString inputFileName, TString outputFileName, float tFinal = -1.)’ -
Use afterburner (GitHub - Troshinvv/mpd_afterburner_stable) to parameterize lambdas with collective flow and global polarization. Important: the last argument(4th) should be 0.
Usage:
3.1 Run generate_dicts.c to generate temporary dictionaries.
root -l -q generate_dicts.c
3.2 Upload in root session these dictionaries and run read_unigen_root.cpp
root -l -b
gSystem->Load(“AutoDict_vector_TVector3__cxx.so”)
gSystem->Load(“AutoDict_vector_UParticle__cxx.so”)
gSystem->Load(“AutoDict_std__vector_ROOT__Math__XYZVector__cxx.so”)
gSystem->Load(“AutoDict_vector_TVector3__cxx.so”)
gSystem->Load(“AutoDict_vector_UParticle__cxx.so”)
gSystem->Load(“AutoDict_std__vector_ROOT__Math__XYZVector__cxx.so”)
.L $read_unigen_root.cpp
read_unigen_root(TString inputFile, TString outputFile, “lamb_config.root”, 0)
.q -
runMC.C file was modified for global polarization transfer. Important: runMC.C should use UniGen generator for reading output from afterburner.
-
reco.C - full detector reconstruction, including the ECAL.
For simulation, please use the attached runMC.C macro.
For reconstruction please use the attached reco.C macro. It produces standard DST files.
To validate the production configuration, we need ~ 1,000,000 events reconstructed for QA.
The output files will have to be copied and stored at NICA cluster.
The ‘mpdroot’ should be tagged for this production.