Spin density matrix of the ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} in the reaction p¯p→ωπ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{p}p}\,\rightarrow \,\omega \pi ^0$$\end{document}

The spin density matrix of the ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} has been determined for the reaction p¯p→ωπ0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{p}p}\,\rightarrow \,\omega \pi ^0$$\end{document} with unpolarized in-flight data measured by the Crystal Barrel LEAR experiment at CERN. The two main decay modes of the ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} into π0γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi ^0 \gamma $$\end{document} and π+π-π0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi ^+ \pi ^- \pi ^0$$\end{document} have been separately analyzed for various p¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{p}}$$\end{document}momenta between 600 and 1940 MeV/c. The results obtained with the usual method by extracting the matrix elements via the ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} decay angular distributions and with the more sophisticated method via a full partial wave analysis are in good agreement. A strong spin alignment of the ω\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\omega $$\end{document} is clearly visible in this energy regime and all individual spin density matrix elements exhibit an oscillatory dependence on the production angle. In addition, the largest contributing orbital angular momentum of the p¯p\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{p}p~}$$\end{document}system has been identified for the different beam momenta. It increases from Lp¯pmax\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^{max}_{{\bar{p}p~}}$$\end{document}=\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$=$$\end{document} 2 at 600 MeV/c to Lp¯pmax\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^{max}_{{\bar{p}p~}}$$\end{document}=\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$=$$\end{document} 5 at 1940 MeV/c.


Introduction
The spin density matrix of particles originating fromp p annihilations provides important information about the underlying production process. The knowledge of this property is quite scarce in the low energy regime forp p in-flight reactions and is, however, very fundamental for high quality and high statistics future experiments like PANDA [1]. One This work comprises part of the thesis of J. Pychy. a e-mail: bertram@ep1.rub.de major physics topic of PANDA is the spectroscopy of exotic and non-exotic states in the charmonium and open charm mass regions inp p production or formation processes. For the identification of such resonances it is very helpful to know which initialp p states contribute and in particular how the corresponding production mechanism can be described in detail. The information about the contributing orbital angular momenta of the initialp p system and about the spin alignment of vector mesons produced in such processes is therefore an excellent key to gain a deeper insight into the production mechanisms. Therefore the investigation of the reaction p p → ωπ 0 with a relatively simple final state and without complex decay trees via intermediate resonances provides an excellent access to these questions. The ωπ 0 state couples only to isospin I = 1 and the C-parity C = −1 of thē p p system.
The data presented here have been measured with the Crystal Barrel experiment at LEAR in the years 1995 and 1996. A partial wave analysis has been performed with the PAWIAN software (Partial Wave Interactive Analysis Software) [2] by making use of the helicity formalism and considering the complete reaction chain. Various beam momenta have been studied between 600 and 1940 MeV/c and for two different ω decay modes, ω → π 0 γ and ω → π + π − π 0 , respectively. For the neutral decay mode the polarization of the radiative photon has not been measured and thus it is needed to average over this property.
Similar studies of this reaction for the charged decay mode of the ω have already been published in [3]. The results presented in the following rely on a more accurate data selection and a refined analysis. First preliminary results for the charged decay mode have already been presented in [2].

Crystal barrel experiment
The Crystal Barrel detector, which has been described in detail elsewhere [4], has been designed with a cylindrical geometry along the beam axis. Thep p annihilation took place in a liquid hydrogen target cell with a length of 4.4 cm and a diameter of 1.6 cm located in the center of the detector. This target was surrounded by a silicon vertex detector. This inner part was surrounded by a jet drift chamber which covered 90 and 64 % of the full solid angle for the inner and outer layer, respectively. These devices together with a solenoid magnet providing a homogeneous 1.5 T magnetic field parallel to the incident beam guaranteed a good vertex reconstruction, tracking and identification for charged particles. For accurate measurements of photons the detector was equipped with a barrel of 1380 CsI(Tl) crystals covering the full azimuthal range of 360 • and polar angles from 12 • to 168 • . With this electromagnetic calorimeter, assembled between the jet drift chamber and the solenoid magnet, an energy resolution of σ E /E ≈ 2.5 % and an angular resolution of 1.2 • in θ and φ each have been achieved.

