Jet-associated resonance spectroscopy

TL;DR

This paper develops a model-independent framework for analyzing potential resonances in jet-photon and jet-Z final states at hadron colliders, aiding the identification of new physics signals.

Contribution

It provides a comprehensive parameterization of resonance amplitudes up to spin two, enabling detailed simulation of kinematic distributions for coupling analysis.

Findings

Framework constrains resonance spin and helicity couplings.

Generalizes existing signal generation strategies for collider experiments.

Offers insights into resonance properties despite limited P/CP determination.

Abstract

We present a model-independent study aimed at characterizing the nature of possible resonances in the jet-photon or jet-$Z$ final state at hadron colliders. Such resonances are expected in many models of compositeness and would be a clear indication of new physics. At leading order, in the narrow width approximation, the matrix elements are parameterized by just a few constants describing the coupling of the various helicities to the resonance. We present the full structure of such amplitudes up to spin two and use them to simulate relevant kinematic distributions that could serve to constrain the coupling structure. This also generalizes the signal generation strategy that is currently pursued by ATLAS and CMS to the most general case in the considered channels. While the determination of the P/CP properties of the interaction seems to be out of reach within this framework, there is a wealth of information to be gained about the spin of the resonance and the relative couplings of the helicities.

Figures

Tables

Table 1. Table 1: Central value cross sections for the γ 𝛾 \gamma +jet scenarios discussed in Sec.  3 for coupling values of 10 − 1 superscript 10 1 10^{-1} . The results for P/CP-violating parameters are identical.

scenario	$\sigma \left(M=500\text{GeV}\right)$ [pb]	${\Gamma }_X$ [GeV]	BR$\left(X\to \gamma j\right)$
$a_0^{gg}\otimes a_0^{\gamma g}$	83.58	0.36	0.55
$a_0^{q\overline{q}}\otimes a_0^{\gamma g}$	48.83	0.60	0.33
$a_{1/2}^{gq}\otimes a_{1/2}^{\gamma q}$	82.54	0.93	0.43
$a_1^{gg}\otimes a_1^{\gamma g}$	127.2	0.16	0.40
$a_1^{q\overline{q}}\otimes a_1^{\gamma g}$	22.22	0.20	0.33
$b_1^{q\overline{q}}\otimes a_1^{\gamma g}$	25.45	0.20	0.33
$a_{3/2}^{gq}\otimes a_{3/2}^{\gamma q}$	214.9	0.46	0.43
$a_2^{gg}\otimes a_2^{\gamma g}$	68.71	0.07	0.55

Table 2. Table 2: Central value cross sections for the Z 𝑍 Z +jet scenarios discussed in Sec.  3 for coupling values of 10 − 1 superscript 10 1 10^{-1} . The results for P/CP-violating parameters are identical.

scenario	$\sigma \left(M=500\text{GeV}\right)$	${\Gamma }_X$ [GeV]	BR$\left(X\to Zj\right)$
$a_0^{gg}\otimes a_0^{Zg}$	66.00	0.36	0.54
$a_{1/2}^{gq}\otimes a_{1/2}^{Zq}$	314.0	0.92	0.43
$a_1^{gg}\otimes a_1^{Zg}$	105.9	0.17	0.40
$a_1^{gg}\otimes c_1^{Zg}$	30.70	0.11	0.08
$a_{3/2}^{gq}\otimes a_{3/2}^{Zq}$	74.03	0.46	0.42
$a_2^{gg}\otimes a_2^{Zg}$	53.97	0.07	0.55

Equations

$$\gamma g\to X:$$

$$gg\to X:$$

$$q\bar{q}\to X:$$

$$\gamma q\to X:$$

$$\gamma \bar{q}\to \bar{X}:$$

$$gq\to X:$$

$$g\bar{q}\to \bar{X}:$$

$$\gamma g\to X:$$

$$gg\to X:$$

$$q\bar{q}\to X:$$

$$\gamma q\to X:$$

$$+\frac{\sqrt{3}}{2}{\tilde{a}}_{3/2}^{\gamma q}({\epsilon }_{\mu }{\bar{U}}_{\nu }{\gamma }^{\mu }{\gamma }^{5}u{q}^{\nu }/M+2{\epsilon }_{\mu }{q}^{\mu }{\bar{U}}_{\nu }{\gamma }^{5}u{q}^{\nu }/M^2)+i{\tilde{b}}_{3/2}^{\gamma q}({\epsilon }^{\mu }{\bar{U}}_{\mu }{\gamma }^{5}u-\frac{1}{2}{\epsilon }_{\mu }{\bar{U}}_{\nu }{\gamma }^{\mu }{\gamma }^{5}u{q}^{\nu }/M)),$$

$$\gamma \bar{q}\to \bar{X}:$$

$$+\frac{\sqrt{3}}{2}{\tilde{a}}_{3/2}^{\gamma \bar{q}}({\epsilon }_{\mu }\bar{v}{\gamma }^{\mu }{\gamma }^5{V}_{\nu }{q}^{\nu }/M+2{\epsilon }_{\mu }{q}^{\mu }\bar{v}{\gamma }^5{V}_{\nu }{q}^{\nu }/M^2)+i{\tilde{b}}_{3/2}^{\gamma \bar{q}}({\epsilon }^{\mu }\bar{v}{\gamma }^5{V}_{\mu }-\frac{1}{2}{\epsilon }_{\mu }\bar{v}{\gamma }^{\mu }{\gamma }^5{V}_{\nu }{q}^{\nu }/M)),$$

$$gq\to X:$$

