University of Birmingham Observation of photon polarization in the b s transition

This Letter presents a study of the flavor-changing neutral current radiative B (cid:1) → K (cid:1) π ∓ π (cid:1) γ decays performed using data collected in proton-proton collisions with the LHCb detector at 7 and 8 TeV center-of-mass energies. In this sample, corresponding to an integrated luminosity of 3 fb − 1 , nearly 14 000 signal events are reconstructed and selected, containing all possible intermediate resonances with a K (cid:1) π ∓ π (cid:1) final state in the ½ 1 . 1 ; 1 . 9 (cid:2) GeV =c 2 mass range. The distribution of the angle of the photon direction with respect to the plane defined by the final-state hadrons in their rest frame is studied in intervals of K (cid:1) π ∓ π (cid:1) mass and the asymmetry between the number of signal events found on each side of the plane is obtained. The first direct observation of the photon polarization in the b → s γ transition is reported with a significance of 5 . 2 σ .

Observation of Photon Polarization in the b → sγ Transition R. Aaij et al. * (LHCb Collaboration) (Received 27 February 2014;published 22 April 2014) This Letter presents a study of the flavor-changing neutral current radiative B AE → K AE π ∓ π AE γ decays performed using data collected in proton-proton collisions with the LHCb detector at 7 and 8 TeV centerof-mass energies. In this sample, corresponding to an integrated luminosity of 3 fb −1 , nearly 14 000 signal events are reconstructed and selected, containing all possible intermediate resonances with a K AE π ∓ π AE final state in the ½1.1; 1.9 GeV=c 2 mass range. The distribution of the angle of the photon direction with respect to the plane defined by the final-state hadrons in their rest frame is studied in intervals of K AE π ∓ π AE mass and the asymmetry between the number of signal events found on each side of the plane is obtained. The first direct observation of the photon polarization in the b → sγ transition is reported with a significance of 5.2σ. DOI: 10.1103/PhysRevLett.112.161801 PACS numbers: 13.20.He, 12.15.Mm, 14.40.Nd The standard model (SM) predicts that the photon emitted from the electroweak penguin loop in b → sγ transitions is predominantly left-handed, since the recoiling s quark that couples to a W boson is left-handed. This implies maximal parity violation up to small corrections of the order m s =m b . While the measured inclusive b → sγ rate [1] agrees with the SM calculations, no direct evidence exists for a nonzero photon polarization in this type of decay. Several extensions of the SM [2], compatible with all current measurements, predict that the photon acquires a significant right-handed component, in particular, due to the exchange of a heavy fermion in the penguin loop [3].
This Letter presents a study of the radiative decay B þ → K þ π − π þ γ , previously observed at the B factories with a measured branching fraction of ð27.6 AE 2.2Þ × 10 −6 [1,4,5]. The inclusion of charge-conjugate processes is implied throughout. Information about the photon polarization is obtained from the angular distribution of the photon direction with respect to the normal to the plane defined by the momenta of the three final-state hadrons in their center-of-mass frame. The shape of this distribution, including the up-down asymmetry between the number of events with the photon on either side of the plane, is determined. This investigation is conceptually similar to the historical experiment that discovered parity violation by measuring the up-down asymmetry of the direction of a particle emitted in a weak decay with respect to an axial vector [6]. In B þ → K þ π − π þ γ decays, the up-down asymmetry is proportional to the photon polarization λ γ [7,8] and therefore measuring a value different from zero corresponds to demonstrating that the photon is polarized. The currently limited knowledge of the structure of the K þ π − π þ mass spectrum, which includes interfering kaon resonances, prevents the translation of a measured asymmetry into an actual value for λ γ .
The differential B þ → K þ π − π þ γ decay rate can be described in terms of θ, defined in the rest frame of the final state hadrons as the angle between the direction opposite to the photon momentump γ and the normal p π;slow ×p π;fast to the K þ π − π þ plane, wherep π;slow and p π;fast correspond to the momenta of the lower and higher momentum pions, respectively. Following the notation and developments of Ref. [7], the differential decay rate of B þ → K þ π − π þ γ can be written as a fourth-order polynomial in cos θ dΓ dsds 13  where s ij ¼ ðp i þ p j Þ 2 and s ¼ ðp 1 þ p 2 þ p 3 Þ 2 , and p 1 , p 2 , and p 3 are the four-momenta of the π − , π þ , and K þ mesons, respectively. The functions a k depend on the resonances present in the K þ π − π þ mass range of interest and their interferences. The up-down asymmetry is defined as which is proportional to λ γ . The LHCb detector [9] is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, designed for the study of particles containing b or c quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon-strip 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 placed downstream. The combined tracking system provides a momentum measurement with relative uncertainty that varies from 0.4% at 5 GeV=c to 0.6% at 100 GeV=c, and impact parameter resolution of 20 μm for tracks with a few GeV=c of transverse momentum (p T ). Different types of charged hadrons are distinguished by information from two ring-imaging Cherenkov detectors. Photon, electron, and hadron candidates are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter, and a hadronic calorimeter. The trigger 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.
Samples of simulated events are used to understand signal and backgrounds. In the simulation, pp collisions are generated using PYTHIA [10] with a specific LHCb configuration [11]. Decays of hadronic particles are described by EVTGEN [12], in which final state radiation is generated using PHOTOS [13]. The interaction of the generated particles with the detector and its response are implemented using the GEANT4 toolkit [14] as described in Ref. [15].
This analysis is based on the LHCb data sample collected from pp collisions at 7 and 8 TeV center-of-mass energies in 2011 and 2012, respectively, corresponding to 3 fb −1 of integrated luminosity. The B þ → K þ π − π þ γ candidates are built from a photon candidate and a hadronic system reconstructed from a kaon and two oppositely charged pions satisfying particle identification requirements. Each of the hadrons is required to have a minimum p T of 0.5 GeV=c and at least one of them needs to have a p T larger than 1.2 GeV=c. The isolation of the K þ π − π þ vertex from other tracks in the event is ensured by requiring that the χ 2 of the three-track vertex fit and the χ 2 of all possible vertices that can be obtained by adding an extra track differ by more than 2. The K þ π − π þ mass is required to be in the ½1.1; 1.9 GeV=c 2 range. High transverse energy (> 3.0 GeV) photon candidates are constructed from energy depositions in the electromagnetic calorimeter. The absence of tracks pointing to the calorimeter is used to distinguish neutral from charged electromagnetic particles. A multivariate algorithm based on the energy cluster shape parameters is used to reject approximately 65% of the π 0 → γγ background in which the two photons are reconstructed as a single cluster, while keeping about 95% of the signal photons. The B þ candidate mass is required to be in the ½4.3; 6.9 GeV=c 2 range. Backgrounds that are expected to peak in this mass range are suppressed by removing all candidates consistent with aD 0 → K þ π − π 0 or ρ þ → π þ π 0 decay when the photon candidate is assumed to be a π 0 .
A boosted decision tree [16,17] is used to further improve the separation between signal and background. Its training is based on the following variables: the impact parameter χ 2 of the B þ meson and of each of the final state hadrons, defined as the difference between the χ 2 of a primary vertex (PV) reconstructed with and without the considered particle; the cosine of the angle between the reconstructed B þ momentum and the vector pointing from the PV to the B þ decay vertex; the flight distance of the B þ meson, and the K þ π − π þ vertex χ 2 .
The mass distribution of the selected B þ → K þ π − π þ γ signal is modeled with a double-tailed Crystal Ball [18] probability density function (PDF), with power-law tails above and below the B mass. The four tail parameters are fixed from simulation; the width of the signal peak is fit separately for the 2011 and 2012 data to account for differences in calorimeter calibration. Combinatorial and partially reconstructed backgrounds are considered in the fit, the former modeled with an exponential PDF, the latter described using an ARGUS PDF [19] convolved with a Gaussian function with the same width as the signal to account for the photon energy resolution. The contribution to the partially reconstructed background from events with only one missing pion is considered separately.
The fit of the mass distribution of the selected B þ → K þ π − π þ γ candidates (Fig. 1) returns a total signal yield of 13876 AE 153 events, the largest sample recorded for this channel to date. Figure 2 shows the background-subtracted K þ π − π þ mass spectrum determined using the technique of Ref. [20], after constraining the B mass to its nominal value. No peak other than that of the K 1 ð1270Þ þ resonance can be clearly identified. Many kaon resonances, with various masses, spins, and angular momenta, are expected to contribute and interfere in the considered mass range [1].
The contributions from single resonances cannot be isolated because of the complicated structure of the K þ π − π þ mass spectrum. The up-down asymmetry is thus studied inclusively in four intervals of K þ π − π þ mass. The ½1.4; 1.6 GeV=c 2 interval, studied in Ref. [7], includes the K 1 ð1400Þ þ , K Ã 2 ð1430Þ þ and K Ã ð1410Þ þ resonances with small contributions from the upper tail of the K 1 ð1270Þ þ . At the time of the writing of Ref. [7], the K 1 ð1400Þ þ was believed to be the dominant 1 þ resonance, so the K 1 ð1270Þ þ was not considered. However, subsequent experimental results [21] demonstrated that the K 1 ð1270Þ þ is more prominent than the K 1 ð1400Þ þ ; hence, the ½1.1; 1.3 GeV=c 2 interval is also studied here. The ½1.3; 1.4 GeV=c 2 mass interval, which contains the overlap region between the two K 1 resonances, and the ½1.6; 1.9 GeV=c 2 high mass interval, which includes spin-2 and spin-3 resonances, are also considered.
In each of the four K þ π − π þ mass intervals, a simultaneous fit to the B -candidate mass spectra in bins of the photon angle is performed in order to determine the background-subtracted angular distribution; the previously described PDF is used to model the mass spectrum in each bin, with all of the fit parameters being shared except for the yields. Since the sign of the photon polarization depends on the sign of the electric charge of the B candidate, the angular variable cosθ ≡ chargeðBÞ cos θ is used. The resulting background-subtracted cosθ distribution, corrected for the selection acceptance and normalized to the inverse of the bin width, is fit with a fourth-order polynomial function normalized to unit area, where L i ðxÞ is the Legendre polynomial of order i and c i is the corresponding coefficient. Using Eqs. (1) and (3) the up-down asymmetry defined in Eq. (2) can be expressed as As a cross-check, the up-down asymmetry in each mass interval is also determined with a counting method, rather than an angular fit, as well as considering separately the B þ and B − candidates. All these checks yield compatible results.
The results obtained from a χ 2 fit of the normalized binned angular distribution, performed taking into account the full covariance matrix of the bin contents and all of the systematic uncertainties, are summarized in Table I. These systematic uncertainties account for the effect of choosing a different fit model, the impact of the limited size of the simulated samples on the fixed parameters, and the possibility of some events migrating from a bin to its neighbor because of the detector resolution, which gives the dominant contribution. The systematic uncertainty associated with the fit model is determined by performing the mass fit using several alternative PDFs, while the other two are estimated with simulated pseudoexperiments. Such uncertainties, despite being of the same size as the statistical uncertainty, do not substantially affect the fit results since they are strongly correlated across all angular bins.
The fitted distributions in the four K þ π − π þ mass intervals of interest are shown in Fig. 3. In order to illustrate the effect of the up-down asymmetry, the results of another fit imposing c 1 ¼ c 3 ¼ 0, hence forbidding the terms that carry the λ γ dependence, are overlaid for comparison.
The combined significance of the observed up-down asymmetries is determined from a χ 2 test where the null hypothesis is defined as λ γ ¼ 0, implying that the up-down asymmetry is expected to be zero in each mass interval. The corresponding χ 2 distribution has 4 degrees of freedom, and the observed value corresponds to a p value of 1.7 × 10 −7 . This translates into a 5.2σ significance for nonzero up-down asymmetry. Up-down asymmetries can be computed also for an alternative definition of the photon Background-subtracted K þ π − π þ mass distribution of the B þ → K þ π − π þ γ signal. The four intervals of interest, separated by dashed lines, are shown.
TABLE I. Legendre coefficients obtained from fits to the normalized background-subtracted cosθ distribution in the four K þ π − π þ mass intervals of interest. The up-down asymmetries are obtained from Eq. (4). The quoted uncertainties contain statistical and systematic contributions. The K þ π − π þ mass ranges are indicated in GeV=c 2 and all the parameters are expressed in units of 10 −2 . The covariance matrices are given in Ref. [22]. angle, obtained using the normalp π − ×p π þ , instead of p π;slow ×p π;fast . The obtained values, along with the relative fit coefficients, are listed in Table II. To summarize, a study of the inclusive flavor-changing neutral current radiative B þ → K þ π − π þ γ decay, with the K þ π − π þ mass in the ½1.1; 1.9 GeV=c 2 range, is performed on a data sample corresponding to an integrated luminosity of 3 fb −1 collected in pp collisions at 7 and 8 TeV centerof-mass energies by the LHCb detector. A total of 13876 AE 153 signal events is observed. The shape of the angular distribution of the photon with respect to the plane defined by the three final-state hadrons in their rest frame is determined in four intervals of interest in the K þ π − π þ mass spectrum. The up-down asymmetry, which is proportional to the photon polarization, is measured for the first time for each of these K þ π − π þ mass intervals. The first observation of a parity-violating photon polarization different from zero at the 5.2σ significance level in b → sγ transitions is reported. The shape of the photon angular distribution in each bin, as well as the values for the up-down asymmetry, may be used, if theoretical predictions become available, to determine for the first time a value for the photon polarization, and thus constrain the effects of physics beyond the SM in the b → sγ sector.
We express our gratitude to our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We thank the technical and administrative staff at the LHCb institutes. We acknowledge support from CERN and from the national agencies: CAPES, CNPq, FAPERJ, and FINEP (Brazil); NSFC (China); CNRS/IN2P3 and Region Auvergne (France); BMBF, DFG, HGF, and MPG (Germany); SFI (Ireland); INFN (Italy); FOM and NWO (The Netherlands); SCSR (Poland); MEN/IFA (Romania); MinES, Rosatom, RFBR, and NRC "Kurchatov Institute" (Russia); MinECo, XuntaGal, and GENCAT (Spain); SNSF and SER (Switzerland); NAS Ukraine (Ukraine); STFC (United Kingdom); NSF (USA). We also acknowledge the support received from the ERC under FP7. The Tier1 computing centres are supported by IN2P3 (France), KIT and BMBF (Germany), INFN (Italy), NWO and SURF (Netherlands), PIC (Spain), GridPP (United Kingdom). We are indebted towards the communities behind the multiple open source software packages we depend on. We are also thankful for the computing resources and the access to software R&D tools provided by Yandex LLC (Russia). TABLE II. Legendre coefficients obtained from fits to the normalized background-subtracted cosθ distribution, using the alternative normalp π − ×p π þ for defining the direction of the photon, in the four K þ π − π þ mass intervals of interest. The up-down asymmetries are obtained from Eq. (4). The quoted uncertainties contain statistical and systematic contributions. The K þ π − π þ mass ranges are indicated in GeV=c 2 and all the parameters are expressed in units of 10 −2 . The covariance matrices are given in Ref. [22].