This white paper presents the current status of the R&D for Analysis Facilities (AFs) and attempts to summarize the views on the future direction of these facilities. These views have been collected through the High Energy Physics (HEP) Software Foundation's (HSF) Analysis Facilities forum, established in March 2022, the Analysis Ecosystems II workshop, that took place in May 2022, and the WLCG/HSF pre-CHEP workshop, that took place in May 2023. The paper attempts to cover all the aspects of an analysis facility.

High-luminosity particle collider experiments such as the ones planned at the High-Luminosity Large Hadron Collider require ever-greater vertexing precision of the tracking detectors, necessitating reductions in the material budget of the detectors. Traditionally, the fractional radiation length (x/X-0) of detectors is either estimated using known properties of the constituent materials, or measured in dedicated runs of the final detector. In this paper, we present a method of direct measurement of the material budget of a CMS prototype module designed for the Phase-2 upgrade of the CMS detector using a 40-65 MeV positron beam. A total of 630 million events were collected at the Paul Scherrer Institut PiE1 experimental area using a three-plane telescope consisting of the prototype module as the central plane, surrounded by two MALTA monolithic pixel detectors. Fractional radiation lengths were extracted from scattering angle distributions using the Highland approximation for multiple scattering. A statistical technique recovered runs suffering from trigger desynchronisation, and several corrections were introduced to compensate for local inefficiencies related to geometric and beam shape constraints. Two regions of the module were surveyed and yielded average x/X-0 values of (0.72 +/- 0.05)% and (0.95 +/- 0.09)%, which are compatible with empirical estimates for these regions computed from known material properties of 0.753% and 0.892%, respectively. Two types of higher-granularity maps of the fractional radiation length were produced, subdivided either into rectangular regions of uniform size, or polygonal-shaped regions of uniform material composition. The results bode well for the CMS Phase-2 upgrade modules, which will play a key role in the minimisation of the material of the upgraded detector.

The Large Hadron Collider at CERN will undergo an upgrade in order to increase its luminosity to 7.5 × 10 34 cm -2 s -1 . The increased luminosity during this High-Luminosity running phase, starting around 2029, means a higher rate of proton-proton interactions, hence a larger ionizing dose and particle fluence for the detectors. The current tracking system of the CMS experiment will be fully replaced in order to cope with the new operating conditions. Prototype planar pixel sensors for the CMS Inner Tracker with square 50 μm × 50 μm and rectangular 100 μm × 25 μm pixels read out by the RD53A chip were characterized in the lab and at the DESY-II testbeam facility in order to identify designs that meet the requirements of CMS during the High-Luminosity running phase. A spatial resolution of approximately 3.4 μm (2 μm) is obtained using the modules with 50 μm × 50 μm (100 μm × 25 μm) pixels at the optimal angle of incidence before irradiation. After irradiation to a 1 MeV neutron equivalent fluence of Φ eq = 5.3 × 10 15 cm -2 , a resolution of 9.4 μm is achieved at a bias voltage of 800 V using a module with 50 μm × 50 μm pixel size. All modules retain a hit efficiency in excess of 99% after irradiation to fluences up to 2.1 × 10 16 cm -2 . Further studies of the electrical properties of the modules, especially crosstalk, are also presented in this paper.

The Large Hadron Collider (LHC) at CERN will undergo major upgrades to increase the instantaneous luminosity up to 5–7.5×10 34 cm -2 s -1 . This High Luminosity upgrade of the LHC (HL-LHC) will deliver a total of 3000–4000 fb -1 of proton-proton collisions at a center-of-mass energy of 13–14 TeV. To cope with these challenging environmental conditions, the strip tracker of the CMS experiment will be upgraded using modules with two closely-spaced silicon sensors to provide information to include tracking in the Level-1 trigger selection. This paper describes the performance, in a test beam experiment, of the first prototype module based on the final version of the CMS Binary Chip front-end ASIC before and after the module was irradiated with neutrons. Results demonstrate that the prototype module satisfies the requirements, providing efficient tracking information, after being irradiated with a total fluence comparable to the one expected through the lifetime of the experiment.

