We study the sensitivity of top pair production with six-fermion decay at the LHC to the presence and nature of an underlying boson, accounting for full tree-level Standard Model interference, with all intermediate particles allowed off-shell. We concentrate on the lepton-plus-jets final state and simulate experimental conditions, including kinematic requirements and top quark pair reconstruction in the presence of missing transverse energy and combinatorial ambiguity in jet-top assignment. We focus on the differential mass spectra of the cross section and asymmetry observables, especially demonstrating the use of the latter in probing the coupling structure of a new neutral resonance, in addition to cases in which the asymmetry forms a complementary discovery observable.

## 1 Introduction

New fundamental, massive, neutral, spin-1 gauge bosons () appear ubiquitously in theories Beyond the Standard Model (BSM). The strongest limits for such a state generally exist for the and signatures, known collectively as Drell-Yan (DY). However, in addition to their importance in extracting the couplings to top quarks, resonance searches in the channel can offer additional handles on the properties of a due to uniquely available asymmetry observables, owing to the fact that (anti)tops decay prior to hadronisation and spin information is effectively transmitted to decay products. Their definition in , however, requires the reconstruction of the top quark pair. In these proceedings we summarise our study of the sensitivity to the presence of a single boson at the Large Hadron Collider (LHC) arising from a number of generationally universal benchmark models (section 2), as presented in our recently submitted paper [1]. We simulate top pair production and six-fermion decay mediated by a with full tree-level SM interference and all intermediate particles allowed off-shell, with analysis focused on the lepton-plus-jets final state, and imitating some experimental conditions at the parton level (section 3). We assess the prospect for an LHC analysis to profile a boson mediating production, using the cross section in combination with asymmetry observables, with results and conclusions in section 4 and 5, respectively.

## 2 Models

There are several candidates for a Grand Unified Theory (GUT), a hypothetical enlarged gauge symmetry, motivated by gauge coupling unification at approximately the GeV energy scale. often arise due to the residual U gauge symmetries after their spontaneous symmetry breaking to the familiar SM gauge structure. We study a number of benchmark examples of such models. These may be classified into three types: inspired models, generalised Left-Right (GLR) symmetric models and General Sequential Models (GSMs) [2].

One may propose that the gauge symmetry group at the GUT scale is E. When recovering the SM, two residual symmetries U and U emerge, which may survive down to the TeV scale. LR symmetric models introduce a new isospin group, SU, perfectly analogous to the SU group of the SM, but which acts on right-handed fields. This symmetry may arise naturally when breaking an SO gauge symmetry. We are particularly interested in the residual U and U symmetries, where the former is related to , and and refer to Baryon and Lepton number, respectively. An SSM has fermionic couplings identical to those of the SM boson, but is generically heavier. In the SM the couplings to fermions are uniquely determined by well defined eigenvalues of the and generators, the third isospin component and the Electro-Magnetic (EM) charge.

For each class we may take a general linear combination of the appropriate operators and fix , varying the angular parameter dictating the relative strengths of the component generators, until we recover interesting limits. These models are all universal, with the same coupling strength to each generation of fermion. Therefore, as with an SSM , the strongest experimental limits come from the DY channel. The limits for these models have been extracted based on DY results, at and TeV with an integrated luminosity of fb, from the CMS collaboration [3] by Accomando et al. [5], with general consensus that such a state is excluded below TeV.

## 3 Method

Measuring as the angle between the top and the incoming quark direction, in the parton centre of mass frame, we define the forward-backward asymmetry:

(1) |

With hadrons in the initial state, the quark direction is indeterminate. However, the is likely to carry a larger partonic momentum fraction than the in . Therefore, to define we choose the axis to lie along the boost direction. The top polarisation asymmetry (), measures the net polarisation of the (anti)top quark by subtracting events with positive and negative helicities:

(2) |

where () denote the eigenvalues under the helicity operator of (). Information about the top spin is preserved in the distribution of . We construct two dimensional histograms in and , and equate the gradient of a fitted straight line to .

In each of the models, the residual U gauge symmetry is broken around the TeV scale, resulting in a massive boson. This leads to an additional term in the low-energy Lagrangian, from which we may calculate the unique coupling structure for each observable:

(3) | ||||

(4) | ||||

(5) | ||||

(6) |

where and are the vector and axial-vector couplings of a specific fermion ().

While a parton-level analysis, we incorporate restraints encountered with reconstructed data, to assess, in a preliminary way, whether these observables survive. The collider signature for our process is a single or produced with at least four jets, in addition to missing transverse energy (). Experimentally, the -tagged jet charge is indeterminate and there is ambiguity in -jet (anti)top assignment. We solely identify with the transverse neutrino momentum. Assuming an on-shell we may find approximate solutions for the longitudinal component of the neutrino momentum as the roots of a quadratic equation. In order to reconstruct the event, we account for bottom-top assignment and solution selection simultaneously, using a chi-square-like test, by minimising the variable :

(7) |

