From 97e1ffd933ef26a6ed26a2d9e74c5c7c68e20db2 Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Thu, 9 Jan 2025 19:12:05 +0100 Subject: [PATCH 1/2] PWGHF: rename config of OmegaC->Omega in OmegaC NoDecay --- ...mega.ini => GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.ini} | 2 +- ...to_Omega.C => GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.C} | 0 ...c_to_Omega.cfg => pythia8_charmhadronic_OmegaC_NoDecay.cfg} | 3 --- 3 files changed, 1 insertion(+), 4 deletions(-) rename MC/config/PWGHF/ini/{GeneratorHF_D2H_bbbar_Mode2_OmegaC_to_Omega.ini => GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.ini} (89%) rename MC/config/PWGHF/ini/tests/{GeneratorHF_D2H_bbbar_Mode2_OmegaC_to_Omega.C => GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.C} (100%) rename MC/config/PWGHF/pythia8/generator/{pythia8_charmhadronic_Omegac_to_Omega.cfg => pythia8_charmhadronic_OmegaC_NoDecay.cfg} (95%) diff --git a/MC/config/PWGHF/ini/GeneratorHF_D2H_bbbar_Mode2_OmegaC_to_Omega.ini b/MC/config/PWGHF/ini/GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.ini similarity index 89% rename from MC/config/PWGHF/ini/GeneratorHF_D2H_bbbar_Mode2_OmegaC_to_Omega.ini rename to MC/config/PWGHF/ini/GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.ini index 60af2005c..2661ee333 100644 --- a/MC/config/PWGHF/ini/GeneratorHF_D2H_bbbar_Mode2_OmegaC_to_Omega.ini +++ b/MC/config/PWGHF/ini/GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.ini @@ -4,5 +4,5 @@ fileName=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/external/generator/generator_py funcName=GeneratorPythia8GapTriggeredBeauty(1, -1.5, 1.5, -1.5, 1.5, {4332}) [GeneratorPythia8] -config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_Omegac_to_Omega.cfg +config=${O2DPG_MC_CONFIG_ROOT}/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_OmegaC_NoDecay.cfg includePartonEvent=true diff --git a/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Mode2_OmegaC_to_Omega.C b/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.C similarity index 100% rename from MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Mode2_OmegaC_to_Omega.C rename to MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.C diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_Omegac_to_Omega.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_OmegaC_NoDecay.cfg similarity index 95% rename from MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_Omegac_to_Omega.cfg rename to MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_OmegaC_NoDecay.cfg index d34c47a34..20d4c2fb6 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_Omegac_to_Omega.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_charmhadronic_OmegaC_NoDecay.cfg @@ -40,6 +40,3 @@ BeamRemnants:saturation 5 ### switch off OmegaC decays 4332:onMode = off - -### Omega_c -> Omega + X -4332:onIfAny = 3334 From 556e374b533ffb32ebe745d38701383dc6b598f9 Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Thu, 9 Jan 2025 19:17:51 +0100 Subject: [PATCH 2/2] Modify test accordingly --- ...GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.C | 17 +++++++---------- 1 file changed, 7 insertions(+), 10 deletions(-) diff --git a/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.C b/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.C index 31f001e08..73788a608 100644 --- a/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.C +++ b/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Mode2_OmegaC_NoDecay.C @@ -4,7 +4,7 @@ int External() { int checkPdgQuarkOne = 5; int checkPdgHadron{4332}; - int checkHadronDecays{3334}; + int checkHadronDecays{0}; TFile file(path.c_str(), "READ"); if (file.IsZombie()) { @@ -19,7 +19,7 @@ int External() { tree->SetBranchAddress("MCEventHeader.", &eventHeader); int nEventsInj{}; - int nQuarks{}, nSignals{}, nSignalGoodDecay{}; + int nQuarks{}, nSignals{}; auto nEvents = tree->GetEntries(); for (int i = 0; i < nEvents; i++) { @@ -44,10 +44,8 @@ int External() { nSignals++; // count signal PDG for (int j{track.getFirstDaughterTrackId()}; j <= track.getLastDaughterTrackId(); ++j) { - auto pdgDau = tracks->at(j).GetPdgCode(); - if (std::abs(pdgDau) == checkHadronDecays) { - nSignalGoodDecay; - break; + if (j >= 0) { + checkHadronDecays += 1; } } } @@ -59,7 +57,7 @@ int External() { std::cout << Form("# events injected with %d quark pair: ", checkPdgQuark) << nEventsInj << "\n"; std::cout << Form("# %d (anti)quarks: ", checkPdgQuark) << nQuarks << "\n"; std::cout <<"# signal hadrons: " << nSignals << "\n"; - std::cout <<"# signal hadrons decaying in the correct channel: " << nSignalGoodDecay << "\n"; + std::cout <<"# signal hadrons decaying : " << checkHadronDecays << "\n"; if (nEventsInj < nEvents) { std::cerr << "Number of generated events with triggered events different than expected\n"; @@ -76,9 +74,8 @@ int External() { return 1; } - float fracForcedDecays = float(nSignalGoodDecay) / nSignals; - if (fracForcedDecays < 0.9) { // we put some tolerance (it should not happen, but to be conservative) - std::cerr << "Fraction of signals decaying into the correct channel " << fracForcedDecays << " lower than expected\n"; + if (checkHadronDecays > 0) { + std::cerr << "Decayed OmegaC, it should never decay\n"; return 1; }