Exclusive Measurement of the Decay
Abstract
An exclusive measurement of the decay has been performed at the WASA facility at COSY. The mesons were produced in the fusion reaction pdHe X at a proton beam momentum of 1.7 GeV/c. Efficiency corrected differential distributions have been extracted based on 13960140 events after background subtraction. The measured pion angular distribution is consistent with a relative p-wave of the two-pion system, whereas the measured photon energy spectrum was found at variance with the simplest gauge invariant matrix element of . A parameterization of the data can be achieved by the additional inclusion of the empirical pion vector form factor multiplied by a first-order polynomial in the squared invariant mass of the system.
keywords:
meson, box anomaly, exclusive measurementPacs:
14.40.Be, 13.20.Jf, 12.39.Fe1 Introduction
The meson plays a special role in understanding low–energy Quantum Chromo Dynamics (QCD). Chiral symmetry, its realization in hadron physics at low energies and the role of explicit chiral symmetry breaking due to the masses of the light quarks can be investigated using decays. This work focuses on the anomalous sector of QCD, which is manifested in the radiative decays of the meson.
The radiative decay is the fourth strongest decay mode of the meson with a branching ratio of 4.60 0.16% Nakamura:2010zzi . Conservation of charge conjugation and angular momentum including parity constrain the dynamics of the decay products. The photon and the meson are eigenstates of the charge symmetry transformation with the eigenvalues C = –1 and C = +1, respectively. Therefore, due to C invariance the system must have C = –1. To ensure C invariance, the orbital angular momentum between the two pions must be odd. All involved particles have negative parity. Consequently, parity invariance demands that the orbital angular momentum between the photon and the two-pion system must also be odd. Finally, total angular momentum conservation incorporating the fact that the intrinsic spin of the photon is unity leads to the requirement . Thus, the lowest partial waves which contribute are p-waves. Presumably, higher partial waves with are not very important.
Radiative decays of the meson are driven by the chiral anomaly of QCD. The effects of the anomaly have been summarized by Wess and Zumino in an effective Lagrangian Wess:1971 . As shown by Witten, this Lagrangian is an essential part of effective field theories, because it is necessary in order to correctly incorporate the parity transformation of QCD Witten:1983 . At the chiral limit of zero momentum and massless quarks the decay is determined by the box anomaly term of the Wess-Zumino-Witten Lagrangian, which describes the direct coupling of three pseudoscalar mesons and a photon. The dynamic range of the decay is limited by two pion rest masses and the mass, , and is, thus, far from the chiral limit. As a consequence, the decay rate calculated from the box anomaly term at the tree level is smaller by a factor of two compared to the measured value. Higher order terms of the chiral Lagrangian have to be taken into account to achieve a correct description of the decay . Calculations at the one-loop level show an improved agreement between experiment and theory Bijnens:1989 . But there remains, however, a significant difference. Several efforts have been made to include final state interactions by unitarized extensions to the box-anomaly term, e.g. a momentum dependent Vector Meson Dominance (VMD) model Picciotto:1991 , the Hidden Local Symmetry model Benayoun:2003 , an Omnes-function which accounts for the p-wave pion scattering phase shift Venugopal:1998fq ; Holstein:2001 , and a Chiral Unitary approach Borasoy:2004 .
To test the validity of the different models, not only the decay rate but also differential distributions of the Dalitz plot variables need to be compared with experimental data. For this purpose it is useful to parameterize the Dalitz plot in terms of the photon energy in the rest frame of the meson and the angle of the relative to the photon in the pion-pion rest frame. E is related to the squared invariant mass of the pion pair according to
(1) |
The photon energy distribution has been subject of only a few measurements forty years ago Crawford:1966 ; Cnops:1968 ; Gormley:1970 ; Layter:1973 . The results of the two statistically most significant publications Gormley:1970 ; Layter:1973 are presented without acceptance corrections. Instead, the models which have been used for the interpretation of the data were folded with the acceptance and the results seem to be inconsistent Benayoun:2003 ; Borasoy:2007 .
Due to the limitations of the currently available experimental data the potential of the decay to provide insight to the anomalous sector of QCD cannot be fully exploited. In order to perform compelling tests of the box anomaly and its higher order terms in Chiral Perturbation Theory as well as of the non-perturbative extensions of the box-anomaly terms a new measurement with higher precision is called for.
2 The Experiment
The results presented in this paper are based on a measurement with the WASA detection system Adam:2004 installed at the Cooler Synchrotron COSY Maier:1997zj at the Forschungszentrum Jülich in Germany. A pellet target system produces small spheres of frozen hydrogen or deuterium which interact with the ion beam of the accelerator. The interaction region is surrounded by a central detector covering scattering angles from 20 to 169 degrees. It consists of a straw tube drift chamber, which is operated in the magnetic field of 0.85 T provided by a superconducting solenoid for the momentum reconstruction of charged particles, an electromagnetic calorimeter to measure energies of charged as well as neutral particles and thin plastic scintillators to discriminate charged and neutral particles already at the trigger level. Energy loss patterns in the calorimeter and the plastic scintillators allow to identify charged particles. For the identification and reconstruction of particles emitted at polar angles from 3 to 18 degrees a forward detector is used. While track coordinates are measured precisely by a straw tube drift chamber, the kinetic energies of the ejectiles are reconstructed from the signals in plastic scintillators of different thickness, using the E-E method.
The mesons have been produced in the fusion reaction pdHe X at a proton beam momentum of 1.7 GeV/c. This corresponds to an excess energy of 60 MeV in the center of mass for the reaction pdHe at a cross section of 0.412 b Bilger:2002 . At the trigger level events with one track in the forward detector with a high energy deposit in the thin scintillator detectors close to the exit window of the scattering chamber have been accepted. This condition selects He ions without bias on the decay system.
Fig. 1 shows the energy loss correlations of all tracks in the forward detector. The structure attributed to He is clearly visible and can be easily selected. It is well separated from the structures of protons and deuterons, which are strongly suppressed due to the trigger conditions. The right panel of Fig. 1 shows the inclusive missing mass distribution of He. A pronounced peak at the meson mass is visible. It contains 1.1 tagged mesons on top of a continuous background which originates from direct multi-pion production. The cross sections of the relevant background processes pdHe and pdHe are unknown at the center-of-mass energy of this measurement. However, they can be estimated, by extrapolation of results at lower energies Bashkanov:2005fh , to be about 5 b and 0.4 b, respectively. The missing mass spectrum in the right panel of Fig. 1 illustrates that the direct reactions contribute to the background for decays only in a range of a few MeV/c around the meson mass. This remaining background is subtracted in a model independent way, bin by bin from the final spectra.
3 Data Analysis
For the selection of the final state, one charged particle, which is identified as He, is required in the forward detector. In the central detector, two charged particle tracks of opposite curvature and one neutral particle track with an energy deposit in the calorimeter of at least 20 MeV are required in addition. Both charged particles are identified as pions. Background events which are picked up by the selection rules stem mainly from charged multi-pion production and other decay modes.
Two-pion production is expected to form the largest background contribution in the selected data due to the large cross section and the possibility of cluster splitoffs in the calorimeter, which would fake the signal of a photon candidate in the actually photon free final state. Splitoffs are characterized by a small energy deposit and a small distance to either of the pion candidates in the calorimeter. Therefore, this background contributes predominantly to the low end of the reconstructed spectrum. In order to increase the signal to background ratio in this region a correlated condition is imposed on the energy deposits of the photon candidates and the distance between the shower positions of photon and pion candidates in the calorimeter. The latter is measured by the opening angle between the reconstructed cluster positions. In Fig. 2 the condition is shown as graphical cut with a dashed curve. A strong enhancement is observed at low photon energies and small opening angles between photon and pion candidates from the contribution of the pdHe reaction. The cut used is the best compromise between a high signal-to-background ratio and a large reconstruction efficiency of the signal channel. The contribution from two-pion production is reduced by more than 55% by rejecting all events below the curve.
Final states with three pions contribute to the background only if one photon of the decay remains undetected. The can be identified from the distribution of the squared missing mass of the He system, as demonstrated in Fig. 3. Here, the experimental missing mass spectrum is reproduced with Monte Carlo distributions of the the signal () and background () contributions from decays as well as direct two- and three-pion production. The cross sections of the direct pion processes were fitted using simultaneously the He missing mass distribution and the He missing mass for the selected data sample. The result is in agreement with the expectations Bashkanov:2005fh . The contributions of the different meson decay channels are fixed by the known branching ratios Nakamura:2010zzi . The discrepancy between Monte Carlo and experiment for negative masses might be attributed to the unknown production mechanisms of the direct processes. In the simulations isotropic phase space population has been assumed. The discrepancy disappears when the contribution of the direct production is subtracted bin-by-bin, as discussed later. By rejecting events with a squared missing mass value larger than 0.0125 GeV/c, as indicated with a vertical dashed line in Fig. 3, of the background from the He final states in the mass region of the meson is removed.
Additional suppression of background is achieved by a kinematic fit of the complete final state to the reaction hypothesis pdHe , using four-momentum conservation as the only fit constraint. The uncertainties of the kinematic variables at the input to the fit have been extracted depending on energy and angle from a GEANT Monte Carlo simulation, tuned to match the experimental resolutions. After the fit, all events with a probability of less than 10% are rejected. The invariant mass distribution of the fitted system is shown for the remaining events in Fig. 4. The condition on the probability distribution additionally suppresses background. In particular a comparison with Monte Carlo shows a reduction of the He background to the level of in the selected events.
The remaining background is subtracted bin by bin from the and distributions. The bin size of 5 MeV for the photon energy and 0.1 for the pion angular distributions is chosen to reflect the resolution achieved in the respective observables. The bin width of the angular distribution is similar to previous measurements Gormley:1970 ; Layter:1973 the bin width of the spectrum is smaller by about a factor of two. The left panel in Fig. 5 shows the correlation between the photon energy and the invariant mass of the system. Two structures can be seen in the plot. One at the mass of the meson containing the signal events and another structure, showing a correlation of photon energy and invariant mass, caused mainly by background from two-pion production. The invariant mass of the system is calculated for the events in each bin of the and distributions. This is illustrated for one bin of the distribution in the right panel of Fig. 5. Since both, cross sections and differential distributions are unknown for multi-pion production at the energy of this measurement, the individual mass spectra are fitted in two ways in order to determine the amount of continuous background in the mass region. In the first approach a function that is the sum of a Lorentzian for the signal and an exponential for the background is used. In the second method the background is fitted with an exponential function without making an assumption on the signal shape by excluding the range of the signal peak from the fit. The excluded region was determined as the region of a Gaussian fit to the peak of the meson. The individual results of the fits using both methods are shown in the right panel of Fig. 5 with solid and dashed curves, respectively. The average number of the background events from both fits is used to calculate the number of events from the meson decay in the corresponding bin of the or distribution.
The remaining background from the decay in the signal region is subtracted using scaled Monte Carlo distributions. The scaling factors are determined by a fit of the experimental spectrum of the squared missing mass of the He system with Monte Carlo distributions of the relevant decay modes, (signal), and (background), after subtraction of the continuous background from multi-pion production. The contribution of determined to be 7%. This background is isotropic in and contributes to the E distribution in the energy region above 50 MeV with a maximum at 120 MeV. The total statistics in the final distributions is events of the decay. This is the largest number of events from an exclusive measurement of this decay mode.
For acceptance corrections the general form of the squared matrix element for the decay
(2) |
with being the pion momentum in the pion-pion rest frame, and the form factor according to the VMD calculations in Ref. Picciotto:1991 has been used. Using instead Monte Carlo distributions based on the simplest gauge invariant matrix element () does not alter the experimental result significantly. Thus, it can be concluded that systematic effects due to the applied form factors are negligible. The acceptance varies smoothly as a function of and . For photon energies less than 10 MeV the acceptance becomes vanishingly small. In case of the angular distribution a reduced acceptance is observed for small opening angles between each of the pions and the photon. The reduction of the acceptance in both variables is found to be correlated. It is caused by the method of two-pion suppression presented in Fig. 3, where a condition on the correlation of photon energy and opening angle between pion and photon candidates is used.
4 Results
Fig. 6 shows the background subtracted and acceptance corrected photon energy and pion angular distributions with the statistical errors. The E distribution is also given numerically in Tab. 2.
In the upper panel of Fig. 6 the final distribution of is shown. It can be described by , as indicated by the dashed curve. Thus, the measurement is consistent with the relative p-wave assumed in Eq. 2.
The photon energy distribution in the rest frame is shown in the lower panel of Fig. 6. The line shape obtained from Eq. 2, normalized to the integrated rate, is given by the dashed curve. It does not describe the experimental data, which in comparison are shifted significantly towards lower energies. The observed disagreement confirms the findings of the previous measurements Crawford:1966 ; Cnops:1968 ; Gormley:1970 ; Layter:1973 .
In order to achieve a correct description of the photon energy spectrum in addition to the already properly described angular distribution, Eq. 2 can be multiplied by an energy dependent form factor . The origin of the deviation in the E distribution is predominantly given by the final state interaction in the vector channel. Unitarity and analyticity dictate that this effect should be given by the pion vector form factor multiplied by a polynomial that parameterizes contributions that do not contain the unitarity cut. For a detailed discussion about the multiplier of Eq. 2 see e.g. Refs. Venugopal:1998fq ; Holstein:2001 ; Truong:2001en ; Stollenwerk:IP . The pion vector form factor is experimentally directly accessible via or may be derived using the Omnes representation from the elastic phase shifts in the vector channel — here the representation for derived in Ref. Stollenwerk:IP ; Guo:2008nc is applied. The uncertainty in this form factor is negligible compared to the experimental uncertainties, which are presented in detail in the next paragraph. Furthermore, following Refs. Venugopal:1998fq ; Holstein:2001 ; Stollenwerk:IP , we parameterize the term as a real polynomial of first order:
(3) |
The parameter can then be determined from a fit to the data, which is shown with the solid curve in the lower panel of Fig. 6. The result for is shown as the dotted curve.
Tests for systematic uncertainties of have been performed by varying in the analysis chain one by one the conditions for the suppression of splitoffs, the cut on the missing mass of the He system, the condition on the probability of the kinematic fit, the method of subtracting the background from direct multi-pion production, the contribution and the model used for the acceptance correction. Additionally, subsets of the data, collected with different experimental settings, allowed to cross check the influences of luminosity variations and different RF settings during the measurement.
Using different methods to subtract the continuous multi-pion background, different models to perform the acceptance correction or using different RF settings does not cause statistically significant deviations from the original result. The contributions to the systematic uncertainty of derived from the other tests are listed in Tab. 1. The overall systematic error is obtained by adding the contributions quadratically.
One of the largest uncertainties results from the fluctuations of the luminosity during data taking. The variation of the final result with the chosen luminosity can be explained by the accuracy of the simulations concerning pile-up effects and beam-target overlap parameters, which have not been included in a systematic way.
Test | |
---|---|
Splitoffs | 0.34 |
MM(He) cut | 0.22 |
P(,ndf) | 0.12 |
bkg. | 0.26 |
Luminosity effects | 0.32 |
Taking into account the systematic studies, the final result for the parameter is:
where the theoretical uncertainty results from the uncertainty of the pion vector form factor due to the input of Ref. GarciaMartin:2011 for the phase shifts and the extrapolation beyond the upper cutoff. For more details, see Ref. Stollenwerk:IP .
In comparison to theory, calculations based on vector meson dominance Picciotto:1991 ; Benayoun:2003 result in a shape of the differential distribution corresponding to an . The shape given by a parameterization of the pion vector form factor combined with a fit to vector meson dominance Venugopal:1998fq ; Holstein:2001 corresponds to an . The E spectrum from one–loop Chiral Perturbation Theory Bijnens:1989 can be described with an . Thus, the available theory descriptions produce distributions of E, which are close to the curve of shown with a dotted line in the right panel of Fig. 6. Within the total error of the measurement, the value of found in this work appears to be only compatible with the works of Refs. Venugopal:1998fq ; Holstein:2001 .
[GeV] | Entries [a.u.] | stat. [a.u.] |
---|---|---|
0.0225 | 97 | 653 |
0.0275 | 407 | 394 |
0.0325 | 530 | 302 |
0.0375 | 1537 | 252 |
0.0425 | 1291 | 220 |
0.0475 | 2331 | 227 |
0.0525 | 2397 | 218 |
0.0575 | 2841 | 217 |
0.0625 | 3604 | 228 |
0.0675 | 4171 | 235 |
0.0725 | 4887 | 247 |
0.0775 | 5057 | 248 |
0.0825 | 5091 | 252 |
0.0875 | 6444 | 273 |
0.0925 | 6673 | 282 |
0.0975 | 6595 | 282 |
0.1025 | 6448 | 284 |
0.1075 | 6961 | 291 |
0.1125 | 7051 | 297 |
0.1175 | 7242 | 304 |
0.1225 | 7138 | 300 |
0.1275 | 7514 | 312 |
0.1325 | 6963 | 309 |
0.1375 | 7425 | 317 |
0.1425 | 7014 | 319 |
0.1475 | 6892 | 323 |
0.1525 | 6600 | 324 |
0.1575 | 5833 | 315 |
0.1625 | 5739 | 320 |
0.1675 | 4200 | 294 |
0.1725 | 3942 | 292 |
0.1775 | 2985 | 270 |
0.1825 | 2334 | 244 |
0.1875 | 1760 | 217 |
0.1925 | 1250 | 207 |
0.1975 | 384 | 125 |
0.2025 | 63 | 79 |
5 Conclusions
For the first time, background subtracted and acceptance corrected differential distributions of the decay have been extracted in the analysis of exclusive data. The distributions clearly show the importance of final state interactions. The shape of the spectrum can be very well described by a parameterization that includes the factors required by gauge invariance and the centrifugal barrier as well as the pion vector form factor times a first-order polynomial written as . A fit to the data gives . In order to shed further light on the anomalous sector of QCD, future theoretical studies will have to explain simultaneously both the value of as well as the branching ratio for .
In recent production runs of the WASA facility at COSY further data on decays have been taken with high statistics. From a preliminary analysis at least an order of magnitude more fully reconstructed events is expected. The analysis of the acquired data will significantly decrease not only the statistical but also the systematic uncertainties by an improved understanding of background contributions.
The data will also be used to determine the branching ratio of the decay . A recent measurement of the CLEO collaboration Lopez:2007 shows a relative branching ratio which differs by more than three standard deviations from the results of previous measurements Gormley:1970 ; Thaler:1973 . Due to the unbiased tagging of mesons in the reaction pdHe it is not only possible to extract relative but also absolute branching ratios at the WASA facility. This will help to resolve the discrepancy.
6 Acknowledgments
We thank F.-K. Guo for providing us with the code to calculate the pion vector form factor.
This work was in part supported by: the Forschungszentrum Jülich including the JCHP-FFE program, the European Commission under the 7th Framework Programme through the ’Research Infrastructures’ action of the ’Capacities’ Programme. Call: FP7-INFRASTRUCTURES-2008-1, Grant Agreement N. 227431, the German BMBF, the German-Indian DAAD-DST exchange program, VIQCD (VH-VI-231), the German Research Foundation (DFG), and the Polish National Science Centre and Foundation for Polish Science - MPD program, co-financed by the European Union within the European Regional Development Fund.
We gratefully acknowledge the financial support given by the Knut and Alice Wallenberg Foundation, the Swedish Research Council, the Göran Gustafsson Foundation, the Polish Ministry of Science and Higher Education under the grant PBS 7P-P6-2/07.
We also want to thank the technical and administrative staff at the Forschungszentrum Jülich, especially at the COoler SYnchrotron COSY and at the participating institutes.
This work is part of the PhD Thesis of C. F. Redmer.
References
- (1) Particle Data Group Collaboration, K. Nakamura J. Phys. G37 (2010) 075021.
- (2) J. Wess and B. Zumino Phys. Lett. B37 (1971) 95.
- (3) E. Witten Nucl. Phys. B223 (1983) 422–432.
- (4) J. Bijnens, A. Bramon, and F. Cornet Phys. Lett. B237 (1990) 488.
- (5) C. Picciotto Phys. Rev. D45 (1992) 1569–1574.
- (6) M. Benayoun, P. David, L. DelBuono, P. Leruste, and H. B. O’Connell Eur. Phys. J. C31 (2003) 525–547, arXiv:nucl-th/0306078.
- (7) E. P. Venugopal and B. R. Holstein Phys. Rev. D57 (1998) 4397–4402, arXiv:hep-ph/9710382.
- (8) B. R. Holstein Phys. Scripta T99 (2002) 55–67, arXiv:hep-ph/0112150.
- (9) B. Borasoy and R. Nißler Nucl. Phys. A740 (2004) 362–382, arXiv:hep-ph/0405039.
- (10) F. S. Crawford and L. R. Price Phys. Rev. Lett. 16 (1966) 333.
- (11) A. M. Cnops et al. Phys. Lett. 26B (1968) 398.
- (12) M. Gormley et al. Phys. Rev. D2 (1970) 501–505.
- (13) J. G. Layter et al. Phys. Rev. D7 (1973) 2565–2568.
- (14) B. Borasoy and R. Nißler Eur. Phys. J. A33 (2007) 95–106, arXiv:0705.0954 [hep-ph].
- (15) WASA-at-COSY Collaboration, H. H. Adam et al. arXiv:nucl-ex/0411038
- (16) R. Maier Nucl. Instrum. Meth. A390 (1997) 1–8.
- (17) R. Bilger et al. Phys. Rev. C65 (2002) 044608.
- (18) M. Bashkanov et al. Phys. Lett. B637 (2006) 223–228 arXiv:nucl-ex/050811.
- (19) T. N. Truong Phys. Rev. D65 (2002) 056004, arXiv:hep-ph/0105123 [hep-ph].
- (20) F. Stollenwerk et al., arXiv:1108.2419v2 [nucl-th].
- (21) R. Garcia-Martin et al., Phys. Rev. D83 (2011) 074004.
- (22) F.-K. Guo, C. Hanhart, F. J. Llanes-Estrada, and U.-G. Meißner Phys.Lett. B678 (2009) 90–96, arXiv:0812.3270 [hep-ph].
- (23) CLEO Collaboration, A. Lopez et al. Phys. Rev. Lett. 99 (2007) 122001, arXiv:0707.1601 [hep-ex].
- (24) J. J. Thaler, J. A. Appel, A. Kotlewski, J. G. Layter, W. Lee, and S. Stein Phys. Rev. D7 (1973) 2569–2571.