Constraining scalar resonances with topquark pair production at the LHC
Abstract
Constraints on models which predict resonant topquark pair production at the LHC are provided via a reinterpretation of the Standard Model (SM) particle level measurement of the topantitop invariant mass distribution, . We make use of stateoftheart Monte Carlo event simulation to perform a direct comparison with measurements of in the semileptonic channels, considering both the boosted and the resolved regime of the hadronic top decays. A simplified model to describe various scalar resonances decaying into topquarks is considered, including CPeven and CPodd, colorsinglet and coloroctet states, and the excluded regions in the respective parameter spaces are provided.
Contents:
I Introduction
With its mass being close to the electroweak scale the top quark is very special. It might intimately be connected to the underlying mechanism of electroweak symmetry breaking (EWSB). Consequently, studying topquark production and decays at colliders might provide a portal to New Physics (NP). The Large Hadron Collider (LHC), providing proton–proton collisions currently at 13 TeV centreofmass energy, can be seen as a topquark factory. It allows to search for anomalous topquark production and decay processes, considered as low energy modifications of the Standard Model (SM) parametrized by effective operators (1); (2); (3); (4); (5); (6), or, as the direct production of intermediate resonances, which have been hunted for a long time at different experiments (7); (8); (9).
Heavy scalar resonances that decay into a pair of top quarks are predicted by several NP scenarios, in particular the Two Higgs Doublet Model (THDM), supersymmetric theories and models of dynamical EWSB. In this paper, we provide a framework to reinterpret the SM differential cross section measurements as exclusion limits for signatures of NP resonances decaying into . The framework relies on the comparison between particlelevel data with stateoftheart event simulation and the interpretation of deviations in terms of NP models. It is based on four main ingredients

A Monte Carlo event generator which allows the precise and realistic description of particlelevel observables.
In order to theoretically describe topquark pair production at the LHC, we make use of stateoftheart event simulations provided by the Sherpa (10) eventgenerator framework. This implies the usage of techniques to match leading and nexttoleading order QCD matrix elements with parton showers and merging different partonmultiplicity final states.

A general parametrization of NP whose predictions for colliders can be computed efficiently. We adopt a Lagrangian which describes scalar resonances that can be CPeven or odd and color singlet or octet. We devise a reweighting method to describe the model prediction in the distribution for a wide range of the parameter space in a fast and efficient manner.

A statistical interpretation to decide what regions of parameter space of the model are ruled out at a given confidence level. We adopt here a simplified analysis.
A similar method to constrain NP with SM measurements in several other channels has recently been presented in Ref. (14). These approaches are complementary to modelspecific searches in the respective final states. They provide systematic methods for the theory community to derive more realistic exclusion limits for any particular model, not relying on the experimentspecific assumptions.
In the rest of the paper we explain these 4 points in detail. In Sec. II we describe the setup of our event simulation. In Sec. III we give details on the analyses used in the boosted and the resolved regime and validate our SM predictions by comparing them to experimental data. In Sec. IV we introduce our simplified model of beyond the SM scalar resonances and describe the implementation in our simulation framework, based on an eventbyevent reweighting. In Sec. V we present a statistical analysis to assess the region in parameter space accessible by the LHC experiments and provide interpretations in terms of some specific models. We finally conclude in Sec. VI.
Ii Simulation framework
When searching for imprints of resonant contributions in topquark pair production at the LHC, a detailed understanding of the SM production process is vital. In particular, as there are nontrivial interference effects between NP signals and SM amplitudes that determine the shape of the resulting toppair invariantmass distribution. In order to obtain realistic and reliable predictions for the toppair production process, we make use of stateoftheart particlelevel simulations, based on higherorder matrix elements matched to partonshower simulations and hadronization.
Our analysis focuses on observables in the semileptonic decay channel of topquark pair production, i.e.
(1) 
where denotes muons or electrons, the corresponding neutrinos, are bottom quarks and light quarks or gluons. These decay products and the associated radiation might be reconstructed as wellseparated objects, i.e. lightflavour jets, jets and a lepton, or, in the boosted regime, as a largearea jet, containing the hadronic decay products, additional jets and a lepton. In either case, to realistically simulate the associated QCD activity, higherorder QCD corrections need to be considered.
To describe the SM toppair production process we use the Sherpa eventgeneration framework (15); (10). We employ the techniques to match LO and NLO QCD matrix elements to Sherpa’s dipole shower (16) and to merge processes of variable partonic multiplicity (17); (18). Leadingorder and realemission correction matrix elements are obtained from Comix (19). Virtual oneloop amplitudes, contributing at NLO QCD, are obtained from the Recola generator (20); (21) that employs the Collier library (22). Topquark decays are modelled at leadingorder accuracy through Sherpa’s decay handler, that implements BreitWigner smearing for the intermediate resonances and preserves spin correlations between production and decay (23). We treat bottomquarks as massive in the topquark decays and the finalstate partonshower evolution (24).
To validate the SM predictions we also consider leadingorder simulations in the MadGraph_aMCNLO framework (25). The hardprocess’ partonic configurations get showered and hadronized through Pythia8 (26). The spincorrelated decays of top quarks are implemented through the MadSpin package (27). Samples of different partonic multiplicity are merged according to the MLM prescription described in (28).
For the topquark and boson, the following mass values are used
(2) 
and the corresponding widths are calculated at leading order, assuming for the remaining electroweak input parameters and . In the following section we present a comparison of our simulated predictions against ATLAS measurements and discuss their systematics. Alongside, we give details on the QCD input parameters and calculational choices used there.
Iii Analysis framework
In what follows we describe the event selections used to identify the topquark pairproduction process, used later on to study the imprint of resonant NP contributions. Thereby, we closely follow the strategies used by the LHC experiments. Our simulated events from Sherpa and MadGraph_aMCNLO are produced in the HepMC output format (29) and passed to Rivet (11) where we implement our particlelevel selections.
We consider two analyses, based on measurements performed using the ATLAS detector of the differential production cross sections in protonproton collisions at TeV with an integrated luminosity of (12); (13). Both analyses select events in the leptons+jets decay channel. The two measurements indicated in the following as Resolved and Boosted are optimized for different regions of phase space. The Boosted analysis, cf. Ref. (12), is designed to enhance the selection and reconstruction efficiency of highlyboosted top quarks with transverse momentum 300 GeV, that might originate from the decay of a heavy resonance with mass . In such events the decay products of the hadronic top overlap, due to the high Lorentz boost. In turn, they cannot be reconstructed as three distinct jets. The Resolved analysis, based on Ref. (13), measures the differential cross section as a function of the full kinematic spectrum of the system and is useful to identify and reconstruct rather light resonances.
The selection requirements are applied on leptons and jets at particle level, i.e. after hadronization. In our simulated data we discard any detector resolution, i.e. smearing effects. All the leptons used in the analyses, i.e. and must not originate from hadrons, neither directly nor through a lepton decay. In this way the leptons are guaranteed to originate from boson decays without a specific matching requirement. The fourmomenta of the charged leptons are modified by adding the fourmomenta of all photons found in a cone of = 0.1 around the leptons’ direction, thus representing dressed leptons. The missing transverse energy of the events () is defined from the fourvector sum of the neutrinos not resulting from hadron decays.
Jets are clustered using the anti algorithm (30) with a radius of for smallR jets and for the largeR jets, using all stable particles, excluding the selected dressed leptons, as input. All smallR jets considered during the selections are required to have 25 GeV and 2.5, while for largeR jets we demand 300 GeV and 2. The smallR jets are considered tagged if a hadron with 5 GeV is associated to the jet through a ghostmatching procedure (31); (32). To remove most of the contribution coming from the interaction of the proton remnants, i.e. the underlying event, and to reduce the dependence on the generator, largeR jets are groomed following a trimming procedure with parameters and = 0.05, for details of the procedure see Ref. (33).
Both the Resolved and the Boosted selections require a single lepton with 25 GeV and 2.5. In the Resolved analysis, apart from the leptons, the events are required to have at least four smallR jets and at least two of them have to be tagged. In the Boosted analysis the events are required to have 20 GeV and , with , the transverse mass of the leptonically decaying boson, where denotes the azimuthal angle between the lepton and the vector. The presence of at least one smallR jet with smallR jet) is required. In case more than one jet fulfills this requirement the jet with higher is considered as the jet originating from the leptonic top decay, dubbed lepjet candidate. Furthermore, it is required the presence of a trimmed largeR jet with mass GeV and 40 GeV, where is the distance (34); (35) between the two subjets in the last step of the jet reclustering, i.e. . If more than one largeR jet fulfills these requirements the one with highest transverse momentum is considered as the hadjet candidate. The hadjet candidate must furthermore satisfy certain kinematic requirements: and . The final requirement in the Boosted selection is that at least one tagged jet in is found or that the lepjet candidate is tagged. The Resolved and Boosted event selections are summarized in Tab. 1.
nolistsep
event selections  