To cope with the challenging environment of the planned high luminosity upgrade of the Large Hadron Collider (HL-LHC), scheduled to start operation in 2029, CMS will replace its entire tracking system. The requirements for the tracker are largely determined by the long operation time of 10 years with an instantaneous peak luminosity of up to 7.5 x 1034 cm-2 s-1 in the ultimate performance scenario. Depending on the radial distance from the interaction point, the silicon sensors will receive a particle fluence corresponding to a non-ionising energy loss of up to ?eq = 3.5 x 1016 cm-2. This paper focuses on planar pixel sensor design and qualification up to a fluence of ?eq = 1.4 x 1016 cm-2. For the development of appropriate planar pixel sensors an R&D program was initiated, which includes n+-p sensors on 150 mm (6") wafers with an active thickness of 150 mu m with pixel sizes of 100 x 25 mu m2 and 50 x 50 mu m2 manufactured by Hamamatsu Photonics K.K. (HPK). Single chip modules with ROC4Sens and RD53A readout chips were made. Irradiation with protons and neutrons, as well was an extensive test beam campaign at DESY were carried out. This paper presents the investigation of various assemblies mainly with ROC4Sens readout chips. It demonstrates that multiple designs fulfil the requirements in terms of breakdown voltage, leakage current and efficiency. The single point resolution for 50 x 50 mu m2 pixels is measured as 4.0 mu m for non-irradiated samples, and 6.3 mu m after irradiation to ?eq = 7.2 x 1015 cm-2.

The Short Strip ASIC (SSA) is one of the four front-end chips designed for the upgrade of the CMS Outer Tracker for the High Luminosity LHC. Together with the Macro-Pixel ASIC (MPA) it will instrument modules containing a strip and a macro-pixel sensor stacked on top of each other. The SSA provides both full readout of the strip hit information when triggered, and, together with the MPA, correlated clusters called stubs from the two sensors for use by the CMS Level-1 (L1) trigger system. Results from the first prototype module consisting of a sensor and two SSA chips are presented. The prototype module has been characterized at the Fermilab Test Beam Facility using a 120 GeV proton beam.

Using a data sample of $$\sqrt{s}=13\,\text {TeV}$$ proton-proton collisions collected by the CMS experiment at the LHC in 2017 and 2018 with an integrated luminosity of $$103\text {~fb}^{-1}$$ , the $$\text {B}^{0}_{\mathrm{s}} \rightarrow \uppsi (\text {2S})\text {K}_\mathrm{S}^{0}$$ and $$\text {B}^{0} \rightarrow \uppsi (\text {2S})\text {K}_\mathrm{S}^{0} \uppi ^+\uppi ^-$$ decays are observed with significances exceeding 5 standard deviations. The resulting branching fraction ratios, measured for the first time, correspond to $${\mathcal {B}}(\text {B}^{0}_{\mathrm{s}} \rightarrow \uppsi (\text {2S})K_\mathrm{S}^{0})/{\mathcal {B}}(\text {B}^{0}\rightarrow \uppsi (\text {2S})K_\mathrm{S}^{0}) = (3.33 \pm 0.69 (\text {stat})\, \pm 0.11\,(\text {syst}) \pm 0.34\,(f_{\mathrm{s}}/f_{\mathrm{d}})) \times 10^{-2}$$ and $${\mathcal {B}}(\text {B}^{0} \rightarrow \uppsi (\text {2S})\text {K}_\mathrm{S}^{0} \uppi ^{+} \uppi ^{-})/ {\mathcal {B}}(\text {B}^{0} \rightarrow \uppsi (\text {2S})\text {K}^{0}_{\mathrm{S}}) = 0.480 \pm 0.013\,(\text {stat}) \pm 0.032\,(\text {syst})$$ , where the last uncertainty in the first ratio is related to the uncertainty in the ratio of production cross sections of $$\hbox {B}^{0}_{\mathrm{s}}$$ and $$\hbox {B}^{0}$$ mesons, $$f_{\mathrm{s}}/f_{\mathrm{d}}$$ .

With the ever growing amount of data collected by the ATLAS and CMS experiments at the CERN LHC, fiducial and differential measurements of the Higgs boson production cross section have become important tools to test the standard model predictions with an unprecedented level of precision, as well as seeking deviations that can manifest the presence of physics beyond the standard model. These measurements are in general designed for being easily comparable to any present or future theoretical prediction, and to achieve this goal it is important to keep the model dependence to a minimum. Nevertheless, the reduction of the model dependence usually comes at the expense of the measurement precision, preventing to exploit the full potential of the signal extraction procedure. In this paper a novel methodology based on the machine learning concept of domain adaptation is proposed, which allows using a complex deep neural network in the signal extraction procedure while ensuring a minimal dependence of the measurements on the theoretical modelling of the signal.

The CMS Inner Tracker, made of silicon pixel modules, will be entirely replaced prior to the start of the High Luminosity LHC period. One of the crucial components of the new Inner Tracker system is the readout chip, being developed by the RD53 Collaboration, and in particular its analogue front-end, which receives the signal from the sensor and digitizes it. Three different analogue front-ends (Synchronous, Linear, and Differential) were designed and implemented in the RD53A demonstrator chip. A dedicated evaluation program was carried out to select the most suitable design to build a radiation tolerant pixel detector able to sustain high particle rates with high efficiency and a small fraction of spurious pixel hits. The test results showed that all three analogue front-ends presented strong points, but also limitations. The Differential front-end demonstrated very low noise, but the threshold tuning became problematic after irradiation. Moreover, a saturation in the preamplifier feedback loop affected the return of the signal to baseline and thus increased the dead time. The Synchronous front-end showed very good timing performance, but also higher noise. For the Linear front-end all of the parameters were within specification, although this design had the largest time walk. This limitation was addressed and mitigated in an improved design. The analysis of the advantages and disadvantages of the three front-ends in the context of the CMS Inner Tracker operation requirements led to the selection of the improved design Linear front-end for integration in the final CMS readout chip.

A combination of searches for top squark pair production using proton-proton collision data at a center-of-mass energy of 13 TeV at the CERN LHC, corresponding to an integrated luminosity of 137 fb(-1) collected by the CMS experiment, is presented. Signatures with at least 2 jets and large missing transverse momentum are categorized into events with 0, 1, or 2 leptons. New results for regions of parameter space where the kinematical properties of top squark pair production and top quark pair production are very similar are presented. Depending on themodel, the combined result excludes a top squarkmass up to 1325 GeV for amassless neutralino, and a neutralinomass up to 700 GeV for a top squarkmass of 1150 GeV. Top squarks with masses from 145 to 295 GeV, for neutralino masses from 0 to 100 GeV, with a mass difference between the top squark and the neutralino in a window of 30 GeV around the mass of the top quark, are excluded for the first time with CMS data. The results of theses searches are also interpreted in an alternative signal model of dark matter production via a spin-0 mediator in association with a top quark pair. Upper limits are set on the cross section for mediator particle masses of up to 420 GeV.

Abstract During the operation of the CMS experiment at the High-Luminosity LHC the silicon sensors of the Phase-2 Outer Tracker will be exposed to radiation levels that could potentially deteriorate their performance. Previous studies had determined that planar float zone silicon with n-doped strips on a p-doped substrate was preferred over p-doped strips on an n-doped substrate. The last step in evaluating the optimal design for the mass production of about 200 m2 of silicon sensors was to compare sensors of baseline thickness (about 300 μm) to thinned sensors (about 240 μm), which promised several benefits at high radiation levels because of the higher electric fields at the same bias voltage. This study provides a direct comparison of these two thicknesses in terms of sensor characteristics as well as charge collection and hit efficiency for fluences up to 1.5 × 1015 neq/cm2. The measurement results demonstrate that sensors with about 300 μm thickness will ensure excellent tracking performance even at the highest considered fluence levels expected for the Phase-2 Outer Tracker.

The second-order Fourier coefficients (v(2)) characterizing the azimuthal distributions of Y(1S) and Y(2S) mesons produced in PbPb collisions at root s(NN) = 5.02 TeV are studied. The Y mesons are reconstructed in their dimuon decay channel, as measured by the CMS detector. The collected data set corresponds to an integrated luminosity of 1.7 nb(-1). The scalar product method is used to extract the v2 coefficients of the azimuthal distributions. Results are reported for the rapidity range vertical bar y vertical bar < 2.4, in the transverse momentum interval 0 < pT < 50 GeV/c, and in three centrality ranges of 10-30%, 30-50% and 50-90%. In contrast to the J/psi mesons, the measured v(2) values for the Y mesons are found to be consistent with zero. (C) 2021 The Author(s). Published by Elsevier B.V.

The first measurement of the dependence of $\\gamma\\gamma \\to \\mu^{+}\\mu^{-}$ production on the multiplicity of neutrons emitted very close to the beam direction in ultraperipheral heavy ion collisions is reported. Data for lead-lead interactions at $\\sqrt{s_\\mathrm{NN}} =$ 5.02 TeV, with an integrated luminosity of approximately 1.5 nb$^{-1}$, were collected using the CMS detector at the LHC. The azimuthal correlations between the two muons in the invariant mass region 8 $\\lt m_{\\mu\\mu} \\lt$ 60 GeV are extracted for events including zero, one, or at least two neutrons detected in the forward pseudorapidity range $|\\eta| \\lt$ 8.3. The back-to-back correlation structure from leading-order photon-photon scattering is found to be significantly broader for events with a larger number of emitted neutrons from each nucleus, corresponding to interactions with a smaller impact parameter. This observation provides a data-driven demonstration that the average transverse momentum of photons emitted from relativistic heavy ions has an impact parameter dependence. These results provide new constraints on models of photon-induced interactions in ultraperipheral collisions. They also provide a baseline to search for possible final-state effects on lepton pairs caused by traversing a quark-gluon plasma produced in hadronic heavy ion collisions.

A measurement is presented of electroweak (EW) production of a W boson in association with two jets in proton-proton collisions at s = 13 Te . The data sample was recorded by the CMS Collaboration at the LHC and corresponds to an integrated luminosity of 35.9 fb - 1 . The measurement is performed for the l ν jj final state (with l ν indicating a lepton-neutrino pair, and j representing the quarks produced in the hard interaction) in a kinematic region defined by invariant mass m jj \u003e 120 Ge and transverse momenta p T j \u003e 25 Ge . The cross section of the process is measured in the electron and muon channels yielding σ EW ( W jj ) = 6.23 ± 0.12 (stat) ± 0.61 (syst) pb per channel, in agreement with leading-order standard model predictions. The additional hadronic activity of events in a signal-enriched region is studied, and the measurements are compared with predictions. The final state is also used to perform a search for anomalous trilinear gauge couplings. Limits on anomalous trilinear gauge couplings associated with dimension-six operators are given in the framework of an effective field theory. The corresponding 95% confidence level intervals are - 2.3 c W W W / Λ 2 2.5 Te - 2 , - 8.8 c W / Λ 2 16 Te - 2 , and - 45 c B / Λ 2 46 Te - 2 . These results are combined with the CMS EW Zjj analysis, yielding the constraint on the c W W W coupling: - 1.8 c W W W / Λ 2 2.0 Te - 2 .

A search is conducted for a low-mass charged Higgs boson produced in a top quark decay and subsequently decaying into a charm and a strange quark. The data sample was recorded in proton-proton collisions at $\\sqrt{s}=$ 13 TeV by the CMS experiment at the LHC and corresponds to an integrated luminosity of 35.9 fb$^{-1}$. The search is performed in the process of top quark pair production, where one top quark decays to a bottom quark and a charged Higgs boson, and the other to a bottom quark and a W boson. With the W boson decaying to a charged lepton (electron or muon) and a neutrino, the final state comprises an isolated lepton, missing transverse momentum, and at least four jets, of which two are tagged as b jets. To enhance the search sensitivity, one of the jets originating from the charged Higgs boson is required to satisfy a charm tagging selection. No significant excess beyond standard model predictions is found in the dijet invariant mass distribution. An upper limit in the range 1.68-0.25% is set on the branching fraction of the top quark decay to the charged Higgs boson and bottom quark for a charged Higgs boson mass between 80 and 160 GeV.

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective anisotropic flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in collisions at nucleon-nucleon center-of-mass energies and 5.02 with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.

Events with no charged particles produced between the two leading jets are studied in proton-proton collisions at

Anisotropies in the initial energy density distribution of the quark-gluon plasma created in high energy heavy ion collisions lead to anisotropies in the azimuthal distributions of the final-state particles known as collective anisotropic flow. Fourier harmonic decomposition is used to quantify these anisotropies. The higher-order harmonics can be induced by the same order anisotropies (linear response) or by the combined influence of several lower order anisotropies (nonlinear response) in the initial state. The mixed higher-order anisotropic flow and nonlinear response coefficients of charged particles are measured as functions of transverse momentum and centrality in PbPb collisions at nucleon-nucleon center-of-mass energies $\\sqrt{S}$$_{NN}$=2.76 and 5.02TeV with the CMS detector. The results are compared with viscous hydrodynamic calculations using several different initial conditions, as well as microscopic transport model calculations. None of the models provides a simultaneous description of the mixed higher-order flow harmonics and nonlinear response coefficients.