$$+\frac{\sqrt{3}}{2}{\tilde{a}}_{3/2}^{\gamma q}({\epsilon }_{\mu }{\bar{U}}_{\nu }{\gamma }^{\mu }{\gamma }^{5}u{q}^{\nu }/M+2{\epsilon }_{\mu }{q}^{\mu }{\bar{U}}_{\nu }{\gamma }^{5}u{q}^{\nu }/M^2)+i{\tilde{b}}_{3/2}^{\gamma q}({\epsilon }^{\mu }{\bar{U}}_{\mu }{\gamma }^{5}u-\frac{1}{2}{\epsilon }_{\mu }{\bar{U}}_{\nu }{\gamma }^{\mu }{\gamma }^{5}u{q}^{\nu }/M)),$$

$$g\bar{q}\to \bar{X}:$$

$$+\frac{\sqrt{3}}{2}{\tilde{a}}_{3/2}^{g\bar{q}}({\epsilon }_{\mu }\bar{v}{\gamma }^{\mu }{\gamma }^5{V}_{\nu }{q}^{\nu }/M+2{\epsilon }_{\mu }{q}^{\mu }\bar{v}{\gamma }^5{V}_{\nu }{q}^{\nu }/M^2)+i{\tilde{b}}_{3/2}^{g\bar{q}}({\epsilon }^{\mu }\bar{v}{\gamma }^5{V}_{\mu }-\frac{1}{2}{\epsilon }_{\mu }\bar{v}{\gamma }^{\mu }{\gamma }^5{V}_{\nu }{q}^{\nu }/M)).$$

$$\gamma g\to X:{\delta }^{AB}$$

$$+b_2^{\gamma g}\left(M{S}^{\mu \nu \ast }{\epsilon }_{\mu }^{(\gamma )}{\epsilon }_{\nu }^{(g)}+{\epsilon }_{\mu }^{(g)}{q}^{\mu }{S}^{\nu \rho \ast }{\epsilon }_{\nu }^{(\gamma )}{q}_{\rho }/M+{\epsilon }_{\mu }^{(\gamma )}{q}^{\mu }{S}^{\nu \rho \ast }{\epsilon }_{\nu }^{(g)}{q}_{\rho }/M-\frac{1}{2}{{\epsilon }^{(\gamma )}}^{\mu }{\epsilon }_{\mu }^{(g)}{S}_{\nu \rho }^{\ast }{q}^{\nu }{q}^{\rho }/M\right)$$

$$+\sqrt{\frac{3}{2}}{\tilde{a}}_2^{\gamma g}{\epsilon }_{\mu }^{(\gamma )}{\epsilon }_{\nu }^{(g)}{S}_{\rho \lambda }^{\ast }{q}^{\rho }{q}^{\lambda }{ϵ}^{\mu \nu \alpha \beta }{q}_{\alpha }{p}_{\beta }/M^3$$

$$+\frac{1}{2}{\tilde{b}}_2^{\gamma g}\left({\epsilon }_{\mu }^{(\gamma )}{\epsilon }_{\nu }^{(g)}{S}_{\gamma \rho }^{\ast }{q}^{\rho }{ϵ}^{\mu \nu \gamma \lambda }{q}_{\lambda }/M+{\epsilon }_{\mu }^{(\gamma )}{{\epsilon }^{(g)}}^{\rho }{S}_{\nu \rho }^{\ast }{ϵ}^{\mu \nu \alpha \beta }{q}_{\alpha }{p}_{\beta }/M+{\epsilon }_{\mu }^{(g)}{{\epsilon }^{(\gamma )}}^{\rho }{S}_{\nu \rho }^{\ast }{ϵ}^{\mu \nu \alpha \beta }{q}_{\alpha }{p}_{\beta }/M\right)),$$

$$gg\to X:d^{AB{B}^{\prime }}$$

$$+b_2^{gg}\left(M{S}^{\mu \nu \ast }{\epsilon }_{\mu }{\epsilon }_{\nu }^{\prime }+{\epsilon }_{\mu }^{\prime }{q}^{\mu }{S}^{\nu \rho \ast }{\epsilon }_{\nu }{q}_{\rho }/M+{\epsilon }_{\mu }{q}^{\mu }{S}^{\nu \rho \ast }{\epsilon }_{\nu }^{\prime }{q}_{\rho }/M-\frac{1}{2}{\epsilon }^{\mu }{\epsilon }_{\mu }^{\prime }{S}_{\nu \rho }^{\ast }{q}^{\nu }{q}^{\rho }/M\right)$$

$$+\sqrt{\frac{3}{2}}{\tilde{a}}_2^{gg}{\epsilon }_{\mu }{\epsilon }_{\nu }^{\prime }{S}_{\rho \lambda }^{\ast }{q}^{\rho }{q}^{\lambda }{ϵ}^{\mu \nu \alpha \beta }{q}_{\alpha }{p}_{\beta }/M^3)$$

$$+f^{AB{B}^{\prime }}$$

$$q\bar{q}\to X:T_{b^{\prime }}^{Ab}$$

$$\displaystyle+\tilde{b}^{q\bar{q}}_{2}S_{\mu\nu}^{*}q^{\nu}\bar{v}\gamma^{\mu}\gamma^{5}u/M\bigg{)}.$$

$${}_{\mathrm{out}}⟨\mathcal{P},\lambda ;{\mathcal{P}}^{\prime },{\lambda }^{\prime }|X,m{⟩}_{\mathrm{in}}=(-{)}^{j-m}{}_{\mathrm{out}}⟨\bar{X},-m|\bar{\mathcal{P}},-\lambda ;{\bar{\mathcal{P}}}^{\prime },-{\lambda }^{\prime }{⟩}_{\mathrm{in}},$$

$$Zg\to X:$$

$$Zq\to X:$$

$$\displaystyle+\frac{1}{\sqrt{2}}\tilde{a}^{Zq}_{1/2}((M^{2}-m_{Z}^{2})\zeta_{\mu}\bar{U}\gamma^{\mu}\gamma^{5}u+2M\zeta_{\mu}q^{\mu}\bar{U}\gamma^{5}u)/M^{2}+2\,\tilde{c}^{Zq}_{1/2}\,m_{Z}^{2}\zeta^{\mu}q_{\mu}\bar{U}\gamma^{5}u/M^{3}\Big{)},$$

