Abstract

We study the impact of Run 2 LHC data on general composite Higgs scenarios, where nonlinear effects, mixing with additional scalars, and new fermionic degrees of freedom could simultaneously contribute to the modification of Higgs properties. We obtain new experimental limits on the scale of compositeness, the mixing with singlets and doublets with the Higgs, and the mass and mixing angle of top-partners. We also show that for scenarios where new fermionic degrees of freedom are involved in electroweak symmetry breaking, there is an interesting interplay among Higgs coupling measurements, boosted Higgs properties, SMEFT global analyses, and direct searches for single and double production of vector-like quarks.

1. Introduction

The true origin of the Higgs mechanism is still an open question in Particle Physics, despite the discovery of its key element, the Higgs particle [1, 2], and the observation of the SM-like nature of its couplings to massive particles [3].

The main reason to doubt a purely SM Higgs sector can be traced back its quantum behaviour, a problem often expressed in the context of naturalness: how can a fundamental scalar be so light, yet so sensitive to UV effects. And this is not the only suspicious aspect of the SM Higgs. The hierarchy among its Yukawa couplings or its inability to produce enough baryon asymmetry during the electroweak phase transition also add to the unsatisfactory aspects of the SM Higgs sector.

However, these shortcomings also open opportunities for new physics. The Higgs, due to its scalar nature, could connect to other sectors via mixing with scalars and participate in phase transitions and inflation. Among the theories bringing the Higgs into a new framework, one particularly appealing proposal is the concept of compositeness, leading to composite Higgs models (CHMs). In CHMs, the Higgs is not a truly fundamental particle but a bound state of other, more fundamental, particles.

Primitive proposals for a composite Higgs [4, 5] were replaced by more realistic realisations based on the idea that a Higgs is not just any composite state from strong dynamics, but a pseudo-Goldstone boson [6, 7]. Pseudo-Goldstone bosons appear when approximate global symmetries are spontaneously broken. The size of the symmetry and what it breaks down to determines the amount of scalar degrees of freedom one will have in CHMs. The partial gauging of some of these global symmetries sets what kinds of quantum numbers the scalars have. The minimal framework to achieve successful electroweak symmetry breaking (EWSB), i.e., leading to a Higgs doublet of with a , was found in Ref. [4, 5]. In this setup, the Higgs originated from the breaking , and its radiative potential had to be supplemented with a new vector-like fermion with the same quantum numbers as the top, called top-partner and denoted by . This minimal set-up and its phenomenological consequences have been thoroughly explored in the literature, see, e.g., Ref. [8] for a review.

Typically, one would assume the main manifestation of the composite nature of the Higgs would be the nonlinear origin of its couplings, leading to very specific types of Higgs coupling deviations. Yet, composite Higgs scenarios could exhibit a richer phenomenology, such as the presence of new scalars or fermions which would also modify the Higgs properties. The focus of this paper is the study of the rich patterns which arise in general composite Higgs scenarios and what the Run 2 LHC data can tell us about them. The data we will use for this analysis is summarised in Table 1.

The paper is organised as follows. In Section 2, we describe the different patterns arising in composite Higgs scenarios: nonlinear effects (Section 2.1), extended scalar sectors (Section 2.2), and new fermionic degrees of freedom (Section 2.3). In these sections, we examine the impact that Higgs measurements and direct top-partner searches have on these patterns and in combinations of them, e.g., how nonlinearities and mixing with new scalars do contribute on the same direction to reduce Higgs couplings to massive particles. Hence, as one switches on more than one of these effects, each individual limit becomes stronger. Section 3 is devoted to conclusions.

2. Patterns of Composite Higgs Scenarios