Exactly one lepton ( or ) with GeV and  
Resolved analysis  Boosted analysis 
GeV and 60 GeV  
4 smallR jet:

1 largeR jet:

1 smallR jet:


2 tagged jets 
1 tagged jet:

For the selected events the system is reconstructed based on the event topology:

Resolved analysis: The leptonic top is reconstructed using the tagged jet nearest in to the lepton and the missingmomentum four vector, the hadronic top is reconstructed using the other tagged jet and the two light jets with invariant mass closest to the mass.

Boosted analysis: The leptonic top is reconstructed using the lepjet candidate, the lepton and the missingmomentum four vector, the hadjet candidate is directly considered as the hadronic top.
In order to validate our simulations of SM topquark pairproduction we compare our predictions against ATLAS data for the Boosted and Resolved selection, supplemented by studies of systematic variations. To begin with, we check the impact of the grooming procedure on the reconstructed hadronictop candidate mass, i.e. the mass of the hadjet candidate in the Boosted event selection. We consider event samples from Sherpa and MadGraph_aMCNLO, based on the leadingorder matrix element for topquark pair production, labelled as . In these calculations, i.e. without mergingin highermultiplicity matrix elements, we set the renormalization () and factorization scale () to
(3) 
with () the transverse momentum of the decaying (anti) top quark.
In Fig. 1 we present the resulting invariantmass distributions obtained from Sherpa and MadGraph_aMCNLO before and after applying the grooming procedure. Comparing the untrimmed distributions (left panel) both samples exhibit a clear peak at the nominal topquark mass. However, due to partonshower radiation and nonperturbative corrections from hadronization and underlying event the peak is rather broad and sizeable differences are observed when comparing the predictions from Sherpa and MadGraph_aMCNLO+Pythia8. Note that the uncertainty bands shown represent the statistical uncertainty of the samples only. When applying the trimming procedure to the hadjet candidates the mass distributions agree to a much better degree, both in the tails of the distribution and the peak region. Therefore, trimming of the largeR jets significantly reduces the dependence on the generator and the details of its partonshower formalism and the modelling of nonperturbative effects.
In Figs. 2 and 3 we compare predictions from Sherpa based on LO and NLO matrix elements against data measured by the ATLAS experiment for the Boosted (left panels) and the Resolved (right panels) event selections. For the MEPS@LO sample we merge LO QCD matrix elements for jet production dressed with the Sherpa dipole parton shower (17). The mergingscale parameter is set to . The MEPS@NLO sample combines QCD matrix elements at NLO for jet and jets at LO according to the methods described in (36); (18), again using a merging scale of . Both methods share the eventwise reconstruction of an underlying core process through consecutive clusterings of the external legs. For this reconstructed core process the renormalization and factorization scales are set to , with
(4) 
For the reconstructed clusterings the strong coupling is evaluated at the respective splitting scale. The scale is furthermore used as the resummation, i.e. partonshower starting scale, denoted . To assess the scale uncertainty of the predictions we perform variations by common factors of and for the core scale and the local splitting scales, using the eventreweighting technique described in (37). In the figures the resulting uncertainty estimate is represented by the red band, while the blue band indicates the statistical uncertainty.
For the boostedtop selection we show the transversemomentum distribution of the hadronictop candidate in the left panels of Figs. 2 and 3, respectively. Notably, both samples, i.e. the MEPS@LO and the MEPS@NLO prediction, describe the ATLAS measurement (12) very well, both in terms of the production rate and in particular concerning the shape of the distribution. For the MEPS@LO result the scale uncertainty is quite significant, reaching up to 50%. However, the dominant effect is a mere rescaling of the total production rate, the shape of the distribution stays almost unaltered. This is also observed for the MEPS@NLO sample, however, the scale uncertainty reduces to .
For the resolveddecay selection we compare the Sherpa MEPS@(N)LO predictions for the reconstructed invariant mass of the system against data from the ATLAS experiment (13), see right panels of Figs. 2 and 3. Note that the data and the theoretical predictions are normalized to their respective fiducial cross section. The MEPS@LO and MEPS@NLO results agree very well with the data. For this normalized distribution the scale uncertainties largely cancel. For the MEPS@LO sample this results in an uncertainty estimate of . For the MEPS@NLO sample the shape modifications induced by the scale variations amount to .
For both observables considered, the MEPS@(N)LO predictions from Sherpa yield a very satisfactory description of the data. No significant alteration of the distributions shape is observed upon inclusion of the QCD oneloop corrections in the MEPS@NLO sample. However, in particular the uncertainty on the production rate reduces significantly. For the normalized toppair invariant mass distribution we consider the more realistic estimate from the MEPS@NLO calculation. By normalizing the distribution to the cross section in a certain mass window, this uncertainty might in fact be reduced further, cf. Ref. (38), where, ultimately, an uncertainty estimate of was quoted for the corresponding NNLO QCD prediction.
In what follows we want to study the imprint of New Physics resonant contributions on the toppair invariant mass distribution. To this end we currently rely on a leadingorder description of the signal, interfering with the corresponding SM amplitudes. However, from the considerations above we can conclude that the MEPS@LO calculation of the SM production process captures the dominant QCD corrections, which are of realradiation type. To illustrate this further, we present in Fig. 4 a comparison of MEPS@LO samples using different partonmultiplicity matrix elements for the mass and the transverse momentum of the system in the Boosted selection. These results get compared to the corresponding MEPS@NLO prediction described above.
For the toppair invariant mass all the predictions with at least one extra hard jet agree within their statistical errors. In particular, even the MEPS@LO sample, based on merging the LO matrix elements for jet only, well reproduces the MEPS@NLO result and greatly improves the 0jet sample. As might be expected, for the transverse momentum of the system, the inclusion of highermultiplicity matrix elements improves the agreement with the MEPS@NLO result. The MEPS@LO calculation based on jet predicts a somewhat softer spectrum, i.e. is lacking configuration corresponding to multiple hard emissions. However, the bulk of the events in the Boosted selection is reasonably modeled by this simple LO merging setup and describes the data presented above very well. We will therefore rely on this setup when invoking New Physics contributions.
In the following we also introduce a simple Parton Analysis, used to quantify the effect of the NP without any smearing due to the reconstruction of the top quarks. In the Parton Analysis no cuts are applied to the events and the two top quarks are identified, before any decay, using truthlevel information from the generator.
Iv Simplified model
Several models of NP predict resonances decaying to topquarks. Scalar resonances in particular have large branching ratios in this decay channel due to the fact that their couplings with fermions are often proportional to the fermion masses. In this case, the resonance is at the LHC dominantly produced via gluon fusion through loops of colored particles. These colored particles can be either light compared to the resonance (like the top quark itself), in which case the structure of the loop is resolved as illustrated in Fig. 5(a), or they can be heavy, in which case a pointlike interaction sketched in Fig. 5(b) can describe the interactions.
It has been shown in (39) that the most general scalar extension of the SM which couples to fermions and maintains naturally small flavour changing neutral currents is provided by scalars with the same quantum numbers of the Higgs doublet or that transform as a color octet under the SU(3)SU(2)U(1) SM gauge group. Color neutral and octet scalars arise also naturally in several models of dynamical EWSB, such as in the seminal FarhiSusskind model (40) and models where the top is partially composite (41). Although the specific origin of the scalartop couplings is important, determining the relation to other couplings and their magnitudes, we here adopt a more phenomenological simplified approach relevant for topquark pair production, in which the lefthanded top is stripped off from its doublet and couples directly to the scalars.
In our simplified model we assume the only light state running in the loop to be the topquark. This is a good approximation if two conditions are fulfilled: (i)  the bottomquark contribution is suppressed; and (ii)  the extra states contributing significantly to the gluon–scalar couplings are heavy (at least as much as the scalar resonance itself). This is a good approximation in many models beyond the SM. In the THDM (42) for example, there is no new particle living at higher scale apart from the new scalar sector. Moreover, the loop of bottomquarks is usually suppressed in the cases relevant for production. Specializations of the THDM such as the Minimal Supersymmetric Standard Model (MSSM) where the superpartners are heavy enough to be integrated out can also be described in this framework. Composite models typically predict relatively degenerate spectra of first excitations, thus they can be usually described by the effective pointlike interaction. Similarly, for the color octet in the model of Manohar and Wise (39) the scalars are produced purely by top and bottom loops. In some other models intermediate states much lighter than the first scalar excitations are present, e.g. top partners and stops may be light in some models of partial compositeness and SUSY – in these cases our approximation is not applicable.
Under this assumption we can describe the scalar sector interactions relevant for production via the following Lagrangian:
(5)  
It contains a CPodd isosinglet scalar , a CPeven isosinglet scalar , a CPodd color octet scalar and a CPeven octet scalar which we collectively call . is the gluon fieldstrength tensor, , are the SU(3) generators and is the fully symmetric SU(3) tensor.
The topquark loops generate form factors that describe the gluonscalar interaction. The loop triangles contribute to the trilinear vertices in the form
(6)  
(7) 
with
(8)  
(9)  
(10) 
Similar expressions for the color octet topquark loop generated form factor can be found e.g. in (43).
As a matter of fact, resonant top pair production is accompanied by other signatures. In particular, diphoton, dijet, , and signatures are generated via diagrams induced by a topquark loop, and in general by highscale physics. Treelevel , decay channels are typically present for a scalar state, while decays into lighter fermions are typically suppressed. Color octets decays into and might give striking signatures. The detailed analysis of these channels is not in the scope of this work, however, we provide some qualitative discussion about the regions in parameter space where they can be competitive in sensitivity to search.
Loop (or anomaly) induced decays are typically suppressed and might be competitive to searches only for small Yukawa couplings . They are often the only possible decay channels for pseudoscalars besides that into . As an example, consider some partial widths of a colorsinglet pseudoscalar
(11)  
(12)  
(13) 
Here we parametrize the photon interaction with by the following gauge invariant operators
(14) 
with . These operators also give rise to decays into weak bosons, but not competitive in sensitivity to diphoton searches (unless there is some cancellation in ). From the above expressions it can be noticed that the partial width is much larger than , however, the corresponding search is not as competitive to the diphoton channel due to the clean signature of the latter.
On the other hand, scalar resonances tend to decay into weak bosons at tree level, with large contributions to their decay width and good sensitivity in the corresponding channels.
The color octets have more unexplored signatures, like e.g. , studied for example in Ref. (44); (45).
iv.1 Model description and simulation
Our goal is to achieve accurate predictions for a wide parameter range of our generic model in an efficient and fast way. For this purpose, the Lagrangian given in eq. (5) has been implemented into the FeynRules (46) package to produce a corresponding UFO model file (47). The required helicity amplitudes have been extracted to C++ codes via the Madgraph (48) program and incorporated in the Rivet analyses in order to perform a reweighting method and reproduce the signal lineshape. To this end, each event of the Sherpa SM event sample is given a weight, , proportional to the ratio of the amplitudes,
(15) 
where is the SM amplitude squared summed and averaged over color and spin. In the numerator the amplitude corresponding to the resonant diagrams depicted in Fig. 5 is added on top of the SM diagrams. The further decay of top quarks is included neglecting nonresonant diagrams. Therefore, the full process in eq. (1) – including possible extra hard radiation – is considered with full spin correlation of the topquark decays.
We note that our signal includes not only the purely resonant contribution. The complete squared amplitude can be split into three contributions:
(16) 
The last term defines the SM background (), the pure signal () and the interference between signal and SM ().
We use as the test observable the distribution of the signal hypothesis normalized binbybin to the SM QCD prediction,
(17) 
The signal hypothesis differential cross section is defined as the total differential cross section subtracted by the SM prediction. Such normalized distribution is less affected by systematic errors, i.e. theoretical uncertainties (49).
In order to assess the importance of the interference we study both the full signal including interference and the pure signal hypothesis neglecting interference . To simplify the notation in the remaining of the text we use the following definitions:
(18) 
Interference between signal (Fig. 5) and QCD diagrams are known to be important in this process. In fact, they can completely change the lineshape of the resonance from a pure BreitWigner peak to a peakdip structure, or even dippeak, pure dip or an enhanced peak (50); (51); (52); (53); (54); (55). QCD corrections to this effect have recently been computed (56); (57); (58) and shown to be important. A pilot experimental analysis investigating such interference effects has been presented recently (59).
The form factors in eq. (7) have been implemented in the helicity amplitudes used in the reweighting step. However, the corresponding box diagram contributing to the fourgluon–scalar coupling was kept as an effective vertex without momentum dependence. For the color octet the form factor is approximated by a fixed momentum flowing through the loop that is equal to the mass of the resonance. The interference between topquark loops and pointlike interactions is also manifest in the calculation.
Higherorder QCD corrections are partially taken into account through the radiation of extra gluons in the MEPS@LO simulation. The contribution from realemission matrix elements also get reweighted with the NP theory hypotheses.
We note however that the method neglects the signals’ colorsinglet color flow contribution when attaching parton showers, which affects the subsequent radiation pattern only. We nevertheless found that these effects are small in the description of the toppair mass distribution. In Fig. 6 we show the distribution of variable defined in eq. (17) for a colorsinglet pseudoscalar of mass TeV in the pure signal hypothesis, comparing the Sherpa reweighted events with a dedicated simulation of the full process with MadGraph_aMCNLO+Pythia8. In the latter, the colorflow contribution corresponding to the signal diagrams are considered as seeds for the subsequent parton shower. The error near the resonance peak is about 10% and the reweighted prediction underestimates the yields. We removed the topquark loop form factor considering only the effective scalargluon coupling for this comparison. The distributions were derived according to the Parton Analysis framework described in sec. (III). In the more realistic boosted analysis we expect the reweighting method to predict a more smeared distribution due to the extra connected color lines that favor extra hard radiation connecting the top quarks with initial gluons. We will neglect these effects and employ the reweighting method in what follows to make predictions for a large region of parameter space of the model, while avoiding massive time and machine consuming event generation and “fake” MC statistical error. Our results are expected to give conservative limits since for coloursinglet resonances the signal color flow induces less smearing of the resonance peak.
V Results
Resonant topquark pair production at the LHC has been analyzed for several of the models mentioned above already. Color neutral resonances decaying into have been studied in several works for a large number of models (50); (51); (53); (54); (55), even including interference effects at NLO in QCD (56); (57); (58). The case of a coloroctet signal has been considered in (52); (60); (61); (43), also considering other production channels, e.g. via initial states, or even double scalar production (62); (63); (64). Our approach differs from previous studies because we adopt the strategy of directly comparing to data which has been shown to agree well with the SM prediction, and therefore, can be used to put direct limits on the model parameters, in the same spirit as (14). Indeed, the recent ATLAS measurement of the topquark pair differential cross section at shows good agreement with various SM Monte Carlo generators (65). However, there are no measurements of the invariant mass in the boosted regime at this energy yet. Moreover, the uncertainties are still quite large, since only the 2015 data, corresponding to , were used, but we expect that an update of the analysis will be available in the near future, with improved systematics and statistical uncertainty (comparable to the ones presented in this paper) allowing to derive real exclusion limits. We assume in what follows that data will be well described by the SM expectation, and take the SM prediction from Sherpa as mock data.
The method proposed allows theorists to derive realistic exclusion limits on a variety of NP scenarios without a dedicated and expensive experimental analysis. It opens a new path to search for NP, with the experiments providing precision measurements of SM processes. With respect to dedicated experimental searches, it can serve as check and as an alternative (lessexpensive) approach to look for more general parametrizations of deviations caused by New Physics. For instance, in the ATLAS and CMS collaborations’ analyses (8); (9); (66); (67); (59), only a leptophobic Z’ bosons (present for instance in topcolor scenarios), a KaluzaKlein excitation of the gluon and heavy states in THDM were searched for. Moreover, interference effects were considered only in Ref. (59). With our technique we are able to provide limits for a whole wealth of models.
In order to assess the possibility to observe the signals described above we perform a simple analysis using the bins of the distribution. We consider the mass window and compute
(19) 
with the number of bins taken into account, according to the assumed resolution of the measurement. is the distribution integrated over bin and is the hypothesis (either or ). is the variance on each bin of the distribution.
The variance is derived according to the rules of propagation of uncertainties and is estimated by
(20) 
We kept the indexes implicit in the expression. The first term accounts for statistical error, the second for systematic uncertainties of experimental sources, and the third for theoretical uncertainties. We assume a flat distribution for theory and systematic uncertainty, and that statistical uncertainties are dominated by the background, with a small ratio signal over background. We take for both and , assuming other errors are strongly correlated and will be canceled when taking the ratio distribution. The experimental uncertainty is more important and we consider three benchmark estimates for :