$$Z\bar{q}\to \bar{X}:$$

$$\displaystyle+\frac{1}{\sqrt{2}}\tilde{a}^{Z\bar{q}}_{1/2}((M^{2}-m_{Z}^{2})\zeta_{\mu}\bar{v}\gamma^{\mu}\gamma^{5}V+2M\zeta_{\mu}q^{\mu}\bar{v}\gamma^{5}V)/M^{2}+2\,\tilde{c}^{Z\bar{q}}_{1/2}\,m_{Z}^{2}\zeta^{\mu}q_{\mu}\bar{v}\gamma^{5}V/M^{3}\Big{)}.$$

$$Zg\to X:{\delta }^{AB}$$

$$+2c_1^{Zg}{m}_Z^2\left((M^2-m_Z^2){\epsilon }^{\mu }{S}_{\mu }^{\ast }{\zeta }^{\nu }{q}_{\nu }+{\zeta }^{\mu }{q}_{\mu }{\epsilon }^{\nu }{q}_{\nu }{S}^{\ast \rho }{q}_{\rho }\right)/M^4$$

$$\displaystyle+\tilde{a}^{Zg}_{1}S^{*}_{\mu}q^{\mu}\zeta_{\nu}\varepsilon_{\rho}q_{\lambda}p_{\sigma}\epsilon^{\nu\rho\lambda\sigma}/M^{2}+2\,\tilde{c}^{Zg}_{1}\,m_{Z}^{2}\zeta^{\mu}q_{\mu}\varepsilon_{\nu}S^{*}_{\rho}q_{\lambda}p_{\sigma}\epsilon^{\nu\rho\lambda\sigma}/M^{4}\bigg{)}.$$

$$Zq\to X:$$

$$+b_{3/2}^{Zq}((M^2-m_Z^2{)}^2{\zeta }^{\mu }{\bar{U}}_{\mu }u-\frac{1}{2}M(M^2-m_Z^2){\zeta }_{\mu }{\bar{U}}_{\nu }{\gamma }^{\mu }u{q}^{\nu }+m_Z^2{\zeta }^{\mu }{q}_{\mu }{q}^{\nu }{\bar{U}}_{\nu }u)/M^4$$

$$+\sqrt{6}{c}_{3/2}^{Zq}{m}_Z^2{\zeta }^{\mu }{q}_{\mu }{q}^{\nu }{\bar{U}}_{\nu }u/M^4$$

$$+\frac{\sqrt{3}}{2}{\tilde{a}}_{3/2}^{Zq}((M^2-m_Z^2){\zeta }_{\mu }{\bar{U}}_{\nu }{\gamma }^{\mu }{\gamma }^{5}u{q}^{\nu }+2M{\zeta }_{\mu }{q}^{\mu }{\bar{U}}_{\nu }{\gamma }^{5}u{q}^{\nu })/M^3$$

$$\displaystyle+i\,\tilde{b}^{Zq}_{3/2}\Big{(}(M^{2}-m_{Z}^{2})^{2}\zeta^{\mu}\bar{U}_{\mu}\gamma^{5}u-\frac{1}{2}M(M^{2}-m_{Z}^{2})\zeta_{\mu}\bar{U}_{\nu}\gamma^{\mu}\gamma^{5}uq^{\nu}+m_{Z}^{2}\zeta^{\mu}q_{\mu}q^{\nu}\bar{U}_{\nu}\gamma^{5}u\Big{)}/M^{4}$$

$$\displaystyle+\sqrt{6}\,\tilde{c}^{Zq}_{3/2}\,m_{Z}^{2}\zeta^{\mu}q_{\mu}q^{\nu}\bar{U}_{\nu}\gamma^{5}u/M^{4}\Bigg{)},$$

$$Z\bar{q}\to \bar{X}:$$

$$+b_{3/2}^{Z\bar{q}}((M^2-m_Z^2{)}^2{\zeta }^{\mu }\bar{v}{V}_{\mu }+\frac{1}{2}M(M^2-m_Z^2){\zeta }_{\mu }\bar{v}{\gamma }^{\mu }{V}_{\nu }{q}^{\nu }+m_Z^2{\zeta }^{\mu }{q}_{\mu }\bar{v}{V}_{\nu }{q}^{\nu })/M^4$$

$$+\sqrt{6}{c}_{3/2}^{Z\bar{q}}{m}_Z^2{\zeta }^{\mu }{q}_{\mu }\bar{v}{V}_{\nu }{q}^{\nu }/M^4$$

$$+\frac{\sqrt{3}}{2}{\tilde{a}}_{3/2}^{Z\bar{q}}((M^2-m_Z^2){\zeta }_{\mu }\bar{v}{\gamma }^{\mu }{\gamma }^5{V}_{\nu }{q}^{\nu }+2M{\zeta }_{\mu }{q}^{\mu }\bar{v}{\gamma }^5{V}_{\nu }{q}^{\nu })/M^3$$

$$\displaystyle+i\,\tilde{b}^{Z\bar{q}}_{3/2}\Big{(}(M^{2}-m_{Z}^{2})^{2}\zeta^{\mu}\bar{v}\gamma^{5}V_{\mu}-\frac{1}{2}M(M^{2}-m_{Z}^{2})\zeta_{\mu}\bar{v}\gamma^{\mu}\gamma^{5}V_{\nu}q^{\nu}+m_{Z}^{2}\zeta^{\mu}q_{\mu}\bar{v}\gamma^{5}V_{\nu}q^{\nu}\Big{)}/M^{4}$$

