Measurement of η meson production in γγ interactions and Γ( η → γγ ) with the KLOE detector The KLOE-2 collaboration

: We present a measurement of η meson production in photon-photon interactions produced by electron-positron beams colliding with √ s = 1 GeV. The measurement is done with the KLOE detector at the φ -factory DAΦNE with an integrated luminosity of 0.24 fb − 1 . The e + e − → e + e − η cross section is measured without detecting the out-going electron and positron, selecting the decays η → π + π − π 0 and η → π 0 π 0 π 0 . The most relevant background is due to e + e − → ηγ when the monochromatic photon escapes detection. The cross section for this process is measured as σ ( e + e − → ηγ ) = (856 ± 8 stat ± 16 syst ) pb. The combined result for the e + e − → e + e − η cross section is σ ( e + e − → e + e − η ) = (32 . 72 ± 1 . 27 stat ± 0 . 70 syst ) pb. From this we derive the partial width Γ( η → γγ ) = (520 ± 20 stat ± 13 syst ) eV. This is in agreement with the world average and is the most precise measurement to date.


Introduction
Photon-photon production of neutral mesons provides basic information on their structure. The strength of the coupling, measured by the partial decay width Γ(X → γγ), is related to the quark content of the meson and gives information on the relations between the hadronic state and its qq representation. For the light pseudoscalar mesons π 0 , η and η ′ , the coupling to real photons is measured in their γγ decays, while the coupling to space-like photons can be measured in γγ interactions. This is of particular interest in evaluating the light-by-light contribution to the anomalous magnetic moment of the muon [1]. Photonphoton interactions in electron-positron colliders were pioneered at the Frascati collider Adone in the ′ 70s [2][3][4] and since then have been used to study the production of hadrons in almost all e + e − colliders in a variety of conditions in low-and high-q 2 processes [5][6][7].
We present a measurement of the cross section e + e − → e + e − η with the KLOE detector at the φ-factory DAΦNE. The cross section σ(e + e − → e + e − η) is a convolution of the differential γγ luminosity and the γγ → η cross section. The η partial decay width Γ(η → γγ) is obtained by extrapolating the value of σ(γγ → η) for real photons.
DAΦNE is an e + e − collider designed to operate at high luminosity at the mass of the φ resonance, 1020 MeV. We analyzed data collected with DAΦNE operating off the φ peak, at √ s = 1 GeV, to reduce the large background from φ decays. The final state e + and e − are not detected, being emitted with high probability in the forward directions outside the acceptance of the detector. The production of the η meson is identified in two decay modes, η → π + π − π 0 and η → π 0 π 0 π 0 , that exploit in a complementary way the tracking system and calorimeter of the detector. The most relevant background is the radiative process e + e − → ηγ and, in both measurements, the yield of η mesons is controlled by the e + e − → ηγ cross section measured in the same data sample with a dedicated analysis. The data sample used in the analyses corresponds to an integrated luminosity of 0.24 fb −1 .

Signal and background model
For electron and positron beams colliding with energy E, the cross section for production of a state X in γγ interactions with photon 4-momenta q 1 and q 2 is where the γγ differential luminosity Φ(q 1 , q 2 ) has been calculated in [14][15][16] using different approximations and is proportional to (α/2π) 2 (ln E/m e ) 2 . For a narrow resonance of spin 0 the formation cross section is where Γ X→γγ is the radiative width, and w 2 = (q 1 + q 2 ) 2 . The transition form factor, F (q 2 1 , q 2 2 ), is equal to one for real photons and is usually parametrized in the form inspired by the Vector Dominance Model [17]. The parameter b for the η meson has been measured at high q 2 values in γγ experiments with single-tagging [18][19][20] and in the η leptonic radiative decays η → ℓ + ℓ − γ [21][22][23] at low q 2 values, closer to those of this measurement. The results do not show appreciable dependence on q 2 and the value assumed in this analysis, b η = (1.94 ± 0.15) GeV −2 , was obtained as an average of the measurements at low q 2 . The detector response for signal and background events is fully simulated with the Monte Carlo (MC) program Geanfi [24]. While Geanfi contains the event generator for all background processes, a new generator for e + e − → e + e − X events is developed and JHEP01(2013)119 interfaced to the detector simulation. Events are generated with exact matrix element according to full 3-body phase space distributions [25]. This results in the production of η mesons with non negligible transverse momentum. The relative error due to high-order radiative corrections to equation (2.1) is estimated to be 1% [26]. All background processes have been extensively studied in other analyses. A source of irreducible background is the reaction e + e − → ηγ when the monochromatic photon is emitted at small angles and is not detected. The cross section for this process is measured in the same data sample with two independent methods and the results agree with each other providing an important consistency check of the analysis. The beam-induced backgrounds were measured during data taking and background events are added to simulated events in the MC on a run-byrun basis.

