University of Birmingham Observation of the B 0 00 decay from an amplitude analysis of B 0 ( +-) ( +-) decays LHCb Collaboration

Citation for published version (Harvard): LHCb Collaboration, Aaij, R, Adeva, B, Adinolfi, M, Affolder, A, Ajaltouni, Z, Akar, S, Albrecht, J, Alessio, F, Alexander, M, Ali, S, Alkhazov, G, Alvarez Cartelle, P, Alves, AA, Amato, S, Amerio, S, Amhis, Y, An, L, Anderlini, L, Anderson, J, Andreotti, M, Andrews, JE, Appleby, RB, Aquines Gutierrez, O, Archilli, F, Artamonov, A, Artuso, M, Aslanides, E, Auriemma, G, Baalouch, M, Bachmann, S, Back, JJ, Badalov, A, Baesso, C, Baldini, W, Barlow, RJ, Barschel, C, Barsuk, S, Barter, W, Batozskaya, V, Battista, V, Bay, A, Bifani, S, Farley, N, Griffith, P, Kenyon, IR, Lazzeroni, C, Mazurov, A, McCarthy, J, Pescatore, L & Watson, NK 2015, 'Observation of the B 0 00 decay from an amplitude analysis of B 0 ( +)( +) decays', Physics Letters, Section B: Nuclear, Elementary Particle and High-Energy Physics, vol. 747, pp. 468-478. https://doi.org/10.1016/j.physletb.2015.06.027


Introduction
The study of B meson decays to ρρ final states provides the most powerful constraint to date for the Cabibbo-Kobayashi-Maskawa (CKM) angle α ≡ arg (V td V * tb )/(V ud V * ub ) [1][2][3]. Most of the physics information is provided by the decay B 0 → ρ + ρ − as measured at the e + e − colliders at the ϒ(4S) resonance [4,5], 1 for which the dominant decay amplitude, involving the emission of a W boson only (tree), exhibits a phase difference that can be interpreted as the sum of the CKM angles β + γ = π − α in the Standard Model. The subleading amplitude associated with the exchange of a W boson and a quark (penguin) must be determined in order to interpret the electroweak phase difference in terms of the angle α. This is realised by means of an isospin analysis involving the companion modes B + → ρ + ρ 0 [6,7] and B 0 → ρ 0 ρ 0 [8,9]. 2 In particular, the smallness of the amplitude of the latter leads to a better constraint on α.
large LHCb data set may shed light on this discrepancy. In addition, LHCb may confirm the hint of B 0 → ρ 0 f 0 (980) decays reported by Belle [9]. Measurements of the B 0 → ρ 0 ρ 0 branching fraction and longitudinal polarisation fraction at LHCb can be used as inputs in the determination of α [2,3].
This work focuses on the search and study of the B 0 → (π + π − )(π + π − ) decay in which the two (π + π − ) pairs are selected in the low invariant mass range (< 1100 MeV/c 2 ). The B 0 → ρ 0 ρ 0 is expected to dominate the (π + π − ) mass spectrum. The (π + π − ) combinations can actually emerge from S-wave nonresonant and resonant contributions or other P-or D-wave resonances interfering with the signal. Hence, the determination of the B 0 → ρ 0 ρ 0 yields requires a two-body mass and angular analysis, from which the fraction of the longitudinally polarised final state can be measured.
The branching fraction is measured relative to the B 0 → φ K * (892) 0 mode. The B 0 → φ K * (892) 0 decay, which results in four light mesons in the final state, is similar to the signal, thus allowing for a cancellation of the uncertainties in the ratio of selection efficiencies.

