In python/pythia8_conf.py (configurerpvsusy), beautyhistograms = ["b_mu", "b_tau", "b0_nu_mu", "b0_nu_tau"] feeds getmaxsumbrrpvsusy, but only the b_tau (B+) and b0_nu_tau (B0) channels are ever configured via addChannel.
The bench3/bench4 data files (branchingratiosrpvsusybench{3,4}.dat) contain only b0_nu_mu. So maxsumBR > 0 passes the exit_if_zero_br guard, yet both getbr_rpvsusy(h, "b_tau", ...) and getbr_rpvsusy(h, "b0_nu_tau", ...) return 0 — no HNL production channel is added at all, and both B mesons get freshly-overwritten empty decay tables.
The fix is presumably to add the muon-mode channels mirroring the tau modes:
- B+ → N μ+ :
521:addChannel 1 {br} 0 9900015 -13 using getbr_rpvsusy(h, "b_mu", ...)
- B0 → N ν_μ :
511:addChannel 1 {br} 22 9900015 14 using getbr_rpvsusy(h, "b0_nu_mu", ...)
(with the corresponding sumBR += accumulation for the 22 22 filler).
An RPV-SUSY expert should confirm the intended channel definitions and PDG assignments before this generator config is changed. Found during a code review of the tracked tree.
In
python/pythia8_conf.py(configurerpvsusy),beautyhistograms = ["b_mu", "b_tau", "b0_nu_mu", "b0_nu_tau"]feedsgetmaxsumbrrpvsusy, but only theb_tau(B+) andb0_nu_tau(B0) channels are ever configured viaaddChannel.The bench3/bench4 data files (
branchingratiosrpvsusybench{3,4}.dat) contain onlyb0_nu_mu. SomaxsumBR > 0passes theexit_if_zero_brguard, yet bothgetbr_rpvsusy(h, "b_tau", ...)andgetbr_rpvsusy(h, "b0_nu_tau", ...)return 0 — no HNL production channel is added at all, and both B mesons get freshly-overwritten empty decay tables.The fix is presumably to add the muon-mode channels mirroring the tau modes:
521:addChannel 1 {br} 0 9900015 -13usinggetbr_rpvsusy(h, "b_mu", ...)511:addChannel 1 {br} 22 9900015 14usinggetbr_rpvsusy(h, "b0_nu_mu", ...)(with the corresponding
sumBR +=accumulation for the22 22filler).An RPV-SUSY expert should confirm the intended channel definitions and PDG assignments before this generator config is changed. Found during a code review of the tracked tree.