A dispersive analysis of the pion vector form factor and $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ decay

TL;DR

This paper performs a dispersive analysis of the pion vector form factor and tau decays to extract properties of the $ ho(1450)$ and $ ho(1700)$ resonances, emphasizing the importance of systematic errors and combining data from different decay channels.

Contribution

It introduces an improved method for resonance parameter extraction using dispersion relations and combined decay data, highlighting the dominance of systematic errors over statistical uncertainties.

Findings

Systematic errors dominate over fit uncertainties in resonance parameters.

Provides detailed form factor modulus and phase with errors as ancillary material.

First analysis of $ au^{-} o K^{-}K_{S} u_{ au}$ decay spectrum with implications for resonance studies.

Abstract

We explore the potential of a combined analysis of the decays $\tau^{-}\to\pi^{-}\pi^{0}\nu_{\tau}$ and $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ in the determination of the $\rho(1450)$ and $\rho(1700)$ resonance properties in the frame of resonance Chiral Theory supplemented by dispersion relations. On the one hand, we take advantage of the very precise data on the modulus squared of the pion vector form factor $|F_{V}^{\pi}|^{2}$ obtained by Belle to carry out a very dedicated analysis of the region where these resonances come up into play. Our study provides an improved treatment of the systematic theoretical errors and, as a most important result, we conclude that they dominate over the fit uncertainties in the determination of the $\rho(1450)$ and $\rho(1700)$ pole parameters and tend to be larger than in other determinations quoted in the literature where these errors were ignored or underestimated. The results of our analysis are summarized in numerical tables for the form factor modulus and phase, including both statistical and systematic errors, that can be found as ancillary material of this paper. As a byproduct, we also determine the low-energy observables of the pion vector form factor and the $\rho$-pole position. On the other hand, we benefit from the recent experimental data for the transition $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ released by BaBar to perform a first analysis of its decay spectrum and discuss the role of these resonances in this decay. We point out that higher-quality data on the $K^{-}K_{S}$ decay channel will allow compete with the $|F_{V}^{\pi}|^{2}$ ones and improve the determination of the $\rho(1450)$ and $\rho(1700)$ resonance parameters as a result of a combined analysis. We hope our study to be of interest for present and future experimental analysis of these decays.

Tables

Table 1. Table 1: Results for the fits obtained with a three-times-subtracted dispersion relation including three vector resonances in F V π ​ ( s ) superscript subscript 𝐹 𝑉 𝜋 𝑠 F_{V}^{\pi}(s) according to Eq. ( 2.19 ) for four representative values of s cut subscript 𝑠 cut s_{\rm{cut}} in the dispersive integral.

	Parameter	$s_{\mathrm{cut}}$ [GeV2]
Fits		$m_{\tau }^2$	$4$ (reference fit)	$10$	$\infty$
Fit 1	${\alpha }_1$ [GeV-2]	$1.87\left(1\right)$	$1.88\left(1\right)$	$1.89\left(1\right)$	$1.89\left(1\right)$
	${\alpha }_2$ [GeV-4]	$4.40\left(1\right)$	$4.34\left(1\right)$	$4.32\left(1\right)$	$4.32\left(1\right)$
	$m_{\rho }$ [MeV]	$=773.6\left(9\right)$	$=773.6\left(9\right)$	$=773.6\left(9\right)$	$=773.6\left(9\right)$
	$M_{\rho }$ [MeV]	$=m_{\rho }$	$=m_{\rho }$	$=m_{\rho }$	$=m_{\rho }$
	$M_{{\rho }^{\prime }}$ [MeV]	$1365\left(15\right)$	$1376\left(6\right)$	$1313\left(15\right)$	$1311\left(5\right)$
	${\Gamma }_{{\rho }^{\prime }}$ [MeV]	$562\left(55\right)$	$603\left(22\right)$	$700\left(6\right)$	$701\left(28\right)$
	$M_{{\rho }^{\prime \prime }}$[MeV]	$1727\left(12\right)$	$1718\left(4\right)$	$1660\left(9\right)$	$1658\left(1\right)$
	${\Gamma }_{{\rho }^{\prime \prime }}$ [MeV]	$278\left(1\right)$	$465\left(9\right)$	$601\left(39\right)$	$602\left(3\right)$
	$\gamma$	$0.12\left(2\right)$	$0.15\left(1\right)$	$0.16\left(1\right)$	$0.16\left(1\right)$
	${\varphi }_1$	$-0.69\left(1\right)$	$-0.66\left(1\right)$	$-1.36\left(10\right)$	$-1.39\left(1\right)$
	$\delta$	$-0.09\left(1\right)$	$-0.13\left(1\right)$	$-0.16\left(1\right)$	$-0.17\left(1\right)$
	${\varphi }_2$	$-0.17\left(5\right)$	$-0.44\left(3\right)$	$-1.01\left(5\right)$	$-1.03\left(2\right)$
	${\chi }^2$/d.o.f	$1.47$	$0.70$	$0.64$	$0.64$
Fit 1-$\rho$	${\alpha }_1$ [GeV-2]	$1.88\left(1\right)$	$1.88\left(1\right)$	$1.89\left(1\right)$	$1.88\left(1\right)$
	${\alpha }_2$ [GeV-4]	$4.37\left(3\right)$	$4.34\left(1\right)$	$4.31\left(3\right)$	$4.34\left(1\right)$
	$m_{\rho }$ [MeV]	$773.9\left(3\right)$	$773.8\left(3\right)$	$773.9\left(3\right)$	$773.9\left(3\right)$
	$M_{\rho }$ [MeV]	$=m_{\rho }$	$=m_{\rho }$	$=m_{\rho }$	$=m_{\rho }$
	$M_{{\rho }^{\prime }}$ [MeV]	$1382\left(71\right)$	$1375\left(11\right)$	$1316\left(9\right)$	$1312\left(8\right)$
	${\Gamma }_{{\rho }^{\prime }}$ [MeV]	$516\left(165\right)$	$608\left(35\right)$	$728\left(92\right)$	$726\left(26\right)$
	$M_{{\rho }^{\prime \prime }}$[MeV]	$1723\left(1\right)$	$1715\left(22\right)$	$1655\left(1\right)$	$1656\left(8\right)$
	${\Gamma }_{{\rho }^{\prime \prime }}$ [MeV]	$315\left(271\right)$	$455\left(16\right)$	$569\left(160\right)$	$571\left(13\right)$
	$\gamma$	$0.12\left(13\right)$	$0.16\left(1\right)$	$0.18\left(2\right)$	$0.17\left(1\right)$
	${\varphi }_1$	$-0.56\left(35\right)$	$-0.69\left(1\right)$	$-1.40\left(19\right)$	$-1.41\left(8\right)$
	$\delta$	$-0.09\left(3\right)$	$-0.13\left(1\right)$	$-0.17\left(4\right)$	$-0.17\left(3\right)$
	${\varphi }_2$	$-0.19\left(69\right)$	$-0.45\left(12\right)$	$-1.06\left(10\right)$	$-1.05\left(11\right)$
	${\chi }^2$/d.o.f	$1.09$	$0.70$	$0.63$	$0.66$

Table 2. Table 2: Results for the fits obtained with a three-times-subtracted dispersion relation including three vector resonances in F V π ​ ( s ) superscript subscript 𝐹 𝑉 𝜋 𝑠 F_{V}^{\pi}(s) according to Eq. ( 2.19 ) with s cut = 4 subscript 𝑠 cut 4 s_{\rm{cut}}=4 GeV 2 in the dispersive integral for four representative values of the matching point.

	Parameter	Matching point [GeV]
Fits		$0.85$	$0.9$	$0.95$	$1$ (reference fit)
Fit I	${\alpha }_1$ [GeV-2]	$1.88\left(1\right)$	$1.88\left(1\right)$	$1.88\left(1\right)$	$1.88\left(1\right)$
	${\alpha }_2$ [GeV-4]	$4.35\left(1\right)$	$4.35\left(1\right)$	$4.34\left(1\right)$	$4.34\left(1\right)$
	$m_{\rho }$ [MeV]	$=773.6\left(9\right)$	$=773.6\left(9\right)$	$=773.6\left(9\right)$	$=773.6\left(9\right)$
	$M_{\rho }$ [MeV]	$=m_{\rho }$	$=m_{\rho }$	$=m_{\rho }$	$=m_{\rho }$
	$M_{{\rho }^{\prime }}$ [MeV]	$1394\left(6\right)$	$1374\left(8\right)$	$1351\left(5\right)$	$1376\left(6\right)$
	${\Gamma }_{{\rho }^{\prime }}$ [MeV]	$592\left(19\right)$	$583\left(27\right)$	$592\left(2\right)$	$603\left(22\right)$
	$M_{{\rho }^{\prime \prime }}$[MeV]	$1733\left(9\right)$	$1715\left(1\right)$	$1697\left(3\right)$	$1718\left(4\right)$
	${\Gamma }_{{\rho }^{\prime \prime }}$ [MeV]	$562\left(3\right)$	$541\left(45\right)$	$486\left(7\right)$	$465\left(9\right)$
	$\gamma$	$0.12\left(1\right)$	$0.12\left(1\right)$	$0.13\left(1\right)$	$0.15\left(1\right)$
	${\varphi }_1$	$-0.44\left(3\right)$	$-0.60\left(1\right)$	$-0.80\left(1\right)$	$-0.66\left(1\right)$
	$\delta$	$-0.13\left(1\right)$	$-0.13\left(1\right)$	$-0.13\left(1\right)$	$-0.13\left(1\right)$
	${\varphi }_2$	$-0.38\left(3\right)$	$-0.51\left(2\right)$	$-0.62\left(1\right)$	$-0.44\left(3\right)$
	${\chi }^2$/d.o.f	$0.75$	$0.74$	$0.68$	$0.70$

Table 3. Table 3: Results for the fits obtained with a three-times-subtracted dispersion relation including three vector resonances in F V π ​ ( s ) superscript subscript 𝐹 𝑉 𝜋 𝑠 F_{V}^{\pi}(s) according to Eq. ( 2.19 ) with s cut = 4 subscript 𝑠 cut 4 s_{\rm{cut}}=4 GeV 2 in the dispersive integral with (second column) and without (last column) the K ​ K ¯ 𝐾 ¯ 𝐾 K\bar{K} channel in the ρ ′ superscript 𝜌 ′ \rho^{\prime} and ρ ′′ superscript 𝜌 ′′ \rho^{\prime\prime} energy-dependent resonance widths, and with the additional inclusion of the ρ ′ → ω ​ π → superscript 𝜌 ′ 𝜔 𝜋 \rho^{\prime}\to\omega\pi contribution into the ρ ′ superscript 𝜌 ′ \rho^{\prime} width (third column). See main text for details.

