Photoproduction of $\omega$ mesons on nuclei near the production threshold

The photoproduction of $\omega$ mesons on LH_2, C and Nb has been studied for incident photon energies from 900 to 1300 MeV using the CB/TAPS detector at ELSA. The $\omega$ line shape does not show any significant difference between the LH_2 and the Nb targets. The data are compared with results of transport calculations that predict a sensitivity of the $\omega$ line shape to in-medium modifications near the production threshold on a free nucleon of $E_{\gamma}^{lab}$ = 1109 MeV.}


Introduction
Modifications of hadron properties in a strongly interacting environment have attracted a lot of attention and have been intensively studied both theoretically and experimentally. These studies were motivated by the expec-Correspondence to: Mariana.Nanova@physik.uni-giessen.de tation that chiral symmetry, a fundamental symmetry of Quantum Chromodynamics (QCD) in the limit of vanishing quark masses, may be restored in a nuclear medium at high temperatures or densities. In vacuum this symmetry is broken as visible in the low mass part of the hadronic spectrum: chiral partners -hadronic states with the same spin but opposite parity -like the ρ and a 1 meson are  different in mass while they should be mass degenerate if chiral symmetry were to hold. It turns out, however, (for details see [1]) that a connection between chiral symmetry restoration and hadronic in-medium spectral functions is much more involved. QCD sum rules provide a link between the quark-gluon sector and hadronic descriptions but do not fix the properties of hadrons in the strongly interacting medium. They only provide constraints for hadronic models which are still needed for calculating the in-medium self-energies of hadrons and their spectral functions.
Experimentally this field has been addressed in heavy ion collisions and reactions with proton and photon beams. Light vector mesons are particularly suited for these investigations since -after production in a nuclear reactionthey decay in the nuclear medium with sizable probability because of their short lifetimes. Experimental results are summarized and critically evaluated in recent reviews [1,2]. Almost all experiments report a softening of the spectral functions of the light vector mesons ρ, ω and φ. Increases in width are observed depending on the density and temperature of the hadronic environment. Mass shifts are only reported by the KEK group [3,4] who studied ρ, ω and φ production in proton nucleus reactions at 12 GeV. The claim of a mass shift of the ω meson in photoproduction on N b [5] has not been confirmed in a re-analysis of the data [6].
In the latter experiment incident photon energies covered the range from 900 -2200 MeV. Because of the increase of the production cross section with photon energy most of the observed ω mesons are produced with photons of energy larger than 1500 MeV. For the energy range of 1500 -2200 MeV transport calculations [7,8] have shown that the ω lineshape is rather insensitive to different inmedium modification scenarios like broadening or broad-ening and mass shift since most of the ω decays occur outside of the nuclear medium, even despite of a cut on the ω momentum (p ω ≤ 500 MeV/c). Furthermore, due to the experimentally observed strong absorption of ω mesons in the nuclear medium [9] ω mesons produced in the interior of the nucleus are largely removed by inelastic reactions and do not reach the detector; information on possible in-medium mass shifts thereby gets lost. The limited sensitivity of the ω lineshape to in-medium modifications has been confirmed experimentally in [6].
Gallmeister et al. [7] argue that a search for medium effects would be much more promising for incident photon energies below or near the photoproduction threshold on a free nucleon of E lab γ = 1109 MeV. New calculations along these lines illustrate in Fig.1 the expected sensitivity of the ω signal to various in-medium changes, such as mass shift with and without collisional broadening for two different energy ranges. It is seen that the lower-energy window indeed leads to a more pronounced -though not dramatic -sensitivity than the higher-energy one. This relatively weak sensitivity is to a large extent simply a consequence of the density profile of the nucleus that spans all densities from 0 to ρ 0 and thus smears any density-dependent signal. Assuming for the sake of the argument a density profile with a constant density of 0.6ρ 0 -roughly corresponding to the average density in nuclei -and a sharp fall off at the surface the dash-dotted line in Fig.1 is obtained; here the in-medium signal is significantly stronger. For a realistic nuclear density profile contributions to the spectral function from the surface dominate, suppressing contributions from higher density regions [10].
In both energy windows a tail towards lower masses is predicted for the scenario of a dropping ω mass. This tail is due to ω mesons which are produced off-shell within the nucleus. In [7,8]  mass tail region was obtained. This calculation used a phenomenological method for the off-shell propagation while the present results are based on the theoretical framework provided by Leupold [11] and Juchem and Cassing [12] who have derived equations of motion for the testparticles that represent the spectral function in the transport simulation starting from the general Kadanoff-Baym equations [13]. While the Leupold and Juchem-Cassing method has a firm theoretical basis it is nevertheless useful to remember that it relies on the assumption of adiabaticity, i.e. small gradients in space and time of all physical proper-ties (potential, density, spectral function). Non-adiabatic effects such as those at level-crossings are not taken into account.