The idea of the Higgs as a composite state could be realised in many ways. Considering that theoretical constraints and precision measurements substantially reduce the possible scenarios, yet there are still many possibilities to contemplate. Nevertheless, one can describe modern composite Higgs scenarios as exhibiting one or more of the following patterns: (1)A composite Higgs sector would be able to trigger electroweak symmetry in a nonlinear fashion. This leads to new possibilities for the light Higgs boson. Its couplings, and even quantum numbers, could be different from the SM expectation. Indeed, in the SM the Higgs is a -doublet, and the way it induces EWSB completely determines how it couples to fermions and bosons, including couplings to more than one Higgs. On the other hand, in CHMs, these couplings would be more general, even opening the possibility of a singlet Higgs [22, 23]. This additional freedom leads to the expectation that the SM-like Higgs, the 125 GeV state observed by the CERN collaborations, should exhibit nonstandard interactions. We will discuss this first composite Higgs pattern in Section 2.1(2)The second set of patterns found in CHMs refers to the origin of the composite Higgs in the spontaneous breaking of an approximate global symmetry. This breaking could lead to just the right amount of degrees of freedom to match one single light Higgs particle and the three missing polarisations of the and bosons. But, more generically, this breaking would lead to new scalar degrees of freedom which, after partially gauging the global symmetry, could be identified as singlets, doublets, or even higher representations. These nonminimal coset embeddings would lead to additional degrees of freedom which could also develop their own potential, participate in EWSB, and typically mix with the SM-like Higgs. As we have not observed more scalars than the Higgs, these additional scalars would be typically required to be heavier than the Higgs. Therefore, their indirect effects could appear through virtual corrections to the Higgs coupings. At leading order, the effects could be at tree level in the effective theory via mixing with the SM-like Higgs. The study of the Run 2 limits on these extensions of the SM scalar sector can be found on Section 2.2. Note that the additional pseudo-Goldstone degrees of freedom would exhibit their own phenomenology in the form of the so-called axion-like particles (ALPs), see Ref. [24] for a brief introduction to ALPS. The same parameter which we will constrain using the Higgs couplings would also influence the ALP phenomenology, see, e.g., Refs. [2530](3)The third set of patterns found in CHMs refers to the presence of new, fermionic degrees of freedom which would assist the SM-like Higgs in its job of triggering EWSB. These fermions, often called top-partners, should be relatively light to efficiently modify the Higgs potential. They can be searched for directly, via their production at the LHC through couplings to the strong or electroweak sector, or indirectly, as they modify the Higgs properties. The discussion on direct and indirect searches for top-partners will be developed in Section 2.3(4)In each section, we also discuss situations where more than one dominant pattern could be at play in the LHC measurements. We will explain how in this case, their effect in the Higgs couplings would go in the same direction, ruling out cancellations which could invalidate the limits from the previous sections

These patterns are summarised in Figure 1.

2.1. CHM Pattern I: Nonlinear Effects

In CHM, one assumes the existence of a global symmetry , spontaneously broken to a smaller subgroup . The Higgs particle and the would-be Goldstone bosons for the and originate from this breaking. A part of is then weakly gauged to provide the Higgs and massive gauge bosons with the correct properties under the SM.

To go beyond the SM structure of EWSB, it is useful to follow the Callan, Coleman, Wess, and Zumino prescription [31]. Within this prescription, one aims to build objects with definite transformation properties under . In CHMs, the object that contains the Goldstones from the global breaking is , where is the Goldstones and generators in the coset and is the scale associated with the spontaneous breaking. At low energies, below , one can write an effective Lagrangian involving

which, after making a choice for weakly gauging , determines the Higgs and vector boson interactions. Assuming the Higgs is a doublet of , similarly as in the SM, the structure that arises from this kinetic term is

One could also allow the Higgs to be a singlet under , but that would lead to the need to impose a tuning in the effective Lagrangian to explain the specific relations between the and masses [32, 33]. We will not follow this path in this paper and embed the CHM global symmetry structure in a way that preserves the custodial symmetry present in the SM.

Expanding around the physical Higgs boson vacuum expectation value (VEV), one finds clear predictions of the structure of the Higgs-vector boson couplings where .

Turning now to the Higgs-fermion couplings, we note that the predictions are not unique. In particular, the fermion mass generation mechanism in CHMs usually relies on the concept of partial compositeness [32, 34], i.e., the idea that SM fermions feel EWSB via mixing with fermionic bound states. To build structures leading to fermion couplings to the Higgs, one needs to specify the embedding of the fermionic degrees of freedom in the global symmetry structure . As an example, in the minimal composite Higgs models, based on coset group, will transform as 5-plet of which could form a invariant either with two 5-plets or pair of 5 and 10-plets.

More general CHMs lead to very different patterns of breaking and types of embeddings for the fermion content of the model. Yes, despite spanning many model building options, it was noted in Ref. [35] that the Yukawa couplings usually fall into two choices,