Data sets and selection requirements
The analysed data correspond to an integrated luminosity of 1.0 fb −1 and 2.0 fb −1 from pp collisions recorded at a centre-ofmass energy of 7 TeV, collected in 2011, and 8 TeV, collected in 2012, by the LHCb experiment at CERN.
The LHCb detector [10,11] is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5 precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region [12], a large-area siliconstrip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes [13] placed downstream of the magnet.
The tracking system provides a measurement of momentum, p, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary vertex, the impact parameter, is measured with a resolution of (15 + 29/p T ) μm, where p T is the component of the momentum transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov (RICH) detectors [14]. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers [15]. The online event selection is performed by a trigger [16], which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction. In this analysis two categories of events that pass the hardware trigger stage are considered: those where the trigger decision is satisfied by the signal b-hadron decay products (TOS) and those where only the other activity in the event determines the trigger decision (TIS). The software trigger requires a two-, three-or fourtrack secondary vertex with large transverse momenta of charged particles and a significant displacement from the primary pp interaction vertices (PVs). At least one charged particle should have p T > 1.7 GeV/c and is required to be inconsistent with originating from any primary interaction. A multivariate algorithm [17] is used for the identification of secondary vertices consistent with the decay of a b hadron.
Further selection criteria are applied offline to reduce the number of background events with respect to the signal. The (π + π − ) candidates must have transverse momentum larger than 600 MeV/c, with at least one charged decay product with p T > 1000 MeV/c. The two (π + π − ) pairs are then combined to form a B 0 candidate with a good vertex quality and transverse momentum larger than 2500 MeV/c. The invariant mass of each pair of opposite-charge pions forming the B 0 candidate is required to be in the range 300-1100 MeV/c 2 . The identification of the final-state particles (PID) is performed with dedicated neural-networks-based discriminating variables that combine information from the RICH detectors and other properties of the event [14]. The combinatorial background is further suppressed with multivariate discriminators based on a boosted decision tree algorithm (BDT) [18,19]. The BDT is trained with simulated B 0 → ρ 0 ρ 0 (where ρ 0 → π + π − ) events as signal sample and candidates reconstructed with fourbody mass in excess of 5420 MeV/c 2 as background sample. The discriminating variables are based on the kinematics of the B decay candidate (B p T and the minimum p T of the two ρ 0 candidates) and on geometrical vertex measurements (quality of the B candidate vertex, impact parameter significances of the daughters, B flight distance significance and B pointing to the primary vertex). The optimal thresholds for the BDT and PID discriminating variables are determined simultaneously by means of a frequentist estimator for which no hypothesis on the signal yield is assumed [20]. The B 0 meson candidates are accepted in the mass range 5050-5500 MeV/c 2 .
The normalisation mode B 0 → φ K * (892) 0 is selected with similar criteria, requiring in addition that the invariant mass of the (K + π − ) candidate is found in a range of ±150 MeV/c 2 around the known value of the K * (892) 0 meson mass [21] and the invariant mass of the (K + K − ) pair is in a range of ±15 MeV/c 2 centred at the known value of the φ meson mass [21]. A sample enriched in B 0 → (K + π − )(π + π − ) events is selected using the same ranges in (π + π − ) and (K + π − ) masses to estimate the background with one misidentified kaon.
The presence of (π + π − ) pairs originating from J /ψ , χ c0 and χ c2 charmonia decays is vetoed by requiring the invariant masses M of all possible (π + π − ) pairs to be |M − M 0 | > 30 MeV/c 2 , where M 0 stands for the corresponding known values of the J /ψ , χ c0 and χ c2 meson masses [21]. Similarly, the decays D 0 → K − π + and D 0 → π + π − are vetoed by requiring the corresponding invariant masses to differ by 25 MeV/c 2 or more from the known D 0 meson mass [21]. To reduce contamination from other charm backgrounds and from the B 0 → a + 1 (→ ρ 0 π + )π − decay, the invariant mass of any three-body combination in the event is required to be larger than 2100 MeV/c 2 .
Simulated B 0 → ρ 0 ρ 0 and B 0 → φ K * (892) 0 decays are also used for determining the relative reconstruction efficiencies. The pp collisions are generated using Pythia [22] with a specific LHCb configuration [23]. Decays of hadronic particles are described by EvtGen [24]. The interaction of the generated particles with the detector and its response are implemented using the Geant4 toolkit [25] as described in Ref. [26].