Data selection and measured angular distributions
The data for this analysis have been taken over various beam times in the years 1995 and 1996 using an unpolarizedpbeam and an unpolarized liquid hydrogen target. In most cases the data samples have been recorded by utilizing a 0prong trigger for the neutral and a 2-prong trigger for the charged decay mode. In addition, a mixed trigger has been used where events with exactly 0 and 2 detected charged particles have been accumulated.
The offline reconstruction and event selection have been performed similarly to thep p annihilation at rest data [5]. In addition neural networks have been applied for the recognition of misleadingly reconstructed photons induced from electromagnetic [6] and hadronic [7] split-offs in the calorimeter. Only exclusive events are considered where all final state particles have been detected. In order to simply reduce the data samples to a more manageable size, preselection cuts have been carried out as follows: exact number of charged particles and photons in the final state and conservation of the total energy ( E tot = |E tot p p − E tot rec | < 500 MeV) and momentum ( p tot = |p tot p p − p tot rec | < 500 MeV/c) for the desired reaction. In addition exactly one π + π − pair must be reconstructed for the charged decay mode originating from a common vertex which is required to be within the target cell. After that kinematic fits with the hypothesesp p → π + π − 4γ, π + π − 2π 0 , ω2γ for the charged decay mode andp p → 5γ, π 0 π 0 γ, ω2γ for the neutral decay mode have been performed. Each individual fit requires the conservation of the momentum and energy of the events (4 constraints) and additional constraints on the π 0 -mass. Due to the fact that even with these fits the width of the reconstructed invariant mass of the ω is still dominated by the detector resolution further improvements of the quality of the data has been achieved by constraining the narrow mass of this vector meson (7-constraint fit:p p → ωπ 0 ). It is required that the fit converges with a confidence level (CL) greater than 10 % for each hypothesis. For all beam momenta the distribution of the confidence level is nearly flat and the distributions of the individual pulls are found to be Gaussian centered at 0 with a width of about σ ≈ 1. This is an indication for a good data quality and for a proper adjusted error matrix. As an example Fig. 1 shows these distributions for the neutral decay mode at 900 MeV/c.

