Eta photoproduction in a combined analysis of pion- and photon-induced reactions

Eta photoproduction in a combined analysis of pion- and photon-induced reactions

D. Rönchen Helmholtz-Institut für Strahlen- und Kernphysik (Theorie) and Bethe Center for Theoretical Physics, Universität Bonn, 53115 Bonn, Germany    M. Döring Institute for Nuclear Studies and Department of Physics, The George Washington University, Washington, DC 20052, USA    H. Haberzettl Institute for Nuclear Studies and Department of Physics, The George Washington University, Washington, DC 20052, USA    J. Haidenbauer Institut für Kernphysik and Jülich Center for Hadron Physics, Forschungszentrum Jülich, 52425 Jülich, Germany Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany    U.-G. Meißner Helmholtz-Institut für Strahlen- und Kernphysik (Theorie) and Bethe Center for Theoretical Physics, Universität Bonn, 53115 Bonn, Germany Institut für Kernphysik and Jülich Center for Hadron Physics, Forschungszentrum Jülich, 52425 Jülich, Germany Institute for Advanced Simulation, Forschungszentrum Jülich, 52425 Jülich, Germany    K. Nakayama Institut für Kernphysik and Jülich Center for Hadron Physics, Forschungszentrum Jülich, 52425 Jülich, Germany Department of Physics and Astronomy, University of Georgia, Athens, Georgia 30602, USA

The final state is isospin-selective and thus provides access to the spectrum of excited nucleons without being affected by excited states. To this end, the world database on eta photoproduction off the proton up to a center-of-mass energy of  GeV is analyzed, including data on differential cross sections, and single and double polarization observables. The resonance spectrum and its properties are determined in a combined analysis of eta and pion photoproduction off the proton together with the reactions , , and . For the analysis, the so-called Jülich coupled-channel framework is used, incorporating unitarity, analyticity, and effective three-body channels. Parameters tied to photoproduction and hadronic interactions are varied simultaneously. The influence of recent MAMI and asymmetry data on the eta photoproduction amplitude is discussed in detail.

11.80.Gw, 13.60.Le, 13.75.Gx.

I Introduction

The determination of the spectrum of excited baryons from experimental data is necessary to understand Quantum Chromodynamics at low and medium energies. In this non-perturbative regime, quark models Capstick:1993kb (); Ronniger:2011td () and lattice calculations Edwards:2012fx (); Lang:2012db (); Engel:2013ig () predict more exited states than found so far in partial-wave analyses of experimental data. This dilemma is known as “missing resonance problem” Koniuk:1979vw (). In the past, the dominant source of information on resonance properties was provided by elastic scattering Cutkosky:1979fy (); Hoehler1 (); Arndt:2006bf (). The analysis of inelastic reactions is, however, essential Ronchen:2012eg (); Burkert:2014wea () when aiming at a reliable extraction of the entire spectrum and, consequently, at an identification of missing states that might couple predominantly to channels other than .

Among the inelastic channels accessible in scattering, the channel plays a crucial role. It couples exclusively to states with isospin allowing for the extraction of states unaffected by contributions from states. Moreover, the channel opens at relatively low energies, in a region which is populated by numerous nucleon resonances. For example, some states like the four-star resonance are known to have a large branching ratio. Less well-established resonances like the , whose parameters are only weakly constrained from elastic scattering Arndt:2006bf (); Burkert:2014wea (); Ceci:2011ae (); Ceci:2006ra (), may show a noticeable signal in their decay Ronchen:2012eg (); Burkert:2014wea (). A narrow structure discovered in eta photoproduction on the neutron Kuznetsov:2006kt (); Miyahara:2007zz (); Werthmuller:2013rba () at around  GeV could also appear in eta production on the proton target.

In this respect it is unfortunate that the database for the reaction is problematic over the whole energy range. Its coverage in scattering angles and energies is too limited to perform a well-founded resonance analysis. Furthermore, for energies around 100 MeV from the threshold and beyond, only some experiments were performed and many of those are known to suffer from systematic uncertainties Clajus:1992 (); Arndt:2006bf (); Ronchen:2012eg ().

An alternative experimental window has opened in recent years with high-precision measurements of cross sections and polarization observables in eta photoproduction at photon-beam facilities like ELSA, JLab or MAMI, see Refs. Krusche:2014ava (); Aznauryan:2011qj (); Klempt:2009pi () for reviews. For example, the first measurements from the JLab FROST target (polarization for ) have appeared only recently Strauch:2015zob (), and many more polarization data are expected. The database for eta photoproduction is not yet as large as the one for pion photoproduction, but it is rapidly growing. In addition, the data already available are of much higher quality than those for the pion-induced reaction . Recently, the first data for the beam-target asymmetry in was presented by the A2 collaboration at MAMI, together with a measurement of the target asymmetry  Akondi:2014ttg ().

A key question is whether eta photoproduction data can be used to access the poorly known branching ratios to the channel. The combined analysis of elastic scattering and pion photoproduction determines, at least in principle, the helicity couplings and branching ratios. With known helicity couplings, then indeed the data on eta photoproduction allow one to pin down the resonance branching ratios, without having to resort to the problematic data Clajus:1992 (); Arndt:2006bf (); Ronchen:2012eg (). This argumentation is, however, limited by the fact that even for pion photoproduction the database is not yet fully complete. The extension of the analysis to eta photoproduction will then improve the knowledge of the branching ratios, but, due to incomplete databases, will also lead to changes both in the helicity couplings and branching ratios, as will be seen.

Given the discussed mismatch in data quality, it is not recommendable to only fit model parameters tied to photoproduction and leave the hadronic interaction unchanged: The hadronic amplitude, poorly fixed from the data, appears as a sub-process in photoproduction. Thus, to avoid any bias, an unconstrained fit is needed. In this way, the higher statistical weight of the eta photoproduction data even provides a better constraint on the hadronic amplitude.

The photoproduction of mesons is also a prime candidate for a “complete experiment”. From a mathematical point of view, a complete experiment Barker:1975bp () consists of a set of eight carefully chosen observables, which resolve all discrete ambiguities up to an overall phase Chiang:1996em (); Keaton:1996pe (). For example, a complete set including and is given by  Chiang:1996em (). For experiments with realistic uncertainties, however, eight observables are not sufficient Sandorfi:2010uv (); Vrancx:2013pza (); Nys:2015kqa (). Less than eight observables are required in a truncated partial-wave analysis Omelaenko (); Wunderlich:2014xya ().