The KLOE detector
The KLOE detector consists of a large volume cylindrical drift chamber, surrounded by a lead-scintillating fibers finely segmented calorimeter. A superconducting coil around the calorimeter provides a 0.52 T axial magnetic field. The beam pipe at the interaction region is spherical in shape with 10 cm radius, it is made of a Beryllium-Aluminum alloy of 0.5 mm thickness. Low-beta quadrupoles are located at ∼ 50 cm distance from the interaction region. Two small lead-scintillating tiles calorimeters (QCAL) [27] are wrapped around the quadrupoles.
The drift chamber (DC) [28], 4 m in diameter and 3.3 m long, has 12,582 drift cells arranged in 58 concentric rings with alternated stereo angles and is filled with a low-density gas mixture of 90% Helium-10% isobutane. The chamber shell is made of carbon fiberepoxy composite with an internal wall of 1.1 mm thickness at 25 cm radius. The spatial resolutions are σ xy ∼ 150 µm and σ z ∼ 2 mm. 1 The momentum resolution for long tracks is σ(p T )/p T ∼ 0.4%. Vertices are reconstructed with a spatial resolution of ∼ 3 mm.
The calorimeter [29] is divided into a barrel and two end-caps and covers 98% of the solid angle. The readout granularity is (4.4 × 4.4) cm 2 , for a total of 2440 cells arranged in five layers. Each cell is read out at both ends by photomultipliers. The energy deposits are obtained from the signal amplitude while the arrival times and the position along the fibers are obtained from the time differences. Cells close in time and space are grouped into energy clusters. The cluster energy E is the sum of the cell energies. The cluster time t and position r are energy-weighted averages. Energy and time resolutions are σ E /E = 0.057/ E (GeV) and σ t = 57 ps/ E (GeV) ⊕ 100 ps, respectively. The cluster space resolution is σ = 1.4 cm/ E (GeV) along the fibers and σ ⊥ = 1.3 cm in the orthogonal direction.
The trigger [30] uses both calorimeter and chamber information. For this analysis the events are selected by the calorimeter trigger, requiring two energy deposits with E > 50 MeV in the barrel or E > 150 MeV in the end-caps. A higher-level cosmic-ray veto rejects events with at least two energy deposits above 30 MeV in the outermost calorimeter layer. Data are then analyzed by an event classification filter [24], which selects and streams various categories of events in different output files.

JHEP01(2013)119 4 Data sample and event preselection
The data were collected at √ s = 1000.1 MeV with electron and positron beams colliding at a small angle with an average transverse momentum of 12.7 MeV in the horizontal plane. The average instantaneous luminosity was 7 × 10 31 cm −2 s −1 and the analysis is based on an integrated luminosity of 242.5 pb −1 measured with a precision of 0.3% recording large angle Bhabha scattering events [31]. Data are selected with a background rejection filter [24] before event reconstruction. A 1/20 sample of unfiltered data, corresponding to about 11 pb −1 , is also reconstructed to define the preselection filter used for the analysis and to evaluate its efficiency for event selection. The preselection filter requires • at least two energy clusters, neutral (not associated to any track) and prompt (with |t − r/c| < 5σ t ); • all prompt neutral clusters are required to have energy E γ > 15 MeV and polar angle 20 • < θ γ < 160 • ; • at least one prompt neutral cluster with energy greater than 50 MeV; • a ratio of the two highest energy neutral prompt clusters to the total calorimeter energy R = (E γ1 + E γ2 )/E tot > 0.3; • 100 MeV < E tot < 900 MeV, to reject low energy background events and the high rate processes e + e − → e + e − (γ), e + e − → γγ.

