[Paper Review] Study of variants for Monte Carlo generators of $ au o 3\pi u$ decays

Z. Wąs, J. Zaremba|arXiv (Cornell University)|Aug 26, 2015

Particle physics theoretical and experimental studies1 references3 citations

TL;DR

This paper evaluates and compares multiple hadronic current parameterizations used in Monte Carlo generators for τ → 3πν decay, focusing on their theoretical foundations, numerical implementations, and compatibility with experimental data. It demonstrates that one-dimensional distributions are insufficient to constrain model parameters accurately, and three-dimensional phase-space distributions—as used by CLEO—are essential for reliable fitting, especially when testing models like Resonance Chiral Lagrangian (RChL) against data with high precision.

ABSTRACT

Low energy QCD (below 2 GeV) is a region of resonance dynamics, sometimes lacking satisfactory description as compared to precision of available experimental data. Hadronic $\ au$ decays offer a probe for such energy regime. In general, predictions for decays are model dependent, with parameters fitted to experimental results. Parameterizations differ by amount of assumptions and theoretical requirements taken into account. Both model distributions and acquired data samples used for fits are results of complex effort. In this paper, we investigate main parameterizations of $\ au$ decay matrix elements for the one- and three-prong channels of three-pion $\ au$ decays. Differences in analytical forms of the currents and resulting distributions used for comparison with the experimental data are studied. We use invariant mass spectra of all possible pion pairs and the whole three-pion system. Also three-dimensional histograms spanned over all distinct squared invariant masses are used to represent results of models and experimental data. We present distributions from {\ t TAUOLA} Monte Carlo generation and semi-analytical calculation. These are necessary steps in development for fitting in as model-independent way as possible, and to explore multi-million event experimental data samples. This includes response of distributions to model variants, and/or numerical values of parameters. Interference effects of currents parts are also studied. For technical purposes, weighted events are introduced. Even though we focus on $3\\pi\ u_\ au$ modes, technical aspects of our study are relevant for all $\ au$ decay modes into three hadrons.

Motivation & Objective

• To evaluate and document the differences among various hadronic current parameterizations used in TAUOLA for τ → 3πν decays.
• To assess the impact of using one-dimensional vs. three-dimensional experimental distributions on model fitting and parameter constraints.
• To investigate the predictive power and consistency of model-independent approaches, such as those based on Resonance Chiral Lagrangian (RChL), with high-precision data.
• To highlight the limitations of current experimental data representations in constraining complex hadronic current structures, especially when phase-space dimensionality is reduced.
• To support the development of flexible, well-documented Monte Carlo models that can be updated efficiently as new high-precision data become available.

Proposed method

• Systematic comparison of multiple hadronic current variants used in TAUOLA, including CLEO, BaBar, and RChL-based models.
• Use of Monte Carlo event generation to compute matrix element squared differences (|wt−1|) across event samples to quantify model discrepancies.
• Adaptation of model-independent fitting techniques from Ref. [4] to relativistic τ decays where ντ momentum is not reconstructible.
• Semi-analytical calculation methods for multidimensional distributions, enabling efficient comparison with experimental histograms.
• Numerical evaluation of interference effects and contributions from different current components to understand model sensitivity.
• Use of one-dimensional histograms (from CLEO and BaBar) and three-dimensional distributions (from CLEO) as input for model comparisons.

Experimental results

Research questions

• RQ1 How do different hadronic current parameterizations in TAUOLA affect the prediction of τ → 3πν decay distributions?
• RQ2 To what extent do one-dimensional experimental distributions fail to constrain model parameters compared to three-dimensional distributions?
• RQ3 What are the numerical consequences of interference terms between different current components in τ → 3πν decays?
• RQ4 How well do RChL-based models agree with experimental data when constrained by one-dimensional vs. three-dimensional distributions?
• RQ5 What are the implications of phase-space constraints and limited data dimensionality for the reliability of model fits in τ decays?

Key findings

• One-dimensional histograms used in CLEO and BaBar fits are insufficient to fully constrain model parameters, as evidenced by discrepancies in predictions for π⁻π⁻ invariant mass.
• Three-dimensional distributions—as used by CLEO—provide significantly better constraints on hadronic currents and are necessary for reliable model fitting.
• A factor of two discrepancy exists between RChL and CLEO model predictions in certain bins of three-dimensional distributions, indicating that 1D data cannot capture full model differences.
• Even though RChL models agree with data at the ~30% level (1/N_C), experimental precision is an order of magnitude better, necessitating higher-dimensional data for accurate fitting.
• Interference terms between current components significantly affect final distributions, and their effects are not captured by low-dimensional data representations.
• Model developers and experimentalists must align on data representation and model flexibility to ensure reliable parameter constraints and future adaptability.