Experimental set up
The experiment was performed at the ELSA electron accelerator facility [14,15] at the University of Bonn, using the Crystal Barrel/TAPS detector set up which provides an almost complete coverage of the full solid angle. The features of this calorimeter system and its capability for the detection of multi-photon final states have been described in detail elsewhere [16,17,18]. Tagged photon beams of 900 -1300 MeV were generated by bremsstrahlung and impinged in subsequent runs on LH 2 , C and N b targets with thicknesses of 53 mm, 20 mm and 1 mm (30 mm diameter), respectively. The running conditions were the same as described in [6].

Analysis
The analysis of the data follows exactly the scheme described in detail in [6]. In fact, the data discussed here represent a subset of the data published in [6]. ω mesons were reconstructed in the reaction γA → (A − 1)pω → (A − 1)pπ 0 γ → (A − 1)pγγγ from events with 3 photons and 1 proton in the final state. Fig. 2a shows the π 0 γ invariant mass spectrum for the carbon target for events with 1 charged and 3 neutral hits in the system. An ω signal is observed on a steeply falling background. As discussed in [6] this background stems dominantly from π 0 π 0 and also π 0 η → 4γ final states where due to shower overlap in the detector or detector inefficiencies 1 out of the 4 photons escapes detection. Fig. 2a also shows the π 0 γ invariant mass spectrum for events with 1 charged and 4 neutral hits where 1 neutral hit has been arbitrarily omitted to simulate the background in the 3 neutral and 1 charged spectrum. The slopes of the two spectra are different, reflecting the energy dependence of the probability to register only 3 out of 4 photons relative to detecting all 4 photons. The ratio of the two spectra in Fig. 2a is given in Fig. 2b. This ratio is a smooth function of the invariant mass as verified by Monte Carlo simulations, also shown in Fig. 2b. Since no strong in-medium effects are expected for a light nucleus like carbon this ratio can be applied to correct the π 0 γ background spectrum derived from events with 1 charged and 4 neutral hits measured for the N b target. Fig. 2c shows the π 0 γ invariant mass spectrum for N b together with this background distribution after applying this correction. The normalization of the background spectrum is done by requesting the same number of events in the mass range from 400 -900 MeV, excluding the counts in the ω peak which account for only 2 % of the total yield in this mass range. It is important to note that this determination of the background in magnitude and shape does not pay any attention to the ω signal region. The background in the (3 neutral +1 charged) signal spectrum is well reproduced by the corrected (4 neutral + 1 charged) spectrum. This is demonstrated by the ratio of the two spectra shown in Fig. 2d. In the mass range from 400 -700 MeV the average deviation from unity is 5 %. This uncertainty determines the systematic error in the background subtraction.

Comparison to the LH 2 reference measurement and to GiBUU simulations
After subtraction of the background spectrum from the π 0 γ signal spectrum the ω lineshape shown in Fig. 3a is obtained for the N b target. The experimental distribution has been fitted using the Novosibirsk function [19]: where Here A is the amplitude of the signal, x 0 is the peak position, σ is FWHM/2.35 and τ is the asymmetry parameter. This function takes into account the tail in the region of lower invariant masses resulting from the energy response of the calorimeters. The fit in the mass range 700 to 820 MeV is compared with the ω signal measured on the LH 2 target and with a Monte Carlo simulation of the ω signal in Fig. 3a. It should be noted that the 53 mm length of the LH 2 target leads to an increase in ω-signal width by less than 10 %. The agreement between the ω signal on the LH 2 target with the Monte Carlo simulation demonstrates that the detector response is under control. The fit parameters are σ = (29.1 ± 2.8) MeV, τ = -0.29 ± 0.16 for N b and σ = (28.1 ± 0.6) MeV, τ = -0.09 ± 0.02 for LH 2 , respectively. The fit of the data with the function of eqs.(1),(2) yields a χ 2 /DoF = 18.95/12 with a χ 2 -probability of 9.0% while a fit of the data with the ω line shape measured on the LH 2 target gives a χ 2 /DoF = 32.8/15, corresponding to a χ 2 -probability of 0.5%. Nevertheless, in view of the systematic and statistical uncertainties no significant deviation from the reference signals is claimed. Higher statistics will be needed to draw any conclusion. Corresponding data have been taken at MAMI-C using the Crystal Ball/TAPS set up. The analysis is ongoing.
In Fig. 3b the measured ω signal is compared to predictions of transport calculations using the GiBUU model [20] for the same scenarios as in Fig. 1. While all the curves seem to underestimate the data slightly on the low mass side of the ω peak, the experimental data obviously do not allow to distinguish between the various theoretical scenarios, in contrast to initial expectations.

Conclusions
The photoproduction of ω mesons on LH 2 , C and N b targets has been measured for incident photon energies from 900 -1300 MeV, i.e. near the photoproduction threshold on a free nucleon of E lab γ = 1109 MeV. The experimentally observed ω signal does not allow to distinguish between various in-medium scenarios which -also near thresholdlead all to only a weak tail at low invariant masses. Access to the in-medium spectral function of vector mesons is thus very limited, mainly due to the dependence of the in-medium properties such as mass and width on the nuclear density [10,21] and the inherent density-smearing caused by the density profile of nuclei. On the other hand, transparency measurements [22,23] can give at least access to the imaginary part of the in-medium self energy of the hadron. Another promising tool could be the measurement of excitation functions [24]. Such experiments are presently being analyzed.