where and are the invariant mass of the leptonic and hadronic (anti)top, respectively.

In order to characterise the sensitivity to each of these models, we test the null hypothesis, which includes only the known processes of the SM, assuming the alternative hypothesis (), which includes the SM processes with the addition of a single , using the profile Likelihood ratio as a test statistic, approximated using the large sample limit, as described in [9]. This method is fully general for any D histogram, and we test both D histograms in , and D in and the defining variable of each asymmetry to assess their combined significance.

## 4 Results

Figure 1 shows plots for the differential cross section, and . The statistical error is quantified for this luminosity assuming Poisson errors. The absent models, including all of the class, only produce an asymmetry via the interference term, which generally gives an undetectable enhancement with respect to the SM yield. The absence of a corresponding peak in either asymmetry offers an additional handle on diagnosing a discovered . The cross section, profiled in , shows a very visible peak for all models. The GSM models feature a greater peak, and width, consistent with their stronger couplings, but the impact on the cross section is otherwise similar for both classes. Mirroring the cross section, the distribution clearly distinguishes between the models and SM, with the difference in width even more readily apparent. The best distinguishing power over all the models investigated comes from the distribution, which features an oppositely signed peak for the GLR and GSM classes.

To evaluate the significance of each asymmetry as a combined discovery observable we bin in both and its defining variable. For , the asymmetry is calculated directly. Therefore, we divide the domain of into just two equal regions. is extracted from the gradient of the fit to for each mass slice, and we calculate the significance directly from this histogram. The final results of the likelihood-based test, as applied to each model, and tested against the SM, are presented in table 1. The models with non-trivial asymmetries consistently show an increased significance for the 2D histograms compared with using alone, illustrating their potential application in gathering evidence to herald the discovery of new physics.

Class | U | Significance () | ||

& | & | |||

E | U | - | - | |

U | - | - | ||

U | - | - | ||

U | - | - | ||

U | - | - | ||

U | - | - | ||

GLR | U | |||

U | - | - | ||

U | ||||

U | ||||

GSM | U | |||

U | ||||

U | - | - |

## 5 Conclusions

We have investigated the scope of the LHC in accessing semileptonic final states produced by pairs emerging from the decay of a heavy state. We tested a variety of BSM scenarios embedding one such a state, and show that asymmetry observables can be used to not only aid the diagnostic capabilities provided by the cross section, in identifying the nature of a possible signal, but also to increase the combined significance for first discovery. While the analysis was performed at the parton level, we have implemented a reconstruction procedure of the (anti)tops that closely mimics experimental conditions. We have, therefore, set the stage for a fully-fledged analysis eventually also to include parton-shower, hadronisation, and detector reconstruction, which will constitute the subject of a forthcoming publication. In short, we believe that our results represent a significant phenomenological advancement in proving that charge and spin asymmetry observables can have a strong impact in accessing and profiling signals during Run 2 of the LHC. This is all the more important in view of the fact that several BSM scenarios, chiefly those assigning a composite nature to the recently discovered Higgs boson, embed one or more state which are strongly coupled to top (anti)quarks [10].

## Acknowledgements

We acknowledge the support of ERC-CoG Horizon 2020, NPTEV-TQP2020 grant no. 648723, European Union. DM is supported by the NExT Institute and an ATLAS PhD Grant, awarded in February 2015. SM is supported in part by the NExT Institute and the STFC Consolidated Grant ST/L000296/1. FS is supported by the STFC Consolidated Grant ST/K001264/1. We would like to thank Ken Mimasu for all his prior work on phenomenology in , as well as his input when creating the generation tools used for this analysis. Thanks also go to Juri Fiaschi for helping us to validate our tools in the case of Drell-Yan production. Additionally, we are very grateful to Glen Cowan for discussions on the statistical procedure, and Lorenzo Moneta for aiding with the implementation.

## References

- [1] L. Cerrito, D. Millar, S. Moretti, and F. Spano, arXiv:1609.05540 (2016).
- [2] E. Accomando, A. Belyaev, L. Fedeli, S. F. King, and C. Shepherd-Themistocleous, Phys. Rev. D83, 075012 (2011).
- [3] CMS collaboration, JHEP 04, 025 (2015).
- [4] ATLAS collaboration, Phys. Rev. D90, 052005 (2014).
- [5] E. Accomando et al., JHEP 01, 127 (2016).
- [6] K. Hagiwara, J. Kanzaki, H. Murayama, and I. Watanabe, KEK-91-11 (1992).
- [7] T. Stelzer and W. F. Long, Comput. Phys. Commun. 81, 357 (1994).
- [8] G. P. Lepage, J. Comput. Phys. 27, 192 (1978).
- [9] G. Cowan, K. Cranmer, E. Gross, and O. Vitells, Eur. Phys. J. C71, 1554 (2011).
- [10] D. Barducci, S. De Curtis, K. Mimasu, and S. Moretti, Phys. Rev. D88, 074024 (2012).