In Ref. (59) the total systematics on the background were estimated as 10% and 11%. As a pessimistic case we consider .

As an optimistic scenario we vary it to lower values considering a future improved understanding of the uncertainties and the reduction in uncertainty associated to normalization. Since we are using a normalized distribution many of the uncertainties estimated in the previous benchmark are strongly correlated and will be canceled out. For this we use .

As the most optimistic case we assume experimental uncertainties can be drastically reduced to the level of theoretical, which according to Ref. (49) results in .
We consider for a bad resolution case, assuming the experiment can resolve only the full window of 400 GeV in , and assuming a mass resolution in of GeV.
We consider as a criterion for exclusion, which corresponds roughly to an exclusion at 95% of confidence level.
This simple analysis is intended to be a first approximation to a full statistical data analysis that will be carried out eventually. In particular we assume the same uncertainty for every bin without correlation between them, and we assume only two cases of resolution independent of the bin. In the following we discuss some benchmark scenarios and the respective results.
v.1 Pseudoscalar color octet
The first scenario we consider is when the resonance represents a pseudoscalar color octet () with total width dominated by the decays to pairs of tops and gluons
(21) 
In Fig. 7 we show the resulting distribution assuming a color octet resonance with mass GeV and the parameters , (left) and (right) at parton level, i.e. using the Parton Analysis described in sec. (III). We show both the full lineshape, which comprises signal and interference with QCD background (S+I), and the pure signal (S) for comparison. The importance of taking into account interference effects can clearly be noticed.
Similarly, in Fig. 8, we present the effect of a resonance with mass and couplings , (left) and (right), reconstructed using the Boosted Analysis. The excess reaches more than 10%, which indicates that even a pessimistic estimate of the uncertainties is sufficient to exclude the existence of this state for values of of order 1. We thus use the most pessimistic value for the systematic error, .
In Fig. 9 the corresponding exclusion limits are shown, assuming a fixed value of . The bands correspond to a systematic uncertainty on the measurement running from to . The limits are evaluated considering the interference effect (dashed lines) or neglecting it (continuous lines). The interference has a significant effect in the low mass region (). The excluded region corresponds to larger values of . We show the exclusion for integrated luminosities of (blue line) and (black). In the leftpanel we use bins of width in the invariantmass distribution to compute while on the rightpanel we use only a single 400 GeV bin centered around the resonance mass, . The comparison between the left and right panel shows the importance of a good resolution and for a lineshape analysis.
We expect striking signatures in other channels, but little has been studied. For instance, in the analysis of +jets in Ref. (45) a color octet has not been considered.
v.2 Pseudoscalar singlet
For the benchmark scenario of a pseudoscalar color singlet we again assume the resonance’ width is dominated by the top and gluon decays, as in eq. (21).
We show in Fig. 10 the distribution of the normalized distribution assuming . In the lefthand (righthand) panel we consider (). The lineshapes of this scenario are highly nontrivial, they strongly depend on the mass and couplings, and can feature pure dips, pure peaks and intermediate peakdip or dippeak structures. A sample of different lineshapes is shown in app. (A).
In Fig. 11 we show the exclusion limits in the parameter space plane for . The band represents the different assumptions for the systematic uncertainty, 5% and 10%. The effect of interference is important for low masses , where also systematics dominate and have a huge impact on the exclusion power. The use of the full lineshape in the statistical analysis improves the exclusion power mostly for low masses where more distinct lineshapes are present. For masses above , higher luminosities than are needed.
In Fig. 12 we show the corresponding exclusion limits in the plane for a fixed mass . The effect of interference is important for large top couplings, , which is directly related to the size of the width. The use of full lineshape gives a mild improvement in the exclusion power.
For very low masses the Resolved analysis can be slightly more powerful than the Boosted. In Fig. 13 on the left we show an example of a lineshape and on the right the exclusion limit provided by the Resolved analysis. Compared to fig. (11) it can be noticed that the low mass region can be better covered by the Resolved selection. We note as well that the case of negative is less excluded due to the fact that larger cancellations between topquark loop and effective vertex happens for these masses.
Diphoton and dijet searches might be relevant in extreme regions of parameter space, i.e. for very small , and large masses, due to the dependence of the and partial widths on as opposed to the linear dependence of the decay width. In fig. (14) we show the 95%CL excluded region derived from the limits provided by the ATLAS collaboration in the dijet search (68). We used the case and assumed an acceptance of 50%. In the same figure we show the 95%CL excluded region in the diphoton channel using the exclusion limits by the ATLAS analysis in Ref. (69). We used the case and the spin0 selection. To derive cross sections we used the NLO result for Higgs production cross section (70) and rescale by the LO decay width,
(22) 
is given in eq. (13) and the form factors in eqs. (8–10). The shaded area in the figure represents the region where BR is larger than the excluded line in the respective references, and BR is the corresponding branching ratios. We can notice that these channels get competitive in sensitivity to analysis at low and large mass, but only if is particularly large. In particular, even for , for the dijet search seems to be more sensitive to New Physics.
Interpretation for Composite Higgs models with Top Partial Compositeness
As an ultraviolet realization of the pseudoscalar scenario we consider the composite models M3, M8 and M9 of Ref. (41). These models are constituted by two additional confining fermions, and , which form several composite states among which a top partner that can generate a mass to the top quarks through the partialcompositeness mechanism. In addition, they present two isosinglet pseudoscalar mass eigenstates and . In general, the observation of such pseudoscalar state decaying into top quarks can shed light on the mechanism of fermion mass generation (71). These models present extra parameters which determine the couplings, given by a pair of integers and the relation between the mixing angle and the ratio of scales and U(1) charges, . We do not enter a discussion of the details of these the models and their parameters here but invite the reader to consult Ref. (41). We choose and the values of which provide the largest couplings to the tops, and . We neglect contributions to the resonance width from the decays into , and , which are subdominant. The relevant couplings are summarized in Tab. 2.
model  

