{'fraction_non_alphanumeric': 0.02976364879382752, 'fraction_numerical': 0.010450239281179803, 'mean_word_length': 4.184430379746836, 'pattern_counts': {'":': 0, '<': 5, '<?xml version=': 0, '>': 11, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 6, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Mass-degenerate Heavy Vector Mesons at Hadron Colliders

We study the LHC phenomenology of a couple of mass-degenerate heavy new gauge bosons with the quantum numbers of the Z and photon. We give a leading-order estimate of the number of events expected in Drell-Yan processes in terms of the parameters of the model (mass and coupling) and of the LHC machine specifications (luminosity and energy). We consider the feasibility of measuring a forward-backward asymmetry for various choices of the parameters and estimate the potential reach. We comment on how the results may affect future collider design and the results for a specific model of electro-weak symmetry breaking by way of example.


Introduction
The main purpose of the LHC is to discover the origin of electro-weak (EW) symmetry breaking and any associated new physics at the TeV scale.
One clean signal for the LHC would be the discovery of new, higher-mass copies of the standard-model (SM) gauge fields. For such a scenario, the LHC reach in observable masses and couplings is of interest. There are theoretical motivations for such searches. Many weakly-coupled extensions of the SM gauge group yield new heavy vector fields [1]. Such new fields also appear in little Higgs models [2]. Even strongly-coupled extensions, such as technicolor (TC), predict the existence of towers of heavy spin-1 fields [3].
In TC models EW symmetry is broken by a condensate, along the lines of chiral symmetry breaking in QCD. Thus these models do not contain a light elementary Higgs particle ‡. As a result the most promising test of these models at the LHC is a search for higher-mass spin-1 resonances. In particular, a copy of the full SU(2)×U(1) realized as heavy gauge bosons would be distinctive. If the new bosons are nearly degenerate in mass then the model will satisfy the precision EW constraints [5,6], without suppressing the coupling to SM fermions. This situation would imply that the new boson resonances overlap at the LHC.
Once the LHC reaches its design specification, proton-proton collisions will scan over the energy range current data allow higher resonances to occur at. In generic models this is in the 1 TeV region, but in many concrete realizations the bound is above the 2 − 3 TeV range [7,8].
In this paper, we consider the scenario in which a complete set of four spin-1 copies of the SM gauge bosons with heavy degenerate masses are produced in pp collisions. We concentrate on detection of the neutral heavy-Z and heavy-γ through their decay into muons. The phenomenology of the degenerate scenario has been largely ignored to date in the literature. This is in contrast to the single Z ′ situation, see for example [9] and references therein. For simplicity, we model the decay width of these particles by assuming it to be dominated by the decay into SM fermions. This avoids complications (and model-dependencies) arising when including the decay into light scalars and gauge bosons. If the underlying dynamics is based on an SU(N c ) gauge theory with large values of N c then the boson to boson decays are suppressed by the N c counting.
At LHC energies the SM background is mostly hadronic, so by looking at the neutral sector of higher resonances we exploit the very low background for muon and electron pair production. Within these assumptions we can study the event rate to leading order. Our aim is to identify what portion of parameter space will be testable at the LHC as a function of integrated luminosity and energy.
In the SM, a forward-backward asymmetry arises from interference terms between Z and γ exchange diagrams. It is a sizable effect, even in the TeV energy range. Measuring such an effect at a hadron collider is onerous but a method and its use at the LHC was surveyed for instance by Dittmar [10]. It was shown that a deviation of the measured asymmetry from the SM prediction is a possible signature of the presence of a new, very broad resonance with the quantum numbers of a Z ′ .
In contrast we will study a relativity narrow resonance. In our case, near the pole of some new resonance well into the TeV region the number of events is well in excess of the SM. So much so that the study of the di-lepton invariant-mass distribution will contain a clear peak. This will be a suitable discovery signature. Furthermore for all practical purposes the SM contribution to this peak can be ignored. In this case, one expects no significant asymmetry to be measurable at the peak from the SM.
If such a narrow peak is discovered in the data, but it results from the existence of a γ ′ , in addition to the Z ′ , then the interference between these two new states yields a large asymmetry. Measuring this asymmetry is hence a clean way to identify the presence of two distinct particles in the peak region, with masses so close together that they cannot be resolved in the invariant-mass distributions [11].
The organization of our paper is the following. We begin by introducing a model that encompasses the phenomenology of many specific models of EWSB. We then summarize some known results about decay widths and cross-sections. Then we proceed to study the Drell-Yan production of spin-1 resonances at the LHC. To illustrate the dependence of the event rate on model parameters we make a series of approximations to give a simple formula for the total number of events near a peak. We provide estimates of the discovery reach for two collider energies, 10 and 24 TeV.
We then derive a simple formula for computing the reach of forward-backward asymmetry measurement at a collider. At a 10 TeV machine the ability to search for a forward-backward asymmetry, of the type we are interested in, is limited to already excluded parameter space. Therefore we present data at 14 TeV (which is close to the border of feasibility) and 24 TeV. To illustrate the results we finish by concentrating on a particular model. Individual models give relations between the free parameters in our model-independent analysis and restrict the parameter space. This allows us to evaluate a typical reach scenario. A few words of explanation about these choices are needed. The 10 TeV centre of mass energy analysis is representative of the LHC experimental program. As we will see, large integrated luminosities are needed for discovery and hence comparing to the 7 TeV -with low luminosity -initial phase of LHC running is not very useful for the cases we study. By 10 TeV energy, however, the reach is quite substantial. This is the reason for our conservative choice of energy. The 24 TeV energy analysis is somewhat speculative, but provides a way to gauge how an energy upgrade compares to a luminosity upgrade of the apparatus, in reference to the specific process we study.

Set-up and Rationale
This section reviews cross-section and asymmetry calculations as well as defining the notations used throughout the paper. Most of the material in this section is well-known, but it is useful to remind the reader about it, and to fix the notation used in later sections. We model the new neutral resonances by adding an extra set of SU(2) × U(1) gauge fields into the SM electro-weak Lagrangian. We assume that some mechanism spontaneously breaks the SU(2) 2 × U(1) 2 gauge symmetry to the SM SU(2) L × U(1) Y above the TeV scale, and that in the process a copy of the W , Z and γ acquire a common mass m. In general, the latter is not rigorously correct: the new heavy bosons would acquire different masses in any realization. However, the idea that the mass splitting is parametrically small is suggested by the precision electro-weak parameters. We will provide a simple illustration of this in the next sub-section. To parametrize the couplings of the new heavy gauge bosons to the SM currents we introduce a parameter R. This is defined so that if the coupling of a given SM boson to a given SM current is g SM then the corresponding new heavy boson has a coupling of √ Rg SM to the same current. For example if the SM Z coupling to an electron current is g Z then the corresponding new boson, the Z ′ , would have a coupling of √ Rg Z to the electron current. The common mass m and the universal ratio R of the couplings are the only new parameters in the model.
It must be noted that in the paper we focus our attention on the neutral vectors only, setting aside the charged W ′ . This choice is due to the fact that for the processes we consider, the kinematics of the leptonic final state can be reconstructed precisely when the intermediate state is neutral, while this is difficult for charged states, in which case one has either to deal with missing energy, or to consider final states with hadronic jets.

Mass Degeneracy
As previously noted, in this work we assume that the new heavy gauge fields are degenerate in mass. As an example, illustrating the rationale behind this assumption, we examine the precision parameterŜ [5,6]. Recall that on general grounds one can write the precision parameterŜ aŝ with g the SU(2) L coupling in the SM, f and m being the decay constant and mass of the vector (v) and axial (a) mesons. In models where vector dominance holds the lightest among the states dominate the sum. Denoting them states as the techni-ρ meson and the techni-a 1 meson we truncate the sum to these two states. Under this truncation, we write the parameters of these two states in terms of a constant a 1 and a mass splitting ∆m 2 , Working with these definitions to leading order in ∆m 2 /m 2 ρ we see that, Experimentally,Ŝ is found to be O(10 −3 ) [6]. This requires to implement some mechanism to obtain an acceptable value from eq.(4). Choosing a 1 = 1 would require some sort of fine-tuning, hence we assume that a 1 be some generic O(1) number. Another possibility, that f 2 ρ /m 2 ρ ≪ 1, is realized in QCD. However, this ratio scales as N c in the large-N c limit (for N c the number of colors of the confining gauge theory). This makes it an unsuitable mechanism in a large-N c theory. Also, allowing for f 2 ρ /m 2 ρ ≪ 1 to be present would mean that the (techni-)rho mesons coupling to the fermions is heavily suppressed (the parameter R is very small), and its main decay is into longitudinally polarized standard-model vector mesons. This is in analogy to what happens in QCD where the rho meson decays into pions. Finally, in this case the rho mesons are actually strongly coupled, in the sense that their width is large. For all of these reasons we prefer not to use this suppression mechanism. Another possibility to satisfy the electro-weak precision constraints is to assume that ∆m 2 /m 2 ρ ≪ 1, small enough thatŜ 10 −3 . Interestingly, this observation extends beyond the vector-dominance approximation, at least at large-N c , in the sense thatŜ is suppressed if the mass splitting between new axial and vector meson states is small. This is the regime which inspires this present work.

Parton-level processes
In this section we begin the computation of the relevant cross-section, following standard textbook procedures. Initially this is done for a lepton machine then the results are extended to a hadron machine. This choice is made to highlight the physics of the model considered.
We find it useful to introduce the function P V , where √ŝ is the total energy, in the centre of mass frame, for the collision whose crosssection we compute, m V is the mass of the particle V and Γ V is the decay width of the particle V for V = γ, Z, γ ′ , Z ′ . This function enters the propagator of the boson V . Note that this requires computing the decay widths for the two new resonances. The new heavy gauge bosons can also decay into SM gauge boson pairs (e.g Z ′ → W + W − ). However the rate of such processes is very model-dependent, and at the same time assuming it to dominate the decay width might imply that the new gauge bosons be strongly coupled, as anticipated. In this study we assume such channels are negligible. This is equivalent to setting the 3-point and 4-point meson vertices to zero. This assumption is justified by the 1/N c counting if the underlying dynamics is similar to that of a large-N c gauge theory and the heavy spin-1 states are mesons of the strongly coupled sector, in analogy with the ρ and a 1 of a (large-N c ) QCD-like theory.
Within these assumptions, the γ ′ and Z ′ have decay widths proportional to those of the γ and Z in the SM. The constant of proportionality is the ratio between the SM and new gauge boson couplings. This is the parameter R introduced earlier on. Denoting the SU (2) coupling of the standard model by g, the Weinberg angle by θ w , the fine structure constant by α (and the new heavy boson mass by m) the decay widths into SM fermions of the new states are, where we summed over all the SM fermions, including the top, neglecting the fermion masses. Typical values of the new heavy masses are O(1) TeV, thus the error introduced by neglecting fermion masses in the phase-space factors is negligible.
As an illustration of the model, we begin by removing the complications of pp collisions and take the case of e + e − → µ + µ − . There are four possible diagrams. In addition to the SM channels of γ and Z exchange there is also γ ′ and Z ′ exchange. Following standard notation, introduce c V and c A for the vector and axial couplings defined as where Q denotes the electric charge and T 3 the third component of weak isospin for the relevant fermion. We denote for Q i and Q f the charges of the initial and final fermions in the Feynman diagram. Using this notation the angular 2 → 2 partonic cross-section can be written out as, forθ the angle of inclination with respect to the beam andΩ the solid angle of the scattering process. The A and B are defined as, where ℜ denotes the real part.
Similarly we can write the forward-backward asymmetry A F B as We obtain, It is useful to combine the above expressions and rewrite the cross-section in terms of an angular function g(θ): with This form of the partonic cross-section, in terms of g(θ), is the most useful form for our discussion of asymmetry. Above an energy of 500 GeV the SM contribution becomes negligible for our purposes, as anticipated in the introduction. Events near the resonance of the new heavy states,ŝ ≃ m 2 , are dominated by the new states as long as mΓ V ′ ≪ŝ, since in this case Therefore unless R is very small the new resonances will dominate. The definitions of Π V given can now be approximated by dropping the P V terms.
We use this approximation to re-write the asymmetry near the new resonance (ŝ ≃ m 2 ) as This approximation is used in computing the numerical values of A F B that we work with later. This method avoids the complications introduced by accounting for the weak dependence of A F B onŝ . Now we quote the correct formulae in the case of a hadron collider. The formulae for asymmetry and cross-section at partonic level can be amended by adding a factor for the number of colors N c , and changing to quark couplings. Denoting the quark/lepton couplings with a q/l superscript, A is modified to Including the appropriate color factor, the partonic cross-section for qq → l + l − iŝ where the angular dependence has been integrated over. Adjusting B to account for quarks gives The asymmetry at the partonic level is hence In summary, our study of event rates focuses around the partonic cross-section eq.(24). Our asymmetry calculations makes use of the angular partonic cross-section The forward-backward asymmetry required in g(θ) is computed from eq.(26) under the approximation of no SM background. Explicitly this means that

Drell-Yan processes
We will follow the standard recipe to model processes at a hadron collider. We consider the incident protons to be composed of partons holding fraction ξ of the total proton momentum p. Numerous sets of parton distribution functions (PDF) are available through Fortran code [14]. For this paper we choose the Fermi02 PDF set [15]. For practical purposes, one can treat the proton and its constituent partons as massless.
The two beams of protons each have energy √ s/2. The energy of a parton-parton collisionŝ is given, in terms of the two fractions ξ 1 and ξ 2 , bŷ If, for a parton x, the PDF is denoted by f x (ξ) we can follow convention and write the pp → l + l − cross-section as, To fix the total energy √ s of the pp collision, we impose the condition of Eq. (29). This means we may define the parton luminosity L q to be and write the full cross-section summed over all quarks as, Finally we can write the total number of events N in some range ofŝ at a machine with integrated luminosity L, We allow √ s to vary between 10 and 28 TeV. The coupling ratio R will be allowed to vary between 0 and 1, where the gauge boson coupling is certainly perturbative. Then the mass range will be taken as lying between 1 and 4 TeV. Analysis will also require the assumption that decay widths are small compared to the resonance mass. Inspection of the decay widths, under the limits on R and m, shows this condition is satisfied. Asymmetry measurements at a hadron collider are made more difficult (in comparison with a leptonic machine) by the unknown directions of the initial quark and anti-quark. However, at high rapidity there is a much greater chance that the quark carries the bulk of the momentum. Thus if data near a resonance in the invariant-mass distribution are cut so that only events with a rapidity Y greater than some value Y c (|Y | > Y c ) are used, then there is a corresponding degree of certainty regarding the quark directions. This method is used and explained more fully by Dittmar [10]. In summary, the magnitude of the rapidity controls the probability of misidentification of the quark/anti-quark directions, while the sign of the rapidity represents whether the quark was from the left or right in the collider. This methodology of inserting cuts, Lorentz boosting with the rapidity modulus to the muon centre of mass frame and finally using the sign of Y to define the angle consistently allows one to infer an asymmetry. Such an asymmetry is the remnant of the partonic asymmetry.

Event Rate
In this section we introduce a set of approximations that are useful in order to study the dependence of observable quantities (such as number of expected events and forwardbackward asymmetry) on the parameters of the model and on the machine specifications. The numerical part of the analysis will, however, use the complete equations and PDF distributions without approximations. From the approximations made we gain insight into how event rates behave as a function of model parameters.
Using standard electro-weak relations the coupling g can be written in terms of α and sin 2 θ w ≃ 0.22. Hence, near the new resonance: Using Eqs. (23)- (24) we can comment on the relative strength of the γ ′ and Z ′ contributions to the total cross-section. Because c l v ≪ 1, the cross term in Eq. (24) is suppressed. Also, in the last term of this equation c f 2 V ≪ c f 2 A ≃ 1/4, and hence the contribution of the γ ′ is a factor of ∼ 16Q 2 q larger than the contribution coming from the Z ′ . All of this, together with the fact that the PDF for the up-quark and downquarks are the most important, and that the coupling of the up quarks to γ ′ is four times stronger than the coupling of down quarks, means that the most significant contribution to the total event rate is given by the subprocess uū → γ ′ → µ + µ − . The next largest contribution to the event rate is from the cross term between the new photon and Zlike states (which is approximately half the value of the leading contribution) and is responsible for the asymmetry we discuss later.
Focusing on the u-quark PDF, we may attempt to approximate it with some form of expansion. The only relevant dimensionless combinations areŝ/s and m 2 /s. As we are looking at the peak we can replace the square mass withŝ so that only one relevant dimensionless combination remains. From the definition of L the overall dimension is 1/s. If the mass of the resonance were to increase one would expect the event rate to decrease. It is also reasonable to think increasing the machine energy s will increase the event rate. We use an approximate form of, and consider various n so as to obey the above criteria and then numerically fit the constant of proportionality. The minimal n would be −2 to fulfill our expectations for the behavior of L . Viewing our expression as the leading term in some expansion a polynomial can be added to increase accuracy. We find that an acceptable approximation is, provided √ŝ ∈ [1,4] TeV. We then go through a numerical fit of the constants a 1 and a 2 . In the range √ s ∈ [10,28] TeV just those two terms are needed.
We are interested in the event rate N in the peak region Approximating L u with its value at the pole (which is allowed since the resonances are very narrow), and using yields where the SM background from γ and Z has been neglected. Finally, the contribution of up-quarks to the total rate is, up to a proportionality constant that contains the electro-weak couplings and the normalization of the PDF distribution.
To correct for the neglected contributions of dd events and Z ′ exchange we can consider the ratio of the found rate N u to a more complete expression involving uū and dd events with γ ′ and Z ′ exchange. For notational simplicity let us absorb the coupling coefficients into the definition of P V briefly and write the ratio η as, An expression similar to Eq. (36) describes, with different choices of the coefficients, L d . Then we find that a very good approximation is given by for some b i and c i that are combinations of the up and down coefficients. The total event rate can be written, up to an overall constant, as We fix the parameters, a 1 , a 2 , b 1 , b 2 , c 1 , c 2 plus an overall normalization with the direct numerical computation using the PDFs. We estimate the error hence introduced between the estimated event rate and the direct computation to be about 5%. Discovery of the new resonance requires detecting a minimum number of events N at a machine integrated luminosity L. The approximate event rate can be solved for Suppose one considers N = 10 events to be necessary for discovery of a new resonance. Setting N = 10, if one has data corresponding to a 40 fb −1 integrated luminosity, then the line N/L = 1/4 fb indicates the minimum coupling R required to see such an effect. In the case of a 10 TeV machine, for a new resonance with mass m = 3 TeV, then the minimum R is around 0.4. Therefore any new boson with coupling ratio R > 0.4 and m < 3 TeV will produce at least 10 events in the data set obtained with 40 fb −1 integrated luminosity.
Changing the machine energy from 10 TeV to 24 TeV means moving from Fig. 1 to Fig. 2. At 10 TeV, by improving the integrated luminosity by a factor of 20, one may move from a reach in masses below 2.2 TeV to roughly 3.3 TeV, provided R = 1 (see the leftmost and rightmost lines in Fig. 1). Whilst improving the collision energy takes the reach to well above 4.0 TeV, even with the no improve in integrated luminosity, for the same R = 1 (see the leftmost curve in Fig. 2).

Forward-Backward Asymmetry
In the present study we limit ourselves to theoretical considerations only, so that we may ascertain how feasible asymmetry measurements are from a theoretical viewpoint. In fact one may argue that, due to large PDF uncertainties, a full scale simulation of the experiments is overkill before LHC running. It would be interesting to see a detailed simulation once LHC data becomes available, and new data about the PDFs is available. We restrict our attention to two main sources of error in the asymmetry extraction algorithm: the potential for misidentification of the quark direction and finite number of events passing the rapidity cut Y > Y c . We set aside experimental errors, coming for instance from the measurement of the energy of the final leptons.
Consider the angular distribution in Eq. (18) as a probability distribution function. This distribution has a mean defined as and a variance defined by Interestingly, A F B turns out to be just proportional to θ . It is possible to estimate the deviation of the mean θ from its true value as σ p / √ n − 1 where n is the number of events used to study the angular distribution. Hence we estimate the difference between the true and measured asymmetry due to statistics to be, Here n should be thought of as representing the number of events remaining after a cut in rapidity is applied to the total number of events N.
Secondly there is a probability of quark misidentification introduced by the methodology. The assumption that the initial quark direction be always the one of the out-going γ ′ or Z ′ , implies a possible misidentificationθ ↔ π −θ. Let the probability of misidentification be p(Y c ), then in a sample of n events np(Y c ) are misidentified. Let n F (n B ) events travel in the forward (backward) direction in the collider. Then n = n F + n B and the forward-backward asymmetry is defined as (n F − n B )/n in this notation. The worst case scenario would be that all np(Y c ) events were placed in n F but should have been placed in n B or vice versa. In this case the difference in true and measured asymmetry would be, Adding the misidentification and the statistical together we can define a 'total misidentification' of, We conservatively add the two linearly. Now we define P (Y c ) to be the fraction of events with |Y | > Y c , i.e. the probability an event survives the chosen rapidity cut. Then we can relate the number of events in the peak, N, to the number used for the asymmetry extraction, n, through P (Y c ) by, Suppose we are interested in measuring an asymmetry at some level of significance. To do this introduce E and require that ∆A F B < E . We must impose E > 2p(Y c ) so that the chosen cut value is capable of producing data of sufficiently low misidentification. This implies that we require, or equivalently, total events. For the model considered, at the partonic level the asymmetry is A F B ≃ −0.5 in the case of uū → µ + µ − , while for dd events the asymmetry would be A F B ≃ −0.15. (For comparison, SM exchange plus only one new resonance the asymmetry vanishes up to 3 decimal places.) The measured asymmetry will then be some average of the two. Using what we showed before, namely that the uū events give a larger contribution to the total rate, we estimate the combined asymmetry to be A F B ∼ −0.42. Therefore seeing an asymmetry at a collider would indicate that there are two particles in the new neutral resonance.
To measure such an effect at a collider one would require a 'total misidentification' E < 0.14 in order to conclude, at the 3σ level, an asymmetry is present. Consider the case that all events happen to fall precisely at the rapidity cut Y c . Then the probability p(Y c ) is calculable, yielding a conservative estimate: with real data having |Y | > Y c one can in principle use a larger value of p(Y c ). The fraction P (Y c ) can also be computed.
Optimal choice of Y c for a minimal N can be achieved by writing a simple minimization algorithm. A high Y c gives high quality data (i.e. contaminated by small misidentification systematics) with low p(Y c ) but few events will pass such the cut, reducing also P (Y c ). The converse will be true for a low cut Y c . We used an optimization algorithm finding a Y c such that the required number of events has been minimized.
We compute the minimum number of events required for a 3σ uncertainty at given mass. Using the approximate event rate given by eqs. (43) and (44)  modest gains in reach, compared to the gain obtained by going to larger energies, as shown in Fig. 4. One can also see from fig. 3 that in order to test the interesting range m > 2 TeV, one requires a very high integrated luminosity of 100 fb −1 . At a 10 TeV machine it is not possible to measure an asymmetry of a neutral heavy gauge boson in this range of most interest, and hence we omit showing the related plot.

PDF Uncertainties
So far the uncertainty on the PDFs themselves has been completely ignored. The Fortran module supplied by [14] has the capability to compute PDF uncertainties. As the computational cost to produce uncertainties with the package is very large, we use an approximate expression which captures the essential features of the PDF uncertainties. Namely, the quark PDFs are well known up to fractions of around 0.5, then they become poorly understood while the anti-quark PDFs are poorly understood across all of the fraction range. Two piecewise functions, one for the u and one for theū, approximate this behavior These uncertainties are approximations to the results given by using the package [14] for our chosen PDF set [15] and the results broadly agree with those published by the CTEQ collaboration for their PDF sets [16]. Even if an event is not identified as being uū or dd in our analysis, the uū dominates events and hence the uncertainty indicated will offer a good guide as to the situation. As the asymmetry extraction methodology involves cutting the data for large Y , it probes the region with a smaller anti-quark fraction and a larger quark fraction. If the quark fraction gets above 0.5 the uncertainty grows rapidly. At a 10 TeV collider this  Figure 4. The minimum ratio R between couplings of the new and SM boson to SM currents, such that at the 3σ level (see text for definition) a forward-backward asymmetry of -0.42 is detectable, plotted for various integrated luminosities at a collision energy of √ s =24 TeV and as a function of the heavy boson mass. The leftmost line is for an integrated luminosity of 1 fb −1 , the middle line for 10 fb −1 and the rightmost 100 fb −1 .
effect will come in faster as high fractions are required to make large invariant masses while the effect will be slower to appear in a 24 TeV collider.
In general, most events in a peak will have low rapidity. As a result we expect far lower uncertainties for event rate measurements compared to asymmetry measurements.

An Example Model
In order to illustrate the results of the model-independent study it is useful to apply them to a specific model. For example in the AdS/TC model in [17] a 5D model is formulated in AdS space and bounds on the (R, m) parameter space are discussed.
In this work, several qualitatively different regions of the (R, m) parameter space are identified. Firstly there is a region where the dominant decay mode of the techniρ particles is into SM gauge bosons, and not SM fermions. This region is marked as invisible in figs. 5 & 6 because the branching ratio into fermions would strongly suppress the Drell-Yan process of interest here. There is then a second region already excluded by precision measurements. Finally the constraint R 0.65 comes from specific details of the model.
Putting this information into a plot and overlaying the search capabilities of a collider produces Figs. 5 and 6. At the very high integrated luminosity of Fig. 5 most of the parameter space will be tested for the existence of new resonances. However, there is no reach in terms of asymmetry. Thus it will be hard to conclude whether there is more than one particle responsible for the peak. In contrast, with the higher collision energy of Fig. 6, a very large fraction of the parameter space presents a testable asymmetry.
Both Figs. 5 & 6 include a 1σ PDF uncertainty marked by thin red lines and a shading between. The lower of these two lines falls below the axis and is barely visible. It is clear that the event rate uncertainty is much smaller in the 24 TeV case with the cone barely visible on the plot. As for the asymmetry, the uncertainty at 14 TeV grows much faster for those masses where the 14 and 24 TeV lines are both able to reach. In the case of the 14 TeV curve, the uncertainty in coupling is 10 times the coupling estimate by 2 TeV. In comparison the coupling uncertainty will not reach a factor of 10 until 3.5 TeV in the 24 TeV machine.
For the representative point m = 2.6 TeV and R = 0.4 we simulate the asymmetry and see the rapidity cut in action. For simplicity, we assume that only uū events contribute at the partonic level. As explained, this is the dominant subprocess, although as a check simulations including dd were run which confirmed expectations or lowering the asymmetry mildly. In this limit, the asymmetry is A F B ≃ −0.5.
First let us look at the events with √ s = 24 TeV. We apply a cut rapidity Y > Y c = 1.0. In this way the misidentification due to the methodology is relatively small (≃ 0.1). Also, the number of events that pass the cut is large enough as to make the statistical element very small (∼ 0.03). As a result, the simulated set plotted in Fig 7 clearly shows a large asymmetry, compatible with the true one, and ∼ 4σ away from zero. By comparison, in Fig. 8 we show a set of simulated events with the same choices of parameters, but √ s = 14 TeV. In this case we apply a lower cut of Y > Y c = 0.75. As a result, the asymmetry of the simulated data is A F B = −0.45 ± 0.07(stat) ± 0.14(sys), lower than the expected true value, and affected by large uncertainties. A final comment is in order. In this numerical analysis we did not take into account the effect of experimental errors, nor any kinematical cuts that might be applied to extract analyzable data. For example, when | cosθ| ≃ 1 the experimental extraction of the kinematics of the two final leptons might be problematic. In order to understand how important this effect is, we extracted the asymmetry from the simulated data by imposing the further constraint | cosθ| < 0.8. This procedure effectively reduces the forward-backward asymmetry to A F B = −0.37 (A F B = −0.28) in the 24 (14) TeV case.
These two examples show how the method works in practice. However caution must be noted that in performing such simulations there is a degree of sensitivity to Drell-Yan production is suppressed or not accessible by the specific model. The (thick) blue line represents the minimum coupling ratio R required for at least 10 events to be observed. (Thin) blue lines follow the contours of an upper and lower 1σ PDF uncertainty . The region between the upper and lower 1σ uncertainties is shaded blue. The (thick) red line represents the minimum coupling ratio R for an asymmetry to be visible at the 3σ level (see text for definition). (Thin) red lines follow the contours of an upper and lower 1σ PDF uncertainty (the lower line largely falls below the axis). The region between the upper and lower 1σ uncertainties is shaded red.
the choice of PDF sets. We illustrate here that the method functions but very different results may be obtained depending on how the large-ξ distributions are modelled, and ultimately this problem will be resolved only once actual data are collected.

Non-degenerate Couplings
The region at low rapidity exhibits an interesting feature. Due to the differences in the up and down PDF distributions, the rapidity distributions of these two types of partonic events have very distinct shapes: up-type events peak at finite Y , while down-type peak at Y = 0. In principle, this might allow to discriminate between different values of the couplings R γ ′ and R Z ′ of the new resonances to the up and down quarks.
To investigate the feasibility of such a measurement we examine the two extreme cases, that R γ ′ = 0 versus R Z ′ = 0. We normalize the distributions to the same area, hence fixing the event number. The result is shown in Fig. 9. We see the largest effect is around Y = 0 where there is a 1 part in 25 difference, or 5%. The chosen situation is the Cos Θ Figure 7. An example simulation of the angular distribution for the predicted number of events at a √ s =24 TeV machine with L = 100 fb −1 , for boson mass m = 2.6 TeV and coupling ratio R = 0.4. The simulated data has been restricted to rapidity Y > Y c = 1.0 and the angular distribution in the parton-parton centre of mass frame prediction based on the detector angular distribution is plotted. The observed asymmetry is A F B = −0.51 ± 0.03(stat) ± 0.1(sys). This is the remnant of initial asymmetry of -0.42 that the data was simulated with. most extreme, a light mass in terms of what the precision constraints allow and a high energy collider. We can conclude that the effect of a non-degenerate R will be so small as not to be measurable, even assuming perfect detector efficiency and no kinematic cuts.

Conclusions
We performed a model-independent study of the Drell-Yan processes at the LHC for two new neutral gauge bosons with the quantum numbers of the photon and Z boson. We assumed the behavior of the new states are analogous to their standard-model cousins, but with couplings to SM fermions rescaled by a common factor of R < 1. We also chose to look at the case where the peaks of the two new states were so close as to be indistinguishable in the invariant mass plots.
We assumed perfect detector performance, omitting experimental errors, such as energy resolution, and without including the effect of kinematical cuts. We did not include any QCD corrections to the tree-level partonic calculation of pp → µ + µ − and pp → e + e − . In relation to the forward-backward asymmetry, we included in the analysis the main sources of theoretical uncertainty: a statistical element related to the number of expected events, a systematic due to the misidentification of the direction of the Cos Θ Figure 8. An example simulation of the angular distribution for the predicted number of events at a √ s =14 TeV machine with L = 100 fb −1 , for boson mass m = 2.6 TeV and coupling ratio R = 0.4. The simulated data has been restricted to rapidity Y > Y c = 0.75 and the angular distribution in the parton-parton centre of mass frame prediction based on the detector angular distribution is plotted. The observed asymmetry is A F B = −0.45 ± 0.07(stat) ± 0.14(sys). This is the remnant of initial asymmetry of -0.42 that the data was simulated with.  . The expected rapidity distribution of a sample of events in the case that a single new Z ′ (blue) and a single γ ′ (red) is the only exchange particle present in a new resonance. The new boson has a mass set at 2.6 TeV and a machine energy of √ s = 24 TeV is used. The same number of events is assumed for both.

Number of Events
initial (anti-)quark, and the uncertainty due to the PDF distributions.
The results can be interpreted as a best case scenario at the LHC, limited only by theoretical considerations. We compared different choices of the integrated luminosity and center-of-mass energy of the pp collisions. For illustration purposes, we also referred to one specific model, in which such mass-degeneracy of the new weakly-coupled gauge bosons is required by imposing the bounds on precision electro-weak parameters.
For a 10 TeV collider, and luminosities in the 2−40 fb −1 , bounds on new resonances no higher than 3 TeV can be obtained. With a 24 TeV collider the region well above 4 TeV can be explored, depending on the specific coupling R. The uncertainty imparted upon this prediction from the PDFs is quite modest for a high energy collider, while at lower energies, like 10 TeV, the PDF uncertainty at large momentum fraction is an important limiting factor. This is in substantial agreement with the literature.
We studied in detail the forward-backward asymmetry as a method of telling a two particle resonance from a one particle peak. The collider ability to measure an asymmetry at the 3σ level is heavily restricted to light masses. Even at very high luminosities measuring the asymmetry for masses of 3 TeV is not possible at a 14 TeV machine. Factoring in the PDF uncertainty, it is hard to conceive of any possible asymmetry measurements in a 14 TeV (or lower) collider, in relation to the scenarios considered here. By contrast, at a 24 TeV machine the PDF uncertainty is markedly lower and the reach far larger. At such a machine it may be possible to measure asymmetries and greatly restrict currently open parameter space.
In many phenomenological models one would wish to measure precisely the physics of heavy gauge bosons. Measuring the forward-backward asymmetry in events with neutral intermediate states provides a theoretically clean signal allowing one to discriminate the presence of more than one boson in a heavy resonance. Such a measurement is very difficult in practice from an experimental viewpoint. On top of the practicalities there are the large PDF uncertainties at large momentum fraction. Together with the scaling of the event rate with s, the situation suggests upgrading the LHC to a higher energy (besides increasing the luminosity) to learn about such new bosons, if the new bosons are heavy as is suggested by the precision tests.