By contrast, for a complete experiment on the reaction only four observables are needed Jackson:2013bba (); Jackson:2015dva (). Given the data situation for that reaction, a re-measurement would greatly advance our understanding of the final state. Physics opportunities with a pion beam are discussed in Ref. hadwp (). In any case, the current database does not contain a complete set of observables, neither for photon- nor pion-induced eta production. One then has to resort to other approaches that often combine data from different initial and final states.

Over the years, a variety of theoretical approaches has been applied to analyze the pion- and photon-induced production of mesons. For example, photoproduction of mesons in the resonance region not too far from threshold was studied in the framework of unitarized chiral perturbation theory in Refs. Kaiser:1996js (); Inoue:2001ip (); Borasoy:2002mt (); Doring:2009qr (); Mai:2012wy (). Considering a broader energy range up to and beyond 2 GeV, -matrix Arndt:2005dg (); Anisovich:2012ct (); Shklyar:2012js (); Shrestha:2012ep (); Batinic:2010zz (); Cutkosky:1979fy (); Gridnev:1999sz () and unitary isobar Chiang:2001as () models are practical tools to perform an analysis of large amounts of data. Sometimes, the real part of the self-energies is neglected and only on-shell intermediate states are maintained, which reduces the complexity of the calculations. For the purpose of a combined analysis of different reactions over a wide energy range, so-called dynamical coupled-channel (DCC) models provide a particularly suited framework. Theoretical constraints of the -matrix, like two- and three-body unitarity, analyticity, left-hand cuts and complex branch points, are manifestly implemented or at least approximated. This enables the reliable determination of the resonance spectrum in terms of pole positions, residues, and helicity couplings in the complex energy plane. The production of mesons in DCC approaches was studied, e.g., in Refs. Kamano:2013iva (); Chen:2007cy ().

Here, we extend the Jülich model, a DCC approach pursued over many years Doring:2009bi (); Doring:2009yv (); Doring:2010ap (); Ronchen:2012eg (); Ronchen:2014cna () starting with Ref. Schutz:1994ue (), to perform a simultaneous analysis of the pion-induced reactions , , , , , , and the photon-induced reactions , , and . We allow the hadronic amplitudes themselves to vary, in addition to the parameters tied to photoproduction. As discussed, this is necessary because the quality of the data in is much inferior to the data in eta photoproduction. In a simultaneous fit to all pion- and photon-induced data, we observe that, indeed, the eta photoproduction data have a strong influence. This influence reaches beyond the electromagnetic resonance properties and affects also resonance pole positions and hadronic branching ratios.

The paper is organized as follows: In Sec. II, a short overview of the applied formalism is given. For a more detailed introduction of the semi-phenomenological approach to meson photoproduction we refer the reader to Ref. Ronchen:2014cna (). In Ref. Ronchen:2012eg (), an extensive description of the hadronic Jülich DCC framework is provided. In Sec. III.1 we describe the data analysis and in Sec. III.2 the fit results are shown. The extracted resonance parameters are presented and discussed in Sec. IV. Technical details about the renormalization of the nucleon mass are summarized in an appendix.

Ii Formalism

In the approach (referred to as “Jülich model”), the hadronic scattering potential is iterated in a Lippmann-Schwinger equation formulated in time-ordered perturbation theory (TOPT) and two-body unitarity is, thus, automatically fulfilled. The three-body states are parameterized through the channels , and . These effective three-body channels are included dynamically, i.e., the and subsystems match the corresponding phase shifts Doring:2009yv (). The analytic structure of the amplitude is given through real and complex branch points Ceci:2011ae () and the real, dispersive contributions of the intermediate states. Moreover, - and -channel exchanges of known mesons and baryons constitute the non-resonant part of the amplitude and serve as “background”. While the -channel diagrams approximate the left-hand cuts, -channel meson exchanges are essential to achieve three-body unitarity Aaron:1969my (). Note that the latter is, at the moment, only approximately satisfied in the Jülich model. By means of this explicit treatment of the background, strong correlations between the different partial waves and a non-trivial energy and angular dependence of the observables are generated. Although - and -channel processes are necessary for analytic structure and unitarity, they do not fully determine the amplitude. Bare resonance states are included as -channel processes. In contrast to previous versions of the approach, here we also allow for additional contact interactions. Such interactions do not spoil the analytic properties ensured by -, - and -channel interactions. They absorb physics beyond the explicit processes and, thus, increase the model-independence of the approach at the cost of a few more parameters. Practically, the changes in the amplitudes induced by the contact terms are comparatively small and the so-called background is still dominated by the - and -channel exchanges. Details are given in the following section. Note that contact terms are also included in Ref. Jackson:2015dva ().

In Ref. Doring:2010ap () the approach was extended to the strangeness sector incorporating the final state in the analysis. In the Jülich2012 model of Ref. Ronchen:2012eg (), the spectrum of nucleon and resonances was extracted from a simultaneous analysis of the reactions , , and . The extension of the Jülich approach to pion photoproduction in a field-theoretical formulation, that respects the generalized off-shell Ward-Takahashi identity, was achieved in Ref. Huang:2011as (). In Ref. Ronchen:2014cna (), by contrast, the photon interaction is approximated in a phenomenological framework and the Jülich2012 analysis serves as final-state interaction. The flexible formulation of Ref. Ronchen:2014cna (), used to study the world data on pion photoproduction on the proton, proved to be capable of analyzing large amounts of data while at the same time maintaining the analytic properties of the Jülich approach. This framework, with the addition of contact terms, will be applied in the present study.

ii.1 Pion-induced reactions

The pion-induced reactions are treated within the Jülich dynamical coupled-channel formalism Ronchen:2012eg (). The -matrix which describes the scattering process of a baryon and a meson can be formulated in the partial-wave basis and reads


In Eq. (1) and in the following, always means the scattering energy in the center-of-mass (c.m.) frame and () is the modulus of the outgoing (incoming) three-momentum that can be on- or off-shell. The channel indices , and represent the incoming, outgoing and intermediate meson-baryon pairs, respectively.

The propagators for channels with stable particles, i.e. , , , or , are given by


with and being the on-mass-shell energies of the intermediate particles and in channel with respective masses and . In case of the channels with unstable particles , and that parameterize the channels, the propagators are more complex; for details see Ref. Doring:2009yv ().

The scattering potentials can be decomposed into a pole and a non-pole part


