arxiv: 2602.18378 · v4 · submitted 2026-02-20 · ✦ hep-ph · hep-ex

Recognition: 1 theorem link

· Lean Theorem

On the simulated kinematic distributions of heavy meson decays

Florian Herren , Raynette van Tonder

Authors on Pith no claims yet

Pith reviewed 2026-05-15 20:25 UTC · model grok-4.3

classification ✦ hep-ph hep-ex

keywords EvtGenphase-space samplingsemileptonic B decaysresonanceskinematic distributionsMonte Carlo simulationreweightingflavour physics

0 comments

The pith

EvtGen's phase-space sampling produces unphysical features in kinematic distributions of semileptonic B decays involving resonances.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines the phase-space sampling algorithm inside EvtGen, a Monte Carlo generator widely used for heavy-hadron decays. It demonstrates that the algorithm omits phase-space factors when resonances appear in the decay chain of semileptonic B mesons. These omissions create artificial structures in the simulated distributions of particle momenta and invariant masses. The authors supply a practical reweighting procedure that adjusts the hadronic invariant-mass spectrum to restore physical shapes in existing samples. Accurate kinematic modeling is required for reliable extraction of parameters from modern flavour-physics data sets.

Core claim

EvtGen's phase-space sampling algorithm generates unphysical features in kinematic distributions of semileptonic B decays involving resonances, originating from neglected phase-space factors. Reweighting the hadronic invariant mass distribution supplies a short-term correction for the affected simulated samples.

What carries the argument

EvtGen phase-space sampling algorithm, which neglects resonance-related phase-space factors during event generation.

If this is right

Kinematic distributions in simulated semileptonic B decays that include resonances will exhibit artificial peaks or distortions.
Analyses that rely on these distributions for background estimation or signal extraction may inherit systematic biases.
Reweighting the hadronic invariant mass restores the expected physical shapes in existing EvtGen samples.
The sampling issue appears whenever resonances are treated in the decay chain of heavy mesons.
Measurements of branching fractions or angular observables in resonance-mediated modes require the correction to reach target precision.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Similar sampling shortcuts may exist in other Monte Carlo packages that generate multi-body decays with intermediate resonances.
Incorporating the full phase-space factors directly into EvtGen's sampling routine would eliminate the need for post-hoc reweighting.
The correction could be tested on additional decay modes such as those involving higher resonances or different lepton flavours.
Validation against analytic phase-space integrals for simple resonance cases would quantify the size of the residual discrepancy.

Load-bearing premise

The unphysical features arise specifically from neglected phase-space factors in EvtGen's algorithm and can be fully removed by reweighting the hadronic invariant mass without side effects.

What would settle it

Generate the same semileptonic B decay with resonances using an independent generator that retains all phase-space factors and compare the resulting kinematic distributions to the uncorrected and reweighted EvtGen output.

read the original abstract

Modern measurements in flavour physics rely on accurate simulations of signal and background processes, provided by a wide range of general-purpose and specialised Monte-Carlo event generators. Due to the inclusion of a larger amount of specialised decays of heavy hadrons, EvtGen is often the tool of choice for many scenarios. We investigate the phase-space sampling algorithm of EvtGen and demonstrate that it generates unphysical features in kinematic distributions of semileptonic $B$ decays involving resonances, originating from neglected phase-space factors. We provide a short-term solution to correct the affected simulated samples through reweighting of the hadronic invariant mass distribution.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

EvtGen's phase-space sampler skips factors in resonant semileptonic B decays and creates unphysical kinks; the authors give a hadronic-mass reweight as a quick patch.

read the letter

The paper flags that EvtGen's sampling algorithm for semileptonic B decays with resonances neglects phase-space factors, which produces artificial structures in the kinematic distributions. They show this and offer reweighting the hadronic invariant mass as a short-term fix for existing samples. This is the core takeaway: a concrete, usable warning for anyone who relies on EvtGen for flavor-physics Monte Carlo. What the work does well is identify a real limitation in a tool that most B-physics analyses actually run, rather than a toy generator. The abstract is direct about the origin of the problem and the proposed correction, which keeps the scope narrow and practical. Groups at LHCb or Belle II who simulate resonant modes will recognize the issue and can apply the reweight immediately. The soft spot is the one the stress test raises. Reweighting only the hadronic mass restores that marginal, but the full phase-space density involves correlated variables such as q² and angles. Unless the sampling error factorizes exactly, other distributions may still be off. The abstract does not detail how they validated the joint distributions against a reference integrator, so the full paper must show that the fix does not introduce new biases elsewhere. If the validation is only one-dimensional, the claim is weaker than it first appears. This paper is for simulation users and analysts who generate or fit EvtGen samples for semileptonic decays. A reader who cares about MC accuracy in current flavor measurements will get immediate value from the warning and the patch. It deserves peer review because the claim is testable with public tools and the affected analyses are common; referees can check the validation plots and decide how complete the correction really is. I would send it rather than desk-reject.

Referee Report

3 major / 2 minor

Summary. The manuscript investigates the phase-space sampling algorithm implemented in EvtGen and shows that it produces unphysical features in the kinematic distributions of semileptonic B decays involving resonances. These features are attributed to the neglect of phase-space factors in the sampling procedure. The authors propose a practical correction by reweighting simulated events according to the hadronic invariant mass distribution.