For example, the models based on coset groups [7, 36], [3739], [40], and [41, 42] have fermions-Higgs couplings modified by . On the other hand, -type couplings could exist in [36, 4345], [46], [47], and [3739] group-based models. In all these models, the Higgs doublet lies within the unbroken subgroup, but in some of those, there could be an extra singlet or Higgs doublet.

A composite Higgs, with nonstandard couplings to massive fermions and bosons, would also exhibit nonstandard loop-level couplings to gluons and photons () [9, 48].

as well as deviations on the Higgs width ()

Note that and do not have to be equal (see Equations (4) and (5)) depending on how fermions are embedded within the global symmetry group.

Putting all this together, we are ready to compare measurements of the Higgs properties with expectations from a nonlinear realisation of EWSB via a composite Higgs. The experimental inputs, described in Tables 2 and 3, are Higgs signal strengths for different production and decay channels. We compare these measurements with theoretical predictions of the signal strength as a function of the modifiers , , , , and . For example, for the gluon fusion () channel, it adopts the form: . Note that all ’s are a function of a single parameter model, , as shown in Equations (3)–(5).

To compare the CHM predictions with Higgs signal strength measurements from CMS and ATLAS experiments, we evaluate the statistic test

Here, , , and denote the model prediction of the signal strength, experimental measurement, and error for the experimental measurement, respectively. The index runs over all the measurements from Run 1 and Run 2. Note that the correlations among the different experimental measurements are not considered as they are subdominant with respect to the model uncertainties introduced by the choices of .

In this section, we only focus on nonlinear effects; hence, the rates are sensitive to one parameter, the ratio, and choices for the fermion couplings from Equations (4) and (5). In Figure 2, we show the fit of different choices for to the combined LHC Run 1 and Run 2 data. Run 1 data is taken from the combined CMS and ATLAS analysis [9]. For Run 2 data, individual measurements from CMS [1216] (see Table 2) and ATLAS [10, 11] (see Table 3) are considered. For decay channel, both CMS and ATLAS have combined Run 1 and Run 2 data measurements so we have not considered the single Run 1 measurements.

In Figure 2, we show the 1- and 2-sigma limits on in green and yellow, respectively. The value is approx 46 with 59 observables, which leads to approximately 0.78. The values of depend on whether the top, bottom, and tau modifiers are of type or . Following Ref. [35], we plot all the lines with all combinations. We also denote which combination is less constrained, corresponding to , , and , and denoted by the symbol ABB. We also denote which combination of coupling choices is constrained more strongly by LHC data, BAA in the same notation. Note that the less constrained option, ABB, corresponds to a smaller deviation of top couplings than bottom and tau, see Equations (4) and (5).

Comparing with Ref. [35], we obtain a stronger lower bound on the scale of compositeness of  GeV ( is a minimum scale at which the considered fermion coupling combination is compatible with the data at 2 level.). The different curves highlight the fact that there is roughly 500 GeV variation in the scale, reaching  TeV for the most tightly bound scenario.

2.2. CHM Pattern II: Extended Higgs Sectors

When the group breaks down to , light degrees of freedom are generated. Only in specific configurations one would expect exactly four Goldstones, to match the SM needs. Hence, generic composite Higgs scenarios would typically exhibit a pattern of extended Higgs sectors, with more light scalars involved in EWSB. In the simplest of such extensions, one would consider an additional singlet, see, e.g., [49], which could mix with the SM-like Higgs doublet. The next simplest iteration would introduce new doublets, as is the case in, e.g., the coset group [39], where the light degrees of freedom organise as a two Higgs doublet model (2HDMs) [50].

Doublets and singlets are not the only model-building possibilities which composite Higgs models offer, e.g., on the coset, one has two Higgs doublets but also a custodial bitriplet [51]. Generically speaking, a composite Higgs scenario leading to one single doublet at low energies is minimal but not typical, and one should consider the effect of more scalars involved in EWSB.

Moreover, this pattern opens new, interesting possibilities for model building beyond EWSB. These extra singlet or doublet pseudo-Goldstones could play a role in Inflation [52, 53] or act as dark matter candidates [51]. In composite Higgs scenarios, these new scalars could help on enhancing the strength of the electroweak phase transition [54] by introducing new sources of CP violation and more interesting phase diagram structures, and even play some role in the QCD CP problem [55].