Parameter	$s_{\mathrm{cut}}=4$ GeV2
	Fit A	Fit B	reference fit
${\alpha }_1$ [GeV-2]	$1.87\left(1\right)$	$1.88\left(1\right)$	$1.88\left(1\right)$
${\alpha }_2$ [GeV-4]	$4.37\left(1\right)$	4.35(1)	$4.34\left(1\right)$
$m_{\rho }$ [MeV]	$=773.6\left(9\right)$	$=773.6\left(9\right)$	$=773.6\left(9\right)$
$M_{\rho }$ [MeV]	$=m_{\rho }$	$=m_{\rho }$	$=m_{\rho }$
$M_{{\rho }^{\prime }}$ [MeV]	$1373\left(5\right)$	$1441\left(3\right)$	$1376\left(6\right)$
${\Gamma }_{{\rho }^{\prime }}$ [MeV]	$462\left(14\right)$	$576\left(33\right)$	$603\left(22\right)$
$M_{{\rho }^{\prime \prime }}$[MeV]	$1775\left(1\right)$	$1733\left(9\right)$	$1718\left(4\right)$
${\Gamma }_{{\rho }^{\prime \prime }}$ [MeV]	$412\left(27\right)$	$349\left(52\right)$	$465\left(9\right)$
$\gamma$	$0.13\left(1\right)$	$0.15\left(3\right)$	$0.15\left(1\right)$
${\varphi }_1$	$-0.80\left(1\right)$	$-0.53\left(5\right)$	$-0.66\left(1\right)$
$\delta$	$-0.14\left(1\right)$	$-0.14\left(1\right)$	$-0.13\left(1\right)$
${\varphi }_2$	$-0.44\left(2\right)$	$-0.46\left(3\right)$	$-0.44\left(3\right)$
${\chi }^2$/d.o.f	$0.93$	$0.70$	$0.70$

Table 4. Table 4: Comparison between different results for the model parameters and corresponding pole positions for the ρ ′ superscript 𝜌 ′ \rho^{\prime} (upper table) and ρ ′′ superscript 𝜌 ′′ \rho^{\prime\prime} (lower table) resonances. The first and second errors refer, respectively, to the statistical and systematic uncertainties.

Reference	Model parameters	Pole parameters	Data
	$M_{{\rho }^{\prime }},{\Gamma }_{{\rho }^{\prime }}$ [MeV]	$M_{{\rho }^{\prime }}^{\mathrm{pole}},{\Gamma }_{{\rho }^{\prime }}^{\mathrm{pole}}$ [MeV]
ALEPH [52]	$1328\pm 15,468\pm 41$	$1268\pm 19,429\pm 31$	$\tau$
ALEPH [52]	$1409\pm 12,501\pm 37$	$1345\pm 15,459\pm 28$
Belle (fixed ${\left|F_V^{\pi }\left(0\right)\right|}^2$) [31]	$1446\pm 7\pm 28,434\pm 16\pm 60$	$1398\pm 8\pm 31,408\pm 13\pm 50$	$\tau$
Belle (all free) [31]	$1428\pm 15\pm 26,413\pm 12\pm 57$	$1384\pm 16\pm 29,390\pm 10\pm 48$	$\tau$
Dumm et. al. [7]	—	$1440\pm 80,320\pm 80$	$\tau$
Celis et. al. [8]	$1497\pm 7,785\pm 51$	$1278\pm 18,525\pm 16$	$\tau$
Bartos et. al. [51]	—	$1342\pm 47,492\pm 138$	$e^{+}{e}^{-}$
Bartos et. al. [51]	—	$1374\pm 11,341\pm 24$	$\tau$
This work	$1376\pm 6_{-73}^{+18},603\pm {22}_{-141}^{+236}$	$1289\pm 8_{-143}^{+52},540\pm {16}_{-111}^{+151}$	$\tau$
Reference	Model parameters	Pole parameters	Data
	$(M_{{\rho }^{\prime \prime }},{\Gamma }_{{\rho }^{\prime \prime }})$ [MeV]	$(M_{{\rho }^{\prime \prime }}^{\mathrm{pole}},{\Gamma }_{{\rho }^{\prime \prime }}^{\mathrm{pole}})$ [MeV]
ALEPH [52]	$=1713,=235$	$1700,232$	$\tau$
ALEPH [52]	$1740\pm 20,=235$	$1728\pm 20,232$
Belle (fixed ${\left|F_V^{\pi }\left(0\right)\right|}^2$) [31]	$1728\pm 17\pm 89,164\pm {21}_{-26}^{+89}$	$1722\pm 18,163\pm {21}_{-27}^{+88}$	$\tau$
Belle (all free) [31]	$1694\pm 41,135\pm {36}_{-26}^{+50}$	$1690\pm 94,134\pm {36}_{-28}^{+49}$	$\tau$
Dumm et. al. [7]	—	$1720\pm 90,180\pm 90$	$\tau$
Celis et. al. [8]	$1685\pm 30,800\pm 31$	$1494\pm 37,600\pm 17$	$\tau$
Bartos et. al. [51]	—	$1719\pm 65,490\pm 17$	$e^{+}{e}^{-}$
Bartos et. al. [51]	—	$1767\pm 52,415\pm 120$	$\tau$
This work	$1718\pm 4_{-94}^{+57},465\pm 9_{-53}^{+137}$	$1673\pm 4_{-125}^{+68},445\pm 8_{-49}^{+117}$	$\tau$

Table 5. Table 5: Low-energy observables of the pion vector form factor up to the quadratic term. Statistical and systematic uncertainties have been added in quadrature. Some of the values of the charged pion radius ⟨ r 2 ⟩ V π superscript subscript delimited-⟨⟩ superscript 𝑟 2 𝑉 𝜋 \langle r^{2}\rangle_{V}^{\pi} given in the table are not quoted in the original literature in units of GeV -2 but rather in fm 2 and the conversion has been evaluated by us.

Reference	${⟨r^2⟩}_V^{\pi }$ [GeV-2]	$c_V^{\pi }$ [GeV-4]
Colangelo et. al. [54]	$11.07\pm 0.66$	$3.2\pm 1.03$
Bijnens et. al. [32]	$11.22\pm 0.41$	$3.85\pm 0.60$
Pich et. al. [6]	$11.04\pm 0.30$	$3.79\pm 0.04$
Bijnens et. al. [33]	$11.61\pm 0.33$	$4.49\pm 0.28$
de Troconiz et. al. [55]	$11.10\pm 0.03$	$3.84\pm 0.02$
Masjuan et. al. [56]	$11.43\pm 0.19$	$3.30\pm 0.33$
Guo et. al. [57]	—	$4.00\pm 0.50$
Lattice [58]	$10.50\pm 1.12$	$3.22\pm 0.40$
Ananthanarayan et. al. [59]	$11.17\pm 0.53$	$[3.75,3.98]$
Ananthanarayan et. al. [60]	$[10.79,11.3]$	$[3.79,4.00]$
Schneider et. al. [47]	$10.6$	$3.84\pm 0.03$
Dumm et. al. [7]	$10.86\pm 0.14$	$3.84\pm 0.03$
Celis et. al. [8]	$11.30\pm 0.07$	$4.11\pm 0.09$
Ananthanarayan et. al. [61]	$11.10\pm 0.11$	—
Hanhart et. al. [62]	$11.34\pm 0.01\pm 0.01$	—
Colangelo et. al. [39]	$11.02\pm 0.10$	—
PDG [42]	$11.61\pm 0.28$	—
This work	$11.28\pm 0.08$	$3.94\pm 0.04$

Table 6. Table 6: Values for the subtraction constants calculated from the sum rule Eq. ( 2.22 ), as for the results of our fits given in Eq. ( 2.35 ), for three different values of s cut subscript 𝑠 cut s_{\rm{cut}} in the dispersive integral.

Sum rule	$s_{\mathrm{cut}}$ [GeV2]
	4	10	$\infty$	Fit Eq. (2.35)
${\alpha }_1$	$1.52$	$1.66$	$1.75$	$1.88\pm 0.01\pm 0.01$
${\alpha }_2$	$4.26$	$4.30$	$4.31$	$4.34\pm 0.01\pm 0.03$

Table 7. Table 7: Results for the fit to the BaBar τ − → K − ​ K S ​ ν τ → superscript 𝜏 superscript 𝐾 subscript 𝐾 𝑆 subscript 𝜈 𝜏 \tau^{-}\to K^{-}K_{S}\nu_{\tau} data [ 16 ] with a three-times-subtracted dispersion relation including two vector resonances in F V K ​ ( s ) superscript subscript 𝐹 𝑉 𝐾 𝑠 F_{V}^{K}(s) according to Eq. ( 3.9 ) with s cut = 4 subscript 𝑠 cut 4 s_{\rm{cut}}=4 GeV 2 in the dispersive integral.

Parameter	$s_{\mathrm{cut}}=4$ [GeV2]
	Fit $i)$	Fit $ii)$	Fit $iii)$	Fit $iv)$
${\tilde{\alpha }}_1$	$=1.88\left(1\right)$	$=1.84$	$=1.88\left(1\right)$	—
${\tilde{\alpha }}_2$	$=4.34\left(1\right)$	$=4.34$	$=4.34\left(1\right)$	—
$M_{{\rho }^{\prime }}$ [MeV]	$1467\left(24\right)$	$1538\left(32\right)$	$1489\left(25\right)$	1411(12)
${\Gamma }_{{\rho }^{\prime }}$ [MeV]	$415\left(48\right)$	$604\left(83\right)$	$297\left(36\right)$	394(35)
$\tilde{\gamma }$	$0.10\left(2\right)$	$0.36\left(11\right)$	$0.10\left(2\right)$	0.09(1)
${\tilde{\varphi }}_1$	$-1.19\left(16\right)$	$-1.48\left(13\right)$	$-1.10\left(15\right)$	$-1.88\left(9\right)$
${\chi }^2$/d.o.f.	$2.9$	$1.9$	$2.9$	3.3

Table 8. Table 8: Simultaneous fit results for different choices regarding resonance mixing and linear slopes parameters obtained with a three-times-subtracted dispersion relation including three vector resonances in F V π ​ ( s ) superscript subscript 𝐹 𝑉 𝜋 𝑠 F_{V}^{\pi}(s) and F V K ​ ( s ) superscript subscript 𝐹 𝑉 𝐾 𝑠 F_{V}^{K}(s) according to Eqs. ( 2.19 ) and ( 3.9 ) with s cut = 4 subscript 𝑠 cut 4 s_{\rm{cut}}=4 GeV 2 in the dispersive integral.

Parameter	$s_{\mathrm{cut}}=4$ [GeV2]
	Fit $a$	Fit $b$	Fit $c$