$$\displaystyle+\sqrt{6}\,\tilde{c}^{Z\bar{q}}_{3/2}\,m_{Z}^{2}\zeta^{\mu}q_{\mu}\bar{v}\gamma^{5}V_{\nu}q^{\nu}/M^{4}\Bigg{)}.$$

$$Zg\to X:{\delta }^{AB}$$

Full text

IPPP/17/49

1 Introduction

Many scenarios of dynamics Beyond the Standard Model (BSM), built with the aim to ameliorate the hierarchy problem, predict the existence of new resonances at the TeV scale. Examples of such states abound in many different contexts such as vector-like confinement [1, 2], compositeness [3, 4, 5, 6, 7, 8, 9, 10], partial compositeness [11, 12, 13, 14, 15, 16] and excited quarks [17, 18, 19, 20, 21, 22].

If such states are indeed observed in the future at the LHC (or at a future hadron collider), the first step in obtaining a detailed picture of the underlying theory is a dedicated spectroscopy program targeting the nature of their couplings, spin and CP properties.

These new resonances do not necessarily have to be QCD singlets and, thus, the general spectroscopy program that has been pursued in conjunction to the Higgs discovery [23, 24, 25, 26, 27, 28, 29, 30], using the pioneering techniques of [31] needs to be augmented by considering jet-inclusive final states.

While di-jet analyses exist and are already used to constrain the presence of BSM physics, analyses of electro-weak bosons in association with jets have received less attention (but do exist as well, see e.g. [32, 33, 34, 35]). This is predominantly due to the fact that these channels are less common in established BSM scenarios and limits are typically dominated by cleaner $ZZ$ or $\gamma\gamma$ channels. However, particularly in the aforementioned scenarios of (partial) compositeness, these channels do provide important information about the couplings of a possible discovery. This motivates searches and a characterization program of jet-$\gamma$ and jet-$Z$ resonances (related by gauge invariance) as an important probe of BSM physics at the LHC. Conversely, the lack of an observation in these channels would help to restrict the parameter space for such models. In either case, a detailed understanding of the possible dynamical scenarios is needed to perform an efficient analysis.

From a phenomenological perspective, jet-associated resonances are appealing as they combine potentially large signals, due to the presence of colored particles, with the precision of having a highly energetic photon or final state leptons, thereby being well-covered by existing trigger requirements. Searches so far have not given any hint for such resonances. However, much more data is being collected and hopefully will soon be analyzed. We expect that studies of jet-$\gamma$ resonances will feature among the many channels in which searches for BSM physics are being carried out.

The purpose of this note is to present a model-independent leading-order analysis of the various possibilities of jet-$\gamma/Z$ interactions, aimed at extracting information about the spin of such resonances from the kinematic distribution of the outgoing photon or reconstructed $Z$ boson. (It is important to stress that EFT-based power-counting arguments might not be valid in case of a strongly interacting nature of such a resonance.)

In Sec. 2, we analyze the cases of spin $j=0,\frac{1}{2},1,\frac{3}{2},2$, and construct the general form for the model-independent amplitude at parton level. In Sec. 3, we bridge these amplitudes to the hadron level for chosen benchmark scenarios and point out the sensitivity that can be expected from a dedicated spectroscopy program in these channels.

2 Model-independent amplitudes

In this section we present the amplitudes relevant for the study of jet-$\gamma$ and jet-$Z$ resonances. In the narrow width approximation, at leading order, we can decompose the amplitude into $2\to 1$ (production) and $1\to 2$ (decay) on-shell processes, each characterized by a handful of coefficients coupling the different helicities.

For the jet-$\gamma$ case, these coefficients are denoted by $a^{\mathcal{P}\mathcal{P}^{\prime}}$ and $b^{\mathcal{P}\mathcal{P}^{\prime}}$, where $\mathcal{P}$ and $\mathcal{P}^{\prime}$ are the relevant “partons” and $b$ refers to the larger helicity component (below we will also label these coefficients with a subscript indicating the spin of the resonance). The P/CP properties are indicated by putting a tilde on those coefficients related to amplitudes containing a $\gamma^{5}$ or a $\epsilon$ tensor.

We retain the same notation for the jet-$Z$ case, so that amplitudes with the same $a$ or $b$ coefficients reduce to the previous ones in the $m_{Z}\to 0$ limit. The additional coefficients arising from the longitudinal modes of the $Z$ are denoted by $c^{\mathcal{P}\mathcal{P}^{\prime}}$ and the corresponding amplitudes vanish in the $m_{Z}\to 0$ limit.

Since all $2\to 1$ and $1\to 2$ amplitudes have the dimension of a mass, we divide by the appropriate power of the resonance mass $M$ to get the right overall dimension so that all coefficients $a,b,c$ are dimensionless. In the case where these coefficients arise from an effective field theory at the scale $\Lambda\gg M$, they will then scale by the appropriate powers of $M/\Lambda$ but we find it unnecessary to introduce an additional scale in the kinematics at this stage.

2.1 The jet$+\gamma$ case

We begin by looking at $2\to 1$ processes involving 2 incoming massless particles of spin $1\text{ or }\frac{1}{2}$ creating a resonance of spin $0,\frac{1}{2},1,\frac{3}{2}$ or $2$. These amplitudes are those of relevance to the production of the resonance but also, by using CPT, to the decay into a parton and a photon. In the next subsection we will consider the inclusion of the $Z$ boson.

For the process $\mathcal{P}_{1}\mathcal{P}_{2}\to X$ we introduce the partons’ on-shell four-momenta, generically denoted by $p_{1}$ and $p_{2}$, ($p_{1}^{2}=p_{2}^{2}=0$) as well as $p=p_{1}+p_{2}$ and $q=p_{1}-p_{2}$, obeying $p^{2}=-q^{2}=2p_{1}\cdot p_{2}=M^{2}$. The gluon and photon four-dimensional polarization vectors are transverse and the quark spinors obey the massless Dirac equation with the associated momenta.