Signal-background separation
The background contamination is caused by a variety of different sources. One scenario is that channels decaying to slightly different combinations of final state particles contribute where one particle remains undetected or energy deposits in the electromagnetic calorimeter originating from split-off effects are misinterpreted as an additional photon. Another possibility for the fulfillment of all selection criteria  Fig. 2 Dalitz plots for the selected π 0 π 0 γ events at thep momentum of 900 MeV/c (a) and 1940 MeV/c (b). The ω signal is visible as strong narrow bands parallel to the horizontal and vertical axis at ≈ 6·10 5 MeV 2 /c 4 . The remaining bands mainly originate from the f 2 (1270) π 0 background channel. The cross regions between the ω and the background bands are marked by white ellipses. The comparison between the two plots clearly demonstrates that the positions of the crossing regions strongly depend on the incident beam momentum is that even channels containing the same final state particles can contribute as background due to misleadingly combined decay products.
For the neutral channel the Dalitz plot of the selected π 0 π 0 γ events sheds light on the most crucial background source (Fig. 2). Besides the clear ω signal, structures from background events are visible whose major origin has been identified as the channelp p → f 2 (1270) π 0 → (π 0 π 0 ) π 0 → 6γ where one photon remains undetected. In this case the most problematic events are those which appear in the crossing regions of the signal and background band. Due to the fact that in this region the events are located in the same phase space volume it is impossible to reject the background by just applying the selection criteria as described above. Moreover these inhomogeneities of the background events along the ω band whose distribution is directly correlated to the one of the ω decay angle would result in huge systematic uncertainties for the determination of the spin density matrix. Since the positions of the crossing regions vary with the incident beam momentum, this situation becomes even more problematic.
In order to separate these non-interfering background sources from the signal events, an elaborated technique has been used where a signal weight factor Q has been assigned to each event. The strategy has been described in detail in [8] and was successfully applied on CLAS data for the reaction γ p → p ω [9,10]. Usual separation methods like the side-band subtraction method are based on the requirement of a binned data set. This exhibits disadvantages due to the complexity in a high dimensional phase space. Instead, the advantages of the technique used here is that it is an event based method and that detailed information about the specific background sources is not needed.
The method takes advantage of the fact that all noninterfering background events cannot reproduce the narrow  Fig. 4 Histogram a shows the normalized transition rate λ/λ max of the ω decay. Histogram b represents the invariant π + π − π 0 mass of all selected π + π − 2π 0 events at the beam momentum of 900 MeV/c. The excellent background separation power can be seen by the shaded areas representing the fraction of the signal events and the dotted red lines illustrating the background content resonance shape of the ω meson in the corresponding invariant mass spectrum. Therefore not the fitted ωπ 0 events but rather all selected and fitted π 0 π 0 γ events for the neutral and π + π − π 0 π 0 events for the charged decay mode appearing within a certain window around the relevant ω-mass shape (see Figs. 3b, 4b) are considered for the determination of the Q-value. The procedure starts with the assignment of the nearest neighbors for each event by defining a metric with the relevant kinematic observables. For the neutral channel the metric has been defined via three observables: the polar angle of the ω production in thep p rest frame and the azimuth and polar angle of the ω decay in its helicity system, in which the y-axis is defined to be parallel to the normal vector of the production plane. A subset of 200 neighbors for each event has been chosen which ensures that the associated events cover only a small region of the phase space. A Q value for each event is then obtained by the determination of the signal to background ratio in the invariant mass spectrum of the corresponding data subset. For this an unbinned fit has been performed with a convolution of a Gaussian and a nonrelativistic Breit-Wigner function for the description of the ω signal and a linear approximation for the background content. This approximation can be justified by the assumption that the background events are homogeneously distributed within the small region of the phase space. One example of this fit procedure is illustrated in Fig. 3a. The invariant π 0 γ spectrum ( Fig. 3b) shows the excellent result for the global signal-background separation obtained for the beam momentum at 900 MeV/c. With the outcome of this approach the final data set for the input of the partial wave analysis has been selected with the fitted ωπ 0 events, each weighted with the corresponding Q-factor.
The same event weight method has been performed for the charged decay mode. Here, the non-interfering background events exhibit as well different shapes in the invariant π + π − π 0 mass distribution in comparison to the ω signal. Potential interfering background sources are channels decaying into the same final state particles and can be estimated from thep p annihilation into the four charged pion final state [11], i.e.
Due to kinematic reasons these events do not overlap with the ωπ 0 events in the phase-space volume and thus do not contribute to the background. Also the small fraction of combinatorial background of these events do not interfere with the ωπ 0 channel and can therefore be eliminated by the Qweight method. For the charged decay mode the metric has been defined with four independent observables: the polar angle of the ω production in thep p rest frame, the azimuth and polar angle of the normal of the ω decay plane in its helicity system and the transition rate λ of the ω decay, which is characterized by the cross product of two pion momenta in the ω helicity frame [9,12,13]: where T π represents the kinetic energy of the individual pions. Figure 4 shows very impressively the obtained background separation power. Especially the shape of the normalized transition rate λ demonstrates the proper distinction between signal and background events. While the signal events follow the expected λ-shape for the ω decay with a linear increase and an intersection at the origin of the Table 1 Used data samples and number of selected events for the chan- axis (0 at λ = 0) the background results in an almost flat distribution.
3.2 Overview of the selected data samples Tables 1 and 2 summarize the numbers of the selected ωπ 0 events without and with the obtained weight factor Q for both decay modes. The number of ωπ 0 signal events is found to be between 1698 at 1525 MeV/c and 12,823 at 900 MeV/c for the charged decay mode and between 1113 at 600 MeV/c and 53,788 at 900 MeV/c for the neutral decay mode, respectively. The large variations of the ratio between the event numbers of the two decay modes for the different beam momenta are mainly caused by the use of different trigger configurations during the data taking. The final data sets consist of sufficient numbers of events for achieving significant results for the partial wave analysis and in particular for the determination of the spin density matrix of the ω. The background contamination estimated by the weight factor (1-Q) depends slightly on the beam momentum and on the decay pattern and varies between 9.2 and 14.6 % for the charged and 13.7 and 21.4 % for the neutral decay mode. Acceptance corrected angular distributions for the channel p p → ωπ 0 → (π 0 γ )π 0 as a function of the production angle (first column) and of the decay angle in cos(θ ω γ ) (second column) and φ ω γ (third column). The production angle is defined in thep p rest frame by the direction of the ω related to the beam axis. The decay angles arespecified by the helicity system of the ω meson. The production angular distribution is given integrated over all ω-decay angles, the decay angular distributions are given integrated over all production angles. While the data are marked with red error bars, the fit results (Sect. 4) are plotted with black lines. Each row represents one specific beam momentum Acceptance corrected angular distributions for the channel p p → ωπ 0 → (π + π − π 0 )π 0 as a function of the production angle cos(θp p ω ) (first column) and of the decay angle in cos(θ ω n ) (second column) and φ ω n (third column). The production angle is defined in thep p rest frame by the direction of the ω related to the beam axis. The decay angles θ ω n and φ ω n are specified by the normal of the decay plane of the ω meson in its helicity system. The production angular distribution is integrated over all ω-decay angles, the decay angular distributions are integrated over all production angles. The fourth column represents the normalized transition rate λ of the ω decay. A detailed description of this property can be found in Sect. 3.1. While the data are marked with red error bars, the fit results (Sect. 4) are plotted with black lines. Each row represents one specific beam momentum background rejection criteria for the neutral and charged decay mode, respectively. The distributions of the ω production angle are integrated over the ω-decay distributions and are characterized by fluctuations of the intensity with a higher number of extrema for increased beam momenta. This is an indication that more waves contribute with the rise of the center of mass energy. The huge error bars and the absence of entries around cos(θp p ω ) = 1 are caused by the acceptance leakage of the detector in the very forward and backward region. These inefficiencies are more distinctive for the charged decay mode due to the limited angular coverage of the tracking devices and become even more apparent with increasing beam momentum. The distributions of the ω-decay angles are integrated over all production angles and exhibit typical shapes for this particle (see Sect. 5).