${\alpha }_1$	$1.88\left(1\right)$	$1.89\left(1\right)$	$1.87\left(1\right)$
${\alpha }_2$	$4.34\left(2\right)$	$4.31\left(2\right)$	$4.38\left(3\right)$
${\tilde{\alpha }}_1$	$={\alpha }_1$	$={\alpha }_1$	$1.88\left(24\right)$
${\tilde{\alpha }}_2$	$={\alpha }_2$	$={\alpha }_2$	$4.38\left(29\right)$
$m_{\rho }$ [MeV]	$=773.6\left(9\right)$	$=773.6\left(9\right)$	$=773.6\left(9\right)$
$M_{\rho }$ [MeV]	$=m_{\rho }$	$=m_{\rho }$	$=m_{\rho }$
$M_{{\rho }^{\prime }}$ [MeV]	$1396\left(19\right)$	$1453\left(19\right)$	$1406\left(61\right)$
${\Gamma }_{{\rho }^{\prime }}$ [MeV]	$507\left(31\right)$	$499\left(51\right)$	$524\left(149\right)$
$M_{{\rho }^{\prime \prime }}$[MeV]	$1724\left(41\right)$	$1712\left(32\right)$	$1746\left(1\right)$
${\Gamma }_{{\rho }^{\prime \prime }}$ [MeV]	$399\left(126\right)$	$284\left(72\right)$	$413\left(362\right)$
$\gamma$	$0.12\left(3\right)$	$0.15\left(3\right)$	$0.11\left(11\right)$
$\tilde{\gamma }$	$0.11\left(2\right)$	$=\gamma$	$0.11\left(5\right)$
${\varphi }_1$	$-0.23\left(26\right)$	$0.29\left(21\right)$	$-0.27\left(42\right)$
${\tilde{\varphi }}_1$	$-1.83\left(14\right)$	$-1.48\left(13\right)$	$-1.90\left(67\right)$
$\delta$	$-0.09\left(2\right)$	$-0.07\left(2\right)$	$-0.10\left(5\right)$
$\tilde{\delta }$	$=0$	$=0$	$-0.01\left(4\right)$
${\varphi }_2$	$-0.20\left(31\right)$	$0.27\left(29\right)$	$-1.15\left(71\right)$
${\tilde{\varphi }}_2$	$=0$	$=0$	$0.40\left(3\right)$
${\chi }^2$/d.o.f	$1.52$	$1.19$	$1.25$

Equations

$$⟨{\pi }^0{\pi }^{-}|\bar{d}{\gamma }^{\mu }u|0⟩=\sqrt{2}{\left(p_{{\pi }^{-}}-p_{{\pi }^0}\right)}^{\mu }{F}_V^{\pi }(s),$$

$$F_V^{\pi }(s){|}_{\mathrm{ChPT}}=1+\frac{2L_9^r(\mu )}{F_{\pi }^2}s-\frac{s}{96{\pi }^2{F}_{\pi }^2}\left[A_{\pi }(s,{\mu }^2)+\frac{1}{2}{A}_K(s,{\mu }^2)\right],$$

$$A_P(s,{\mu }^2)=\log \frac{m_P^2}{{\mu }^2}+\frac{8m_P^2}{s}-\frac{5}{3}+{\sigma }_P^3(s)\log \left(\frac{{\sigma }_P(s)+1}{{\sigma }_P(s)-1}\right),$$

$${\sigma }_P(s)=\sqrt{1-\frac{4m_P^2}{s}}.$$

$$F_V^{\pi }(s)=1+\frac{F_V{G}_V}{F_{\pi }^2}\frac{s}{M_{\rho }^2-s},$$

$$F_V{G}_V=F_{\pi }^2,$$

$$F_V^{\pi }(s)=\frac{M_{\rho }^2}{M_{\rho }^2-s}.$$

$$L_9^r(M_{\rho })=\frac{F_V{G}_V}{2M_{\rho }^2}=\frac{F_{\pi }^2}{2M_{\rho }^2}\backsimeq 7.2\cdot 10^{-3},$$

$$F_V^{\pi }(s)=\frac{M_{\rho }^2}{M_{\rho }^2-s}-\frac{s}{96{\pi }^2{F}_{\pi }^2}\left[A_{\pi }(s,{\mu }^2)+\frac{1}{2}{A}_K(s,{\mu }^2)\right],$$

$$F_{V}^{\pi}(s)=\frac{M_{\rho}^{2}}{M_{\rho}^{2}-s}\exp\Big{\{}-\frac{s}{96\pi^{2}F_{\pi}^{2}}\left[A_{\pi}(s,\mu^{2})+\frac{1}{2}A_{K}(s,\mu^{2})\right]\Big{\}}\,.$$

$${\Gamma }_{\rho }(s)$$

$$F_{V}^{\pi}(s)=\frac{M_{\rho}^{2}}{M_{\rho}^{2}-s-iM_{\rho}\Gamma_{\rho}(s)}\exp\Big{\{}-\frac{s}{96\pi^{2}F_{\pi}^{2}}{\rm{Re}}\left[A_{\pi}(s,\mu^{2})+\frac{1}{2}A_{K}(s,\mu^{2})\right]\Big{\}}\,,$$

$$F_V^{\pi }(s)$$

$${\Gamma }_{{\rho }^{\prime },{\rho }^{\prime \prime }}(s)$$

$$M_{\rho }=775.2(4)\mathrm{MeV},\gamma =0.15(4),{\varphi }_1=-0.36(24),$$

$$M_{{\rho }^{\prime }}=1438(39)\mathrm{MeV},{\mathrm{\Gamma }}_{{\rho }^{\prime }}=535(63)\mathrm{MeV},\delta =-0.12(4),{\varphi }_2,=-0.02(45),$$

$$M_{{\rho }^{\prime \prime }}=1754(91)\mathrm{MeV},{\mathrm{\Gamma }}_{{\rho }^{\prime \prime }}=412(102)\mathrm{MeV},$$

$$M_{\rho }^{\mathrm{pole}}=762.0(3)\mathrm{MeV},{\mathrm{\Gamma }}_{\rho }^{\mathrm{pole}}=143.0(2)\mathrm{MeV},$$

$$M_{{\rho }^{\prime }}^{\mathrm{pole}}=1366(38)\mathrm{MeV},{\mathrm{\Gamma }}_{{\rho }^{\prime }}^{\mathrm{pole}}=488(48)\mathrm{MeV},$$

$$M_{{\rho }^{\prime \prime }}^{\mathrm{pole}}=1718(82)\mathrm{MeV},{\mathrm{\Gamma }}_{{\rho }^{\prime \prime }}^{\mathrm{pole}}=397(88)\mathrm{MeV},$$

$$F_V^{\pi }(s)=\Omega (s)=\exp \left(\frac{s}{\pi }{\int }_{4m_{\pi }^2}^{\infty }d{s}^{\prime }\frac{{\delta }_1^1(s^{\prime })}{s^{\prime }\left(s^{\prime }-s\right)}\right).$$

$$\tan {\delta }_1^1(s)=\frac{\mathrm{Im}{F}_V^{\pi }(s)}{\mathrm{Re}{F}_V^{\pi }(s)},$$

$$F_V^{\pi }(s)=\exp \left[{\alpha }_{1}s+\frac{{\alpha }_2}{2}{s}^2+\frac{s^3}{\pi }{\int }_{4m_{\pi }^2}^{s_{\mathrm{cut}}}d{s}^{\prime }\frac{{\delta }_1^1(s^{\prime })}{(s^{\prime }{)}^3(s^{\prime }-s-i0)}\right],$$

$$F_V^{\pi }(s)=1+\frac{1}{6}⟨r^2{⟩}_V^{\pi }s+c_V^{\pi }{s}^2+d_V^{\pi }{s}^3+\cdots .$$

$$⟨r^2{⟩}_V^{\pi }=6{\alpha }_1,c_V^{\pi }=\frac{1}{2}\left({\alpha }_2+{\alpha }_1^2\right).$$

$${\alpha }_k=\frac{k!}{\pi }{\int }_{4m_{\pi }^2}^{s_{\mathrm{cut}}}d{s}^{\prime }\frac{{\delta }_1^1(s^{\prime })}{s^{\prime k+1}}.$$

$$d_V^{\pi }=\frac{1}{6}\left({\alpha }_3+3{\alpha }_1{\alpha }_2+{\alpha }_1^3\right).$$

$${\Gamma }_{{\rho }^{\prime },{\rho }^{\prime \prime }}(s)$$

$${\Gamma }_{{\rho }^{\prime }}(s)=$$

$${\sigma }_{\omega \pi }(s)=\frac{1}{s}\sqrt{(s-(m_{\omega }-m_{\pi }{)}^2)(s-(m_{\omega }+m_{\pi }{)}^2)},$$

$$F_V^{\pi }(s)=f(s)\Sigma (s),$$

$$f(s)=\exp \left[\frac{s}{\pi }{\int }_{4m_{\pi }^2}^{s_{\mathrm{cut}}}d{s}^{\prime }\frac{{\delta }_1^1(s^{\prime })}{s^{\prime }(s^{\prime }-s)}+\frac{s}{\pi }{\int }_{s_{\mathrm{cut}}}^{\infty }d{s}^{\prime }\frac{{\bar{\delta }}_1^1(s^{\prime })}{s^{\prime }(s^{\prime }-s)}\right].$$

Full text

1 Introduction

Tau lepton decays into a neutrino and hadrons provide a privileged scenario to investigate the non-perturbative regime of QCD under rather clean conditions since half of the transition is purely electroweak and can be computed straightforwardly. Such advantageous framework is used to improve our understanding of the hadronization of QCD currents as well as for determining the physical parameters, mass and width, of the intermediate resonances produced in the decay [1]. The strong dynamics is encoded within the hadronic matrix element which, in turn, is given in terms of form factors. As it is well-known, Chiral Perturbation Theory (ChPT) [2, 3] provides a successful description of these form factors valid at very low energies $E<<M_{\rho}$, where $M_{\rho}$ is the $\rho(770)$ resonance mass. However, as one approaches the resonance region $E\sim M_{\rho}$, ChPT ceases to provide a good description of the physics and the resonance fields shall be explicitly incorporated into the description as new degrees of freedom. This is the aim of Resonance Chiral Theory (RChT) [4], an effective field theory with resonance fields built in. RChT has been used by different groups as the initial setup approach to describe two meson tau decays providing a good description of the experimental measurements after being supplemented by arguments of analyticity and unitarity through dispersion relations. For example, the analyses of the $\pi^{-}\pi^{0}$ [5, 6, 7, 8] and $K\pi$ [9, 10, 11, 12] decay channels were found to be in a nice agreement with the rich experimental data provided by experiments.

