A simple alternative to the Breit-Wigner distribution

First, we discuss the conditions under which the non-relativistic and relativistic types of the Breit-Wigner energy distributions are obtained. Then, upon insisting on the correct normalization of the energy distribution, we introduce a Flatt\'{e}-like relativistic distribution -- denominated as Sill distribution -- that (i) contains left-threshold effects, (ii) is properly normalized for any decay width, (iii) can be obtained as an appropriate limit in which the decay width is a constant, (iv) is easily generalized to the multi-channel case (v) as well as to a convoluted form in case of a decay chain and -- last but not least -- (vi) is simple to deal with. We compare the Sill distribution to spectral functions derived within specific QFT models and show that it fairs well in concrete examples that involve a fit to experimental data for the $\rho$, $a_1(1260)$, and $K^*(982)$ mesons as well as the $\Delta(1232)$ baryon. We also present a study of the $f_2(1270)$ which has more than one possible decay channels. Finally, we discuss the limitations of the Sill distribution using the $a_0(980)$-$a_0(1450)$ and the $K_0^\ast(700)$-$K_0^\ast(1430)$ resonances as examples.


Introduction
The parametrization of the energy distribution (or spectral functions) of resonances is an important topic of both Quantum Mechanics (QM) and Quantum Field Theory (QFT). Namely, new resonances are often obtained by fits involving the line-shape of unstable states, see e.g. the Particle Data Group (PDG) [1], and the activities in various ongoing and planned experiments [2,3,4,5,6,7,8,9,10].
The very well known Breit-Wigner (BW) distribution [11,12]: where M is the energy (or mass) of the unstable state and Γ the decay width (the lifetime being simply 1/Γ in natural units), has been widely used in countless applications. Quite remarkably, the BW distribution has no low-energy threshold and the corresponding decay law is exactly exponential. It is correctly normalized, +∞ −∞ dEd BW S (E) = 1, allowing for the standard interpretation of d S (E)dE as the probability that the unstable state possesses an energy (or mass) between E and E + dE. It may be regarded as a very good approximation in most cases, especially in the context of QM when Γ M. In general, Γ is itself a function of energy and is interpreted as the imaginary part of the self-energy of the one-particle unstable state; the BW function naturally arises as the limiting case in which the coupling of the unstable state to the final decay products is energy independent.
In the relativistic case, the Breit-Wigner formula (rBW) for an unstable state S reads [1,13,14] (see later on for details): where, for relativistic particles, the variable E represents the particle running mass and must be non-negative. The normalization of the spectral function, necessary for its probabilistic energy interpretation, is not fulfilled since ∞ 0 dEd rBW S (E) < 1. The loss of normalization, even if small for narrow states, is due to the fact the limiting process that works in the QM case does not apply in the relativistic case. Technically speaking, when the imaginary part of the self-energy is a constant for E ≥ 0, the real part of it is not a simple constant that can be subtracted.
In this work we aim to show that a rather simple spectral function can be obtained as a natural limit of the relativistic context. First, we introduce the threshold E th associated with a certain unstable state (for a two-body decay it is simply E th = m 1 + m 2 , where m 1 and m 2 are the rest masses of the produced particle). We then introduce a Flatté-like spectral function [15,16] that we call "Sill"(Sill stands for threshold and for the form of the function close to it): where the rescaled widthΓ is given byΓ = ΓM/ M 2 − E 2 th (for a nonrelativistic Flatté-like parameterization, see Refs. [17,18,19]).
The function d Sill S (E) is always normalized to 1 ( ∞ E th dEd Sill S (E) = 1) for any M and Γ. The form d Sill S (E) is very similar to the relativistic BW distribution(s) when Γ M, but is utterly different when Γ/M is sizable. Moreover, the constraints of complex analysis are correctly taken into account (details of its derivation and properties are in the paper): the function can be seen as the proper limit of a distribution in which the decay width is a constant (the real part of the self-energy vanishes). Quite importantly, d Sill S (E) is a continuous function for any value of E.
We therefore argue that a fit to data using the function in Eq. (3) may be advantageous. If BW or rBW work well, then the Sill also does: it simply delivers similar results. Yet, whenever a threshold is close enough to the mass of the unstable state, d Sill S (E) fairs better, as we shall show with some concrete examples.
Indeed, Eq. (3) is closely reminiscent of the Flatté distribution [15] employed for the resonance a 0 (980) meson decaying into two channels,KK and πη, where however typically only the KK channel is treated as energy dependent, while the πη decay channel is treated as a constant, see e.g. Ref. [16]. In the context of the resonance a 0 (980), theKK channel is very close to the nominal mass, hence the peculiar shape of the spectral function can be described.
Our claim is that Eq. (3) -as well as multi-channel and /or convoluted extensions of it-can be also used independently on the location of threshold(s). In addition, the Sill distribution can be applied as it stands to the case in which the decay product consists of particles of different masses without any loss of normalization of the spectral function 1 . Summarizing, it may be regarded as a useful phenomenological approximation able to capture in first approximation some essential (but of course not all) features needed to properly describe various physical resonances. In this respect, it must be seen as a valid alternative to the rBW function because it contains threshold effect(s) and is mathematically consistent. Surely, it is not an alternative to a detailed study of a given resonance, in which the proper coupling to the decay products is taken into account. Yet, it can be used for both standard and non-standard mesons (such as the a 0 (980) mentioned above), since it is not dependent on the internal structure of hadrons In Fig. 1 we compare the plots of the BW, rBW, and Sill distributions for the ρ and a 1 (1260) mesons using the parameters of the PDG. This picture serves a first illustration of the differences of these functions, which are evident for broad resonances (right plot). In the following, for these two resonances as well as for the narrow K * (892) we shall perform a fit to data. We shall see that the Sill fairs better than BW and rBW for the ρ and especially for the a 1 (1260), for which the threshold effects are important, but there is basically no difference for the narrow K * (892), since then Sill reduces to rBW in the small-width limit.
The paper is organized as follows: in Sec. 2 we discuss the non-relativistic case by using the propagator formalism and derive the standard BW spectral function; in Sec. 3 we introduce the relativistic propagator, the standard relativistic approximations as well as the decay of a scalar QFT (note, Secs. 2 and 3 contain known results which are useful for the comparison with the Sill distribution); in Sec. 4 we introduce the Sill distribution, compare it to the results of Sec. 3 and extend it to the multi-channel and convoluted cases; in Sec. 5 we discuss various numerical examples, that include fit of the BW, rBW, and the Sill to experimental data for the resonances ρ, a 1 (1260), and K * (982) mesons and the ∆(1232) baryon. Moreover, we shall study the f 2 (1270) as an example of multi-channel decays. In Sec. 6, we present the a 0 (980) and k as examples of non-conventional states, and present the limitations of our approach. In Sec. 7, we discuss the domain fo applicability of the Sill and give a empirical condition when the use of Sill will be favorable. Finally, in Sec. 8 we present our conclusions.
2 Brief review of the non-relativistic case In this section we recall the derivation of the non-relativistic BW function. To this end, we start from the propagator of an unstable quantum state |S as function of the energy E: where the self-energy function Π(E) for complex E can be expressed via the dispersion relation with E th being the lowest energy threshold for the decay of the unstable state. (For a well defined theoretical setup in which the above equations can be defined, we refer to the Lee (or Friedrichs) model, see Ref. [20,21,22,23,24,25,26] and refs. therein). The function is the so-called decay width function, whose specific form depends on the details of the interaction. A possible definition of the mass of the unstable state S is obtained via which shows that quantum fluctuations may shift the bare mass M . Of course, with a suitable subtraction, one can impose that Re Π(M ) = 0, thus M r = M. We shall adopt this choice, thus M can be regarded as the nominal energy/mass of the state. The on-shell decay width is obtained by setting Γ = Γ(M ). Alternatively, a commonly used definition is the pole energy/mass z pole obtained by the equation where II refers to the second Riemann sheet. Then: If more poles are present, typically the one(s) closest to the real axis is (are) the most important, but exceptions are possible. The position of the pole, being independent of the particular reaction needed to form the unstable state, is regarded as a stable feature of a given resonance. The spectral function is obtained as the imaginary part of the propagator or conversely (11) in agreement with the interpretation of dE d S (E) as a mass/energy probability between E and E + dE [14,27]. Namely, the full propagator has been expressed as the 'sum' of free propagtors with mass/energy E , each one weighted by d S (E ). The function d S (E) is correctly normalized to unity as long as lim E→∞ |Π(E)| /E = 0 (see Refs. [24,28]): In fact, one may consider the state |S as the eigenstate of the unperturbed Hamiltonian H 0 which fulfills the normalization condition S|S = 1. The full set of eigenstates of the Hamiltonian H reads {|E } with H |E = E |E and E ≥ E th . Expressing |S in terms of |E implies The quantity d S (E) = |a S (E)| 2 is the 'spectral function' which is evaluated in this work as the imaginary part of the propagator (see also Ref. [25] for a more detailed discussion of these relations). It naturally follows that The survival probability amplitude as function of time reads: In general, the decay is not exponential, see Refs. [22,29,30,31] for a description of this issue. Quite interestingly, also in presence of deviations from the exponential decay, the pole contribution implies that where the exponential contribution turns out to be dominant for intermediate times.
The Breit-Wigner limit is obtained in two steps. First, let us consider Γ(E) = Γ = const and introduce a maximal energy Λ: It follows that Next, one considers the limit E th → −∞ and Λ → ∞ together with the constraint E th /Λ = −1.
Clearly, Im Π(E) = Γ/2 for each E. On the other hand: (another choice for the E th /Λ ratio could be re-absorbed in the definition of M.) Hence, the Breit-Wigner propagator reads whose pole is simply given by The spectral function in this case takes the very well known Breit-Wigner form The survival probability amplitude is in this case exactly an exponential: The normalization is clearly fulfilled. Yet, the Breit-Wigner distribution can be only considered as a useful but unphysical limit, the most clear drawback being the lack of the left energy threshold (the energy is unbounded from below). Yet, it is very good close to the peak for sufficiently narrow (alias long-lived) quantum states.

