Search

Abstract We present a detailed mathematical study of the Monte Carlo replica method as applied in the global fitting literature from the high-energy physics theory community. For the first time, we provide a rigorous derivation of the parameter distributions implied by the method, and show that, whilst they agree with Bayesian posteriors for linear models, they disagree otherwise. We proceed to numerically quantify the disagreement between the Monte Carlo replica method and the Bayesian method in the context of two phenomenologically relevant scenarios: fits of the SMEFT Wilson coefficients, and fits of PDFs (albeit in a toy scenario). In both scenarios, we find that uncertainty estimates of the quantities of interest are discrepant between the two approaches when non-linearity is relevant. Our findings motivate future investigation of Bayesian methodologies for global PDF fits, especially in the context of simultaneous determination of PDFs and SMEFT Wilson coefficients.

Abstract Due to the pseudo-scalar nature of the axion-like particle (ALP), the CP-conserving production of two Higgs bosons via the ALP necessarily involves an additional Z or γ boson. We examine the existing constraints from di-Higgs searches at Run 2 of the LHC and find that, despite the presence of extra objects in the final state, these searches are sensitive to a combination of ALP couplings to gluons and three-bosons in the TeV scale range. Additionally, we propose a specialized search strategy incorporating an energetic leptonic Z boson. This refined ALP-induced production process would allow for the identification of the h h → 4 b-jet final state and could potentially probe the TeV scale using data from Run 2 of the LHC. This production process can also occur through a coupling between the top quark and the ALP. We translate the current constraints on di-Higgs production into new limits on the ALP-top coupling.

Abstract We present the open-source SIMUnet code, designed to fit Standard Model Effective Field Theory (SMEFT) Wilson coefficient alongside Parton Distribution Functions (PDFs) of the proton. SIMUnet can perform SMEFT global fits, as well as simultaneous fits of the PDFs and of an arbitrarily large number of SMEFT degrees of freedom, by including both PDF-dependent and PDF-independent observables. SIMUnet can also be used to determine whether the effects of any New Physics models can be fitted away in a global fit of PDFs. SIMUnet is built upon the open-source NNPDF code and is released together with documentation, and tutorials. To illustrate the functionalities of the new tool, we present a new global analysis of the SMEFT Wilson coefficients accounting for their interplay with the PDFs. We increment our previous analysis of the LHC Run II top quark data with both (i) the Higgs production and decay rates data from the LHC, and (ii) the precision electroweak and diboson measurements from LEP and the LHC.

Abstract The interpretation of LHC data, and the assessment of possible hints of new physics, require the precise knowledge of the proton structure in terms of parton distribution functions (PDFs). We present a systematic methodology designed to determine whether and how global PDF fits might inadvertently ‘fit away’ signs of new physics in the high-energy tails of the distributions. We showcase a scenario for the High-Luminosity LHC, in which the PDFs may completely absorb such signs of new physics, thus biasing theoretical predictions and interpretations. We discuss strategies to single out the effects in this scenario, and disentangle the inconsistencies that stem from them. Our study brings to light the synergy between the high luminosity programme at the LHC and future low-energy non-LHC measurements of large-x sea quark distributions. The analysis code used in this work is made public so that any users can test the robustness of the signal associated to a given BSM model against absorption by the PDFs.

Abstract In this paper we explore the coupling of a light axion-like particle (ALP) to top quarks. We use high-energy LHC probes, and examine both the direct probe to this coupling in associated production of a top-pair with an ALP, and the indirect probe through loop-induced gluon fusion to an ALP leading to top pairs. Using the latest LHC Run II data, we provide the best limit on this coupling. We also compare these limits with those obtained from loop-induced couplings in diboson final states, finding that the t t ¯ $$ t\overline{t} $$ +MET channel is the best current handle on this coupling.

Abstract 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 cal- culations, resulting in bounds on 25 Wilson coefficients modifying top quark interactions. Subsequently, we integrate the two analyses within the SIMUnet approach to realise a simul- taneous 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.

Abstract Theoretical interpretations of particle physics data, such as the determination of the Wilson coefficients of the Standard Model Effective Field Theory (SMEFT), often involve the inference of multiple parameters from a global dataset. Optimizing such interpretations requires the identification of observables that exhibit the highest possible sensitivity to the underlying theory parameters. In this work we develop a flexible open source frame-work, ML4EFT, enabling the integration of unbinned multivariate observables into global SMEFT fits. As compared to traditional measurements, such observables enhance the sensitivity to the theory parameters by preventing the information loss incurred when binning in a subset of final-state kinematic variables. Our strategy combines machine learning regression and classification techniques to parameterize high-dimensional likelihood ratios, using the Monte Carlo replica method to estimate and propagate methodological uncertainties. As a proof of concept we construct unbinned multivariate observables for top-quark pair and Higgs+Z production at the LHC, demonstrate their impact on the SMEFT parameter space as compared to binned measurements, and study the improved constraints associated to multivariate inputs. Since the number of neural networks to be trained scales quadratically with the number of parameters and can be fully parallelized, the ML4EFT framework is well-suited to construct unbinned multivariate observables which depend on up to tens of EFT coefficients, as required in global fits.

