Non-linear regression models for Approximate Bayesian Computation

Approximate Bayesian inference on the basis of summary statistics is well-suited to complex problems for which the likelihood is either mathematically or computationally intractable. However the methods that use rejection suffer from the curse of dimensionality when the number of summary statistics is increased. Here we propose a machine-learning approach to the estimation of the posterior density by introducing two innovations. The new method fits a nonlinear conditional heteroscedastic regression of the parameter on the summary statistics, and then adaptively improves estimation using importance sampling. The new algorithm is compared to the state-of-the-art approximate Bayesian methods, and achieves considerable reduction of the computational burden in two examples of inference in statistical genetics and in a queueing model.Comment: 4 figures; version 3 minor changes; to appear in Statistics and Computin

Finite one dimensional impenetrable Bose systems: Occupation numbers

Bosons in the form of ultra cold alkali atoms can be confined to a one dimensional (1d) domain by the use of harmonic traps. This motivates the study of the ground state occupations $\lambda_i$ of effective single particle states $\phi_i$, in the theoretical 1d impenetrable Bose gas. Both the system on a circle and the harmonically trapped system are considered. The $\lambda_i$ and $\phi_i$ are the eigenvalues and eigenfunctions respectively of the one body density matrix. We present a detailed numerical and analytic study of this problem. Our main results are the explicit scaled forms of the density matrices, from which it is deduced that for fixed $i$ the occupations $\lambda_i$ are asymptotically proportional to $\sqrt{N}$ in both the circular and harmonically trapped cases.Comment: 22 pages, 8 figures (.eps), uses REVTeX

Search For Trapped Antihydrogen

We present the results of an experiment to search for trapped antihydrogen atoms with the ALPHA antihydrogen trap at the CERN Antiproton Decelerator. Sensitive diagnostics of the temperatures, sizes, and densities of the trapped antiproton and positron plasmas have been developed, which in turn permitted development of techniques to precisely and reproducibly control the initial experimental parameters. The use of a position-sensitive annihilation vertex detector, together with the capability of controllably quenching the superconducting magnetic minimum trap, enabled us to carry out a high-sensitivity and low-background search for trapped synthesised antihydrogen atoms. We aim to identify the annihilations of antihydrogen atoms held for at least 130 ms in the trap before being released over ~30 ms. After a three-week experimental run in 2009 involving mixing of 10^7 antiprotons with 1.3 10^9 positrons to produce 6 10^5 antihydrogen atoms, we have identified six antiproton annihilation events that are consistent with the release of trapped antihydrogen. The cosmic ray background, estimated to contribute 0.14 counts, is incompatible with this observation at a significance of 5.6 sigma. Extensive simulations predict that an alternative source of annihilations, the escape of mirror-trapped antiprotons, is highly unlikely, though this possibility has not yet been ruled out experimentally.Comment: 12 pages, 7 figure

Model-independent search for CP violation in D0→K−K+π−π+ and D0→π−π+π+π− decays

A search for CP violation in the phase-space structures of D0 and View the MathML source decays to the final states K−K+π−π+ and π−π+π+π− is presented. The search is carried out with a data set corresponding to an integrated luminosity of 1.0 fb−1 collected in 2011 by the LHCb experiment in pp collisions at a centre-of-mass energy of 7 TeV. For the K−K+π−π+ final state, the four-body phase space is divided into 32 bins, each bin with approximately 1800 decays. The p-value under the hypothesis of no CP violation is 9.1%, and in no bin is a CP asymmetry greater than 6.5% observed. The phase space of the π−π+π+π− final state is partitioned into 128 bins, each bin with approximately 2500 decays. The p-value under the hypothesis of no CP violation is 41%, and in no bin is a CP asymmetry greater than 5.5% observed. All results are consistent with the hypothesis of no CP violation at the current sensitivity

Search for the lepton-flavor-violating decays Bs0→e±μ∓ and B0→e±μ∓