Brief review of the relativistic case
In this section we present how the previous formalism can be generalized to the relativistic case. We shall describe the series of approximations leading to the relativistic BW distribution. Moreover, we shall also discuss a scalar QFT as a notable example.

Propagator and spectral function
The relativistic case can be easily obtained by replacing the variable E with the variable s = E 2 into Eq. (4), hence the propagator (as function of s) takes the form (see e.g. [32,33]) It is clear from the very beginning that s ≥ 0, hence in a relativistic framework a left threshold is built in. The self-energy reads [34] with s th = E 2 th ≥ 0 is the low energy threshold. In this case, the imaginary part takes the form The mass is obtained by Also in this case, a suitable subtraction can be made to assure that Re Π(M 2 ) = 0, and hence, M r = M. The corresponding on-shell decay width is then Γ = Γ(M ), see below.
Alternatively, the pole position s pole is such that and leads to the pole mass and decay width defined as: The spectral function as function of s is obtained as the imaginary part of the propagator then is the Källen-Lehmann representation in QFT, e.g. Ref. [32]. As function of the energy E one obtains 2 Just as in the non-relativstic QM case as long as lim s→∞ |Π(s)| /s = 0. In general, the function Π(s = E 2 ) depends on the employed model. The imaginary part Im Π(s = E 2 ) can be calculated in the framework of a given QFT. In the easiest case, it is evaluated at tree-level, see e.g. Ref. [27] and Section 3.3 for a concrete example. Higher orders can also be evaluated but are technically more involved, even though in general the correction with respect to the tree-level result is small [35].
The time evolution of the state in QFT reads [25,29] Note, the first integral extends from −∞ to +∞, i.e. it includes also unphysical values of the variable s.

Standard relativistic approximation
In order to obtain possible expressions for the relativistic BW distribution, let us consider which closely resembles the BW case of Sec. 2 by including a lowest threshold at zero. Of course, Im Π(s) = M Γθ(s) is itself not physical, but -as we shall see -additional approximations are needed to obtain the standard relativistic BW distribution.
The loop function takes the form (detailed studies of the structure of the loops have been carried out before on a case by case basis, see, e.g, [36]): By including a subtraction constant C = − M Γ π ln M 2 Λ 2 and by taking Λ → ∞ we get: which assures that Re Π(s = M 2 ) = 0, thus the mass is still M even after including the oneparticle self-energy. The important point is that the real part of the loop Re Π(s) = M Γ π ln M 2 s is not a simple constant. The propagator reads Note, the pole has no simple algebraic solution. The corresponding spectral function takes the form which is correctly normalized to 1, as it should. The propagator for the relativistic BW approximation is obtained by approximating the previous propagator upon setting the real part of the loop artificially to zero, thus finding: whose corresponding pole is Indeed, the form of Eq. (42) is directly used in many applications, e.g. Refs. [37,38,39,40] and refs. therein. The rBW spectral function as function of s reads to be understood for s ≥ 0, hence the step function has been added. The related spectral function, as a function of the energy E takes the form that contains an unphysical jump at E = 0. We refer e.g. Refs. [41,42,43,44] for a practical use of this function 3 This distribution function is in fact an approximation where, the factor of EΓ in the numerator has been replaced 4 by M Γ.
The normalization to 1 is lost due to the employed approximation: the real part is not a constant and therefore cannot be completely subtracted. Of course, one may normalize the distribution of Eq. (46) by "brute force "upon considering but this step is not mathematically consistent. wit Curiously, as we shall show in Sec. 5, the rBW fits to the data in a similar way as the BW. As shown in the numerical examples, the two plots are not discernible to the naked eye. However, the pole positions and the decay widths are marginally different.

