# What is the right formalism to search for resonances?

###### Abstract

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 and 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.

###### pacs:

11.55.Bq, 11.80.Cr, 11.80.Et^{†}

^{†}preprint: JLAB-THY-17-2606

Joint Physics Analysis Center

## I 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 study the hadron spectrum from first-principles lattice QCD calculations Wilson *et al.* (2015); Briceño *et al.* (2017a, b, c); Hu *et al.* (2016).

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 Lebed *et al.* (2017); Esposito *et al.* (2016); Olsen *et al.* (2017), 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 Anisovich *et al.* (2012); Workman *et al.* (2012),
or meson spectroscopy from pion or photon beam fragmentation Abbon *et al.* (2015); Shepherd (2009); Glazier (2015).

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 Chung (1993); Anisovich *et al.* (2012); Adolph *et al.* (2017).
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 Collins (2009).
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 Cohen-Tannoudji *et al.* (1968a), 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) Chung (1993); Chung and Friedrich (2008); Filippini *et al.* (1995); Anisovich and Sarantsev (2006),
that starts from a LS partial wave model (or equivalently the Cartesian, aka Zemach amplitudes Zemach (1965)) but writes them
in a covariant fashion.
The method has a drawback, which is related to the behavior under crossing (see Section III.1). The alternative, which we refer to as
the canonical approach Jacob and Wick (1959); Chung (1971); Collins (2009); Kutsckhe , is to use the well known relation between the helicity amplitudes and the helicity partial waves Collins (2009) 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 , with .
The analyses by Belle and LHCb show nontrivial structures appearing in the Mizuk *et al.* (2009); Chilikin *et al.* (2013); Aaij *et al.* (2015a, b), and in the channel Chilikin *et al.* (2014).
These are of particular interest, because a resonance in these channels would require an exotic interpretation Lebed *et al.* (2017); Esposito *et al.* (2016); Olsen *et al.* (2017).
The rest of the paper is organized as follows. In Sec. II we discuss the canonical approach on the example of the 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 Matevosyan *et al.* ; Neubert .
In Sec. III we examine the CPM approach and compare this model with the findings from Sec. II. We mention the crossing symmetry properties of CPM using, as an example, , which was recently analyzed by LHCb within this formalism Aaij *et al.* (2015c).
Summary and conclusions are given in Section IV.

## Ii Analyticity constraints for

In Fig. 1 we show a diagram representing the kinematics of the decay . The spinless particles , , 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 , to label the momenta of , , and respectively. The momentum of the will be denoted by , for a reason which we will explain below. The helicity amplitude for the decay process is denoted by , being the helicity of , i.e. . The amplitude depends on the standard Mandelstam variables , , and with .

The meson decays weakly, so 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 (, , ). 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 setting Trueman and Wick (1964). 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 -channel. The -channel resonances correspond to the ’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 -channel. Once we have constructed the -channel amplitudes, the -channel ones can be treated similarly (cf. Appendix B).

In the center of mass of the -channel scattering process, the momentum defines the -axis, the momenta and lie in the -plane. We call () to the magnitude of the incoming (outgoing) three momentum. The scattering angle is a polar angle of the pion (see Fig. 2). The quantities depend on the Mandelstam invariants through

(1) |

with . The function is a polynomial in . To incorporate resonances in the system with a certain spin , we expand the amplitude in partial waves,

(2) |

where are the helicity partial wave amplitudes in the -channel. In Eq. (2) the entire dependence enters though the functions. The functions have singularities in which lead to kinematical singularities in of the helicity amplitudes . The dynamical singularities in , 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 - or -channel resonances (singularities) are accounted for explicitly through - or -channel partial waves, and to avoid double counting each series is truncated at a finite number of terms. This defines the so-called isobar model in which

(3) |

with,

(4) |