One of the purposes of this work is to extend our series of dedicated analyses of two meson tau decays based on the framework of RChT supplemented by dispersion relations i.e. $\tau^{-}\to\pi^{-}\pi^{0}\nu_{\tau}$ [7], $\tau^{-}\to K_{S}\pi^{-}\nu_{\tau}$ and $\tau^{-}\to K^{-}\eta^{(\prime)}\nu_{\tau}$ [13, 14], and $\tau^{-}\to\pi^{-}\eta^{(\prime)}\nu_{\tau}$ [15], to the $K^{-}K_{S}$ final state meson system. The topic is of timely interest due to the recent measurement of the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ decay spectrum released by the BaBar collaboration [16]. This measurement is based on a sample of 223741 events and significantly improves the mass spectrum measured by CLEO in 1996 [17] were only 100 events in the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ were selected. The threshold for the decay $\tau\to K^{-}K_{S}\nu_{\tau}$ opens around 1000 MeV which is $\sim$100 MeV larger than $M_{\rho}+\Gamma_{\rho}$, a characteristic energy scale for the $\rho(770)$-dominance region. This implies that the $K^{-}K_{S}$ decay mode is not sensitive to the $\rho(770)$ peak, and consequently not useful to study its properties, but rather enhances its sensitivity to the properties of the heavier copies $\rho(1450)$ and $\rho(1700)$.

Also, a precise theoretical determination of the two-pion vector form factor within the SM (with a robust error band) is needed to increase the accuracy of the search for non-standard interactions in semileptonic weak charged currents [18, 19, 20]. This information is included as supplementary material of this paper (see Appendix A).

These facts motivate the present work where we intend to demonstrate that a reanalysis of the $\pi^{-}\pi^{0}$ and $K^{-}K_{S}$ decay spectra at present B-factories such Belle-II [21] could help improve notably the knowledge of the $\rho(1450)$ and $\rho(1700)$ resonance properties.

First, we reexamine the pion vector form factor focusing our effort on the improvement of the description of the energy region where the $\rho(1450)$ and $\rho(1700)$ come up into play. Our analysis is based on a three-times-subtracted dispersion relation. For the required input of the form factor phase entering the dispersive integral, we rely on Watson’s theorem [22]111Watson’s theorem applied to the pion vector form factor tell us that the form factor phase equals that of the two-pion scattering within the elastic region. and take advantage of the well-known parametrization of the $\pi\pi$ scattering phase shift existent in the literature, to drive the form factor phase up to 1 GeV. Above 1 GeV, we get a model for the phase from the exponential Omnès representation that we explicate in detail in section 2. This parametrization establishes the framework of our form factor description. However, we will also consider a number of variants to this approach that will be used to assess the (important) role of the systematic uncertainties in the extraction of the $\rho(1450)$ and $\rho(1700)$ resonance parameters that have been usually ignored or underestimated in the literature so far. Altogether results, to the best of our knowledge, in the most dedicated analysis of the intermediate- and high-energy region of the pion vector form factor experimental data to date. Also, the impact on the low-energy observables of the pion vector form factor is addressed and discussed as a byproduct of our approach.

Second, we built a parametrization for the kaon vector form factor in a similar fashion as for the pion one and perform a first analysis of the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ BaBar experimental measurement. The role of each participating resonance is discussed and the corresponding parameters are extracted from fits to data.

Finally, in view of the findings obtained from our analyses of the individual $\pi^{-}\pi^{0}$ and $K^{-}K_{S}$ channels, we perform a combined analysis to both data sets to see what can be learned. From our study, we can anticipate that although the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ BaBar data supersede the old CLEO data, still the precision is not sufficiently good enough to compete with the pion vector form factor modulus squared data from Belle. In all, we hope our study to be of interest for present and future experimental analysis of these decays.

This article is structured as follows. In section 2 we provide a bottom-up review of the pion vector form factor organized according to the fulfillment of unitarity and analyticity constraints. For illustrative purposes, we start with the ChPT calculation at order $\mathcal{O}{(p^{4})}$ and follow by the explicit inclusion of vector resonance states. For our study, we will consider three resonances i.e. $\rho(770),\rho(1450)$ and the $\rho(1700)$, and then submit the form factor to a unitarization procedure through the Omnès integral that leads to the Omnès exponentiation of the form factor. This parametrization allows us to get a model for the phase of the form factor valid up to the mass of the $\tau$. This phase is then inserted into a three-times-subtracted dispersion relation that completes our representation of the form factor, and the corresponding model parameters are fitted to the Belle measurement of the modulus squared of the pion vector form factor $|F_{V}^{\pi}(s)|^{2}$. Predictions and fits to the BaBar $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ decay spectrum measurement are discussed in section 3. In section 4, we perform joint fits to the decays $\tau^{-}\to\pi^{-}\pi^{0}\nu_{\tau}$ and $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ and finally, our conclusions are presented in section 5.

2 The pion vector form factor

The pion vector form factor has been measured in $e^{+}e^{-}\to\pi^{+}\pi^{-}$ [23, 24, 25, 26, 27, 28, 29] and in $\tau\to\pi^{-}\pi^{0}\nu_{\tau}$ [30, 31] and widely studied in the literature [5, 6, 7, 8, 32, 33, 34, 35, 36, 37, 38, 39] since it enters the description of many physical processes. As it is the main object concerning our analysis, we will thus provide in the following a brief, but detailed review of its calculation following a bottom-up approach according to the fulfillment of unitarity and analyticity constraints. We will start with the Chiral Perturbation Theory calculation at $\mathcal{O}(p^{4})$ to follow with the explicit inclusion of resonances as degrees of freedom. The $\pi\pi$ and $KK$ final-state interactions will be then resummed to all order through an Omnès exponentiation. Finally, a three-times-subtracted dispersion relation completes our form factor representation.

2.1 Exponential Omnès representation

The pion vector form factor $F_{V}^{\pi}(s)$ is defined through the hadronic matrix element of the vector current between the vacuum and the pion pair

[TABLE]

where $s=(p_{\pi^{-}}+p_{\pi^{0}})^{2}$.

At very low-energies, the pion vector form factor is well described by ChPT. The resulting calculation at $\mathcal{O}(p^{4})$ reads [3] (in the limit of exact isospin symmetry)

[TABLE]

where $L_{9}^{r}(\mu)$ is one of the renormalized low-energy couplings constants in the $\mathcal{O}(p^{4})$ chiral Lagrangian, and $A_{P}(s,\mu^{2})$ are two-pseudoscalars loop-functions (whose renormalization-scale dependence cancels out with the one in $L_{9}^{r}$) accounting for the unitary corrections given by

[TABLE]

where

[TABLE]

The validity of ChPT is restricted to very low energies, and as one approaches the region where the influence of new degrees of freedom, the lightest meson resonances, becomes important, ChPT ceases to provide a good description. Resonance Chiral Theory partly cures this limitation incorporating such resonances explicitly. For the case concerning us, the $\rho(770)$ resonance dominates the form factor. At leading order in powers of $1/N_{C}$, which it is $\mathcal{O}(p^{4})$ in the chiral expansion, the result is given by [4]

[TABLE]

where $F_{V}$ and $G_{V}$ measure, respectively, the strength of the $\rho V^{\mu}$ and $\rho\pi\pi$ couplings, with $V^{\mu}$ being the quark vector current. Assuming that the form factor vanishes when $s\to\infty$ one gets the condition

[TABLE]

that yields the usual Vector Meson Dominance (VMD) in the zero-width approximation

[TABLE]

Re-expanding Eq. (2.7) in $s$ and comparing with the polynomial part of its ChPT counterpart in Eq. (2.2) one gets an estimate for the $\mathcal{O}(p^{4})$ chiral coupling $L_{9}^{r}$

[TABLE]

which is in very good agreement with the value extracted from phenomenology. This result shows explicitly that the $\rho(770)$ contribution is indeed the dominant physical effect in the pion vector form factor.

An improved realization of the pion vector form factors stems from combining Eqs. (2.2) and (2.7), which yields:

[TABLE]

where the first term (VMD) is the dominant one in $1/N_{C}$ and resums an infinite number of local contributions in ChPT to all orders, while the second term includes the loop contributions that are next order in $1/N_{C}$.

In the spirit of Refs. [5, 9, 10], one can do better and perform a resummation of the $\pi\pi$ and $KK$ final-state interactions to all orders relying on unitarity and analyticity constraints. This leads to the Omnès exponentiation of the full loop function

[TABLE]

However, the previous expression still has an obvious defect, it cannot describe the energy region of the peak of the $\rho$-meson. For that, it is necessary to incorporate its width. The energy dependent width of the $\rho$ is related to the imaginary part of the loop function and is given by [40]

[TABLE]

In order to account for the $\rho$-resonance width, we insert $\Gamma_{\rho}(s)$ in the propagator of the $\rho$ in Eq. (2.10) arriving at:

[TABLE]

where, in order to avoid double counting of the imaginary part of the loop functions, only the real part of the loops is kept in the exponential. We would like to point out here that strict analyticity and unitarity is only maintained if the real part of the loop integral function is resummed in the propagator together with the imaginary part [11]. Up to $\mathcal{O}(p^{4})$ in the chiral expansion, resumming the real part in the propagator or in the exponential is fully equivalent and differences between the two approaches start to appear at $\mathcal{O}(p^{6})$. This effect, however, is seen to be numerically negligible in the $\pi\pi$ and $K\pi$ systems.

In all, equation (2.12) provides a suitable description of the $\rho$-dominance region. However, the precise measurement of the $\tau^{-}\to\pi^{-}\pi^{0}\nu_{\tau}$ decay spectrum by the Belle collaboration [31] indicates that heavier resonance excitations i.e. $\rho^{\prime}\equiv\rho(1450)\,,\rho^{\prime\prime}\equiv\rho(1700)$, contribute and cannot be simply neglected. We incorporate these higher excited resonances into the description in an analogous fashion as we have done for the $\rho$. The resulting expression takes the form:

[TABLE]

where the coefficients $\gamma$ and $\delta$ measure the relative weight between the contributions of the different resonances while the phases $\phi_{1}$ and $\phi_{2}$ stand for the corresponding interference, and with

[TABLE]

Notice that, for the energy-dependent width of the $\rho^{\prime}$ and $\rho^{\prime\prime}$ given in the previous equation, we have assumed that these resonances only decay into $\pi\pi$ since in our resonance chiral framework there is no warranty that other intermediate states contribute in the proportion given in Eq. (2.11) for the $\rho$-meson width. However, one can still explore a model including other intermediate states in a similar fashion and see what can be learned. We shall come back to this in section 2.2.

At this point, we would like to anticipate that the exponential Omnès representation of the pion vector form factor given by Eq. (2.13) will be used as input for the parametrization of the form factor phase entering the dispersive approach described in section 2.2 which, in turn, will be used to get the central results of this work. In particular, the extracted phase will be employed to describe the energy region from 1 GeV to $m_{\tau}$, and matched smoothly to the $\pi\pi$ scattering phase-shift solution of the Roy equations of [41] at 1 GeV. By doing this matching, sensitivity is lost to whether or not the real part of the loop function is resummed into the propagator denominator or kept into the exponential as discussed few lines above since the differing numerical results are tiny.