The quantity denotes the sum of all - and -channel exchange diagrams, while comprises the -channel resonance graphs. The functions () correspond to the bare creation (annihilation) vertices of a resonance with bare mass and are constructed from an effective Lagrangian which can be found in Table 8 of Ref. Doring:2010ap (). Explicit expressions of and are also given in Ref. Doring:2010ap (), cf. also Appendix A of Ref. Ronchen:2012eg (). The exchange potentials that constitute , also derived from effective Lagrangians, are compiled in Appendix B of Ref. Ronchen:2012eg (). The decomposition of Eq. (3) is slightly modified in the present approach. We implement here contact terms that do not introduce any singularities and that are used to absorb physics not explicitly contained in the parameterization through -, - and -channel processes. The contact terms are introduced in a separable form and separately for every partial wave,


where the () have the same functional form as the resonance vertices () in Eq. (3). The associated couplings in the are now new free parameters to be adapted in the fit to the data.

Formally, the numerator structure of Eq. (4) is the same as the one of -channel pole terms, to ensure the correct threshold behavior. Also, the contact terms carry channel indices. Formally, we can treat the contact terms as bare -channel processes and absorb the contributions in the definition of :


For a compact notation, we will no longer distinguish between bare resonance vertices and structures of the contact term in the following. Then, the index can refer to a bare -channel resonance vertex or one of the terms in the numerator of Eq. (4).

Similar to the potential , the scattering matrix can also be written as the sum of a pole and a non-pole part,


with the unitary defined through


Here and in the following, we do not display function arguments and integration symbols in favor of a more compact notation. Besides the resonance poles arising from -channel diagrams in , a dynamical generation of poles in is also possible, as explained in Refs. Doring:2009bi (); Ronchen:2012eg ().

Note that the pole part can be evaluated from the non-pole part . To this purpose one defines the dressed creation (annihilation) vertex () via


where is the self-energy. The indices and label the -channel states or a contact diagram in a given partial wave.

The pole part is then given by


with the resonance propagator . For example, if there are two -channel resonances with masses and (indices ) plus one contact term (indices ), the expression reads

The decomposition into pole and non-pole part performed here has mostly technical reasons: The numerical evaluation of the non-pole part is much more time-consuming than the evaluation of the pole part. This leads to an effective, nested fitting workflow as discussed in detail in Refs. Ronchen:2012eg (); Ronchen:2014cna (). Here, we have introduced the contact terms technically on the same footing as resonances, which allows to fit those terms computationally more effectively.

In Sec. 2.2 of Ref. Ronchen:2012eg (), the renormalization of the nucleon mass and coupling in the presence of two bare -channel states in the partial wave was derived. With the implementation of contact diagrams, this procedure has to be extended to two resonances and one contact diagram. This is addressed in Appendix A.

ii.2 Photon-induced reactions

A field-theoretical description of the photoproduction amplitude within a gauge-invariant framework that respects the generalized off-shell Ward-Takahashi identity Haberzettl:2011zr (); Haberzettl:2006bn (); Haberzettl:1997jg () was successfully applied in the analysis of pion photoproduction in Ref. Huang:2011as (). In this study an earlier version of the Jülich model was utilized to provide the hadronic final-state interaction. This field-theoretical method allows one to gain insight into the microscopic reaction dynamics of the photo-interaction. By contrast, the semi-phenomenological approach to pseudoscalar meson photoproduction, developed in Ref. Ronchen:2014cna () and used here, is more flexible and facilitates the analysis of a large amount of data. Here, the photo-interaction kernel is approximated by energy-dependent polynomials, while the hadronic final-state interaction is provided by the Jülich DCC model described in Sec. II.1. The formalism is inspired by the GW-SAID CM12 parameterization Workman:2012jf () and will be applied in the present study. Nonetheless, we consider the present analysis as an intermediate step towards an expansion of the field-theoretical framework of Ref. Huang:2011as (). A detailed introduction to the semi-phenomenological approach was given in Sec. 2.2 of Ref. Ronchen:2014cna (). In the following, we recapitulate the basic elements.

The multipole amplitude of the photoproduction process is given by


where the index denotes the initial state. is the hadronic half-off-shell -matrix introduced in Sec. II.1 with the intermediate (final) meson-baryon channel () and the corresponding off-shell momentum (on-shell momentum ). Integration over the intermediate off-shell momentum similar to Eq. (1) is suppressed here in the notation, following the convention of Eq. (7). In the present analysis the photon is allowed to couple to the intermediate channels , and , while we have and as final states .

The photoproduction kernel is written as


Here, stands for the photon coupling to the non-pole part of the photoproduction kernel. The tree-level coupling of the channel to the nucleon and resonances is represented by the vertex function , where denotes the resonance number in a given partial wave. The hadronic resonance annihilation vertex is exactly the same as in Eq. (3), which results in the cancellation of the explicit singularity at . Both quantities, and , are approximated by energy-dependent polynomials ,


In Eq. (13), the vertex function has the same form as in Eq. (12) but without any dependence on the resonance number . The polynomials , for a given multipole, are parameterized as

with and being free parameters that are fitted to the data. The upper limits of the summation indices , and , are chosen so as to permit a good data description, but are restricted to be less than 4. In order to fulfill the decoupling theorem, which states that resonance contributions are parametrically suppressed at threshold, the summation for starts with . The expansion point is chosen to be close to the threshold,  MeV. In this way, the factor absorbs the potentially strong energy dependence at the threshold, which is not too far from the threshold. Moreover, this factor guarantees a well-behaved multipole amplitude in the high-energy limit, although a more quantitative matching to Regge amplitudes remains to be done Huang:2008nr ().

In order to achieve a good description of the high-precision data for pion photoproduction close to threshold, we take into account some isospin breaking effects, i.e. we apply different threshold energies for the and the channels, as explained in Sec. 2.3 in Ref. Ronchen:2014cna (). Similarly, we take the physical threshold of the final state when calculating observables. Note that, in general, isospin-averaged masses are used in the Jülich model.

A multipole decomposition of the photoproduction amplitude of pseudoscalar mesons can be found in Appendix A of Ref. Ronchen:2014cna ().

Iii Results

iii.1 Database and free parameters

The database of the present study comprises the hadronic data used in Ref. Ronchen:2012eg (). This represents the world database on the reactions , and up to  GeV, plus the WI08 energy-dependent solution of the GWU/INS SAID group Workman:2012hx (). In addition, we include almost all published data on pion photoproduction off the proton up to  GeV, with some forward regions at high energy excluded Ronchen:2014cna (). Third, in this study we add the world database for the reaction , again up to  GeV. These data were taken from the GW-SAID database Workman:2012hx ().

