Skip to content

Commit 6e01146

Browse files
hscheidhscheid
authored
add param for dielectron OO cocktail based on pythia (#2124)
* add param for OO based on pythia * add ini for OO LMee cocktail * Add test for OO cocktail * BR with max 6 digits after decimal point
1 parent f535c06 commit 6e01146

File tree

4 files changed

+134
-54
lines changed

4 files changed

+134
-54
lines changed

‎MC/config/PWGEM/decaytables/decaytable_LMee.dat‎

Lines changed: 54 additions & 54 deletions
Original file line numberDiff line numberDiff line change
@@ -280,74 +280,74 @@
280280
1 0 0.000000 -11 11 -11 11 0
281281
1 0 0.000000 -11 11 0 0 0
282282
221 eta 0 0 0 0.54786 0.001308 0.00000 1.50886E-07 0 1
283-
1 0 0.3935878 22 22 0 0 0
284-
1 0 0.3256175 111 111 111 0 0
285-
1 0 0.0002554 111 22 22 0 0
286-
1 0 0.2301877 211 -211 111 0 0
287-
1 0 0.0427689 22 211 -211 0 0
288-
1 2 0.0070044 22 11 -11 0 0
289-
1 0 0.0003101 22 13 -13 0 0
290-
1 0 0.0000058 -13 13 0 0 0
291-
1 0 0.0000007 -11 11 0 0 0
292-
1 0 0.0000000 -11 11 -11 11 0
293-
1 0 0.0000000 211 -211 11 -11 0
283+
1 0 0.393588 22 22 0 0 0
284+
1 0 0.325618 111 111 111 0 0
285+
1 0 0.000255 111 22 22 0 0
286+
1 0 0.230188 211 -211 111 0 0
287+
1 0 0.042769 22 211 -211 0 0
288+
1 2 0.007004 22 11 -11 0 0
289+
1 0 0.000310 22 13 -13 0 0
290+
1 0 0.000006 -13 13 0 0 0
291+
1 0 0.000001 -11 11 0 0 0
292+
1 0 0.000000 -11 11 -11 11 0
293+
1 0 0.000000 211 -211 11 -11 0
294294
331 eta' 0 0 0 0.95778 0.000188 0.002 1.00336E-09 0 1
295-
1 0 0.4246864 211 -211 221 0 0
296-
1 0 0.2948085 22 113 0 0 0
297-
1 0 0.2237140 111 111 221 0 0
298-
1 0 0.0251900 22 223 0 0 0
299-
1 0 0.0001970 223 -11 11 0 0
300-
1 0 0.0230680 22 22 0 0 0
301-
1 0 0.0024994 111 111 111 0 0
302-
1 0 0.0001130 22 -13 13 0 0
303-
1 0 0.0036058 211 -211 111 0 0
304-
1 0 0.0000010 211 -211 211 -211 0
305-
1 0 0.0001800 211 -211 111 111 0
306-
1 0 0.0000000 211 -211 -11 11 0
307-
1 2 0.0004803 22 11 -11 0 0
308-
1 0 0.0000000 -11 11 -11 11 0
295+
1 0 0.424687 211 -211 221 0 0
296+
1 0 0.294809 22 113 0 0 0
297+
1 0 0.223714 111 111 221 0 0
298+
1 0 0.025190 22 223 0 0 0
299+
1 0 0.000197 223 -11 11 0 0
300+
1 0 0.023068 22 22 0 0 0
301+
1 0 0.002499 111 111 111 0 0
302+
1 0 0.000113 22 -13 13 0 0
303+
1 0 0.003606 211 -211 111 0 0
304+
1 0 0.000001 211 -211 211 -211 0
305+
1 0 0.000180 211 -211 111 111 0
306+
1 0 0.000000 211 -211 -11 11 0
307+
1 2 0.000480 22 11 -11 0 0
308+
1 0 0.000000 -11 11 -11 11 0
309309
113 rho0 0 0 0 0.77526 0.14739 0.40000 1.33000E-12 0 1
310310
1 3 0.988927 211 -211 0 0 0
311311
1 0 0.009900 211 -211 22 0 0
312-
1 0 0.0004669 111 22 0 0 0
313-
1 0 0.0002997 221 22 0 0 0
314-
1 0 0.0000445 111 111 22 0 0
315-
1 0 0.0000472 11 -11 0 0 0
316-
1 0 0.0000455 13 -13 0 0 0
312+
1 0 0.000467 111 22 0 0 0
313+
1 0 0.000300 221 22 0 0 0
314+
1 0 0.000045 111 111 22 0 0
315+
1 0 0.000047 11 -11 0 0 0
316+
1 0 0.000046 13 -13 0 0 0
317317
1 0 0.000101 211 -211 111 0 0
318318
1 0 0.000018 -211 211 -211 211 0
319319
1 0 0.000016 -211 211 111 111 0
320320
223 omega 0 0 0 0.78266 0.00868 0.10000 2.33000E-11 0 1
321-
1 1 0.8924576 211 -211 111 0 0
322-
1 0 0.0832642 22 111 0 0 0
323-
1 3 0.0153487 211 -211 0 0 0
324-
1 0 0.0004478 221 22 0 0 0
325-
1 0 0.0001336 111 13 -13 0 0
326-
1 0 0.0007673 111 11 -11 0 0
327-
1 0 0.0000741 11 -11 0 0 0
328-
1 0 0.0000666 111 111 22 0 0
329-
1 0 0.0000740 13 -13 0 0 0
321+
1 1 0.892458 211 -211 111 0 0
322+
1 0 0.083264 22 111 0 0 0
323+
1 3 0.015349 211 -211 0 0 0
324+
1 0 0.000448 221 22 0 0 0
325+
1 0 0.000134 111 13 -13 0 0
326+
1 0 0.000767 111 11 -11 0 0
327+
1 0 0.000074 11 -11 0 0 0
328+
1 0 0.000067 111 111 22 0 0
329+
1 0 0.000074 13 -13 0 0 0
330330
333 phi 0 0 0 1.01946 0.004249 0.01500 4.63000E-11 0 1
331-
1 3 0.4988842 321 -321 0 0 0
332-
1 3 0.3364070 130 310 0 0 0
331+
1 3 0.498884 321 -321 0 0 0
332+
1 3 0.336407 130 310 0 0 0
333333
1 0 0.042259 213 -211 0 0 0
334334
1 0 0.042259 -213 211 0 0 0
335335
1 0 0.042259 113 111 0 0 0
336336
1 1 0.026534 211 -211 111 0 0
337-
1 0 0.0130621 221 22 0 0 0
338-
1 0 0.0013280 111 22 0 0 0
339-
1 0 0.0002963 11 -11 0 0 0
340-
1 0 0.0002862 13 -13 0 0 0
341-
1 0 0.0001084 221 -11 11 0 0
342-
1 0 0.0000950 211 -211 0 0 0
343-
1 0 0.0000473 223 111 0 0 0
344-
1 0 0.0000410 -211 211 22 0 0
345-
1 0 0.0001126 111 111 22 0 0
346-
1 0 0.0000039 -211 211 -211 211 0
347-
1 0 0.0000133 111 -11 11 0 0
348-
1 0 0.0000726 111 221 22 0 0
337+
1 0 0.013062 221 22 0 0 0
338+
1 0 0.001328 111 22 0 0 0
339+
1 0 0.000296 11 -11 0 0 0
340+
1 0 0.000286 13 -13 0 0 0
341+
1 0 0.000108 221 -11 11 0 0
342+
1 0 0.000095 211 -211 0 0 0
343+
1 0 0.000047 223 111 0 0 0
344+
1 0 0.000041 -211 211 22 0 0
345+
1 0 0.000113 111 111 22 0 0
346+
1 0 0.000004 -211 211 -211 211 0
347+
1 0 0.000013 111 -11 11 0 0
348+
1 0 0.000073 111 221 22 0 0
349349
1 0 0.000063 331 22 0 0 0
350-
1 0 0.0000143 -13 13 22 0 0
350+
1 0 0.000014 -13 13 22 0 0
351351
443 J/Psi 0 0 0 3.09690 0.000093 0.00000 2.13000E-09 0 1
352352
1 12 0.881200 82 -82 0 0 0
353353
1 0 0.059710 11 -11 0 0 0

‎MC/config/PWGEM/ini/GeneratorEMCocktail_Run3_OO.ini‎

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
### The setup uses an external event generator
2+
### This part sets the path of the file and the function call to retrieve it
3+
4+
[GeneratorExternal]
5+
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/GeneratorEMCocktailV2.C
6+
funcName=GenerateEMCocktail(200,0,3,63,"${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/parametrizations/OO_536TeV_param.json","536TeV_MB_pythia",200,0.0,30.0,10000,1,1,0,0,"",0,1.1,"${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/decaytables/decaytable_LMee.dat",1)

‎MC/config/PWGEM/ini/tests/GeneratorEMCocktail_Run3_OO.C‎

Lines changed: 64 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,64 @@
1+
int External()
2+
{
3+
4+
int checkPdgDecay = -11;
5+
std::string path{"o2sim_Kine.root"};
6+
TFile file(path.c_str(), "READ");
7+
if (file.IsZombie()) {
8+
std::cerr << "Cannot open ROOT file " << path << "\n";
9+
return 1;
10+
}
11+
12+
auto tree = (TTree*)file.Get("o2sim");
13+
std::vector<o2::MCTrack>* tracks{};
14+
tree->SetBranchAddress("MCTrack", &tracks);
15+
16+
int nMesons{};
17+
int nMesonsDiElectronDecay{};
18+
auto nEvents = tree->GetEntries();
19+
20+
for (int i = 0; i < nEvents; i++) {
21+
tree->GetEntry(i);
22+
for (auto& track : *tracks) {
23+
auto pdg = track.GetPdgCode();
24+
auto y = track.GetRapidity();
25+
if ((pdg == 111) || (pdg == 221) || (pdg == 331) || (pdg == 223) || (pdg == 113) || (pdg == 333)) {
26+
if ((y>-1.2) && (y<1.2)) {
27+
nMesons++;
28+
Int_t counterel = 0;
29+
Int_t counterpos = 0;
30+
int k1 = track.getFirstDaughterTrackId();
31+
int k2 = track.getLastDaughterTrackId();
32+
// k1 < k2 and no -1 for k2
33+
for (int d=k1; d <= k2; d++) {
34+
if (d>0) {
35+
auto decay = (*tracks)[d];
36+
int pdgdecay = decay.GetPdgCode();
37+
if (pdgdecay == 11) {
38+
counterel++;
39+
}
40+
if (pdgdecay == -11) {
41+
counterpos++;
42+
}
43+
}
44+
}
45+
if ((counterel>0) && (counterpos>0)) nMesonsDiElectronDecay++;
46+
}
47+
}
48+
}
49+
}
50+
51+
std::cout << "#events: " << nEvents << "\n"
52+
<< "#mesons: " << nMesons << "\n"
53+
<< "#mesons which decay semi-electronicly: " << nMesonsDiElectronDecay << "\n";
54+
if (nMesonsDiElectronDecay < nEvents) {
55+
std::cerr << "One should have at least one meson that decays into dielectrons per event.\n";
56+
return 1;
57+
}
58+
if (nMesons < nEvents) {
59+
std::cerr << "One meson per event should be produced.\n";
60+
return 1;
61+
}
62+
63+
return 0;
64+
}

‎MC/config/PWGEM/parametrizations/OO_536TeV_param.json‎

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,10 @@
1+
{
2+
"536TeV_MB_pythia": {
3+
"111_pt": "56.6336*x*pow(1+(sqrt(x*x+0.138918*0.138918)-0.138918)/(6.8304*0.139501),-6.8304)",
4+
"221_pt": "1939.44*x*pow(1+(sqrt(x*x+-0.424653*-0.424653)--0.424653)/(6.84504*0.0745758),-6.84504)",
5+
"331_pt": "2.89402*x*pow(1+(sqrt(x*x+-0.322689*-0.322689)--0.322689)/(8.32826*0.211736),-8.32826)",
6+
"113_pt": "75.4035*x*pow(1+(sqrt(x*x+-0.360389*-0.360389)--0.360389)/(7.74615*0.16301),-7.74615)",
7+
"223_pt": "2.29*x*pow(1+(sqrt(x*x+0.405007*0.405007)-0.405007)/(7.61366*0.243978),-7.61366)",
8+
"333_pt": "0.646087*x*pow(1+(sqrt(x*x+-0.257169*-0.257169)--0.257169)/(8.17199*0.261382),-8.17199)"
9+
}
10+
}

0 commit comments

Comments
 (0)