From a2e65c977d12fcd2842c0a1e22453db12c75a7a9 Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Thu, 24 Sep 2026 14:22:22 +0200 Subject: [PATCH 1/3] Update decay table for b-hadron to open charm MC --- ...thia8_beautyhadronic_with_decays_Mode2.cfg | 79 +++++-------------- 1 file changed, 19 insertions(+), 60 deletions(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg index 3ac775eb6..c862a114b 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg @@ -48,24 +48,12 @@ BeamRemnants:saturation 5 ### add D+ decays absent in PYTHIA8 decay table and set BRs from PDG for other 411:oneChannel = 1 0.0752 0 -321 211 211 411:addChannel = 1 0.0104 0 -313 211 -411:addChannel = 1 0.0156 0 311 211 -411:addChannel = 1 0.0752 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi +411:addChannel = 1 0.0100 0 333 211 # to have the same amount of D+->KKpi and D+->Kpipi ## add Lc decays absent in PYTHIA8 decay table and set BRs from PDG for other 4122:oneChannel = 1 0.0196 100 2212 -313 4122:addChannel = 1 0.0108 100 2224 -321 4122:addChannel = 1 0.022 100 102134 211 4122:addChannel = 1 0.035 0 2212 -321 211 -4122:addChannel = 1 0.0159 0 2212 311 -### add Xic+ decays absent in PYTHIA8 decay table -4232:addChannel = 1 0.2 0 2212 -313 -4232:addChannel = 1 0.2 0 2212 -321 211 -4232:addChannel = 1 0.2 0 3324 211 -4232:addChannel = 1 0.2 0 3312 211 211 -### add Xic0 decays absent in PYTHIA8 decay table -4132:addChannel = 1 0.0143 0 3312 211 -### add OmegaC decays absent in PYTHIA8 decay table -4332:addChannel = 1 0.5 0 3334 211 -4332:addChannel = 1 0.5 0 3312 211 ### K* -> K pi 313:onMode = off @@ -82,25 +70,12 @@ BeamRemnants:saturation 5 ### for Lambda_c -> Lambda(1520) K- 102134:onMode = off 102134:onIfAll = 2212 321 -### for Xic0 -> pi Xi -> pi pi Lambda -> pi pi pi p -### and Omega_c -> pi Xi -> pi pi Lambda -> pi pi pi p -3312:onMode = off -3312:onIfAll = 3122 -211 -3122:onMode = off -3122:onIfAll = 2212 -211 -### for Omega_c -> pi Omega -> pi K Lambda -> pi K pi p -3334:onMode = off -3334:onIfAll = 3122 -321 ### switch off all decay channels 411:onMode = off 421:onMode = off 431:onMode = off 4122:onMode = off -4232:onMode = off -4132:onMode = off -443:onMode = off -4332:onMode = off 511:onMode = off 521:onMode = off 531:onMode = off @@ -131,60 +106,44 @@ BeamRemnants:saturation 5 4122:onIfMatch = 102134 211 ### Lambda_c -> p K pi 4122:onIfMatch = 2212 321 211 -### Lambda_c -> pK0s -4122:onIfMatch = 2212 311 - -### Xic+ -> pK*0 -4232:onIfMatch = 2212 313 -### Xic+ -> p K- pi+ -4232:onIfMatch = 2212 321 211 -### Xic+ -> Xi*0 pi+, Xi*->Xi- pi+ -4232:onIfMatch = 3324 211 -### Xic+ -> Xi- pi+ pi+ -4232:onIfMatch = 3312 211 211 - -### Xic0 -> Xi- pi+ -4132:onIfMatch = 3312 211 - -### Omega_c -> Omega pi -4332:onIfMatch = 3334 211 -### Omega_c -> Xi pi -4332:onIfMatch = 3312 211 ### Force also golden beauty hadrons decay modes for D2H studies ### add B0 decays 511:oneChannel = 1 0.4 0 -411 211 -511:addChannel = 1 0.3 0 -413 211 -511:addChannel = 1 0.2 0 431 -211 +511:addChannel = 1 0.4 0 -413 211 ### add B+ decays -521:oneChannel = 1 0.8 0 -421 211 +521:oneChannel = 1 0.7 0 -421 211 ### add Bs0 decays -531:oneChannel = 1 0.8 0 -431 211 +531:oneChannel = 1 0.7 0 -431 211 ### add Lb decays -5122:oneChannel = 1 0.8 0 4122 -211 +5122:oneChannel = 1 0.6 0 4122 -211 ### we also add channels useful for studies of partly reconstructed decays / correlated backgrounds ### add B0 decays 511:addChannel = 1 0.1 0 -411 213 +511:addChannel = 1 0.1 0 -411 321 ### add B+ decays 521:addChannel = 1 0.1 0 -421 213 521:addChannel = 1 0.1 0 -423 211 +521:addChannel = 1 0.1 0 -421 321 ### add Bs0 decays 531:addChannel = 1 0.1 0 -431 213 531:addChannel = 1 0.1 0 -433 211 +531:addChannel = 1 0.1 0 -433 321 ### add Lb decays -5122:addChannel = 1 0.1 0 4122 -213 -5122:addChannel = 1 0.05 0 4122 211 -211 -211 -5122:addChannel = 1 0.05 0 4212 -211 +5122:addChannel = 1 0.10 0 14122 -211 +5122:addChannel = 1 0.10 0 4124 -211 +5122:addChannel = 1 0.10 0 4212 -211 +5122:addChannel = 1 0.10 0 4122 -321 ### B0 -> D pi 511:onIfMatch = 411 211 ### B0 -> D* pi 511:onIfMatch = 413 211 -### B0 -> Ds pi -511:onIfMatch = 431 211 ### B0 -> D rho 511:onIfMatch = 411 213 +### B0 -> D rho +511:onIfMatch = 411 321 ### B+ -> D0 pi 521:onIfMatch = 421 211 @@ -202,9 +161,9 @@ BeamRemnants:saturation 5 ### Lb -> Lc pi 5122:onIfMatch = 4122 211 -### Lb -> Lc rho -5122:onIfMatch = 4122 213 -### Lb -> Lc pi pi pi -5122:onIfMatch = 4122 211 211 211 -### Lb -> Sc pi +### Lb -> Lc(2595)+ pi +5122:onIfMatch = 14122 211 +### Lb -> Lc(2625)+ pi +5122:onIfMatch = 4124 211 +### Lb -> Sc(2455)+ pi 5122:onIfMatch = 4212 211 From 6093c1faec616427818e2ef3f435212140e29716 Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Thu, 24 Sep 2026 14:37:03 +0200 Subject: [PATCH 2/3] Update test macro --- .../GeneratorHF_D2H_bbbar_Bforced_gap5_Mode2.C | 13 +++++-------- .../pythia8_beautyhadronic_with_decays_Mode2.cfg | 8 +++++--- 2 files changed, 10 insertions(+), 11 deletions(-) diff --git a/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Bforced_gap5_Mode2.C b/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Bforced_gap5_Mode2.C index 629aa3ae1..7a5fdab70 100644 --- a/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Bforced_gap5_Mode2.C +++ b/MC/config/PWGHF/ini/tests/GeneratorHF_D2H_bbbar_Bforced_gap5_Mode2.C @@ -6,17 +6,14 @@ int External() { std::vector checkPdgHadron{411, 421, 431, 4122, 4132, 4232, 4332, 511, 521, 531, 5122}; std::map>> checkHadronDecays{ // sorted pdg of daughters - {411, {{-321, 211, 211}, {-313, 211}, {211, 311}, {211, 333}}}, // D+ + {411, {{-321, 211, 211}, {-313, 211}, {211, 333}}}, // D+ {421, {{-321, 211}, {-321, 111, 211}}}, // D0 {431, {{211, 333}, {-313, 321}}}, // Ds+ - {4122, {{-313, 2212}, {-321, 2224}, {211, 102134}, {-321, 211, 2212}, {311, 2212}}}, // Lc+ - {4132, {{211, 3312}}}, // Xic0 - {4232, {{-313, 2212}, {-321, 3324}, {211, 211, 3312}, {-321, 211, 2212}}}, // Xic+ - {4332, {{211, 3334}}}, // Omegac+ - {511, {{-411, 211}, {-413, 211}, {-411, 213}, {431, -211}}}, // B0 - {521, {{-421, 211}, {-423, 211}, {-421, 213}}}, // B+ + {4122, {{-313, 2212}, {-321, 2224}, {211, 102134}, {-321, 211, 2212}}}, // Lc+ + {511, {{-411, 211}, {-413, 211}, {-411, 213}, {-411, 321}, {-413, 213}, {-413, 321}}}, // B0 + {521, {{-421, 211}, {-423, 211}, {-421, 213}, {-421, 321}}}, // B+ {531, {{-431, 211}, {-433, 211}, {-431, 213}}}, // Bs0 - {5122, {{4122, -211}, {4122, -213}, {4122, 211, -211, -211}, {4212, -211}}} // Lb0 + {5122, {{4122, -211}, {4122, -213}, {14122, -211}, {4124, -211}, {4212, -211}}} // Lb0 }; TFile file(path.c_str(), "READ"); diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg index c862a114b..75815f5af 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg @@ -109,8 +109,8 @@ BeamRemnants:saturation 5 ### Force also golden beauty hadrons decay modes for D2H studies ### add B0 decays -511:oneChannel = 1 0.4 0 -411 211 -511:addChannel = 1 0.4 0 -413 211 +511:oneChannel = 1 0.35 0 -411 211 +511:addChannel = 1 0.35 0 -413 211 ### add B+ decays 521:oneChannel = 1 0.7 0 -421 211 ### add Bs0 decays @@ -121,7 +121,9 @@ BeamRemnants:saturation 5 ### we also add channels useful for studies of partly reconstructed decays / correlated backgrounds ### add B0 decays 511:addChannel = 1 0.1 0 -411 213 -511:addChannel = 1 0.1 0 -411 321 +511:addChannel = 1 0.1 0 -413 213 +511:addChannel = 1 0.05 0 -411 321 +511:addChannel = 1 0.05 0 -411 321 ### add B+ decays 521:addChannel = 1 0.1 0 -421 213 521:addChannel = 1 0.1 0 -423 211 From f1cb404f2e4a6a6e782c88827dca75dc3c0748d8 Mon Sep 17 00:00:00 2001 From: Fabrizio Grosa Date: Thu, 24 Sep 2026 14:37:55 +0200 Subject: [PATCH 3/3] Fix typo --- .../generator/pythia8_beautyhadronic_with_decays_Mode2.cfg | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg b/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg index 75815f5af..0cb1bbaf5 100644 --- a/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg +++ b/MC/config/PWGHF/pythia8/generator/pythia8_beautyhadronic_with_decays_Mode2.cfg @@ -123,7 +123,7 @@ BeamRemnants:saturation 5 511:addChannel = 1 0.1 0 -411 213 511:addChannel = 1 0.1 0 -413 213 511:addChannel = 1 0.05 0 -411 321 -511:addChannel = 1 0.05 0 -411 321 +511:addChannel = 1 0.05 0 -413 321 ### add B+ decays 521:addChannel = 1 0.1 0 -421 213 521:addChannel = 1 0.1 0 -423 211