As far as the polarizations of the resonance $X$ are concerned, we have $\not\!pU=MU$ and $\not\!pV=-MV$ for spin $\frac{1}{2}$; $p^{\mu}S_{\mu}=0$ for spin 1; $\not\!pU_{\mu}=MU_{\mu}$ and $\not\!pV_{\mu}=-MV_{\mu}$, $\gamma^{\mu}U_{\mu}=\gamma^{\mu}V_{\mu}=p^{\mu}U_{\mu}=p^{\mu}V_{\mu}=0$ for spin $\frac{3}{2}$; $S_{\mu\nu}=S_{\nu\mu}$, $S_{\mu}^{\mu}=0$ and $p^{\mu}S_{\mu\nu}=0$ for spin 2.

It is then straightforward to construct all combinations that are Lorentz invariant and obey the Ward identities. What is a bit more tedious is to eliminate all the linearly dependent combinations, particularly for the case containing the $\epsilon$-tensor where one needs to use Schouten’s identity. This computation is greatly simplified by going to the center-of-mass frame (cms) of the resonance.

We always assume flavor conservation, so the quark polarizations $u,v$ refer to the same quark flavors. The color indexes of the resonance are denoted by $A$ (octet) or $a$ (triplet/anti-triplet), those of the partons by $B$, $B^{\prime}$, $b$ or $b^{\prime}$. In the case of two gluons we denote by $\varepsilon$ and $\varepsilon^{\prime}$ their polarizations and we can use either the totally anti-symmetric $f^{ABB^{\prime}}$ or the totally symmetric $d^{ABB^{\prime}}$ to respect the overall Bose symmetry, rendering selections rules a-la’ Landau-Yang irrelevant in this case (a fact also mentioned in [36]). The amplitudes read:

Spin 0:

[TABLE]

Spin $\frac{1}{2}$:

[TABLE]

Spin 1:

[TABLE]

Spin $\frac{3}{2}$:

[TABLE]

Spin 2:

[TABLE]

The amplitudes involving two incoming gluons respect Bose symmetry under the exchange: $B\leftrightarrow B^{\prime}$, $\varepsilon\leftrightarrow\varepsilon^{\prime}$, $q\leftrightarrow-q$.

It is then straightforward to write down the non-zero production amplitudes ${}_{\mathrm{out}}\langle X,m|\mathcal{P},\lambda;\mathcal{P}^{\prime},\lambda^{\prime}\rangle_{\mathrm{in}}$ for the process where parton $\mathcal{P}$ with helicity $\lambda$ coming along the positive $\hat{z}$ axis ($\theta=0$) and parton $\mathcal{P}^{\prime}$ with helicity $\lambda^{\prime}$ coming along the negative $\hat{z}$ axis ($\theta=\pi$) create a resonance of spin $j$ and spin projection along $\hat{z}$ equal to $m=\lambda-\lambda^{\prime}$. As already mentioned, the coefficients $a$ refer to the pairs of lowest helicity (RR or LL) while $b$ refers to RL or LR, if present.

The CPT theorem

[TABLE]

allows then one to write down the $1\to 2$ amplitudes for the decay of the resonance $X$ in its rest frame into one parton $\mathcal{P}$ and a photon.

If one assumes the hermiticity of the effective interaction giving rise to the couplings (e.g. when the interaction arises by integrating out heavy mediators, as also discussed in [37]), one can assume ${}_{\mathrm{out}}\langle X,m|\mathcal{P},\lambda;\mathcal{P}^{\prime},\lambda^{\prime}\rangle_{\mathrm{in}}=\;_{\mathrm{out}}\langle\mathcal{P},\lambda;\mathcal{P}^{\prime},\lambda^{\prime}|X,m\rangle^{*}_{\mathrm{in}}$. The above assumption, combined with the CPT theorem, forces the coefficients of the integer-spin amplitudes to be real and those of the half-integer ones to be conjugate of each other. In the simulations of Sec. 3 we will always assume this to be the case.

2.2 The jet+$Z$ boson case

For the decay into a jet and a $Z$ boson we need to generalize some of the amplitudes in the previous section to include a massive vector boson. The notation is as before, but now $p_{1}^{\mu}$ is the momentum of the $Z$ boson, with $p_{1}^{2}=m_{Z}^{2}$. The other parton is still massless and $q=p_{1}-p_{2}$ and $p=p_{1}+p_{2}$ now obey $q^{2}=2m_{Z}^{2}-M^{2}$, $p^{2}=M^{2}$ and $p\cdot q=m_{Z}^{2}$. The energy of the massless parton in the center-of-mass frame is now $(M^{2}-m_{Z}^{2})/2M$, so the on-shell normalization of the quark wave-function has changed. All remaining polarizations are the same and, for the $Z$, we have an additional longitudinal polarization. In order not to confuse the gluon and the $Z$ polarization, we refer to the gluon as $\varepsilon^{\mu}$ and denote the $Z$ polarizations by $\zeta^{\mu}$.

We write down the production amplitudes for the process $Z\,\mathcal{P}\to X$. Of course, in this case one is really only interested in the conjugate process $X\to\,Z\mathcal{P}$, but we stick with this notation for ease of comparison with the previous formulas. The coefficients are chosen in such a way that for $m_{Z}\to 0$ these amplitudes reduce to the previous one with a photon.

Spin 0:

[TABLE]

Spin $\frac{1}{2}$:

[TABLE]

Spin 1:

[TABLE]

Spin $\frac{3}{2}$:

[TABLE]

Spin 2:

[TABLE]