In the following, however, we would like first to prove this parametrization against experimental data as a warm-up. In total, we have nine unknown parameters, $\{M_{\rho},\gamma,\phi_{1},M_{\rho^{\prime}},\\ \Gamma_{\rho^{\prime}},\delta,\phi_{2},M_{\rho^{\prime\prime}},\Gamma_{\rho^{\prime\prime}}\}$, that can be inferred from fits to the measured modulus squared of the pion vector form factor extracted from the Belle $\tau^{-}\to\pi^{-}\pi^{0}\nu_{\tau}$ measurement [31]. The resulting fit parameters are222In all our fits throughout the paper, whenever the pion form factor is involved, we employ $F_{\pi}=92.316$ MeV, which is the central value using the restriction $\sqrt{2}|V_{ud}|F_{\pi}=(127.13\pm 0.02\pm 0.13)$ MeV from the 2018 PDG edition [42] with the 2006 PDG reported value $|V_{ud}|=0.97377\pm 0.00027$ used by Belle.

[TABLE]

with a $\chi^{2}$/d.o.f $=48.9/53\sim 0.92$, and where the associated uncertainty is the statistical error resulting from the fit. Notice, however, that $M_{\rho,\rho^{\prime},\rho^{\prime\prime}}$ and $\Gamma_{\rho^{\prime},\rho^{\prime\prime}}$ are model input parameters and do not correspond to the physical resonance mass and width. In order to extract the physical resonance pole parameters one should compute the pole position in the complex $s_{R}$ plane according to $\sqrt{s_{R}}=M_{R}-\frac{i}{2}\Gamma_{R}$. By doing so, the pole parameters associated to the $\rho$, $\rho^{\prime}$ and $\rho^{\prime\prime}$ resonances are then found to be

[TABLE]

where the uncertainties are calculated by assuming a Gaussian error propagation while simultaneously varying the corresponding unphysical masses and widths in Eq. (2.1).

The resulting form factor corresponding to our fit is displayed in Fig. 1 (solid green line) confronted to Belle data [31]. As can be seen from the plot and the $\chi^{2}$/d.o.f, the agreement with data is very satisfactory. In the figure, the ChPT calculation at $\mathcal{O}(p^{4})$ (cf. Eq. (2.2)) is also shown (dashed gray line) for illustrative purposes333For simplicity, in the figure we do not represent the ChPT calculation at $\mathcal{O}(p^{6})$ [32, 33] since including higher-order chiral corrections increases very little the energy region where the data is described well..

2.2 Dispersive representation

A Cauchy dispersion relation representation of the pion vector form factor is fully determined by the discontinuity across the cut along the positive real axis. Contributions to the discontinuity arise every time an intermediate state production threshold opens starting at $s_{\rm{th}}=4m_{\pi}^{2}$, the lightest possible contribution. The elastic approximation is confined to the two-pion contribution to the discontinuity and neglects heavier intermediate state contributions. In this limit, Watson’s theorem [22] states that the phase of the form factor equals that of the elastic $\pi\pi$ scattering phase and the form factor admits an analytic solution given in terms of the phase shift. Such a solution is the well-known Omnès equation [43] that, in terms of the $I=1$ $P$-wave $\pi\pi$ scattering phase shift $\delta_{1}^{1}(s)$ concerning us, reads

[TABLE]

The phase $\delta_{1}^{1}(s)$ entering the dispersive integral encodes the physics of the $\rho$-meson and it is known with an excellent precision in the elastic region $s\leq 1$ GeV2 from the solutions of Roy equations [41, 44] that are valid roughly up to $s_{0}=1.3$ GeV. However, uncertainties associated to the input of $\delta_{1}^{1}(s)$ can be estimated between the region of 1.3 GeV and the mass of the $\tau$ lepton $\sim 1.8$ GeV. The precise Belle measurement of the $\tau\to\pi^{-}\pi^{0}\nu_{\tau}$ decay distribution [31] indicates relevant contributions from the excited resonances $\rho^{\prime}$ and $\rho^{\prime\prime}$ that cannot be simply neglected. To include them, we adopt the form factor representation given in Eq. (2.13) to get a model for the form factor phase including the three participating resonances. This phase can be extracted from the relation

[TABLE]

that is only valid in the $\tau$ decay region $(s<m_{\tau}^{2})$. For the high-energy region $(s>m_{\tau}^{2})$ we guide smoothly the phase to $\pi$ [45] to ensure the correct asymptotic $1/s$ fall-off of the form factor [46].

The phase that we will use for our analysis consists in matching smoothly at 1 GeV the phase as extracted in Eq. (2.18) to the phase-shift solution of the Roy equations of Ref. [41]444Another successful parametrization of the phase shift is given in Ref. [44]. We only consider one of these parametrizations [41] since both agree rather well.. This procedure is inspired by, but slightly different than, the ones followed in Ref. [47] where a fit to data is first performed with the exponential Omnès representation of the pion vector form factor without kaon loops (cf. Eq. (2.13)) and then matched to the phase-shift solution of the Roy equations. For our fits, we follow the representation of the pion vector form factor outlined in Refs. [7, 8] and write a thrice-subtracted dispersion relation

[TABLE]

where $\alpha_{1}$ and $\alpha_{2}$ are two subtraction constants that can be related to chiral low-energy observables, namely the squared charged pion radius $\langle r^{2}\rangle_{V}^{\pi}$ and the coefficients of $\mathcal{O}(s^{2})$ and $\mathcal{O}(s^{3})$ terms in the chiral expansion, $c_{V}^{\pi}$ and $d_{V}^{\pi}$, respectively, appearing in the low-energy expansion of the form factor:

[TABLE]

Explicitly, the relations for the linear and quadratic slope parameters $\langle r^{2}\rangle_{V}^{\pi}$ and $c_{V}^{\pi}$ read

[TABLE]

These subtraction constants can be calculated theoretically through the sum rule

[TABLE]

For our analysis, however, we treat them as free parameters to be determined from fits to data. This has the advantage that they turn out to be less model dependent. Higher slope parameters can be computed from the previous sum rule. For example, the cubic slope parameters $d_{V}^{\pi}$ can be obtained through

[TABLE]

The use of a three-times dispersion relation in Eq. (2.19) makes the fit less sensitive to the higher-energy region of the dispersive integral where the phase is less well-known. In total, we have ten free parameters555The parametrization for the phase shift $\delta_{1}^{1}(s)$ of Ref. [41] contains a parameter for the $\rho$-meson mass, that we name $m_{\rho}$, that denotes the energy at which the phase shift passes through $\pi/2$ (and therefore it shall not be confused with real part of the pole of the $\rho$) and its quoted value is $m_{\rho}=773.6(9)$ MeV. For our study, in a first approximation we fix the model input parameter for the $\rho$-meson mass, $M_{\rho}$ in Eq. (2.13), to this value. However, will also test the sensitivity of our fits to this parameter by allowing it to float. entering $F_{V}^{\pi}(s)$ to be determined by a fit to the Belle data. Regarding the integral cutoff $s_{\rm{cut}}$, one should take a value as large as possible so as not to spoil the a priori infinite interval of integration and to avoid the effects of the spurious singularities generated due to the upper limit being finite, but low enough that the phase is well known within the interval. The parameters resulting from the fits are given in Table 1 as Fit 1 for four representative values of $s_{\rm{cut}}$, namely $m_{\tau}^{2}$ (third column), $4$ GeV2 (fourth column), $10$ GeV2 (fifth column) and finally the $s_{\rm{cut}}\to\infty$ limit (last column). The choice of $s_{\rm{cut}}=m_{\tau}^{2}$ is motivated by the fact that the model used to get the phase, Eq. (2.18), is only valid within the $\tau$ decay region $s\leq s_{\rm{cut}}\sim m_{\tau}^{2}$ and beyond that point the dispersive integral has no physical content. The resulting form factor corresponding to this cutoff generates, as mentioned above, a singularity at $s=m_{\tau}^{2}$ after bending the form factor shape in the preceding neighborhood region. As a consequence, the high-energy data points are not described well and, in turn, the values for the resonance parameters should be taken with great care. In fact, this fit is seen very sensitive to such singularities and the resulting parameters are found to be strongly correlated with unstable associated uncertainties. Therefore, we consider the results with $s_{\rm{cut}}=m_{\tau}^{2}$ only for illustrative purposes throughout the paper. Our reference fit corresponds to $s_{\rm{cut}}=4$ GeV2 (fourth column in Table 1) since this value of the cutoff deals well with the imbalance mentioned above. This does not mean that we know the phase shift up to that point of the integration interval but rather that the chosen cutoff is large enough that avoids, to large extent, the effects caused by the spurious singularity that arises at $s=s_{\rm{cut}}$. The results obtained by varying $s_{\rm{cut}}$ in Table 1 will be used to assess the systematic uncertainties of our fit results obtained with $s_{\rm{cut}}=4$ GeV2 666We would like to note the slightly low $\chi^{2}$/d.o.f. that in general, and in line with Ref. [8], we find along the fits of this section. This may indicate that there are too many free parameters to fit eventually, but as each of them has a physical meaning, it is reasonable to keep them all.. We shall return to a discussion on the integral cutoff below.

In Fig. 2, we show the resulting phase shifts for the chosen $s_{\rm{cut}}$. The phase shift solution of the Roy equations is given by the solid black curve while the variations due to $s_{\rm{cut}}$ are given by the dot-dashed blue $(s_{\rm{cut}}=m_{\tau}^{2}$), solid red $(s_{\rm{cut}}=4$ GeV${}^{2})$, dashed green $(s_{\rm{cut}}=10$ GeV${}^{2})$ and dotted black $(s_{\rm{cut}}\to\infty)$ curves, respectively. In the figure, the statistical uncertainty associated to our reference fit $(s_{\rm{cut}}=4$ GeV${}^{2})$ is also shown by the light red error band. Also, there are two brown dashed vertical lines shown in the figure. They are placed at 1 GeV and $\sqrt{s}=m_{\tau}$ and denote, respectively, the range where the phase shift from $\pi\pi$ scattering is used and the validity of the form factor phase shift parametrization as extracted through Eq. (2.18). In Fig. 3, we provide a graphical account of the resulting form factor for $s_{\rm{cut}}=4$ GeV2 (solid red curve and light red error band) and $s_{\rm{cut}}\to\infty$ (dotted black curve) compared to Belle data [31] . As can be seen from the figure and the corresponding $\chi^{2}$/d.o.f, excellent agreement with experimental data is seen with all data points.

