From a76f58ecf3b6e2e49c5fd06b254ae4eb0baeb98a Mon Sep 17 00:00:00 2001 From: DanielSamitz Date: Tue, 23 Jul 2024 16:08:13 +0200 Subject: [PATCH 1/2] no forced b-decays --- .../PWGEM/ini/Pythia8_Beauty_Cocktail.ini | 4 +-- .../PWGEM/ini/Pythia8_Charm_Cocktail.ini | 2 +- ...onic.cfg => force_semileptonic_beauty.cfg} | 22 ---------------- .../decayer/force_semileptonic_charm.cfg | 26 +++++++++++++++++++ 4 files changed, 29 insertions(+), 25 deletions(-) rename MC/config/PWGEM/pythia8/decayer/{force_semileptonic.cfg => force_semileptonic_beauty.cfg} (75%) create mode 100644 MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg diff --git a/MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail.ini b/MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail.ini index ffc32dca6..a1b41e14a 100644 --- a/MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail.ini +++ b/MC/config/PWGEM/ini/Pythia8_Beauty_Cocktail.ini @@ -1,6 +1,6 @@ [GeneratorExternal] fileName = ${O2DPG_ROOT}/MC/config/PWGEM/external/generator/Generator_pythia8_forcedDecays.C -funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;511;521;531;5122;5132;5232;5332") +funcName=GeneratePythia8ForcedDecays("411;421;431;4122;4232;4132;4332;") [GeneratorPythia8] config = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/generator/pythia8_hf_cocktail.cfg @@ -10,4 +10,4 @@ includePartonEvent=true [DecayerPythia8] config[0] = ${O2DPG_ROOT}/MC/config/common/pythia8/decayer/base.cfg -config[1] = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg +config[1] = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg diff --git a/MC/config/PWGEM/ini/Pythia8_Charm_Cocktail.ini b/MC/config/PWGEM/ini/Pythia8_Charm_Cocktail.ini index 57f79129d..23f36b30c 100644 --- a/MC/config/PWGEM/ini/Pythia8_Charm_Cocktail.ini +++ b/MC/config/PWGEM/ini/Pythia8_Charm_Cocktail.ini @@ -10,4 +10,4 @@ includePartonEvent=true [DecayerPythia8] config[0] = ${O2DPG_ROOT}/MC/config/common/pythia8/decayer/base.cfg -config[1] = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg +config[1] = ${O2DPG_ROOT}/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg diff --git a/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg b/MC/config/PWGEM/pythia8/decayer/force_semileptonic_beauty.cfg similarity index 75% rename from MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg rename to MC/config/PWGEM/pythia8/decayer/force_semileptonic_beauty.cfg index 8cf853b1e..a9f15fff9 100644 --- a/MC/config/PWGEM/pythia8/decayer/force_semileptonic.cfg +++ b/MC/config/PWGEM/pythia8/decayer/force_semileptonic_beauty.cfg @@ -1,25 +1,3 @@ -### only semileptonic decays for charm -### D+ -411:oneChannel = 1 0.087 0 -311 -11 12 -411:addChannel = 1 0.040 0 -321 211 -11 12 -411:addChannel = 1 0.037 0 -313 -11 12 -### D0 -421:oneChannel = 1 0.035 0 -321 -11 12 -421:addChannel = 1 0.022 0 -323 -11 12 -421:addChannel = 1 0.016 0 -321 111 -11 12 -421:addChannel = 1 0.014 0 -311 -211 -11 12 -### Ds -431:oneChannel = 1 0.025 0 333 -11 12 -431:addChannel = 1 0.027 0 221 -11 12 -### Lambdac -4122:oneChannel = 1 0.036 0 3122 -11 12 -### chi_{c}^{+} -4232:oneChannel = 1 0.07 0 3322 -11 12 -### chi_{c}^{0} -4132:oneChannel = 1 0.014 0 3312 -11 12 -### Omega_{c} -4332:oneChannel = 1 0.01224 0 3334 -11 12 - ### only semileptonic decays for beauty ### B0 511:oneChannel = 1 0.0207000 0 12 -11 -411 diff --git a/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg b/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg new file mode 100644 index 000000000..d605ab18d --- /dev/null +++ b/MC/config/PWGEM/pythia8/decayer/force_semileptonic_charm.cfg @@ -0,0 +1,26 @@ +### only semileptonic decays for charm +### D+ +411:oneChannel = 1 0.087 0 -311 -11 12 +411:addChannel = 1 0.040 0 -321 211 -11 12 +411:addChannel = 1 0.037 0 -313 -11 12 +### D0 +421:oneChannel = 1 0.035 0 -321 -11 12 +421:addChannel = 1 0.022 0 -323 -11 12 +421:addChannel = 1 0.016 0 -321 111 -11 12 +421:addChannel = 1 0.014 0 -311 -211 -11 12 +### Ds +431:oneChannel = 1 0.025 0 333 -11 12 +431:addChannel = 1 0.027 0 221 -11 12 +### Lambdac +4122:oneChannel = 1 0.036 0 3122 -11 12 +### chi_{c}^{+} +4232:oneChannel = 1 0.07 0 3322 -11 12 +### chi_{c}^{0} +4132:oneChannel = 1 0.014 0 3312 -11 12 +### Omega_{c} +4332:oneChannel = 1 0.01224 0 3334 -11 12 + +# Correct OmegaC decay length (wrong in PYTHIA8 decay table) (mm/c) +4332:tau0 = 0.08000000000 +# Correct Lb decay length (wrong in PYTHIA8 decay table) +5122:tau0 = 4.41000e-01 \ No newline at end of file From b1352c9112176aada84d907c8f644682ffcec3ad Mon Sep 17 00:00:00 2001 From: DanielSamitz Date: Wed, 31 Jul 2024 14:11:26 +0200 Subject: [PATCH 2/2] check for forced charm decays --- MC/config/PWGEM/ini/tests/Pythia8_Beauty_Cocktail.C | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/MC/config/PWGEM/ini/tests/Pythia8_Beauty_Cocktail.C b/MC/config/PWGEM/ini/tests/Pythia8_Beauty_Cocktail.C index 42efb8fd8..ef3c3929e 100644 --- a/MC/config/PWGEM/ini/tests/Pythia8_Beauty_Cocktail.C +++ b/MC/config/PWGEM/ini/tests/Pythia8_Beauty_Cocktail.C @@ -1,6 +1,6 @@ int External() { - + std::string path{"o2sim_Kine.root"}; TFile file(path.c_str(), "READ"); if (file.IsZombie()) { @@ -24,11 +24,11 @@ int External() if (pdg == -5){ found_bbar = true; } - if ( abs(pdg) == 511 || abs(pdg) == 521 || abs(pdg) == 531 || abs(pdg) == 5122 || abs(pdg) == 5132 || abs(pdg) == 5232 || abs(pdg) == 5332){ + if ( abs(pdg) == 411 || abs(pdg) == 421 || abs(pdg) == 431 || abs(pdg) == 4122 || abs(pdg) == 4232 || abs(pdg) == 4132 || abs(pdg) == 4332){ int ifirstdaughter = track.getFirstDaughterTrackId(); int ilastdaughter = track.getLastDaughterTrackId(); if (ifirstdaughter == -1 || ilastdaughter == -1){ - std::cerr << "Found beauty hadron that does not have daughters" << "\n"; + std::cerr << "Found charm hadron that does not have daughters" << "\n"; return 1; } bool found_electron = false; @@ -47,7 +47,7 @@ int External() } } if (!found_electron){ - std::cerr << "Found beauty hadron that does not decay to electron" << "\n"; + std::cerr << "Found charm hadron that does not decay to electron" << "\n"; return 1; } } @@ -56,6 +56,6 @@ int External() std::cerr << "Found event without b-bbar pair" << "\n"; return 1; } - } + } return 0; }