Fit A Fit B
PWA GW-SAID WI08 Workman:2012hx ()
, ,
, ,
6000 data points
, , , , , ,
, , , , , ,
, , , , , ,
29,392 data points 29,680 data points
Table 1: Data included in fits A and B. In contrast to fit A, fit B contains the recent MAMI measurements of and  Akondi:2014ttg ().

In a first fit, called fit A in the following, we include all data except the recent eta photoproduction measurement of the transverse target asymmetry and the beam-target asymmetry by the A2 collaboration at MAMI Akondi:2014ttg (). These data are added in a second fit, fit B. Performing these two fits allows us to estimate the influence of the new polarization data on the extracted resonance spectrum. An overview of the data included in the fits can be found in Table 1.

Compared to the high-precision data nowadays available in case of pseudoscalar meson photoproduction, the data situation for the pion-induced reactions is difficult in large parts, the lack of polarization measurements being one of the major issues. In Sec. 3 of Ref. Ronchen:2012eg () the situation for the individual hadronic channels was discussed in detail. In the present study, we adopt the systematic errors and mainly also the special weights applied to certain data sets.

Moreover, we continue along the lines of Ref. Ronchen:2014cna () and apply an additional systematic error of 5 % to all photon-induced data in order to account for discrepancies in the data. The amount of 5% is an estimate at this point. More advanced techniques have been applied, e.g., by the GWU/SAID group allowing for normalization corrections Arndt:2006bf (); Workman:2012jf (). We plan to improve our analysis along these lines in the future. To compensate for the smaller number of data points for and in eta photoproduction, those data are weighted in fit B with generic factors around 10 as found necessary for obtaining satisfactory fit results. To achieve a good description of the data at higher energies, additional weights have to be applied. The situation is similar in case of the beam asymmetry in . Compared to the number of data points for the differential cross section, 5680, only a few are available for , namely 189. Furthermore, we did not attempt to achieve a good description of the recoil polarization for as only seven data points are available and their influence on the fit is very limited. The data are, however, included in the fit but no special weights were applied.

In the Jülich approach the free parameters tied to the hadronic interaction are the bare coupling constants and masses of the -channel resonances, the strengths of the contact terms, and the cut-off parameters in the - and -channel diagrams constituting . In addition, there are certain couplings in some of the latter diagrams that cannot be connected to other coupling constants via SU(3) flavor symmetry and have to be fitted to data as well. In the present study we do not alter the parameters but employ the values found in the Jülich2012 analysis of pion-induced reactions Ronchen:2012eg ().

We do, however, vary the parameters tied to hadronic resonances and contact terms. This is a necessity in the current situation, given that the data of pion-induced eta production is of less quality than the eta photoproduction data. The higher-quality photoproduction data can, thus, help to constrain the hadronic amplitude. In both fits A and B the number of -channel resonance parameters amounts to 128. For each of the 11 genuine and 10 genuine resonances (i.e., those resonances included in form of a pole diagram in the potential Doring:2010ap (); Ronchen:2012eg ()) the parameters are given by one bare mass and the couplings to the channels , , , , and as allowed by isospin. Contact terms introduced in Eq. (4) are switched on in both fits in the and partial waves for the and channel, and in the partial wave for the , , , and channel, giving rise to 9 additional fit parameters. The free parameters that are used to tune the interaction of the photon with hadrons are the resonance parameters and and the non-pole couplings and with , or , cf. Eq. (LABEL:polys). Formally, all fit parameters up to are implemented in the computer code. However, the actual number of fit parameters is chosen as required by data and the remaining parameters are set to zero. In order to achieve a good description of and in fit B, additional fit parameters that were set to zero in fit A had to be released. Thus, we have 443 parameters tied to the photon interaction in fit A and 456 in fit B. In total, the number of fit parameters adds up to 580 in fit A and 593 in fit B.

The free parameters are adjusted to the data in simultaneous fits of all pion- and photon-induced reactions using MINUIT on the JUROPA supercomputer at the Forschungszentrum Jülich.

iii.2 Fit results

In the following, only selected fit results for the reactions and are shown. Data sets with energies that differ by less than 5.5 MeV are sometimes displayed in the same graph. The full fit results for all pion- and photon-induced reactions included in this analysis can be found online Juelichmodel:online ().

The definition of the various photoproduction observables in terms of the multipole amplitudes is given in Appendix B of Ref. Ronchen:2014cna (). The convention agrees with the one of the SAID group Sandorfi:2011nv (); Arndt:2002xv ().


In Figs. 1 and 2 we show selected results for the differential cross section. In the threshold region only very small differences can be observed between fits A and B. Starting at  MeV, the two fits show some discrepancies at extreme angles, fit A providing a slightly better description of the data at medium energies, cf.  MeV in Fig. 1. At higher energies, the differences are most apparent at very forward angles. However, the data situation does not allow for an assessment regarding which of the fits is best. This can be seen, e.g., at  MeV in Fig. 2 where both fits seem to describe the data equally well.

Figure 1: Differential cross section of the reaction . Dashed (blue) line: fit A; solid (red) line: fit B; data: MC10 McNicoll:2010qk (), BA07 Bartalini:2007fg (), KR95 Krusche:1995nv (), NA06 Nakabayashi:2006ut (), DY95 Dytman:1995vm (), CR05 Crede:2003ax (), DU02 Dugger:2002ft (), CR09 Crede:2009zzb (), WI09 Williams:2009yj (). In this and all following figures, the numbers in the plots specify the pertinent center-of-mass energy [MeV].
Figure 2: Differential cross section of the reaction . Dashed (blue) line: fit A; solid (red) line: fit B; data: cf. Fig. 1 and SU09 Sumihama:2009gf ().
Figure 3: Beam asymmetry of the reaction . Dashed (blue) line: fit A; solid (red) line: fit B; data: BA07 Bartalini:2007fg (), ER07 Elsner:2007hm ().

The situation is similar in case of the beam asymmetry in Fig. 3. Disagreements of the fits A and B are visible predominantly at higher energies and at forward angles where there are no data.

Figure 4: Recoil polarization of the reaction . Dashed (blue) line: fit A; solid (red) line: fit B; data: HE70 Heusch:1970tr (), HO71 Hongoh:1971pv ().

Only seven data points are available for the recoil polarization , cf. Fig. 4. Although the results of our two fits are very different, especially at higher energies, both describe the data more or less well but with some deficiencies. A larger database of this observable may help constrain the partial wave content.