where is finite. The expressions for the and isobars are similar to Eq. (4).
Note, that due to the superscript the amplitudes are not identical to the helicity partial waves, of Eq. (2).
This is because the other two terms on the right hand side of Eq. (3) also contribute to the -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 Khuri and Treiman (1960); Niecknig *et al.* (2012); Niecknig and Kubis (2015); Danilkin *et al.* (2015); Guo *et al.* (2015); Pilloni *et al.* (2017); Albaladejo and Moussallam (2017).
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 Hara (1964); Wang (1966); Jackson and Hite (1968); Cohen-Tannoudji *et al.* (1968b); Martin and Spearman (1970); Collins (2009).
We recall that , where is the so-called half angle factor that contains all the kinematical singularities in . The reduced rotational function is a polynomial in and of order divided by the factor .
The helicity partial waves have singularities in . These have both dynamical and kinematical origin. The former arise, for example, from -channel resonances.
The kinematical singularities, just like the -dependent kinematical singularities, arise because of external particle spin.
We explicitly isolate the kinematic factors in ,
and denote the kinematical singularity-free helicity partial wave amplitudes by .^{1}^{1}1For fermion-boson scattering, the can still have a branch point at , as discussed in Cohen-Tannoudji *et al.* (1968b).
First, the term is factorized out from the helicity amplitude . This factor is there
to cancel the threshold and pseudothreshold singularities in that appear in .
Second, we follow Collins (2009) and introduce the additional kinematic factor
(‘’ is short for ). These factors are required to account for a mismatch between the and dependence in the angular momentum barrier factors in presence of particles with spin.
Finally, the kinematical singularity-free helicity partial wave amplitudes are defined by

(5a) | ||||

(5b) | ||||

(5c) | ||||

with and given by | ||||

(5d) | ||||

(5e) |

Specifically, it is expected that at threshold, where and are the lowest possible orbital angular momenta in the given helicity and parity combination. This explains why requires a special treatment Martin and Spearman (1970), since for we have , but for the lowest is . In addition, the -factors have powers of as required to ensure factorization of the vertices of Regge poles. Similarly, as explained before, is dynamical and thus the kinematical singularity-free amplitudes are not expected to contain singularities in , and as will be seen below, the dependence of the -factor takes care of that. The 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.

We now seek a representation of in terms of the scalar functions, as discussed in Sec. I. For the PV amplitude this is given by

(6) |

Although the second term in brackets may look like an extra pole, it cancels when multiplied by . 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 -channel center of mass frame the polarization vectors are given by for the transverse polarizations and for the longitudinal polarization. The energies are calculated from the momenta and are fully determined by . The functions and are the kinematical singularity free scalar amplitudes discussed in the Sec. I.

We can match Eqs. (2) and (6), and express the scalar functions as a sum over kinematical singularity free helicity partial waves. The ratio , computed using Eq. (6), is compared to the same ratios computed using the helicity partial waves from Eq. (2). This yields

(7) | ||||

(8) |

from and , respectively, which can be combined to obtain

(9) |

Neither nor can have kinematical singularities in or . In Eqs. (7)-(9), is regular in , and the singularities at (pseudo)thresholds are canceled by the factor .
The latter factor contains a high-order pole at . Such pole
is a feature of the dynamical model, and specifically arises because at the little group is not 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 Toller (1965); Collins (2009).^{2}^{2}2For example, in Regge theory these poles are canceled by the daughter Regge trajectories Freedman and Wang (1967); Collins (2009). This construction hence does not constrain the poles at .

For the same reason the sum in Eq. (9) has no kinematical singularities in and , however the factor in front of the sum generates two poles at , unless the expression in brackets vanishes at those points. This means that the with different cannot be independent functions at the (pseudo)threshold. Explicitly, in the limit at fixed one has and using Jackson and Hite (1968),

(10) |

one finds that the expression within the brackets in Eq. (9) behaves as

(11) |

This combination has to vanish to cancel the , thus one finds (for )

(12a) | ||||

(12b) |

where , , , and are regular functions at , and . 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 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 Sec. III. 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.

Upon restoration of the kinematic factors, the original helicity partial wave amplitudes read ()

(13a) | ||||

(13b) |

and , where is regular at (pseudo)threshold. A particular choice of the functions , , and constitutes a given hadronic model. A specific example is given in Appendix B.

### ii.1 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 , and corresponding helicities , , the relation between a two-particle state in the helicity and LS basis is

