Modern global analyses of the structure of the proton include collider measurements which probe energies well above the electroweak scale. While these provide powerful constraints on the parton distribution functions (PDFs), they are also sensitive to beyond the Standard Model (BSM) dynamics if these affect the fitted distributions. Here we present a first simultaneous determination of the PDFs and BSM effects from deep-inelastic structure function data by means of the NNPDF framework. We consider representative four-fermion operators from the SM Effective Field Theory (SMEFT), quantify to which extent their effects modify the fitted PDFs, and assess how the resulting bounds on the SMEFT degrees of freedom are modified. Our results demonstrate how BSM effects that might otherwise be reabsorbed into the PDFs can be systematically disentangled.

Stefano Carrazza Celine Degrande Shayan Iranipour Juan Rojo Maria Ubiali

Processes involving heavy quarks are a crucial component of the LHC physics program, both by themselves and as backgrounds for Higgs physics and new physics searches. In this work, we critically reconsider the validity of the widely-adopted approximation in which heavy quarks are generated at the matrix-element level, with special emphasis on the impact of the collinear logarithms associated with final-state heavy quark and gluon splittings. Our study, based on a perturbative fragmentation-function approach, explicitly shows that neglecting the resummation of collinear logarithms may yield inaccurate predictions, in particular when observables exclusive in the heavy quark degrees of freedom are considered. Our findings motivate the use of schemes which encompass the resummation of final-state collinear logarithms.

Giovanni Ridolfi Maria Ubiali Marco Zaro

In most extensions of the Standard Model, heavy charged Higgs bosons at the LHC are dominantly produced in association with heavy quarks. An up-to-date determination of the next-to-leading order total cross section in a type-II two-Higgs-doublet model is presented, including a thorough estimate of the theoretical uncertainties due to missing higher-order corrections, parton distribution functions and physical input parameters. Predictions in the four- and five-flavour schemes are compared and reconciled through a recently proposed scale-setting prescription. A four- and five-flavour scheme matched prediction is provided for the interpretation of current and future experimental searches for heavy charged Higgs bosons at the LHC.

Martin Flechl Richard Klees Michael Kramer Michael Spira Maria Ubiali

We present the software framework underlying the NNPDF4.0 global determination of parton distribution functions (PDFs). The code is released under an open source licence and is accompanied by extensive documentation and examples. The code base is composed by a PDF fitting package, tools to handle experimental data and to efficiently compare it to theoretical predictions, and a versatile analysis framework. In addition to ensuring the reproducibility of the NNPDF4.0 (and subsequent) determination, the public release of the NNPDF fitting framework enables a number of phenomenological applications and the production of PDF fits under user-defined data and theory assumptions.

Richard D. Ball Stefano Carrazza Juan Cruz-Martinez Luigi Del Debbio Stefano Forte Tommaso Giani Shayan Iranipour Zahari Kassabov Jose I. Latorre Emanuele R. Nocera Rosalyn L. Pearson Juan Rojo Roy Stegeman Christopher Schwan Maria Ubiali Cameron Voisey Michael Wilson

We compute the perturbative component of the fragmentation function of the $b$ quark to the best of the present theoretical knowledge. The fixed-order calculation to order $\alpha_s^2$ of the fragmentation function at the initial scale is matched with soft-emission logarithm resummation to next-to-next-to-leading logarithmic accuracy, so that order-$\alpha_s^2$ corrections are accounted for exactly, and logarithmically enhanced contributions from loops of $b$ quarks are included. This requires the calculation of the Mellin transform of the order-$\alpha_s^2$ result in the whole complex plane for the Mellin variable, which we provide for the first time for all the fragmenting partons. Evolution is performed to next-to-next- to-leading log accuracy, and mixing with the gluon fragmentation function is taken into account. The perturbative fragmentation functions are made available via LHAPDF grids.

Fabio Maltoni Giovanni Ridolfi Maria Ubiali Marco Zaro