Four-body mass fit
The four-body mass spectrum M(π + π − )(π + π − ) is fit with an unbinned extended likelihood. The fit is performed simultaneously for the two data taking periods together with the normalisation channel M(K + K − )(K + π − ) and PID misidentification control channel M(K + π − )(π + π − ) mass spectra. The four-body invariant mass models account for B 0 and possible B 0 s signals, combinatorial backgrounds, signal cross-feeds and background contributions arising from partially reconstructed b-hadron decays in which one or more particles are not reconstructed.
The B 0 and B 0 s meson shapes are modelled with a modified Crystal Ball distribution [27]. A second power-law tail is added on the high-mass side of the signal shape to account for imperfections of the tracking system. The model parameters are determined from a simultaneous fit of simulated signal events that fulfil the trigger, reconstruction and selection chain, for each data taking period. The values of the tail parameters are identical for the B 0 and B 0 s mesons. Their mass difference is constrained to the value from Ref. [21]. The mean and width of the modified Crystal Ball function are free parameters of the fit to the data.
The combinatorial background in each four-body spectrum is described by an exponential function where the slope is allowed to vary in the fit.
The misidentification of one or more final-state hadrons may result in a fully reconstructed background contribution to the corresponding signal spectrum, denoted signal cross-feed. The magnitude of the branching fractions of the signal and control modes as well as the two-body mass selection criteria make these signal cross-feeds negligible, with one exception: the misidentification of the kaon of the decay B 0 → (K + π − )(π + π − ) as a pion yields a significant contribution in the M(π + π − )(π + π − ) mass spectrum. The mass shape of Ball function, whose parameters are determined from simulated events. The yield of this signal cross-feed is allowed to vary in the fit. The measurement of the actual number of reconstructed B 0 → (K + π − )(π + π − ) events multiplied by the data-driven estimate of the misidentification efficiency is consistent with the measured yield.
The partially reconstructed background is modelled by an AR-GUS function [28] convolved with a Gaussian function accounting   Table 1 Yields from the simultaneous fit for the 2011 and 2012 data sets. The first and second uncertainties are the statistical and systematic contributions, respectively.

Decay mode
Signal yields 2011 Signal yields 2012 for resolution effects. Various mass shape parameterisations are examined. The best fit is obtained when the endpoint of the AR-GUS function is fixed to the value expected when one pion is not attributed to the decay. The other shape parameters of the ARGUS function are free parameters of the fit, common to the two data taking periods. The floating width parameter of the signal mass shape is constrained to be equal to the width of the Gaussian function used in the convolution.
spectra with the fit results overlaid. The signal event yields are shown in Table 1. Aside from the prominent signal of the B 0 → is observed with a statistical significance of more than 10 standard deviations. The statistical significance is evaluated by taking the ratio of the likelihood of the nominal fit and of the fit with the signal yield fixed to zero. A systematic uncertainty due to the fit model is associated to the measured yields. The dominant uncertainties arise from the knowledge of the signal and signal cross-feed shape parameters determined from simulated events. Several pseudoexperiments are generated while varying the shape parameters within their uncertainties, and the systematic uncertainties on the yields are estimated from the differences in results with respect to the nominal fit.