(14) |

where is the projection of the total angular momentum . For the amplitude, it implies the following relation between the LS amplitudes and the helicity amplitudes,

(15) |

The amplitudes with and differ by parity. Equation (15) can be inverted to relate the helicity partial wave amplitudes with the LS amplitudes ,

(16) |

In Eq. (16) we denoted the LS partial wave amplitudes with the threshold factors explicitly factored out by , i.e. . We now compare the general expression for the helicity partial waves with the spin-orbit LS partial waves. We find that Eq. (16) matches the general form in Eq. (12) when

(17a) | ||||

(17b) | ||||

(17c) | ||||

(17d) |

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 Eqs. (17) are free from kinematical singularities. For example, if one takes the functions and to be proportional to Breit-Wigner functions with constant couplings, the amplitudes and would end up having a pole at , and/or a branch point at unexpected for boson-boson scattering. On the other hand, as discussed in Sec. II, the pole at is part of the dynamical model. It is clear that using Breit-Wigner parameterizations, or any other model for helicity amplitudes, i.e. the left-hand sides of Eq. (17), instead of the LS amplitudes helps prevent unwanted singularities.

## Iii Comparison with the Covariant Projection Method

We consider now the CPM approach of Chung (1993); Chung and Friedrich (2008); Filippini *et al.* (1995); Anisovich and Sarantsev (2006). As said, the method is based on the construction of explicitly covariant expressions. To describe the decay , we first consider the polarization tensor of each particle with index and spin , . Using the decay momentum and the total momentum , we build a tensor to represent the orbital angular momentum of the system.
In order to find total angular momentum tensor, we first combine the polarizations of and into a “total spin” tensor
(orthogonal to the momentum ).
Then, we combine the tensor with the orbital tensor
and finally contract the result with the polarization of , thus mimicking the LS construction.
The tensors and 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, .
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 . 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 with . We start with the tensor amplitude for the scattering process . The orbital angular momentum of the decay in -wave is given by . The tensor is constructed from a four-vector of the relative momentum and the total momentum of the system . For the PV amplitudes, the initial process is described by two waves. The corresponding orbital tensors are the unit rank- tensor for the -wave and rank-2 tensor , with , for the -wave. Hence

(18) |

where is the momentum. The final -wave orbital tensor is . The -wave orbital tensor , with , and . Explicitly,

(19) |

and matching with Eq. (12) gives

(20a) | ||||

(20b) | ||||

(20c) |

The threshold conditions are satisfied, and the functions and are regular at the thresholds. Finally, we show the relation between the CPM and the LS amplitudes. The comparison with Eq. (16) leads to

(21a) | ||||

(21b) |

Although the and of the CPM formalism, see Eq. (18), are typically interpreted as the and partial wave amplitudes, we see that this is the case only at (pseudo)threshold , where the factor vanishes. In Sec. III.2 we discuss a specific example to show the numerical difference between the various approaches.

### iii.1 Crossing symmetry; and the decay

An issue with the CPM formalism is the explicit violation of crossing symmetry.
The recursive procedure explained in Chung (1993); Chung and Friedrich (2008); Filippini *et al.* (1995); Anisovich and Sarantsev (2006) 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 . The tensor describes the -wave decay as before.
- and -waves are still possible for the decay .
The same symbolic expression in Eq. (18) holds for the decay kinematics,
but is now constructed from the relative momentum between and in the rest frame (we restored the for the momentum of the in the decay kinematics), and orthogonalized with respect to the momentum ,

(22) |

where depends on , and . As mentioned before, crossing symmetry requires the helicity amplitude to be the same up to a phase for the decay and the scattering process. The expressions for the helicity amplitudes read

(23a) | ||||

(23b) |

where is the boost factor of in the rest frame discussed in Filippini *et al.* (1995); Chung (1993).
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 . This reaction has been analyzed by LHCb using the CPM formalism Aaij *et al.* (2015c). 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 -channel scattering kinematics . We use the indices , , and for the , , and momenta. Therefore, is the center of mass momentum, is the invariant mass of the system, and and are the breakup momenta for the initial and final states. For simplicity, we restrict this discussion to the case of a spin- isobar, .
The CPM amplitude is given by