In this section, we explore what the Run 2 LHC data can tell us about these extensions. We will discuss the modification of the gauge boson and fermions couplings to the SM-like Higgs boson when the additional scalars mix with the SM-like Higgs. Their effect could also be felt at loop level, even in the absence of a mixing, but these contributions would be suppressed by loop factors respect to mixing. Both possibilities, mixing and loop effects, were computed and explored in Ref. [56], and here, we make use of these theoretical calculations and update the limits in the context of CHMs and the possible interplay between mixing and nonlinearities.

2.2.1. Singlet Scalar

The presence of an extra singlet modifies the by a factor of , where denotes the mixing angle between the neutral scalar and the extra singlet scalar field. This effect is simply due to linear mixing terms when both the singlet and the SM-like Higgs get their VEVs. Nonlinearities due to the origin of the SM-like Higgs as a composite Higgs would still be present, in exactly the same way we discussed in Section 2.1 (One could consider situations where the additional singlets or doublets do participate directly in the mechanism for EWSB, as shown in the seesaw composite Higgs [57], where this assumption would fail). Therefore, the Higgs couplings to vector bosons would be doubly modified as where we have expanded for small values of nonlinearities and mixing to show how these two effects work in the same direction, i.e., to reduce the coupling value from the SM expectation.

The nonlinear part of the modification of the fermion couplings depends on the fermion embedding in representations of the global symmetry. As discussed before, we find two main choices for , namely, or where again we expand for small modifications and to show explicitly the cooperative effort of both effects to lower the coupling.

We perform the fit to the Higgs data as in Section 2.1 but now including the effect of mixing, . As shown in Figure 2, we vary options for the fermion couplings and find the . We perform this fit for different values of the mixing angle, always finding the minimum corresponding to a set of options for . The result of this procedure is shown in Figure 3, where we plot scale as a function of the singlet mixing. On the right side of the plot, we see the effect of small mixing, where we recover the result from the previous section. As we move towards the left in the plot, the mixing becomes more important and at some point , we obtain the pure mixing limit of 0.96.

Note that the mixing angle is determined by singlet mass (), singlet vev (), and quadratic coupling of light Higgs doublet and singlet scalar (), as given below where  GeV. Therefore, the limit 0.96 would be equivalent to 450 GeV.

In between these two asymptotic limits, the constraints on the scales and become stronger, as both mixing and nonlinearities work together to reduce the Higgs couplings.

2.2.2. Two Higgs Doublet Models (2HDM)

As discussed in the previous section, the mixing effect due the additional Higgs in a composite Higgs model does not typically couple with the nonlinear effects at leading order. Hence, in composite two Higgs doublet models, the vector and femion couplings to the Higgs, , would still adopt this factorisable form , where has been discussed in Section 2.1 and correspond to the modifications of the SM-like Higgs couplings in 2HDM models.

The explicit form of is given in Table 4 for various types of 2HDM models. In these scenarios, vector and fermion modifiers are determined by two parameters, the ratio of symmetry breaking VEVs of Higgs doublets (), and neutral Higgs mixing angle ().

In our analysis, we only consider the main scenarios, Type I and Type II. The statistical test for a composite Higgs with possible mixing with another doublet becomes a function of three parameters: , , and .

First, we discuss the information that Run 2 LHC data provides on type I and II 2HDMs, without considering the nonlinear effects. This information is shown in Figure 4, where the coloured regions green, yellow, and red represent 68%, 95%, and 98% C.L. limits in the plane, respectively. As expected, the SM-like Higgs prefers a region with , the alignment limit [58].

The alignment limit can be understood in terms of the potential parameters in the 2HDM, couplings between the two doublets and masses , see Ref. [56] for more details. In terms of these parameters, the deviation from the alignment limit can be parametrised as .

Next, we introduce the additional nonlinear effects in Figure 5(a), where we show how the 95% C.L. regions are modified with nonlinear effects. The nonlinear effects, as in the case of the singlet, add up to the doublet mixing. Larger nonlinearities restrict further the parameter space of 2HDMs.