Amplitude analysis
An amplitude analysis is used to determine the vector-vector (VV) contribution B 0 → ρ 0 ρ 0 by using two-body mass spectra and angular variables. The four-body mass spectrum is first analysed with the sPlot technique [29] to subtract statistically the background under the B 0 → (π + π − )(π + π − ) signal.  [31,32], and a low-mass component. The latter includes the broad low-mass resonance f 0 (500) and a non-resonant contributions, which are jointly modelled in the framework of the K -matrix formalism [33] and referred as M (ππ) 0 (m i ). Following the K -matrix formalism, the amplitude for the low-mass π + π − S-wave can be written as where κ is measured to be −0.07 ± 0.24 from a fit to the inclusive π + π − mass distribution and m 0 and are the nominal mass and mass-dependent width of the f 0 (500), as determined in Ref. [34]. The functions ρ(m) and q(m), defined in Ref. [33], are the phase space factor and the relative momentum of a pion in the ρ 0 centre-of-mass system. By convention, the phase of the M (ππ) 0 (m i ) mass propagator is set to zero at the ρ 0 nominal mass. as defined in Ref. [35]. The contribution of the decays The differential decay rate for B 0 → (π + π − )(π + π − ) decays at the B 0 production time t = 0 is given by where the variables θ 1 , θ 2 and ϕ are the helicity angles, described in Fig. 2, and 4 is the four-body phase space factor. The notations of the complex amplitudes, A i , and the expressions of their related angular distributions, f i , are displayed in Table 2. The mass propagators included in the f i functions are normalised to unity in the fit range.
For the CP conjugated mode, B 0 → (π + π − )(π + π − ), the decay rate is obtained under the transformation A i → η i A i , where η i is the CP eigenvalue of the CP eigenstate i, shown in Table 2.
The untagged time-integrated decay rate of B 0 and B 0 to four pions, assuming no CP violation, can be written as d 5 ( + ) d cos θ 1 d cos θ 2 dϕ dm 2 (5) where δ ij = 1 when i = j and δ ij = 0 otherwise.
The efficiency of the selection of the final state B 0 → (π + π − )(π + π − ) varies as a function of the helicity angles and the two-body invariant masses. To take into account variations in the efficiencies, four event categories k are defined according to their hardware trigger decisions (TIS or TOS) and data taking period (2011 and 2012).
The acceptance is accounted for through the complex integrals where f i are the functions given in Table 2 and the overall efficiency. The integrals are computed with simulated events of each of the four considered categories, selected with the same criteria as those applied to data, following the method described in Ref. [38].
The coefficients ω k ij are used to determine the efficiency and to build a probability density function for each category, which is defined as The four event categories are used in the simultaneous unbinned maximum likelihood fit which depends on the 19 free parameters indicated in Table 3.
Systematic effects are estimated by fitting with the angular model and ensemble of 1000 pseudoexperiments generated with the same number of events as observed in data. The biases are for the parameters of interest consistent with zero. A systematic uncertainty is assigned by taking 50% of the fit bias or the uncertainty on the rms when the latter is bigger in order to account for possible statistical fluctuations.
Several model related uncertainties are envisaged. The B 0 → a ± 1 π ∓ angular model requires knowledge of the lineshape of the a ± 1 meson. The a ± 1 natural width is chosen to be 400 MeV/c 2 . The difference to the fit results obtained by varying the width from 250 to 600 MeV/c 2 is taken as the corresponding systematic uncertainty. In addition, a systematic uncertainty is obtained by introducing the CP-odd component in the fit model of the decay amplitude B 0 → a ± 1 π ∓ by fixing the relative amplitudes of B 0 → a + 1 π − and B 0 → a − 1 π + components to the values measured in Ref. [39]. Table 2 Amplitudes, A i , CP eigenvalues, η i , and mass-angle distributions, f i , of the B 0 → (π + π − )(π + π − ) model. The indices ijkl indicate the eight possible combinations of pairs of opposite-charge pions. The angles α kl , β ij and kl are defined in Ref. [37].

Table 3
Results of the unbinned maximum likelihood fit to the angular and two-body invariant mass distributions. The first uncertainty is statistical, the second systematic.
Parameter Definition Fit result 0.29 +0.12 −0.09 ± 0.08 Another source of uncertainty originates in the modelling of the low mass (π + π − ) S-wave lineshape. The f 0 (500) mass and natural width uncertainties from Ref. [34] and the uncertainty on the parameter that quantifies the non-resonant contribution are propagated to the angular analysis parameters by generating and fitting 1000 pseudoexperiments in which these input values are varied according to a Gaussian distribution having their uncertainties as widths. The root mean square of the distribution of the results is assigned as a systematic uncertainty. The same strategy is followed to estimate the systematic uncertainties originating from the ρ 0 , f 0 (500) and ω lineshape parameters.
The uncertainty related to the background subtraction method is estimated by varying within their uncertainties the fixed parameters of the mass fit model and studying the resulting angular distributions and two-body mass spectra. The difference to the fit results is taken as a systematic uncertainty. An alternative subtraction of the background estimated from the high-mass sideband is performed, yielding compatible results.
The knowledge of the acceptance model described in Eq. (6) comes from a finite sample of simulated events. An ensemble of pseudoexperiments is generated by varying the acceptance weights according to their covariance matrix. The root mean square of the distribution of the results is assigned as a systematic uncertainty.
The resolution on the helicity angles is evaluated with pseudoexperiments resulting in a negligible systematic uncertainty. The systematic uncertainty related to the (π + π − ) mass resolution is estimated with pseudoexperiments by introducing a smearing of the (π + π − ) mass. Differences in the parameters between the fit with and without smearing are taken as a systematic uncertainty. Table 4 details the contributions to the systematic uncertainty in the measurement of the fraction of B 0 → ρ 0 ρ 0 signal decays in the B 0 → (π + π − )(π + π − ) and its longitudinal polarisation fraction.
The final results of the combined two-body mass and angular analysis are shown in Fig. 3 and Table 3. The fit also allows for Table 4 Relative systematic uncertainties on the longitudinal polarisation parameter, f L , and the fraction of B 0 → ρ 0 ρ 0 decays in the B 0 → (π + π − )(π + π − ) sample. The model uncertainty includes the three uncertainties below.