M3  0.934/1.09/2.65  5.44 
M8  0.926 /1.54/2.16  1.54 
M9  0.293/0.195/1.37  8.6 
In Fig. 15 we show the value of and for each model together with the exclusion region (above the black curve) for a fixed mass . We consider an integrated luminosity of for the exclusion limit and a systematic error . The different line colors in the figure refer to the different models: red is M8, yellow M9 and brown M3. The styles of the lines represent the fermionic charges: (solid line), (4,2) (dashed) and (4,2) (dotdashed). Each line scans the values of from (most external and largest couplings) to (most internal and smallest couplings), the dots represent the values , with an integer between 1 and 8 included. Also shown for reference in the upper region the couplings of the state of the FahriSusskind onefamily model (40).
From the figure we can get the minimal value of the compositeness scale for which state would still not have been observed for different scenarios. For instance, for model M8 (red lines), depending on the values of . The model M3 is more constrained, and for (4,2) and for (2,0) or (4,2). Model M9 has low values of but values can be excluded for the case (4,2), while the other scenarios are hard to access in the search.
Other decay channels have been analyzed in Ref. (41).
v.3 Broad scalar color singlet
In this benchmark scenario we assume a CPeven colorsinglet scalar that can, apart from top quarks and gluons, also decay into other particles and is thus much broader than the previous scenarios. We choose a total width of 20% of the resonance mass . The rationale for choosing a larger width is the fact that the scalar tends to decay also to weak bosons. Indeed, we expect a large sensitivity in this decay channel which might be competitive w.r.t. top pair production.
In this scenario the signal is very weak and thus hard to be observed unless the systematic uncertainty is improved to values below 5% or higher values of are considered. In Fig. 16 on the left we show the lineshape for , . It can be noticed that the yields are always below 5%. On the right panel we show the contours in the parameter space plane for . Varying the assumed systematic uncertainties between determines the band of the exclusion limit. The integrated luminosities are (blue line) and (black). Limits are given considering interference (dashed lines) and neglecting it (solid lines). A large interference effect can be noticed, which is in fact larger than the pure signal.
Vi Conclusion
In this work we have provided a framework to reinterpret the SM differential cross section measurements in terms of exclusion limits for signatures of NP scalar resonances decaying into . The method relies on the detailed simulation of the SM prediction at particle level with the Sherpa Monte Carlo, the subsequent analysis in the Rivet framework, which can be directly compared with the measured distributions provided by the experimental collaborations, a modeling of the NP scenarios efficient enough to allow a scan over a large range in parameter space, and finally a statistical analysis to determine the excluded regions.
In the simulation of toppair production we take into account higherorder QCD corrections through matching LO or NLO matrix elements to parton showers and merging partonic processes of varying multiplicity. To validate our simulation we compare to data from the ATLAS collaboration, finding very good agreement. As New Physics contributions we consider CPeven and CPodd scalar resonances, being either colorsinglets or octets. To model the signal we devise an efficient and fast reweighting method allowing to scan large regions of parameter space without the need of full resimulation and reanalysis for each parameter point. For our simplified model we have derived exclusion limits based on a simple analysis, that can subsequently be used to set limits on other specific models, and we consider a model of partial compositness as an example. We showed the importance of properly accounting for interference between the New Physics signal and the SM background in setting the exclusion limit, as well as of using a full lineshape analysis which is not necessarily a simple BreitWigner shape due to the interference effects.
By confronting SM precision measurements with hypotheses for New Physics models stringent exclusion limits on the parameters of the latter can be obtained, providing complementary sensitivity to direct searches. The methodology laid out here can be readily applied to other observables than the toppair invariant mass considered here. It relies on a solid understanding of the respective SM expectation and the uncertainties related to the theoretical predictions and the experimental data.
Acknowledgments
This work was supported by the European Union through the FP7 training network MCnetITN (PITNGA2012315877) and the Horizon2020 Marie SkłodowskaCurie network MCnetITN3 (722104). Federica Fabbri especially wants to thank MCnetITN for the opportunity to hold a shortterm studentship at the II. Institute for Physics at Göttingen University.
Appendix A Lineshapes samples
Footnotes
 preprint: MCnet1717