(24) |

In the center of mass frame of the -channel, and 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 is replaced by to be orthogonalized to the four-momentum . As a result, a factor appears, and the breakup momentum from the orbital tensor is evaluated in the rest frame of . The amplitudes for the decay process crossed to the scattering kinematics are

(25) |

The two amplitudes differ by a factor . While this factor is analytic in 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.

### iii.2 -mass distribution in different approaches

To explore the differences between the various approaches, we consider the example of two intermediate vectors in the channel: the , with mass and width MeV, MeV, and the , with MeV, MeV. The differential width is given by the expression,

(26) |

where , and is a normalization constant.

In Fig. 3 we show the results for five different scenarios. We consider the CPM formalisms discussed in Eq. (19) and Eq. (23) (for the scattering and decay kinematics, respectively), setting and , with

(27) |

For the LS formalism, we choose the couplings in Eq. (16) to be , . The LS amplitude in the decay kinematics differs from the one in the scattering kinematics only because of the breakup momentum of , calculated in the rest frame or in the rest frame, respectively. Finally, we draw a line for our proposal given by Eq. (35), the only nonzero term in the sum is .

The partial wave amplitudes for two-to-two scattering processes are proportional to , where () are the initial (final) state break up momentum and the particles in ()-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}^{3}3i.e. having left-hand singularities only.
The most popular approach is based on a nonrelativistic model introduced by Blatt and Weisskopf Blatt and Weisskopf (1952); Von Hippel and Quigg (1972).
For the plot on the right in Fig. 3 we multiply the amplitudes by the Blatt-Weisskopf barrier factors

(28a) | ||||

(28b) |

for the initial - and final -waves, respectively. The couplings are set as and for the corresponding formalisms in Eqs. (18), (19) and (35). The constant is chosen to be GeV, which corresponds fm, i.e. the scale of the strong interaction.

We see in Fig. 3 that the 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. (35) and (36), the threshold factor in the function in Eq. (35) is set to , in contrast to the CPM and LS formalisms where the factor is used. This makes the differential width distribution different by the factor . Another difference originates from the factor for the 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 at the amplitude level, resulting in difference in the the differential width.

## Iv 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 -matrix constraints Chung and Friedrich (2008); Filippini *et al.* (1995). Specifically, we have compared the
canonical helicity formalism Jacob and Wick (1959); Chung (1971); Collins (2009); Kutsckhe and
the covariant projection method Chung (1993); Chung and Friedrich (2008); Filippini *et al.* (1995); Anisovich and Sarantsev (2006).
We used analyticity as a guiding principle to examine these approaches.
Using as example the decay , 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 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 Chung (1993); Chung and Friedrich (2008); Filippini *et al.* (1995); Anisovich and Sarantsev (2006),
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 one-to-three decay, when the mother particle has an (infinitely) narrow width. The extension to other reactions requires dedicated studies.

###### Acknowledgements.

We are grateful for the inspiring atmosphere at the PWA9/ATHOS4 workshop where the idea for this work was born. AP thanks Greig Cowan and Jonas Rademacker for useful discussions about the application of the different formalisms in LHCb. MM and AP would like to thank Suh-Urk Chung and Andrey Sarantsev for several useful discussion about the CPM approach, Dmitry Ryabchikov for sharing his experience with partial wave analysis using the Zemach and LS formalisms, and Anton Ivashin for several useful discussions. This work was supported by BMBF, the U.S. Department of Energy under grants No. DE-AC05-06OR23177 and No. DE-FG02-87ER40365, PAPIIT-DGAPA (UNAM, Mexico) grant No. IA101717, CONACYT (Mexico) grant No. 251817, Research Foundation – Flanders (FWO), U.S. National Science Foundation under award numbers PHY-1507572, PHY-1415459 and PHY-1205019, and Ministerio de Economía y Competitividad (Spain) through grant No. FPA2016-77313-P.## Appendix A Parity-conserving amplitudes in the channel

We consider the most general scattering amplitude for with parity conservation enforced. The only tensor structure allowed is