Fast and flexible data structures for the LHCb Run 3 software trigger

Starting in 2022, the upgraded LHCb detector will collect data with a pure software trigger. In its first stage, reducing the rate from 30MHz to about 1MHz, GPUs are used to reconstruct and trigger on B and D meson topologies and high-pT objects in the event. In its second stage, a CPU farm is used to reconstruct the full event and perform candidate selections, which are persisted for offline use with an output rate of about 10 GB/s. Fast data processing, flexible and custom-designed data structures tailored for SIMD architectures and efficient storage of the intermediate data at various steps of the processing pipeline onto persistent media, e.g. tapes is essential to guarantee the full physics program of LHCb. In this talk, we will present the event model and data persistency developments for the trigger of LHCb in Run 3. Particular emphasize will be given to the novel software-design aspects with respect to the Run 1+2 data taking, the performance improvements which can be achieved and the experience of restructuring a major part of the reconstruction software in a large HEP experiment.Comment: Connecting The Dots (CTD 2022), Princeton, United States Of America, 31 May - 2 Jun 202

Multidifferential study of identified charged hadron distributions in ZZ-tagged jets in proton-proton collisions at s=\sqrt{s}=13 TeV

Jet fragmentation functions are measured for the first time in proton-proton collisions for charged pions, kaons, and protons within jets recoiling against a $Z$ boson. The charged-hadron distributions are studied longitudinally and transversely to the jet direction for jets with transverse momentum 20 $< p_{\textrm{T}} < 100$ GeV and in the pseudorapidity range $2.5 < \eta < 4$. The data sample was collected with the LHCb experiment at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 1.64 fb$^{-1}$. Triple differential distributions as a function of the hadron longitudinal momentum fraction, hadron transverse momentum, and jet transverse momentum are also measured for the first time. This helps constrain transverse-momentum-dependent fragmentation functions. Differences in the shapes and magnitudes of the measured distributions for the different hadron species provide insights into the hadronization process for jets predominantly initiated by light quarks.Comment: All figures and tables, along with machine-readable versions and any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-PAPER-2022-013.html (LHCb public pages

Study of the B−→Λc+Λˉc−K−B^{-} \to \Lambda_{c}^{+} \bar{\Lambda}_{c}^{-} K^{-} decay

The decay $B^{-} \to \Lambda_{c}^{+} \bar{\Lambda}_{c}^{-} K^{-}$ is studied in proton-proton collisions at a center-of-mass energy of $\sqrt{s}=13$ TeV using data corresponding to an integrated luminosity of 5 $\mathrm{fb}^{-1}$ collected by the LHCb experiment. In the $\Lambda_{c}^+ K^{-}$ system, the $\Xi_{c}(2930)^{0}$ state observed at the BaBar and Belle experiments is resolved into two narrower states, $\Xi_{c}(2923)^{0}$ and $\Xi_{c}(2939)^{0}$, whose masses and widths are measured to be $$m(\Xi_{c}(2923)^{0}) = 2924.5 \pm 0.4 \pm 1.1 \,\mathrm{MeV}, \\ m(\Xi_{c}(2939)^{0}) = 2938.5 \pm 0.9 \pm 2.3 \,\mathrm{MeV}, \\ \Gamma(\Xi_{c}(2923)^{0}) = \phantom{000}4.8 \pm 0.9 \pm 1.5 \,\mathrm{MeV},\\ \Gamma(\Xi_{c}(2939)^{0}) = \phantom{00}11.0 \pm 1.9 \pm 7.5 \,\mathrm{MeV},$$ where the first uncertainties are statistical and the second systematic. The results are consistent with a previous LHCb measurement using a prompt $\Lambda_{c}^{+} K^{-}$ sample. Evidence of a new $\Xi_{c}(2880)^{0}$ state is found with a local significance of $3.8\,\sigma$, whose mass and width are measured to be $2881.8 \pm 3.1 \pm 8.5\,\mathrm{MeV}$ and $12.4 \pm 5.3 \pm 5.8 \,\mathrm{MeV}$, respectively. In addition, evidence of a new decay mode $\Xi_{c}(2790)^{0} \to \Lambda_{c}^{+} K^{-}$ is found with a significance of $3.7\,\sigma$. The relative branching fraction of $B^{-} \to \Lambda_{c}^{+} \bar{\Lambda}_{c}^{-} K^{-}$ with respect to the $B^{-} \to D^{+} D^{-} K^{-}$ decay is measured to be $2.36 \pm 0.11 \pm 0.22 \pm 0.25$, where the first uncertainty is statistical, the second systematic and the third originates from the branching fractions of charm hadron decays.Comment: All figures and tables, along with any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-PAPER-2022-028.html (LHCb public pages

Measurement of the ratios of branching fractions R(D∗)\mathcal{R}(D^{*}) and R(D0)\mathcal{R}(D^{0})

The ratios of branching fractions $\mathcal{R}(D^{*})\equiv\mathcal{B}(\bar{B}\to D^{*}\tau^{-}\bar{\nu}_{\tau})/\mathcal{B}(\bar{B}\to D^{*}\mu^{-}\bar{\nu}_{\mu})$ and $\mathcal{R}(D^{0})\equiv\mathcal{B}(B^{-}\to D^{0}\tau^{-}\bar{\nu}_{\tau})/\mathcal{B}(B^{-}\to D^{0}\mu^{-}\bar{\nu}_{\mu})$ are measured, assuming isospin symmetry, using a sample of proton-proton collision data corresponding to 3.0 fb${ }^{-1}$ of integrated luminosity recorded by the LHCb experiment during 2011 and 2012. The tau lepton is identified in the decay mode $\tau^{-}\to\mu^{-}\nu_{\tau}\bar{\nu}_{\mu}$. The measured values are $\mathcal{R}(D^{*})=0.281\pm0.018\pm0.024$ and $\mathcal{R}(D^{0})=0.441\pm0.060\pm0.066$, where the first uncertainty is statistical and the second is systematic. The correlation between these measurements is $\rho=-0.43$. Results are consistent with the current average of these quantities and are at a combined 1.9 standard deviations from the predictions based on lepton flavor universality in the Standard Model.Comment: All figures and tables, along with any supplementary material and additional information, are available at https://cern.ch/lhcbproject/Publications/p/LHCb-PAPER-2022-039.html (LHCb public pages

Performance optimization for the LHCb experiment

The LHCb experiment, at CERN, is preparing a major upgrade of its detector and a change from an hardware-based to a fully software-based trigger system. It is now facing the challenge of being able to process incoming events at a rate of 30 million events per second. To cope with this massive data input, the software must be optimized to use the processing power of the filtering farm more efficiently. This thesis focus on the first algorithm of LHCb’s High Level Trigger software: the Vertex Locator (VELO) reconstruction algorithm. The VELO is the first detector encountered by particles, directly surrounding the interaction region. Its goal is to find the initial track candidate that are then followed through the other layers of the LHCb detector with a good enough resolution that they could also be used to locate the origin of the collisions. The first step of this algorithm is to prepare the data by grouping pixels of the silicon sensors into hits; this process is called connected component analysis (CCA). This thesis presents multiple new CCA algorithms for both CPU and GPU architectures. The first algorithm, HA4, was developed at the very start of this thesis and improved the state-of-the-art in connected component labeling on GPUs, as well as being the first efficient implementation of connected component analysis on GPUs. The second algorithm is a GPU port of the FLSL SIMD CPU algorithm, inspired by the LSL algorithm. FLSL on GPUs improved upon HA4 by reducing the memory accesses conflicts that are especially presents on new hardware with a lot of cores. Along with FLSL, two other optimisations aimed at further reducing conflicts are presented and evaluated. On CPU, two new algorithms were made for this thesis. The first one is a modification of the classic Rosenfeld algorithm to use SIMD. The second one is a new algorithm, named SparseCCL, which takes advantage of the sparsity of the input images. A new VELO reconstruction algorithm using SIMD is presented, that enable LHCb to process events in real time and improve the quality of the reconstruction. The SIMDWrapper library, developed for the new VELO algorithm, is now part of LHCb’s software and is used in other algorithms

Fast and flexible data structures for the LHCb Run 3 software trigger

Starting in 2022, the upgraded LHCb detector will collect data with a pure software trigger. In its first stage, reducing the rate from 30MHz to about 1MHz, GPUs are used to reconstruct and trigger on B and D meson topologies and high-$p_T$ objects in the event. In its second stage, a CPU farm is used to reconstruct the full event and perform candidate selections, which are persisted for offline use with an output rate of about 10 GB/s. Fast data processing, flexible and custom-designed data structures tailored for SIMD architectures and efficient storage of the intermediate data at various steps of the processing pipeline onto persistent media, e.g. tapes is essential to guarantee the full physics program of LHCb. In this talk, we will present the event model and data persistency developments for the trigger of LHCb in Run 3. Particular emphasize will be given to the novel software-design aspects with respect to the Run 1+2 data taking, the performance improvements which can be achieved and the experience of restructuring a major part of the reconstruction software in a large HEP experiment

Improved measurement of CPCP violation parameters in Bs0→J/ψK+K−B_s^0\to J/\psi K^+K^- decays in the vicinity of the ϕ(1020)\phi(1020) resonance

The decay-time-dependent $CP$ asymmetry in $B_s^0\to J/\psi(\to \mu^+\mu^-) K^+ K^-$ decays is measured using proton-proton collision data, corresponding to an integrated luminosity of $6 {\rm fb}^{-1}$, collected with the LHCb detector at a center-of-mass energy of 13 TeV. Using a sample of approximately 349 000 $B_s^0$ signal decays with an invariant $K^+ K^-$ mass in the vicinity of the $\phi(1020)$ resonance, the $CP$-violating phase $\phi_s$ is measured, along with the difference in decay widths of the light and heavy mass eigenstates of the $B_s^0$-$\overline{B}_s^0$ system, $\Delta\Gamma_s$, and the difference of the average $B_s^0$ and $B^0$ meson decay widths, $\Gamma_s-\Gamma_d$. The values obtained are $\phi_s = -0.039 \pm 0.022 \pm 0.006$ rad, $\Delta\Gamma_s = 0.0845 \pm 0.0044 \pm 0.0024 ~{\rm ps}^{-1}$ and $\Gamma_s-\Gamma_d = -0.056^{\:+\:0.0013}_{\:-\:0.0015} \pm 0.0014 ~{\rm ps}^{-1}$, where the first uncertainty is statistical and the second systematic. These are the most precise single measurements to date and are consistent with expectations based on the Standard Model and with the previous LHCb analyses of this decay. These results are combined with previous independent LHCb measurements. The phase $\phi_s$ is also measured independently for each polarization state of the $K^+K^-$ system and shows no evidence for polarization dependence.The decay-time-dependent $CP$ asymmetry in $B^0_s\to J/\psi(\to \mu^{+}\mu^{-}) K^{+}K^{-}$ decays is measured using proton-proton collision data, corresponding to an integrated luminosity of 6 $fb^{-1}$, collected with the LHCb detector at a center-of-mass energy of 13 TeV. Using a sample of approximately 349 000 $B^{0}_{s}$ signal decays with an invariant $K^{+}K^{-}$ mass in the vicinity of the $\phi(1020)$ resonance, the $CP$-violating phase $\phi_s$ is measured, along with the difference in decay widths of the light and heavy mass eigenstates of the $B^0_s$-$\bar{B}^0_s$ system, $\Delta\Gamma_s$, and the difference of the average $B^0_s$ and $B^0$ meson decay widths, $\Gamma_s-\Gamma_d$. The values obtained are $\phi_s = \ -0.039 \pm 0.022 \pm 0.006$ rad, $\Delta\Gamma_s = 0.0845 \pm 0.0044 \pm 0.0024$ ps$^{-1}$ and $\Gamma_s-\Gamma_d = -0.0056 ^{+ 0.0013}_{-0.0015} \pm 0.0014$ ps$^{-1}$, where the first uncertainty is statistical and the second systematic. These are the most precise single measurements to date and are consistent with expectations based on the Standard Model and with the previous LHCb analyses of this decay. These results are combined with previous independent LHCb measurements. The phase $\phi_s$ is also measured independently for each polarization state of the $K^{+}K^{-}$ system and shows no evidence for polarization dependence

Search for Bc+→π+μ+μ−B_c^+\to\pi^+\mu^+\mu^- decays and measurement of the branching fraction ratio B(Bc+→ψ(2S)π+)/B(Bc+→J/ψπ+){\cal B}(B_c^+\to\psi(2S)\pi^+)/{\cal B}(B_c^+\to J/\psi \pi^+)

International audienceThe first search for nonresonant $B_c^+\to\pi^+\mu^+\mu^-$ decays is reported. The analysis uses proton-proton collision data collected with the LHCb detector between 2011 and 2018, corresponding to an integrated luminosity of 9 fb$^{-1}$. No evidence for an excess of signal events over background is observed and an upper limit is set on the branching fraction ratio ${\cal B}(B_c^+\to\pi^+\mu^+\mu^-)/{\cal B}(B_c^+\to J/\psi \pi^+) < 2.1\times 10^{-4}$ at $90\%$ confidence level. Additionally, an updated measurement of the ratio of the $B_c^+\to\psi(2S)\pi^+$ and $B_c^+\to J/\psi \pi^+$ branching fractions is reported. The ratio ${\cal B}(B_c^+\to\psi(2S)\pi^+)/{\cal B}(B_c^+\to J/\psi \pi^+)$ is measured to be $0.254\pm 0.018 \pm 0.003 \pm 0.005$, where the first uncertainty is statistical, the second systematic, and the third is due to the uncertainties on the branching fractions of the leptonic $J/\psi$ and $\psi(2S)$ decays. This measurement is the most precise to date and is consistent with previous LHCb results

Search for Bc+→π+μ+μ−B_c^+\to\pi^+\mu^+\mu^- decays and measurement of the branching fraction ratio B(Bc+→ψ(2S)π+)/B(Bc+→J/ψπ+){\cal B}(B_c^+\to\psi(2S)\pi^+)/{\cal B}(B_c^+\to J/\psi \pi^+)

International audienceThe first search for nonresonant $B_c^+\to\pi^+\mu^+\mu^-$ decays is reported. The analysis uses proton-proton collision data collected with the LHCb detector between 2011 and 2018, corresponding to an integrated luminosity of 9 fb$^{-1}$. No evidence for an excess of signal events over background is observed and an upper limit is set on the branching fraction ratio ${\cal B}(B_c^+\to\pi^+\mu^+\mu^-)/{\cal B}(B_c^+\to J/\psi \pi^+) < 2.1\times 10^{-4}$ at $90\%$ confidence level. Additionally, an updated measurement of the ratio of the $B_c^+\to\psi(2S)\pi^+$ and $B_c^+\to J/\psi \pi^+$ branching fractions is reported. The ratio ${\cal B}(B_c^+\to\psi(2S)\pi^+)/{\cal B}(B_c^+\to J/\psi \pi^+)$ is measured to be $0.254\pm 0.018 \pm 0.003 \pm 0.005$, where the first uncertainty is statistical, the second systematic, and the third is due to the uncertainties on the branching fractions of the leptonic $J/\psi$ and $\psi(2S)$ decays. This measurement is the most precise to date and is consistent with previous LHCb results

Measurement of CP violation in B0→ψ(→ℓ+ℓ−)KS0(→π+π−)B^0\to\psi(\to\ell^+\ell^-)K^0_S(\to\pi^+\pi^-) decays

International audienceA measurement of time-dependent CP violation in the decays of $B^0$ and $\overline{B}^0$ mesons to the final states $J/\psi(\to\mu^+\mu^-)K^0_S$, $\psi(2S)(\to\mu^+\mu^-)K^0_S$ and $J/\psi(\to e^+e^-)K^0_S$ with $K^0_S\to\pi^+\pi^-$ is presented. The data correspond to an integrated luminosity of 6 fb${}^{-1}$ collected at a centre-of-mass energy of $\sqrt{s}=13$ TeV with the LHCb detector. The CP-violation parameters are measured to be \begin{align*} S_{\psi K^0_S} &= 0.717 \pm 0.013 (\text{stat}) \pm 0.008 (\text{syst}), \\ C_{\psi K^0_S} &= 0.008 \pm 0.012 (\text{stat}) \pm 0.003 (\text{syst}). \end{align*} This measurement of $S_{\psi K^0_S}$ represents the most precise single measurement of the CKM angle $\beta$ to date and is more precise than the current world average. In addition, measurements of the CP-violation parameters of the individual channels are reported and a combination with the LHCb Run 1 measurements is performed