Figure 5: Transverse target asymmetry and beam-target asymmetry in the reaction . Dashed (blue) line: prediction of fit A; solid (red) line: fit B; data: Ref. Akondi:2014ttg ().

The fit results for the transverse target asymmetry and the beam-target asymmetry can be found in Fig. 5. These data were only included in fit B, meaning that the blue dashed line in Fig. 5 (fit A) represents a prediction for these observables. Note that older data for from Ref. Bock:1998rk (), that are partially in conflict with the MAMI data, were not fitted. The prediction of fit A at lower energies  MeV is acceptable for and good for , while the shortcomings of the prediction are considerable at higher energies. However, once the data are included (fit B) a good description over the whole energy range is achieved.


In Fig. 6 we show selected fit results for the differential cross section of the reaction . In addition to fits A and B, we display results from the old fit A of the Jülich2012 analysis Ronchen:2012eg (), in which only pion-induced reactions were analyzed. We call the latter fit A in the following, where the subscript serves as a reminder that the fit in Ref. Ronchen:2012eg () included only hadronic data. Since the resonance vertex function appears in the construction of the hadronic amplitude in Eq. (3) and also in the photoproduction kernel in Eq. (12), the bare parameters of have influence on the pion- as well as on the photon-induced production of the final state. As much less data are available for compared to , the photon-induced eta production does impose constraints on the fit results of the hadronic channel.

As can be seen in Fig. 6, starting already at energies of  GeV, data sets from different experimental groups show conflicting behavior. Moreover, at energies  MeV all data for the differential cross section of stem from a measurement Brown:1979ii () deemed problematic due to a miscalibration of the beam momentum, see Ref. Clajus:1992 () for details. Those data enter the fit with a much reduced weight and all three fits yield different results. This is most apparent at higher energies, where fits A and B are dominated by photon-induced data. Here, fits A and B mostly differ at extreme forward angles while for other angles they coincide reasonably well. This discrepancy has been also noted for the differential cross section and beam asymmetry in fits A and B in eta photoproduction, mostly due to less precise data at forward angles.

In Fig. 7 we present results of the recoil polarization in . In principle, the only published data for this observable Baker:1979aw () exhibit the same problems as the differential cross sections of Ref. Brown:1979ii () because the same experimental set up was used. Those data are fitted with a very low weight, too. As in the case of the differential cross section, the differences in the three fits are obvious. However, in view of the discussed quality issues with the available data for , we see no reason to enforce a better data description.

The comparison made here demonstrates the need for a pion beam to re-measure the reaction . Note that only four observables are needed for a complete experiment, as discussed in the Introduction.

Figure 6: Differential cross section of the reaction . Dashed (blue) line: fit A; solid (red) line: fit B; dash-dotted (green) line: fit A of Ref. Ronchen:2012eg (). Data: filled circles from Ref. Bayadilov:2008zz (); filled squares from Ref. Prakhov:2005qb (); empty triangles up from Ref. Kozlenko:2003hu (); stars from Ref. Debenham:1975bj (); empty squares from Ref. Richards:1970cy (); empty diamonds from Ref. Brown:1979ii (). Note that the data situation for this reaction is problematic, see text.
Figure 7: Polarization of the reaction . Dashed (blue) line: fit A; solid (red) line: fit B; dash-dotted (green) line: fit A of Ref. Ronchen:2012eg (). Data: Ref. Baker:1979aw (). Note that the data situation for this reaction is problematic, see text.

iii.3 Multipoles

In Fig. 8 the electric and magnetic multipoles from fits A and B for the reaction are shown. Additionally, we display the multipoles of the Bonn-Gatchina BG2014-02 solution Gutz:2014wit (). Note that those amplitudes were not included in the fit.

Generally speaking, the multipoles extracted in our analysis and the ones of the Bonn-Gatchina group exhibit large differences. An exception is the and, to a certain degree, also the multipole. While for differences between our fits A and B are hardly noticeable at all, the multipole features clearly visible deviations in the two fits at energies higher than 1650 MeV. As will be discussed in Sec. IV.1, this is the energy regime where the pole of the resonance is located. Among the lower multipoles, the shows the most striking discrepancies between fits A and B, see also the discussion on the influence of variations in the partial wave on the description of and in Sec. IV.1. For higher multipoles fit B sometimes shows a stronger energy dependence than fit A, cf. , and , . In summary, the new MAMI data for and have a large impact on the multipole amplitudes.

We observe that for lower partial waves, the eta photoproduction multipoles of the present study exhibit less agreement with the Bonn-Gatchina multipoles than our pion photoproduction multipoles with the Bonn-Gatchina BG2014-02 Gutz:2014wit () or the GW-SAID CM12 solution Workman:2012jf (). This suggests that the multipole content of the reaction is much less established than in the case of pion photoproduction where the various analyses agree better. Figures showing the pion photoproduction multipoles can be found online Juelichmodel:online ().

Figure 8: Electric and magnetic multipoles for the reaction . (Black) dash-dotted line: BG 2014-02 solution Gutz:2014wit (). Dashed (blue) line: fit A; solid (red) line: fit B.

Iv Resonance spectrum

A resonance state is uniquely defined by its pole position in the complex energy plane, the residues associated with the channel transitions, and the Riemann sheet the pole is located on. With the exception of the physical sheet of the lowest lying channel, the poles can appear on various Riemann sheets, but not all of them are of physical interest. Usually, only the poles on the sheet which is closest to the physical axis are considered. We select this sheet by rotating the right-hand cuts of all channels in the direction of the negative imaginary axis. In this way, we define the second sheet where all poles extracted in the present study lie. See Ref. Doring:2009yv () for a detailed discussion.

In order to determine the pole positions, the scattering amplitude has to be continued to the second Riemann sheet. For this purpose we apply the method of analytic continuation following Ref. Doring:2009yv () where the amplitude on the second sheet is accessed via a contour deformation of the momentum integration. The calculation of the residues proceeds via the formalism illustrated in Appendix C of Ref. Doring:2010ap (). Definitions of the normalized residue and the branching ratio into a specific channel are given in Sec. 4.1 of Ref. Ronchen:2012eg (). In case of the latter two quantities we use the same definitions as the Particle Data Group Agashe:2014kda (). For a reliable extraction of the resonance parameters, the correct structure of branch points, including the complex branch points of the channels including unstable particles , and , is crucial. In Ref. Ceci:2011ae () it was shown that the absence of the latter might lead to false resonance signals.