A scalar QFT
In this context it is worth reminding the results for the QFT case in which a scalar particle decays into two (pseudo)scalars. The corresponding interacting part of the Lagrangian takes the form [27] out of which (at one-loop): The on-shell decay width is given by The loop takes the form Let us then take the limit Λ → ∞ and introduce a subtraction constant in such a way that Re[Π(M 2 )] = 0. Upon employing arctan (z) = 1 2i log i−iz i+iz ) one has: Then, the propagator is perfectly defined and the spectral function is correctly normalized. On the other hand, the real part of the loop is a nontrivial function with a cusp at 4m 2 , which is nonzero above threshold [27].
The present model has a great advantage with respect to the BW and rBW cases, since it correctly implements the lowest energy threshold and it arises from a -albeit simple -QFT approach. Yet, it would be awkward to use such a function in fitting experimental data for various reasons: 1) the QFT model of Eq. (48) can be only considered as a first approximation for physical resonances.
2) Even in the simplest decays of a scalar particles into two (pseudo)scalar ones there can be additional terms, such as S∂ µ ϕ∂ µ ϕ [34,45,46,47,48]. For instance, this is the case for the spectral function of the resonances a 0 (980), a 0 (1450), K * 0 (700), and K * 0 (700). 3) Decaying particles (especially hadrons) involve also particle with spins. For the use of spectral functions for various vector particles, see for instance Ref. [49] for the vector field ψ(3770), Ref. [50,51,52,53,54,55] for the X(3872), Ref. [56] for the state ψ(4040), Ref. [57] for Y (4260), and Ref. [58] for the ψ(4660). 4) Typically, when derivative and/or higher spin are considered, the divergence of the loop integral is not just logarithmic, hence additional subtractions or cutoff functions are needed (note, one may also use explicit form factors, which however corresponds to nonlocal Lagrangians, see e.g. Ref. [59,60] and refs. therein). 5) Quite interestingly, the form of Eq. (49) has been used in a variety of applications, see for instance the recent paper by the ALICE collaboration [61] (but the coprresponding real part has been neglected).

A useful relativistic Flatté-like spectral function: the Sill
In this section we present a relativistic spectral function that properly includes the threshold (with no unphysical jump) and is always correctly normalized for both equal or unequal masses in the final state. Moreover, it is easily extendable to the multichannel case and to decay chains.

Definition and properties
Let us consider a certain decay whose threshold is s th = E 2 th . For a two-body decay, one has E th = m 1 + m 2 = √ s th . We then assume the following simple imaginary part of the loop whereΓ is a rescaled width linked to the standard decay width Γ via the relation ΓM = Γ M 2 − E 2 th , see below. Notice that in this context the function Γ(s) is actually not a constant. In fact, from Im Π(s) = √ sΓ(s) it follows that (for s ≥ s t h) which on-shell reduces to Yet, even if not a constant, the width Γ(s) saturates for large s to a constant equal toΓ. For M sufficiently larger than E th , it follows that Γ Γ . The loop takes the form Distribution  Table 1: Parameters used to calculate the spectral function as given in eq.(94) of the text.
By taking the limit Λ → ∞ together with one subtraction, the loop function reduces to which is then only imaginary for s ≥ s th (conversely, it is only real below threshold). Note, the function Π(s) for complex values of the variable s is a multivalued function with two Riemann sheets (RS); the branch cut is taken for s ≥ s th ≥ 0. The expression above is therefore on the I RS, while the expression on the II RS is simply given by Π II (s) = −Π(s). For real s < s th the function Π(s) is real, thus the Schwarz reflection principle is fulfilled. The spectral representation of Π(s) (obtained with one subtraction) takes the form: valid for any s besides the cut (s th , ∞).
The propagator takes the form whose relevant pole (on the II RS)is Note, forΓ 2 sufficiently smaller than M 2 − s th , the pole of s can be approximated as just as in the rBW case. Next, the spectral function reads (65) or, as function of the energy E, as: The normalization  770). Experimental data corresponds to the vector channel of the τ → ππν τ decay reported in [67].  Table 2: Mass and widths of ρ(770) fitted using the three distributions discussed in the text, their error estimates, and the poles.
for any E th , M , andΓ is a consequence of the proper treatment of the real part of the loop according to the general proof reported in Refs. [24,28]. Of course, it can be also directly proven for the specific spectral function d Sill S (E). Some comments are in order: 1) The Sill distribution correctly reduces to the rBW, d sill both in the numerator and denominator. This is allowed when the width is small. It is also interesting to notice that d Sill S (E) scales for large E as E −2 , just as in the BW case (and not as E −3 as in the rBW of Eq. (2)).
2) A comparison with the QFT case of Sec. 3.2 is useful. In that case, one had hence the difference is clear: the additional energy dependence in the denominator does not allow for a constant real part. Yet, the numerator is the same of the Sill distribution in the case of equal masses. The QFT case of Section 3.3 translates into the Sill form upon the substitution a step which is necessary to get a vanishing real part of the self-energy. Here, it must be stressed that the modification of the high-energy tail does not constitute a drawback for the Sill 5 , since the high-energy region has a small influence on the behavior of the distribution close to the peak.
On the other hand, the fact that for the Sill Γ(s → ∞) =Γ (i.e. it is a constant for high energies) implies that in this respect the Sill closely resembles (r)BW. Also, for the Sill the real part of the self-energy vanishes, in agreement with BW (but not with rBW). Colloquially speaking, we might say that the Sill is "closer" to the BW than rBW itself, and it achieves that by properly taking into account the threshold.
3) Next, let us consider the case of the QFT decays in two distinct particles with different mass. The corresponding scalar QFT Lagrangianes is L int = gSϕ 1 ϕ 2 , then the imaginary part reads with Thus, the Sill distribution proposed in Eq. (3) amounts to: Indeed, the behaviors of the functions k(E) and E 2 − E 2 th are quite similar for E > E th , but the latter is easier to implement. Moreover, in the case of unequal masses, the choice of Im Π(s) as proportional to k(E) does not lead to a properly normalized spectral function (even when substituting the factor E in the denominator with a constant M ). Since this requirement is central to our approach, the choice implemented in the Sill is both simpler and theoretically favoured.
4) Finally, it should be noticed that in some cases also the generalized form involving the angular momentum l is introduced as yet for l > 0 the divergence of the loop function is too severe and the normalization of the spectral function is lost. This aspect shows that some sort of cutoff function is needed for a proper treatment of the problem and for a meaningful probabilistic interpretation of the mass/energy distribution. This is why we do not attempt to include l into the proposed phenomenological description introduced in Eq. (54) and to the spectral function d Sill S (E) of Eq. (3). As we shall see later on in concrete examples, the Sill fairs well also when higher values of the angular momentum l are actually involved.