Measured angular and λ-distributions
For all beam momenta the normalized λ-distributions for the ω-decay to π + π − π 0 (Fig. 6) are in excellent agreement with the expected shape. This illustrates again the high purity of the ωπ 0 data for the individual beam momenta.

Amplitudes
p p in-flight reactions where mesons and photons are exclusively involved are dominated by the s-channel process. Therefore the partial wave analyses for those reactions have been started usually with the J PC system initiated from thē p p annihilation. One difficulty of this method is that additional Clebsch-Gordan coefficients for the coupling of thē p p system with the J PC intermediate state are not considered correctly. In order to avoid such error-prone procedure, the analysis performed on the data here is based on the description of the complete reaction chain starting from thep p coupling up to the final states. This new method is summarized in detail in [14] and can also be applied to other p p reactions in flight.
The starting point is the description of the differential cross section of the whole reaction chain where the transition amplitude depending on the helicities of the involved particles is divided into the ω-production and the ω-decay amplitude. For the neutral channel this cross section is expressed by where dτ represents the infinitesimal volume element of the phase-space, w the transition probability, λ the helicities of all involved particles, T the production and A the decay amplitude in the helicity frame. The two neutral pions are distinguished by the notation π 0 r for the recoil particle and π 0 d for the ω decay particle. Due to the fact that a mass constraint for the ω has been used for the kinematic fit the dynamics for this meson (e.g. a Breit-Wigner distribution) has not been taken into account. It is noteworthy to mention that the components of the transition amplitude are added coherently over the helicities of the intermediate ω-resonance and incoherently over the helicities of all initial and final state particles. Equation 4 is expanded into states with definite J PC -values defining the partial wave helicity amplitudes T J PC λp,λ p ,0,λ ω and A 1 − λ ω λ γ . These partial wave amplitudes are further expanded in states with definite J PC , L, S-values where L , S are the respective orbital angular momenta and total spins of thep p , ωπ 0 and π 0 γ -system (Lp p , Sp p , L ωπ 0 , S ωπ 0 (= 1), L π 0 γ (= 1), S π 0 γ (= 1)), defining the amplitudes T Jp p Lp p ,Sp p ,L ωπ 0 and A 1 L π 0 γ ,S π 0 γ = A 1 11 . Here, the quantum number J represents the total angular momentum, L the orbital angular momentum and S the total spin of the related system composed of two particles. The underlying formalism for theses expansions can be found in detail elsewhere [15]. With the requirement that the parity, charge conjugation and total angular momentum are conserved for strong and electromagnetic interactions the differential cross section can be described by incoherent sums over thep p singlet and triplet states and over the helicity of the radiative photon of the final state system [14]. In terms of L S-amplitudes Eq. 4 reads [14]: with Mp p = λp −λ p . The summation over λp and λ p can be arranged in such a way, that one incoherent term for singlet states (Sp p = 0, Mp p = 0) and three incoherent terms for triplet states (Sp p = 1, Mp p = 0, ±1) appear. The direction of thep beam is chosen as the quantization axis which results in the restriction of the z-component of Jp p to Mp p = 0, ±1. The ωπ 0 system is fully characterized by L ωπ 0 , S ωπ 0 = 1, the helicity λ ω and the production angle θp p ω of the ω in thep p rest frame. Due to the fact that thē p p system is unpolarized the angle φp p ω is not defined. The ω decay system is characterized by the angular momentum L π 0 d γ = 1, the total spin S π 0 d γ = 1, the helicity λ γ and the decay angles θ ω γ and φ ω γ of the γ in the helicity frame of the ω meson. The d Jp p Mp p ,λ ω denotes the Wigner-d function for the decay of thep p system, D * J ω =1 λ ω ,λ γ the complex conjugate of the Wigner-D function for the ω decay and L , 0, S, λ 1 − λ 2 |J, λ 1 − λ 2 and j 1 , λ 1 , j 2 , −λ 2 |S, λ 1 − λ 2 the Clebsch-Gordan coefficients for the L S-and j j-coupling respectively. As A 1 11 for a givenp p -energy is a fixed complex number, the product T Jp p Lp p ,Sp p ,L ωπ 0 · A 1 11 is handled as one complex parameter α J PC Lp p L ωπ 0 . For the reactionp p → ωπ 0 → (π + π − π 0 ) π 0 Eq. 5 has to be modified [15]. The incoherent sum over λ γ vanishes and the ω decay amplitude A 1 − λ ω λ γ has to be replaced by where ω n , ω n , γ ω n are the Euler angles of the normal of the 3π -decay plane (n) in the ω-helicity system with μ = J ω ·n . In general μ takes the values ±1, 0, but in the ω → π + π − π 0 case, only μ = 0 is allowed. A 1 μ (E π + , E π − ) describes the amplitude in the Dalitz plot, which is proportional to |P π + × P π − | [12].
By making use of the conservation principles and the selection rules one can easily extract the specific combinations of the relevant quantum numbers allowed for the reactionp p → ωπ 0 (Table 3).