References
 V. Barger, W.Y. Keung, and B. Yencho, Azimuthal Correlations in Top Pair Decays and The Effects of New Heavy Scalars, Phys. Rev. D85 (2012) 034016, [arXiv:1112.5173].
 S. Choi and H. S. Lee, Azimuthal decorrelation in production at hadron colliders, Phys. Rev. D87 (2013), no. 3 034012, [arXiv:1207.1484].
 D. Buarque Franzosi and C. Zhang, Probing the topquark chromomagnetic dipole moment at nexttoleading order in QCD, Phys. Rev. D91 (2015), no. 11 114010, [arXiv:1503.08841].
 C. Zhang, Single Top Production at NexttoLeading Order in the Standard Model Effective Field Theory, Phys. Rev. Lett. 116 (2016), no. 16 162002, [arXiv:1601.06163].
 C. Englert, L. Moore, K. Nordstr\unichar246m, and M. Russell, Giving top quark effective operators a boost, Phys. Lett. B763 (2016) 9–15, [arXiv:1607.04304].
 V. Cirigliano, W. Dekens, J. de Vries, and E. Mereghetti, Constraining the topHiggs sector of the Standard Model Effective Field Theory, Phys. Rev. D94 (2016), no. 3 034031, [arXiv:1605.04311].
 CDF Collaboration, T. Aaltonen et al., Search for New ColorOctet Vector Particle Decaying to t antit in p antip Collisions at s**(1/2) = 1.96TeV, Phys. Lett. B691 (2010) 183–190, [arXiv:0911.3112].
 CMS Collaboration, V. Khachatryan et al., Search for resonant production in protonproton collisions at 8\unichar8201\unichar8201TeV, Phys. Rev. D93 (2016), no. 1 012001, [arXiv:1506.03062].
 CMS Collaboration, S. Chatrchyan et al., Searches for new physics using the invariant mass distribution in pp collisions at =8\unichar8201\unichar8201TeV, Phys. Rev. Lett. 111 (2013), no. 21 211804, [arXiv:1309.2030]. [Erratum: Phys. Rev. Lett.112,no.11,119903(2014)].
 T. Gleisberg, S. Höche, F. Krauss, M. Schönherr, S. Schumann, F. Siegert, and J. Winter, Event generation with SHERPA 1.1, JHEP 02 (2009) 007, [arXiv:0811.4622].
 A. Buckley, J. Butterworth, L. Lönnblad, D. Grellscheid, H. Hoeth, J. Monk, H. Schulz, and F. Siegert, Rivet user manual, Comput. Phys. Commun. 184 (2013) 2803–2819, [arXiv:1003.0694].
 ATLAS Collaboration, G. Aad et al., Measurement of the differential crosssection of highly boosted top quarks as a function of their transverse momentum in = 8 TeV protonproton collisions using the ATLAS detector, Phys. Rev. D93 (2016), no. 3 032009, [arXiv:1510.03818].
 ATLAS Collaboration, G. Aad et al., Measurements of topquark pair differential crosssections in the lepton+jets channel in collisions at TeV using the ATLAS detector, Eur. Phys. J. C76 (2016), no. 10 538, [arXiv:1511.04716].
 J. M. Butterworth, D. Grellscheid, M. Kr\unichar228mer, B. Sarrazin, and D. Yallup, Constraining new physics with collider measurements of Standard Model signatures, JHEP 03 (2017) 078, [arXiv:1606.05296].
 T. Gleisberg, S. Höche, F. Krauss, A. Schälicke, S. Schumann, and J.C. Winter, SHERPA 1. alpha: A Proof of concept version, JHEP 02 (2004) 056, [hepph/0311263].
 S. Schumann and F. Krauss, A Parton shower algorithm based on CataniSeymour dipole factorisation, JHEP 03 (2008) 038, [arXiv:0709.1027].
 S. Höche, F. Krauss, S. Schumann, and F. Siegert, QCD matrix elements and truncated showers, JHEP 05 (2009) 053, [arXiv:0903.1219].
 S. Höche, F. Krauss, M. Schönherr, and F. Siegert, QCD matrix elements + parton showers: The NLO case, JHEP 04 (2013) 027, [arXiv:1207.5030].
 T. Gleisberg and S. Höche, Comix, a new matrix element generator, JHEP 12 (2008) 039, [arXiv:0808.3674].
 S. Actis, A. Denner, L. Hofer, J.N. Lang, A. Scharf, and S. Uccirati, RECOLA: REcursive Computation of OneLoop Amplitudes, Comput. Phys. Commun. 214 (2017) 140–173, [arXiv:1605.01090].
 B. Biedermann, S. Br\unichar228uer, A. Denner, M. Pellen, S. Schumann, and J. M. Thompson, Automation of NLO QCD and EW corrections with Sherpa and Recola, Eur. Phys. J. C77 (2017) 492, [arXiv:1704.05783].
 A. Denner, S. Dittmaier, and L. Hofer, Collier: a fortranbased Complex OneLoop LIbrary in Extended Regularizations, Comput. Phys. Commun. 212 (2017) 220–238, [arXiv:1604.06792].
 S. Höche, S. Kuttimalai, S. Schumann, and F. Siegert, Beyond Standard Model calculations with Sherpa, Eur. Phys. J. C75 (2015), no. 3 135, [arXiv:1412.6478].
 F. Krauss, D. Napoletano, and S. Schumann, Simulating associated production of and Higgs bosons with the SHERPA event generator, Phys. Rev. D95 (2017), no. 3 036012, [arXiv:1612.04640].
 J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, O. Mattelaer, H. S. Shao, T. Stelzer, P. Torrielli, and M. Zaro, The automated computation of treelevel and nexttoleading order differential cross sections, and their matching to parton shower simulations, JHEP 07 (2014) 079, [arXiv:1405.0301].
 T. Sj\unichar246strand, S. Ask, J. R. Christiansen, R. Corke, N. Desai, P. Ilten, S. Mrenna, S. Prestel, C. O. Rasmussen, and P. Z. Skands, An Introduction to PYTHIA 8.2, Comput. Phys. Commun. 191 (2015) 159–177, [arXiv:1410.3012].
 P. Artoisenet, R. Frederix, O. Mattelaer, and R. Rietkerk, Automatic spinentangled decays of heavy resonances in Monte Carlo simulations, JHEP 03 (2013) 015, [arXiv:1212.3460].
 J. Alwall et al., Comparative study of various algorithms for the merging of parton showers and matrix elements in hadronic collisions, Eur. Phys. J. C53 (2008) 473–500, [arXiv:0706.2569].
 M. Dobbs and J. B. Hansen, The HepMC C++ Monte Carlo event record for High Energy Physics, Comput. Phys. Commun. 134 (2001) 41–46.
 M. Cacciari, G. P. Salam, and G. Soyez, The Antik(t) jet clustering algorithm, JHEP 04 (2008) 063, [arXiv:0802.1189].
 M. Cacciari, G. P. Salam, and G. Soyez, The Catchment Area of Jets, JHEP 04 (2008) 005, [arXiv:0802.1188].
 M. Cacciari and G. P. Salam, Pileup subtraction using jet areas, Phys. Lett. B659 (2008) 119–126, [arXiv:0707.1378].
 D. Krohn, J. Thaler, and L.T. Wang, Jet Trimming, JHEP 02 (2010) 084, [arXiv:0912.1342].
 ATLAS Collaboration, G. Aad et al., Performance of jet substructure techniques for large jets in protonproton collisions at = 7 TeV using the ATLAS detector, JHEP 09 (2013) 076, [arXiv:1306.4945].
 J. M. Butterworth, B. E. Cox, and J. R. Forshaw, scattering at the CERN LHC, Phys. Rev. D65 (2002) 096014, [hepph/0201098].
 S. Höche, F. Krauss, M. Schönherr, and F. Siegert, NLO matrix elements and truncated showers, JHEP 08 (2011) 123, [arXiv:1009.1127].
 E. Bothmann, M. Sch\unichar246nherr, and S. Schumann, Reweighting QCD matrixelement and partonshower calculations, Eur. Phys. J. C76 (2016), no. 11 590, [arXiv:1606.08753].
 M. Czakon, D. Heymes, and A. Mitov, Bump hunting in LHC events, Phys. Rev. D94 (2016), no. 11 114033, [arXiv:1608.00765].
 A. V. Manohar and M. B. Wise, Flavor changing neutral currents, an extended scalar sector, and the Higgs production rate at the CERN LHC, Phys. Rev. D74 (2006) 035009, [hepph/0606172].
 E. Farhi and L. Susskind, Technicolor, Phys. Rept. 74 (1981) 277.
 A. Belyaev, G. Cacciapaglia, H. Cai, G. Ferretti, T. Flacke, A. Parolini, and H. Serodio, Diboson signatures as Standard Candles for Partial Compositeness, JHEP 01 (2017) 094, [arXiv:1610.06591].
 G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Rebelo, M. Sher, and J. P. Silva, Theory and phenomenology of twoHiggsdoublet models, Phys. Rept. 516 (2012) 1–102, [arXiv:1106.0034].
 A. Hayreter and G. Valencia, LHC constraints on color octet scalars, arXiv:1703.04164.
 ATLAS Collaboration, G. Aad et al., Search for new phenomena with photon+jet events in protonproton collisions at TeV with the ATLAS detector, JHEP 03 (2016) 041, [arXiv:1512.05910].
 ATLAS Collaboration, M. Aaboud et al., Search for new phenomena in highmass final states with a photon and a jet from collisions at = 13 TeV with the ATLAS detector, arXiv:1709.10440.
 A. Alloul, N. D. Christensen, C. Degrande, C. Duhr, and B. Fuks, FeynRules 2.0  A complete toolbox for treelevel phenomenology, Comput. Phys. Commun. 185 (2014) 2250–2300, [arXiv:1310.1921].
 C. Degrande, C. Duhr, B. Fuks, D. Grellscheid, O. Mattelaer, et al., UFO  The Universal FeynRules Output, Comput.Phys.Commun. 183 (2012) 1201–1214, [arXiv:1108.2040].
 J. Alwall, M. Herquet, F. Maltoni, O. Mattelaer, and T. Stelzer, MadGraph 5 : Going Beyond, JHEP 1106 (2011) 128, [arXiv:1106.0522].
 M. Czakon, D. Heymes, A. Mitov, D. Pagani, I. Tsinikos, and M. Zaro, Toppair production at the LHC through NNLO QCD and NLO EW, arXiv:1705.04105.
 K. J. F. Gaemers and F. Hoogeveen, Higgs Production and Decay Into Heavy Flavors With the Gluon Fusion Mechanism, Phys. Lett. B146 (1984) 347.
 D. Dicus, A. Stange, and S. Willenbrock, Higgs decay to top quarks at hadron colliders, Phys. Lett. B333 (1994) 126–131, [hepph/9404359].
 R. Frederix and F. Maltoni, Top pair invariant mass distribution: A Window on new physics, JHEP 01 (2009) 047, [arXiv:0712.2355].
 S. Gori, I.W. Kim, N. R. Shah, and K. M. Zurek, Closing the Wedge: Search Strategies for Extended Higgs Sectors with Heavy Flavor Final States, Phys. Rev. D93 (2016), no. 7 075038, [arXiv:1602.02782].
 S. Jung, J. Song, and Y. W. Yoon, Dip or nothingness of a Higgs resonance from the interference with a complex phase, Phys. Rev. D92 (2015), no. 5 055009, [arXiv:1505.00291].
 A. Djouadi, J. Ellis, and J. Quevillon, Interference effects in the decays of spinzero resonances into and , JHEP 07 (2016) 105, [arXiv:1605.00542].
 W. Bernreuther, P. Galler, C. Mellein, Z. G. Si, and P. Uwer, Production of heavy Higgs bosons and decay into top quarks at the LHC, Phys. Rev. D93 (2016), no. 3 034032, [arXiv:1511.05584].
 B. Hespel, F. Maltoni, and E. Vryonidou, Signal background interference effects in heavy scalar production and decay to a topantitop pair, JHEP 10 (2016) 016, [arXiv:1606.04149].
 D. Buarque Franzosi, E. Vryonidou, and C. Zhang, Scalar production and decay to top quarks including interference effects at NLO in QCD in an EFT approach, JHEP 10 (2017) 096, [arXiv:1707.06760].
 ATLAS Collaboration, M. Aaboud et al., Search for heavy Higgs bosons decaying to a top quark pair in collisions at TeV with the ATLAS detector, arXiv:1707.06025.
 P. Fileviez Perez, R. Gavin, T. McElmurry, and F. Petriello, Grand Unification and Light ColorOctet Scalars at the LHC, Phys. Rev. D78 (2008) 115017, [arXiv:0809.2106].
 I. V. Frolov, M. V. Martynov, and A. D. Smirnov, Resonance contribution of scalar color octet to production at the LHC in the minimal fourcolor quark\unichar8211lepton symmetry model, Mod. Phys. Lett. A31 (2016), no. 39 1650224, [arXiv:1610.08409].
 M. Gerbush, T. J. Khoo, D. J. Phalen, A. Pierce, and D. TuckerSmith, Coloroctet scalars at the CERN LHC, Phys. Rev. D77 (2008) 095003, [arXiv:0710.3133].
 C. Kim and T. Mehen, Color Octet Scalar Bound States at the LHC, Phys. Rev. D79 (2009) 035011, [arXiv:0812.0307].
 S. Schumann, A. Renaud, and D. Zerwas, Hadronically decaying coloradjoint scalars at the LHC, JHEP 09 (2011) 074, [arXiv:1108.2957].
 ATLAS Collaboration, M. Aaboud et al., Measurements of topquark pair differential crosssections in the lepton+jets channel in collisions at =13 TeV using the ATLAS detector, arXiv:1708.00727.
 CMS Collaboration, S. Chatrchyan et al., Search for resonant production in lepton+jets events in collisions at TeV, JHEP 12 (2012) 015, [arXiv:1209.4397].
 CMS Collaboration, A. M. Sirunyan et al., Search for resonances in highlyboosted lepton+jets and fully hadronic final states in protonproton collisions at 13 TeV, arXiv:1704.03366.
 ATLAS Collaboration, M. Aaboud et al., Search for new phenomena in dijet events using 37 fb of collision data collected at 13 TeV with the ATLAS detector, Phys. Rev. D96 (2017), no. 5 052004, [arXiv:1703.09127].
 ATLAS Collaboration, M. Aaboud et al., Search for new phenomena in highmass diphoton final states using 37 fb of proton–proton collisions collected at TeV with the ATLAS detector, Phys. Lett. B775 (2017) 105–125, [arXiv:1707.04147].
 C. Anastasiou, C. Duhr, F. Dulat, E. Furlan, T. Gehrmann, F. Herzog, A. Lazopoulos, and B. Mistlberger, CPeven scalar boson production via gluon fusion at the LHC, arXiv:1605.05761.
 T. Alanne, M. T. Frandsen, and D. Buarque Franzosi, Testing a dynamical origin of Standard Model fermion masses, Phys. Rev. D94 (2016) 071703, [arXiv:1607.01440].