Extension to the multi-channel case
In many cases unstable resonances decay in more than a single decay channel [1].
We consider as an illustrative example the case of an unstable state S with two open decay channels with rescaled partial decay widthsΓ 1 andΓ 2 (linked to the partial decay width via k,th /M ) as well as the two thresholds s 1,th and s 2,th ≥ s 1,th . The propagator is: Care is needed when deriving the spectral function, since different regions emerge: Clearly, the normalization still holds from the lowest threshold: This is a consequence of unitarity, which is a necessary feature of the unstable state S. Thus, the present generalization of the Sill is a simple, straightforward extension to the case of multiple decay channels that is both consistent and sufficient for our purposes. However, we stress that the present multichannel extension refers only to the shape of the spectral function of the unstable state in the framework of the used approximation(s), but not to the whole multichannel scattering problem that involve scattering amplitudes, see for instance [55].
The extension to the N channels is straightforward: The spectral function is which takes the explicit form for e.g. s Q,th ≤ s ≤ s Q+1,th : with Q = 1, ..., N . The normalization is fulfilled also in this general case. Note, in other intervals of s one should carefully check which pieces appear in the numerator and which ones affect the mass. Within this context it is also useful to introduce the partial i-th spectral function as The latter reads explicitly (for s Q,th ≤ s ≤ s Q+1,th ): where only the i-th channel is present in the numerator (the denominator being unchanged). The quantity d (i) S (s)ds may be interpreted as the probability that the unstable state has a squared energy between s and s + ds and decays subsequently in the i-th channel. As a consequence the integral ∞ s i,th d (i) S (s)ds (84) can be seen as the branching ratio for the decay in the i-th channel. Of course, the sum of all branching ratios is one. In Sec. 5 we shall apply the multichannel Sill to the resonances f 2 (1270) and a 0 (980).

Convolution of the Sill
Quite often, decay chains take place. As an example, let us consider the unstable state S decaying into where ϕ 1 is stable but the field ϕ 2 is itself unstable and decays into The spectral function of the state ϕ 2 reads (using the Sill for consistency, for other approaches see Refs. [62,63,64]) then the convoluted spectral function of the state S takes this aspect into account. Next, we considerd in which the mass of the field ϕ 2 is taken as the running quantity √ s . Then the convoluted Sill reads d conv-sill Intuitively, we are taking the average of all possible masses of the unstable field ϕ 2 weighted by d ϕ 2 (s ). Note, forΓ 2 → 0 one has d ϕ 2 (s) = δ(s − m 2 2 ), thus one reobtains the expected spectral function. The normalization of the convoluted Sill is also straightforward: where we used that Of course, one can also generalize the study of decay chains to the case of multiple decays. For specific examples, see the discussion about the resonances a 1 (1260) and f 2 (1270) in Sec. 5.
In general, the convoluted Sill distribution is more involved and requires more efforts to be used in a fit than the plain Sill and it can be usually omitted in first approximation. Nevertheless, in presence of precise data it can be an additional tool in some selected cases since it represents a straightforward, self-contained, and normalized extension of the Sill distribution.

Numerical examples
The spectral functions studied in the previous sections (BW, rBW, and Sill) enter the relation where F th represents the theoretical quantity that can be compared to data, g is a scaling factor (typically a coupling constant, but that depends on the considered experiment), and b is the background. Since g is often an important physical quantity, it is evident that the correct normalization of the spectral function d S (E) is relevant. The values of these coefficients used in the present work are listed in Table. 1. Indeed, for the ρ and the a 1 cases, the coupling constants can be seen as physical couplings, as explained in Ref. [67]. On the other hand, for the K * , the experimental value refer to row counts. The parameters listed in the subsequent tables were estimated by minimizing the χ 2 defined as for the respective distributions. The errors in the parameters were calculated using the Hesse matrix method as explained in Appendix B.
We stress that the examples that we shall present are illustrative and aim to show the utility of the Sill distribution as a useful function to fit data in first approximation. Yet, for each of the presented resonances, more detailed and advanced studies have been performed (e.g, see [65,66] for a review of various methods as well as the PDG review on Resonances [1]). The Sill is not a substitute of such advanced studies, but a simple useful tool to study resonances by taking into account threshold effects.

The ρ meson
The three functions F for the ρ(770), fitted using the three distribution functions (BW, rBW, Sill), along with the experimental data [67] are given in Fig. 2. The mass and the width used to fit the data, their error estimates, and the pole position are listed in Table 2.
We see that the BW and rBW distributions fit the data poorly (as indicated by the higher value of the reduced χ 2 ), compared to the Sill. There is also a marked difference between the mass and the width of the ρ(770) estimated by the Sill distribution compared to the other distributions. The difference between the three distributions is not so large close to the peak because the 2π threshold is far away. Yet, as expected, the Sill fits the left tail clearly better. This behavior plays an important role in the case of resonances with the poles closer to the threshold (e.g. the a 1 (1260) discussed next).
We notice that the Sill can describe the ρ spectral function relatively well even though the correct angular momentum is not properly taken into account. Namely, for the ρ, one has actually the scaling k(E) 3 . As the detailed study of Ref. [68] shows, a proper implementation of the angular momentum is important for an accurate description of data. Yet, the effects of the threshold are -in the first approximation-more relevant than the angular momentum ones.
The study in Ref. [68] shows another important point. In that work, the real part of the loop was evaluated by including two subtractions to guarantee convergence of the one-loop function. The treatment is more involved and definitely preferable whenever accurate data are available. On the other hand, the Sill distribution discussed here is much simpler since the real part can be (consistently) set to zero.
In conclusions, the Sill does not represent a substitute to a detailed study of a given resonance, but represents a simple function to start with, especially if some thresholds are close. It seems therefore that the role of the left energy threshold is enough for a fair description of data. These are good news, since the neglect (in first approximation) of the angular momentum renders the treatment of line shapes easier and systematically possible within the Sill function presented here.  Table 3: Mass and width of a 1 (1260) fitted using the three distributions discussed in the text, their error estimates, and the poles (as described in the text).