In order to optimize the phase shift in the fit to the pion form factor measurement from $\tau^{-}\to\pi^{-}\pi^{0}\nu_{\tau}$ data, we have also run fits allowing to float the parameter for the $\rho$-meson mass entering $\delta_{1}^{1}(s)$. The corresponding fit results are contained in Table 1 as Fit 1-$\rho$ and the resulting parameters are found to be in general accordance with those of Fit 1.

For our reference fit, we have also probed the theoretical uncertainty associated to the choice of the matching point with the phase-shift solution of the Roy equations. The dependence of the fitted parameters on the matching point is explored through the fits that we collect in Table 2 as Fit I for a fixed $\rho$ mass parameter. The matching point values $0.85$ GeV (third column), $0.90$ GeV (fourth column), $0.95$ GeV (fifth column) and 1 GeV (last column) are used (the matching point $1$ GeV corresponds to our reference fit in the fourth column of Table 1 and is repeated here for ease of comparison). A look at the results of this table reveals that, in general, all the parameters are rather stable against variations of the matching point, although the $\rho^{\prime\prime}$ width is sensitive to these variations with a tendency of becoming larger as it decreases.

We shall now return to the discussion mentioned in section 2.1 about the inclusion of intermediate states other than $\pi\pi$ into the $\rho^{\prime}$ and $\rho^{\prime\prime}$ decay widths. We next allow such resonances to decay not only in $\pi\pi$ but also into $K\bar{K}$. In this case, the energy-dependent widths read

[TABLE]

and these are incorporated into the corresponding resonance propagators of Eq. (2.13), while the real part of the kaon loop is resumed in the exponential in a similar fashion as the pion ones. The corresponding fit results are gathered in Table 3 as Fit A (second column) and are compared to our reference fit (last column) again repeated here for ease of comparison. In this case, the $\chi^{2}$ increases by 12 and, as seen, there are two or three parameters that are affected by the fact of including kaons into the decay widths, the rest are seen rather stable upon comparison. The $\rho^{\prime}$ width decreases by $\sim 140$ MeV while the $\rho^{\prime\prime}$ mass increases up to $1775$ MeV. The $\rho^{\prime\prime}$ width is slightly shifted downwards but the associated error is enlarged.

We have also explored a variant of Eq. (2.24) that includes, moreover, the contribution $\rho^{\prime}\to\omega\pi$ into the $\rho^{\prime}$ decay width. In spirit of [48], we write the energy-dependent width as

[TABLE]

with

[TABLE]

and where $\mathcal{BR}(\rho^{\prime}\to(\pi\pi+KK))$ and $\mathcal{BR}(\rho^{\prime}\to\omega\pi)$ are relative branching ratios normalized with their sum equal to one. We use the PDG estimate $\mathcal{BR}(\rho^{\prime}\to\omega\pi)\equiv\Gamma(\rho^{\prime}\to\omega\pi)/\Gamma^{\rm{total}}_{\rho^{\prime}}\sim 0.21$ [42], which implies $\mathcal{BR}(\rho^{\prime}\to(\pi\pi+KK))\sim 0.79$, and the resulting fit results are presented in Table 3 as Fit B (third column). In this case, the $\rho^{\prime}$ mass(width) is shifted upwards(downwards) by $65(27)$ MeV with respect to our reference fit, while the $\rho^{\prime\prime}$ width is seen decreased by $116$ MeV. The other parameters remain rather stable. This last exercise serves to have an idea of the potential impact of the channel $\rho^{\prime}\to 4\pi$ in $\Gamma_{\rho^{\prime}}(s)$.

Finally, we also come back to the discussion on the integral cutoff $s_{\rm{cut}}$ and the corresponding generated singularities. We found that $s_{\rm{cut}}=4$ GeV2 is a suitable value for the integral cutoff in the dispersive integral. However, we would also like to consider a parametrization that allows both to cut the integral and avoid such singularities. In the case at hand, we have a parametrization for the form factor phase $\delta_{1}^{1}(s)$ which is valid up to $s=m_{\tau}^{2}$ and the idea is to extend it to the full region in an appropriate way. We follow Refs. [49, 50] and write

[TABLE]

where the function $f(s)$ is given by a once subtracted dispersion relation, that ensures $f(0)=1$, defined by

[TABLE]

The value of the phase $\bar{\delta}_{1}^{1}(s)$ in Eq. (2.28) should be such that it avoids the generation of spurious singularities and ensures the $1/s$ behavior of the form factor for $s\to\infty$. In order to fulfill these properties, we choose a smooth interpolation in $s$ for $\delta_{1}^{1}(s)$ above $s_{\rm{cut}}$ as simply as

[TABLE]

so that $\bar{\delta}_{1}^{1}(s_{\rm{cut}})=\delta_{1}^{1}(s_{\rm{cut}})$ and $\delta_{1}^{1}(s)\to\pi$ for large $s$ recovering the $1/s$ fall-off of $F_{V}^{\pi}(s)$. In this case, the integral going from $s_{\rm{cut}}$ to $\infty$ in Eq. (2.28) can be calculated explicitly and we arrive at [49, 50]

[TABLE]

Regarding the function $\Sigma(s)$, it contains the (inelastic) contributions beyond $s_{\rm{cut}}$ and may be understood as giving the correction to the linear continuation of the phase $\delta_{1}^{1}(s)$ above $s_{\rm{cut}}$ as we did in Eq. (2.29). It is described by an analytical function on the $s$-plane with a cut from $s_{\rm{cut}}$ to $\infty$ and should be obtained from a model or fitted to experiment since it is largely unknown. Often, it is parametrized by a conformal transformation that maps the right-hand cut in the complex $s$-plane into the unit circle through

[TABLE]

with the variable $\omega(s)$ given by

[TABLE]

For our analysis, we take the condition $a_{0}=1$ to ensure $F_{V}^{\pi}(0)=1$.

We next probe the application of Eq. (2.30) against data for $s_{\rm{cut}}=4$ GeV2 with one and two parameters in the expansion of Eq. (2.31). The resulting fit parameters are found to be

[TABLE]

with a $\chi^{2}$/d.o.f $=32.3/53\sim 0.61$ for the one-parameter fit, and

[TABLE]

with a $\chi^{2}$/d.o.f $=35.6/52\sim 0.63$ for the two-parameter fit. The large uncertainty associated to $a_{2}$ suggests not continuing the expansion to higher orders.

We are now in the position to combine all the results from the different fits that we have obtained from our dedicated analysis discussed above and that we graphically compare in Figs. 4 and 5 for the form factor phase shift and modulus squared, respectively. In particular, we show Fit 1 (reference fit) and Fit 1-$\rho$ with $s_{\rm{cut}}=4$ GeV2 from Table 1, Fit I at the matching point of 0.85 GeV from Table 2, Fit A of Table 3 and the fit results of Eq. (2.34) named as Fit singularities. In the figures, the statistical uncertainty associated to Fit 1 (reference fit) is also displayed by the light red error band. Tables with the corresponding numerical values including both statistical and systematic errors are given as ancillary material (see Appendix A). As seen from these figures, both the phase and the form factor absolute value squared are rather stable and only small differences are seen in the dip region $\sim 2.5$ GeV2 caused by the destructive interference between the $\rho^{\prime}$ and $\rho^{\prime\prime}$ resonances. To present our central results, we quote the values of our reference fit in Table 1 with $s_{\rm{cut}}=4$ GeV2 (Fit 1) and ascribe a conservative systematic uncertainty coming from the largest variations of central values with respect to the differing results shown in Tables 1 and 2 while changing $s_{\rm{cut}}$ and the matching point with the Roy equations, respectively, in Table 3 due to the inclusion of the $K\bar{K}$ decay channel into the $\rho^{\prime}$ and $\rho^{\prime\prime}$ widths and in Eq. (2.34) due to the parametrization presented in Eq. (2.30) that avoids the aforementioned singularities. We then obtain

[TABLE]

where the first uncertainty is the statistical fit error while the second is our estimated systematic uncertainty.

As has been already stated in section 2.1, the resonance mass and width parameters of Eq. (2.35) are unphysical fit parameters. To obtain the physical resonance mass and width, we calculate the pole positions in the complex $s$-plane. This yields:

[TABLE]

where the systematic uncertainties are calculated by assuming a Gaussian error propagation while simultaneously varying the corresponding unphysical mass and width given in Eq. (2.35). The results given in Eq. (2.36) constitute one of the fundamental results of the article, we show that the extraction of the pole mass and width of the $\rho^{\prime}$ and $\rho^{\prime\prime}$ resonances is limited by theoretical errors that, as we will see in the following, have been usually ignored or underestimated in the literature.

In Table 4, we show our results for the $\rho(1450)$ and $\rho(1700)$ resonance parameters compared to other determinations quoted in the literature. We consider the pole mass and width as the relevant resonance properties since one expects the pole parameters to be essentially model independent. To shed further light on the comparison, we have calculated the pole mass and width of the $\rho(1450)$ and $\rho(1700)$ parameters appearing in phenomenological amplitudes where resonances are introduced through Breit-Wigner (BW) type functions, thus being strongly model dependent. This is the case of the Gounaris and Sakurai parametrization [53] used in the fits by ALEPH [52] and Belle [31], and also in [8] where a BW function supplemented by a dispersion relation was used to fit Belle data. For the $\rho(1450)$, we obtain a pole mass(width) on the lower(upper) side, but in general agreement with respect to previous determinations. Our results are found to be in a remarkable agreement with [8], especially after the pole position is computed, while they are seen in a mild tension with respect to the PDG educated guess reported values $M_{\rho^{\prime}}^{\rm{PDG}}=1450\pm 25$ MeV and $\Gamma_{\rho^{\prime}}^{\rm{PDG}}=400\pm 60$ MeV [42]. For the $\rho(1700)$ we obtain, on the one hand, a pole mass slightly lower than, but in agreement with, previous determinations and the PDG reported value $M_{\rho^{\prime\prime}}^{\rm{PDG}}=1720\pm 20$ MeV, and $\sim 300$ MeV higher than [8]. On the other hand, as seen from the results quoted in the table, the values for the $\rho(1700)$ width show some scatter. Our value is found to be 250-300 MeV higher than the Belle value and than the findings of [7], $\sim$200 MeV higher than the PDG value $\Gamma_{\rho^{\prime}}^{\rm{PDG}}=250\pm 100$ MeV, in line with [51] and $\sim$150 MeV lower than [8]. Due to the large uncertainties associated to the $\rho(1450)$ and $\rho(1700)$ resonances parameters, however, we cannot draw any more definite conclusion. For that, more precise data in the $\rho(1450)$ and $\rho(1700)$ region would be highly desirable.