Event selection
In addition to the preselection, candidate decays η → π + π − π 0 should fulfill the following requirements • two and only two neutral prompt clusters with |t − r/c| < 3σ t and polar angle 23 • < θ γ < 157 • ; • at least two tracks with opposite curvature that are extrapolated inside a cylinder ρ = x 2 + y 2 < 8 cm and |z| < 8 cm centered around the average beam collision point; • the distance of the first DC hit to the average beam collision point to be less than 50 cm for both tracks (in case of two or more tracks with the same curvature, the track with best quality parameters is chosen); • sum of the two tracks momenta | p 1 | + | p 2 | < 700 MeV. To minimize any selection bias and to optimize the selection efficiency, there is no requirement for the tracks to be associated to clusters in the calorimeter nor that they form a vertex. The number of selected events is 3.9 × 10 6 . A small fraction of fully neutral final states can survive the two tracks requirement in case of photon conversion γN → e + e − N or π 0 Dalitz decay.
Many background contributions have been considered. The e + e − → ηγ process is a source of irreducible background when η decays to π + π − π 0 and the monochromatic photon, E γ = 350 MeV, is emitted at small polar angles and is not detected. However, the correlation of the squared missing mass, m 2 mis and the η longitudinal momentum p Lη can be used to separate the signal from the background. For the background p Lη = E γ cos θ ≃ 350 MeV and m 2 mis ≃ 0 while the signal, for small values of p T η , is characterized by Lη , as shown in figure 1. The process e + e − → ωπ 0 , with ω → π + π − π 0 , has four photons in the final state and therefore produces the same final state as the signal when two photons are not detected. The cross section has been measured with data from the same run [32], σ(e + e − → ωπ 0 → π + π − π 0 π 0 ) = (5.72 ± 0.05) nb. The e + e − → K L K S events can mimic the signal either when the K L decays to π ± ℓ ∓ ν close to the collision point and K S decays to π 0 π 0 , or when the K L escapes detection and K S → π 0 π 0 is followed by photon conversion in an e + e − pair or by a π 0 Dalitz decay. The e + e − → K + K − events can mimic the signal when both kaons decay close to the collision point, either K ± → π ± π 0 or K ± → π 0 ℓ ± ν in coincidence with K ∓ → µ ∓ ν. Also Bhabha radiative events, e + e − → e + e − γ, given the large cross section, can be a source of background in case of accidental or split clusters.

