Skip to content

add param for dielectron OO cocktail based on pythia #2124

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
dsekihat merged 4 commits into from
Sep 24, 2025
Merged

add param for dielectron OO cocktail based on pythia #2124

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
e12274c
add param for OO based on pythia
hscheid Sep 19, 2025
b532e0f
add ini for OO LMee cocktail
hscheid Sep 19, 2025
8fc1ebb
Add test for OO cocktail
hscheid Sep 20, 2025
22b1fb7
BR with max 6 digits after decimal point
hscheid Sep 23, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
108 changes: 54 additions & 54 deletions MC/config/PWGEM/decaytables/decaytable_LMee.dat
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -280,74 +280,74 @@
1 0 0.000000 -11 11 -11 11 0
1 0 0.000000 -11 11 0 0 0
221 eta 0 0 0 0.54786 0.001308 0.00000 1.50886E-07 0 1
1 0 0.3935878 22 22 0 0 0
1 0 0.3256175 111 111 111 0 0
1 0 0.0002554 111 22 22 0 0
1 0 0.2301877 211 -211 111 0 0
1 0 0.0427689 22 211 -211 0 0
1 2 0.0070044 22 11 -11 0 0
1 0 0.0003101 22 13 -13 0 0
1 0 0.0000058 -13 13 0 0 0
1 0 0.0000007 -11 11 0 0 0
1 0 0.0000000 -11 11 -11 11 0
1 0 0.0000000 211 -211 11 -11 0
1 0 0.393588 22 22 0 0 0
1 0 0.325618 111 111 111 0 0
1 0 0.000255 111 22 22 0 0
1 0 0.230188 211 -211 111 0 0
1 0 0.042769 22 211 -211 0 0
1 2 0.007004 22 11 -11 0 0
1 0 0.000310 22 13 -13 0 0
1 0 0.000006 -13 13 0 0 0
1 0 0.000001 -11 11 0 0 0
1 0 0.000000 -11 11 -11 11 0
1 0 0.000000 211 -211 11 -11 0
331 eta' 0 0 0 0.95778 0.000188 0.002 1.00336E-09 0 1
1 0 0.4246864 211 -211 221 0 0
1 0 0.2948085 22 113 0 0 0
1 0 0.2237140 111 111 221 0 0
1 0 0.0251900 22 223 0 0 0
1 0 0.0001970 223 -11 11 0 0
1 0 0.0230680 22 22 0 0 0
1 0 0.0024994 111 111 111 0 0
1 0 0.0001130 22 -13 13 0 0
1 0 0.0036058 211 -211 111 0 0
1 0 0.0000010 211 -211 211 -211 0
1 0 0.0001800 211 -211 111 111 0
1 0 0.0000000 211 -211 -11 11 0
1 2 0.0004803 22 11 -11 0 0
1 0 0.0000000 -11 11 -11 11 0
1 0 0.424687 211 -211 221 0 0
1 0 0.294809 22 113 0 0 0
1 0 0.223714 111 111 221 0 0
1 0 0.025190 22 223 0 0 0
1 0 0.000197 223 -11 11 0 0
1 0 0.023068 22 22 0 0 0
1 0 0.002499 111 111 111 0 0
1 0 0.000113 22 -13 13 0 0
1 0 0.003606 211 -211 111 0 0
1 0 0.000001 211 -211 211 -211 0
1 0 0.000180 211 -211 111 111 0
1 0 0.000000 211 -211 -11 11 0
1 2 0.000480 22 11 -11 0 0
1 0 0.000000 -11 11 -11 11 0
113 rho0 0 0 0 0.77526 0.14739 0.40000 1.33000E-12 0 1
1 3 0.988927 211 -211 0 0 0
1 0 0.009900 211 -211 22 0 0
1 0 0.0004669 111 22 0 0 0
1 0 0.0002997 221 22 0 0 0
1 0 0.0000445 111 111 22 0 0
1 0 0.0000472 11 -11 0 0 0
1 0 0.0000455 13 -13 0 0 0
1 0 0.000467 111 22 0 0 0
1 0 0.000300 221 22 0 0 0
1 0 0.000045 111 111 22 0 0
1 0 0.000047 11 -11 0 0 0
1 0 0.000046 13 -13 0 0 0
1 0 0.000101 211 -211 111 0 0
1 0 0.000018 -211 211 -211 211 0
1 0 0.000016 -211 211 111 111 0
223 omega 0 0 0 0.78266 0.00868 0.10000 2.33000E-11 0 1
1 1 0.8924576 211 -211 111 0 0
1 0 0.0832642 22 111 0 0 0
1 3 0.0153487 211 -211 0 0 0
1 0 0.0004478 221 22 0 0 0
1 0 0.0001336 111 13 -13 0 0
1 0 0.0007673 111 11 -11 0 0
1 0 0.0000741 11 -11 0 0 0
1 0 0.0000666 111 111 22 0 0
1 0 0.0000740 13 -13 0 0 0
1 1 0.892458 211 -211 111 0 0
1 0 0.083264 22 111 0 0 0
1 3 0.015349 211 -211 0 0 0
1 0 0.000448 221 22 0 0 0
1 0 0.000134 111 13 -13 0 0
1 0 0.000767 111 11 -11 0 0
1 0 0.000074 11 -11 0 0 0
1 0 0.000067 111 111 22 0 0
1 0 0.000074 13 -13 0 0 0
333 phi 0 0 0 1.01946 0.004249 0.01500 4.63000E-11 0 1
1 3 0.4988842 321 -321 0 0 0
1 3 0.3364070 130 310 0 0 0
1 3 0.498884 321 -321 0 0 0
1 3 0.336407 130 310 0 0 0
1 0 0.042259 213 -211 0 0 0
1 0 0.042259 -213 211 0 0 0
1 0 0.042259 113 111 0 0 0
1 1 0.026534 211 -211 111 0 0
1 0 0.0130621 221 22 0 0 0
1 0 0.0013280 111 22 0 0 0
1 0 0.0002963 11 -11 0 0 0
1 0 0.0002862 13 -13 0 0 0
1 0 0.0001084 221 -11 11 0 0
1 0 0.0000950 211 -211 0 0 0
1 0 0.0000473 223 111 0 0 0
1 0 0.0000410 -211 211 22 0 0
1 0 0.0001126 111 111 22 0 0
1 0 0.0000039 -211 211 -211 211 0
1 0 0.0000133 111 -11 11 0 0
1 0 0.0000726 111 221 22 0 0
1 0 0.013062 221 22 0 0 0
1 0 0.001328 111 22 0 0 0
1 0 0.000296 11 -11 0 0 0
1 0 0.000286 13 -13 0 0 0
1 0 0.000108 221 -11 11 0 0
1 0 0.000095 211 -211 0 0 0
1 0 0.000047 223 111 0 0 0
1 0 0.000041 -211 211 22 0 0
1 0 0.000113 111 111 22 0 0
1 0 0.000004 -211 211 -211 211 0
1 0 0.000013 111 -11 11 0 0
1 0 0.000073 111 221 22 0 0
1 0 0.000063 331 22 0 0 0
1 0 0.0000143 -13 13 22 0 0
1 0 0.000014 -13 13 22 0 0
443 J/Psi 0 0 0 3.09690 0.000093 0.00000 2.13000E-09 0 1
1 12 0.881200 82 -82 0 0 0
1 0 0.059710 11 -11 0 0 0
Expand Down
6 changes: 6 additions & 0 deletions MC/config/PWGEM/ini/GeneratorEMCocktail_Run3_OO.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,6 @@
### The setup uses an external event generator
### This part sets the path of the file and the function call to retrieve it