The definition of the photocouplings at the pole ,


can be found in Appendix C of Ref. Ronchen:2014cna () and is identical to the definition given in Ref. Workman:2013rca (). The photocoupling at the pole characterizes the coupling of the channel to a resonance, independently of the final state in the reaction under consideration. Note that, in general, the complex cannot be compared to the real-valued helicity amplitudes , see Sec. D of Ref. Ronchen:2014cna () for further remarks.

In Tables 2 to 4 we list the pole positions, residues and the photocouplings at the pole of the present study. In addition to the values extracted in the current fits A and B we list the pole positions and residues found in fit A of the Jülich2012 analysis Ronchen:2012eg (), called fit A in the present study, and the photocouplings of fit 2 from the Jülich2013 analysis Ronchen:2014cna (). Note that in the latter study, the parameters of the hadronic -matrix were not altered, i.e. the resonance pole positions and hadronic residues are the same as in fit A of Ref. Ronchen:2012eg (). An overview of the pole positions of fit A and B is also given in Fig. 9.

In Table 3 the channel labeled (6) corresponds to the case where and the one labeled (7) to . Also the channel can couple to a resonance with a given in multiple ways, cf. Table 11 in Ref. Ronchen:2012eg (). Here, we only quote normalized residues for , since at energies well above the threshold this channel can be regarded as being composed of the two stable particles and . In general, the resonance coupling at the pole to a channel like is a function of the center-of-mass momentum of the stable particle (that equals the summed momenta of the decay products of the unstable particle). Here, we do not quote this function of but choose as on-shell three-momentum of a stable of mass  MeV and a pion. Obviously, this prescription does not lead to meaningful results for the very broad in the channel, or the channel. In the latter channel, most resonances are not far above the threshold that is situated around  GeV, and the cannot be considered a stable particle.

Figure 9: Pole positions of the isospin =1/2 (above) and =3/2 (below) resonances extracted from fit A (filled squares) and from fit B (empty diamonds). For a better differentiation, the names and pole positions of some resonances are marked with different colors. The green crosses denote the branch points of the amplitude. Note that all cuts, starting at the branch points, are chosen in the negative Im direction.

iv.1 Discussion of specific resonances

Compared to the earlier analysis of Ref. Ronchen:2012eg (), the extension of the model to eta photoproduction did not require the inclusion of additional bare -channel states, and we find no new dynamically generated resonances either. In particular, there is no need to include a narrow state at around  GeV. The narrow structure discovered in eta photoproduction on the neutron Kuznetsov:2006kt (); Miyahara:2007zz (); Werthmuller:2013rba () is absent in the present analysis of eta photoproduction on the proton.

In the following, when discussing selected resonance states, we always refer to the values quoted in Tables 2 to 4.

: While the real part of the pole position of the 1/2 is very stable throughout all three fits A, B and A, the width is by 30 MeV larger in the new fits that include eta photoproduction data. The new value is close to the one found in a recent analysis Svarc:2014aga () of the GW-SAID WI08 solution Workman:2012hx (), where elastic and data were fitted. Also the mass in the latter analysis is very similar to our fits. We obtain the same values for the normalized residue in all three fits. While the elastic residue is larger in the new fits, the coupling to this channel is still considerably smaller than the coupling to . In both current fits A and B, the magnitude of the photocoupling is more than twice as large as in the previous fit 2 of Ref. Ronchen:2014cna () where only pion photoproduction data were considered. This change is related to the increase of the width for the , because the resonance width and size of the photocoupling at the pole are strongly correlated. On the whole, comparing the resonance parameters of the of the older fits to the new ones, the inclusion of eta photoproduction seems to have noticeable impact for this resonance. The influence of the new and data Akondi:2014ttg (), on the other hand, is rather limited as the parameters in fit A and B do not exhibit major differences. This observation is in agreement with the similarity of fit A and B for and at  GeV (cf. Fig. 5 ).

By contrast, the pole of the second resonance in the partial wave, , is located in an energy region where the deficiencies of the prediction of fit A for and become more apparent. Accordingly, slightly larger variations are found in the pole positions of fits A and B, and also compared to fit A. Although the width is smaller in the new fits A and B, the resonance is still broader and has a higher mass than in Ref. Svarc:2014aga (). In Ref. Svarc:2014sqa (), however, the results from a Laurent-Pietarinen (L+P) expansion of the GW-SAID CM12 solution Workman:2012jf () for pion photoproduction give a pole position of  MeV, which is closer to our values. As in fit A, the current fits reveal a strong coupling to the channels. The magnitude of the photocoupling to the is more than twice as large in the fits including eta photoproduction compared to the older fit 2 of Ref. Ronchen:2014cna (), where only pion photoproduction was considered. While fit 2 yielded a value smaller than the one of Ref. Svarc:2014sqa (), the photocoupling is now larger in fit A and B.

In summary, the pole positions of the and the come closer to the values of the GW-SAID analysis. As resonance widths and sizes of helicity couplings are correlated, the values of the latter also approach the values of that analysis. The eta photoproduction data have a strong influence on the resonance properties.

: Besides the nucleon pole, we find two resonances in the partial wave. One of them, the Roper resonance 1/2 is dynamically generated from the interplay of the - and -channel diagrams. As the fit parameters corresponding to these diagrams are not altered in the present study the extracted resonance parameters do not change much. The situation is different for the third state, the 1/2. This explicit -channel state was introduced in the Jülich model in Ref. Ronchen:2012eg () mainly to improve the description of the pion-induced and channels. Since it couples only weakly to the channel its resonance parameters are poorly constrained from elastic scattering. As can be seen in Table 2, the extension of the fit to new, inelastic reaction channels results in noticeable changes in the pole position. Moreover, also the inclusion of new observables for a specific reaction, here and in in fit B, leads to significant variations not only in the pole position but also for the residues and photocouplings. The latter observation suggests that additional information from inelastic channels, e.g. in form of new polarization measurements, might help to fix the parameters of the 1/2.

In all our fits, the 1/2 has a lower mass and is narrower than in recent analyses by the ANL-Osaka ( MeV) Kamano:2013iva () and the Bonn-Gatchina groups ( MeV) Anisovich:2011fc (). In the L+P analysis of the GW-SAID CM12 solution in Ref. Svarc:2014sqa (), a broad state with a higher mass associated with this resonance is found that can be alternatively explained as the complex branch point. However, the authors state that additional information from other decay channels beside is needed to distinguish between the two options. See also Ref. Svarc:2014aga () by the same authors where the same conclusion was drawn. Note that in the present study the complex branch point is included explicitly.