Reconstruction of
To identify the π 0 meson, clusters are paired choosing the combination that minimizes the difference between the two-cluster invariant mass and the π 0 mass. This is performed using a pseudo-χ 2 variable The energy resolution function is given in section 3, the γγ invariant mass resolution is dominated by the calorimeter energy resolution while the angle measurement gives a negligible contribution. Figure 2 shows the distribution of the χ 2 γγ variable for MC signal events and for data. In the following analysis we select events with χ 2 γγ < 8. The two tracks momenta are combined with the π 0 to identify η → π + π − π 0 decay candidates, assigning the charged pion mass to the tracks. A kinematic fit is done requiring the invariant mass of ππγγ equal to the η mass. In the fit the energies, E i , the times, t i , and the coordinates of the cluster centroid position x i , y i , z i , for the two clusters are varied. The track momenta are not varied in the minimization since they are measured with much better precision than the cluster energies. There are four constraints: the promptness of the two clusters, t i − r i /c = 0, and the mass values m γγ = m π 0 , m ππγγ = m η . Figure 3 shows the distribution of χ 2 η from the kinematic fit for MC signal events and for data. We require χ 2 η < 20 to reduce the η(→ π + π − π 0 )γ background. This process has a long tail in the χ 2 η distribution due to events with the monochromatic photon in the detector acceptance and one photon from the π 0 decay undetected, that are not rejected by the χ 2 γγ < 8 requirement. At this stage of the selection, radiative Bhabha scattering, e + e − → e + e − γ, and e + e − → γγ annihilation followed by photon conversion are still a source of background. Separation of charged pion from electron/positron tracks is done using a likelihood method when a cluster is associated to the track [33]. A cluster is associated if the distance between the centroid and the extrapolation of the track to the calorimeter wall is less than 50 cm. The e-π likelihood is based on three variables: i) the difference of time of flight, ii) the energy of the cluster; iii) the fractions of energy deposited in the first and in the fifth calorimeter layers. In this analysis events with a cluster associated to each track and a value of the likelihood estimator log L π /L e < 0 for both clusters are rejected. The background from e + e − → γγ annihilation is reduced requiring that the most energetic cluster satisfies the conditions E γ1 < 230 MeV and 27.5 • < θ γ1 < 152.5 • .
Opposite curvature track pairs can originate from split tracks. This is due to the track finding algorithm that looks for secondary vertices of kaon decays. Background of split tracks is reduced applying a topological cut based on the correlation between the tracks opening angle, α ππ , and the distance between the first DC hits associated to the two tracks by the tracking algorithm. The background from kaon decays in e + e − → K L K S and e + e − → K + K − is reduced applying a cut on the two tracks opening angle α ππ > 50 • .
Kaon decays are characterized by non prompt energy clusters, thus both the time and the energy assigned by the kinematical fit to the neutral clusters are modified by the fit constraints. This effect is observed in the time and energy pulls built with the two neutral clusters Table 1. Selection efficiency, in %, for the signal and the most relevant backgrounds. The column Selection includes the efficiency of the trigger, the background filter and the data filters described in sections 4 and 5.1. The global efficiency for e + e − γ is derived before the cut on the e-π likelihood.
where the superscript meas and fit indicate the values measured and returned by the fit, respectively. The pulls are required to satisfy χ 2 t < 7 and χ 2 E < 8. The selection efficiencies are evaluated with the MC simulation and are listed in table 1 for the signal and the most relevant background sources. The column Selection includes the efficiency of the trigger, the background filter and the data filters described in sections 4 and 5.1. The trigger efficiency is controlled by comparison of the calorimeter trigger with a complementary trigger based on the drift chamber hit patterns [30]. A sample of unfiltered data is used to control the filter efficiency.
The signal is simulated with different values of the b η parameter 2 of the form factor in equation (2.3) and the fit to derive the signal yield is repeated for each value. The values of efficiencies shown in table 1 correspond to b η = 1.94 GeV −2 .

Cross section evaluation
The analysis cuts described in section 5.2 select 2977 events. The number of signal events is derived with a 2-dimensional fit to the data. The variables used to discriminate the signal from background are the squared missing mass and the η transverse momentum in the interval −0.15 GeV 2 < m 2 mis < 0.25 GeV 2 and p T η < 300 MeV that contains 2720 events. The fit to the data is done using the simulated shapes for the signal and backgrounds, the weights are left free except for η(→ π + π − π 0 )γ whose cross section and error, measured in the same data sample (see section 8), is used as a constraint in the fit. The fit returns the fraction of data events f i = n i /n tot with the constraint i f i = 1.
The projections of the m 2 mis × p T η distribution are shown in figure 4 for the data and the backgrounds weighted by their fractions f i , and the p Lη distribution is shown in figure 5. The most relevant background is e + e − → ηγ characterized by m 2 mis ≃ 0 and p Lη ≃ ±350 MeV. Table 2 Figure 5. Distribution of the longitudinal π + π − γγ momentum. The contribution of the signal is blue, e + e − → ηγ is red, e + e − → ωπ 0 is black, e + e − → e + e − γ is green, e + e − → K + K − is light blue and e + e − → K S K L is purple.
are compared for data and MC simulation, weighted by the fractions f i returned by the fit, and good agreement is observed. The fit finds 394 ± 29 signal events.  Table 3. Systematic errors determined varying the cuts for each variable for the σ(e + e − → e + e − η → e + e − π + π − π 0 ) measurement.