The a 1 (1260) meson
As discussed before, the non-relativistic BW distribution is ill-equipped to handle resonances lying close to the threshold. On the other hand, even though the relativistic BW distribution can be modified to include a threshold, the latter has to be artificially imposed. This limitation of the BW distributions appear stark in case of resonances like the a 1 (1260) as well as for the later discussed a 0 (980).
The plots of the fits of the three distributions (BW, rBW, and Sill) discussed in the text for the ALEPH data for the a 1 (1260) are shown in Fig. 3 along with the experimental data [67]. The fitted mass and the width of the a 1 (1260), their error estimates, and the pole position are listed in Table 3. The table shows that both the BW distributions -relativistic and nonrelativistic, over-estimate the position of the peak compared to the experimental data. As indicated by the reduced χ 2 , the Sill distribution estimates the peak and the width of the spectral function better than the BW distributions.
The spectral function of a 1 (1260) poses a special problem as the peak is close to the ρπ threshold. In the region close to the threshold, the BW and the rBW distributions overestimate the spectral function. As in the case of the ρ, the Sill can describe data fairly well even if it neglects the proper angular momentum. This fact confirms that the the inclusion of this effect, while surely relevant for a detailed description that goes beyond the simple Sill, is subleading w.r.t the threshold effect.
The Sill distribution cannot reproduce the data close to threshold, because of the imposed threshold m ρ + m π . This drawback can be overcome by convoluting the spectral function of the a 1 (1260) with the spectral function of the ρ(770), see Sec. 4.3 for details. The plot comparing the convoluted Sill and the plain Sill with the experimental data are shown in Fig. 4. It is visible that the convoluted Sill can qualitatively reproduce the behavior of the spectral function near the threshold while still roughly fitting the data elsewhere. Note, however, that in this case the convoluted Sill has not been fitted to the experimental data and instead the parameters have been adjusted to fit the data points. The new values of the relevant parameters are: g a 1 = 1.48 MeV, Γ a 1 = 625 MeV, Γ ρ = 130 MeV. The rest of the parameters take the values given in Tables 1, 2 and 3. In the future, it would be advisable to fit directly the convoluted Sill to data, as this would allow to describe the data points close to the m ρ + m π as well as close to the peak. Here, we did not include a possible direct 3-pion decay of the a 1 meson.  Figure 3: The spectral function for the a 1 (1260). Experimental data corresponds to the spectral function of the τ → πππν τ decay reported in [67]. It is visible that the Sill fares better than the (r)BW distribution.

The K * (892)
The plot involving the spectral functions BW, rBW and Sill for the meson K * (892), fitted using the experimental data [69] (similar experimental data can be found in [70,71]) is given in Fig.  5. The mass and the width used to fit the data, their error estimates, and the pole position are listed in Table 4. All three distributions fit the data with acceptable (and similar) accuracy. As mentioned earlier, the difference between Sill and BW and rBW distributions is negligible when Γ M . This feature is evident in the case of K * (892), in which the Kπ threshold is sufficiently far away from the peak to have a sizable impact on it.  Table 4: Mass and width of K * (892) fitted using the three distributions discussed in the text, their error estimates, and the poles (as described in the text). It must be stressed that the convoluted Sill is not fitted to data (see text for details) and hence, the plain Sill fairs better in the high energy region. Yet, the convoluted Sill shows that events below the ρπ threshold are possible.

Baryonic resonances: The ∆(1232)
We now provide the example of the ∆(1232) resonance to demonstrate the applicability of the Sill distribution to the baryonic resonances. The ∆(1232) is the lightest baryonic resonances with isospin 3/2. It has been studied extensively in the photoproduction (N γ) and π-p scattering experiments. Its dominant mode of decay is into a nucleon and a pion. The PDG lists the mass and width of these resonance as 1232 MeV and 117 MeV respectively, and the Breit-Wigner pole at √ s = (1210 − i50) MeV. The BW, rBW, and Sill distributions were fitted to the spectral function of the ∆(1232) extracted from the πp scattering cross section [72]. The plots of the distributions are shown in Fig. 6, and the values of the parameters are listed in Table 5. The experimental data plotted in Fig. 6 corresponds to sin 2 (δ(E)), where δ(E) is the πp scattering phase shift and this quantity is proportional to the spectral function. The error involved in this calculation has not been reported, hence we have taken the lowest value to be the estimate of the error. In addition, some remarks of caution are in order: the πp scattering process in this energy region and angular momentum -isospin channel (P 33 ) gets a significant background from the ∆(1600) which has a mass of around 1510 MeV, and width 270 MeV for which N π is a possible decay channel. The presence of this background can cause the spectral function to become asymmetric. We have not modelled this non-trivial background. Further, there exist non-resonant background effects which we have not taken into account (see Refs. [73,74,75] for some details).
The mass of the ∆(1232) is close to the πp threshold. However, since the width of the ∆(1232) resonance is small, the threshold does not overlap significantly with the resonance, contrary to the case of the a 1 (1260). Thus, the effects of the threshold are less pronounced on the ∆(1232) resonance compared to the a 1 (1260). Yet, tt can be seen from the Table 5 that s GeV th Κ Figure 5: The spectral function for the K * (892). Experimental data (counts) from [70]. All three functions fair similarly well.
the Sill fits better than the BW and the rBW to the data. The mass, width and the position of the pole estimated by the BW, rBW, and the Sill distributions are nearly the same. However, the Sill estimation is more accurate than the (r)BW estimate as shown by the χ 2 /d.o.f. Also, an inspection of the plots given in Fig. 6 shows that although the three distributions fare identically around the peak, the Sill fits the data better closer to the threshold.  Table 5: Mass and width of ∆(1232) fitted using the three distributions discussed in the text, their error estimates, and the poles (as described in the text).