Fits to data and determination of the parameters α
Unbinned maximum likelihood fits were performed for each beam momentum and decay mode individually in order to Table 3 Combinations of the allowed quantum numbers for the reaction p p → ωπ 0 . The J PC combinations even +− and odd −+ are forbidden for thep p system and the combinations even −+ , even ++ and odd ++ are not allowed for the ωπ 0 coupling. The quantum numbers for the ω decay to π 0 γ (L π 0 γ = 1, S π 0 γ = 1) and to π + π − π 0 (L π + π − π 0 = 1, S π + π − π 0 = 0) and S ωπ 0 = 1 for the ωπ 0 coupling are not given explicitly determine the best hypothesis with the resulting fit parameters α J PC Lp p L ωπ 0 . Input for this method are the selected data with the obtained event weights Q i as well as phasespace distributed Monte Carlo events. For properly taking into account the detector resolution and acceptance the GEANT3 transport code has been used. To considering also the correct reconstruction efficiency these Monte Carlo events were then undergoing the same reconstruction and selection criteria as applied for data events and described in Sect. 3. The extended likelihood function L is defined as [16]: where n data denotes the number of data events, τ the phase-space coordinates, α the complex fit parameter, (τ ) the acceptance and reconstruction efficiency at the position τ and n = n data · w(τ , α) (τ ) dτ/ (τ ) dτ . The w(τ , α) represents the transition probability given by Eq. 5. By logarithmizing Eq. 7, approximating the integrals with Monte Carlo events and introducing the weight Q i for each event, the final function to be minimized is then given by: where n MC represents the number of selected Monte Carlo events. Table 4 Best fit results for L max p p for the channelp p → ωπ 0 → (π + π − π 0 )π 0 . The significant improvement in comparison to the hypothesis with L max p p −1 and the marginal improvement of the assumption with L max p p + 1 is a good indication for the unambiguousness of the fit result. The significance is denoted in units of the standard deviation σ  Table 5 Best fit results for L max p p for the channelp p → ωπ 0 → (π 0 γ )π 0 . For further explanations see the caption of Table 4 Momentum To obtain the best hypothesis for the description of the data a strategy has been carried out where fits with successive increase of the maximal contributing orbital angular momentum L max p p have been performed. For each of those fits all allowed waves with Lp p ≤ L max p p have been taken into account. The fit results have been compared using the likelihood ratio. With this strategy it was feasible to determine unambiguously the best hypothesis and thus the largest contributing orbital angular momentum L max p p for all data samples. Summaries of the obtained results are listed in Tables 4 and 5, respectively. Except for the beam momentum of 1525 MeV/c the results for the charged and neutral decay modes are consistent. The slight discrepancy for only one beam momentum is likely caused by the limited acceptance of the detector for the charged decay mode and thus the results for the neutral decay mode are more reliable. For the initial states J PC = even −− and J PC = odd +− two different orbital angular momenta L ωπ 0 = J −1 and L ωπ 0 = J +1 for the ωπ 0 -system are possible (see Table 3). It turned out that for all fits both waves for this system contribute. As an example the fit result with the obtained parameter values are Table 6 Fitted parameters for the reactionp p → ωπ 0 → (π + π − π 0 ) π 0 at 900 MeV/c beam momentum and L max p p = 4. The L,S correspond to the p p, ωπ 0 or ω system, respectively. For technical reasons the α J PC Lp p L ωπ 0 were split into two parameters where some can be fixed due to the linear dependency. Furthermore, two additional phases can be fixed as only relative phases are relevant. The errors are statistical only and originated from the covariance matrix obtained by the fit  Table 6 for the charged decay mode at the beam momentum of 900 MeV/c. The maximal contributing orbital angular momentum L max p p increases continuously from 2 at the lowest beam momentum of 600 MeV/c up to 5 at the highest beam momentum of 1940 MeV/c. These values are in good agreement with a former analysis [17]. Partial wave annihilation cross sections as a function of thep beam momentum for several Lp p values have been estimated [18,19]. Figure 7 shows the outcome of these model calculations for a typical hadronic radius of the baryon core of r 2 B 1/2 = 0.6 f m. Under this assumption the minimump-beam momentum for the production of Lp p = 3 states is expected to be roughly 0.7 GeV, for Lp p = 4, 5 states it is expected to be 1.0 and 1.5 GeV, respectively. The results presented here are in good agreement with these model calculations and only differ in slightly lower momentum thresholds.  [18,19]. The main input parameter for this model is the nucleon radius, which is assumed to be r 2 B 1/2 = 0.6 f m. This figure is extracted from [18]

