What is the right formalism to search for resonances?

Hadron decay chains constitute one of the main sources of information on the QCD spectrum. We discuss the differences between several partial wave analysis formalisms used in the literature to build the amplitudes. We match the helicity amplitudes to the covariant tensor basis. Hereby, we pay attention to the analytical properties of the amplitudes and separate singularities of kinematical and dynamical nature. We study the analytical properties of the spin-orbit (LS) formalism, and some of the covariant tensor approaches. In particular, we explicitly build the amplitudes for the B→ψπK\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B\rightarrow \psi \pi K$$\end{document} and B→D¯ππ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B\rightarrow \bar{D}\pi \pi $$\end{document} decays, and show that the energy dependence of the covariant approach is model dependent. We also show that the usual recursive construction of covariant tensors explicitly violates crossing symmetry, which would lead to different resonance parameters extracted from scattering and decay processes.


Introduction
The high quality data on hadron production and decays that are or will be available from BaBar, BelleII, BESIII, CMS, CLAS12, COMPASS, GlueX, LHCb, and other experiments, necessitate rigorous amplitude analysis. This is particularly true for the extraction of resonance parameters that are based on analytical partial waves. Moreover, analytical reaction amplitudes are needed in conjunction with lattice data to a e-mail: mikhail.mikhasenko@hiskp.uni-bonn.de b e-mail: pillaus@jlab.org study the hadron spectrum from first-principles lattice QCD calculations [1][2][3][4][5].
In this paper, we focus on three-body decays, aka 1-to-3 processes. In recent years such reactions have led to ample data that resulted in the observation of new exotic phenomena, e.g. the so-called XYZ states in heavy meson decays [6][7][8], and that are also used in studies of excited mesons and baryons. The issues we address and the methodology we present are, however, of relevance to other analyses as well, for example to baryon resonance studies in photoproduction [9,10], or meson spectroscopy from pion or photon beam fragmentation [11][12][13].
In the modern literature, there seems to be a lot of confusion regarding properties of the reaction amplitudes employed in analyses of such processes. This is often stated in the context of a potentially nonrelativistic character of certain approaches [9,14,15]. As we explain below, however, rather than arising from relativistic kinematics, the differences between the various formalisms have a dynamical origin. Reaction amplitudes are given by the scattering matrix elements between initial and final states that represent asymptotically free particles. Such states belong to a unitary, noncovariant representation of the Lorentz group. Since the scattering operator is a Lorentz scalar, reaction amplitudes share the transformation properties of the free particle states. A typical three particle decay process is dominated by two-body resonances, and can be well approximated by a finite number of partial waves. The latter can be given by the helicity partial waves or the Russell-Saunders, aka LS amplitudes [16]. For the LS amplitudes, one couples particle states in the canonical representation. The relation between the helicity and LS basis is a straightforward orthogonal transformation. Because of the noncovariant transformation properties of the reaction amplitude, partial waves transform in a nontrivial way as well, e.g. helicity amplitudes mix under Lorentz boosts through Wigner rotations. Nevertheless, all of the amplitudes referred to above (the helicity amplitudes, the helicity partial waves, the LS partial wave amplitudes) are relativistic, i.e. have well defined behavior under Lorentz transformations.
Since the helicity amplitudes involve asymptotically free particle states, they must be proportional to free particle wave functions, e.g. Dirac spinors or polarization tensors. These wave functions have mixed transformation properties, i.e. have both covariant (Lorentz or Dirac), and noncovariant (helicity) indices. The Lorentz and Dirac indices need to be contracted with covariant tensors built from particle four-vectors and Dirac gamma matrices to yield the noncovariant helicity amplitudes. Helicity amplitudes can therefore be expressed as linear combinations of products of covariant tensors and wave functions with coefficients that are scalar functions of the Mandelstam invariants. It can be shown that these scalar functions have only dynamical singularities as demanded by unitarity [17], and for this reason are useful when analyzing singularities of the partial waves. Furthermore, these scalar functions are invariant under crossing which makes them convenient to relate amplitudes in the decay and scattering kinematics.
There exist an approach for constructing the scalar functions from an assumed model for the partial waves, hereafter referred to as the covariant projection method (CPM) [14,[18][19][20], that starts from a LS partial wave model (or equivalently the Cartesian, aka Zemach amplitudes [21]) but writes them in a covariant fashion. The method has a drawback, which is related to the behavior under crossing (see Sect. 3.1). The alternative, which we refer to as the canonical approach [16,[22][23][24][25], is to use the well known relation between the helicity amplitudes and the helicity partial waves [16] to determine the scalar functions in terms of the partial wave models. The differences between these two approaches to relate partial waves and scalar functions result in factors which are confusingly referred in the literature as "relativistic corrections". These are actually Lorentz invariant functions and therefore can be absorbed into the scalar functions. In both the CPM and canonical approaches, the relativistic kinematics is properly taken into account. Thus, the differences in these approaches are dynamical in nature.
In what follows, we present a detailed comparison of these two approaches, paying specific attention to the analytical properties, which are among the few constraints that can be imposed in a model independent way. Instead of presenting results for a general case, we find it more pedagogical to compare these constructions in a few concrete examples. The examples we discuss are of special interest to various ongoing analyses, and are complex enough to illustrate the general principles. The first example is the parity violating (PV) three-body decay B 0 → ψπ − K + , with ψ = J/ψ, ψ(2S). The analyses by Belle and LHCb show nontrivial structures appearing in the ψ(2S) π [25][26][27][28], and in the J/ψ π channel [29]. These are of particular interest, because a resonance in these channels would require an exotic interpretation [6][7][8]. The rest of the paper is organized as follows. In Sect. 2 we discuss the canonical approach on the example of the B → ψπ K decay. By relating the helicity partial waves to the scalar amplitudes via the partial wave expansion, we derive constraints and isolate the kinematical singularities. We also discuss implication of these constraints for the LS partial wave amplitudes. The details of the amplitude parameterizations are given in the Appendices and are presented in a way that can be implemented in the standard data analysis tools [30,31]. In Sect. 3 we examine the CPM approach and compare this model with the findings from Sect. 2. We mention the crossing symmetry properties of CPM using, as an example, B 0 →D 0 π + π − , which was recently analyzed by LHCb within this formalism [32]. Summary and conclusions are given in Sect. 4.

Analyticity constraints for B → ψπ K
In Fig. 1 we show a diagram representing the kinematics of the decay B → ψ(→ μ + μ − )π K . The spinless particles B, π , K are stable against the strong interaction. The ψ is narrow enough to completely factorize its decay dynamics. Thus, we construct the amplitude considering ψ to be stable. More details, including the dilepton decay of the ψ, are given in Appendices A and B. We use p i , i = 2 . . . 4 to label the momenta of B, π , and K respectively. The momentum of the ψ will be denoted byp 1 , for a reason which we will explain below. The helicity amplitude for the decay process B → ψπ K is denoted by A λ (s, t), λ being the helicity of ψ, i.e. ψπ K , out|B, in = (2π) 4 The amplitude depends on the standard Mandelstam variables s = ( p 3 + p 4 ) 2 , t = (p 1 + p 3 ) 2 , and u = (p 1 + p 4 ) 2 with s + t + u = i m 2 i . The B meson decays weakly, so A λ is given by the sum of a PV and parity conserving (PC) amplitude. The difficulty treating the decay channel directly is that the mass of the decaying particle should be considered on the same footing as the other dynamical variables (s, t, u). This is demanded by unitarity, which implies that above a threshold, the amplitude is a singular function of the corresponding dynamical variable. It is therefore simpler to study singularities in a scattering channel and cross to the other channels by analytical continuation in the momentum of the ψ, i.e. by settinḡ p 1 = −p 1 [33]. In general, under crossing, helicity amplitudes are mixed by Wigner rotations. In our case, however, since crossing can be realized through a (unphysical) boost in the direction of motion of the ψ, there is no change in helicity.
We begin with the discussion of the PV amplitudes in the s-channel. The s-channel resonances correspond to the K * 's and dominate the reaction. As discussed in the previous section, the analysis of the experimental data indicates a possible signal of resonances in the exotic ψπ spectrum, which in our notation correspond to the t-channel. Once we have constructed the s-channel amplitudes, the t-channel ones can be treated similarly (cf. Appendix B).
In the center of mass of the s-channel scattering process, the ψ momentum defines the z-axis, the momenta p 3 and p 4 lie in the xz-plane. We call p (q) to the magnitude of the incoming (outgoing) three momentum. The scattering angle θ s is a polar angle of the pion (see Fig. 2). The quantities depend on the Mandelstam invariants through

Fig. 2
Scattering kinematics in the s-channel rest frame. In the decay kinematics, the momentum and the spin of the ψ is reversed, so to keep the same helicity To incorporate resonances in the π K system with a certain spin j, we expand the amplitude in partial waves, where A j λ (s) are the helicity partial wave amplitudes in the schannel. In Eq. (2) the entire t dependence enters though the d functions. The d functions have singularities in z s which lead to kinematical singularities in t of the helicity amplitudes A λ . The dynamical singularities in t, related to, for example, the possible resonances in the ψπ channel, can only be reproduced if the the sum contains the infinite number of partial waves. In practice the t-or u-channel resonances (singularities) are accounted for explicitly through t-or u-channel partial waves, and to avoid double counting each series is truncated at a finite number of terms. This defines the socalled isobar model in which where J max is finite. However, we remark that the analysis of kinematical singularities has general validity, and might be applied to the original untruncated series. The expressions for the (t) and (u) isobars are similar to Eq. (4). Note, that due to the superscript (s) the amplitudes A (s) j λ (s) are not identical to the helicity partial waves, A j λ (s) of Eq. (2). This is because the other two terms on the right hand side of Eq. (3) also contribute to the s-channel partial wave expansion. We refer to the former as the isobar partial waves or simply, isobars. The difference between the partial waves, which are defined in a model independent way, and isobars, which appear in the specific model as in Eq. (3), has important consequences when establishing the relation between phases of the isobar amplitudes and those of the partial waves [34][35][36][37][38][39][40]. This issue, however, is not directly related to the topic of this paper and we do not discuss it any further.
We return to the partial wave expansion, and proceed with the analysis of kinematical singularities. An extensive discussion and the full characterization of these singularities can be found in [16,[41][42][43][44][45]. We recall that d j is the so-called half angle factor that contains all the kinematical singularities in t. The reduced rotational function d j λ0 (z s ) is a polynomial in s and t of order j − |λ| divided by the factor λ . The helicity partial waves A j λ (s) have singularities in s. These have both dynamical and kinematical origin. The former arise, for example, from schannel resonances. The kinematical singularities, just like the t-dependent kinematical singularities, arise because of external particle spin. We explicitly isolate the kinematic factors in s, and denote the kinematical singularity-free helicity partial wave amplitudes byÂ j λ (s). 1 First, the term ( pq) j−|λ| is factorized out from the helicity amplitude A j λ (s). This factor is there to cancel the threshold and pseudothreshold singularities in s that appear ind j λ0 (z s ). Second, we follow [16] and introduce the additional kinematic factor K λ0 ('±' is short for λ = ±1). These factors are required to account for a mismatch between the j and L dependence in the angular momentum barrier factors in presence of particles with spin. Finally, the kinematical singularity-free helicity partial wave amplitudesÂ j λ (s) are defined by with K 00 and K ±0 given by Specifically, it is expected that A j λ (s) ∼ p L 1 q L 2 at threshold, where L 1 and L 2 are the lowest possible orbital angular momenta in the given helicity and parity combination. This explains why j = 0 requires a special treatment [45], since for j ≥ 1 we have L 1 = j − 1, but for j = 0 the lowest is L 1 = j + 1. In addition, the K -factors have powers of √ s as required to ensure factorization of the vertices of Regge poles. Similarly, as explained before, m 1 is dynamical and thus the kinematical singularity-free amplitudes are not expected to contain singularities in m 2 1 , and as will be seen below, the m 1 dependence of the K -factor takes care of that. TheÂ j λ (s) are left as dynamic functions, which are unknown in general and cannot be calculated from first principles. Usually they are parameterized in terms of a sum of Breit-Wigner amplitudes with Blatt-Weisskopf barrier factors. 1 For fermion-boson scattering, theÂ j λ (s) can still have a branch point at s = 0, as discussed in [44].
We now seek a representation of A λ (s, t) in terms of the scalar functions, as discussed in Sect. 1. For the PV amplitude this is given by Although the second term in brackets may look like an extra 1/s pole, it cancels when multiplied by ( p 3 + p 4 ) μ . This choice simplifies the final expressions, but we remark that any other choice of independent tensor structures would lead to the same results. In the s-channel center of mass frame the ψ polarization vectors are given by The energies E i are calculated from the momenta and are fully determined by s. The functions B(s, t) and C(s, t) are the kinematical singularity free scalar amplitudes discussed in the Sect. 1.
We can match Eqs. (2) and (6), and express the scalar functions as a sum over kinematical singularity free helicity partial waves. The ratio A λ (s, t)/ (K λ0 ξ λ0 (z s )), computed using Eq. (6), is compared to the same ratios computed using the helicity partial waves from Eq. (2). This yields from λ = ± and λ = 0, respectively, which can be combined to obtain Neither B(s, t) nor C(s, t) can have kinematical singularities in s or t. In Eqs. (7)-(9),d j 10 (z s ) is regular in t, and the s singularities at (pseudo)thresholds are canceled by the factor ( pq) j−1 . The latter factor contains a high-order pole at s = 0. Such pole is a feature of the dynamical model, and specifically arises because at s = 0 the little group is not SO(3) anymore. The latter motivates the partial wave expansion, thus it is not surprising that the truncation of the partial wave series results in such singularities [16,46]. 2 This construction hence does not constrain the poles at s = 0.
For the same reason the sum in Eq. (9) has no kinematical singularities in s and t, however the 1/λ 12 factor in front of the sum generates two poles at s ± = (m 1 ± m 2 ) 2 , unless the expression in brackets vanishes at those points. This means that theÂ j λ (s) with different λ cannot be independent functions at the (pseudo)threshold. Explicitly, in the limit s → s ± at fixed t one has z s → ∞ and using [43], one finds that the expression within the brackets in Eq. (9) behaves aŝ This combination has to vanish to cancel the 1/λ 12 , thus one finds (for j > 0) where g j (s), f j (s), g j (s), and f j (s) are regular functions at s = s ± , and g j (s ± ) = g j (s ± ). Note that, while the functional form considered in Eq. (12) complies with the general requirements we are imposing, it actually implements more freedom than required by the former. For instance, one could take f j (s) = 0 without any loss of generality. The particular choice taken in Eq. (12), however, turns out to be useful for the comparisons with other parameterizations (LS and CPM) which we will discuss in Sect. 3. Together with Eq. (12), the expressions in Eqs. (6), (8) and (9) provide the most general parameterization of the amplitude that incorporates the minimal kinematic dependence that generates the correct kinematical singularities as required by analyticity.

Implications for the LS partial wave amplitudes
The advantage of the LS basis is that the identification of the correct threshold factors is straightforward. For a given system of two particles with spins j 1 , j 2 and corresponding helicities λ 1 , λ 2 , the relation between a two-particle state in the helicity and LS basis is where is the projection of the total angular momentum j. For the B → ψπ K amplitude, it implies the following relation between the LS amplitudes G and the helicity amplitudes, The amplitudes with L = j ± 1 and L = j differ by parity. Equation (15) can be inverted to relate the helicity partial wave amplitudes with the LS amplitudes G j l (s), In Eq. (16) we denoted the LS partial wave amplitudes with the threshold factors explicitly factored out byĜ j l (s), i.e. G j j±1 (s) = p j±1 q jĜ j j±1 (s). We now compare the general expression for the helicity partial waves with the spinorbit LS partial waves. We find that Eq. (16) matches the general form in Eq. (12) when The common lore is that the LS formalism is intrinsically nonrelativistic. However, the matching in Eq. (17) proves that the formalism is fully relativistic, but care should be taken when choosing a parameterization of the LS amplitude so that the expressions in Eq. (17) are free from kinematical singularities. For example, if one takes the functionsĜ

Comparison with the covariant projection method
We consider now the CPM approach of [14,[18][19][20]. As said, the method is based on the construction of explicitly covariant expressions. To describe the decay a → bc, we first consider the polarization tensor of each particle with index i and spin j i , i μ 1 ...μ j i ( p i ). Using the decay momentum p bc = ( p b − p c )/2 and the total momentum P bc = p b + p c , we build a tensor X μ 1 ...μ L ( p bc , P bc ) to represent the orbital angular momentum of the bc system. In order to find total angular momentum tensor, we first combine the polarizations of b and c into a "total spin" tensor S μ 1 ...μ S ( b , c ) (orthogonal to the momentum p b + p c ). Then, we combine the tensor S μ 1 ...μ S ( b , c ) with the orbital tensor X μ 1 ...μ L and finally contract the result with the polarization of a, thus mimicking the LS construction. The tensors S and X have definite spin and parity, i.e. are in an irreducible representation of the rotation group in the particle rest frame. Thus they must be symmetric, traceless, and orthogonal to the total momentum in the particles system, p b + p c . If one of the daughters is unstable, we can implement its decay in a similar way. The procedure is recursive, and relatively simple for low spins.
Together with the explicit covariance, it makes the formalism very attractive.
We use the CPM to build the amplitude for B → ψπ K . The construction of an amplitude for an arbitrary spin of the intermediate state is cumbersome, and we limit ourselves to the special case of an intermediate K * with j = 1. We start with the tensor amplitude for the scattering process ψ B → K * → π K . The orbital angular momentum of the decay K * → π K in P-wave is given by X ρ (q, P). The tensor is constructed from a four-vector of the relative momentum q μ = ( p μ 3 − p μ 4 )/2 and the total momentum of the system For the PV amplitudes, the initial process ψ B → K * is described by two waves. The corresponding orbital tensors are the unit rank-0 tensor for the S-wave and rank-2 tensor X ρμ ( p, P), where P is the K * momentum. The final P-wave orbital tensor is X ν (q, and matching with Eq. (12) gives The threshold conditions g 1 (s ± ) = g 1 (s ± ) are satisfied, and the functions f 1 (s) and f 1 (s) are regular at the thresholds. Finally, we show the relation between the CPM and the LS amplitudes. The comparison with Eq. (16) leads to Although the g S (s) and g D (s) of the CPM formalism, see Eq. (18), are typically interpreted as the S and D partial wave amplitudes, we see that this is the case only at (pseudo)threshold s = s ± , where the factor E 1 /m 1 − 1 vanishes. In Sect. 3.2 we discuss a specific example to show the numerical difference between the various approaches.

Crossing symmetry; and the decay B →Dππ
An issue with the CPM formalism is the explicit violation of crossing symmetry. The recursive procedure explained in [14,[18][19][20] produces different scalar amplitudes if applied in the scattering or in the decay kinematics. For the decay kinematics, the CPM amplitude is constructed according to a chain B → ψ K * (→ K π). The tensor X ν (q, P) describes the P-wave decay K * → K π as before. S-and D-waves are still possible for the decay B → ψ K * . The same symbolic expression in Eq. (18) holds for the decay kinematics, but X μν is now constructed from the relative momentum p = (P −p 1 )/2 between K * and ψ in the B rest frame (we restored thep 1 for the momentum of the ψ in the decay kinematics), and orthogonalized with respect to the B momentum p 2 , A λ (s, t) = * μ (λ,p 1 ) −g μν + P μ P ν s X ν (q, P)g S (s) where X ρμ depends onp μ ⊥ =p μ − p μ 2 p 2 ·p/m 2 2 , and g ρμ ⊥ = g ρμ − p ρ 2 p μ 2 /m 2 2 . As mentioned before, crossing symmetry requires the helicity amplitude A λ (s, t) to be the same up to a phase for the decay and the scattering process. The expressions for the helicity amplitudes read where γ (s) = (s − m 2 1 + m 2 2 )/(2m 2 √ s) is the boost factor of K * in the B rest frame discussed in [14,19]. The matching to the general form in Eq. (12) is analogous. Although Eqs. (19) and (23) agree at threshold, the dynamical models differ in general, and the additional factors appearing in Eq. (23) are part of the model.
The issue with the crossing symmetry is particularly interesting, and we want to illustrate it further on a simpler example. We consider the decay B →Dππ. This reaction has been analyzed by LHCb using the CPM formalism [32]. Since none of the external particles have spin, the reaction is described by a single scalar amplitude. Because of crossing symmetry, the amplitude is the same scalar function of the Mandelstam variables in both scattering and decay kinematics. We consider the s-channel scattering kinematics B D → π + π − . We use the indices 1, 2, 3 and 4 for the D, B, π + and π − momenta. Therefore, P = p 1 + p 2 is the center of mass momentum, s = P 2 is the invariant mass of the ππ system, and p and q are the breakup momenta for the initial and final states. For simplicity, we restrict this discussion to the case of a spin-1 isobar, B D → ρ → π + π − . The CPM amplitude is given by In the center of mass frame of the s-channel, X μ ( p, P) = (0, 0, 0, p) and X ν (q, P) = (0, q sin θ s , 0, q cos θ s ) are purely spacelike vectors proportional to the breakup momenta. Therefore, the amplitude in Eq. (24) matches the expectations. For the decay process, the orbital tensor X μ ( p, P) is replaced by X μ (p, p 2 ) to be orthogonalized to the fourmomentum p 2 . As a result, a factor γ (s) = (s − m 2 1 + m 2 2 )/(2m 2 √ s) appears, and the breakup momentum from the B →Dρ orbital tensor is evaluated in the rest frame of B. The amplitudes for the decay process crossed to the scattering kinematics are A [B D→ππ] = pq cos θ s g P (s), The two amplitudes differ by a factor γ (s) √ s/m 2 = (s − m 2 1 + m 2 2 )/(2m 2 2 ). While this factor is analytic in s and does not spoil the counting of kinematical singularities discussed in the previous section, its appearance breaks crossing symmetry and this shows the drawback of the CPM formalism.
The issues arise from the construction of an amplitude as subsequent one-to-two decays. At first sight this appears as a natural choice. However, a well defined amplitude should have only asymptotic states on the external legs. This would exclude any decay into a resonance. One needs to take a step back to the definition of a resonance, i.e. a pole in the scattering amplitude. Therefore, the consistent procedure would be to write the amplitude in the scattering kinematics and then use crossing symmetry to analytically continue the amplitude into the decay region.

K π -mass distribution in different approaches
To explore the differences between the various approaches, we consider the example of two intermediate vectors in the π K channel: the K * (892), with mass and width M K * = 892 MeV, K * = 50 MeV, and the K * (1410), with M K * = 1414 MeV, K * = 232 MeV. The differential width is given by the expression, where ρ(s) = λ 1/2 12 λ 1/2 34 /s, and N j is a normalization constant. In Fig. 3 we show the results for five different scenarios. We consider the CPM formalisms discussed in Eqs. (19) and (23) (for the scattering and decay kinematics, respectively), setting g S (s) = 0 and g D (s) = T K * (s), with For the LS formalism, we choose the couplings in Eq. (16) to beĜ 1 0 (s) = 0,Ĝ 1 2 (s) = T K * (s). The LS amplitude in the decay kinematics differs from the one in the scattering kinematics only because of the breakup momentum of B → ψ K * , calculated in the B rest frame or in the K * rest frame, respectively. Finally, we draw a line for our proposal given by Eq. (B3), the only nonzero term in the sum is F 1 2 (s) = T K * (s).
The partial wave amplitudes for two-to-two scattering processes are proportional to p L 1 q L 2 , where p (q) are the initial (final) state break up momentum and the particles in L 1 (L 2 )-waves. This behavior comes from the expansion of the amplitude at threshold, and generates an unphysical growth at higher energies. This behavior is customarily modified by model-dependent form factors. 3 The most popular approach is based on a nonrelativistic model introduced by Blatt and Weisskopf [48,49]. For the plot on the right in Fig. 3 we multiply the amplitudes by the Blatt-Weisskopf barrier factors for the initial P-and final D-waves, respectively. The couplings are set as g S (s) =Ĝ 1 0 (s) = F 1 0 (s) = 0 and for the corresponding formalisms in Eqs. (18), (19) and (B3). The constant R is chosen to be 5 GeV −1 , which corresponds 1 fm, i.e. the scale of the strong interaction.
We see in Fig. 3 that the K π invariant mass squared distribution is distorted differently in all models. It is straightforward to track down where the differences come from. In the JPAC amplitude of Eqs. (B3) and (B4), the threshold factor in the F 1 2 (s) function in Eq. (B3) is set to λ 12 , in contrast to the CPM and LS formalisms where the factor p 2 = λ 12 /(4s) is used. This makes the differential width distribution different by the factor 1/s 2 . Another difference originates from the factor E 1 /m 1 for the A 0 (s) amplitude which was required by analyticity. We showed that it is not present in the LS approach, as one can see in Eq. (16). In the physical domain this factor behaves as 1/ √ s at the amplitude level, resulting in 1/s difference in the the differential width.

Summary and conclusions
We considered different approaches for constructing amplitudes for scattering and decay processes. Although the problem might be viewed as standard exercise, there seems to be confusion among amplitude analysis practitioners, as to which formalism best represents S-matrix constraints [18,19]. Specifically, we have compared the canonical helicity formalism [16,[22][23][24][25] and the covariant projection method [14,[18][19][20]. We used analyticity as a guiding principle to examine these approaches. Using as example the decay B → ψπ K , and the helicity formalism, we separated the kinematical factors from the dynamical functions. We then matched the helicity amplitudes with the most general covariant expression. In this process we identified kinematical constraints on the helicity amplitudes. We have shown that the naïve parameterization of the LS couplings fails to satisfy all the constraints required by analyticity. We found that, in contrast to LS parameterization, an extra factor (s + m 2 1 − m 2 2 )/(2m 1 √ s) naturally appears in the tensor formalism when written for the scattering kinematics. More interestingly, the customary recipes in the CPM approach explicitly violate crossing symmetry. In particular, we showed that the tensor approach discussed in [14,[18][19][20], when applied to the decay kinematics directly, introduces a peculiar energy dependence which has no clear physical motivation.
To address the issue of the relativistic corrections, we recall the relation between the helicity and the LS amplitudes. This relation is valid for any energy. The concept of the spin-orbit decomposition is fully relativistic. However, analyticity prevents the LS couplings to be parameterized as simple constants. We remark that our observations and conclusions are strictly valid only when asymptotic states are considered. We performed extensive studies for the four-legs process which describes two-to-two scattering, or the oneto-three decay, when the mother particle has an (infinitely) narrow width. The extension to other reactions requires dedicated studies.
Regarding the question raised in the title, we conclude that the helicity formalism is the right guide. Specifically, kinematical singularities have been classified for helicity partial waves. However, crossing to the decay channel is more natural for covariant amplitudes, while the LS amplitudes have the correct threshold behavior built in. Since a model for an amplitude is always required, examining the behavior of all of these amplitudes as we did here can maximize the consistency of the model with S-matrix principles.
where the Levi-Civita tensor is defined by 0123 = 1 and D(s, t) is the singularity free scalar amplitude.
As follows from the derivation in [16], the kinematical factor is K ±0 = √ s pq = λ 1/2 12 λ 1/2 34 /(4 √ s), i.e. the minimal L in the final and initial states matches j. Removing also the half-angle factor ξ 10 (z s ) = sin θ s , we get where the partial wave helicity amplitudes can be written as The decay ψ → μ + (q 2 )μ − (q 1 ) can be attached by contracting the tensor amplitude A μ from Eq. (B5) with the tensor given by the fermion vector current.