JHEP01(2013)119
The contributions to the systematic error are evaluated by varying the analysis cuts by the r.m.s. width of the distributions of each variable: χ 2 γγ , χ 2 η , E γ1 , θ γ1 , α ππ , χ 2 t , χ 2 E , accounting for their correlation. This results in a systematic relative error of −2.4% to +2.6%. The contributions are listed in table 3.

Event selection
In addition to the preselection described in section 4, candidate decays η → 3π 0 should fulfill the following requirements • six and only six neutral prompt clusters with E γ > 15 MeV, |t − r/c| < 3σ t and polar angle 23 • < θ γ < 157 • ; • no tracks in the drift chamber.
The number of selected events is 9857. Many background contributions have been considered. As in the charged decay analysis, the e + e − → ηγ process is a source of irreducible background when η decays to 3π 0 and the recoil photon is not detected. The process e + e − → ωπ 0 with ω → π 0 γ produces 5 photons in the final state and is important in case of accidental or split clusters. The cross section has been measured with the same data set [32]: σ(e + e − → ωπ 0 → π 0 π 0 γ) = (0.550 ± 0.005) nb. The process e + e − → a 0 (980)γ → ηπ 0 γ can mimic the signal when η decays to 3π 0 and three photons are not detected, or it decays to γγ with split or accidental clusters. Similarly for e + e − → f 0 (980)γ → π 0 π 0 γ and e + e − → η ′ γ when η ′ decays to neutrals. Also the process e + e − → K L K S with K S → π 0 π 0 and undetected K L can mimic the signal in case of split or accidental clusters.

Reconstruction of η → 3π 0 decay
The six photons are paired choosing the combination that minimizes the difference between the γγ invariant mass of the pairs and the mass of the π 0 as described in section 5.2. In the following analysis we select events with χ 2 6γ < 14. A kinematic fit is done requiring the 6γ invariant mass to be equal to the η mass. In the fit the energies, E i , the times, t i , and the coordinates of the centroid positions x i , y i , z i , for the six clusters are varied. There are seven constraints: the promptness of the six clusters t i − r i /c = 0 and m 6γ = m η . Figure 6 shows the distribution of χ 2 η from the kinematic fit for MC signal events and for data, we require χ 2 η < 20 to reduce the background. MC simulation shows that e + e − → ηγ gives a large contribution to the tail of the distribution when the monochromatic photon is in the acceptance and is wrongly paired with a photon from η decay. In this case it also produces an enhancement at large values of the 6γ invariant mass distribution. To reduce the background we require the highest energy neutral cluster to have E γ1 < 260 MeV and the six-photon invariant mass m 6γ < 630 MeV.