The multi-channel sill distribution: The f 2 (1270)
The multi-channel Sill distribution was discussed in sec. 4.2. Here, we demonstrate the usefulness of the Sill distribution taking the case of f 2 (1270) as an example. The well-known resonance f 2 (1270) decays mostly into ππ,KK, and 4π states. Here, we aim to study it using the Sill in order to show how to deal with multichannel decay as well as with the convolution leading to the four-pion decay mode. This resonance is often understood as a standard quark-antiquark state belonging to a nonet of ground-state tensor mesons (yet other interpretations are possible, see Ref. [76]).
The ππ channel is the dominant channel and accounts for nearly 84.2% of the decays [1]. We assume that the 4π decay proceeds through intermediate ρρ i.e, f 2 (1270) → ρρ → 4π. This is a suppressed process because the sum of the masses of the intermediate ρ mesons is larger than the mass of the f 2 (1270).  Figure 6: The spectral function for the ∆(1232). Experimental data from [72]. It is visible that the Sill fairs marginally better than the (r)BW distributions.
Out of the two observed 4π decays, the 2π + 2π − , which within our framework arises from the ρ 0 ρ 0 decay, is a factor two smaller than the π + π − 2π 0 decay, which comes from 2ρ + ρ − (this is a consequence of isospin counting). We do not distinguish between the two channels, but we note that the experimental result is consistent with our assumptions. The fraction of the 4π decay is the sum of the fractions of the two possible 4π decay channels and is ∼ 10.5%. With this simplification, the 4π channel becomes the second most favored channel for the decay of f 2 (1270).
The last channel we have considered is theKK channel. This channel is also suppressed because of the highKK threshold of ∼ 980MeV. The plot of the (convoluted) Sill distribution and the different channels is shown in Fig. 7. The list of the employed parameters is given in Table 6. The values ofΓ i are calculated using Eq. (57) with the partial decay widths calculated from the branching ratios. For the 4π channel, the 2ρ threshold is higher than the mass of the f 2 (1270). Hence, Eq. (57) cannot be used. Instead, we use a value such that the branching ratio after convolution matches the value given in PDG [1]. The convoluted 4π channel (inset, Fig. 7) stands out because of the difference in the position of the maximum. The peak of the ρρ mediated 4π channel appears at ∼ 1.62 GeV, in contrast with the mass of f 2 which is ∼ 1.275GeV. However, the effect of this channel is less pronounced as its contribution is small. Further, the convolution pushes the peak of the f 2 to a slightly higher value of 1.283 GeV. The pole and the branching ratios extracted from the convoluted sill are listed in Table 7. The values branching ratios listed here are in good agreement with the PDG.
In conclusion, the multi-channel treatment of the f 2 (1270) via the Sill offers an interesting, albeit simple, description which is valid in first approximation. Of course, a more detailed and correct treatment needs to make use of the correct QFT involving tensor mesons. The four channels considered are the ππ channel, theKK channel, and the two 4π channels. The convoluted 4π channel, whose peak is shifted to higher energies, is shown in the inset.   Table 7: The pole and the branching ratios of the three decay channels of f 2 (1270) extracted from the Sill distribution.

Further examples and limitations
In this subsection, we present some additional examples which show the limitations of the Sill.