Finally, as the dependence is rather mild for 2, we fix the value of 6 to explore the interplay between and . In Figure 5(b), we plot the one-, two-, and three- contours in this plane. In the alignment limit, , we recover the pure nonlinear limit at 95% C.L. As we move from the alignment limit and allow a larger amount of doublet mixing, the limit on becomes milder, leading to a stronger bound on the compositeness scale . On the left side of the plot where , we find the limit on non-decoupling effects .25 at 95% C.L., which would correspond to the scale for the second doublet 500 GeV for .

2.3. CHM Pattern III: Extended Fermionic Sector

We finish this section on patterns, analysing a typical building block in composite Higgs scenarios: the presence of new fermionic degrees of freedom, not too far from the electroweak scale. These fermionic degrees of freedom are usually coming along as top partners, new vector-like fermions which mix with the SM top quark and modify its couplings with the Higgs.

We first focus on the mixing and present results for the -singlet top-partner scenario . We calculate the indirect bound on these states from measurements on the differential distributions of the the boosted Higgs in association with an energetic jet. The methodology used here was developed in Ref. [59] and further expanded in Refs. [60, 61]. We will then evaluate other sources of experimental constraints, including direct searches for the new state.

To describe the effect of the top-partner in the Higgs behaviour, we need to specify the mixing parameters between the top-partner and the SM top quark. For the singlet case, the mass matrix of the SM top quark and Dirac fermion top-partner can be written as where describes the mixing between and and is the Dirac mass of the top partner. After diagonalising this matrix by a biunitary transformation, we can calculate the mass eigenstates and and the mixing angle

Note that and are related in a simple way,

2.3.1. Top-Partner Indirect Searches

Contrary to the previous sections, total Higgs rates and their respective limits on vector and fermion couplings would not be sensitive to the presence of the top-partner. On the other hand, high- probes would access the top-partner mixing and mass scale indirectly [59, 60].

To analyse the effect of in the Higgs production in association with radiation, it is useful to define a quantity which depends on and based on semidifferential measurements,

In particular, we define a quantity based on these differential measurements, which exhibits good properties from the point of view of systematic and statistical fluctuations, see Ref. [59] for more details.

With this observable as an indirect probe for fermionic top-partners, we calculated and for 13 TeV by varying and . The results can be seen in Figure 6, where we plot the relative change in the Higgs+jet cross-section as a function of for different combinations of the mixing angle and . For this plot, is varied from 500 GeV to 1500 GeV, and one can see that has a weak dependence on . We can then use the experimental measurements in this channel to put a bound on the mixing angle. In particular, we used the ATLAS Higgs differential cross-section analysis [17], which reports an 8% deviation at 95% C.L. in the 250-350 GeV (diphoton) bin from the SM cross-section. We can then look back at our Figure 6 and note the variation of as a function of for =250 GeV.

In Figure 7, we show the variation of as a function of for a fixed  GeV. The green band corresponds to varying from 500 GeV to 1500 GeV. The 8% limit is noted with a hashed black line, which allows us to set a limit on . Using Equation (15), the corresponding limit in terms of the parameter is also obtained, see the lower panel of Figure 7.

Another source of constraints from a singlet VLQ top-partner would come from a SM effective field theory (SMEFT) global analysis of Run 2 LHC and LEP measurements. The singlet top-partner would produce a SMEFT pattern characterised by relations among some SMEFT operators, whereas all the other operators would be zero at tree level, see Ref. [62] for a dictionary between many extensions of the SM and their SMEFT matching, see also the discussion in Ref. [63] on the top-partners and the SMEFT framework. Using this dictionary and a global SMEFT analysis, in Ref. [64], the parameters of the singlet top-partner were bound to for  TeV.

2.3.2. Top-Partner Direct Searches

Boosted Higgs measurements provide one handle on top-partners. The LHC New Physics search programme includes very mature analyses looking for single and pair production of vector-like Quarks (VLQ). Pair production of VLQs is a dominantly strong interaction process, whereas single production relies on electroweak couplings.

A singlet up-type VLQ (T) is excluded up to 1.31 TeV by the ATLAS  TeV (36.1 fb) search for pair produced VLQs [65] after considering all decay modes. This limit is a combined limit from several searches for the pair production of VLQs.

From the CMS side, there are separate limits from each analysis of pair production of VLQs. The exclusion limit on singlet from dilepton final state of VLQs decaying to boson is 1280 GeV [66]. The CMS analysis for all leptonic final states excludes in the range 1140-1300 GeV corresponding to different branching fractions [67], and the analysis of the final state excludes singlet up to 1295 GeV [68]. The analysis targeting fully hadronic final states provides a bound up to 1.3 TeV for a specific combination of the branching fractions of VLQs [69].

