2. Cross-Sections
We consider multiple-final-state mesons produced by linearly polarized photons diffractively colliding off a proton target at rest. The outline of the reaction is shown in
Figure 1.
Let
represent the complete set of variables needed to describe the decay of the resonance. In the case of two final state mesons, only two angles will be needed. We use the (
) angles of one meson in the Gottfried–Jackson (GJ) frame (see reference [
3], Appendix A, for frame definitions). In the case of more than two mesons in the final state, at least two more angles for each extra meson will be required.
The cross-section for the reaction
, where
will be written as [
3]
where
is the Lorentz-invariant (transition or scattering) amplitude and
is the Lorentz-invariant phase-space element (LIPS). The spin’s incoming and outgoing degrees of freedom are included in the sum over spins. The LIPS includes the kinematical constraints and
M includes the spin and production/decay-internal (transition) degrees of freedom. We can write
, where
i runs over all the incoming and outgoing particles.
To measure cross-sections experimentally, we normally “bin” or divide data into small ranges of one variable such that the dependence of the cross-section on that variable is suppressed. For example, if this division is performed with the mass, energy, and t-Mandelstam, only the angular dependencies for two produced particles will remain (more are needed for more final state particles, i.e., perhaps isobar properties; see
Section 5.1.2). All the “external” (normalization) dependencies can be taken into an overall constant,
(which we will just drop out afterward from the formulas, as they will not affect the overall behavior). Therefore, in a data bin (
,
t, and
M), we define an intensity
M is a representation of the scattering operator or transition operator,
T, and it can be written as
and then
and, further,
In the case that we have spin information from the incoming photon beam, we can include the beam polarization in our amplitude by defining the matrix
, and the
photon spin density matrix operator,
(see reference [
4] for details), as
Suppose that we prepare the polarization of the incoming photons or measure their states of polarization. The average over photon polarization will be completely described by this spin density matrix. In the case of a beam of linearly polarized photons, any polarized state can be written as a linear combination of two pure polarization states. Therefore, the general structure of this
matrix (for example, in the helicity basis defined by
and
) will be
The
are now the target (
) and recoil (
) helicities. We will assume, as is tradition in the study of meson production [
3,
5] and discussed in reference [
6] (Section 16.6) that the transition operator can be factorized into two parts: the production (of X) and the decay operators (of X), such that
Furthermore, we can take a completely orthogonal set of states,
, such that
, and include them in the previous relation such that
The set of states,
, a full set of intermediate states, we will call the
partial waves. Each of these states can be described by a set of quantum numbers, for example,
l,
m, isobars (mass, width). The total angular momentum by
(
and the total spin
, where l = 0, 1, 2... (S, P, D...) and m
will define the “waves” of the expansion. In practice, this expansion is truncated (to a very few states). We will refer to these quantum numbers as (l,m,I), where the I include all other parameters needed for a more extended model, for example, the isobar model parameters (see
Section 5.1.2).
The production amplitudes describe the strong interaction production mechanism that we are not able to calculate (without a phenomenological model). In a mass-independent PWA fit, the production amplitudes, in a given bin, will be considered constant, independent of the decay properties (for example, final particle angles). They function as
weights on each partial decay amplitude, and will be extracted (fitted) from the data. We will rewrite
and
being the production amplitudes and
the decay amplitudes. Note that the
A’s and
T’s are both complex numbers and that
A depends of the resonance quantum numbers and angles, while
T depends on the resonance and beam (photon) quantum numbers. The photon spin density matrix depends on the partial polarization
P and the polarization angle
(see reference [
4] for definitions).
Therefore, in the helicity basis [
7],
being the helicities of the incoming photon, and
the helicities of the target (outgoing) nucleons:
For example, in a two-meson final state, we have 2(from 2(from 2(from unknown parameters () for each wave to be fitted to the data.
3. Reflectivity
The effect of parity is defined as the inversion of the spatial coordinates with respect to the origin of coordinates. Most reactions in High-Energy Physics (HEP) are unchanged under this operation (as only weak interactions violate parity).
In our case, assuming vector meson dominance [
8] for the photon and diffractive scattering from the nucleon, and since the strong interaction conserves parity, the parity operator commutes with the scattering matrix (or transition operator). Helicity states, however, are not eigenstates of the parity operator and therefore, they are not directly related to the parity exchanged in the reaction.
The parity operation is equivalent to a “mirror reflection” with respect to an arbitrary plane, followed by a
rotation with respect to an axis orthogonal to that plane. Let us call
the parity operator. Since the parity operation acting on rotations only changes the direction (sign), in the canonical representation (and in the rest frame of the particle), we have
where
are its eigenvalues. Let us consider a particle moving with momentum
in the
z direction. We can obtain this state by boosting (
L is a Lorentz transformation) the state at rest:
Applying the parity operator,
To get back from
to
, we need a rotation of modulo
around the
y axis
and we know that
Therefore, we finally have
Since any other direction can be constructed by rotation, and the parity operator commutes with rotations (in the x-z plane), we can express the former formula in the rest frame of the resonance, with
y perpendicular to the production plane (GJ/HEL frames), with the spin quantization in the z-axis given by
mIt is useful to define the reflection operator [
9]
which involves parity and a
angular rotation around the
y axis either in the GJ or HEL frames. It represents a mirror
reflection through the production plane (x,z). This operator commutes with the transition operator. The
y axis in the GJ/HEL frame is perpendicular to the production plane; therefore, the transition matrix is independent of
y, and only the
coordinates participate in the parity transformation. Reflection commutes with the Hamiltonian. The reflection operator acting on the resonance states produces
where
P are the parity eigenvalues
. We can build the following eigenstates of
(since the reflection changes signs on the z-projection quantum numbers,
m, we will create eigenstates that are a linear combination of both (m) signs’ states with adequate coefficients):
The sign between both terms in Equation (
26) is arbitrary. We use the sign definition in reference [
3] and define
and
It can be shown (see ref. [
3]) that the
’s are the real (for mesons) eigenvalues of the reflectivity operator. We define a
resonance reflectivity = as
In our previous notation,
Notice that since each state defined in the reflectivity basis includes a combination of
m and
, the projections of the spin on the quantization axis,
m, are replaced by
(a kind of “absolute value”). We can think of the reflectivity
“carrying” the sign of
m. When we sum over possible quantum numbers for each wave (
l), we have (2l + 1) terms in this sum; we have
for two reflectivities for each
plus one
for m = 0 [
3].
In pion beam experiments [
5] (spinless beam) or past photo production experiments (CLAS) [
10], where no information on the beam polarization was available, the spin density matrix is (or is considered) a constant (see [
3]), and therefore, it can be factored out from the intensity expression. The past CLAS formalism [
3] includes the helicity of the photon in the rank of the matrices (in the external sum of spins). Invoking parity conservation, we still reduced the number of degrees of freedom from eight to four, and the reflectivity was only defined for the resonance. Again, this was done for unpolarized photons or when no information about the photon polarization was available.
In the case we are considering, having information about the photon polarization, it will be proper to also define a reflectivity state for the photon. For a real photon
,
and
; therefore, we define a
photon reflectivity =
from
then (the reflectivity eigenvalues for a photon are
),
Equation (
14), in this new (two) reflectivity basis, is then
The photon spin density matrix in the photon reflectivity basis has the following form (see
Appendix A or reference [
4]):
We now write the expression for the intensity, where are defined positive in the reflectivity basis, and we include the “resonance” and “photon” reflectivities.
There are only two degrees of freedom associated with
target spins; we will call them
(spin-flop and no spin-flop).
We have organized the indices such that
k are the external or non-interfering indices. Expanding the sum over the photon reflectivities (using the photon spin density matrix), we have (just for clarity, we drop the
indexes in the next expression)
We define the
resonance spin density matrices as
We can write (see reference [
4])
where
are the Pauli matrices and
the photon polarization vector. Therefore,
which can be written as
with
being the polarized SDME (Spin Density Matrix Elements).
4. Naturality and Reflectivity
A state is said to have natural parity if
, while is said to have unnatural parity if
. We can recast this definition by introducing the
naturality of the particle,
, as
Naturality is (natural) for (i.e., ) and (unnatural) for (i.e., ). Determining (or constraining) the naturality of the production (exchange particle) will give us extra information on the produced resonances.
The reflectivities are defined for the resonance decay
and for the incoming photon. Reflection is a conserved quantum number since both rotation and parity are conserved. Therefore, at least at higher energies (see Appendix C of reference [
11]), the product of the initial photon reflectivity and the exchange particle reflectivity must equal the reflectivity of the resonance:
or
And since
, then
The photon spin density matrix, in the photon reflectivity basis, represents a mix of photon states, as seen in Equation (34). We can also see that the resonance spin density matrix (Equation (
37)) will not be diagonal in this formalism. Only for full polarization
might there be defined reflectivity configurations that contribute to the reaction. These are when
only (
) contribute and when
only (
) contribute. Using linearly polarized photons at those explicit configurations, we could then constrain the naturality of the exchange and particles produced. For example, in the case of pion exchange (or other Regge unnatural trajectory particles), the reflectivity of the resonance (
) will be opposite to that of the photon (
). In the case of
exchange (or other Regge natural trajectory particles), the reflectivity of the resonance and the photon will be the same. For unpolarized beams, the reflectivity is only defined for the resonance and the spin density matrix of the reaction becomes diagonal. For polarized beams, we can still use similar methods if we include the beam polarization in the rank of the sum (added to the external spin). The JPAC collaboration [
11] defined a reflectivity for the case of two pseudo-scalar final states, taking into account combined photon resonance parity conservation. In that case, there is only one reflectivity, and the spin density matrix becomes diagonal. In the JPAC definition,
or
and the reflectivity coincides directly (by construction) with the naturality of the resonance. It has been shown [
12] that the JPAC definition and the two reflectivity scenarios defined in this paper are equivalent for the case of two final state pseudo-scalars.
6. Phenomenological Models
After performing mass-independent fits in each bin of
M (or
M and
t) for a given
, we obtained the predicted mass distribution of
for each partial wave included in the fit. Nevertheless, merely identifying peaks in the mass spectrum falls short of substantiating the existence of a resonance. In the past, the mass dependence of those partial waves has been described by a coherent sum of Breit–Wigner amplitudes and, if needed, a phenomenological model of the background or other effects (i.e., Deck mechanisms) [
3,
5]. Such a procedure can produce a good fit to the data; however (especially using the isobar approximation), it violates fundamental principles such as probability conservation and causality. Therefore, in order to obtain more physically grounded amplitudes, models that fulfill the principles of unitarity and analyticity (which originate from probability conservation and causality) are to be used. Unitarity is especially important when we deal with resonances since it controls resonance widths and pole positions in the complex energy plane. One will first look for regions of enhancement (peaks or valleys) in the distributions and fit a theoretically based distribution to obtain the resonance properties (mass and width). However, interference and overlapping can greatly disturb the appearance of the spectrum. The properties (and positions) of the resonances should be obtained from the poles on the complex amplitudes of the S-Matrix expansion [
17]. These poles (and thresholds) had been studied using the Regge treatment of the S-matrix [
8]. Resonances are poles in the complex plane (Riemann surfaces) and only their projected real axis values can then be evaluated experimentally. In the case of multiple poles with the same quantum numbers and/or poles far from the real axis, the axis projections can deviate considerably from the BW distribution. The shape of these distributions is also influenced by the QCD dynamics. Effective field theories, i.e., Chiral Perturbation Theory, has been combined with the dispersion relations to obtain better parameterization of the mass distributions [
18]. Recent studies (i.e., reference [
19,
20]) have used those approaches to obtain mass and width values for several resonances. A comprehensive description of recent efforts by the Joint Physics Analysis Center (JPAC) in this direction can be found in reference [
21].