The a 0 (980)
The resonance a 0 (980) is peculiar, since it was used to introduce for the first time the well-known Flatté distribution, which is generalized by the Sill in this work. Moreover, the kaon-kaon mode of this resonance opens just at the mass of this resonance, thus making it an interesting case for a multichannel study. On the other hand, care is needed, since the resonance a 0 (980) has been studied in a variety of works as a exotic meson candidate (e.g. in the form of a molecular state or a companion pole) (e.g, [65,77,78,79,80], for a review on the status of exotic hadrons, see [81]). We thus proceed as follows. First, we discuss the a 0 (980) as a single normalized state, ignoring its nontrivial nature, and later on we include modification/extension due to its nonconventional origin.
The a 0 (980) decays predominantly into the ηπ and theKK channels. As discussed in Ref. [82], the ηπ channel is the dominant decay channel for a 0 (980). TheKK channel is special here, as the mass of the resonance lies very close to theKK threshold. The branching ratios of the decay B 0 → ψa 0 (980) →KK is very much sensitive to the mass of the a 0 (980) [83].
The multi-channel Sill is capable of explaining some of the peculiar features of the shape of its spectral function. As numerical inputs, we use the branching ratios and the width at half-maximum given in the PDG [1] to estimate the parameters listed in Table 8 employed for plotting Fig. 8. While this study should be only regarded only as heuristic, it represents a useful demonstration of a channel opening right at the threshold. In fact, this threshold effect is visible in the spectral function as well as in both the decay channels at √ s = 2m K . For the parameters given in Tabe Table 8: Parameters of the Sill for a 0 (980): the mass and the rescaled partial widths of the ηπ andKK channels used in plotting Fig. 8.
Yet, as mentioned above, the a 0 (980) is nowadays regarded as a non-conventional dynamically generated mesonic state. In this context, the state associated with a 0 (980) does not need to be normalized to unity (as naively done in Fig. 8), as approaches in which this state is interpreted a kaon-antikaon bound state show [79]. Moreover, a similar conclusion is obtained in Refs. [34,85,86], in which a 0 (980) is described as a companion pole of the predominantly quark-antiquark state a 0 (1450). In both pictures, the normalization of a 0 (980) amounts to a number Z < 1.
In this respect, the question is how to modify the Sill in order to take into account the more complicated nature of both resonances a 0 (980) and a 0 (1450). We thus write down a normalized two-Sill as a simple model for two-state system in which one is a companion pole of the other as: where the superscripts low and high refer to the a 0 (980) and the a 0 (1450) respectively and where Z (with 0 ≤ Z ≤ 1) is a weight factor. Thus, the two-Sill is a simple parametrization aimed to capture some of the features of rather complex phenomenon that involve hadronic loops. More specifically, in the large-N c limit one expects a unique quark-antiquark state |qq , whose mass is about 1.3 GeV and corresponds (roughly) with a 0 (1450). Then, the interaction of the original seed state with other mesons induces a mixing that adds new components to the Fock space that are proportional to K K and |πη . Thus, the original state is a mixture of both quark-antiquark and meson-meson component. This mechanism is quite general and takes place also for standard resonances, such as the already mentioned ρ-meson. In the specific case of the a 0 -system, the mixing can be strong enough to generate an additional pole that is identified to a 0 (980) (the κ − K * 0 (1430) system is also similar, as discussed in the next subsection), which, as a consequence, can be seen as predominantly a meson-meson state.
The values of these poles for a 0 (980) are: and for the a 0 (1450): Note, BP provide the pole mass and pole width as M pole = 1082 MeV and Γ pole = 309 MeV. For the a 0 (1450), we have used the BW mass and width given in the PDG [1].
The constant Z is taken to be 0.6, 0.5, and 0.7 for the TR, WGR, and the BP curves respectively, in such a way to mimic those corresponding calculations (see the plots for all 3 cases in Ref. [85]).
The two-Sill distributions in Fig. 9 are actually pretty similar to those of TR, BP, and WGR, showing that a qualitative agreement is offered by our simple approach. Clearly, it must be stressed again that the two-Sill does not represent a substitute of a full loop calculation, but it shows that the presented modification of the Sill offers a potentially useful way to study, in the first approximation, some non-conventional resonances.
As a final remark, we underline that the companion pole scenario presented above is only one framework among the possibilities to understand the dynamically generated states. Other relevant approaches exist, for instance the already mentioned molecular picture for the a 0 (980) of Ref. [79], which is based on the compositeness condition [87,88]. Quite interestingly, in Ref. [79] the integral over the spectral function between 2m K − 0.05 GeV and 2m K + 0.05 GeV ranges between 0.2-0.5. For comparison, a standard BW spectral function with a width of 50 MeV would deliver (2/π) arctan(2) ∼ 0.7. The fact that obtained value (0.2-0.5) is smaller than 0.7 was interpreted as a sign of an important molecular component in this state (an even smaller value is found for the f 0 (980)). A similar integral for the three curves of Fig. 9 delivers the values 0.05 (TR), 0.1 (BP), and 0.25 (WGR), the last of which is consistent with [79]. Surely, this result is built-in by the form of Eq. (96), which is inspired by previous works on the subject. Nevertheless, as a general remark, the lack of normalization of a certain resonance around its peak can be seen as a hint of its non-standard nature.
5.6.2 The κ and K * 0 (1430) The next question is if the description mentioned above for the a 0 (980) can be applied also to the resonance κ, which does not even show a peak in the cross-section data.
Indeed, also here the normalization of κ is not a requirement, since this state is also (predominantly) not quark-antiquark mesons, see e.g. Refs. [47,89].
As an exercise, we try to apply the two-Sill to this case as well. We thus consider κ as the companion pole of K * 0 (1430) and write down: where in this case low refers to κ and high to K * 0 (1430). The pole of the κ is listed by the PDG as √ s pole = (630 − 730) − i(260 − 340) MeV. We adjust the mass and the width of the d low s (E) to get the pole at √ s pole ≈ 700 − i295 MeV. The PDG does not list a value for the pole of K * 0 (1430). We thus adjusted the mass and width of the K * 0 (1430) to obtain the peak seen in the data, which leads to a pole at √ s pole = 1345 − i88 MeV. The weight factor Z takes the value 0.5. The masses and widths so obtained were used to calculate the effective distribution similar to in Eq. 96 and the resultant curve is plotted in Fig. 10, along with the data from the πK scattering experiment [90]. The plotted experimental values correspond to sin 2 (δ πK (E)), where  Table 9: Illustrative list of some mesonic resonances near the respective thresholds. ∆m is the difference between the masses of the initial and final states.
δ πK (E) is the πK scattering phase shift. Also plotted in the Fig. 10 is the spectral function of the κ. It is visible that the two-Sill is peaked at about 1.3 GeV and the light κ represents only an enhancement, as more advanced model show. Yet, it is not possible to get a qualitative agreement if we impose that the pole of the light κ agrees with the PDG range. (A better agreement could be obtained at the price of a pole at higher values, but this is in conflict with many works on the subject).
We thus encounter a clear quantitative limitation of the double-Sill for this specific case. This is somewhat expected: the light κ is a pole that is caused by pion-kaon loop effects, and the (double-)Sill is not capable of reproducing the details of this system.
For what concerns the f 0 (500), the Sill that corresponds to the pole quoted by the PDG has a similar form as that of the light κ reported in Fig. 10, but is even more elongated. This is in agreement with the non-conventional nature of this state. A more detailed study of this sector would necessitate even more than two resonances, but the situation is at least as complex as the light κ: one does not expect to get a correct description of data for such a non-perturbative system using the Sill.

Applicability of the Sill
In this subsection, we discuss under which conditions one may expect the proposed Sill distribution to be -in first approximation-useful.   An important feature of the Sill distribution is its (simple but consistent) treatment of the thresholds. As the various examples in Sec. 5 have shown, the Sill works better than (r)BW when dealing with the states/resonances lying close to the threshold. In particular, denoting by ∆m as the distance to a certain threshold, we expect some advantage of using the Sill when Γ ∼ ∆m. This condition is indeed fulfilled for the a 1 -meson with Γ ∼ 300 MeV and ∆m ∼ 300 MeV and for the ∆-baryon with Γ = 117 MeV and ∆m = 157 MeV. Even in the case of the ρmeson, for which Γ ∼ 150 MeV is three times smaller than ∆m ∼ 500 MeV (thus Γ/∆m ∼ 1/3) the fit with the Sill is favored. On the other hand, for K * (892) (for which Γ ∼ 50 MeV and ∆m ∼ 250 MeV, resulting in a value Γ/∆m ∼ 1/5) the Sill and the (r)BW deliver very similar results, thus in this case the use of the Sill does not provide additional advantages. Out of these considerations, it seems that a ratio Γ/∆m 1/3 is required for the Sill to be advantageous. For values smaller than this ratio we do not expect the Sill to be better than (r)BW.
Interestingly, states that have masses lying close to significant thresholds appear quite often. Some of such states have been listed in Table 9 and Table 10 along with their masses, decay widths, prominent decay channels whose threshold is close to the mass of the parent, the distance to the threshold, branching ratios (when known), and the 3-momentum carried by the decay products.
In the meson sector, notable examples are the four heavy mesons ψ(3770), ψ(4040), χ c1 (4140), and Υ(4S). For instance, the ψ(3770) and the Υ(4S) are very close to the DD and the BB threshold respectively [1]. As one can see from the Table 9, the difference between the masses of the parent and the decay products is less than or very close to the decay widths of the parent states. In this Table, we include also the hybrid candidate π 1 (1600) (note, as an hybrid, this state is not dynamically generated, but a confined state that survives in the large-N c limit). For this resonance, the lattice calculations show that the b 1 π channel is the dominant channel [91].
Also among baryons there are resonances which have one of the dominant modes of decay as sub-threshold. One such example is the N (1720). The N (1720) is a 3/2 + baryon with the two dominant decay channels: N ππ and N ω. Of these, the N ω decay channel is sub-threshold with the mass difference between −42 MeV and 28 MeV. Other resonances, like the ∆(1600), the Λ(1600), and the N (1535) lie very close to the thresholds of their dominant decay channels, as displayed in Table 10.
The closeness of these states to their respective thresholds is expected to alter their line shapes, making the Sill as a valid alternative to the conventional relativistic Breit-Wigner distribution functions. Of course, the Sill is a quite general parametrization with some specific interesting features, that however does not represent a substitute to more complete studies that take into account the properties related to the specific quantum numbers as well a more advanced inclusion of loop effects of a given state.
The resonances listed in the Tables can be interpreted as conventional quark-antiquark mesons or three-quark baryons (with the exception of the already mentioned hybrid π 1 (1600) resonance). In case of dynamically generated states, care is needed, as we discuss in detail in Sec. 5.6. One may still use the Sill -e.g. in the framework of the companion pole scenario-but the results can be at best qualitative.
Summarizing, the Tables 9 and 10 should be seen as illustrative examples (beyond Sec. 5) that show that threshold effects are pretty common. Yet, for those states, that are analogous to the resonances discussed in Sec. 5, we may expect the Sill to be useful. In the future, whenever a certain resonance sufficiently close to a threshold is found, together with the usual (r)BW, we think that a fitting via the Sill could be advisable.
As a final remark, we mention also that the Sill can be employed as toy model. The fact that the normalization of the spectral function is fulfilled also for broad states close to the thresholds, make it an interesting tool for testing its properties that appear in certain hadronic models, for instance, when calculating decays by taking into account finite width effects or as an intermediate state in three-body decay.