The single VLQ production process is mediated by its coupling with SM particles; hence, it provides a bound on the mixing as a function of the mass of VLQs [1820]. Both the production cross-section and decay branching ratios are sensitive to mixing effects. The communication of the single-production results is a bit more cumbersome than for double production. In some analyses, this mixing is directly parametrised by the mixing angle () and sometimes by the coupling where [70] is just a rescaling of the mixing angle. ATLAS analysis cross-section limits are interpreted in terms of either , or both and . Note that both the angles and are related by Equation (16).

On the other hand, CMS exclusion limits on production cross-section times branching fraction are provided only as a function of . Additionally, the theoretical cross-section is provided for the fixed value of and . We use this additional information to calculate the approximate bound on as a function of .

We translate CMS bound in terms of mixing angle by parametrizing the production cross-section as , where is calculated from the theoretical cross-section with fixed . Then, using the bound on cross-section (including BR), we calculate the limit on which we further translate to using the relation . Note that we have not considered the effect of variation on the width of . We also consider other decay channels of the singlet vector-like top to and . Note that CMS analyses for [20] and [21] final states do assume 100% branching ratio, and the ATLAS collaboration has considered and top+MET final states [19] assuming BR from the singlet model.

In Figure 8, we show all these direct and indirect limits together in the plane of top-partner mass versus mixing angle. Double production of VLQs (dark-green vertical line) sets a limit on the mass 1.3 TeV. The SMEFT limit (light blue horizontal line) is very strong and independent of the top-partner mass as long as . The current sensitivity of the other indirect probe for top-partners, the +jet channel, is below the SMEFT fit, although it is worth noticing that our +jet analysis is not at the same level of sophistication as the SMEFT’s [64]. Note that for the ATLAS top+MET analysis (black line) the values covered by the closed contour are excluded.

3. Conclusions

In this work, we have studied the impact of the LHC Run 2 measurements on general composite Higgs scenarios.

The dataset we used includes Higgs signal strength measurements from CMS and ATLAS, differential properties of boosted Higgs production, several searches for single and double production of vector-like quarks and a specific result from a global fit to SMEFT properties obtained in Ref. [64].

We have described the most significant deviations expected from composite Higgs scenarios, going from the simplest effects, nonlinear couplings, to the most complex interplay between direct and indirect searches for additional fermionic degrees of freedom.

We have shown how in general composite Higgs scenarios one should expect simultaneous effects from more than one source. Typically, the composite nature of the Higgs is tied to the non-linear realisation of EWSB and the scale of compositeness . But generic composite Higgs scenarios exhibit a richer phenomenology, in particular new light scalars or fermions are also typical predictions of these scenarios.

These new degrees of freedom, if close enough to the electroweak scale, will also modify the Higgs couplings. For example, mixing of additional scalars participating in EWSB will add to nonlinearities to reduce even further the couplings of the Higgs to massive particles, leading then to even stronger limits on the scale of compositeness . Note that this scale is tied to the mass of new vector resonances like and and to the degree of tuning of Higgs potential. The Run 2 LHC data has pushed this scale into the TeV range, and we have shown that more complex scenarios would push this tuning even further. Hence, Run 3 and future LHC runs should have a good handle at natural composite Higgs scenarios.

Finally, we would like to mention that the same scale studied in this work could be linked to the the phenomenology of additional pseudo-Goldstone degrees of freedom which could appear in generic CHMs. The study of the correlations between the Higgs phenomenology presented here and the axion-like particles (ALPs) would be very interesting but goes beyond the scope of the paper.

Data Availability

This paper does not use data sources; it simply uses public cited results on limits from experimental searches. Therefore, no datasets have been generated during this study.

Disclosure

The manuscript is uploaded on arxiv with the number (arXiv:2102.13429 [hep-ph]) [71].

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work of V.S. is supported by the Science Technology and Facilities Council (STFC) under grant number ST/P000819/1. We want to thank Stephan Huber and Jack Setford for many fruitful discussions during the early stages of this work. We thank Andrea Banfi for the help regarding running the Higgs+jet process in his code implementation.