<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Big Data</journal-id>
<journal-title>Frontiers in Big Data</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Big Data</abbrev-journal-title>
<issn pub-type="epub">2624-909X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fdata.2022.1008737</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Big Data</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Effect of the high-level trigger for detecting long-lived particles at LHCb</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Calefice</surname> <given-names>Lukas</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Hennequin</surname> <given-names>Arthur</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Henry</surname> <given-names>Louis</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1974290/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Jashal</surname> <given-names>Brij</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2058328/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Mendoza</surname> <given-names>Diego</given-names></name>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2036716/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Oyanguren</surname> <given-names>Arantza</given-names></name>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1601843/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Sanderswood</surname> <given-names>Izaac</given-names></name>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Sierra</surname> <given-names>Carlos V&#x000E1;zquez</given-names></name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhuo</surname> <given-names>Jiahui</given-names></name>
<xref ref-type="aff" rid="aff5"><sup>5</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1994727/overview"/>
</contrib>
<contrib contrib-type="author">
<collab>Part of LHCb-RTA Collaboration</collab>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Laboratoire de Physique Nucl&#x000E9;aire et de Hautes Energies, Sorbonne Universit&#x000E9;, CNRS/IN2P3</institution>, <addr-line>Paris</addr-line>, <country>France</country></aff>
<aff id="aff2"><sup>2</sup><institution>Fakult&#x000E4;t Physik, Technical University Dortmund</institution>, <addr-line>Dortmund</addr-line>, <country>Germany</country></aff>
<aff id="aff3"><sup>3</sup><institution>Laboratory for Nuclear Science, Massachusetts Institute of Technology</institution>, <addr-line>Cambridge, MA</addr-line>, <country>United States</country></aff>
<aff id="aff4"><sup>4</sup><institution>European Organization for Nuclear Research (CERN)</institution>, <addr-line>Geneva</addr-line>, <country>Switzerland</country></aff>
<aff id="aff5"><sup>5</sup><institution>Instituto de F&#x000ED;sica Corpuscular (IFIC), Consejo Superior de Investigaciones Cient&#x000ED;ficas, University of Valencia</institution>, <addr-line>Valencia</addr-line>, <country>Spain</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Michael Spannowsky, Durham University, United Kingdom</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Letizia Lusito, Inside Climate Service, Italy; Alexander Kish, Fermi National Accelerator Laboratory (DOE), United States</p></fn>
<corresp id="c001">&#x0002A;Correspondence: Arantza Oyanguren <email>arantza.oyanguren&#x00040;ific.uv.es</email></corresp>
<fn fn-type="other" id="fn001"><p>This article was submitted to Big Data and AI in High Energy Physics, a section of the journal Frontiers in Big Data</p></fn></author-notes>
<pub-date pub-type="epub">
<day>07</day>
<month>11</month>
<year>2022</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>5</volume>
<elocation-id>1008737</elocation-id>
<history>
<date date-type="received">
<day>01</day>
<month>08</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>03</day>
<month>10</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2022 Calefice, Hennequin, Henry, Jashal, Mendoza, Oyanguren, Sanderswood, Sierra, Zhuo and part of LHCb-RTA Collaboration.</copyright-statement>
<copyright-year>2022</copyright-year>
<copyright-holder>Calefice, Hennequin, Henry, Jashal, Mendoza, Oyanguren, Sanderswood, Sierra, Zhuo and part of LHCb-RTA Collaboration</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license> </permissions>
<abstract>
<p>Long-lived particles (LLPs) show up in many extensions of the Standard Model, but they are challenging to search for with current detectors, due to their very displaced vertices. This study evaluated the ability of the trigger algorithms used in the Large Hadron Collider beauty (LHCb) experiment to detect long-lived particles and attempted to adapt them to enhance the sensitivity of this experiment to undiscovered long-lived particles. A model with a Higgs portal to a dark sector is tested, and the sensitivity reach is discussed. In the LHCb tracking system, the farthest tracking station from the collision point is the scintillating fiber tracker, the SciFi detector. One of the challenges in the track reconstruction is to deal with the large amount of and combinatorics of hits in the LHCb detector. A dedicated algorithm has been developed to cope with the large data output. When fully implemented, this algorithm would greatly increase the available statistics for any long-lived particle search in the forward region and would additionally improve the sensitivity of analyses dealing with Standard Model particles of large lifetime, such as <inline-formula><mml:math id="M1"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> or &#x0039B;<sup>0</sup> hadrons.</p></abstract>
<kwd-group>
<kwd>LHCb</kwd>
<kwd>trigger</kwd>
<kwd>real time analysis</kwd>
<kwd>long-lived particles</kwd>
<kwd>GPU</kwd>
<kwd>SciFi</kwd>
<kwd>beyond standard physics</kwd>
</kwd-group>
<counts>
<fig-count count="12"/>
<table-count count="0"/>
<equation-count count="3"/>
<ref-count count="16"/>
<page-count count="13"/>
<word-count count="5341"/>
</counts>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1. Introduction</title>
<p>The Large Hadron Collider beauty (LHCb) (Alves et al., <xref ref-type="bibr" rid="B5">2008</xref>) is one of the experiments based at the CERN Large Hadron Collider (LHC) that aims to explore the limitations of the well-tested Standard Model of particle physics. The imbalance of matter and antimatter in the Universe, the presence of neutrino masses, and the nature of dark matter are some examples of open questions for models of undiscovered physics beyond the Standard Model (BSM) that has been proposed.</p>
<p>The observation of the new processes that may be currently indiscernible is directly related to the accelerator luminosity. To meet the new challenges posed by Run 3 data taking at the LHC, the LHCb detector has undergone a major upgrade (LHCb, <xref ref-type="bibr" rid="B10">2014a</xref>) to run at a higher luminosity of <inline-formula><mml:math id="M2"><mml:msub><mml:mrow><mml:mrow><mml:mi mathvariant="script">L</mml:mi></mml:mrow></mml:mrow><mml:mrow><mml:mstyle class="mbox"><mml:mtext>inst</mml:mtext></mml:mstyle></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mn>2</mml:mn><mml:mo>&#x000D7;</mml:mo><mml:mn>1</mml:mn><mml:msup><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mn>33</mml:mn></mml:mrow></mml:msup><mml:mtext>c</mml:mtext><mml:msup><mml:mrow><mml:mtext>m</mml:mtext></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mtext>s</mml:mtext></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:math></inline-formula>. This luminosity corresponds to an averaged bunch collision frequency of 30 MHz, from where the full detector readout produces around 4 TB/s of raw data. To handle such detector occupancy and to search for rare processes among well-known SM events are computing challenges managed by the <italic>trigger</italic> system. The trigger system decides what information is recorded for later analysis. It must use simplified criteria to maintain a minimal latency. The previously used L0 hardware system is not suitable for such a high luminosity. Therefore, the LHCb trigger has been redesigned to remove the L0 hardware system and rely on a fully-software based architecture, which consists of two stages: HLT1 and HLT2 (LHCb, <xref ref-type="bibr" rid="B12">2020</xref>), explained in the following sections.</p>
</sec>
<sec sec-type="materials and methods" id="s2">
<title>2. Materials and methods</title>
<sec>
<title>2.1. The LHCb detector</title>
<p>The upgraded LHCb detector, which will be operating during Run 3, detailed by Bediaga et al. (<xref ref-type="bibr" rid="B6">2012</xref>), is summarized below. The most important detector improvements as compared to the previous detector (Alves et al., <xref ref-type="bibr" rid="B5">2008</xref>) are the new tracking system and the readout architecture. The readout architecture uses a purely software-based trigger that enables data to be recorded at five times the previous rate, and increases the acceptance of purely hadronic b decays by up to a factor of 2. LHCb is comprised of a tracking system (Vertex Locator, Upstream Tracker, and SciFi tracker), as well as two Ring Imaging Cherenkov detectors, a hadronic calorimeter, and an electromagnetic and muon chambers. These are summarized below.</p>
<p>The VErtex LOcator (VELO) is based on pixelated-silicon sensors and is critical to determine the decay vertices of b and c flavored hadrons. The Upstream Tracker, UT, contains vertically-segmented silicon strips and continues the tracking upstream of the VELO. It is also used to determine the momentum of charged particles with &#x00394;<italic>p</italic>/<italic>p</italic> &#x02248; 10% and is useful for removing low-momentum tracks from being extrapolated downstream, thus speeding up the software trigger by about a factor of 3. Tracking after the magnet is handled by the new scintillating fiber-based Scintillating Fiber detector (SciFi). Two Ring Imaging CHerenkov (RICH) detectors supply particle identification. RICH1 is mainly for lower momentum particles and RICH2 is for higher momentum ones. The Electromagnetic Calorimeter (ECAL) identifies electrons and reconstructs photons and neutral pions. The Hadronic Calorimeter (HCAL) measures the energy deposits of hadrons, and the muon chambers M2-M5 are mostly used for muon identification. <xref ref-type="fig" rid="F1">Figure 1</xref> shows the LHCb upgrade detector.</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p>Sketch of the new LHCb detector.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0001.tif"/>
</fig>
</sec>
<sec>
<title>2.2. The LHCb event-builder</title>
<p>The main goal of the upgraded LHCb data-acquisition system is to enable the event filter to efficiently select events based on the full data readout from the detectors. The heart of the LHCb online system is the event builder, which is capable of collecting and aggregating data fragments coming from the FPGA front-ends of the various sub-detectors at the full collision rate of 30 MHz (LHCb, <xref ref-type="bibr" rid="B11">2014b</xref>). The data size of a full event during the Run 3 is expected to be as high as 150 kB. Therefore, taking into account the data packed and the transfer protocol overheads, the new event builder has to sustain a throughput of up to 40 Tb/s. The event builder is essentially made up of three components: the readout units, the builder units, and the network interconnecting them. Each readout unit receives the event fragments from the detector point-to-point links and makes them available on the EB network. For each event, one of the builder units gathers all the event fragments from all the readout units, assembles them into a complete event, and writes it into HLT1 input buffer for event filtering. The present event-buider architecture relies on 163 EB servers interconnected through a 200 Gbits/s bi-directional network, as shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. The EB nodes have free PCIe slots that host GPU (Graphic Processing Units) cards to process the HLT1. This has the advantage of reducing the amount of data sent to the Event Filter Farm (EFF) and therefore reduces network costs.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p>The architecture of the new LHCb data readout and event-buider system.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0002.tif"/>
</fig>
</sec>
<sec>
<title>2.3. Track types at LHCb</title>
<p>In LHCb, the tracking system for charged particles consists of three subsystems: VELO, UT, and SciFi. The magnet is also necessary to bend particles to estimate their momentum. Several track types are defined depending on the subdetectors involved in the reconstruction, as shown in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p>Track types in LHCb, as defined in LHCb (<xref ref-type="bibr" rid="B10">2014a</xref>).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0003.tif"/>
</fig>
<p>The main track types considered for physics analyses are as follows:</p>
<list list-type="bullet">
<list-item><p>Long tracks, which have information from at least the VELO and the SciFi, and possibly the UT, are the main tracks used in physics analyses and at all stages of the trigger;</p></list-item>
<list-item><p>Downstream tracks, which have information from the UT and the SciFi, but not VELO, typically correspond to decay products of <inline-formula><mml:math id="M3"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> and &#x0039B;<sup>0</sup> hadron decays; and</p></list-item>
<list-item><p>T-tracks, which only have hits from the SciFi, are typically not included either in the trigger or in analysis. Nevertheless, their potential for physics has been recently outlined (LHCb, <xref ref-type="bibr" rid="B13">2022</xref>).</p></list-item>
</list>
<p>When simulating collision data, tracks meeting certain thresholds are defined to be reconstructible and have an assigned type based on the subdetector reconstructibility. This is based on the existence of reconstructed detector digits or clusters in the emulated detector, which are matched to the simulated particles if the detector hits they originated from are properly linked (Li et al., <xref ref-type="bibr" rid="B14">2021</xref>). One can distinguish reconstructibility in the following subdetectors:</p>
<list list-type="bullet">
<list-item><p>VELO: at least 3 pixel sensors and at least 1 digit each.</p></list-item>
<list-item><p>UT: at least 2 clusters where 1 cluster has to be in layer one or two and 1 in layer three or four. The clusters can be x or stereo clusters.</p></list-item>
<list-item><p>SciFi: at least 1 x cluster and 1 stereo cluster in each of the 3 SciFi stations.</p></list-item>
</list>
<p>Requirements for reconstructible Long tracks implies VELO and SciFi reconstructibility, Downstream tracks must satisfy the UT and SciFi reconstructibility, and T-tracks only requires SciFi reconstrucibility.</p>
</sec>
<sec>
<title>2.4. The high-level trigger</title>
<p>The trigger of the LHCb detector readout in Run 3 and later is fully software based and comprises two levels: HLT1 and HLT2, described in LHCb (<xref ref-type="bibr" rid="B12">2020</xref>). Most notably, the HLT1 level has to be executed at a 30 MHz rate, which, thus, suffers from heavy constraints on timing for event reconstruction.</p>
<p>The HLT1 is the first trigger stage and performs partial event reconstruction to reduce the data rate. Tracking algorithms play a key role on fast event decisions, and the fact that they are inherently parallelizable processes suggests a way to increase the trigger performance. Thus, the HLT1 has been implemented on a number of GPUs through the Allen project (Aaij et al., <xref ref-type="bibr" rid="B2">2020</xref>), which manages 4 TB/s and reduces the data rate by a factor of 30. After this initial selection, data are passed to a buffer system that allows real-time calibration and alignment of the detector. This is used for the full and improved event reconstruction carried out by HLT2. <xref ref-type="fig" rid="F4">Figure 4</xref> shows the dataflow at LHCb in the first trigger level.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p>Dataflow for the first level of the trigger in the upgraded LHCb detector.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0004.tif"/>
</fig>
<p>Due to timing constraints, the current LHCb implementation in this stage is based on partial reconstruction and focuses solely on <italic>long</italic> tracks, i.e., tracks that have hits in the VELO. This trigger thus significantly affects the identification of particles with long lifetimes, in particular LLP searches in LHCb, where some of the final-state particles are created further than roughly a meter away from the IP and thus outside of the VELO acceptance.</p>
<p>Implementing new fast tracking algorithms in HLT1, allowing direct access to Downstream and T tracks from HLT1 selections, could solve this actual LHCb inability in LLPs reconstruction.</p>
</sec>
<sec>
<title>2.5. The SciFi seeding algorithm</title>
<p>Hybrid seeding, detailed in Aiola et al. (<xref ref-type="bibr" rid="B4">2021</xref>), is an iterative reconstruction algorithm that uses information from the SciFi subdetector of LHCb to reconstruct track segments, also called T-tracks. It exploits the four-layered structure of the SciFi stations, as shown in <xref ref-type="fig" rid="F5">Figure 5</xref>, which have two vertical outer (X) layers and two inclined inner (U, V) layers, which are inclined by &#x000B1;5&#x000B0; to determine the <italic>y</italic>-position of the trajectory.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p>Representation of a typical T-track trajectory in the SciFi. The x-z trajectory is modeled as a parabola with a fixed cubic correction, while the y-z trajectory is a line.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0005.tif"/>
</fig>
<p>This algorithm operates in iterations or cases, starting from different layers to cover for hit inefficiency. Between each iteration, hits that were used to form a good-quality track are removed from the container, which allows the loosening of tolerance windows and track quality requirements as the environment is cleaned up by the algorithm.</p>
<p>Each iteration can be separated into two parts: the <italic>x</italic> &#x02212; <italic>z</italic> track finding and the U/V hit addition. The trajectory of a track, as shown in <xref ref-type="fig" rid="F5">Figure 5</xref>, can be approximated as a parabola with a fixed cubic correction in the <italic>x</italic> &#x02212; <italic>z</italic> plane and a straight line in the <italic>y</italic> &#x02212; <italic>z</italic> plane. Due to the geometry of the SciFi detector, y information can only be accessed using the <italic>x</italic>(<italic>z</italic>) equation.</p>
<p><xref ref-type="fig" rid="F6">Figure 6</xref> shows the principle of the triplet finding in the <italic>x</italic> &#x02212; <italic>z</italic> plane, which is the most computationally expensive part of the algorithm. Each hit in the first considered X layer is propagated onto a second X layer, assuming the track has infinite momentum and comes from the origin. A tolerance window is opened around that projection, assuming a minimum momentum, and all hits within this window are collected. The difference between the position of the second hit and the infinite momentum projection is a measurement of the momentum, which allows to estimate the bending within SciFi and to look for the third hit. Each triplet then provides the full track equation, and remaining hits are searched for. Tracklets are then selected according to quality criteria such as &#x003C7;<sup>2</sup> of the fit to the trajectory and the number of hits.</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p>Formation of the triplets in x-z search. The position of a hit in a first layer is projected onto a second layer assuming infinite momentum, and all hits around that projection, given a momentum-dependent tolerance, are used to form 2-hit combination. A third hit is then looked for in a third layer using the estimated momentum from the 2-hits combination.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0006.tif"/>
</fig>
<p>Given the <italic>x</italic>(<italic>z</italic>) track equation, hits on tilted U and V layers can be converted to a <italic>y</italic> measurement. Considering that we expect tracks to be a straight line coming from the origin, these <italic>y</italic> measurements are converted to a slope measurement, and the algorithm looks for the accumulation of similar slope measurements in different U/V layers. This is called Hough transformation and clustering. Each Hough cluster which has enough hits is then tested by adding it to the <italic>x</italic> &#x02212; <italic>z</italic> track and performing a global fit to the trajectory. Tracks that do not pass quality requirements on that fit are rejected. The seeds reconstructed by the Hybrid seeding are then used to form more complex tracks. They are matched with VELO segments to form Long tracks or with hits in the UT to form Downstream tracks. The remaining tracklets are then called T-tracks.</p>
<p>Recently, an optimization and high-throughput configuration of this algorithm has been implemented in the new GPU architecture of the LHCb trigger system. The main features of this new seeding is the high level of parallelism and the reimplementation of the U/V hit addition as a sequential hit search, similar to that performed in the <italic>x</italic> &#x02212; <italic>z</italic> part. Additionally, some quality criteria have been tightened. <xref ref-type="fig" rid="F7">Figure 7</xref> shows the preliminary physics performance of the Hybrid seeding algorithm as implemented on GPUs on Long tracks coming from the <italic>B</italic><sub><italic>s</italic></sub> &#x02192; <italic>&#x003D5;&#x003D5;</italic> decay channel, where the <italic>&#x003D5;</italic> meson is reconstructed from two kaons of opposite signs that traverse the SciFi detector.</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p>Efficiency of the Hybrid seeding reconstruction on GPUs as a function of the track momentum of Long tracks produced from <italic>B</italic><sub><italic>s</italic></sub> &#x02192; <italic>&#x003D5;&#x003D5;</italic> decays [LHCb-FIGURE-2022-010].</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0007.tif"/>
</fig>
</sec>
<sec>
<title>2.6. Higgs portal to dark matter</title>
<p>Some models of physics beyond the Standard Model (BSM) predict the Standard Model (SM) Higgs field serving as the portal to a dark sector, which could accommodate dark matter candidates (Kachanovich, Aliaksei et al., <xref ref-type="bibr" rid="B9">2020</xref>). One of the simplest and well-known models in this regard predicts the existence of a mixed state between a new scalar low-mass boson (H&#x00027;) and the SM Higgs (H), regulated by the mixing strength &#x003B8;:</p>
<disp-formula id="E1"><label>(1)</label><mml:math id="M4"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>h</mml:mi><mml:mo>=</mml:mo><mml:mi>H</mml:mi><mml:mtext>cos</mml:mtext><mml:mi>&#x003B8;</mml:mi><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:mtext>sin</mml:mtext><mml:mi>&#x003B8;</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>In this model, H&#x00027; can be interpreted as a mediator to a dark sector of unknown mass and lifetime. This model could be validated through the experimental signature of the decay <italic>B</italic>&#x02192;<italic>H</italic>&#x02032;<italic>K</italic>, with the H&#x00027; decaying into &#x003C0;<sup>&#x0002B;</sup>&#x003C0;<sup>&#x02212;</sup>, <italic>K</italic><sup>&#x0002B;</sup><italic>K</italic><sup>&#x02212;</sup>, <italic>&#x003BC;</italic><sup>&#x0002B;</sup><italic>&#x003BC;</italic><sup>&#x02212;</sup>, or&#x003C4;<sup>&#x0002B;</sup>&#x003C4;<sup>&#x02212;</sup>, depending on its mass. A displaced vertex could be determined, allowing to reconstruct the H&#x00027; mass from the kinematics and the identification of the two decay particles. The sensitivity to this model depends nevertheless on the H&#x00027; mass and lifetime, which could lead the new scalar to decay outside the detector fiducial volume. In particular, if the H&#x00027; has long lifetime, the two final decay particles would not be selected by the first level of the LHCb trigger, thus escaping detection. <xref ref-type="fig" rid="F8">Figure 8</xref> shows the parameter regions and the sensitivity of different (present and future) experiments to the H&#x00027; (Kachanovich, Aliaksei et al., <xref ref-type="bibr" rid="B9">2020</xref>). One should note that the mixing angle is directly related to the decay width, and thus to the lifetime, with small mixing angles implying long lifetimes.</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p>Parameter space regions for the <italic>B</italic> &#x02192; <italic>H</italic>&#x02032;<italic>K</italic> and <italic>B</italic> &#x02192; <italic>H</italic>&#x02032;<italic>K</italic><sup>&#x0002A;</sup> decay channels <bold>(top)</bold> and sensitivity to the H&#x00027; displaced vertex by different present and future experiments <bold>(bottom)</bold>, (Figures adapted from Kachanovich, Aliaksei et al., <xref ref-type="bibr" rid="B9">2020</xref>).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0008.tif"/>
</fig>
<p>Considering the leptonic decay mode H&#x00027; &#x02192; <italic>&#x003BC;</italic><sup>&#x0002B;</sup><italic>&#x003BC;</italic><sup>&#x02212;</sup>, the decay rate can be expressed by:</p>
<disp-formula id="E2"><label>(2)</label><mml:math id="M5"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mo>&#x00393;</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:mo>&#x02192;</mml:mo><mml:mi>&#x02113;</mml:mi><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mo class="qopname">sin</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mi>&#x003B8;</mml:mi><mml:mfrac><mml:mrow><mml:msub><mml:mrow><mml:mi>G</mml:mi></mml:mrow><mml:mrow><mml:mi>F</mml:mi></mml:mrow></mml:msub><mml:msub><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:msub><mml:msubsup><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mrow><mml:mrow><mml:mn>4</mml:mn><mml:msqrt><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msqrt><mml:mi>&#x003C0;</mml:mi></mml:mrow></mml:mfrac><mml:msup><mml:mrow><mml:mrow><mml:mo stretchy="true">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>-</mml:mo><mml:mfrac><mml:mrow><mml:mn>4</mml:mn><mml:msubsup><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02113;</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mrow><mml:mrow><mml:msubsup><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup></mml:mrow></mml:mfrac></mml:mrow><mml:mo stretchy="true">)</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mn>3</mml:mn><mml:mo>/</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>where <italic>G</italic><sub><italic>F</italic></sub> is the Fermi constant and <italic>m</italic><sub><italic>l</italic></sub> the lepton mass. The Higgs lifetime can be computed:</p>
<disp-formula id="E3"><label>(3)</label><mml:math id="M6"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>&#x003C4;</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mfrac><mml:mrow><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mo>&#x00393;</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:mo>&#x02192;</mml:mo><mml:msup><mml:mrow><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x0002B;</mml:mo></mml:mrow></mml:msup><mml:msup><mml:mrow><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow></mml:mfrac><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>The following sections discuss the sensitivity of the LHCb experiment at the trigger level to the H&#x00027; decay into the <italic>&#x003BC;</italic><sup>&#x0002B;</sup><italic>&#x003BC;</italic><sup>&#x02212;</sup> final state. This decay channel binds the Higgs mass to be <italic>m</italic><sub><italic>H</italic></sub> &#x0003E; 2<italic>m</italic><sub><italic>&#x003BC;</italic></sub> &#x02248; 212 MeV/c<sup>2</sup>. The H&#x00027; mass is also constrained to <inline-formula><mml:math id="M7"><mml:msub><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>H</mml:mi></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup></mml:mrow></mml:msub><mml:mo>&#x0003C;</mml:mo><mml:msub><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>B</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x0002B;</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>m</mml:mi></mml:mrow><mml:mrow><mml:msup><mml:mrow><mml:mi>K</mml:mi></mml:mrow><mml:mrow><mml:mo>&#x0002B;</mml:mo></mml:mrow></mml:msup></mml:mrow></mml:msub><mml:mo>&#x02248;</mml:mo><mml:mn>4700</mml:mn></mml:math></inline-formula> MeV/c<sup>2</sup>.</p>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3. Results</title>
<sec>
<title>3.1. Sensitivity to <italic>B</italic>&#x02192;<italic>H</italic>&#x02032;(&#x02192; <italic>&#x003BC;</italic><sup>&#x0002B;</sup><italic>&#x003BC;</italic><sup>&#x02212;</sup>)<italic>K</italic> decays</title>
<p>Using the upgraded LHCb simulation, 99 MC samples of 7,000 events each are simulated from proton-proton collisions. The samples are generated using Pythia8 and assuming Run 3 beam conditions. The <italic>B</italic>&#x02192;<italic>H</italic>&#x02032;(&#x02192; <italic>&#x003BC;</italic><sup>&#x0002B;</sup><italic>&#x003BC;</italic><sup>&#x02212;</sup>)<italic>K</italic> decay channel is generated considering H&#x00027; masses in the range of 500&#x02013;4,500 MeV and lifetimes from 1 to 2,000 ps.</p>
<p>The decay vertex of the H&#x00027; is expected to be displaced with a dependence on these variables, and they will be labeled in the following according to the track type of the two muons (two Long tracks = LL, two Downstream tracks = DD and two T-tracks = TT). LL vertices are thus expected to be produced in the Velo detector, DD between the Velo and the UT, and TT vertices between the UT and SciFi. The reconstructibility of these vertices is defined according to the track reconstructibility criteria defined in Section 2.3, and imposing that the two muons are coming from the same vertex. <xref ref-type="fig" rid="F9">Figure 9</xref> shows the reconstructibility of the decay vertex of the H&#x00027; particle as a function of its mass and lifetime. For lifetimes below 10 ps, a large proportion of LL vertex topologies is found, as expected, where the H&#x00027; decays in the VELO acceptance and both muons can be reconstructed as Long tracks. Nevertheless, for H&#x00027; lifetimes larger than 100ps (and small-mixing angle), most of the decays are produced downstream the VELO, resulting in a large proportion of the DD and TT topologies. These fractions are very similar in the case that the H&#x00027; decays into two hadrons<xref ref-type="fn" rid="fn0001"><sup>1</sup></xref>.</p>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p>Reconstructibility of the decay vertex of the H&#x00027; particle as a function of its mass and lifetime. Decay topologies are shown, from top to bottom, in the order: LL, DD, and TT, corresponding to the track types of the two muons.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0009.tif"/>
</fig>
<p><xref ref-type="fig" rid="F10">Figure 10</xref> shows the LHCb HLT1 effect when triggering on the H&#x00027; decay products (Trigger on Signal, TOS). Since only Long tracks are reconstructed at the HLT1 level, a high inefficiency can be observed for large H&#x00027; lifetimes, going down 10% for lifetimes larger than 500ps. A loss in sensitivity for small H&#x00027; masses is also observed, and since the H&#x00027; is experiencing larger boosts, muons, are escaping from detection in the VELO.</p>
<fig id="F10" position="float">
<label>Figure 10</label>
<caption><p>Proportion of events triggered by the HLT1 decision on the H&#x00027; decay products (Trigger on Signal, TOS).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0010.tif"/>
</fig>
</sec>
<sec>
<title>3.2. Study of long-lived SM particles: &#x0039B;<sup>0</sup> and <inline-formula><mml:math id="M8"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula></title>
<p>Some strange SM particles, such as <inline-formula><mml:math id="M9"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> and &#x0039B;<sup>0</sup> baryons, also have long lifetimes of the order of 100ps. These particles are reconstructed in many physics analysis to measure observables which can be sensitive to new physics models. One example is the rare <inline-formula><mml:math id="M10"><mml:msub><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02192;</mml:mo><mml:msup><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msup><mml:mi>&#x003B3;</mml:mi></mml:math></inline-formula> decay, where the &#x0039B;<sup>0</sup> baryon decays into a proton and a pion. Measurements of the branching fraction and the angular distribution of the decay particles are sensitive to models with non-standard right-handed currents (Aaij et al., <xref ref-type="bibr" rid="B3">2019</xref>, <xref ref-type="bibr" rid="B1">2022</xref>; Garc&#x00301;&#x00131;a Mart&#x00301;&#x00131;n et al., <xref ref-type="bibr" rid="B8">2019</xref>). Another example is the very rare <inline-formula><mml:math id="M11"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> &#x02192; <italic>&#x003BC;</italic><sup>&#x0002B;</sup><italic>&#x003BC;</italic><sup>&#x02212;</sup> decay, very suppressed in SM, not observed experimentally at present, which is sensitive to different BSM scenarios such as SUSY (Zhu, <xref ref-type="bibr" rid="B16">2012</xref>) or presence of leptoquarks (Bobeth and Buras, <xref ref-type="bibr" rid="B7">2018</xref>). <xref ref-type="fig" rid="F11">Figure 11</xref> shows the sensitivity to the <inline-formula><mml:math id="M12"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> &#x02192; <italic>&#x003BC;</italic><sup>&#x0002B;</sup><italic>&#x003BC;</italic><sup>&#x02212;</sup> decay channel (in terms of branching fraction limit) as a function of the product of the trigger efficiency and luminosity.</p>
<fig id="F11" position="float">
<label>Figure 11</label>
<caption><p>Sensitivity (green band) to the <inline-formula><mml:math id="M13"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> &#x02192; <italic>&#x003BC;</italic><sup>&#x0002B;</sup><italic>&#x003BC;</italic><sup>&#x02212;</sup> branching fraction as function of the trigger efficiency &#x000D7; luminosity. The SM prediction is shown in black (Santos, <xref ref-type="bibr" rid="B15">2022</xref>).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0011.tif"/>
</fig>
<p>Using 10,000 simulated events, the HLT1 trigger effect has been studied on these two decay channels involving &#x0039B;<sup>0</sup> and <inline-formula><mml:math id="M14"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> particles. <xref ref-type="fig" rid="F12">Figure 12</xref> shows the normalized number of reconstructible events (LL&#x0002B;DD&#x0002B;TT) as a function of the end decay vertex of the &#x0039B;<sup>0</sup> and of the <inline-formula><mml:math id="M15"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula>.</p>
<fig id="F12" position="float">
<label>Figure 12</label>
<caption><p>&#x0039B;<sup>0</sup> <bold>(top)</bold> and <inline-formula><mml:math id="M16"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> <bold>(bottom)</bold> reconstructible candidates as function of the decay vertex position. The &#x0039B;<sup>0</sup> candidates decay into a pion and a proton. The <inline-formula><mml:math id="M17"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> candidates decay into two muons of opposite charges. Vertical color lines indicate the positions of the Velo, UT, and SciFi detectors in the Z-axis.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fdata-05-1008737-g0012.tif"/>
</fig>
<p>The relative proportions of track types are 12% (LL), 51% (DD), and 37% (TT) for the <inline-formula><mml:math id="M18"><mml:msub><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02192;</mml:mo><mml:msup><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msup><mml:mi>&#x003B3;</mml:mi></mml:math></inline-formula> decay channel. For the case of prompt <inline-formula><mml:math id="M19"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula>, the amount of reconstructible tracks is 46% (LL), 38% (DD), and 16% (TT), One should notice that a large amount of decays occurs at the end of the VELO and they do not let enough hits to be reconstructed as Long tracks, being reconstructed thus in the DD category. That is, the amount of DD&#x0002B;TT tracks is about 88% for the <inline-formula><mml:math id="M20"><mml:msub><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02192;</mml:mo><mml:msup><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msup><mml:mi>&#x003B3;</mml:mi></mml:math></inline-formula> decay channel and 54% for prompt <inline-formula><mml:math id="M21"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula>.</p>
<p>One can apply the HLT1 conditions on the reconstructible events, and in particular, to check how many events are selected by inclusive trigger lines, such as the OneTrackMVA or TwoTrackMVA. They require tracks with minimum transverse momentum and minimum Impact Parameter (IP) significance with respect to the Primary Vertex, and for the latter, they form a vertex with minimum requirements. In the case of <inline-formula><mml:math id="M22"><mml:msub><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02192;</mml:mo><mml:msup><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msup><mml:mi>&#x003B3;</mml:mi></mml:math></inline-formula>, the HLT1 signal efficiency for the proton and pion coming from the &#x0039B;<sup>0</sup> is found to be less than 10%. In the case of <inline-formula><mml:math id="M23"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula>, the HLT1 efficiency on the muons, adding some inclusive muon lines and dedicated lines for <inline-formula><mml:math id="M24"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> selection, is less than 25%. We normalize the efficiency to the sum of LL&#x0002B;DD&#x0002B;TT reconstructible events. The <inline-formula><mml:math id="M25"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> candidates in this study are prompt, produced at the interaction point. In the case that <inline-formula><mml:math id="M26"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> candidates are coming from the decays of b or c-hadrons, the amount of reconstructible LL candidates is expected to decrease, thus increasing the HLT1 inefficiency.</p>
</sec>
</sec>
<sec id="s4">
<title>4. Prospects</title>
<p>A downstream algorithm, making use of the high performance of the Seeding algorithm described earlier, is currently being developed within the LHCb-RTA Collaboration. Preliminary studies showed that, in the case that DD tracks could be reconstructed and selected at the HLT1 stage, efficiencies for LLPs could increase to 60% or more. A lot of effort is also being devoted to the development of inclusive trigger lines based on SciFi seeds, even if the high occupancy and track combinatorics in this detector makes it quite challenging in terms of allowed rates. The use of machine learning algorithms is also being investigated, which could be applied at an early stage to improve the T-track classification. With the new HLT1 framework based on the GPU architectures, this task seems feasible, now or in a near future, by using faster and better performing GPU cards.</p>
</sec>
<sec sec-type="conclusions" id="s5">
<title>5. Conclusion</title>
<p>In summary, during the past years, the LHCb detector has undergone a major upgrade, which will allow the data for Run 3 to be taken at an averaged proton-proton collision rate of 30 MHz. The first trigger stage of the LHCb will be based for the first time entirely on GPU cards following the Allen framework, a complete high-throughput trigger algorithm design. This will largely increase the capabilities of the trigger system for selecting the data of interest. Tracking reconstruction at early stages of the data chain is becoming an arduous task due to the increase of the output data. This is even more challenging if the reconstruction involves particle trajectories that do not leave signals in the first tracker of the LHCb detector, hitting only the UT and/or the SciFi, which are bigger detectors with high occupancy. Even if the main program of the LHCb experiment is based on the reconstruction of b- and c-hadron decay particles, which are usually Long tracks, there is an important physics case for Downstream and T-tracks to be participating in the first stage of the trigger. The trigger efficiency for long-lived particles is at present very low. This is true for both the cases of BSM and in the Standard Model formalism. In this study, we studied the impact of the HLT1 on a model based on a Higgs portal to a dark sector, where the Higgs can have several masses and lifetimes. For long lifetimes and small masses of the Higgs, the HLT1 inefficiency is found to be quite large. We also tested how the HLT1 decisions affect &#x0039B;<sup>0</sup> and <inline-formula><mml:math id="M27"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> decays, as observed in many Standard Model processes. For these cases, the selected events are one-fourth or less of the possible reconstructible candidates. To improve the HLT1 efficiency, an algorithm is being developed in the Allen framework, which makes the use of the SciFi seeds to reconstruct trajectories of particles with long lifetimes. The inclusion of this algorithm will allow to largely increase the efficiency for &#x0039B;<sup>0</sup> and <inline-formula><mml:math id="M28"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula> and would open a large window for the searches of many particles beyond the Standard Model.</p>
</sec>
<sec id="s6">
<title>6. LHCb-RTA collaboration</title>
<p>The authors of this article want to acknowledge the contribution of Christina Agapopoulou and Lorenzo Pica to the development of the seeding and matching trigger algorithms on GPUs. The contribution of Renato Quagliani has also been decisive for the development of the algorithm. We acknowledge the theoretical discussions with Xabier Cid Vidal. The study of Fernando Martinez Vidal and Nicola Neri, as well as the Universit&#x000E0; di Milano LHCb group on the usability of T-tracks in LHCb, is also very relevant in this context. Similarly, the study of Giovanni Punzi and the LHCb group of the Universit&#x000E0; di Pisa related to other hardware architectures will be crucial for this study in the near future. We also express our gratitude to our colleagues of the LHCb Collaboration, in particular to the LHCb&#x00027;s computing, simulation, DPA, and online projects. This study has been developed in the framework of the RTA project, so we thank our colleagues of RTA for their support and dedication.</p>
</sec>
<sec sec-type="data-availability" id="s7">
<title>Data availability statement</title>
<p>The datasets used in this work have been produced and analyzed inside the LHCb computing framework and are subject to the LHCb agreements and rules. They could be made available upon request, ensuring compliance with privacy agreements and other requirements.</p>
</sec>
<sec id="s8">
<title>Author contributions</title>
<p>LC, LH, AH, BJ, and JZ worked on the Hybrid seeding algorithm and achieved its implementation on GPUs. CVS developed the framework for MC simulations BSM. AO, LH, and DM performed the sensitivity studies including different decay channels and physics models. JZ also contributed to get the track reconstructibility for SM samples. IS worked on the reconstruction of T-tracks from the SciFi seeds and validated their use for physics analyses. AO is the main editor of the paper. All authors have contributed in the edition of the paper and the creation of the figures.</p>
</sec>
<sec sec-type="funding-information" id="s9">
<title>Funding</title>
<p>This work is supported and has been developed in the framework of the RTA project inside the LHCb collaboration. This work has been partially funded by the CERN and the MICINN national agency in Spain. AH is supported by the NSF award 1904160. DM is supported by the CONEXION AIHUB-CSIC. LC is funded by the ERC Refs: 724777-RECEPT and 714536-PRECISION, as well as the PUNCH4NFDI consortium supported by the DFG fund NFDI 39/1.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x00027;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
</body>
<back>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aaij</surname> <given-names>R.</given-names></name> <name><surname>Abdelmotteleb</surname> <given-names>A. S. W.</given-names></name> <name><surname>Beteta</surname> <given-names>C. A.</given-names></name> <name><surname>Abudin&#x000E9;n</surname> <given-names>F. J.</given-names></name> <name><surname>Ackernley</surname> <given-names>T.</given-names></name> <name><surname>Adeva</surname> <given-names>B.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Measurement of the photon polarization in <inline-formula><mml:math id="M29"><mml:msubsup><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msubsup></mml:math></inline-formula> &#x02192; &#x0039B; &#x003B3; decays</article-title>. <source>Phys. Rev. D</source> <volume>105</volume>, <fpage>L051104</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevD.105.L051104</pub-id></citation>
</ref>
<ref id="B2">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aaij</surname> <given-names>R.</given-names></name> <name><surname>Albrecht</surname> <given-names>J.</given-names></name> <name><surname>Belous</surname> <given-names>M.</given-names></name> <name><surname>Billoir</surname> <given-names>P.</given-names></name> <name><surname>Boettcher</surname> <given-names>T.</given-names></name> <name><surname>Rodr&#x000ED;guez</surname> <given-names>A. B.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Allen: a high-level trigger on GPUs for LHCb</article-title>. <source>Comput. Software Big Sci</source>. <volume>4</volume>, <fpage>7</fpage>. <pub-id pub-id-type="doi">10.1007/s41781-020-00039-7</pub-id><pub-id pub-id-type="pmid">33385105</pub-id></citation></ref>
<ref id="B3">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aaij</surname> <given-names>R.</given-names></name> <name><surname>Beteta</surname> <given-names>C. A.</given-names></name> <name><surname>Adeva</surname> <given-names>B.</given-names></name> <name><surname>Adinolfi</surname> <given-names>M.</given-names></name> <name><surname>Aidala</surname> <given-names>C. A.</given-names></name> <name><surname>Ajaltouni</surname> <given-names>Z.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>First observation of the radiative decay <inline-formula><mml:math id="M30"><mml:msubsup><mml:mrow><mml:mo>&#x0039B;</mml:mo></mml:mrow><mml:mrow><mml:mi>b</mml:mi></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow></mml:msubsup><mml:mo>&#x02192;</mml:mo><mml:mo>&#x0039B;</mml:mo><mml:mi>&#x003B3;</mml:mi></mml:math></inline-formula></article-title>. <source>Phys. Rev. Lett</source>. <volume>123</volume>, <fpage>031801</fpage>. <pub-id pub-id-type="doi">10.1103/PhysRevLett.123.031801</pub-id><pub-id pub-id-type="pmid">31386473</pub-id></citation></ref>
<ref id="B4">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Aiola</surname> <given-names>S.</given-names></name> <name><surname>Amhis</surname> <given-names>Y.</given-names></name> <name><surname>Billoir</surname> <given-names>P.</given-names></name> <name><surname>Jashal</surname> <given-names>B. K.</given-names></name> <name><surname>Henry</surname> <given-names>L.</given-names></name> <name><surname>Oyanguren</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Hybrid seeding: a standalone track reconstruction algorithm for scintillating fibre tracker at lhcb</article-title>. <source>Comput. Phys. Commun</source>. <volume>260</volume>, <fpage>107713</fpage>. <pub-id pub-id-type="doi">10.1016/j.cpc.2020.107713</pub-id></citation>
</ref>
<ref id="B5">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Alves Jr</surname> <given-names>A. A.</given-names></name> <name><surname>Barbosa</surname> <given-names>A. F.</given-names></name> <name><surname>Bediaga</surname> <given-names>I.</given-names></name> <name><surname>Cernicchiaro</surname> <given-names>G.</given-names></name> <name><surname>Guerrer</surname> <given-names>G.</given-names></name> <name><surname>Lima Jr</surname> <given-names>H. P.</given-names></name> <etal/></person-group>. (<year>2008</year>). <article-title>The LHCb detector at the LHC</article-title>. <source>JINST</source> <volume>3</volume>, <fpage>S08005</fpage>. <pub-id pub-id-type="doi">10.1088/1748-0221/3/08/S08005</pub-id></citation>
</ref>
<ref id="B6">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bediaga</surname> <given-names>I.</given-names></name> <name><surname>De Miranda</surname> <given-names>J. M.</given-names></name> <name><surname>Rodrigues</surname> <given-names>F.</given-names></name> <etal/></person-group>. (<year>2012</year>). <source>Framework TDR for the LHCb Upgrade: Technical Design Report, CERN-LHCC-2012-007</source>. Technical report.</citation>
</ref>
<ref id="B7">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bobeth</surname> <given-names>C.</given-names></name> <name><surname>Buras</surname> <given-names>A. J.</given-names></name></person-group> (<year>2018</year>). <article-title>Leptoquarks meet &#x003B5;&#x02032;/&#x003B5; and rare Kaon processes</article-title>. <source>JHEP</source> <volume>02</volume>, <fpage>101</fpage>. <pub-id pub-id-type="doi">10.1007/JHEP02(2018)101</pub-id></citation>
</ref>
<ref id="B8">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Garc&#x000ED;a Mart&#x000ED;n</surname> <given-names>L. M.</given-names></name> <name><surname>Jashal</surname> <given-names>B.</given-names></name> <name><surname>Mart&#x000ED;nez Vidal</surname> <given-names>F.</given-names></name> <name><surname>Oyanguren</surname> <given-names>A.</given-names></name> <name><surname>Roy</surname> <given-names>S.</given-names></name> <name><surname>Sain</surname> <given-names>R.</given-names></name> <etal/></person-group>. (<year>2019</year>). <article-title>Radiative <italic>b</italic>-baryon decays to measure the photon and <italic>b</italic>-baryon polarization</article-title>. <source>Eur. Phys. J. C</source> <volume>79</volume>, <fpage>634</fpage>. <pub-id pub-id-type="doi">10.1140/epjc/s10052-019-7123-7</pub-id></citation>
</ref>
<ref id="B9">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kachanovich</surname> <given-names>A.</given-names></name> <name><surname>Nierste</surname> <given-names>U</given-names></name> <name><surname>Nisandzi&#x00107;</surname> <given-names>I.</given-names></name></person-group> (<year>2020</year>). <article-title>Higgs portal to dark matter and <italic>B</italic>&#x02192;<italic>K</italic><sup>(</sup>&#x0002A;) decays</article-title>. <source>Eur. Phys. J. C</source> <volume>80</volume>, <fpage>669</fpage>. <pub-id pub-id-type="doi">10.1140/epjc/s10052-020-8240-z</pub-id></citation>
</ref>
<ref id="B10">
<citation citation-type="journal"><person-group person-group-type="author"><collab>LHCb</collab></person-group> (<year>2014a</year>). <source>LHCb Tracker Upgrade Technical Design Report, LHCB-TDR-015</source>, Technical report.</citation>
</ref>
<ref id="B11">
<citation citation-type="journal"><person-group person-group-type="author"><collab>LHCb</collab></person-group> (<year>2014b</year>). <source>LHCb Trigger and Online Upgrade Technical Design Report, LHCB-TDR-016</source>, Technical report.</citation>
</ref>
<ref id="B12">
<citation citation-type="journal"><person-group person-group-type="author"><collab>LHCb</collab></person-group> (<year>2020</year>). <source>LHCb Upgrade GPU High Level Trigger Technical Design Report, LHCB-TDR-021</source>, Technical report, CERN, Geneva.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal"><person-group person-group-type="author"><collab>LHCb</collab></person-group> (<year>2022</year>). <source>Long-lived particle reconstruction downstream of the LHCb magnet, CERN-LHCb-DP-2022-001</source>, Technical report.</citation>
</ref>
<ref id="B14">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>P.</given-names></name> <name><surname>Rodrigues</surname> <given-names>E.</given-names></name> <name><surname>Stahl</surname> <given-names>S.</given-names></name></person-group> (<year>2021</year>). <article-title>Tracking definitions and conventions for run 3 and beyond</article-title>. <source>LHCb-PUB-2021-B-2005</source>.</citation>
</ref>
<ref id="B15">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Santos</surname> <given-names>D. M.</given-names></name></person-group> (<year>2022</year>). <article-title>Studies performed by Diego Mart&#x000ED;nez Santos (IGFAE-USC), private communication</article-title>. <source>EPJ Web of Conferences.</source> <pub-id pub-id-type="doi">10.1051/epjconf/201817901013</pub-id></citation>
</ref>
<ref id="B16">
<citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname> <given-names>X.</given-names></name></person-group> (<year>2012</year>). <article-title><inline-formula><mml:math id="M31"><mml:mrow><mml:msubsup><mml:mi>K</mml:mi><mml:mtext>S</mml:mtext><mml:mn>0</mml:mn></mml:msubsup></mml:mrow></mml:math></inline-formula>, &#x0039B; and &#x0039E; production at intermediate to high <italic>p</italic><sub><italic>T</italic></sub> from au&#x0002B;au collisions at <inline-formula><mml:math id="M32"><mml:msqrt><mml:mrow><mml:msub><mml:mrow><mml:mi>s</mml:mi></mml:mrow><mml:mrow><mml:mi>N</mml:mi><mml:mi>N</mml:mi></mml:mrow></mml:msub></mml:mrow></mml:msqrt></mml:math></inline-formula>= 39, 11.5 and 7.7 gev</article-title>. <source>Open Phys</source>. <volume>10</volume>, <fpage>1345</fpage>&#x02013;<lpage>1348</lpage>. <pub-id pub-id-type="doi">10.2478/s11534-012-0102-3</pub-id></citation>
</ref>
</ref-list>
<fn-group>
<fn id="fn0001"><p><sup>1</sup>One should note that the tracking reconstruction at this level does not include particle identification.</p></fn>
</fn-group>
</back>
</article>