Comparison of data and fits
The fitted ω-production and ω-decay angles and the normalized λ-value (in case of the charged ω-decay) are compared with the data in Figs. 5 and 6. Apart from minor systematic discrepancies the agreement is good. The reasonable description of the data can also be seen in the fit quality summarized in Table 7. The goodness-of-fit has been estimated with the Pearson χ 2 test based on the histograms for the relevant kinematic variables (Figs. 5, 6) by calculating where n represents the number of the relevant kinematic variables, N bins,i the number of bins for the histogram i, ν i j,data/ f it the number of data/fit entries within bin j for the histogram i and N params the number of fit parameters. The χ 2 values divided by the number of degrees of freedom vary between 0.82 and 1.36 which are reasonable results. However, one has to remark that Eq. 9 does not consider the correlations between the different kinematic variables and thus only serves as a rough estimate for the fit quality.

Spin density matrix of the ω
In addition to the contributing orbital momenta, the polarization observables of the ω meson exhibit important information about its production process. These properties are in general defined by spherical momentum tensors or alternatively by the spin density matrix ρ, which is used in the following. Since the ω is a particle with spin 1 its spin density matrix contains 3 × 3 complex elements ρ λ i λ j , where λ i and λ j represent the helicities of the ω-particle. The ρmatrix is hermitian with a trace of 1 by definition. Polarization means ρ 11 = ρ −1−1 and alignment is defined as ρ 11 = ρ −1−1 = ρ 00 . For measurements with unpolarized protons and antiprotons for channels where the parity is conserved and by choosing the quantization axis to be in the production plane, the number of independent ρ-elements is reduced to four real quantities. The ω spin density matrix for the reactionp p → ωπ 0 is given by [20]: The ρ-matrix elements are dependent on the quantization axis which is here chosen to be the one of the ω-helicity system defined by the ω flight direction in thep p center of mass system. The helicity system is the most suitable one to use for this kind ofp p reactions which is strongly dominated by the s-channel process. In addition, the elements are dependent on the center of mass energy and on the production angle. The determination of the ω-matrix elements has been performed by two different methods: (1) by using the results of the partial wave analysis and (2) solely via the angular decay distributions of the ω-meson. The first method is very rarely used and has already been applied successfully for the reaction γ p → pω [10]. It uses the fitted production amplitude, here defined as T λpλ p λ π 0 r λ ω (p p → ωπ 0 ) (Sect. 4), which contains the information of the ω spin density matrix. The individual ρ-elements can be extracted from the production amplitude by [21]: where N is the normalization factor: According to Eq. 11 the ρ-matrix elements have been projected out from the production amplitude obtained from the partial wave fit of the full reaction chain. In our case (unpolarized initial states) the differential cross section is only dependent on ρ 00 , ρ 1−1 and ρ 10 , so that only these matrix elements can be extracted. The results as function of the center of mass energy and the ω-production angle are summarized in Figs. 8 and 9. The statistical errors have been calculated by propagating the covariance matrix obtained by the likelihood fit. Additionally, a much more time consuming bootstrap approach as described in [22] has been tested which yielded results that are in full agreement to the first calculation.
The second and more traditional method, also called Schilling method, does not make use of the results of the partial wave analysis and uses only the distribution of the ω decay angles θ and φ [20]. The angular distribution for the charged decay mode is given by: and for the neutral decay mode by: As can be seen from Eqs. 13 and 14 only the elements of the real part of the matrix are sensitive to the ω decay angular distribution, which are ρ 00 , ρ 1−1 and ρ 10 . The imaginary part ρ 10 related to an eventual ω-polarization perpendicular to the scattering plane is not accessible. The matrix elements have been extracted separately for different bins in the production angle by fitting the two dimensional ω decay distribution according to Eq. 13 for the charged decay mode and Eq. 14 for the neutral decay mode with a maximum likelihood fit procedure analogous to the one described before in Sect. 4.2. While the two methods rely on different approaches, Spin density matrix elements ρ 00 (first column), ρ 1−1 (second column) and ρ 10 (third column) of the ω in its helicity frame as function of the production angle for the reactionp p → ωπ 0 → (π + π − π 0 )π 0 . While the results obtained via the ω-decay angles are marked with red error bars, the outcome via the partial wave analysis is plotted with continuous black lines. The dashed black lines represent the statistical errors of the partial wave result. Each row represents one specific beam momentum both, however, should yield the same results. Due to the fact that binned data in the production angle are needed for the Schilling method the determination of the ρ-matrix elements with this method is not as accurate as for the first one which naturally imposes all the physical constraints and correlations.
The good agreement between the results obtained with the two different methods can clearly be seen in Fig. 8 for the neutral decay mode and in Fig. 9 for the charged decay mode. It is noticeable that in case of the PWA method the statistical errors are smaller in comparison to the Schilling method. The ρ-matrix elements show a strong oscillatory dependence on the ω-production angle cos(θp p ω ). The ρ 00 and ρ 1−1 have minima and maxima for cos(θp p ω ) = 0 and | cos(θp p ω )| = 0.4, respectively. The ρ 00 -values averaged over the production angle are listed in Table 8. These values show a clear spin alignment effect (ρ 00 = 1/3 would correspond to no spin-alignment). Table 8 ρ 00 -values of the ω meson averaged over the production angle. The averaging is limited on the range of the production angle with a reasonable detector acceptance, which is between −0.85 ≤ cos θp The results for the charged and the corresponding neutral decay mode are in an overall good agreement for all beam momenta. However, differences are visible which are in particular strongly depending on the production angle. These inconsistencies are more significant for the results obtained with the PWA method due to the relatively small statistical errors and might be caused by systematic uncertainties in the simulation and reconstruction procedure.
Similar dependencies on the production angle have already been observed for the tensor polarisation observables of the ω in the same reactionp p → ωπ 0 [3]. In addition the values obtained in the analysis here can be compared with earlier vector meson production experiments inp p -interactions at higher energies. Also there explicit alignment effects for the ρ-meson have been observed [23]. This is in contrast to ppreactions, where negligible alignment for ρ 0 is reported [24]. This trend is also observed in low energy pp-reactions for the orientation of the ω-spin [25].

Summary
The reactionp p → ωπ 0 with unpolarized in-flight data has been analyzed in detail. The ω meson with the neutral decay to π 0 γ as well as with the charged decay to π + π − π 0 has been investigated separately in the low energy regime for variousp beam momenta between 600 and 1940 MeV/c. An excellent background rejection power has been achieved by determining an event based signal weight factor. The performed partial wave analysis has taken into account the complete reaction chain starting from thep p coupling up to the final state particles. It described the data with high precision. The maximal contributing orbital angular momentum L max p p increases continuously from 2 at the lowest beam momentum of 600 MeV/c up to 5 at the highest beam momentum of 1940 MeV/c. The elements of the spin density matrix have been determined with two different methods. The results based on the outcome of the partial wave analysis and those based on the ω decay distributions are in excellent agreement. The first method via the production amplitudes of the PWA was only used in a few cases up to now. The individual elements exhibit a strong dependency on the ω-production angle. A clear spin alignment with ρ 00 values between 0 and 25 % over the whole angular range within | cos(θ )| < 0.9 is visible.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited. Funded by SCOAP 3 / License Version CC BY 4.0.