Although the purpose of this work is to be as model-independent as possible, it is worth commenting on possible scenarios in which such couplings could arise. Historically, bosonic resonances of this type were considered in the context of Technicolor models [4] while fermionic resonances arose considering models of quark compositeness [20]. These particles can be pair produced with ordinary QCD strength and this puts model independent bounds on the low mass region of the spectrum, typically below one TeV. The single production modes, via gluon-gluon or quark-gluon fusion, of interest to us have a higher mass reach but are more model-dependent.

As far as bosonic resonances are concerned, while the original motivation from technicolor has greatly diminished due to the phenomenological difficulties of these models, recently there has been interest in the search for colored pseudo-Nambu-Goldstone bosons that arise in more recent models of partial compositeness [7, 16]. The production of such objects occurs via the anomaly and the cross section scales like $(K/f)^{2}$, where $K$ is the anomaly coefficient and $f$ the decay constant. For $K/f=1/\mbox{TeV}$ the production cross-section ranges from 0.1 pb for $M=500$ GeV to 0.1 fb for $M=3$ TeV at $\sqrt{s}=13$ TeV. Here the bounds are still attractive for LHC searches, allowing, for some models, masses around one TeV for values of the decay constant of 1 TeV.

As far fermionic resonances are concerned, the recent CMS search [35] using the model [22] sets a bound of 5.5 TeV on first and second family excited quarks for amplitude coefficient $a_{1/2}^{gq}$ $a_{1/2}^{\gamma q}$ of order $g_{s}$, $e$ respectively, and 1.8 TeV for excited $b$-quarks. The production cross section time branching ratio into jet-$\gamma$ are about 0.2 fb and 7 fb at the upper limit of the excluded mass range. For amplitude coefficients reduced by a factor 10, the mass reach is below 2 TeV for light quarks while the search lacks sufficient sensitivity to set a bound for the $b$-quark.

3 Elements of Hadron Collider Phenomenology

To gain a quantitative understanding of the phenomenology that we can expect at the LHC, we first consider the most motivated cases of scalars, spin half fermions and vector resonances. We later comment on the qualitative differences compared to the higher spin modes for representative examples, thus generalizing the analysis of [32, 33, 34, 35] to all allowed coupling structures in light of CPT and Lorentz invariance.

We have implemented the couplings of the previous section into the MadEvent [38, 39] event generator with purpose-built Helas routines [40], and have checked these against implementations derived from the FeynRules [41] and Ufo [42] toolkits. As already mentioned, we only consider flavor symmetric cases and treat all quarks in the four-flavor scheme (and hence the involved pdfs) on an equal footing ($q=u,d,c,s$ in the following). We interface the generated parton-level events with Herwig++ [43, 44] for showering and hadronization, and choose as benchmark $M=500~{}{\text{GeV}}$ for demonstration purposes. We will point out the influence of the mass scale $M$ on our analysis below, where we will also comment on detector resolution effects. Throughout, we focus on 13 TeV collisions.

Due to the very character of the above amplitudes being decompositions of physical scattering amplitudes in the narrow width approximation, we choose a small reference width $\Gamma_{X}/M=10^{-4}$. In this work we will not discuss the extraction of the coupling sizes comprehensively, but we note that our analysis of the signal events is entirely insensitive to the exact choice as long as $\Gamma_{X}/M\ll 1$. Since the amplitudes of Sec. 2 are not valid in the off-shell regime of any involved legs, we can expect measurements of the CP character of the resonance along the lines of [45, 46, 47, 48, 49, 50, 51], when relevant, to be significantly limited already at this stage. The non-validity of the above amplitudes for off-shell momenta also does not allow us to perform multi-jet matching for the signal.

We have simulated the contributing backgrounds in an identical way and specifically focus on QCD-induced $\gamma/Z$+jet production which are the largest contributing backgrounds in the SM. The LHC experiments typically estimate these with data-driven methods (see e.g. [32]) and this part of our analysis solely serves as numerical guide to highlight the potential sensitivity for our benchmark. To allow us to compare signal and background on an equal footing, we do not include jet-matching effects for the backgrounds.

Since the resonances we study in this work are motivated from general amplitude Lorentz structures, and since these carry colour charges, such a discovery can be established not only $\gamma/Z+\hbox{jet}$ production, but also in multi-jet final states. Analyses of the latter have been carried out in two jet [52, 53], as well as in four jet final states [54, 55]. A discovery will crucially depend on the sizes of the coefficients quoted in Secs. 2.1 and 2.2, as these can be chosen to avoid discovery in the multi-jet channels.

The LHC already performs bump hunts for resonances in the $\gamma+\hbox{jet}$ channel as detailed in [32, 33, 34, 35]. These are performed by a data-driven estimate of the sensitive invariant mass distribution with the aim to reveal excesses in model-independent approach. If such a search is successful, the question of the precise coupling structure of the new discovery arises. The $Z$ boson takes a special role in this case due to gauge invariance (however at a much smaller rate due to leptonic branching ratios which deliver the cleanest signatures for subsequent analyses). The imminent spectroscopy programme after such a discovery will then need to be informed by a range of searches, in particular because the exclusive rate of decays into the final states we discuss in this work will be influenced by the coefficients of the di-jet channels. For the purpose of this work we will assume a discovery in the jet+$\gamma$ channel at 300/fb (this might be accompanied by a similar observation in the multi-jet channels) and we outline a spectroscopy follow-up programme in the jet+$\gamma/Z$ channels. We will come back to the importance of multi-jet final states at the end of this section.

Different spin expectations can be discriminated by characteristic angular distributions [56, 57, 58, 18, 59, 60, 61]

(see also [29, 24, 62]). The angles that are sensitive to the spin structure of the interactions which are relevant to our final states are (Fig. 1)

[TABLE]

The subscripts $X,Z$ refer to the rest frames in which these angles are defined. The momenta are defined from the decay products, i.e.

[TABLE]

with

[TABLE]