In addition to the 1/2 and the 1/2 we find non-conclusive indications for another, very broad and dynamically generated pole at  GeV.

: We include one bare -channel state in the partial wave, the 3/2. Although we observed a noticeable sensitivity of the description of the channel on variations in the partial wave, the pole position of the 3/2 is very similar in fits A and B. The impact of the new and data from MAMI can be seen in the photocouplings of this state (fit B vs. fit A). This is reflected in the discrepancies observed in the multipole in Fig. 8.

In different GW-SAID solutions Svarc:2014sqa (); Svarc:2014aga () and in the Bonn-Gatchina analysis of Ref. Anisovich:2011fc () the 3/2 has a pole position with a real part 20 to 80 MeV lower than in our fits and an imaginary part more than 50 MeV larger. By contrast, the ANL-Osaka group Kamano:2013iva () finds values closer to ours.

We tested the influence of a second explicit resonance state in the partial wave but observed no significant improvement of the fit results. In this partial wave, the Bonn-Gatchina group finds strong evidence for a state named in the photoproduction of and  Nikonov:2007br (); Anisovich:2011fc (). It has also been confirmed in in an effective Lagrangian model Maxwell:2012zz () and in a covariant isobar-model single channel analysis Mart:2012fa (). It remains to be seen, whether this state is also needed in the Jülich approach once the analysis is extended to kaon photoproduction. Note that the is also included in the ANL-Osaka analysis Kamano:2013iva () and in the Gießen model Penner:2002ma ().

: While the real part of the pole positions of the 3/2 is unchanged in fits A and B, the imaginary part is about 10 MeV smaller in fit B. In the previous fit A the real and the imaginary part were similar to fit A, -2Im =110 MeV. In the GWU-SAID solutions analyzed in Refs. Svarc:2014sqa (); Svarc:2014aga () and the Bonn-Gatchina analysis Anisovich:2011fc (), widths of about 110 MeV are also found. With -2Im =78 MeV, a smaller width was extracted in the ANL-Osaka analysis Kamano:2013iva ().

Moderate changes in our three fits can also be observed for the values of the residues and photocouplings. Although the 3/2 is well determined from elastic scattering and no new information from this channel was included in the new fits, certain changes in the resonance parameters are not surprising. Due to the well-known -wave interference in the pion-induced production resulting in a -shape form of the differential cross section (cf. Fig. 6), the 3/2 shows some sensitivity to the parameterization of the channel. As can be seen in Fig. 6 at  MeV and especially at  MeV, the description of the data differs in all three fits. The energy bin at  MeV is, on the other hand, a prime example for the systematic problems in the data. Data at backward angles are available with small error bars, but not in agreement with other data spanning the entire angular region. Underestimated normalization problems can obviously change the angular dependence significantly and have a large impact on the partial-wave content. Better data are called for.

, : Although the poles of the 5/2 and the 5/2 are located in an energy region where the prediction of fit A for and becomes worse (cf. Fig. 5), the pole positions and residues exhibit only minor differences in fit A and B. Still, as can be seen in Fig. 10, in the current fit B, the and are important to achieve a good description of the new data in eta photoproduction. Whereas for a qualitative description of the data is feasible with the , , and partial waves alone, in case of all -, -, -, and -waves are needed at medium and higher energies.

In the current form of the approach, only one bare -channel state is incorporated in the partial wave. For a discussion of a possible second explicit resonance we refer the reader to Sec. 4.3 of Ref. Ronchen:2012eg ().

Figure 10: Partial wave contribution to (upper row) and (lower row) in fit B. Solid (red) line: fit B; dashed (green) line: , , wave only; dash-dotted (indigo) line: all -, -, -, and -waves.

: The normalized residue of the 7/2 is small. Still, the inclusion of the new polarization data for the reaction results in a pole position with the real part 53 MeV and the imaginary part almost 60 MeV smaller in fit B than in fit A. The magnitudes of photocouplings, on the other hand, show much less variations in the two new fits compared to the previous fit 2. As remarked in Sec. 4.2 of Ref. Ronchen:2012eg (), from an analysis of elastic scattering not much evidence can be claimed for this resonance. However, in our current fit B the partial wave seems to play a certain role in the parameterization of in , cf. Fig. 10. At energies in the range of the pole position of the 7/2, the alone plus the -, -, and -waves does not yield a qualitative description of the data. However, evidence for this resonance from the current database is weak in general.

, , : We include bare -channel states identified with the 7/2, 9/2 and 9/2 resonance. The parameters of those broad states are less stable, as was already observed in the Jülich2012 analysis Ronchen:2012eg ().

Isospin resonances: Since the channel does not couple to resonances with , the pole positions and hadronic residues are very similar in fits A and B. The inclusion of eta photoproduction data can change the resonances only indirectly through the mixed-isospin channels in pion photoproduction, and the mixed-isospin channels in the reaction . Nonetheless, the mass of the 5/2 is about 40 MeV higher and the widths about 70 MeV larger in fit B. Smaller differences in the parameters from fits A and B can also be seen in the 7/2 and the 9/2. The photocouplings at the pole are marginally less stable than the pole positions.

Comparing fit A and fit A of Ref. Ronchen:2012eg (), the influence of pion photoproduction data included in the former fit but not in the latter is visible in the results for some states, as e.g. the 1/2. Note that also the 3/2 changes its pole position slightly.

In the analysis of pion photoproduction within the Jülich framework Ronchen:2014cna () the uncertainties of the extracted photocouplings were estimated from re-fits based on different re-weighted data sets. In Ref. Ronchen:2014cna (), all data included in the fit entered with a universal weight of one. In the present study, however, the situation is different. As described in Sec. III.1, the quality of the hadronic data requires a specific weighting of the various data sets. Moreover, in case of the elastic channel we fit to the energy-dependent partial-wave amplitudes of the GW-SAID group for which no errors are provided. As a side remark, it should be noted that the error bars of the corresponding single-energy solutions do not provide enough information for correlated fits. Furthermore, also for eta photoproduction certain data sets were included with a higher weight. This renders an error estimation as performed in Ref. Ronchen:2014cna () impracticable for the present analysis.