[GeneratorExternal]
fileName = ${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGEM/external/generator/GeneratorEMCocktailV2.C
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)
64 changes: 64 additions & 0 deletions MC/config/PWGEM/ini/tests/GeneratorEMCocktail_Run3_OO.C
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,64 @@
int External()
{

int checkPdgDecay = -11;
std::string path{"o2sim_Kine.root"};
TFile file(path.c_str(), "READ");
if (file.IsZombie()) {
std::cerr << "Cannot open ROOT file " << path << "\n";
return 1;
}

auto tree = (TTree*)file.Get("o2sim");
std::vector<o2::MCTrack>* tracks{};
tree->SetBranchAddress("MCTrack", &tracks);

int nMesons{};
int nMesonsDiElectronDecay{};
auto nEvents = tree->GetEntries();

for (int i = 0; i < nEvents; i++) {
tree->GetEntry(i);
for (auto& track : *tracks) {
auto pdg = track.GetPdgCode();
auto y = track.GetRapidity();
if ((pdg == 111) || (pdg == 221) || (pdg == 331) || (pdg == 223) || (pdg == 113) || (pdg == 333)) {
if ((y>-1.2) && (y<1.2)) {
nMesons++;
Int_t counterel = 0;
Int_t counterpos = 0;
int k1 = track.getFirstDaughterTrackId();
int k2 = track.getLastDaughterTrackId();
// k1 < k2 and no -1 for k2
for (int d=k1; d <= k2; d++) {
if (d>0) {
auto decay = (*tracks)[d];
int pdgdecay = decay.GetPdgCode();
if (pdgdecay == 11) {
counterel++;
}
if (pdgdecay == -11) {
counterpos++;
}
}
}
if ((counterel>0) && (counterpos>0)) nMesonsDiElectronDecay++;
}
}
}
}

std::cout << "#events: " << nEvents << "\n"
<< "#mesons: " << nMesons << "\n"
<< "#mesons which decay semi-electronicly: " << nMesonsDiElectronDecay << "\n";
if (nMesonsDiElectronDecay < nEvents) {
std::cerr << "One should have at least one meson that decays into dielectrons per event.\n";
return 1;
}
if (nMesons < nEvents) {
std::cerr << "One meson per event should be produced.\n";
return 1;
}

return 0;
}
10 changes: 10 additions & 0 deletions MC/config/PWGEM/parametrizations/OO_536TeV_param.json
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,10 @@
{
"536TeV_MB_pythia": {
"111_pt": "56.6336*x*pow(1+(sqrt(x*x+0.138918*0.138918)-0.138918)/(6.8304*0.139501),-6.8304)",
"221_pt": "1939.44*x*pow(1+(sqrt(x*x+-0.424653*-0.424653)--0.424653)/(6.84504*0.0745758),-6.84504)",
"331_pt": "2.89402*x*pow(1+(sqrt(x*x+-0.322689*-0.322689)--0.322689)/(8.32826*0.211736),-8.32826)",
"113_pt": "75.4035*x*pow(1+(sqrt(x*x+-0.360389*-0.360389)--0.360389)/(7.74615*0.16301),-7.74615)",
"223_pt": "2.29*x*pow(1+(sqrt(x*x+0.405007*0.405007)-0.405007)/(7.61366*0.243978),-7.61366)",
"333_pt": "0.646087*x*pow(1+(sqrt(x*x+-0.257169*-0.257169)--0.257169)/(8.17199*0.261382),-8.17199)"
}
}