in case of $\hbox{jet}+Z$ production. Note that the helicity angle $\theta_{h}$ and the azimuthal angle $\phi$ are not observable for decays $X\to\gamma j$, hence limiting the available range of sensitive observables to $\theta^{\ast}$. We have introduced $\phi$ for completeness but we find that it contains no discriminative power for the scenarios we study in this work, and we will not further consider this angle in the following.

The discriminating power of these angles (we will discuss the contributing backgrounds further below) lies in the fact that the boosts into respective rest frames remove some kinematic dependence on the final states and the mass of $X$ in particular. For the jet-$\gamma$ case, in the lab frame, the scattering is fully described through a combination of transverse momentum $p_{T}$ and pseudorapidity $\eta$ of the photon.

While removing this energy dependence has the benefit of projecting out the helicity decomposition of the interactions of Sec. 2 (sculpted to some extent by finite detector coverage), the distribution of the signal events according to energy-sensitive observables such as the transverse photon momentum can provide evidence of the dominant production mechanism via pdf effects. Another avenue we will discuss in the following is the prospects of quark-gluon tagging [63, 64, 65], which can in principle further discriminate the $X$ decay phenomenology, thus providing important information in discriminating the different spin hypotheses. This strategy is particularly motivated, as the threshold induced by the expected large mass of $X$ helps to choose working points that are particularly attractive to quark-gluon tagging [66].

Our event selection is performed fairly inclusively only reflecting the basic trigger thresholds for the signal events to be recorded. Representative generator level ccross sections for the coupling combinations that we study in this paper are tabulated in Tabs. 1 and 2.

For the case of photon-associated jet production we require an isolated photon (defined as isolated if the hadronic energy deposit in an area $R<0.3$ around the photon candidate is less than 5% of the photon’s transverse momentum) with

[TABLE]

to guarantee the event to be triggered [67]. In case of final state leptons, we require two isolated leptons (hadronic energy deposit less than 10% of the lepton candidate’s transverse momentum in $R<0.3$) with opposite charge and

[TABLE]

Again these criteria reflect the standard trigger thresholds [67]. On top of these thresholds we require

the leptons to be compatible with the $Z$ pole mass

[TABLE]

in the leptonic final state case.

The jets are clustered with the anti-kT algorithm [68] with resolution parameter $D=0.4$ using FastJet [69] and we define jets from the thresholds

[TABLE]

We require the leading jet in $p_{T}$ to be inside $|\eta_{j_{1}}|<2.33$ in a back-to-back configuration with the reconstructed photon or $Z$ boson ($p_{Z}=p_{\ell^{-}}+p_{\ell^{+}}$) in the azimuthal angle–pseudo-rapidity plane $R(j,Z/\gamma)>2.5$. Finally we require consistency of the reconstructed resonance with our mass hypothesis within 50 GeV

[TABLE]

This latter criterion, while not relevant for our signal distributions, is crucial for a comparison of the signal with the expected background. For comparison, ATLAS searches are sensitive to width/mass ratios of 2% in Ref. [32], which is well-covered by our representative invariant mass window cut in our signal-like selection, where our approximations can be trusted.

For the jet$+\gamma$ case, where we essentially only have a single angle at our disposal to discriminate the various hypotheses, we can already identify the qualitative overall behavior of the final state. The $p_{T}$ distribution of the photon (Figs. 2 and 3), while giving some indication of the dominant partonic subprocess through the parton distribution functions as well as spin and coupling character, is largely dominated by the threshold of the particle $X$, whose mass gets equally distributed into transverse momentum for central production. The differen-

tial measurement of $\theta^{\ast}$ (Figs. 3 and 4), on the other hand, allows us to formidably discriminate between different spin hypotheses. Some of the remaining qualitative degeneracies can be lifted (see below).111Note that the acceptance selections (which are determined by trigger thresholds and mostly finite detector coverage) bias the constructed distribution, see [62] for a detailed discussion in the context of QCD singlets.

The background distribution is fundamentally different from the signal distribution. As alluded

to above, none of the observables of the hard $2\to 2$ scattering reflects the P/CP character of the couplings and more involved processes that access off-shell information need to be considered (see Fig. 4).

All results in Figs. 2 – 4 do not include the expected effects of energy mis-measurements. As the sensitivity in these angle depends on the correct reconstruction of the $X$ rest frame, detector resolution effects can in principle limit the sensitivity. In practice, when the mass scales are large, the expected calibration of photons and jets has shown to be under extremely good control and is likely to improve even further in the future [70]. However, to provide a more quantitative understanding of how detector effects modify the angular distributions we show a comparison of $\cos\theta^{\star}$ including detector effects in Fig. 6 (we have adopted the energy parameterisation of Delphes3 [71]).

This phenomenological situation remains largely unchanged when considering the $Z$ boson case (Figs. 6 and  7) as far as the $p_{T}$ and $\eta$ distributions are concerned and we choose not to show these distributions for this reason. However, $\cos\theta^{\ast}$ remains an important additional handle to constrain and discriminate different scenarios (albeit at a lower rate due to the leptonic decays of the $Z$ boson that we consider here). The distribution of $\cos\theta_{h}$ essentially discriminates the decays of the transversely polarized $Z$ (coefficients $a$ and $b$) from the longitudinal one (coefficient $c$). See some examples in Fig. 7.

It is important to add that the background distribution is highly sculpted towards the bulk of signal hypotheses, making further discrimination beyond the aforementioned cases increasingly difficult.