Cross section evaluation
The number of signal events is derived with a 2-dimensional fit to the data. The distributions used to discriminate the signal from background are the squared missing mass and the η longitudinal momentum in the interval -0.15 GeV 2 < m 2 mis < 0.35 GeV 2 and -450 MeV < p Lη < 450 MeV that contains 2166 events. The fit to the data is done using the simulated shapes for the signal and backgrounds and the fit returns the fraction of data events f i = n i /n tot with the constraint i f i = 1.
The contribution of all backgrounds, except ηγ production, is very small, below the statistical sensitivity of the fit. The contribution of ωπ 0 derived from the value of the e + e − → ωπ 0 cross section [32] is f ωπ 0 = 0.47%. The contributions of a 0 (980)γ, η ′ γ and K L K S expected extrapolating the measurements at the φ peak are negligible.
The fit with two components gives f eeη = (33.4 ± 1.5)% and f ηγ = (66.6 ± 1.9)% using the signal efficiency evaluated with b η = 1.94 GeV −2 . The fit is repeated for all the other values. The projections of the m 2 mis × p Lη distribution are shown in figure 7 for the data and the background weighted by their relative factors f i , and the p T η distribution is shown in figure 8.
The contributions to the systematic error are evaluated by varying the analysis cuts by the r.m.s. width of the distributions of each variable: χ 2 γγ , χ 2 η , E γ1 , m 6γ , accounting for their correlation. This results in a systematic relative error of −1.5% to +2.6%. The contributions are listed in table 5.
The MC simulation statistical error of 1.0% (table 4)  610 MeV -650 MeV -1.33 +2.38 Table 5. Systematic errors determined varying the cuts for each variable for the σ(e + e − → e + e − η → e + e − 3π 0 ) measurement. Varying the cut on E γ1 gives a negligible contribution.
The two values of the cross section in equations (5.3) and (6.1) are combined accounting for the following sources of correlation: • systematic uncertainties are correlated due to the requirements on the neutral prompt clusters, the photon energy, time and position resolutions common to both selections and fit procedures; • the determination of the signal efficiencies for the two measurements that share the same transition form factor; • the systematic error in the measurement of the luminosity [31]; • the correlation between the η → π + π − π 0 and η → 3π 0 branching ratios [34].
The partial width of the η meson, Γ(η → γγ), can be determined from equations (2.1) and (2.2). The γγ differential luminosity is calculated following reference [25], the program computes also the transition form factor as parametrized in equation (2.3), for the same values of the b η parameter used in evaluating the e + e − → e + e − η cross section. Since the values of the 4-momenta q 1 and q 2 sampled in the two decay modes analyzed in sections 5 and 6 can be slightly different, the partial width is determined separately for the two decays. The theoretical error in evaluating σ(γγ → η) has been added to the systematic error due to the form factor. From the two values of the e + e − → e + e − η cross section, (5.3) and (6.1), we derive The two measurements are combined accounting for their correlations to derive Γ(η → γγ) = (520 ± 20 stat ± 13 syst ) eV . The most relevant background in the measurement of the e + e − → e + e − η cross section is due to the radiative process e + e − → ηγ. The value of the cross section has been used as a constraint in the fit in case of the η → π + π − π 0 decay while it has been derived as a by-product of the analysis of the η → 3π 0 decay. The cross section has been measured by the SND experiment [35] at VEPP-2M in the range √ s = (0.6 − 1.38) GeV, but with less precision than needed to control the analysis of e + e − → e + e − η. The cross section for e + e − → ηγ is measured exploiting the η → π + π − π 0 decay using the same data sample and the same preselection procedure described in sections 4 and 5.1 with the only difference that in this case events with three and only three neutral prompt clusters are selected. The event selection aims at finding two tracks of opposite curvature, compatible with being due to π ± , two neutral prompt clusters compatible with being originated by a π 0 decay, and a third neutral prompt cluster compatible with the photon recoiling against the π + π − π 0 system. Several background processes have been considered. e + e − → ωπ 0 with ω → π + π − π 0 is characterized by two tracks and four photons and can simulate the signal if one photon is not detected. e + e − → π + π − π 0 γ has the same configuration as the signal. e + e − → K L K S can mimic the signal when K L decays to π ± ℓ ∓ ν close to the collision point and K S decays to π 0 π 0 but one photon is not detected. e + e − → K + K − can mimic the signal when both kaons decay close to the collision point to π ± π 0 , π ∓ π 0 and one photon is not detected, or decay to π ± π 0 , µ ∓ ν and the additional photon originates from split or accidental clusters. e + e − → π + π − π 0 and e + e − → π + π − γ can mimic the signal in case of one or two accidental or split clusters. e + e − → e + e − γ has a very large cross section and can be an important background if the electron (positron) is misidentified as a pion and the two additional photons originate from split or accidental clusters. e + e − → γγ has also a large cross section and may originate background in case of photon conversions and there are split or accidental clusters. Beside these, ηγ production with η decaying to π + π − γ or to 3π 0 should be discriminated from the π + π − π 0 signal by the number of prompt neutral clusters.
8.1 Reconstruction of ηγ → π + π − π 0 γ events The identification of the π 0 meson follows the procedure described in section 5.2. No cut is applied to the value of χ 2 γγ . A kinematic fit is applied to the selected combination of three neutral prompt clusters and two tracks, with the assignment of the charged pion mass. The fit uses 15 variables, the energy E i , time t i and cluster coordinates x i , y i , z i of the three clusters, and has 7 constraints, promptness of three clusters t i − r i /c = 0, energy and momentum conservation: i E γi + E π + + E π − = √ s and i p γi + p π + + p π − = p e + e − .
The track momenta are not varied in the minimization procedure. Figure 9 shows the distribution of the χ 2 of the kinematic fit for MC signal events and for data. In the following analysis we select events with χ 2 < 50. The background of e + e − → e + e − γ and e + e − → γγ is reduced using the e-π likelihood estimator as described in section 5.2, and requiring the angle between the two tracks to be α ππ < 160 • , and the angle between any photon pair to be α γγ > 20 of e + e − → ηγ → 3π 0 γ with photon conversion is reduced requiring the sum of the photon energies i E γi < 660 MeV. At this stage of the analysis, the residual background is dominated by the processes e + e − → π + π − γ and e + e − → π + π − π 0 with split clusters, characterized by a neutral energy smaller than for the signal, and e + e − → ωγ → π + π − π 0 γ characterized by the same final state as the signal. These backgrounds are reduced by requiring for the sum of the track momenta | p + | + | p − | < 440 MeV. The effect of these cuts is controlled by the distribution of the energy of the unpaired photon shown in figure 10 where E γ3 is the value returned by the fit and has a resolution greatly improved by the good time and position resolution of the calorimeter. The peaks at the energies of the photon recoiling against the ω and the η are clearly visible over a small background at E γ3 = 194 MeV and E γ3 = 350 MeV, respectively.
The selection efficiencies are evaluated with the MC simulation described in section 2 and are listed in table 6 for the signal and the most relevant background sources.