A comprehensive statistical error analysis is complicated by the large number of data points and free parameters, typically inherent in the kind of analysis at hand. Such an analysis is, to our knowledge, not pursued in any of the current DCC approaches and we postpone a rigorous error analysis to future work. Without such an error analysis, the assessment of the significance of certain less well-determined states, like a potential , is not possible. Concerning the significance of resonance signals, the systematic elimination of states as pursued in Ref. Arndt:2006bf () is also a necessary task which we postpone to future work.

V Summary

Over the last years, measurements of pseudoscalar meson photoproduction reactions with unprecedented quality at facilities like ELSA, MAMI, and JLab have opened a path towards a more complete picture of the baryon spectrum. The photoproduction of mesons is isospin selective and allows for an analysis of states unaffected by contributions from states. Furthermore, the final state is physically open for all resonances in the second resonance region and beyond. Eta photoproduction is, thus, a prime reaction for resonance analysis and future complete experiments. Recently, polarization observables with large angular coverage and high statistics have emerged. Among them are the target asymmetry and the beam-target asymmetry . The latter observable has been measured at MAMI for the very first time in .

However, even with the measurement of more observables and an improved coverage of the data in angles and energy, a reliable determination of resonance properties requires a combined analysis of reactions with different initial and final states. One of these reactions is where, however, the data situation is known to be problematic. To avoid bias, this requires a refit not only of parameters tied to photoproduction, but also of hadronic parameters. Baryon resonance analyses would greatly benefit from a re-measurement of the reaction.

Dynamical coupled-channel (DCC) approaches provide an especially suited tool to combine different reaction channels in a global analysis. In the present study, we extended the Jülich DCC framework to eta photoproduction. Based on a simultaneous analysis of nearly 30,000 data points for pion and eta photoproduction off the proton and the world database on the pion-induced reactions and , we extracted the spectrum of nucleon and resonances in terms of pole positions, residues and photocouplings at the pole in an energy regime from threshold up to 2.3 GeV. In the current approach, unitarity and analyticity are respected which is a prerequisite for a reliable determination of the resonance properties.

Poles and residues were compared to the preceding Jülich2012 analysis Ronchen:2012eg () in which only hadronic data were considered. The effect of the photoproduction data is most apparent for higher resonances, but also noticeable in case of well established states like the two resonances and whose widths change when photoproduction data are included. Also, some photocouplings at the pole changed in the present analysis compared to the Jülich2013 solution Ronchen:2014cna () in which only pion photoproduction data, but no eta photoproduction data were considered.

In order to estimate the influence of the recent MAMI and measurements, two different fits were performed, including the new data only in the second fit. Changes in the resonance parameters are predominantly observed for less well established states like the or higher lying resonances. Smaller but significant changes appear also for well reputed states and particularly for the photocouplings at the pole. Moreover, the new data on and have a major influence on the multipoles.

In general, the multipole content of eta photoproduction is less well established than for pion photoproduction. This calls for further measurements of single- and double polarization observables. Upcoming experiments on polarization observables will have significant impact on the resonance spectrum and will help to identify so-called missing states and determine their resonance parameters.

The authors gratefully acknowledge the computing time granted on the supercomputer JUROPA at Jülich Supercomputing Centre (JSC). This work is also partially supported by the EU Integrated Infrastructure Initiative HadronPhysics3 (contract number 283286), by the DFG (Deutsche Forschungsgemeinschaft, SFB/TR 16, “Subnuclear Structure of Matter”), by the DFG and NSFC through the Sino-German CRC 110, by the National Science Foundation, PIF grant No. PHY 1415459, and by the National Research Foundation of Korea, Grant No. NRF-2011-220-C00011.

Appendix A Renormalization of the nucleon mass and coupling

The nucleon is included as an -channel state in the partial wave. In contrast to the other explicit states in this partial wave the bare mass and the bare coupling of the nucleon are not free parameters but undergo a renormalization process such that the nucleon pole position and residue to the channel correspond to the physical values, i.e.  MeV and  Baru:2011bw (). Note that in the present study the nucleon is only allowed to couple to the channel. Effects of the coupling to other channels with significantly higher threshold energies are small and can be absorbed in the renormalization process.

In Ref. Ronchen:2012eg () the renormalization of the nucleon in the presence of two -channel states was illustrated. In the present study, we introduce an additional contact term in the partial wave. Hence, the renormalization procedure has to be modified.

For this purpose we define the following reduced self-energies


with and the bare coupling constant ,


We also define the reduced nucleon resonance vertices via


The nucleon pole position at is given by a zero of the determinant of of Eq. (LABEL:3res). Using Eqs. (16) and (18) we obtain an expression for the bare nucleon mass :

Here, we introduced the auxiliary quantities and :


In Eq. (A) all quantities , , and are evaluated at .

To determine and thus the bare nucleon coupling , we exploit that at the nucleon pole the physical residue has to agree with the residue of from Eq. (9). The physical residue is given by


where () are the bare nucleon vertices calculated at with the physical nucleon coupling instead of the bare coupling (cf. Appendix B.1. of Ref. Doring:2010ap ()),


where is the particle momentum in the center-of-mass frame. The residue of of Eq. (9) can be calculated as


with . We obtain


Both sides of Eq. (24) are evaluated at . Solving Eq. (24) we arrive at an expression for which depends only on known or fitted quantities. We can calculate the bare mass and coupling of the nucleon by inserting this expression for in Eqs. (17) and (A). Setting we recover the two resonance case, cf. Sec. 2.2 of Ref. Ronchen:2012eg ().

The renormalization procedure is performed for each step in the fitting process.

Re -2Im
[MeV] [MeV] [MeV] [deg] [%] [deg] [%] [deg] [%] [deg]
1/2 A 1497 105 23 51 110 5.6 8.9
B 1499 104 22 51 112 5.0 32 5.0
A 1498 74 17 51 120 7.7 68 15

A 1664 126 31 20 47 21 26
B 1672 137 37 21 48 20 26
A 1677 146 45 15 57 25 26

A 1349 221 32 5.0 11 4.5 1.1 140
B 1355 215 62 7.8 16 145 2.7 113
A 1353 212 59 2 2 156 1 67

A 1611 140 2.7 6.1 175 3.7 0.9
B 1651 121 3.2    55 16 12 0.4
A 1637 97 4 24 130 9.4 3.9

A 1746 317 10 0.7 0.3 127 1.8
B 1747 323 14 0.2 138 0.4 86 1.6
A 1742 318 8    161 0.5 0.8 2.2 4

A 1711 209 5.3 1.0 132 0.3 2.2 96
B 1710 219 4.2 0.7 106 1.1 0.2 79