Abstract We use the Fitmaker tool to incorporate the recent CDF measurement of m W in a global fit to electroweak, Higgs, and diboson data in the Standard Model Effective Field Theory (SMEFT) including dimension-6 operators at linear order. We find that including any one of the SMEFT operators O $$ \mathcal{O} $$ HWB , O $$ \mathcal{O} $$ HD , O $$ \mathcal{O} $$ ℓℓ or O H ℓ 3 $$ {\mathcal{O}}_{H\ell}^{(3)} $$ with a non-zero coefficient could provide a better fit than the Standard Model, with the strongest pull for O $$ \mathcal{O} $$ HD and no tension with other electroweak precision data. We then analyse which tree-level single-field extensions of the Standard Model could generate such operator coefficients with the appropriate sign, and discuss the masses and couplings of these fields that best fit the CDF measurement and other data. In particular, the global fit favours either a singlet Z′ vector boson, a scalar electroweak triplet with zero hypercharge, or a vector electroweak triplet with unit hypercharge, followed by a singlet heavy neutral lepton, all with masses in the multi-TeV range for unit coupling.

Abstract Extensions of the minimal supersymmetric standard model (MSSM) gauge group abound in the literature. Several of these include an additional U(1) X gauge group. Chiral fermions’ charge assignments under U(1) X are constrained to cancel local anomalies in the extension and they determine the structure and phenomenology of it. We provide all anomaly-free charge assignments up to a maximum absolute charge of Q max = 10, assuming that the chiral superfield content of the model is that of the MSSM plus up to three Standard Model (SM) singlet superfields. The fermionic components of these SM singlets may play the rôle of right-handed neutrinos, whereas one of the scalar components may play the rôle of the flavon, spontaneously breaking U(1) X . Easily scanned lists of the charge assignments are made publicly available on Zenodo. For the case where no restriction is placed upon Q max, we also provide an analytic parameterisation of the general solution using simple techniques from algebraic geometry.

Abstract The high-energy tails of charged- and neutral-current Drell-Yan processes provide important constraints on the light quark and anti-quark parton distribution functions (PDFs) in the large-x region. At the same time, short-distance new physics effects such as those encoded by the Standard Model Effective Field Theory (SMEFT) would induce smooth distortions to the same high-energy Drell-Yan tails. In this work, we assess for the first time the interplay between PDFs and EFT effects for high-mass Drell-Yan processes at the LHC and quantify the impact that the consistent joint determination of PDFs and Wilson coefficients has on the bounds derived for the latter. We consider two well-motivated new physics scenarios: 1) electroweak oblique corrections ( W ̂ , Y ̂ $$ \hat{W},\hat{Y} $$ ) and 2) four-fermion interactions potentially related to the LHCb anomalies in R(K (*)). We account for available Drell-Yan data, both from unfolded cross sections and from searches, and carry out dedicated projections for the High-Luminosity LHC. Our main finding is that, while the interplay between PDFs and EFT effects remains moderate for the current dataset, it will become a significant challenge for EFT analyses at the HL-LHC.

Abstract The search for effective field theory deformations of the Standard Model (SM) is a major goal of particle physics that can benefit from a global approach in the framework of the Standard Model Effective Field Theory (SMEFT). For the first time, we include LHC data on top production and differential distributions together with Higgs production and decay rates and Simplified Template Cross-Section (STXS) measurements in a global fit, as well as precision electroweak and diboson measurements from LEP and the LHC, in a global analysis with SMEFT operators of dimension 6 included linearly. We present the constraints on the coefficients of these operators, both individually and when marginalised, in flavour-universal and top-specific scenarios, studying the interplay of these datasets and the correlations they induce in the SMEFT. We then explore the constraints that our linear SMEFT analysis imposes on specific ultra-violet completions of the Standard Model, including those with single additional fields and low-mass stop squarks. We also present a model-independent search for deformations of the SM that contribute to between two and five SMEFT operator coefficients. In no case do we find any significant evidence for physics beyond the SM. Our underlying Fitmaker public code provides a framework for future generalisations of our analysis, including a quadratic treatment of dimension-6 operators.

Abstract We estimate the future sensitivity of the high luminosity (HL-) and high energy (HE-) modes of the Large Hadron Collider (LHC) and of a 100 TeV future circular collider (FCC-hh) to leptoquark (LQ) pair production in the muon-plus-jet decay mode of each LQ. Such LQs are motivated by the fact that they provide an explanation for the neutral current B-anomalies. For each future collider, Standard Model (SM) backgrounds and detector effects are simulated. From these, sensitivities of each collider are found. Our measures of sensitivity are based upon a Run II ATLAS search, which we also use for validation. We illustrate with a narrow scalar (‘$$S_3$$ S3 ’) LQ and find that, in our channel, the HL-LHC has exclusion sensitivity to LQ masses up to 1.8 TeV, the HE-LHC up to 4.8 TeV and the FCC-hh up to 13.5 TeV.

The statistical models used to derive the results of experimental analyses are of incredible scientific value and are essential information for analysis preservation and reuse. In this paper, we make the scientific case for systematically publishing the full statistical models and discuss the technical developments that make this practical. By means of a variety of physics cases -- including parton distribution functions, Higgs boson measurements, effective field theory interpretations, direct searches for new physics, heavy flavor physics, direct dark matter detection, world averages, and beyond the Standard Model global fits -- we illustrate how detailed information on the statistical modelling can enhance the short- and long-term impact of experimental results.