A search for the lepton-flavor-violating decays Bs0→e±μ∓ and B0→e±μ∓ is performed with a data sample, corresponding to an integrated luminosity of 1.0  fb-1 of pp collisions at √s=7  TeV, collected by the LHCb experiment. The observed number of Bs0→e±μ∓ and B0→e±μ∓ candidates is consistent with background expectations. Upper limits on the branching fractions of both decays are determined to be B(Bs0→e±μ∓)101  TeV/c2 and MLQ(B0→e±μ∓)>126  TeV/c2 at 95% C.L., and are a factor of 2 higher than the previous bounds

Measurements of long-range near-side angular correlations in sNN=5\sqrt{s_{\text{NN}}}=5TeV proton-lead collisions in the forward region

Two-particle angular correlations are studied in proton-lead collisions at a nucleon-nucleon centre-of-mass energy of $\sqrt{s_{\text{NN}}}=5$TeV, collected with the LHCb detector at the LHC. The analysis is based on data recorded in two beam configurations, in which either the direction of the proton or that of the lead ion is analysed. The correlations are measured in the laboratory system as a function of relative pseudorapidity, $\Delta\eta$, and relative azimuthal angle, $\Delta\phi$, for events in different classes of event activity and for different bins of particle transverse momentum. In high-activity events a long-range correlation on the near side, $\Delta\phi \approx 0$, is observed in the pseudorapidity range $2.0<\eta<4.9$. This measurement of long-range correlations on the near side in proton-lead collisions extends previous observations into the forward region up to $\eta=4.9$. The correlation increases with growing event activity and is found to be more pronounced in the direction of the lead beam. However, the correlation in the direction of the lead and proton beams are found to be compatible when comparing events with similar absolute activity in the direction analysed.Comment: All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2015-040.htm

Study of the production of Λb0\Lambda_b^0 and B‾0\overline{B}^0 hadrons in pppp collisions and first measurement of the Λb0→J/ψpK−\Lambda_b^0\rightarrow J/\psi pK^- branching fraction

The product of the $\Lambda_b^0$ ($\overline{B}^0$) differential production cross-section and the branching fraction of the decay $\Lambda_b^0\rightarrow J/\psi pK^-$ ($\overline{B}^0\rightarrow J/\psi\overline{K}^*(892)^0$) is measured as a function of the beauty hadron transverse momentum, $p_{\rm T}$, and rapidity, $y$. The kinematic region of the measurements is $p_{\rm T}<20~{\rm GeV}/c$ and $2.0<y<4.5$. The measurements use a data sample corresponding to an integrated luminosity of $3~{\rm fb}^{-1}$ collected by the LHCb detector in $pp$ collisions at centre-of-mass energies $\sqrt{s}=7~{\rm TeV}$ in 2011 and $\sqrt{s}=8~{\rm TeV}$ in 2012. Based on previous LHCb results of the fragmentation fraction ratio, $f_{\Lambda_B^0}/f_d$, the branching fraction of the decay $\Lambda_b^0\rightarrow J/\psi pK^-$ is measured to be \begin{equation*} \mathcal{B}(\Lambda_b^0\rightarrow J/\psi pK^-)= (3.17\pm0.04\pm0.07\pm0.34^{+0.45}_{-0.28})\times10^{-4}, \end{equation*} where the first uncertainty is statistical, the second is systematic, the third is due to the uncertainty on the branching fraction of the decay $\overline{B}^0\rightarrow J/\psi\overline{K}^*(892)^0$, and the fourth is due to the knowledge of $f_{\Lambda_b^0}/f_d$. The sum of the asymmetries in the production and decay between $\Lambda_b^0$ and $\overline{\Lambda}_b^0$ is also measured as a function of $p_{\rm T}$ and $y$. The previously published branching fraction of $\Lambda_b^0\rightarrow J/\psi p\pi^-$, relative to that of $\Lambda_b^0\rightarrow J/\psi pK^-$, is updated. The branching fractions of $\Lambda_b^0\rightarrow P_c^+(\rightarrow J/\psi p)K^-$ are determined.Comment: 29 pages, 19figures. All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2015-032.htm