Regarding the chiral observables associated to the low-energy expansion of the form factor (cf. Eqs. (2.20) and (2.21)), taking into account the results quoted in Eq. (2.35), we obtain the values shown in the last row of Table 5 where we have added the systematic error in quadrature to the statistical uncertainty. In this table, we also display previous determinations of these quantities for comparison.

As can be seen, our results are found to be in good agreement with, but in general more precise than, all previous determinations.

It is opportune to mention again that we have treated the subtraction constants $\alpha_{1}$ and $\alpha_{2}$ as free parameters that capture our ignorance of the higher energy part of the integral. However, in order to check the consistency of our approach, we have also calculated these constants through the sum rule given in Eq. (2.22), as for the central values of our analysis presented in Eq. (2.35), for three different values of $s_{\rm{cut}}$ i.e. 4 GeV2, 10 GeV2 and $\infty$. The values we get are collected in Table 6 and are seen in a reasonable good agreement with the results of our fits that we show in the last column for ease of comparison. In other words, this tell us that the content of the phase is such that saturates rather well the dispersive integral, otherwise the differing results between the sum rules and the fitted subtraction constants would be larger.

The next shape parameter in the expansion, the cubic slope $d_{V}^{\pi}$, is much less known. To the best of our knowledge, there are neither theoretical results from ChPT nor calculations on the Lattice. We obtain

[TABLE]

a value which is seen slightly larger than previous estimates $d_{V}^{\pi}=9.70\pm 0.40$ GeV*-6* [63], $d_{V}^{\pi}=9.84\pm 0.05$ GeV*-6* [7] and $d_{V}^{\pi}=[10.14,10.56]$ GeV*-6* [60].

3 Predictions and fits to $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ BaBar data

The theoretical expression for the differential decay distribution for the transition $\tau^{-}\to K^{-}K^{0}\nu_{\tau}$ in terms of the $K^{-}K^{0}$ invariant mass can be written as [35, 64, 65]

[TABLE]

and it is related to the normalized invariant mass spectrum through

[TABLE]

where $N_{\rm{events}}$ is the total number of measured events, $\Gamma_{\tau}$ is the inverse $\tau$ lifetime and $\Delta^{\rm{bin}}_{K^{-}K_{S}}$ is the bin width. $\bar{B}(\tau^{-}\to K^{-}K_{S}\nu_{\tau})\equiv\bar{B}$ is a normalization constant that, for a perfect description of the spectrum, would equal the branching ratio. For our analysis, we fix this normalization to the BaBar measured branching fraction $\bar{B}=0.739(11)_{\rm{stat}}(20)_{\rm{syst}}\times 10^{-3}$ [16] 777Another possibility would be to let this constant float and infer its value from fits to the data [11, 12, 13]. However, in order to reduce the number of free parameters to fit, we prefer to fix this constant to the branching ratio measured by BaBar.. The factor $1/2$ is due to the $K^{-}K_{S}$ decay channel is analyzed. The corresponding number of events measured by BaBar is $223741\pm 3461$ [16] and the bin width is $0.04$ GeV.

$F_{V}^{K}(s)$ in Eq. (3.1) denotes the participant kaon vector form factor that we will describe in the following. Similar to Eq. (2.1), it can be defined via the matrix element of the vector current between the vacuum and the $K^{-}K^{0}$ pair as

[TABLE]

where, as in the case of the pion vector form factor, the two-kaon final state corresponds to a $I=J=1$ configuration. In order to obtain the expression for $F_{V}^{K}(s)$ at $\mathcal{O}(p^{4})$ in ChPT we need the expressions of the $K^{+}K^{-}$ and $K^{0}\bar{K}^{0}$ form factors. These can be found in the literature and read [3]

[TABLE]

The $I=1$ component corresponding to the $K^{-}K^{0}$ state can be extracted from Eqs. (3.4) and (3.5) and yields

[TABLE]

We would like to note that the ChPT calculation at $\mathcal{O}(p^{4})$ of the pion (cf. Eq. (2.2)) and kaon (cf. Eq. (3.6)) vector form factors are the same. For our study, we will consider exact $SU(3)$ flavor symmetry and assume that, in a first approximation, both form factors are also the same at energies higher than the chiral region. This allows us to predict the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ decay spectrum from the description of the pion vector form factor carried out in the previous section.

In Fig. 6, we show such a prediction (dotted red curve) based on our central results of the pion vector form factor analysis presented in Eq. (2.35), confronted to the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ spectrum measured by BaBar (black circles). A look at the figure shows a clear disagreement between our prediction and the BaBar data. The shape of the decay spectrum is not followed by this approach, as also indicated by the low value of the corresponding branching ratio, ${\rm{BR}}(\tau^{-}\to K^{-}K_{S}\nu_{\tau})=0.545(32)\times 10^{-3}$, which is seen $\sim 5\sigma$ away the BaBar measurement ${\rm{BR}}(\tau^{-}\to K^{-}K_{S}\nu_{\tau})=0.739(11)_{\rm{stat}}(20)_{\rm{syst}}\times 10^{-3}$. From these results, we conclude that using the pion vector form factor to describe the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ decay channel is a too rough approach. This is due to the $K^{-}K^{0}$ production threshold is around $1$ GeV, and thus one should expect the $\rho^{\prime}$ and $\rho^{\prime\prime}$ vector resonances to play a significantly different role than in $\tau^{-}\to\pi^{-}\pi^{0}\nu_{\tau}$ as also noted in [66, 67]. Thus, the weight of each resonance contribution, represented by the coefficients $\gamma$ and $\delta$ and the phases $\phi_{1}$ and $\phi_{2}$ in Eq. (2.13), can vary with respect to those entering the form factor of the pion. We will check this in the following by performing individual fits to the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ spectrum. To this end, we adapt the Omnès representation of the pion vector form factor (cf. Eq. (2.13)) to the kaon vector form factor ones and write:

[TABLE]

where the coefficients $\tilde{\gamma}$ and $\tilde{\delta}$ and the phases $\tilde{\phi}_{1}$ and $\tilde{\phi}_{2}$ account, respectively, for the relative importance between the contributions of the different resonances and the corresponding interference in the $K^{-}K^{0}$ system.

From the kaon vector form factor in Eq. (3.7), we extract its phase through

[TABLE]

and this is inserted888The $\delta_{1}^{KK}$ phase as extracted from Eq. (3.8) is also matched to the $\pi\pi$ scattering at 1 GeV as explained along the lines of section 2. into a three-times-subtracted dispersive representation of the form factor

[TABLE]

where $\tilde{\alpha}_{1}$ and $\tilde{\alpha}_{2}$ are two subtraction constants corresponding to the slope and curvature of the form factor of the kaon.

We have considered different fits to the measured $m_{K^{-}K_{S}}$ invariant mass distribution999We would like to notice here that the last two data points of the BaBar paper’s Table II [16] have been rescaled to match the paper’s figure 12. and found on the one hand that, in full generality, the data is not sensitive either to the low-energy region or to the $\rho(770)$ peak region. This is expected due to the $K^{-}K_{S}$ production threshold opens around 1000 MeV which is some 100 MeV larger than $M_{\rho}+\Gamma_{\rho}$, the energy region more influenced by the presence of the $\rho(770)$. This implies first that the slopes of the kaon vector form factor, which encode the physics immediately above threshold, cannot be fitted with $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ data and second that the fits lead unrealistic parameters under floating the $\rho$-meson parameters. On the other hand, the data is scarce in the $\rho(1700)$ resonance region and thus not suitable for extracting the corresponding resonance parameters. We have therefore fixed the slopes associated to the kaon vector form factor and the $\rho$-meson mass ($773.6(9)$ MeV), and considered fits varying only the $\rho(1450)$-resonance mass and width, and $\tilde{\gamma}$ and $\tilde{\phi}_{1}$, while neglecting the contribution of the $\rho(1700)$ resonance to the decay i.e. $\tilde{\delta}_{1}=0$. In Table 7, we show the results of our fits using different settings. Fit $i)$ corresponds to fixing the slopes $\tilde{\alpha}_{1}$ and $\tilde{\alpha}_{2}$, to the slopes $\alpha_{1}$ and $\alpha_{2}$ given in Eq. (2.35) obtained from the analysis of the pion vector form factor, while Fits $ii)$ and $iii)$ are variants of it. In particular, Fit $ii)$ is the result of fixing the slope to $\tilde{\alpha}_{1}=1.84(30)$ obtained from the kaon radius of Ref. [33] while Fit $iii)$ includes the intermediate $K\bar{K}$ state into the $\rho^{\prime}$ decay width (cf. Eq. (2.24)). Finally, Fit $iv)$ is the result of the direct application of the exponential vector form factor in Eq. (3.7) to fit experimental data.

The fit with best $\chi^{2}$/d.o.f is seen for Fit $ii)$ but the resulting fit parameters carry the larger error, and the (large) uncertainty associated to $\tilde{\alpha}_{1}$ has not been taken into account. Because of that, the results from Fit $ii)$ should be taken with a word of caution and we only consider them as illustration of the potential effects due to the low-energy parameters of the kaon form factor. The values appearing in Table 7 can be translated to pole values along the lines discussed in the previous section. This yields $M^{\rm{pole}}_{\rho^{\prime}}\,=\,1422\pm 22$ MeV and $\Gamma^{\rm{pole}}_{\rho^{\prime}}\,=\,393\pm 41$ MeV (Fit $i)$), $M^{\rm{pole}}_{\rho^{\prime}}\,=\,1453\pm 29$ MeV and $\Gamma^{\rm{pole}}_{\rho^{\prime}}\,=\,546\pm 70$ MeV (Fit $ii)$), and $M^{\rm{pole}}_{\rho^{\prime}}\,=\,1466\pm 23$ MeV and $\Gamma^{\rm{pole}}_{\rho^{\prime}}\,=\,289\pm 33$ MeV (Fit $iii)$) for the dispersive approaches, and $M^{\rm{pole}}_{\rho^{\prime}}\,=\,1370\pm 15$ MeV and $\Gamma^{\rm{pole}}_{\rho^{\prime}}\,=\,373\pm 30$ MeV for the exponential representation (Fit $iv)$). From these results we conclude that while the pole mass of the $\rho(1450)$ resonance as extracted from the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ decay tends to be larger than the values obtained in the previous section from the analysis of the pion vector form factor, the width tends to be smaller and the associated fit uncertainties are in both cases larger. Also, we would like to note that the relative weight $\tilde{\gamma}$ as extracted from the $K^{-}K_{S}$ channel is found to be in accordance with the values determined in the previous section.