An additional handle for disentangling the spin hypothesis is through discriminating the production modes. While the transverse momentum distributions of the $Z$ boson or the photon give some understanding of the dominant perturbative partonic subprocess, another avenue to further access this information is via identifying, at least approximately (see e.g. the recent Ref. [72]), the quark- or gluon-like character of the leading final state jet. This analysis step, which is entirely complementary to the analysis of angles $\theta_{h}$ and/or $\theta^{\ast}$, needs to be understood as an additional criterion that is invoked on a final selection that separates signal from background. The drop in signal rate before and after quark/gluon tagging is applied will provide additional power separating integer spin from the ${1\over 2}$ and ${3\over 2}$ hypotheses. This is shown representatively in Fig. 8 for the concrete example of separating the $b^{q\bar{q}}_{1}\otimes a^{\gamma g}_{1}$ coupling combination as null-hypothesis from the $a^{gq}_{1/2}\otimes a^{\gamma q}_{1/2}$ (alternative) hypothesis. The exclusion limits of Fig. 8 are calculated using the CLs method [73, 74, 75] for a binned log-likelihood based on $\cos\theta^{\ast}$ observable in the jet+$\gamma$ channel of Fig. 4, assuming “even” coupling structures.

This particular choice for a representative example is motivated from the overall similarities in the $\cos\theta^{\ast}$ observable, however, with a clear separation of quark- vs. gluon-initiated hard jet. As representative working point for quark-tagging and gluon-rejection, we use efficiencies $(\epsilon^{t}_{q},\epsilon^{r}_{g})=(0.5,0.13)$, which have been obtained in Ref. [66] under similar kinematical conditions. Throughout, we include the expected background distributions and the relative reduction of the SM $\mathrm{jet}+\gamma$ production after quark/gluon tagging. For the considered mass range of Eq. (11), the SM continuum production is dominated by processes with a final state quark and the reduction of background cross section is mostly determined by the tagger’s working point. Compared to the jet$+\gamma$ production cross section after selection cuts of $\sim 6.6~{}\text{pb}$, quark/gluon tagging reduces the cross section by $\sim 47\%$. In Fig. 8(b), we add the tagged distributions as a separate category to the likelihood, which can be compared to the “raw” $\cos\theta^{\ast}$ discrimination in Fig. 8(a). As can be seen, some of the competing distributions can be excluded to support statistical preference for one particular model for moderate luminosities. This simple hypothesis test which indicates statistical preference between two well-defined hypotheses does not constitute a coupling measurement, but will be the first step in this direction (experimental results related to the Higgs can be found in e.g. [76]).

Continuing with this particular example, at higher luminosity, once the $b^{q\bar{q}}_{1}\otimes a^{\gamma g}_{1}$ character is established, the overall rate in excess the background measurement can be used to constrain the couplings of the model. We show this representatively in Fig. 9 for our $M=500~{}\text{GeV}$ benchmark, showing all parameter combinations that yield a maximum signal cross section of the estimated 5$\sigma$ discovery cross section of $\sim 23.4~{}\text{fb}$. This MC-based toy extraction also demonstrates that additional information from di-jet measurements is necessary to avoid blind directions, which arise from fitting the narrow width approximation: For large values of e.g. $b^{q\bar{q}}_{1}$, the jet$+\gamma$ signal cross section scales as function of $a^{\gamma g}_{1}$ alone and the limit is saturated by a constraint on this single coupling.

A di-jet constraint, on the other hand, will close the blind direction for large values of $b^{q\bar{q}}_{1}$. A remaining question is whether the estimated 5$\sigma$ discovery cross section of $\sim 23.4~{}\text{fb}$ corresponds to a choice of parameters for which our main assumption, i.e. the narrow width approximation is still valid. Scanning the couplings inputing this discovery threshold as a constraint, we can obtain $b^{q\bar{q}}_{1}$ as a function of $a^{\gamma g}_{1}$, which then allows to express $\Gamma_{X}$ as function of $b^{q\bar{q}}_{1}$ alone. For our particular benchmark we obtain $\Gamma_{X}/M\lesssim 0.03$ for $|b^{q\bar{q}}_{1}|<1$, highlighting the good validity of the narrow width approximation in this case. The branching ratios for the parameter choices that are allowed this way vary between dominant $j\gamma$-like final states for $|b^{q\bar{q}}_{1}|\ll 1$ to di–jet like final states for $|b^{q\bar{q}}_{1}|\sim 1$, which again shows the necessity of studying complementary channels in case of a discovery in final states as described in this work.

4 Summary

The absence of conclusive hints for new interactions beyond the Standard Model motivates a wider approach to searches for new states that fall within the energy coverage of the LHC or future hadron colliders. Scenarios of QCD-charged new states that could arise in a range of composite Higgs models have been less investigated in the past, and this work presents a detailed investigation of possible $2\to 2$ scattering processes of such states with jet$+\gamma$ and jet$+Z$ production. This provides the theoretical underpinning for future searches, generalizing the current signal modelling of ATLAS and CMS [32, 33, 34, 35]. We have based our analysis around a decomposition of the scattering amplitudes into irreducible categories, thereby widening the phenomenological range beyond the constraints of effective field theory. This comes at the price of a strict on-shell formulation of the hard scattering amplitude, which removes a straightforward application of CP-discriminating techniques. However, the phenomenology of the scattering amplitudes proves rich and is observable at the LHC over a broad particle mass range. Naturally, the lack of any observation so far that could be interpreted along the lines of this paper leaves the parameter space vastly bigger than the humble set of angular and kinematic distributions of $2\to 2$ scattering is capable to constrain. However, we have shown that by adapting angles from the Higgs characterization program, at least partial insights can be gained into the spin of such a produced state. In particular, we have also demonstrated how advances in jet technology (specifically quark/gluon tagging) can assist the resonance’s discovery or spectroscopy in the future.

Acknowledgements

C.E. is supported by the IPPP Associateship scheme and by the UK Science and Technology Facilities Council (STFC) under grant ST/P000746/1. G.F. would like to thank C. Petersson for discussions. M.S. is supported in part by the European Commission through the “HiggsTools” Initial Training Network PITN-GA-2012-316704.