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

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.


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 ρ , where M ρ is the ρ(770) resonance mass. However, as one approaches the resonance region E ∼ M ρ , 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 π − π 0 [5,6,7,8] and Kπ [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. τ − → π − π 0 ν τ [7], τ − → K S π − ν τ and τ − → K − η ( ) ν τ [13,14], and τ − → π − η ( ) ν τ [15], to the K − K S final state meson system. The topic is of timely interest due to the recent measurement of the τ − → K − K S ν τ 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 τ − → K − K S ν τ were selected. The threshold for the decay τ → K − K S ν τ opens around 1000 MeV which is ∼100 MeV larger than M ρ + Γ ρ , a characteristic energy scale for the ρ(770)-dominance region. This implies that the K − K S decay mode is not sensitive to the ρ(770) peak, and consequently not useful to study its properties, but rather enhances its sensitivity to the properties of the heavier copies ρ(1450) and ρ(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 nonstandard 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 π − π 0 and K − K S decay spectra at present B-factories such Belle-II [21] could help improve notably the knowledge of the ρ(1450) and ρ(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 ρ(1450) and ρ(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] 1 and take advantage of the well-known parametrization of the ππ 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 ρ(1450) and ρ(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 τ − → K − K S ν τ 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 π − π 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 τ − → K − K S ν τ 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 O(p 4 ) and follow by the explicit inclusion of vector resonance states. For our study, we will consider three resonances i.e. ρ(770), ρ(1450) and the ρ(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 τ . 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 (s)| 2 . Predictions and fits to the BaBar τ − → K − K S ν τ decay spectrum measurement are discussed in section 3. In section 4, we perform joint fits to the decays τ − → π − π 0 ν τ and τ − → K − K S ν τ and finally, our conclusions are presented in section 5.

The pion vector form factor
The pion vector form factor has been measured in e + e − → π + π − [23,24,25,26,27,28,29] and in τ → π − π 0 ν τ [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 O(p 4 ) to follow with the explicit inclusion of resonances as degrees of freedom. The ππ 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.

Exponential Omnès representation
The pion vector form factor F π V (s) is defined through the hadronic matrix element of the vector current between the vacuum and the pion pair At very low-energies, the pion vector form factor is well described by ChPT. The resulting calculation at O(p 4 ) reads [3] (in the limit of exact isospin symmetry) where L r 9 (µ) is one of the renormalized low-energy couplings constants in the O(p 4 ) chiral Lagrangian, and A P (s, µ 2 ) are two-pseudoscalars loop-functions (whose renormalizationscale dependence cancels out with the one in L r 9 ) accounting for the unitary corrections given by 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 ρ(770) resonance dominates the form factor. At leading order in powers of 1/N C , which it is O(p 4 ) in the chiral expansion, the result is given by [4] F π V (s) = 1 + where F V and G V measure, respectively, the strength of the ρV µ and ρππ couplings, with V µ being the quark vector current. Assuming that the form factor vanishes when s → ∞ one gets the condition that yields the usual Vector Meson Dominance (VMD) in the zero-width approximation 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 O(p 4 ) chiral coupling L r 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 ππ 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 However, the previous expression still has an obvious defect, it cannot describe the energy region of the peak of the ρ-meson. For that, it is necessary to incorporate its width. The energy dependent width of the ρ is related to the imaginary part of the loop function and is given by [40] Γ In order to account for the ρ-resonance width, we insert Γ ρ (s) in the propagator of the ρ in Eq. (2.10) arriving at: 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 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 O(p 6 ). This effect, however, is seen to be numerically negligible in the ππ and Kπ systems. In all, equation (2.12) provides a suitable description of the ρ-dominance region. However, the precise measurement of the τ − → π − π 0 ν τ decay spectrum by the Belle collaboration [31] indicates that heavier resonance excitations i.e. ρ ≡ ρ(1450) , ρ ≡ ρ(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 ρ. The resulting expression takes the form: where the coefficients γ and δ measure the relative weight between the contributions of the different resonances while the phases φ 1 and φ 2 stand for the corresponding interference, and with (2.14) Notice that, for the energy-dependent width of the ρ and ρ given in the previous equation 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 τ , and matched smoothly to the ππ 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 ρ , γ, φ 1 , M ρ , Γ ρ , δ, φ 2 , M ρ , Γ ρ }, that can be inferred from fits to the measured modulus squared of the pion vector form factor extracted from the Belle τ − → π − π 0 ν τ measurement [31]. 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 χ 2 /d.o.f, the agreement with data is very satisfactory. In the figure, the ChPT calculation at O(p 4 ) (cf. Eq. (2.2)) is also shown (dashed gray line) for illustrative purposes 3 .
2 In all our fits throughout the paper, whenever the pion form factor is involved, we employ F π = 92.316 MeV, which is the central value using the restriction √ 2|V ud |F π = (127.13 ± 0.02 ± 0.13) MeV from the 2018 PDG edition [42] with the 2006 PDG reported value |V ud | = 0.97377 ± 0.00027 used by Belle. 3 For simplicity, in the figure we do not represent the ChPT calculation at O(p 6 ) [32,33] since including higher-order chiral corrections increases very little the energy region where the data is described well.

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 th = 4m 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 ππ 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 ππ scattering phase shift δ 1 1 (s) concerning us, reads The phase δ 1 1 (s) entering the dispersive integral encodes the physics of the ρ-meson and it is known with an excellent precision in the elastic region s ≤ 1 GeV 2 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 δ 1 1 (s) can be estimated between the region of 1.3 GeV and the mass of the τ lepton ∼ 1.8 GeV. The precise Belle measurement of the τ → π − π 0 ν τ decay distribution [31] indicates relevant contributions from the excited resonances ρ and ρ 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 that is only valid in the τ decay region (s < m 2 τ ). For the high-energy region (s > m 2 τ ) we guide smoothly the phase to π [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] 4 . 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 , (2.19) where α 1 and α 2 are two subtraction constants that can be related to chiral low-energy observables, namely the squared charged pion radius r 2 π V and the coefficients of O(s 2 ) and O(s 3 ) terms in the chiral expansion, c π V and d π V , respectively, appearing in the low-energy expansion of the form factor: Explicitly, the relations for the linear and quadratic slope parameters r 2 π V and c π V read These subtraction constants can be calculated theoretically through the sum rule 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 can be obtained through 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 parameters 5 entering F π V (s) to be determined by a fit to the Belle data. Regarding the integral cutoff s 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 cut , namely m 2 τ (third column), 4 GeV 2 (fourth column), 10 GeV 2 (fifth column) and finally the s cut → ∞ limit (last column). The choice of s cut = m 2 τ is motivated by the fact that the model used to get the phase, Eq. (2.18), is only valid within the τ decay region s ≤ s cut ∼ m 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 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 cut = m 2 τ only for illustrative purposes throughout the paper. Our reference fit corresponds to s cut = 4 GeV 2 (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 cut . The results obtained by varying s cut in Table 1 will be used to assess the systematic uncertainties of our fit results obtained with s cut = 4 GeV 2 6 . We shall return to a discussion on the integral cutoff below.
In Fig. 2, we show the resulting phase shifts for the chosen s cut . The phase shift solution of the Roy equations is given by the solid black curve while the variations due to s cut are given by the dot-dashed blue (s cut = m 2 τ ), solid red (s cut = 4 GeV 2 ), dashed green (s cut = 10 GeV 2 ) and dotted black (s cut → ∞) curves, respectively. In the figure, the statistical uncertainty associated to our reference fit (s cut = 4 GeV 2 ) is also shown by 5 The parametrization for the phase shift δ 1 1 (s) of Ref. [41] contains a parameter for the ρ-meson mass, that we name m ρ , that denotes the energy at which the phase shift passes through π/2 (and therefore it shall not be confused with real part of the pole of the ρ) and its quoted value is m ρ = 773.6(9) MeV. For our study, in a first approximation we fix the model input parameter for the ρ-meson mass, M ρ in Eq. (2.13), to this value. However, will also test the sensitivity of our fits to this parameter by allowing it to float. the light red error band. Also, there are two brown dashed vertical lines shown in the figure. They are placed at 1 GeV and √ s = m τ and denote, respectively, the range where the phase shift from ππ 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 cut = 4 GeV 2 (solid red curve and light red error band) and s cut → ∞ (dotted black curve) compared to Belle data [31] . As can be seen from the figure and the corresponding χ 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 τ − → π − π 0 ν τ data, we have also run fits allowing to float the parameter for the ρ-meson mass entering δ 1 1 (s). The corresponding fit results are contained in Table 1 as Fit 1-ρ 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 ρ 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 ρ 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 ππ into the ρ and ρ decay widths. We next allow such resonances to decay not only in ππ but also into KK. In this case, the energy-dependent widths read (2.24) 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 χ 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 ρ width decreases by ∼ 140 MeV while the ρ mass increases up to 1775 MeV. The ρ width is slightly shifted downwards but the associated error is enlarged.
Finally, we also come back to the discussion on the integral cutoff s cut and the corresponding generated singularities. We found that s cut = 4 GeV 2 is a suitable value for the integral cutoff in the dispersive integral. However, we would also like to consider a This work s cut Figure 3: Belle measurement of the absolute value squared of the pion vector form factor |F π V | 2 (black filled circles) [31] as compared to our fit results as presented in Table 1 for s cut = 4 GeV 2 (red solid line) and for s cut → ∞ (black dotted curve).
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 δ 1 1 (s) which is valid up to s = m 2 τ and the idea is to extend it to the full region in an appropriate way. We follow Refs. [49,50] and write where the function f (s) is given by a once subtracted dispersion relation, that ensures . (2.28) The value of the phaseδ 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 → ∞. In order to fulfill these properties, we choose a smooth interpolation in s for δ 1 1 (s) above s cut as simply asδ = 773.6(9) = 773.6(9) = 773.6(9) = 773.6(9) (6) 1374 (8)   = 773.6(9) = 773.6(9) = 773.6(9)  ) with s cut = 4 GeV 2 in the dispersive integral with (second column) and without (last column) the KK channel in the ρ and ρ energydependent resonance widths, and with the additional inclusion of the ρ → ωπ contribution into the ρ width (third column). See main text for details. so thatδ 1 1 (s cut ) = δ 1 1 (s cut ) and δ 1 1 (s) → π for large s recovering the 1/s fall-off of F π V (s). In this case, the integral going from s cut to ∞ in Eq. (2.28) can be calculated explicitly and we arrive at [49,50] Regarding the function Σ(s), it contains the (inelastic) contributions beyond s cut and may be understood as giving the correction to the linear continuation of the phase δ 1 1 (s) above s cut as we did in Eq. (2.29). It is described by an analytical function on the s-plane with a cut from s cut to ∞ 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 with the variable ω(s) given by For our analysis, we take the condition a 0 = 1 to ensure F π V (0) = 1. We next probe the application of Eq. 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-ρ with s cut = 4 GeV 2 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 ∼ 2.5 GeV 2 caused by the destructive interference between the ρ and ρ resonances. To present our central results, we quote the values of our reference fit in Table 1 with s cut = 4 GeV 2 (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 cut and the matching point with the Roy equations, respectively, in Table 3 due to the inclusion of the KK decay channel into the ρ and ρ widths and in Eq. (2.34) due to the parametrization presented in Eq. (2.30) that avoids the aforementioned singularities. We then obtain    Table 1, Fit 1 (solid red curve), and from Fit 1-ρ at the matching point of 0.85 GeV (dotdashed blue curve), Fit I (short dashed purple curve), Fit A (long dashed green curve) and with the fit that avoids singularities (dotted black curve). The vertical dashed brown line is placed at m τ and denotes the validity of the parametrization of the form factor phase shift. All phases are smoothly guided to π for s > m τ . See main text for details.
In Table 4, we show our results for the ρ(1450) and ρ(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 ρ(1450) and ρ(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   Table 1, Fit 1 (solid red curve), and to our Fit 1-ρ at the matching point of 0.85 GeV (dot-dashed blue curve), Fit I (short dashed purple curve), Fit A (long dashed green curve) and with the fit that avoids singularities (dotted black curve). See main text for details. data. For the ρ(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 PDG ρ = 1450 ± 25 MeV and Γ PDG ρ = 400 ± 60 MeV [42]. For the ρ(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 PDG ρ = 1720 ± 20 MeV, and ∼ 300 MeV higher than [8].
On the other hand, as seen from the results quoted in the table, the values for the ρ(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], ∼200 MeV higher than the PDG value Γ PDG ρ = 250 ± 100 MeV, in line with [51] and ∼150 MeV lower than [8]. Due to the large uncertainties associated to the ρ(1450) and ρ(1700) resonances parameters, however, we cannot draw any more definite conclusion. For that, more precise data in the ρ(1450) and ρ(1700) region would be highly desirable.
Regarding the chiral observables associated to the low-energy expansion of the form   ) resonances. The first and second errors refer, respectively, to the statistical and systematic uncertainties.
It is opportune to mention again that we have treated the subtraction constants α 1 and α 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.  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 , is much less known. To the best of our knowledge, there are neither theoretical results from ChPT nor calculations

Predictions and fits to τ − → K − K S ν τ BaBar data
The theoretical expression for the differential decay distribution for the transition τ − → K − K 0 ν τ in terms of the K − K 0 invariant mass can be written as [35,64,65] and it is related to the normalized invariant mass spectrum through 2) where N events is the total number of measured events, Γ τ is the inverse τ lifetime and ∆ bin ≡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 fractionB = 0.739(11) stat (20) syst × 10 −3 [16] 7 . 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 ± 3461 [16] and the bin width is 0.04 GeV. F K V (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 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 K V (s) at O(p 4 ) in ChPT we need the expressions of the K + K − and K 0K 0 form factors. These can be found in the literature and read [3] (3.5) The I = 1 component corresponding to the K − K 0 state can be extracted from Eqs. (3.4) and (3.5) and yields We would like to note that the ChPT calculation at 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 τ − → K − K S ν τ 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 τ − → K − K S ν τ 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, BR(τ − → K − K S ν τ ) = 0.545(32) × 10 −3 , which is seen ∼ 5σ away the BaBar measurement BR(τ − → K − K S ν τ ) = 0.739(11) stat (20) syst ×10 −3 . From these results, we conclude that using the pion vector form factor to describe the τ − → K − K S ν τ 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 ρ and ρ vector resonances to play a significantly different role than in τ − → π − π 0 ν τ as also noted in [66,67]. Thus, the weight of each resonance contribution, represented by the coefficients γ and δ and the phases φ 1 and φ 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 τ − → K − K S ν τ 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: where the coefficientsγ andδ and the phasesφ 1 andφ 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 , (3.8) and this is inserted 8 into a three-times-subtracted dispersive representation of the form factor , (3.9) whereα 1 andα 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 distribution 9 and found on the one hand that, in full generality, the data is not sensitive either to the low-energy region or to the ρ(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 ρ + Γ ρ , the energy region more influenced by the presence of the ρ(770). This implies first that the slopes of the kaon vector form factor, which encode the physics immediately above threshold, cannot be fitted with τ − → K − K S ν τ data and second that the fits lead unrealistic parameters under floating the ρ-meson parameters. On the other hand, the data is scarce in the ρ(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 ρ-meson mass (773.6(9) MeV), and considered fits varying only the ρ(1450)-resonance mass and width, andγ andφ 1 , while neglecting the contribution of the ρ(1700) resonance to the decay i.e.δ 1 = 0. In Table 7, we show the results of our fits using different settings. Fit i) corresponds to fixing the slopesα 1 andα 2 , to the slopes α 1 and α 2 given in Eq. (2.35) obtained from the analysis of the pion vector 8 The δ KK 1 phase as extracted from Eq. (3.8) is also matched to the ππ scattering at 1 GeV as explained along the lines of section 2. 9 We 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.   [16] for the decay τ − → K − K S ν τ (black solid circles) confronted to our prediction (dotted red curve), as for our central results of the pion vector form factor analysis presented in Eq. (2.35), and fits with the exponential representation (dashed blue curve) and the dispersive approach with s cut = 4 GeV 2 (solid green curve).
form factor, while Fits ii) and iii) are variants of it. In particular, Fit ii) is the result of fixing the slope toα 1 = 1.84 (30) obtained from the kaon radius of Ref. [33] while Fit iii) includes the intermediate KK state into the ρ 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 χ 2 /d.o.f is seen for Fit ii) but the resulting fit parameters carry the larger error, and the (large) uncertainty associated toα 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    approaches, and M pole ρ = 1370 ± 15 MeV and Γ pole ρ = 373 ± 30 MeV for the exponential representation (Fit iv)). From these results we conclude that while the pole mass of the ρ(1450) resonance as extracted from the τ − → K − K S ν τ 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γ 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 τ − → K − K S ν τ 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(τ − → K − K S ν τ ) = 0.749(93)×10 −3 (Fit i)) and BR(τ − → K − K S ν τ ) = 0.744(89) × 10 −3 (Fit iv)), the quality of the fit as indicated by the χ 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 | 2 and the BaBar data of the decay τ − → K − K S ν τ . Such analysis shall allow us to determine the ρ(1450) and ρ(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 | 2 Belle and τ − → K − K S ν τ BaBar data The χ 2 minimised in our simultaneous fit is where the first and second terms correspond, respectively, to the Belle pion vector form factor data [31] and to the BaBar τ − → K − K S ν τ measurement [16]. For the later, N exp j and σ N exp j are the experimental normalized number of events for τ − → K − K S ν τ 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, α 1,2 andα 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 ρ and ρ resonances, M ρ ,ρ and Γ ρ ,ρ , used to model the phase entering the dispersive integral. The parameter for the ρ-meson mass, M ρ , is taken equal to that entering the phase shift, m ρ , and its quoted value is fixed to 773.6(9) MeV as discussed in section 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 cut = 4 GeV 2 in the dispersive integral is employed. Fit a (second column) corresponds to fixingα 1 = α 1 andα 2 = α 2 and takingδ = 0 because, as we have discussed in section 3, the BaBar measurement of the decay τ − → K − K S ν τ is still not sensitive to ρ(1700) resonance properties. The corresponding fit results supports the relative weights γ andγ to be the same for the π − π 0 and K − K S channels. This feature is proven in Fit b (third column) by enforcing γ =γ. By doing this, the χ 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 ρ -width, whose central value is shifted downwards by ∼ 100 MeV, and to less extent for the ρ -mass, which suffers a variation of ∼ 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α 1 = α 1 and α 2 = α 2 made in Fits a and b. As a side result, we extract the charge kaon radius through r 2 K V = 6α 1 (cf. Eq. (2.21)). Our value, r 2 K V = (11.28 ± 1.44) GeV −2 , lies in the ballpark of results for this quantity r 2 K V = (9.09 ± 1.82) GeV −2 and r 2 K V = (11.07 ± 1.82) GeV = 773.6(9) = 773.6(9) = 773.6(9) (19) 1453 (  The results of our fits are confronted to the Belle |F π V | 2 form factor measurement and to the BaBar τ − → K − K S ν τ distribution in Figs. 7 and 8, respectively. Satisfactory agreement with the experimental data is seen in accord with the observed χ 2 /d.o.f. However, the potential of a combined analysis of the |F π V | 2 and τ − → K − K S ν τ data is presently limited by the fact that the errors associated to the latter are still relatively large and the BaBar measurement of the τ − → K − K S ν τ spectrum is yet not sensitive to the ρ(1700) resonance properties. This presents a limitation in determining the ρ(1450) and ρ(1700) resonance parameters with improved precision with respect to the individual analysis of |F π V | 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.  Figure 7: Belle measurement of |F π V | 2 (black filled circles) [31] as compared to our fits obtained from a simultaneous analysis of |F π V | 2 and τ − → K − K S ν τ (see Table 8).  [16] for τ − → K − K S ν τ (black solid circles) as compared to our fits obtained from a simultaneous analysis of |F π V | 2 and τ − → K − K S ν τ (see Table 8).

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 τ → π − π 0 ν τ , and provided a parametrization for the kaon vector form factor that describes the decay τ − → K − K S ν τ , 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 ππ 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 ππ interactions capturing the physics of the ρ-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 ππ scattering phase shift existent in the literature. Above 1 GeV, the ρ(1450) and ρ(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 ππ 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 ρ(1450) and ρ(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 τ − → K − K S ν τ 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 ρ(770)-dominance region, the role of the ρ(1450)-resonance is different than in the π − π 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 ρ(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 τ − → K − K S ν τ will allow one to determine the ρ(1450) and ρ(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.

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 FormFactor-Phase.dat and FormFactorModulusSquared.dat, respectively. In the FormFactorPhase.dat file, the first column is √ 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 GeV 2 ) and the second column is the value for the form factor modulus squared for each fit, and is given until s = m 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.