Evaluation of the cross section
The number of signal events is derived with a 2-dimensional fit to the data. The distributions used to discriminate the signal from background are the energy of the unpaired photon and the invariant mass of the two charged pions in the interval 50 MeV < E γ3 < 400 MeV and 280 MeV < m ππ < 520 MeV that contains 55150 events. The fit to the data is done using the simulated shapes for the signal and backgrounds and the weights are left free. The projections of the E γ3 × m ππ distribution are shown in figure 11 for the data and the backgrounds weighted by their relative factors returned by the fit. The result of the fit gives 13536 ± 121 signal events resulting in a cross section σ(e + e − → ηγ → π + π − π 0 γ) = (194.7 ± 1.8 stat ) pb. Final state Selection Global efficiency (%) efficiency (%) η(→ π + π − π 0 )γ 36.71 ± 0.02 28.68 ± 0.02 π + π − π 0 γ 6.08 1.19 ωπ 0 19.80 1.07 η(→ π + π − γ)γ 0.723 0.069 η(→ neutral)γ 0.111 0.002 Table 6. Selection efficiency, in %, for the signal and the most relevant backgrounds. The column Selection includes the efficiency of the trigger, the background filter and the data filter described in section 4.
This value, obtained from a direct measurement, agrees well with the value (6.2) obtained from the analysis of γγ → η → 3π 0 . The result interpolates well with the measurements of the SND experiment [35] and has a better precision.

Summary
The cross section σ(e + e − → e + e − η) has been measured at √ s = 1 GeV with the KLOE detector based on an integrated luminosity of 0.24 fb −1 . The η mesons are selected using the two decays η → π + π − π 0 and η → π 0 π 0 π 0 that exploit in a complementary way the tracking and the calorimeter measurements. Many background processes are considered, the most relevant being e + e − → ηγ when the photon is emitted at small polar angles and escapes detection. As a consistency check, we have measured the cross section for e + e − → ηγ in two independent ways, the two values agree well with each other and we derive σ(e + e − → ηγ) = (856 ± 8 stat ± 16 syst ) pb. This value interpolates well previous measurements by the SND experiment and is more precise. The cross section for e + e − → e + e − η is obtained independently for the two η decay modes with a 2-dimensional fit to the squared missing mass and the η momentum projections. Combining the two measurements we obtain σ(e + e − → e + e − η) = (32.72 ± 1.27 stat ± 0.70 syst ) pb. This value is used to extract the partial width Γ(η → γγ) = (520 ± 20 stat ± 13 syst ) eV. This is in agreement with the world average of (510 ± 26) eV and is the most precise measurement to date.