Systematic effect
Uncertainty on f L (%) Uncertainty on the extraction of the fraction of B 0 → ρ 0 ρ 0 decays in the B 0 → (π + π − )(π + π − ) sample, defined as  and B 0 → a ± 1 π ∓ (light blue) are also displayed. The B 0 → ρ 0 ρ 0 contribution is split into longitudinal (dashed red) and transverse (dotted red) components. Interference contributions are only plotted for the total (black) model. The efficiency for longitudinally polarised B 0 → ρ 0 ρ 0 events is ∼5 times smaller than for the transverse component.
(For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)
The product λ f L · P (B 0 → ρ 0 ρ 0 ) is determined from the amplitude analysis to be 1.13 ± 0.19 (stat) ± 0.10 (syst). This quantity is mainly related to the modelling of the S-wave component, and dominates the systematic uncertainty of the parameters of interest. [40]. A 1% systematic uncertainty is added, accounting for differences in the selection acceptance for P-and S-wave contributions.
The amounts of B 0 → (π + π − )(π + π − ) and B 0 → (K + K − )(K + π − ) candidates are determined from the four-body mass spectra analysis and their associated statistical and systematical uncertainties are propagated quadratically to the branching fraction uncertainty estimate.
The limited size of the simulated events samples that meet all selection criteria result in a systematic uncertainty of 1.7% (2.6%) on the measurement of the relative branching fraction for the 2011 (2012) data-taking period. The impact of the discrepancies between experimental and simulated data related to the B 0 meson kinematical properties is 0.6% (1.2%). The efficiencies of the particle-identification requirements are determined from control samples of data with a systematic uncertainty of 0.5%, mostly originating from the limited size of the calibration samples. An additional 1% systematic uncertainty on the tracking efficiency is added accounting for different interaction lengths between π and K .
The relative branching fraction is measured to be = 0.094 ± 0.017 (stat) ± 0.009 (syst). (10) The agreement between the results obtained in the two datataking periods is tested with the best linear estimator technique [41] yielding compatible results.

Conclusions
The full data set collected by the LHCb experiment in 2011 and 2012, corresponding to an integrated luminosity of 3.0 fb −1 , is analysed to search for the B 0 → ρ 0 ρ 0 decay. A yield of 634 ± 28 ± 8 B 0 → (π + π − )(π + π − ) signal decays with π + π − pairs in the 300-1100 MeV/c 2 mass range is obtained. An amplitude analysis is conducted to determine the contribution from B 0 → ρ 0 ρ 0 decays. This decay mode is observed for the first time with a significance of 7.1 standard deviations. In the same π + π − pairs mass range, B 0 s → (π + π − )(π + π − ) decays are also observed with a statistical significance of more than 10 standard deviations. where the last uncertainty is due to the normalisation channel. These results are the most precise to date and will improve the precision of the determination of the CKM angle α.
The measured longitudinal polarisation fraction is consistent with the measured value from BaBar [8] while it differs by 2.3 standard deviations from the value obtained by Belle [9]. The branching fraction measurement is in agreement with the values measured by both BaBar [8] and Belle [9] Collaborations.