We assess the impact of top quark production at the LHC on global analyses of parton distributions (PDFs) and of Wilson coefficients in the SMEFT, both separately and in the framework of a joint interpretation. We consider the broadest top quark dataset to date containing all available measurements based on the full Run II luminosity. First, we determine the constraints that this dataset provides on the large-x gluon PDF and study its consistency with other gluon-sensitive measurements. Second, we carry out a SMEFT interpretation of the same dataset using state-of-the-art SM and EFT theory calculations, resulting in bounds on 25 Wilson coefficients modifying top quark interactions. Subsequently, we integrate the two analyses within the SIMUnet approach to realise a simultaneous determination of the SMEFT PDFs and the EFT coefficients and identify regions in the parameter space where their interplay is most phenomenologically relevant. We also demonstrate how to separate eventual BSM signals from QCD effects in the interpretation of top quark measurements at the LHC.

Zahari Kassabov Maeve Madigan Luca Mantani James Moore Manuel Morales Alvarado Juan Rojo Maria Ubiali

The unprecedented precision of experimental measurements at the Large Hadron Collider (LHC) and the increased statistics that will be reached in the High-Luminosity phase of the LHC (HL-LHC) are pushing the phenomenology community to a new precision frontier, in which new challenges present themselves and new questions arise. A key ingredients of theoretical predictions at hadron colliders are the Parton Distribution Functions (PDFs) of the proton. This contribution highlights some of the new developments in the determination of PDFs from a global set of experimental data, from approximate N3LO PDFs and the inclusion of theory uncertainties in PDF fits, to the realisation of the non trivial interplay between parton densities at large-x and possible signals of New Physics in high energy tails of the distributions, which highlights the synergy between high energy and low energy experimental programs.

Maria Ubiali

We present a systematic comparison of theoretical predictions and various high-precision experimental measurements, specifically of differential cross sections performed by the LHC run II for Drell-Yan gauge boson, top-quark pair, single-inclusive jet and di-jet production, and by HERA for single-inclusive jet and di-jet production. Theoretical predictions are computed at next-to-next-to-leading order (NNLO) accuracy in perturbative Quantum Chromodynamics. The most widely employed sets of Parton Distribution Functions (PDFs) are used, and PDF, strong coupling, and missing higher order uncertainties are taken into account. We quantitatively assess the predictive power of each PDF set and the contribution of the different sources of experimental and theoretical uncertainty to the agreement between data and predictions. We show that control over all of these aspects is crucial to precision physics studies, such as the determination of Standard Model parameters at the LHC.

Amedeo Chiefa Mark N. Costantini Juan Cruz-Martinez Emanuele R. Nocera Tanjona R. Rabemananjara Juan Rojo Tanishq Sharma Roy Stegeman Maria Ubiali

We present Colibri, an open-source Python code that provides a general and flexible tool for PDF fits. The code is built so that users can implement their own PDF model, and use the built-in functionalities of Colibri for a fast computation of observables. It grants easy access to experimental data, several error propagation methodologies, including the Hessian method, the Monte Carlo replica method, and an efficient numerical Bayesian sampling algorithm. To demonstrate the capabilities of Colibri, we consider its simplest application: a polynomial PDF parametrisation. We perform closure tests using a full set of DIS data and compare the results of Hessian and Monte Carlo fits with those from a Bayesian fit. We further discuss how the functionalities illustrated in this example can be extended to more complex PDF parametrisations. In particular, the Bayesian framework in Colibri provides a principled approach to model selection and model averaging, making it a valuable tool for benchmarking and combining different PDF parametrisations on solid statistical grounds.

Mark N. Costantini Luca Mantani James M. Moore Valentina Schutze Sanchez Maria Ubiali

We show that the well-known divergence of the perturbative expansion of resummed results for processes such as deep-inelastic scattering and Drell-Yan in the soft limit can be treated by Borel resummation. The divergence in the Borel inversion can be removed by the inclusion of suitable higher twist terms. This provides us with an alternative to the standard 'minimal prescription' for the asymptotic summation of the perturbative expansion, and it gives us some handle on the role of higher twist corrections in the soft resummation region.

Stefano Forte Giovanni Ridolfi Joan Rojo Maria Ubiali