As a matter of example, in Fig. 6 we provide a graphical account of the dispersive Fit $i)$ (solid green curve) and of the exponential Fit $iv)$ (dashed blue curve). Notice that, as occurs in [68], the second and third data points are difficult to accommodate in any case. These results show that, although the fits to the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ decay spectrum have considerably improved with respect to predictions discussed above (dotted red curve) and seem to agree rather well with BaBar data, yielding BR$(\tau^{-}\to K^{-}K_{S}\nu_{\tau})=0.749(93)\times 10^{-3}$ (Fit $i)$) and BR$(\tau^{-}\to K^{-}K_{S}\nu_{\tau})=0.744(89)\times 10^{-3}$ (Fit $iv)$), the quality of the fit as indicated by the $\chi^{2}$/d.o.f is not satisfactory enough. This fact motivates the combined analysis, detailed in the next section, of the Belle data of the pion vector form factor modulus squared $|F_{V}^{\pi}|^{2}$ and the BaBar data of the decay $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$. Such analysis shall allow us to determine the $\rho(1450)$ and $\rho(1700)$ resonance parameters with improved precision and obtain a good description of the measured $K^{-}K_{S}$ decay spectrum.

4 Joint fits to $|F_{V}^{\pi}|^{2}$ Belle and $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ BaBar data

The $\chi^{2}$ minimised in our simultaneous fit is

[TABLE]

where the first and second terms correspond, respectively, to the Belle pion vector form factor data [31] and to the BaBar $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ measurement [16]. For the later, $\mathcal{N}_{j}^{\rm{exp}}$ and $\sigma_{\mathcal{N}_{j}^{\rm{exp}}}$ are the experimental normalized number of events for $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ and the associated uncertainties in the $i$-th bin, respectively.

The parameters entering the dispersive representation of the form factors of Eqs. (2.19) and (3.9) are therefore determined by a simultaneous fit to both data sets and include:

• •

The two subtraction constants, $\alpha_{1,2}$ and $\tilde{\alpha}_{1,2}$, corresponding to the slope and curvature parameters associated to the low-energy expansion of the pion and kaon form factors.

• •

The masses and decay widths of the participating $\rho^{\prime}$ and $\rho^{\prime\prime}$ resonances, $M_{\rho^{\prime},\rho^{\prime\prime}}$ and $\Gamma_{\rho^{\prime},\rho^{\prime\prime}}$, used to model the phase entering the dispersive integral. The parameter for the $\rho$-meson mass, $M_{\rho}$, is taken equal to that entering the phase shift, $m_{\rho}$, and its quoted value is fixed to $773.6(9)$ MeV as discussed in section 2.

• •

The resonance mixing parameters, $\gamma,\delta$ and $\tilde{\gamma},\tilde{\delta}$, and their phases, $\phi_{1,2}$ and $\tilde{\phi}_{1,2}$.

In Table 8, we show the results of our simultaneous fits using slightly different settings, though in all of them a three-times-subtracted dispersion relation according to Eqs. (2.19) and (3.9) with $s_{\rm{cut}}=4$ GeV2 in the dispersive integral is employed. Fit $a$ (second column) corresponds to fixing $\tilde{\alpha_{1}}=\alpha_{1}$ and $\tilde{\alpha_{2}}=\alpha_{2}$ and taking $\tilde{\delta}=0$ because, as we have discussed in section 3, the BaBar measurement of the decay $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ is still not sensitive to $\rho(1700)$ resonance properties. The corresponding fit results supports the relative weights $\gamma$ and $\tilde{\gamma}$ to be the same for the $\pi^{-}\pi^{0}$ and $K^{-}K_{S}$ channels. This feature is proven in Fit $b$ (third column) by enforcing $\gamma=\tilde{\gamma}$. By doing this, the $\chi^{2}$/d.o.f is reduced from $1.52$ to $1.19$ and the values of the fitted parameters remain basically the same but for the $\rho^{\prime\prime}$-width, whose central value is shifted downwards by $\sim 100$ MeV, and to less extent for the $\rho^{\prime}$-mass, which suffers a variation of $\sim 50$ MeV upwards, but still compatible within errors. Finally, Fit $c$ (last column) is the result of letting all parameters to float independently and the corresponding fit parameters are found to be compatible with Fits $a$ and $b$, though with larger uncertainties. This fit also yields results that supports the assumption $\tilde{\alpha_{1}}=\alpha_{1}$ and $\tilde{\alpha_{2}}=\alpha_{2}$ made in Fits $a$ and $b$. As a side result, we extract the charge kaon radius through $\langle r^{2}\rangle_{V}^{K}=6\tilde{\alpha}_{1}$ (cf. Eq. (2.21)). Our value, $\langle r^{2}\rangle_{V}^{K}=(11.28\pm 1.44)$ GeV*-2*, lies in the ballpark of results for this quantity $\langle r^{2}\rangle_{V}^{K}=(9.09\pm 1.82)$ GeV*-2* and $\langle r^{2}\rangle_{V}^{K}=(11.07\pm 1.82)$ GeV*-2* [33], $\langle r^{2}\rangle_{V}^{K}=(9.76\pm 0.85)$ GeV*-2* [69] and $\langle r^{2}\rangle_{V}^{K}=[10.02,10.79]$ GeV*-2* [70].

The results of our fits are confronted to the Belle $|F_{V}^{\pi}|^{2}$ form factor measurement and to the BaBar $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ distribution in Figs. 7 and 8, respectively. Satisfactory agreement with the experimental data is seen in accord with the observed $\chi^{2}$/d.o.f. However, the potential of a combined analysis of the $|F_{V}^{\pi}|^{2}$ and $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ data is presently limited by the fact that the errors associated to the latter are still relatively large and the BaBar measurement of the $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ spectrum is yet not sensitive to the $\rho(1700)$ resonance properties. This presents a limitation in determining the $\rho(1450)$ and $\rho(1700)$ resonance parameters with improved precision with respect to the individual analysis of $|F_{V}^{\pi}|^{2}$. Because of that, we postpone a dedicated study of the systematic uncertainties as we did in section 2 for the future, when new and more precise measurements become available.

5 Conclusions

An ideal roadmap for describing meson form factors would require a model-independent approach demanding a full knowledge of QCD in both its perturbative and non-perturbative regimes, knowledge not yet unraveled. An alternative to such enterprise would pursuit a synergy between the formal theoretical calculations and experimental data. In this respect, dispersion relations are a powerful tool to direct oneself towards a model-independent description of form factors. In this paper, we have revisited the pion vector form factor as extracted from $\tau\to\pi^{-}\pi^{0}\nu_{\tau}$, and provided a parametrization for the kaon vector form factor that describes the decay $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$, by exploiting the synergy between dispersion relations and Chiral Perturbation Theory.

The pion vector form factor is a classic object in low-energy QCD that provides a privileged benchmark to study the effects of $\pi\pi$ interaction under rather clean conditions. These pion-pion interactions are universal and enter the description of many physical observables, hence the importance of having good control of them. For our analysis, we have used a three-times-subtractred dispersion relation and exploited Watson’s theorem, and the fact that the elastic $P$-wave $\pi\pi$ interactions capturing the physics of the $\rho$-resonance are encoded in the phase shift, to drive the form factor phase entering the dispersive integral up to 1 GeV from the well-known parametrization of the $\pi\pi$ scattering phase shift existent in the literature. Above 1 GeV, the $\rho(1450)$ and $\rho(1700)$ resonance effects show up and to obtain an improved description of the energy region where these resonances come up into play is one of the purposes of this work. For that, we have used a model for the phase extracted from the exponential Omnès representation of the form factor (cf. Eq. (2.13)), whose direct application to the pion vector form factor experimental data is seen very satisfactory (see Fig. 1), that we match smoothly at 1 GeV to the $\pi\pi$ scattering phase. Armed with this parametrization, we have carried out a very dedicated analysis of the high-statistics Belle experimental data and assessed the role of the theoretical systematic uncertainties in the determination of the $\rho(1450)$ and $\rho(1700)$ physical resonance parameters by considering a number of variants of it. Tables with the corresponding numerical values including both statistical and systematic errors are given as ancillary material of this paper (see Appendix A). From our study, we conclude that the determination of the pole mass and width of these resonances (cf. Eq. (2.36)) is limited by theoretical errors that have been usually ignored or underestimated in the literature so far.

On a second stage, we have performed a first analysis of the recent BaBar measurement experimental data on $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ based on a parametrization of the participant kaon vector form factor that is built in a similar fashion to that of the pion. We have shown that while the production threshold of this decay channel sits around 1000 MeV and therefore it is out of the $\rho(770)$-dominance region, the role of the $\rho(1450)$-resonance is different than in the $\pi^{-}\pi^{0}$ mode and indeed dominates. As a result of our fits (see Table 7 and Fig. 6), we have extracted its associated pole parameters. Regarding the $\rho(1700)$, the data is scarce in this region and thus it is not-yet suitable for extracting the corresponding resonance parameters.

Finally, we have pointed out that high-quality data on the decay $\tau^{-}\to K^{-}K_{S}\nu_{\tau}$ will allow one to determine the $\rho(1450)$ and $\rho(1700)$ resonance pole parameters with improved precision from a combined analysis with the pion vector form factor. In summary, we hope our analysis to be of interest for present and future $Z$, tau-charm and $B$-factories where new measurements should be possible.

Acknowledgements

The authors acknowledge fruitful discussions with Feng-Kun Guo and Bing-Son Zou. The work of S.G-S has been supported in part by the CAS President’s International Fellowship Initiative for Young International Scientists (Grant No. 2018DM0034), by the Sino-German Collaborative Research Center “Symmetries and the Emergence of Structure in QCD” (NSFC Grant No. 11621131001, DFG Grant No. TRR110), by NSFC (Grant No. 11747601) and by the National Science Foundation (PHY-1714253). P. R. acknowledges support from Conacyt through projects FOINS-296-2016 (Fronteras de la Ciencia), 250628 (Ciencia Básica) and Fondo SEP-Cinvestav 2018.

Appendix A Supplementary material

As supplementary material of this paper, we provide our fit results shown in Figs. 4 and 5 for the pion form factor phase and modulus squared in two data files named FormFactorPhase.dat and FormFactorModulusSquared.dat, respectively.

In the FormFactorPhase.dat file, the first column is $\sqrt{s}$ (in GeV) and the other two columns correspond to the numerical value of the phase and the statistical error for each fit, and they are given until 2 GeV. For our reference fit (Fit 1), moreover, we provide a fourth column with our conservative estimate for the systematic error. We calculate this after symmetrizing the asymmetric systematic uncertainties associated to the fit parameters presented in Eq. (2.35) and by assuming a Gaussian error propagation while simultaneously varying the corresponding parameters.

In the FormFactorModulusSquared.dat file, the first column is $s$ (in GeV2) and the second column is the value for the form factor modulus squared for each fit, and is given until $s=m_{\tau}^{2}$. The associated upper and lower statistical error bands are collected in the third and fourth columns of the file, respectively. As before, for our reference fit (Fit 1), we also provide an estimated systematic error, and the resulting upper and lower bands are gathered, respectively, in columns fifth and sixth of the file.