]>PLB34358S0370-2693(19)30001-210.1016/j.physletb.2018.12.064The AuthorsPhenomenologyFig. 1Signal cross section as a function of ma at 13 TeV and fa = 1 TeV for N = 1, 5, 10, computed with MadGraph 5 with matching up to one extra jet after imposing pTleading-jet>300GeV at the generation level. The cross section includes an R-axion mass dependent K-factor correction [16].Fig. 1Fig. 2Comparison of the mSD distribution and the invariant mass distribution corrected by a factor of 1.5 for the SM processes W + jets and Z + jets.Fig. 2Fig. 3Left panel: 95% C.L. exclusion limits for CMS data in the [ma,N] plane. Right panel: 95% C.L. exclusion limits of the signal cross section for CMS data in terms of ma. This cross section has been computed with MadGraph 5 with matching up to one extra jet after imposing pTleading-jet>300GeV at the generation level and it includes an R-axion mass dependent K-factor correction [16]. In both cases the value of fa is set to 1 TeV.Fig. 3Fig. 495% C.L. exclusion limits for CMS data in the [ma,fa] plane for N = 1, 5, 10.Fig. 4Fig. 595% C.L. exclusion limits for CMS data in the [N,fa] plane for ma = 45, 165, 305 GeV.Fig. 5Fig. 6p0 value for CMS data as a function of ma for fa = 1 TeV.Fig. 6Fig. 795% C.L. exclusion limits for the gluon-fusion production cross section at the LHC of an ALP decaying into a jet pair as a function of its mass.Fig. 7Constraining R-axion models through dijet searches at the LHCErnestoArgandaab⁎ernesto.arganda@fisica.unlp.edu.arAníbal D.Medinaa⁎anibal.medina@fisica.unlp.edu.arNicolás I.Mileoa⁎mileo@fisica.unlp.edu.arRoberto A.Moralesa⁎roberto.morales@fisica.unlp.edu.arAlejandroSzynkmana⁎szynkman@fisica.unlp.edu.araIFLP, CONICET, Dpto. de Física, Universidad Nacional de La Plata, C.C. 67, 1900 La Plata, ArgentinaIFLPCONICETDpto. de FísicaUniversidad Nacional de La PlataC.C. 67La Plata1900ArgentinabInstituto de Física Teórica UAM/CSIC, Calle Nicolás Cabrera 13-15, Cantoblanco E-28049 Madrid, SpainInstituto de Física Teórica UAM/CSICCalle Nicolás Cabrera 13-15CantoblancoMadridE-28049Spain⁎Corresponding authors.Editor: A. RingwaldAbstractThe search at hadron colliders for new massive resonances of a few 100 GeVs that couple effectively to colored states is an extremely challenging issue, due principally to the presence of large QCD multijet backgrounds at this energy, rendering the searches at the LHC particularly difficult. Recently, it was realized that these large backgrounds could be overcome by demanding one high-pT jet from initial-state radiation (ISR) and by means of novel jet-reconstruction techniques through which the resulting hadronized products of the massive resonances are reconstructed as a fat-jet, a unique large-radius jet. The ATLAS and CMS Collaborations have recently reported searches for the experimental signature of a single fat-jet in association with an ISR jet. Models of dynamical supersymmetry breaking with an spontaneously broken R-symmetry give rise to the appearance of a pseudo-Nambu–Goldstone boson called the R-axion, which naturally tends to be light. In the parameter space regions where the anomalous R-axion coupling to gluons is boosted, these models can be tested against these new LHC dijet searches. Taking into account the CMS search, we apply the qμ statistical method to the signal events against the background-only expectation and obtain the 95% C.L. exclusion limits on the most relevant model parameters for a particular messenger sector, namely, the R-axion mass ma, the decay constant fa, and the number of color messengers N, being these limits suitable to be applied to more general models with axion-like particles.1IntroductionModels of dynamical supersymmetry (SUSY) breaking provide, via dimensional transmutation mweak∼e−O(1)8π2/g2MPlanck, an elegant solution to the problem of hierarchies between the a priori unrelated scales of weak (mweak≈O(100) GeV) and gravitational interactions (MPlanck≈O(1019) GeV). It was shown in [1] that in a generic class of dynamical models in which SUSY is broken via F-terms, an spontaneously broken R-symmetry is a sufficient condition for dynamical SUSY breaking, thus making the study of SUSY models with a U(1)R symmetry very appealing. The spontaneous breaking of the global R-symmetry leads to the appearance of a Nambu–Goldstone boson called the R-axion. There have been several studies about the nature and phenomenology of the R-axion [2–4]. The cancellation of the cosmological constant provided by the tuning of a constant term in the supergravity (SUGRA) superpotential leads to an explicit breaking of the U(1)R and thus to an unavoidable SUGRA contribution to the R-axion mass, making the R-axion a pseudo-Nambu–Goldstone boson (pNGB). There can be however other explicit sources of U(1)R breaking in the hidden sector that add up to the SUGRA contribution, but keep it light with respect to the other SUSY particles of the theory, respecting the pNGB nature and therefore possibly making it the first sign of SUSY. The couplings of the R-axion to other particles depend on the U(1)R charge assignments one makes on the different sectors of the theory (messenger and visible sectors). Nevertheless, if the sources of explicit U(1)R breaking are kept small, the R-axion couplings still display its Goldstone nature and thus are suppressed by its decay constant fa, which is related to the scale of U(1)R breaking.In [4] a phenomenological study of the R-axion was made where its decay constant and mass were kept as free parameters. The study focused on the possible collider signatures, and assumed that there is a strongly coupled hidden sector where SUSY is broken and from which SUSY breaking is transmitted to the visible sector à la gauge mediation. The U(1)R is realized non-linearly by the R-axion and SUSY is realized non-linearly in the constrained superfield formalism to capture the low-energy behavior, providing the interactions between the R-axion and the effective Gravitino Lagrangian. In particular, the coupling of the R-axion to gravitino pairs11The coupling is dominated by the Goldstinos which are longitudinal modes of the gravitino. can be quite sizable in certain regions of parameter space. Given that Majorana gaugino masses break explicitly the U(1)R, one expects a SUSY spectrum where gaugino-like neutralinos (and possibly Higgsinos) tend to be lighter than the other soft-breaking masses. Thus it is possible to obtain a somewhat natural spectrum at low energies, where the only light particles in the spectrum related to SUSY are the gravitino, the R-axion and possibly the neutralinos.In this work we propose to test a class of R-axion models against the dijet searches recorded recently by the CMS Collaboration [5] at the Large Hadron Collider (LHC) with a center-of-mass energy of s=13TeV and a total integrated luminosity of L=35.9fb−1. The experimental signature consists of a single massive large-radius jet in association with a jet from initial-state radiation, and the collaboration reports a slight deviation with respect the background-only expectation of 2.9σ of local significance and 2.2σ of global significance. The ATLAS Collaboration has performed a more recent search [6] for new resonances identified as massive large-radius jets consistent with a particle decaying into quark pairs, obtaining lower values for the local and global significances of the observed deviations above background. We will make use of some of these experimental results to probe the mentioned models of dynamical SUSY breaking in which the anomalous R-axion coupling to gluons is boosted, imposing constraints on the most relevant model parameters, namely, the R-axion mass ma, the decay constant fa, and the number of color messengers N. In addition, these limits should not be difficult to generalized to more generic axion-like particle (ALP) models that posses a similar coupling structure to the R-axion models exhibited here. Analyses of this type have been carried out in [7] recasting the limits imposed on Z′ models reported by CMS in [5].This paper is organized as follows: In Section 2 we provide a brief theoretical review on the R-axion and its most relevant interactions. In Section 3 we comment on the recent dijet searches performed by the ATLAS and CMS collaborations, while Section 4 is devoted to determine the constraints imposed by these experimental searches on the R-axion model. Finally, we conclude in Section 5.2Brief theoretical review on R-axionWe refer the reader to [3] and [4] for a more thorough description of the framework. We define R-symmetry as the largest subgroup of the automorphism group of the supersymmetry algebra that commutes with the Lorentz group. We focus on the Minimal Supersymmetric extension of the SM (MSSM) and as we mentioned in the Introduction, we expect a natural spectrum where squarks and sleptons of all families are decoupled from the low-energy effective theory. In fact, we consider that the only SUSY particles that remain in the low-energy spectrum are the R-axion and possibly the gravitino.22The gravitino tends to be light in the case that the F-term responsible for SUSY breaking is small. However, it has no impact on the phenomenology we are interested in. The relevant interactions of the R-axion for the phenomenology we wish to describe are given by its coupling to the MSSM gauge sector and in particular its anomalous coupling to gluons and photons. We assume that there exist two classes of messenger fields that transmit SUSY breaking to the visible sector. One is a single 5+5¯ of SU(5) which contributes to the masses of gauginos and sfermions in the usual general gauge-mediated way. The other class consists of N-copies of messengers q+q¯ which are 3+3¯ under SU(3)c and singlets under SU(2)L and U(1)Y weak gauge groups. This messenger sector upsets the unification of gauge couplings at high energies but that is something we are not concerned with in this work. It also implies that gluinos and squarks naturally tend to be heavier than weak gauginos and sleptons. This kind of spectrum is where the current searches at the LHC seem to be leading us towards given the lack of evidence of sparticles, in particular colored ones like gluinos and squarks which should be easily produced at a hadron collider. Another reason to consider this messenger sector lies in the fact that current diphoton resonant searches [8,9] already put strong constraints in the R-axion parameter space if one decides to consider the somewhat standard choice of N copies of 5+5¯ under SU(5) [4]. Allowing to split the messenger sector as we do implies that R-axions now become basically insensitive to diphoton searches but nonetheless, as we will show in later sections, they can still be probed via the novel fat-jet techniques that profit from the anomalous coupling to gluons. Following [3], we can write the relevant couplings of the R-axion to massless gauge bosons as(1)ΔLaggeff=gs232π2afa×[−N−1+3−2×cos2β2]GμνG˜μν,(2)ΔLaγγeff=e232π2afa×[−2+2−3×49×2cos2β]FμνF˜μν, where Gμν and Fμν are the gluon and photon field strengths, G˜μν=ϵμνρσGρσ/2 is the dual gluon field strength and similarly for the photon, and we have explicitly separated the contribution of the single 5+5¯ (the −1 and −2 in the couplings to gluons and photons, respectively) from that of the N-copies of 3+3¯ of SU(3)c which only enter in the coupling to gluons. Furthermore, we have assumed that all of the gauginos are heavy enough such that their loop contribution can be taken as an anomalous contribution as well (the +3 in the coupling to gluons and the +2 in the coupling to photons). The only SM fermion contribution that we are taking into account (given that we plan to consider an R-axion with a mass ma≈O(100) GeV) comes from the top quark, which enters in an indirect way given that we assume no direct coupling of the R-axion to fermions. In fact, we consider an R-symmetry consistent with the μ-term and Yukawa interactions, rHu+rHd=2, rHu+rQ+rU=2, rHd+rQ+rD=2, rHd+rL+rE=2, and in which the R-charges of Hu and Hd are fixed by the condition that the Goldstone boson associated to the Z-gauge boson is invariant under U(1)R, (rHu=2cos2β, rHd=2sin2β with tanβ=〈Hu〉/〈Hd〉). It is clear then that there is no direct coupling of the R-axion to fermions. Instead what happens is that the R-axion mixes with the CP-odd Higgs via the Bμ term in the Higgs potential, and this dictates the way in which the R-axion effectively couples to fermions, via the mixing and the CP-odd Higgs coupling to fermions. We will follow [3], and assume that the Bμ term is generated via Renormalization Group running from gaugino-loops. In that case, the effective coupling of the R-axion to SM fermions takes the form,(3)ΔLaff¯eff=iafa(cos2βmuu¯γ5u+sin2βmdd¯γ5d+sin2βmll¯γ5l). Thus we understand the appearance of the top quark in the anomalous coupling of the R-axion to gluons and photons as it enters in the top quark fermion loop triangle diagram. For example, in the anomalous contribution to gluons, the factor of 1/2 stems from the trace over the SU(3)c generators and the additional factor of 2 comes from the different directions in the fermion loop. In the anomalous coupling to photons, there is a color factor Nc=3, the top-quark electromagnetic charge squared and once again a factor of two from the two possible directions of the fermion loop. The factor cos2β/fa which is common in both diagrams comes from the effective coupling to up-type fermions, Eq. (3).We would like to stress that our results could be easily accommodated to any other different set of R-charges, for example one in which the μ-term carries non-vanishing R-charge when generated by the hidden sector dynamics, which would then modify rHu+rHd and with it all the other choices for the fermion R-charges. Similarly, as was done in [4], one could assume that the hidden dynamics is responsible for the generation of the Bμ term, in which case the effective coupling of the R-axion to fermions would be modified. Given that we plan to consider values of tanβ≫1 rendering the top-quark contribution marginal, and that we also consider R-axion masses such that the light fermion contributions to the anomalies are negligible, we expect that this latter modification on the mixing to be completely negligible.Therefore, the phenomenology in which we are interested in is described by the anomalous coupling to gluons and photons in Eqs. (1) and (2), as well as the coupling to fermions in Eq. (3). The relevant parameters are the number of color messengers, N, the R-axion mass, ma, the decay constant, fa, and tanβ. Throughout all the analysis performed in this paper we will vary N and ma in the ranges [1,10] and [45 GeV, 305 GeV], respectively, and consider fa=1 TeV as a benchmark. Finally, we will fix the value of tanβ to 10, though our results do not depend on this particular choice. In fact, we expect that our conclusions will not change for 4≲tanβ≲30.3Low mass dijet resonances at the LHCThe search at hadron colliders for low-lying-massive resonances that couple to colored states is extremely challenging due to the presence of large SM backgrounds, principally the so-called QCD multijet background. This issue can be faced by demanding one high-pT jet from initial-state radiation, which allows to satisfy the energy trigger requirements. Under these conditions, the CMS Collaboration has recently considered this class of search with an enough high pT for the new resonances [5], in such a way that the resulting hadronized products are reconstructed as an unique large-radius jet, usually called fat-jet. CMS has looked for narrow resonance peaks in the continuous distributions of fat-jet masses, taking advantage of a jet-mass and pT-decorrelated substructure variable which keeps the jet-mass distribution shapes. This novel variable is the jet-mass distribution groomed with the soft-drop algorithm [10,11], mSD, which reduces the distribution of jet masses of the QCD background arising from soft ISR gluons whilst keeps practically unchanged the jet masses coming from the new massive resonances and W/Z backgrounds. The fat-jet candidate in each event is the most energetic jet that fulfills the requirement to be an AK8 jet (reconstructed using the anti-kT algorithm with a radius parameter of 0.8) with pT>500GeV and |η|<2.5. The distribution for the collected data, corresponding to L=35.9 fb−1 at s=13TeV, and the simulated SM background of the leading-jet soft-drop mass mSD is displayed in Figure 1 of [5], while in Fig. 6 the mSD distribution is shown for data and background in five different pT ranges (500, 600, 700, 800, 900, and 1000 GeV) considered by CMS for a signal consisting of a leptophobic vector boson Z′ with a mass of 135 GeV. After computing the upper limits on the production cross section by means of the CLS method [12–14], a maximum local (global) excess of 2.9 (2.2) standard deviations is observed at mSD≃115GeV.More recently, the ATLAS Collaboration has reported a similar search [6], corresponding to L=36.1fb−1 at also s=13TeV. In this case, ATLAS uses two non-exclusive jet categories defined by the radius parameter R of the anti-kT algorithm: large-R jets (J) with R=1.0 and |η|<2.0 and narrow jets (j) with R=0.4 and |η|<2.4. The event selection implies the requirement of having at least one large-R jet with pTJ>450GeV (the resonance candidate) and at least one narrow jet (the ISR jet) with pTj>420GeV and azimuthal angular separation of Δϕ>π/2 with respect to the resonance candidate. The large-R jet mass distributions of data and background, not displayed in pT bins as in the CMS analysis, are presented in Figure 3 and 4 of [6], where it can be seen that the estimated background contributions reproduce very well the observed distributions. After a signal-plus-background fit to the large-R jet mass distribution with a Z′ model assumption, a local (global) excess of 2.5σ (1.1σ) is observed for masses around 150 GeV.The R-axion model described in Section 2 can actually lead to a final state signature with a fat-jet plus an ISR jet via the process pp→a(→jj)j. Although both the ATLAS and CMS searches can lead to constraints on the R-axion model, from now on we will only focus on the latter. On the one hand, the mass range explored by the CMS search (50–300 GeV) is wider than the one considered in the ATLAS analysis (100–220 GeV). On the other hand, the reconstruction of the two non-exclusive jet categories used by ATLAS is more difficult to deal with by means of a fast simulation of the detector response, and this would make our recasting of the ATLAS results less reliable than in the case of CMS. In addition, in the mass range where ATLAS and CMS analyses overlap, the exclusion limits on the signal cross sections are of the same order and therefore, we do not expect that our results would be significantly modified by the inclusion of the ATLAS analysis.In Fig. 1 we show the signal cross section σ(pp→a(→jj)j)=σ(pp→aj)×BR(a→jj) at 13 TeV for three values of N computed with MadGraph 5 [15] with matching up to one extra jet after imposing a cut on the transverse momentum of the leading jet pTleading-jet>300GeV at the generator level.33This cut is imposed in order to make a more efficient generation while considering a possible fraction of events generated with pTleading-jet<500GeV passing the cut pTleading-jet>500GeV applied in the CMS analysis at the reconstruction level. In addition, an R-axion mass dependent K-factor, which includes full NNLO, approximate N3LO and a threshold resummation at N3LL' [16], has been applied to the signal cross section according to Fig. 7 of that reference. Note that this figure provides the K-factor down to masses of 100 GeV and then, in order to cover the region between 45 GeV and 100 GeV, we have performed an extrapolation by adding a K-factor of 3.7 for ma=45GeV in consonance with [7]. As can be seen from Fig. 1, within the considered values for N and ma the cross section spans the range of ∼10−3−10pb. A comment is in order about the particular case in which N=2. From Eq. (1) we see that for this value there is a cancellation that results in an anomalous coupling to gluons proportional to rt=cos2β∼10−2 (since we use tanβ=10). This suppression in the coupling gives rise to a signal cross section many orders of magnitude smaller than the one obtained for the remaining values of N.4ResultsIn order to determine the constraints imposed by the CMS searches on the R-axion model, it is necessary to obtain the mSD distribution corresponding to our signal. However, the soft-drop algorithm required to reconstruct this variable is out of the scope of this work. Instead we used directly the invariant mass distribution of the leading jet corrected by a factor ϵ obtained from the comparison between the number of events at the peak of the mSD and the invariant mass distributions corresponding to the W+jets background. Here we are assuming that the resulting correction factor can also be used to reproduce the resonance peak in the mSD distribution of the signal. As it is stated in [5], the jet mass for merged W→qq¯ (as well as for Z→qq¯ and Z′→qq¯) jets is mostly unchanged by the soft-drop grooming since it arises from the kinematics of the decay process as opposed to QCD background jets. Therefore, expecting the same to happen with a→jj, it is not unlikely to assume a similar correction factor for our signal.We generated 4.5×105 events of the W+jets background at LO using MadGraph 5 [15], with the W boson decaying hadronically. The parton shower and hadronization were carried out with PYTHIA 8 [17] and the detector response was implemented with Delphes 3 [18]. Also, we matched the sample up to three additional jets. Regarding the jet reconstruction, we use the anti-kT algorithm with the R parameter set to 0.8. Hence, all the jets in the event are reconstructed as fat-jets and we associate the fat-jet candidate with the most energetic AK8 jet.In order to obtain the correction factor mentioned above, we plot with MadAnalysis5 [19–22] the invariant mass distribution of the leading jet after applying the cuts pT>500 GeV and |η|<2.5 used in [5]. From this simulated distribution we found that the number of events at the peak exceeds that of the mSD distribution by a factor of ϵ=1.5. In order to test the robustness of this value, we also performed the estimation of the correction factor by using the Z+jets background. In this case we obtain a value ∼1.6, whose difference with the factor derived from the W+jets background is negligible in terms of its impact on the exclusion limits. Based on this check, we used the same correction factor regardless of the resonance mass. In Fig. 2 we show the mSD distribution for the W+jets and the Z+jets backgrounds provided by CMS in [5] along with our estimation through the invariant mass distribution and a correction factor of 1.5 in the vicinity of the resonance peak. It is important to stress that even when ten bins of mSD are displayed in this figure, we concentrate, when estimating the correction factor, on the two bins with the highest number of events. By looking at Fig. 2 one may think that a correction factor ∼3 would be more appropriate than 1.5 for the region below the resonance peak. However, even overestimating this possible change in the correction factor outside the resonance peak and using the value 3 instead of 1.5 to correct the whole distribution, the exclusion limits are slightly weakened in half of the considered mass range. Finally, potential variations of the correction factor in the high invariant mass region above the resonance peak do not affect our results since the invariant mass distribution of the signal falls abruptly there.Regarding the simulation of the signal, we use the same setup as for the W/Z+jets backgrounds and consider values of N in the range [1,10] and resonance masses between 45 GeV and 305 GeV (as in [5]) in steps of 20 GeV. For each point of this grid we generate 5×104 events and build the leading jet invariant mass distribution in the five different pT ranges considered in the CMS analysis (see Figure 6 of [5]). These invariant mass distributions are in turn mapped into mSD distributions by means of the correction factor ϵ=1.5. The mSD distributions corresponding to the different backgrounds were taken from the Figure 6 of [5]. In addition to the correction factor, we need to apply the acceptance corresponding to the cut N21,DDT<0 that is included in the CMS analysis to discriminate two-prong signal jets from multijet QCD background jets preserving the shape of the soft-drop jet mass distribution. From Figure 3 of [5], the acceptance is found to be ∼0.3 both for the W+jets background and for the Z′ model tested in that reference.For purposes of establishing the exclusion limits on the parameter space of the R-axion model, we use the test statistic qμ, which is based on the profile likelihood ratio, and its corresponding p-value, denoted as pμ [14]. Notice that the five mutually exclusive pT regions were taken into account in the computation of the likelihood function. The region in the [ma,N] plane excluded at 95% C.L. by the CMS data is displayed on the left panel of Fig. 3 for fa=1TeV. We show the same limits translated to the signal cross section on the right panel. As can be seen from this figure, values of the number of color messengers N equal or greater than 7 are forbidden at 95% C.L. within the considered range for ma. In the region around ma=115GeV, where CMS has reported a slight deviation from the SM expectation, values of N≥6 are excluded. In order to quantify the impact of a potential systematic uncertainty in our estimate of the factor used to correct the invariant mass distributions, we recast the exclusion limits by considering different values of ϵ both above and below the estimated value (1.5). We found that the exclusion limits are not affected at all by variations of about 10%, while deviations of up to 30% shift the contour only in three mass bins out of the 14 considered in the analysis.It is interesting to reinterpret the exclusion limits in the [ma,fa] plane with the number of color messengers fixed. In Fig. 4 we show the region excluded at 95% C.L. by the CMS data for three different values of N. In contrast to Fig. 3, now the lower region of the plane [ma,fa] is excluded due to the fact that the couplings are inversely proportional to the decay constant. Notice that the exclusion limits become stronger as the value of N increases. In addition, only the range fa>4TeV is allowed for any of the values of ma and N considered along this work.Finally, we display in Fig. 5 the results in the [N,fa] plane for ma=45,165, and 305 GeV. Due to the particular suppression that occurs in the coupling of the R-axion to gluons for N=2, no constraint can be put in this case and, consequently, all the three exclusion regions depicted in Fig. 5 collapse to fa≃0 at N=2. On the other hand, the linear behavior of the exclusion region boundary obtained in the [N,fa] plane can be read off directly from the dependence of the anomalous coupling to gluons on the parameters fa and N. The fact that the exclusion region for ma=165 GeV is wider than for ma=305 GeV is not surprising but consistent with what we previously shown in Fig. 3, where it can be seen that the limit set to the signal cross section is more stringent in the former case than in the latter.As described in Section 3, CMS has observed slight local and global excesses from the background-only hypothesis at mSD≃115 GeV. Although the reported values are not statistically significant, we studied the level of agreement between the data and the background-only hypothesis in the case of the R-axion model. To quantify this, we used the statistic q0 and computed the corresponding p-value, denoted usually as p0 [14]. In Fig. 6 we display the p0 values as a function of the R-axion mass for the decay constant fa=1 TeV. The minimum value of p0, ∼0.09, is obtained at ma=117.5 GeV (N=5) and corresponds to a significance of ∼1.35 standard deviations, which is below the significances obtained within the context of the Z′ model considered in the CMS search. This is mainly caused by the discrete nature, through the parameter N, of the R-axion model, which constrains the optimization of the likelihood function.44Notice that if one allows free variations of fa, larger significances are possible.Lastly, we would like to comment that though we worked in the particular context of the R-axion model with a set of R-charges, the constraints derived from the specific CMS collider search can be easily generalized to other ALP models in which the ALP coupling to photons and SM fermions are suppressed, while the coupling to gluons is enhanced. In fact, this particular CMS search has been considered in [7] in order to put bounds on ALP models in the mass range between 10 GeV and 65 GeV. As a matter of fact, if one were to consider as main motivation for an axion particle the strong CP problem, there is in principle no reason to expect an anomalous coupling to photons given that there is no CP problem for QED. This coupling to photons is usually induced by the fermion electroweak charges to which the axion couples to. However, there is not necessarily in principle a reason for which these fermions should have EW charges in order to solve this QCD anomalous coupling. In that sense, an ALP like the one we are considering could be related more in the spirit to the strong CP problem. In Fig. 7 we reinterpret the upper bounds derived previously in terms of the R-axion production cross section via gluon fusion,55This cross section has been computed with MadGraph with only basic cuts. in such a way that they can be easily generalized to ALP models. For mass values below ∼150 GeV, the cross sections allowed by data varies from a few of tens of pb to ∼103 pb, whilst for larger mass values, the largest allowed cross section is ∼40 pb.5ConclusionsIn this paper we have made use of the results reported by the CMS collaboration in the search for dijet resonances at the LHC with the purpose of imposing constraints on a class of models of dynamical SUSY breaking with a boosted anomalous coupling of the R-axion to gluons and a negligible coupling to photons. From a phenomenological point of view, the most relevant parameters of this R-axion model are the number of color messengers, the mass of the R-axion and its constant decay. In order to establish exclusion limits on these parameters we have used the test statistic qμ and its corresponding p-value pμ, based on the profile likelihood ratio. Setting the decay constant to fa=1TeV, values of the number of color messengers N greater or equal to 7 are forbidden at 95% C.L. for R-axion masses within the considered range, [45 GeV, 305 GeV]. In the vicinity of ma=115GeV, where CMS has reported a slight excess from the background-only expectation, values of N≥ 6 are excluded. On the other hand, only values of the R-axion decay constant fa above 4 TeV are allowed by the CMS data for any N and the full ma range considered here. Finally, to quantify the level of agreement between the CMS data and the background-only hypothesis in the case of the R-axion model, we have calculated the corresponding p-value p0, obtaining a minimum of ∼0.09 at ma=117.5GeV that corresponds to a ∼1.35σ significance.Last but not least, it is important to mention that the results obtained with this analysis can be easily generalized to more generic ALP models with a similar coupling structure to the R-axion models studied here.AcknowledgementsWe thank Aurelio Juste and Carlos Wagner for useful discussions. E.A. warmly thanks IFT of Madrid for its hospitality hosting him during the completion of this work. This work has been partially supported by CONICET and ANPCyT project no. PICT 2016-0164.References[1]A.E.NelsonN.SeibergNucl. Phys. B41619944610.1016/0550-3213(94)90577-0arXiv:hep-ph/9309299[2]J.BaggerE.PoppitzL.RandallNucl. Phys. B4261994310.1016/0550-3213(94)90123-6arXiv:hep-ph/9405345[3]H.S.GohM.IbeJ. High Energy Phys.0903200904910.1088/1126-6708/2009/03/049arXiv:0810.5773 [hep-ph][4]B.BellazziniA.MariottiD.RedigoloF.SalaJ.SerraPhys. Rev. Lett.11914201714180410.1103/PhysRevLett.119.141804arXiv:1702.02152 [hep-ph][5]A.M.SirunyanCMS CollaborationJ. High Energy Phys.1801201809710.1007/JHEP01(2018)097arXiv:1710.00159 [hep-ex][6]M.AaboudATLAS CollaborationarXiv:1801.08769 [hep-ex][7]A.MariottiD.RedigoloF.SalaK.Tobiokahttps://doi.org/10.1016/j.physletb.2018.06.039arXiv:1710.01743 [hep-ph][8]G.AadATLAS CollaborationPhys. Rev. Lett.11317201417180110.1103/PhysRevLett.113.171801arXiv:1407.6583 [hep-ex][9]V.KhachatryanCMS CollaborationPhys. Lett. B767201714710.1016/j.physletb.2017.01.027arXiv:1609.02507 [hep-ex][10]M.DasguptaA.FregosoS.MarzaniG.P.SalamJ. High Energy Phys.1309201302910.1007/JHEP09(2013)029arXiv:1307.0007 [hep-ph][11]A.J.LarkoskiS.MarzaniG.SoyezJ.ThalerJ. High Energy Phys.1405201414610.1007/JHEP05(2014)146arXiv:1402.2657 [hep-ph][12]T.JunkNucl. Instrum. Methods A434199943510.1016/S0168-9002(99)00498-2arXiv:hep-ex/9902006[13]A.L.ReadJ. Phys. G282002269310.1088/0954-3899/28/10/313[14]G.CowanK.CranmerE.GrossO.VitellsEur. Phys. J. C712011155410.1140/epjc/s10052-011-1554-0Erratum:Eur. Phys. J. C732013250110.1140/epjc/s10052-013-2501-zarXiv:1007.1727 [physics.data-an][15]J.AlwallJ. High Energy Phys.1407201407910.1007/JHEP07(2014)079arXiv:1405.0301 [hep-ph][16]T.AhmedM.BonviniM.C.KumarP.MathewsN.RanaV.RavindranL.RottoliEur. Phys. J. C7612201666310.1140/epjc/s10052-016-4510-1arXiv:1606.00837 [hep-ph][17]T.SjöstrandComput. Phys. Commun.191201515910.1016/j.cpc.2015.01.024arXiv:1410.3012 [hep-ph][18]J.de FavereauDELPHES 3 CollaborationJ. High Energy Phys.1402201405710.1007/JHEP02(2014)057arXiv:1307.6346 [hep-ex][19]E.ConteB.FuksG.SerretComput. Phys. Commun.184201322210.1016/j.cpc.2012.09.009arXiv:1206.1599 [hep-ph][20]E.ConteB.DumontB.FuksC.WymantEur. Phys. J. C74102014310310.1140/epjc/s10052-014-3103-0arXiv:1405.3982 [hep-ph][21]B.DumontEur. Phys. J. C75220155610.1140/epjc/s10052-014-3242-3arXiv:1407.3278 [hep-ph][22]E.ConteB.FuksarXiv:1808.00480 [hep-ph]