Geant-4 based PWG4 request for Bi (83/209) +Bi (83/209) collisions at 9.2 GeV with the latest version of dE/dx calculations in the TPC and v.4 of the ECAL geometry.
Please find below the request details:
inputfile (444 Bytes)
runMC.C (10.9 KB)
reco.C (8.7 KB)
MpdDecayConfig.txt (357 Bytes)
MpdStack.cxx (21.5 KB)
This is a general-purpose project. Minbias events are generated by UrQMD, light vector meson BR for dielectron decays are enhanced by x20. Decays of short-lived resonances (to have realistic peak widths) and pi0/eta-meson are managed by Pythia8.
The simulation project consists of a few stages:
Events are generated by UrQMD. It should be run with an attached inputfile, which was set up to generate 500 minbias BiBi@9.2 collisions. Some of the particles are requested to be stable in the input file to be further decayed in the mpdroot. Important: The value of ‘rsd’ parameter should be set as a unique integer for each UrQMD run.
Make a fresh check-out of the -dev branch of the MpdRoot, replace:
simulation/mcStack/MpdStack.cxx
gconfig/MpdDecayConfig.txt
files with the attached versions and then compile.
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 ~ 200,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.
QA done:
/eos/nica/mpd/sim/data/exp/dst-BiBi-09.2GeV-mp04-22-500ev
I copied the data, but did not check whether everything was copied. The NIKA cluster has become unavailable. I’ll check data later when it’s turned on.
I’m not sure everything is correct with QA. This production worked too fast.
mpdroot git tag: mprequest24
Dear Andrey,
Sorry for the late response. There is something wrong with the QA DSTs. They are 5-10 times smaller than I would expect them to be and the impact parameter is always set to zero. I suspect that something could have gone wrong on the level of UrQMD event generation. Could you please have a look? Thanks.
With best regards, Victor.
Sorry, I was out of Dubna at the weekend. I wrote that the files are small and the tasks run very fast. It’s unusual.
I checked the input files, everything is the same, only the rsd parameter changes.
Just in case, I am attaching a couple of input files directly in the text. Looks basically like in all productions.
first:
pro 209 83
tar 209 83
nev 500
imp -15
ecm 9.2
tim 200 200
eos 0
rsd 1009
#do not decay pi
stb 101
#do not decay eta
stb 102
#do not decay omega
stb 103
#do not decay rho
stb 104
#do not decay eta-prime
stb 107
#do not decay phi
stb 109
#do not decay Delta
#stb 17
#stb -17
f13
#f14
f15
f16
f19
f20
xxx
second:
pro 209 83
tar 209 83
nev 500
imp -15
ecm 9.2
tim 200 200
eos 0
rsd 30992
#do not decay pi
stb 101
#do not decay eta
stb 102
#do not decay omega
stb 103
#do not decay rho
stb 104
#do not decay eta-prime
stb 107
#do not decay phi
stb 109
#do not decay Delta
#stb 17
#stb -17
f13
#f14
f15
f16
f19
f20
xxx
In what there can be a problem I do not understand.
Dear Andrey,
Thanks for update. I do not see problems with the input files. Can you download somewhere all files for one job: UrQMD, simulated file, DST. The I will try to localize a problem. Thanks. Victor.
Good afternoon.
A little delayed, the Dirac root certificate expired. I had to wait until admins make and install the new one. All files in:
/eos/nica/mpd/sim/qa
.f14 - urqmd file
.dst.root - simulation
.reco.root reconstruction
I also added the task log file: 1674438_OutputSandbox
Andrey, could you please give me permision to read the file?
Done. Again, the notification from the forum got into spam, the whitelist did not work, I’ll try to solve the problem.
Dear Andrey,
I found a problem. You generate a good UrQMD output file but then for some reason in runMC.C file you do not read it but use the BOX generator instead. So the lines below confirm that runMC.C processes BOX-generated events with 100 particles each:
==================================================================
Run 0 start.
[INFO] FairBoxGenerator: Event, vertex = (0,0,0) cm, multiplicity 100
Event 0
$$$X Primary track 90
$$$X Primary track 80
$$$X Primary track 70
$$$X Primary track 60
$$$X Primary track 50
$$$X Primary track 40
$$$X Primary track 30
$$$X Primary track 20
$$$X Primary track 10
-I- MpdStack: Filling MCTrack array…
-I- MpdStack: Number of primaries = 100
Total number of particles = 46579
Number of tracks in output = 411
-I- MpdStack: Updating track indizes…stack and 6 collections updated.
-I- MpdFfd: 6 points registered in this event.
-I- MpdTpc: 7932 points registered in this event.
User CPU time: 78.4362 Memory: resident 492828, virtual 907004
-I- MpdTof: 246 points registered in this event.
-I- MpdEmcKI: 323 points registered in this event.
-I- MpdZdc: 194 points registered in this event.
[INFO] FairBoxGenerator: Event, vertex = (0,0,0) cm, multiplicity 100
Event 1
$$$X Primary track 90
$$$X Primary track 80
$$$X Primary track 70
$$$X Primary track 60
$$$X Primary track 50
$$$X Primary track 40
$$$X Primary track 30
$$$X Primary track 20
$$$X Primary track 10
-I- MpdStack: Filling MCTrack array…
-I- MpdStack: Number of primaries = 100
Total number of particles = 48513
Number of tracks in output = 559
-I- MpdStack: Updating track indizes…stack and 6 collections updated.
-I- MpdFfd: 7 points registered in this event.
-I- MpdTpc: 8181 points registered in this event.
User CPU time: 80.8105 Memory: resident 497300, virtual 911444
-I- MpdTof: 255 points registered in this event.
-I- MpdEmcKI: 374 points registered in this event.
-I- MpdZdc: 382 points registered in this event.
[INFO] FairBoxGenerator: Event, vertex = (0,0,0) cm, multiplicity 100
======
while these should be UrMQD events like the ones shown below:
URQMD VERSION USED = 30400
-I MpdUrqmdGenerator: Event 2, b = 10.9499998 fm, multiplicity 572, ekin: 43.2400017
Event 1
$$$X Primary track 670
$$$X Primary track 660
$$$X Primary track 650
$$$X Primary track 640
$$$X Primary track 630
$$$X Primary track 620
$$$X Primary track 50
$$$X Primary track 40
$$$X Primary track 30
$$$X Primary track 20
$$$X Primary track 10
-I- MpdStack: Filling MCTrack array…
-I- MpdStack: Number of primaries = 726
Total number of particles = 1370558
Number of tracks in output = 12284
-I- MpdStack: Updating track indizes…stack and 6 collections updated.
-I- MpdFfd: 318 points registered in this event.
-I- MpdTpc: 36660 points registered in this event.
User CPU time: 151.6519012 Memory: resident 1053060, virtual 1419824
-I- MpdTof: 821 points registered in this event.
-I- MpdEmcKI: 12407 points registered in this event.
-I- MpdZdc: 3682 points registered in this event.
So I guess you somehow used a wrong runMC.C file. The one attached to the first message in the thread requests to use UrQMD events:
#define URQMD // Choose generator: URQMD VHLLE FLUID PART ION BOX HSD LAQGSM SMASH
I also run it on UrQMD file and it works OK. So it should be something wrong in your configuration.
With best regards, Victor.
P.S. I also get all notifications from mpdforum in SPAM (I use gmail). I do not know why and how to fix the problem.
Sorry, this is my mistake. In the task template, the basic script from the mpdroot distribution remained, and not the script attached to the production request. Launched qa again.
At last, the problems of the cluster were resolved and the administrators gave access to eos. Posted a new qa in:
/eos/nica/mpd/sim/data/exp/dst-BiBi-09.2GeV-mp05-22-500ev/BiBi/09.0GeV-mb/urqmd/BiBi-09.2GeV-mp05-22-500ev
Dear Andrey,
Thanks, the QA is fine.
I would like to ask you to replace inputfile and runMC.C files with the attached ones and then start a mass production:
inputfile (306 Bytes)
runMC.C (10.9 KB)
Thanks.
With best regards, Victor.
Production completed. It took more time, since SPD production was working at the same time, I had to share the resources of TIER clusters with them.
Files in:
/eos/nica/mpd/sim/data/exp/dst-BiBi-09.2GeV-mp05-22-500ev/BiBi/09.0GeV-mb/urqmd/BiBi-09.2GeV-mp05-22v2-500ev/
While the NIСA cluster was not working, and I couldn’t copy the results of QA sample to the NIСA cluster,
I made a production with the first version of the runMC.C and the inputfile(resources for production have already been allocated, they shouldn’t be idle),
Files in:
/eos/nica/mpd/sim/data/exp/dst-BiBi-09.2GeV-mp05-22-500ev/BiBi/09.0GeV-mb/urqmd/BiBi-09.2GeV-mp05-22-500ev/
Maybe these files will also be useful for analysis