Conclusions
In this work, after a brief recall of the non-relativistic and relativistic Breit-Wigner functions and their justifications as certain limiting cases, we have presented a relativistic and properly normalized distribution -called Sill distribution-which models the presence of energy threshold(s). It can be in principle applied to any unstable state, but the realm of hadrons seems the best place to use such a spectral function, since hadrons are often quite broad and various thresholds are present. Indeed, the form of the Sill is inspired by the way the kaon-kaon decay mode was included in the so-called Flatté distribution for the resonance a 0 (980). The form of the self-energy is applied to any decay channel (even if the masses of the decay products are different), regardless of the particular angular momenta involved in the decay.
The determination of the pole is also straightforward within the Sill distribution. Moreover, the extension to a the multichannel case and to a convolution suitable for decay chains have been also developed.
A direct comparison with a scalar QFT shows under which conditions the Sill function emerges. Moreover, a direct application of the Sill to the well-known resonances ρ and a 1 (1260) shows that, despite its simplicity, it works pretty well. This is especially evident for the latter case, due to the fact that a 1 (1260) is quite broad. A similar, even if less pronounced, improvement has been achieved for the baryonic state ∆(1232), showing that the Sill can be also applied to baryonic resonances. On the other hand, for the relatively narrow resonance such as K * (892), the Sill distribution fairs as well as the standard BW and rBW functions. We have also shown a direct implementation of the multichannel and convoluted sill for the well established resonances f 2 (1270).
Next, we have also discussed the example of non-conventional mesons, such as the resonance a 0 (980). A possible simple description is to extend the Sill to describe the a 0 (980)-a 0 (1450) system. Yet, in doing so we have encountered also limitations of the Sill, that -as expectedbecome evident also in the case of the very broad and unconventional light κ meson.
In conclusion, we have presented a simple and general distribution that can be applied for various resonances, including broad ones affected by threshold effect(s). Since the expected effort of using such a distribution in fits to data w.r.t. the standard relativistic Breit-Wigner seems comparable and in consideration of the relative simplicity and theoretical foundation of behind it, which includes the normalization necessary for a probability distribution, we regard the Sill as an additional useful tool to describe some resonances.
Acknowledgments: The authors thank M. Rybczynski for useful discussions. F. G. and V. S. acknowledge financial support through the Polish National Science Centre (NCN) via the OPUS project 2019/33/B/ST2/00613.

A Alternative relativistic Breit-Wigner distribution
In this appendix we discuss a relativistic by keeping in which we keep the energy dependence in the numerator. The corresponding function, denoted as rBWE, reads: This expression differs from the one given in Eq. (45) in the numerator and the denominator. This alters the form of the distribution slightly. For the purpose of demonstration, we fit Eq.
(98) to the spectral function of the a 1 (1260). The χ 2 fit returns us the following parameters:  Table 3. We see that, the mass of the a 1 (1260) estimated by rBWE is significantly larger than that estimated by rBW and smaller than the value estimated by the Sill. The real part of the pole is larger than that estimated by rBW and the Sill, and the imaginary part lies in between the values estimated by the rBW and the Sill. We display in Fig. 11 the plots of the fits of rBWE along with all the other forms the distribution functions discussed in the text and the experimental data from ALEPH [67]. From the plots and the parameters, we see that the two versions of the rBW distributions provide very similar fits to the data, and estimated values of the mass and width are close to those estimated by the nonrelativistic BW. The fits themselves are better than that of the nonrelativistic BW as evidenced by the χ 2 /d.o.f value but worse than the fit of the Sill distribution.
On the other hand, the spectral function of the ρ(770), as estimated by the 4 distributions appear similar, as shown in Fig. 12. However, the mass, width, and the pole positions estimated by the rBWE are signifcantly different from that estimated by the rBW or the Sill. The rBWE estimates for these quantities are: M = 756.1 ± 0.30 MeV, Γ = 139.0 ± 1.2 MeV, and √ s pole = 759.2 − i69.2 MeV, with a χ 2 /d.o.f of 9.42. The mass and width are closer to the values estimated by the Sill, but the position of the pole is closer to the value estimated by rBW. This estimate is as accurate as the rBW estimate as per the value of χ 2 /d.o.f. Since the ρ(770) peak is far away from the 2π threshold, one cannot expect the Sill to show a large deviation from the rBW or the rBWE. However, according to the statistics, the Sill fits the resonance better than the rBW (or the rBWE).
In general, we may notice that there is not a unique rBW form.

B Hesse matrix and error estimation
The errors in the parameters have been calculated using the Hesse matrix. Here, we provide a brief description of the Hesse matrix and error estimation. Consider a set of data that is expected to be described by some theoretical functions f i , each of them depending on the N parameters {x 1 , x 2 , . . . , x N }. The goodness of the fit can be measured using the χ 2 defined as: where the f exp,i are the measured values of the function f i and the δf exp,i are the error in the measurements. The parameters can be estimated by demanding that the value of χ 2 is minimum. The Hesse matrix (H) is defined as the matrix of the second derivatives of the χ 2 , H ij = 1 2 ∂ 2 χ 2 (x 1 , x 2 , . . . , x N ) ∂x i ∂x j P (100) where P represents the values of the parameters for which χ 2 is at its minimum. With this definition of the Hesse matrix, the errors in the parameters are given by