Evidence for the strangeness-changing weak decay Ξb−→Λb0π−\Xi_b^-\to\Lambda_b^0\pi^-

Using a $pp$ collision data sample corresponding to an integrated luminosity of 3.0~fb$^{-1}$, collected by the LHCb detector, we present the first search for the strangeness-changing weak decay $\Xi_b^-\to\Lambda_b^0\pi^-$. No $b$ hadron decay of this type has been seen before. A signal for this decay, corresponding to a significance of 3.2 standard deviations, is reported. The relative rate is measured to be ${{f_{\Xi_b^-}}\over{f_{\Lambda_b^0}}}{\cal{B}}(\Xi_b^-\to\Lambda_b^0\pi^-) = (5.7\pm1.8^{+0.8}_{-0.9})\times10^{-4}$, where $f_{\Xi_b^-}$ and $f_{\Lambda_b^0}$ are the $b\to\Xi_b^-$ and $b\to\Lambda_b^0$ fragmentation fractions, and ${\cal{B}}(\Xi_b^-\to\Lambda_b^0\pi^-)$ is the branching fraction. Assuming $f_{\Xi_b^-}/f_{\Lambda_b^0}$ is bounded between 0.1 and 0.3, the branching fraction ${\cal{B}}(\Xi_b^-\to\Lambda_b^0\pi^-)$ would lie in the range from $(0.57\pm0.21)\%$ to $(0.19\pm0.07)\%$.Comment: 7 pages, 2 figures, All figures and tables, along with any supplementary material and additional information, are available at https://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2015-047.htm

BB flavour tagging using charm decays at the LHCb experiment

An algorithm is described for tagging the flavour content at production of neutral $B$ mesons in the LHCb experiment. The algorithm exploits the correlation of the flavour of a $B$ meson with the charge of a reconstructed secondary charm hadron from the decay of the other $b$ hadron produced in the proton-proton collision. Charm hadron candidates are identified in a number of fully or partially reconstructed Cabibbo-favoured decay modes. The algorithm is calibrated on the self-tagged decay modes $B^+ \to J/\psi \, K^+$ and $B^0 \to J/\psi \, K^{*0}$ using $3.0\mathrm{\,fb}^{-1}$ of data collected by the LHCb experiment at $pp$ centre-of-mass energies of $7\mathrm{\,TeV}$ and $8\mathrm{\,TeV}$. Its tagging power on these samples of $B \to J/\psi \, X$ decays is $(0.30 \pm 0.01 \pm 0.01) \%$.Comment: All figures and tables, along with any supplementary material and additional information, are available at http://lhcbproject.web.cern.ch/lhcbproject/Publications/LHCbProjectPublic/LHCb-PAPER-2015-027.htm

Branching fraction and CP asymmetry of the decays B+→K0Sπ+ and B+→K0SK+

An analysis of B+ → K0 Sπ+ and B+ → K0 S K+ decays is performed with the LHCb experiment. The pp collision data used correspond to integrated luminosities of 1 fb−1 and 2 fb−1 collected at centre-ofmass energies of √ s = 7 TeV and √ s = 8 TeV, respectively. The ratio of branching fractions and the direct CP asymmetries are measured to be B(B+ → K0 S K+ )/B(B+ → K0 Sπ+ ) = 0.064 ± 0.009 (stat.) ± 0.004 (syst.), ACP(B+ → K0 Sπ+ ) = −0.022 ± 0.025 (stat.) ± 0.010 (syst.) and ACP(B+ → K0 S K+ ) = −0.21 ± 0.14 (stat.) ± 0.01 (syst.). The data sample taken at √ s = 7 TeV is used to search for B+ c → K0 S K+ decays and results in the upper limit ( fc · B(B+ c → K0 S K+ ))/( fu · B(B+ → K0 Sπ+ )) < 5.8 × 10−2 at 90% confidence level, where fc and fu denote the hadronisation fractions of a ¯b quark into a B+ c or a B+ meson, respectively