diff --git a/MC/config/PWGEM/decaytables/decaytable_LMee.dat b/MC/config/PWGEM/decaytables/decaytable_LMee.dat index 5b3953b6f..03f88f65a 100644 --- a/MC/config/PWGEM/decaytables/decaytable_LMee.dat +++ b/MC/config/PWGEM/decaytables/decaytable_LMee.dat @@ -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 diff --git a/MC/config/PWGEM/ini/GeneratorEMCocktail_Run3_OO.ini b/MC/config/PWGEM/ini/GeneratorEMCocktail_Run3_OO.ini new file mode 100644 index 000000000..83898beae --- /dev/null +++ b/MC/config/PWGEM/ini/GeneratorEMCocktail_Run3_OO.ini @@ -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) diff --git a/MC/config/PWGEM/ini/tests/GeneratorEMCocktail_Run3_OO.C b/MC/config/PWGEM/ini/tests/GeneratorEMCocktail_Run3_OO.C new file mode 100644 index 000000000..a78a0cc8e --- /dev/null +++ b/MC/config/PWGEM/ini/tests/GeneratorEMCocktail_Run3_OO.C @@ -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* 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; +} diff --git a/MC/config/PWGEM/parametrizations/OO_536TeV_param.json b/MC/config/PWGEM/parametrizations/OO_536TeV_param.json new file mode 100644 index 000000000..bad021d49 --- /dev/null +++ b/MC/config/PWGEM/parametrizations/OO_536TeV_param.json @@ -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)" + } +} \ No newline at end of file