Significance. If the diagnosis and the reweighting prescription are validated, the result is significant for the flavour-physics community because EvtGen is the default generator for many B-physics analyses at LHCb, Belle II and other experiments. Unphysical distortions in resonance decays can bias acceptance corrections, form-factor extractions and background estimates. The provision of an immediately applicable reweighting recipe is a concrete, short-term contribution that does not require modification of the generator source code.

major comments (3)

[§4] §4 (reweighting procedure): the assertion that one-dimensional reweighting of the hadronic invariant mass fully restores the correct kinematics is not demonstrated. In semileptonic decays the phase-space density depends on correlated variables (m_had, q², lepton angles); no evidence is shown that the joint distributions or the marginals in lepton energy and angular variables match those obtained from a reference integrator after reweighting.
[§3.1] §3.1 (algorithm description): the precise location where the phase-space factor is omitted is not identified with an equation or pseudocode excerpt. Without this, it is impossible to verify that the observed distortions are caused exactly by the claimed omission rather than by another implementation detail.
[§5] §5 (validation): no quantitative comparison is presented against an independent phase-space generator (e.g., a flat or importance-sampled integrator that includes the full Dalitz-plot density). Such a benchmark is required to confirm both the unphysical nature of the features and the absence of side-effects from the reweighting.

minor comments (2)

[Abstract] The abstract states that the features appear in 'semileptonic B decays involving resonances' but does not list the specific decay modes or resonance channels examined; this information should be added for reproducibility.
[Figures] Figure captions should explicitly state whether the plotted distributions are before or after reweighting and what reference distribution (if any) is overlaid.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive comments, which we believe will improve the clarity and robustness of the presented results. We address each major comment below and will incorporate the suggested changes in the revised version.

read point-by-point responses

Referee: §4 (reweighting procedure): the assertion that one-dimensional reweighting of the hadronic invariant mass fully restores the correct kinematics is not demonstrated. In semileptonic decays the phase-space density depends on correlated variables (m_had, q², lepton angles); no evidence is shown that the joint distributions or the marginals in lepton energy and angular variables match those obtained from a reference integrator after reweighting.

Authors: We agree that the current manuscript does not explicitly demonstrate restoration of the full joint distributions. Although the reweighting is motivated by the phase-space density and targets the variable where the sampling omission occurs, additional validation is warranted to confirm no residual biases in correlated observables. In the revised manuscript we will add plots and quantitative comparisons of the q², lepton energy, and angular distributions before and after reweighting, together with direct comparisons against a reference phase-space integrator for the joint (m_had, q²) plane. revision: yes
Referee: §3.1 (algorithm description): the precise location where the phase-space factor is omitted is not identified with an equation or pseudocode excerpt. Without this, it is impossible to verify that the observed distortions are caused exactly by the claimed omission rather than by another implementation detail.

Authors: We acknowledge that a more explicit identification of the omission would strengthen the paper. The phase-space factor is omitted inside the hadronic-mass sampling step of EvtGen’s semileptonic decay classes (specifically, the routine that generates m_had without the appropriate kinematic weight before the subsequent lepton and neutrino sampling). In the revised manuscript we will insert a short pseudocode excerpt of the relevant sampling loop together with the line number reference to the EvtGen source to make the precise location unambiguous. revision: yes
Referee: §5 (validation): no quantitative comparison is presented against an independent phase-space generator (e.g., a flat or importance-sampled integrator that includes the full Dalitz-plot density). Such a benchmark is required to confirm both the unphysical nature of the features and the absence of side-effects from the reweighting.

Authors: We agree that an independent benchmark is necessary for full validation. In the revised manuscript we will include a quantitative comparison against a standalone phase-space integrator that samples the full Dalitz-plot density. This will comprise overlaid histograms, Kolmogorov-Smirnov distances, and pull distributions for all relevant kinematic variables both before and after reweighting, thereby confirming the origin of the artifacts and the absence of side-effects. revision: yes

Circularity Check

0 steps flagged

No significant circularity; external analysis of public tool

full rationale

The paper examines the phase-space sampling in the external EvtGen package for semileptonic B decays with resonances, identifies discrepancies with expected physical distributions, and offers a practical reweighting fix in hadronic mass. No derivation chain reduces a claimed prediction or uniqueness result to the paper's own fitted inputs or self-citations. The central claim rests on direct inspection of the algorithm and comparison to reference integrators rather than self-referential definitions or ansatze smuggled via prior work by the same authors. This is the normal case of an independent diagnostic study.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract provides no information on free parameters, axioms, or invented entities; all fields left empty due to lack of detail.

pith-pipeline@v0.9.0 · 5391 in / 1062 out tokens · 32747 ms · 2026-05-15T20:25:19.198811+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Foundation/RealityFromDistinction.lean reality_from_one_distinction unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We investigate the phase-space sampling algorithm of EvtGen and demonstrate that it generates unphysical features in kinematic distributions of semileptonic B decays involving resonances, originating from neglected phase-space factors. We provide a short-term solution to correct the affected simulated samples through reweighting of the hadronic invariant mass distribution.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Extracting production fractions of $b$ hadrons from exclusive semi-leptonic decays
hep-ph 2026-04 unverdicted novelty 7.0

A method to extract fs/fd and similar ratios from exclusive semi-leptonic B decays yields a 7% uncertainty constraint using existing measurements.
Closing the knowledge gap in semileptonic $B\rightarrow X_c\ell\nu$ decays
hep-ph 2026-05 unverdicted novelty 4.0

The unmeasured fraction of inclusive semileptonic B to charm decays is dominated by final states without D mesons, likely from baryons and Ds mesons.