Phenomenological Study of the Interplay between IRImproved DGLAPCS Theory and the Precision of an NLO ME Matched Parton Shower MC
S.K. Majhi ^{1}^{1}1Work supported by grant Pool No. 8545A, CSIR, IN.
Indian Association for the Cultivation of Science, Kolkata, India
Email: tpskm@iacs.res.in
A. Mukhopadhyay
Baylor University, Waco, TX, USA
Email: aditi_mukhopadhyay@baylor.edu
B.F.L. Ward
Baylor University, Waco, TX, USA
Email: bfl_ward@baylor.edu
S.A. Yost ^{2}^{2}2Work supported in part by U.S. D.o.E. grant DEFG0210ER41694 and grants from The Citadel Foundation.
The Citadel, Charleston, SC, USA
Email: scott.yost@citadel.edu
Abstract
We present a phenomenological study of the current status of the application of our approach of exact amplitudebased resummation in quantum field theory to precision QCD calculations, by realistic MC event generator methods, as needed for precision LHC physics. We discuss recent results as they relate to the interplay of the attendant IRImproved DGLAPCS theory of one of us and the precision of exact NLO matrixelement matched parton shower MC’s in the Herwig6.5 environment as determined by comparison to recent
LHC experimental observations on single heavy gauge boson production and decay. The level of agreement between the new theory and the data continues to be a reason for optimism. In the spirit of completeness, we discuss as well other approaches to the same theoretical predictions that we make here from the standpoint of physical precision with an eye toward the (sub)1% QCDEW
total theoretical precision regime for LHC physics.
BUHEPP1302, Apr., 2013
1 Introduction
With the recent announcement [1] of an BroutEnglertHiggs (BEH) [2] candidate boson after the startup and successful running of the LHC for three years, we have entered the era of precision QCD, by which we mean predictions for QCD processes at the total precision tag of or better. Given the expected role of precision comparison between theory and experiment in determining the detailed properties of the newly discovered EBH boson candidate, the attendant need for exact, amplitudebased resummation of large higher order effects is now more paramount. We have argued in Refs. [3, 4] that such resummation allows one to have better than 1% theoretical precision as a realistic goal in such comparisons, so that one can indeed distinguish new physics(NP) from higher order SM processes and can distinguish different models of new physics from one another as well. In what follows, with an eye toward its interplay with NLO exact, matrix element matched parton shower MC precision issues, we present the status of this approach to precision QCD for the LHC in connection with its attendant IRimproved DGLAPCS [5, 6] theory [7, 8] realization via HERWIRI1.031 [9] in the HERWIG6.5 [10] environment. We realize here the attendant exact, NLO matrix element matched parton shower MC realizations for both HERWIRI1.031 and HERWIG6.5 in the MC@NLO [11] methodology, wherewith we make the corresponding comparisons with recent LHC data that we present herein.
In this way, we continue the strategy of building on existing platforms to develop and realize a path toward precision QCD for the physics of the LHC. What we exhibit is an explicit union of the new IRimproved DGLAPCS theory and the MC@NLO realization of exact NLO matrix element(ME) matched parton shower MC theory. We note that we are also pursuing the implementation [12] of the new IRimproved DGLAPCS theory for HERWIG++ [13], HERWIRI++, for PYTHIA8 [14] and for SHERPA [15], as well as the corresponding NLO ME/parton shower matching realizations in the POWHEG [16] framework. For, our ultimate goal is a provable precision tag on our theoretical predictions and one of the strongest cross checks on theoretical precision is the difference between two independent realizations of the attendant theoretical calculation. We will present such cross checks elsewhere [12].
The paper is organized as follows. We motivate the development of the new IRimproved DGLAPCS theory in our approach to precision QCD theory by showing in the next section how it arises naturally in the effort to realize a provable precision from our approach [4] to precision LHC physics. This will allow us to expose properly the interplay between the NLO ME matched parton shower MC precision and the new IRimproved DGLAPCS theory. In the interest of completeness, we also review in the next section the relevant aspects of the approach in Ref. [4], which is an amplitudebased QEDQCD() exact resummation theory realized by MC methods. This paper is an extension of a shorter report of our results and analysis which is given in Ref. [17]. Accordingly, also in the interest of completeness, we discuss as well in the next section other possible approaches to the predictions we make herein from the standpoint of their physical precision. We believe this discussion will aid the reader in putting the results we present in a proper context. Section 3 contains the applications to the recent LHC data on single heavy gauge boson production with an eye toward the analysis in Refs. [9] of the analogous processes at the Tevatron. In this paper we will focus on the single Z/ production and decay to lepton pairs for definiteness. The other heavy gauge boson processes will be taken up elsewhere [12]. Section 4 contains our concluding observations and our outlook.
2 Review of Our Approach to Precision LHC Physics
We take the starting point for what we discuss here as the following fully differential representation of a hard LHC scattering process:
(1) 
where the and are the respective parton densities and reduced hard differential cross section and where we use the subscript on the latter to indicate that it has been resummed for all large EW and QCD higher order corrections in a manner consistent with achieving a total precision tag of 1% or better for the total theoretical precision of (1). The key theoretical issue for precision QCD for the LHC is then The proof of the correctness of the value of the total theoretical precision of (1) is then the key theoretical issue for the realization of precision QCD for the LHC. The theoretical precision in question can be represented as follows:
(2) 
where is the contribution of the uncertainty on the quantity to . Here, we stress that we discuss the situation in which the two errors in (2) are independent for definiteness; (2) has to be modified accordingly when they are not. We further stress that the proof of the correctness of the value of the total theoretical precision is indispensible in order to validate the application of a given theoretical prediction to precision experimental observations, for the discussion of the signals and the backgrounds for both Standard Model(SM) and new physics (NP) studies, and more specifically for the overall normalization of the cross sections in such studies. If one uses a calculation with an unknown value of for the attendant studies, the NP itself can be missed. We do feel a need to continue to state here that this point simply cannot be emphasized too much.
In the interest of completeness, we note that, by our definition, which follows the discussion in Ref. [18], is the total theoretical uncertainty that derives from the physical and technical precision contributions: the physical precision contribution, , arises from such sources as missing graphs, approximations to graphs, truncations,….; the technical precision contribution, , arises from such sources as bugs in codes, numerical rounding errors, convergence issues, etc. We want to observe that here, when we reference bugs in codes we have in mind that all gross errors such as those that give obviously wrong results, as determined by cross checks, are eliminated and we have left programming errors such as those in the logic: suppose for programming error reasons a DOloop ends at 999 steps instead of the intended 1000 steps, resulting in a per mille level error, that could alternate in sign from event to event. As per mille level accuracy is good enough in many applications, the program would remain reliable, but it would have what we call a technical precision error at the per mille level. With this understanding, the total theoretical error is then given by
(3) 
Although the desired value for depends on the specific requirements of the observations, as a general rule, it should fulfill , where is the respective experimental error and , so that the theoretical uncertainty does not significantly adversely affect the attendant physics analysis of the data.
In the interest of completeness, we will now recapitulate the theory we have developed to realize such precision rigorously. This theory is the resummation theory in Refs. [4] for the reduced cross section in (1) and for the resummation of the evolution of the parton densities therein as well. More precisely, the master formula, which applies both to the resummation of the reduced cross section and to that of the evolution of the parton densities, may be identified as
(4)  
(5)  
(6) 
where is either the reduced cross section or the differential rate associated to a DGLAPCS [5, 6] kernel involved in the evolution of the and where the new (YFSstyle [19, 20]) nonAbelian residuals have hard gluons and hard photons and we show the final state with two hard final partons with momenta specified for a generic final state for definiteness. The infrared functions are defined in Refs. [4, 7, 8] as follows:
(7)  
(8) 
where the dummy parameter is such that nothing depends on it and where we have introduced
(9)  
(10)  
(11) 
Here, the superscript denotes that the infrared functions are DGLAPCS synthesized as explained in Refs. [21, 4, 7, 8] and the infrared functions are given in Refs. [19, 20, 4, 7, 8]. We stress that simultaneous resummation of QED and QCD large IR effects is exact here.
We may describe the physical meanings of the key components in the master formula (6) as follows. In the language of Ref. [22], the exponent sums up to the infinite order the maximal leading IR singular terms in the cross section for soft emission below a dummy parameter and the exponent does the same for the regime above so that (6) is independent of – it cancels between and ^{3}^{3}3If we want to include more of the maximal exponentiating terms from the formalism of Ref. [22] in the two exponents , we may do so with a consequent change in the attendant residuals .. In order to maintain exactness order by order in perturbation theory in both and in the presence of these resummed terms, as explained in Refs. [4, 7, 8] we generate the residuals by iterative computation to match the attendant exact results to all orders in and . We need to emphasize the following most important point. In our formulation in (6) the entire soft gluon phase space is included in the representation – no part of it is dropped.
As it is explained in Refs. [4], the new nonAbelian residuals allow rigorous shower/ME matching via their shower subtracted analogs: in (6) we make the replacements
(12) 
where the have had all effects in the showers associated to the removed from them. To see how we make contact between the and the differential distributions in MC@NLO we proceed as follows. We represent the MC@NLO differential cross section via [11]
where is Born distribution, is the regularized virtual contribution, is the corresponding counterterm required at exact NLO, is the respective exact real emission distribution for exact NLO, is the parton shower real emission distribution so that the Sudakov form factor is
, where as usual it describes the respective noemission probability. The respective Born and real emission differential phase spaces are denoted by , respectively. We may note further that the representation of the differential distribution for MC@NLO in (2) is an explicit realization of the compensation between real and virtual divergent soft effects discussed in the Appendices of Refs. [7, 8] in establishing the validity of (6) for QCD – all of the terms on the RHS of (2) are infrared finite. Indeed, from comparison with (6) restricted to its QCD aspect we get the identifications, accurate to ,
(13) 
where we defined [11]
and we understand here that the DGLAPCS kernels in are to be taken as the IRimproved ones as we exhibit below [7, 8]. Here we have written the QCD virtual and real infrared functions and respectively without the superscript for simplicity of notation and they are understood to be DGLAPCS synthesized as explained in Refs. [4, 7, 8] so that we avoid doubling counting of effects. We also reemphasize that we do not drop any effects here in (13). We observe further that, in view of (13), the way to the extension of frameworks such as MC@NLO to exact higher orders in is therefore open via our and will be taken up elsewhere [12].
We point out that in Refs. [7, 8, 9] our methods for resummation of the QCD theory have been shown to be fully consistent with the methods in Refs. [23, 24]. Specifically, it is shown in Refs. [7, 8, 9] that the methods in Refs. [23, 24] give approximations to our hard gluon residuals ; for, unlike the master formula in (6), the methods in Refs. [23, 24] are not exact results. To see this, observe that the thresholdresummation methods in Refs. [23], using the result that, for any function ,
drop nonsingular contributions to the cross section at in resumming the logs in Mellin space. The SCET theory in Refs. [24] drops terms of at the level of the amplitude, where for a process with the hard scale with so that, for , we have . From the known equivalence of the two approaches, the errors in the threshold resummation must be similar. Evidently, we can only use these approaches as a guide to our new nonAbelian residuals as we develop results for the (sub)1% precision regime.
In view of the specific processes which we consider in the next Section, we will continue here with our consistency discussion as it relates to the theory of QCD resummation for the specific heavy gauge boson productiontype processes offered in Refs. [25, 26], which is again an approximate formalism that could be used to make approximations to our hard gluon residuals . The theory in Refs. [25, 26] is in wide use at the LHC and in the data analyses for the Tevatron – see for example the recent analyses in Refs. [27] where this theory, as it is implemented in the MC integration program RESBOS [28, 29, 30] is compared to recent LHC data and to recent analyses of Tevatron data. We need to note that the theory in Refs. [25, 26] builds on and extends beyond the considerable literature in Refs. [31] aimed at the analogous processes to those under study here. Could we perhaps employ this formalism in Refs. [25, 26] as well to reach physical precision QCD predictions? Let us note that the authors in Ref. [26] give us a hint to the answer to our question in their footnote on the journal page with number 215 for that paper, wherein they equate as ’negligible’ 20% in discussing possible nonperturbative contributions in their formalism. Let us keep this footnote in mind in what follows.
We first note that a defining formula for the approach in Refs. [25, 26] is that for the differential cross section for the distribution for the heavy gauge boson production in hadronhadron collisions, which we specialize here to the case of the DrellYan production for definiteness^{4}^{4}4The analogous results for the and are obtained by straightforward substitutional manipulations of the EW aspects of the formula we record here as described in Ref. [26] so that we omit such manipulations here without loss of content of the QCD aspects our discussion., where we record the result in Eq.(2.2) from Ref. [26]:
(14) 
where we have the usual kinematics so that is the transverse momentum, A,B are protons at the LHC, is the cms squared energy of the protons, is the 4momentum so that is the mass squared, and is the rapidity so that and . We have in mind that is near here. In (14), the term involving carries the effects from QCD resummation as developed in Refs. [25, 26] and the term includes those contributions which are ’regular’ at in the sense of Refs. [25, 26], i.e., order by order in perturbation theory they are derived from the parts of the attendant hard scattering coefficients that are less singular than or as . Our question concerns the physical precision of the term; for, the term is perturbative and can be computed in principle to the required accuracy by the standard methods.
The result for given in Refs. [25, 26] is as follows. When we have where is a typical QCD hadronic mass scale such as the inverse of the proton radius, the result for is [25, 26]
(15) 
(16) 
Here we show explicitly the dependences of the parton density functions (note that our corresponds to in Refs. [25, 26]), and we show as well the dependences of the perturbatively calculable exponentiation and scattering coefficient functions and as they are given in Refs. [25, 26]. The scale is usually set to and and are order 1 constants chosen to optimize the resultant perturbation expansions – see Refs. [25, 26] for more discussion on this point. The authors in Ref [25, 26] note that this result for , when , is accurate up to terms where represents the quark masses. In addition, the latter authors argue that in the result (14) they have dropped terms of in the regime where . These last two statements are seen explicitly in Sect. 9 of Ref. [25] in eq.(9.1) and the equation immediately preceeding it, for reference.
When we have to consider as well the regime , the authors in Ref. [25, 26] argue that we have to replace the functions and according to
(17) 
where we see that we can no longer neglect the quark masses and that we no longer have the convolutions : the functions reduce to the latter form when we take and drop terms , according to the arguments in Refs. [25, 26].
The authors in Refs. [25, 26] then join these two results, that for the regime and that for the regime , via an ansatz as follows. Defining
(18) 
which we see from (15) and (17) is welldefined for both regimes of interest in , the authors in Ref. [26] use
(19) 
where is a parameter still in the perturbative regime where (15) holds, to write the joining ansatz
(20) 
where the functions have nonperturbative content and must be determined from data. With this ansatz, the authors in Ref. [26] arrive at the representation
(21) 
We pointout for completeness and illustration that in RESBOS in Refs. [29, 30] three realizations of the nonperturbative functions in (21) were considered as follows:
(22) 
We observe that the are parameters on the RHS of (22) and that the best fit to the data considered in Ref. [30] was obtained from the BLNY parametrization with the values
(23) 
when and .
On the question of the physical precision of (21), we first observe that, with the errors shown on the constants for the best BLNY parmetrization, the error on represents already a 1.5% uncertainty at the saddle point position for the integration over in the respective production analogue of (14) found in Ref. [26] and recapitulated below in (24) – this alone exceeds the theory error budget in precision QCD theory that we advocate here so that one would need these errors reduced considerably for them to be useful for our purposes. But, what is more important is that, on the LHS of (22) the which do not (do) multiply are unspecified functions of that are required to vanish at while on the RHS these dependences are simplified to second order polynomials in and either no or linear dependence on : these simplifications are generically ad hoc and can not be considered as a rigorous platform for testing the fundamental QCD theory. They(the simplifications) do not follow from the formalism of Ref. [25, 26]. As the authors in Ref. [32] have emphasized, the latter formalism could only be used to give a prediction for perturbative QCD in the regime where the details of the nonperturbative parametrization are not important. We can not emphasize this point too much.
On the question of the physical precision of (21) we additionally note that in the regime it has an error of according to the authors in Ref. [25, 26], as such terms have been dropped, order by order in perturbation theory in the construction of – we make reference again specifically to Eq.(9.1) and the equation immediately preceding it in Ref. [25]. For and GeV, this is an error of and its twice this size at GeV. This shows that the formalism in Ref. [25, 26] can not be used for (sub)1% theory predictions for the heavy gauge boson spectrum at the LHC. One might think that as , where these dropped terms should vanish, the theory in Refs. [25, 26] would be accurate enough for (sub)1% precision QCD predictions. This is not true as well because at , following Parisi and Petronzio in Refs. [31], the authors in Ref. [26] show that the integral over in (14) is dominated by a saddle point at
(24) 
where we have and where is the number of effective quark flavors, which we take as here. Thus, we have [25, 26] so that
For , we get when we follow Ref. [26] and take . This result for is in the perturbative regime. Moreover, we know from Refs. [25, 26] that the error on the production analogue of at this point is – see again Eq.(9.1) in Ref. [25]. Thus, even when the perturbative regime obtains so that the result (15) should be reliable, the error in it is too large for use as anything but a guide in constructing our precision theory residuals in our master formula (6) for precision LHC physics. The results reported in Refs. [27] on the comparison of RESBOS with recent ATLAS and Tevatron data related to the spectrum in single are consistent with our estimates on the physical precision of (21): RESBOS misses the data by 2% for the near 0 and it misses the data by for the regime of .
The nonperturbative factor in (21) compromises the predictive power of the formalism in Refs. [25, 26]. The authors in Ref. [32] have argued that for W, Z production as considered here the prediction in (21) in the regime is insensitive to the nonperturbative parametrization in (22). This leaves the effects we have discussed above as the main obstacles to using the theory in Refs. [25, 26] for precision QCD physics predictions at the LHC in this latter regime when that precision is at or below 1%.
We note that in Ref. [33], another version of (14) is presented which has the same physical precision limits as that realized in RESBOS with a different treatment of the nonperturbative regime. As the comparison with the data shown in Refs. [27] shows in the perturbative regime, our estimates of this physical precision for this regime also apply to the realization of (14) in Ref. [33]: it cannot be used for 1% precision LHC physics studies.
Here, let us note that the authors in Ref. [33] discuss the uncertainty in their results as a function of the variation of the (perturbative scale) values associated with their renormalization, factorization and resummation scales. For example, they estimate that this is at the level of 10% near the peak of the spectrum in the single production at the LHC. These uncertainties can in principle be reduced by computing the perturbative terms in their results to higher and higher orders. They also estimate an additional PDF error at this region at the level of 2%, which in principle can be reduced by improving the determination of the respective PDF’s. We stress that the error in the defining result (14) derived in Refs. [25, 26] on which their results are based that we discuss here is an error of that applies order by order in the perturbation theory – it is separate from the scale and PDF errors discussed in Ref. [33].
Indeed, in Ref. [34], the SCET approach is used to recover (14) and it is shown explicitly that the improvement of the perturbative scale errors indeed occurs when higher order corrections are included in the calculation of the perturbative terms in the respective SCET realization of (14). We stress again that the error that we discuss here of due to the approximations in the defining derivation of (14) in Refs. [25, 26] also applies to the orderbyorder results in Refs. [34]. As we noted above, SCET involves for single production at the LHC the defining error for the typical hadronic transverse size and this is consistent with the approximations made in Refs. [25, 26]. Such approximations cannot be used for the 1% precision QCD theory that we have as our goal for (6) here.
In summary, this last remark is beginning to be manifest as we see in Refs. [27] in the comparison between the recent Tevatron and LHC data on the spectra in the production and the predictions of RESBOS and of Ref. [33]. Indeed, even though a new related variable [36] is used in some of the comparisons, , where is the azimuthal angle between the two leptons which have transverse momenta and is the scattering angle of the dilepton system relative to the beam direction when one boosts to the frame along the beam direction such that the leptons are back to back, one sees that these comparisons also show the underlying physical precision error associated with the defining formula in (14); as expected, the comparisons are somewhat better than the spectra comparisons because this quantity is more inclusive – two different values of with correspondingly compensating differences in the attendant can have the same value of . We want to encourage the LHC experimentalists to continue to produce the real spectra in the production because we believe, even if there are some irreducible experimental systematic errors, these spectra can be very useful in determining which theoretical approach is actually correct as we argue in the next section. In particular, given the intrinsic physical precision error in (14), one has to account for that error when the data are normalized to any prediction which is based on the formula therein. We would like to stress this latter point. We await the more complete data sets accordingly.
From the discussions just completed, we see that, in order to have a strict control on the theoretical precision in (1), we need both the resummation of the reduced cross section and that of the attendant evolution of the . We turn now to the latter.
More specifically, we apply the QCD restriction of the formula in (6) to the calculation of the kernels, , in the DGLAPCS theory itself and thereby get an improvement of the IR limit of these kernels, an IRimproved DGLAPCS theory [7, 8] in which large IR effects are resummed for the kernels themselves. The attendant new resummed kernels, are given in Refs. [7, 8, 9]. We reproduce the new kernels here for completeness:
(25) 
where the superscript “” indicates that the kernel has been resummed as predicted by Eq. (6) when it is restricted to QCD alone. Here () is the quadratic Casimir invariant for the quark(gluon) color representation respectively, and the YFS [19] infrared factor is given by
where is Euler’s gamma function and is Euler’s constant. The definitions of the respective resummation functions are as follows [7, 8] ^{5}^{5}5The improvement in Eq. (25) should be distinguished from the resummation in parton density evolution for the “” Regge regime – see for example Refs. [35, 37]. This latter improvement must also be taken into account for precision LHC predictions.:
(26) 
where is Euler’s gamma function and is Euler’s constant. We use a oneloop formula for in (26), so that
where is the number of active quark flavors. These new kernels provide us with a new resummed scheme for the parton density functions (PDF’s) and the reduced cross section with the same value of in (1):
(27) 
This new scheme has improved MC stability as discussed in Refs. [9] – in the attendant parton shower MC based on the new kernels there is no need for an IR cutoff ‘’ parameter. It is important to note that, while the degrees of freedom below the IR cutoffs in the usual showers are dropped in those showers, in the showers in HERWIRI1.031, as one can see from (6), these degrees of freedom are included in the calculation and are integrated over in the process of generating the GribovLipatov exponents in (25). We note also that the new kernels agree with the usual kernels at as the differences between them start in . This means that the NLO matching formulas in the MC@NLO and POWHEG frameworks apply directly to the new kernels for the realization of exact NLO ME/shower matching.
In the interest of pedagogy, in Fig. 1 we illustrate the basic physical idea, already discussed by Bloch and Nordsieck [38], which underlies the new kernels:
an accelerated charge generates a coherent state of very soft massless quanta of the respective gauge field so that one cannot know which of the infinity of possible states one has made in the splitting process illustrated in Fig. 1. In the new kernels this effect is taken into account by resumming the terms when is the IR limit. From (27) and (1), we see that when the usual kernels are used these terms are generated orderbyorder in the solution for the cross section in (1) so that our resumming them enhances the convergence of the representation in (1) for a given order of exactness in the input perturbative components therein. In the next Section, we illustrate this last remark in the context of the comparison of recent LHC data to NLO parton shower/matrix element matched predictions.
3 Interplay of IRImproved DGLAPCS Theory and NLO Shower/ME Precision: Comparison with LHC Data
In the new MC HERWIRI1.031 [9] we have the first realization of the new IRimproved kernels in the HERWIG6.5 [10] environment. Here, using recent LHC data as our baseline, we compare it with HERWIG6.510, both with and without the MC@NLO [11] exact correction to illustrate the interplay between the attendant precision in NLO ME matched parton shower MC’s and the new IRimprovement for the kernels.
More specifically, in Fig. 2 in panel (a) we show for the single production at the LHC the comparison between the CMS rapidity data [39] and the MC theory predictions and in panel (b) in the same figure we show the analogous comparison with the ATLAS data. Here, the rapidity data are the combined results and the data are those for the bare case; for, the theoretical framework of our simulations corresponds to these data. We do not as yet have complete realization of all the corrections involved in the other ATLAS data in Ref. [40].
These results are better appreciated if they are considered from the perspective of our analysis in Ref. [9] of the FNAL data on the single production in collisions at 1.96 TeV.
More precisely, we direct the reader to the results in Fig. 11 of the second paper in Ref. [9]. In that figure, we showed that, when the intrinsic rms parameter is set to 0 in HERWIG6.5, the MC@NLO/HERWIG6.510 simulations give a good fit to the CDF rapidity distribution data [42] therein but they do not give a satisfactory fit to the D0 distribution data [43] therein. In contrast, the corresponding simulations for MC@NLO/HERWIRI1.031 give good fits to both sets of data with the . Here corresponds to rms value for an intrinsic Gaussian distribution in . The authors of HERWIG [41] already have observed that, to get good fits to both sets of data, one may set GeV. Accordingly, in analyzing the new LHC data, we have set GeV in our HERWIG6.510 simulations while we continue to set PRTMS=0 in our HERWIRI simulations.
We turn now with this perspective to the results in Fig. 2, where we see a confirmation of the finding of the HERWIG authors. One needs to set [44] in the MC@NLO/HERWIG6510 simulations
to get a good fit to both the CMS rapidity data and the ATLAS data. We again see that the MC@NLO/HERWIRI1.031 simulations with at LHC give a good fit to the data for both the rapidity and the spectra.
Quantitatively, we use the as a measure of the goodness
of the respective fits. From the results in Fig. 2 we compute
that
the for the rapidity data and the for the
data are (.72,.72)((.70,1.37)) for the
MC@NLO/HERWIRI1.031(MC@NLO/HERWIG
6510(=2.2GeV)) simulations.
The corresponding results are (.70,2.23) for the
MC@NLO/HERWIG6510(=0) simulations.
To reproduce the LHC data on the distribution of the in the pp collision the usual DGLAPCS kernels require the introduction of a hard intrinsic Gaussian distribution in inside the proton whereas the IRimproved kernels give in fact a better fit to the data without the introduction of such a hard intrinsic component to the motion of the proton’s constituents. The hardness of this intrinsic is the issue as this quality of it is entirely ad hoc; it is in disagreement with the results of all successful models of the proton wavefunction [45], wherein the scale of the corresponding intrinsic is found to be GeV. More significantly, it contradicts the wellknown experimental observation of precocious Bjorken scaling [46, 47]; for, the famous SLACMIT experiments on the deep inelastic electronproton scattering process show that Bjorken scaling occurs already at GeV for with q the 4momentum transfer from the electron to the proton. If the proton constituents really had a Gaussian intrinsic distribution with GeV, these pioneering SLACMIT observations would not be possible. What we advocate now is that the ad hoc “hardness” of the GeV value is really just a phenomenological representation of the more fundamental dynamics described by the IRimproved DGLAPCS theory. This raises the following question: “Is possible to tell the difference between the two representations of the data in Fig. 2?”
One expects physically that more detailed observations should be able to distinguish the two representations of the data in
Fig. 2. In this connection, in Fig. 3 we show for the mass spectrum the MC@NLO/HERWIRI1.031(blue squares) and MC@NLO/HER
WIG6510(=2.2GeV) (green squares) predictions when the decay lepton pairs
satisfy the LHC type requirement that their transverse momenta exceed GeV.
From the results in Fig. 3, wherein the peaks differ by 2.2% for example, we see that the high precision data such as the LHC ATLAS and CMS experiments will have (each already has over lepton pairs) would allow one to distinguish between the two sets of theoretical predictions.
Continuing in this direction with the discussion, we see that the main differences between the three predictions in Fig. 2 (b) occur in the regime below GeV/c. To better probe this latter regime, we make a more detailed snapshot of it in which we plot in Fig. 4 the three respective featured theory predictions with the finer binning of GeV/c instead of the GeV/c binning used in Fig. 2 (b).
The results in Fig. 4 show that the three theoretical predictions have significant differences in the shapes that are testable with the precise data that will be available to the ATLAS and CMS experiments. We would again encourage experimentalists to pursue the measurements of both and spectra as these will both be very useful in establishing the correct theoretical approach to the respective LHC observations. We will pursue elsewhere [12] other such detailed observations that may also reveal the differences between the two descriptions of parton shower physics. In this connection especially, we continue to await the release of the entire data sets from ATLAS and CMS.
4 Conclusions
What we have shown is the following. The realization of IRimproved DGLAPCS theory in HERWIRI1.031, when used in the MC@NLO/HERWIRI1.031 exact ME matched parton shower framework, affords one the opportunity to explain, on an eventbyevent basis, both the rapidity and the spectra of the in pp collisions in the recent LHC data from CMS and ATLAS, respectively, without the need of an unexpectedly hard intrinsic Gaussian distribution with rms value of GeV in the proton’s wave function. Our view is that this can be interpreted as providing a rigorous basis for the phenomenological correctness of such unexpectedly hard distributions insofar as describing these data using the usual unimproved DGLAPCS showers is concerned. Accordingly, we have proposed that comparison of other distributions such as the invariant mass distribution with the appropriate cuts and the more detailed spectra in the regime below GeV be used to differentiate between these phenomenological representations of parton shower physics in MC@NLO/HERWIG6510 and the fundamental description of the parton shower physics in MC@NLO/HERWIRI1.031. We have further emphasized that the precociousness of Bjorken scaling argues against the fundamental correctness of the hard scale intrinsic ansatz with the unexpectedly hard value of GeV, as do the successful models [45] of the proton’s wave function, which would predict this value to be GeV. As an added bonus, we have pointedout that the fundamental description in MC@NLO/HERWIRI1.031 can be systematically improved to the NNLO parton shower/ME matched level – a level which we anticipate is a key ingredient in achieving the (sub)1% precision tag for such processes as single heavy gauge boson production at the LHC.
Evidently, relative to what one could achieve from the fundamental representation of the corresponding physics via IRimproved DGLAPCS theory as it is realized in HERWIRI1.031 when employed in MC@NLO/HERWIRI1.031 simulations, the use of ad hoc hard scales in models would compromise any discussion of the attendant theoretical precision. We are pursuing additional cross checks of the MC@NLO/HERWIRI1.031 simulations against the LHC data.
In our discussion, we have also spent some amount of time discussing alternative approaches to the type of resummation embodied in the IRimproved DGLAPCS theory; for, some of these approaches are in wide use. What we conclude is that the physical precisions of these approaches(see Sect. 2), which are based on Refs. [23, 24, 25, 26], are above the 1% precision tag that we aspire, even though there is no contradiction between our exact approach and these more approximate methods.
In closing, two of us (A.M. and B.F.L.W.) thank Prof. Ignatios Antoniadis for the support and kind hospitality of the CERN TH Unit while part of this work was completed.
References
 [1] F. Gianotti, in Proc. ICHEP2012, in press; J. Incandela, ibid., 2012, in press; G. Aad et al., Phys. Lett. B716 (2012) 1, arXiv:1207.7214; D. Abbaneo et al., ibid.716 (2012) 30, arXiv:1207.7235.
 [2] F. Englert and R. Brout, Phys. Rev. Lett. 13 (1964) 312; P.W. Higgs, Phys. Lett. 12 (1964) 132; Phys. Rev. Lett. 13 (1964) 508; G.S. Guralnik, C.R. Hagen and T.W.B. Kibble, ibid. 13 (1964) 585.
 [3] B.F.L. Ward, S.K. Majhi and S.A. Yost, in PoS(RADCOR2011) (2012) 022.
 [4] C. Glosser, S. Jadach, B.F.L. Ward and S.A. Yost, Mod. Phys. Lett. A 19(2004) 2113; B.F.L. Ward, C. Glosser, S. Jadach and S.A. Yost, in Proc. DPF 2004, Int. J. Mod. Phys. A 20 (2005) 3735; in Proc. ICHEP04, vol. 1, eds. H. Chen et al.,(World. Sci. Publ. Co., Singapore, 2005) p. 588; B.F.L. Ward and S. Yost, preprint BUHEPP0505, in Proc. HERALHC Workshop, CERN2005014; in Moscow 2006, ICHEP, vol. 1, p. 505; Acta Phys. Polon. B 38 (2007) 2395; arXiv:0802.0724, PoS(RADCOR2007)(2007) 038; B.F.L. Ward et al., arXiv:0810.0723, in Proc. ICHEP08; arXiv:0808.3133, in Proc. 2008 HERALHC Workshop,DESYPROC200902, eds. H. Jung and A. De Roeck, (DESY, Hamburg, 2009)pp. 180186, and references therein.
 [5] G. Altarelli and G. Parisi, Nucl. Phys. B126 (1977) 298; Yu. L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641; L. N. Lipatov, Yad. Fiz. 20 (1974) 181; V. Gribov and L. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 675, 938; see also J.C. Collins and J. Qiu, Phys. Rev. D39 (1989) 1398.
 [6] C.G. Callan, Jr., Phys. Rev. D2 (1970) 1541; K. Symanzik, Commun. Math. Phys. 18 (1970) 227, and in Springer Tracts in Modern Physics, 57, ed. G. Hoehler (Springer, Berlin, 1971) p. 222; see also S. Weinberg, Phys. Rev. D8 (1973) 3497.
 [7] B.F.L. Ward, Adv. High Energy Phys. 2008 (2008) 682312.
 [8] B.F.L. Ward, Ann. Phys. 323 (2008) 2147.
 [9] S. Joseph et al., Phys. Lett. B685 (2010) 283; Phys. Rev. D81 (2010) 076008.
 [10] G. Corcella et al., hepph/0210213; J. High Energy Phys. 0101 (2001) 010; G. Marchesini et al., Comput. Phys. Commun.67 (1992) 465.
 [11] S. Frixione and B.Webber, J. High Energy Phys. 0206 (2002) 029; S. Frixione et al., arXiv:1010.0568; B. Webber, talk at CERN, 03/30/2011; S. Frixione, talk at CERN, 05/04/2011.
 [12] A. Mukhopadhyay et al., to appear.
 [13] M. Bahr et al., arXiv:0812.0529 and references therein.
 [14] T. Sjostrand, S. Mrenna and P. Z. Skands, Comput. Phys. Commun. 178 (2008) 852867.
 [15] T. Gleisberg et al., J.High Energy Phys. 0902 (2009) 007.
 [16] P. Nason, J. High Energy Phys. 0411 (2004) 040.
 [17] S. Majhi et al., Phys. Lett. B 719 (2013) 367.
 [18] See for example S. Jadach et al., in Physics at LEP2, vol. 2, (CERN, Geneva, 1995) pp. 229298.
 [19] D. R. Yennie, S. C. Frautschi, and H. Suura, Ann. Phys. 13 (1961) 379; see also K. T. Mahanthappa, Phys. Rev. 126 (1962) 329, for a related analysis.
 [20] See also S. Jadach and B.F.L. Ward, Comput. Phys. Commun. 56(1990) 351; Phys.Lett. B 274 (1992) 470; S. Jadach et al., Comput. Phys. Commun. 102 (1997) 229; S. Jadach, W. Placzek and B.F.L Ward, Phys. Lett. B 390 (1997) 298; S. Jadach, M. Skrzypek and B.F.L. Ward,Phys. Rev. D 55 (1997) 1206; S. Jadach, W. Placzek and B.F.L. Ward, Phys. Rev. D 56 (1997) 6939; S. Jadach, B.F.L. Ward and Z. Was,Phys. Rev. D 63 (2001) 113009; Comp. Phys. Commun. 130 (2000) 260; ibid.124 (2000) 233; ibid.79 (1994) 503; ibid.66 (1991) 276; S. Jadach et al., ibid.140 (2001) 432, 475.
 [21] B.F.L. Ward and S. Jadach, Mod. Phys. Lett. A14 (1999) 491.
 [22] J.G.M. Gatheral, Phys. Lett. B 133 (1983) 90.
 [23] G. Sterman, Nucl. Phys. B 281, 310 (1987); S. Catani and L. Trentadue, Nucl. Phys. B 327, 323 (1989); ibid. 353, 183 (1991).
 [24] See for example C. W. Bauer, A.V. Manohar and M.B. Wise, Phys. Rev. Lett. 91 (2003) 122001; Phys. Rev. D 70 (2004) 034014; C. Lee and G. Sterman, Phys. Rev. D 75 (2007) 014022.
 [25] J.C. Collins and D.E. Soper, Nucl. Phys. B193 (1981) 381; ibid.213(1983) 545; ibid. 197(1982) 446.
 [26] J.C. Collins, D.E. Soper and G. Sterman, Nucl. Phys. B250(1985) 199; in Les Arcs 1985, Proceedings, QCD and Beyond, pp. 133136.
 [27] S. Hassani, in Proc. Recntres de Moriond EW, 2013, in press; H. Yin, ibid., 2013, in press; G. Aad et al., arXiv:1211.6899, and references therein.
 [28] C. Balazs and C. Yuan, Phys. Rev. D56 (1997) 5558.
 [29] G. Ladinsky and C. Yuan, Phys. Rev. D50 (1994) 4239.
 [30] F. Landry et al., Phys. Rev. D63 (2003) 073016.
 [31] Yu.L. Dokshitser, D.I. D’yakonov and S.I. Troyan, Phys. Lett. B78 (1978) 290; Proc. 13th Winter School of the LPNI, Leneingrad(1978); Phys. Rept. 58 (1980) 271; G. Parisi and R. Petronzio, Nucl. Phys. B 154 (1979) 427; G. Curci, M. Greco and R. Srivastava, Phys. Rev. Lett. 43 (1979) 834; Nucl. Phys. B 159 (1979) 451; P. Chiappetta and M. Greco, Nucl.Phys. B199 (1982) 77; ibid. 221 (1983) 268; Phys. Lett. B135 (1984) 187; F. Halzen, A.D. Martin and D.M. Scott, Phys. Rev. D25 (1982) 754; P. Aurenche and R. Kinnunen, Phys. Lett. B 135 (1984) 493; A. Nakamura, G. Pancheri and Y. Srivastava, Z. Phys. 21 (1984) 243; G. Altarelli, R.K. Ellis, M. Greco and G. Martinelli, Nucl. Phys. B 246 (1984) 12.
 [32] C. Davies, B. Webber and W. Stirling, Nucl. Phys. B256 (1985) 413.
 [33] A. Banfi et al., Phys. Lett. B715(2012) 152 and references therein.
 [34] T. Becher, M. Neubert and D. Wilhelm, arXiv:1109.6027; T. Becher, M. Neubert, Eur. Phys. J. C71 (2011) 1665, and references therein.
 [35] B.I. Ermolaev, M. Greco and S.I. Troyan, PoS DIFF2006 (2006) 036, and references therein.
 [36] A. Banfi et al., Eur. Phys. J. C71 (2011) 1600.
 [37] G. Altarelli, R.D. Ball and S. Forte, PoS RADCOR2007 (2007) 028.
 [38] F. Bloch and A. Nordsieck, Phys. Rev. 52 (1937) 54.
 [39] S. Chatrchyan et al., arXiv:1110.4973; Phys. Rev. D85 (2012) 032002.
 [40] G. Aad et al., arXiv:1107.2381; Phys. Lett. B705 (2011) 415.
 [41] M. Seymour, “Event Generator Physics for the LHC”, CERN Seminar, 2011.

[42]
C. Galea, in Proc. DIS 2008, London, 2008,
http://dx.doi.org/10.3360/dis.2008.55
.  [43] V.M. Abasov et al., Phys. Rev. Lett. 100, 102002 (2008).
 [44] P. Skands, private communication, 2011, finds a similar behavior in PYTHIA8 simulations.
 [45] R.P. Feynman, M. Kislinger and F. Ravndal, Phys. Rev. D3 (1971) 2706; R. Lipes, ibid.5 (1972) 2849; F.K. Diakonas, N.K. Kaplis and X.N. Mawita, ibid. 78 (2008) 054023; K. Johnson, Proc. Scottish Summer School Phys. 17 (1976) p. 245; A. Chodos et al., Phys. Rev. D9 (1974) 3471; ibid. 10 (1974) 2599; T. DeGrand et al., ibid. 12 (1975) 2060.
 [46] See for example R.E. Taylor, Phil. Trans. Roc. Soc. Lond. A359 (2001) 225, and references therein.
 [47] J. Bjorken, in Proc. 3rd International Symposium on the History of Particle Physics: The Rise of the Standard Model, Stanford, CA, 1992, eds. L. Hoddeson et al. (Cambridge Univ. Press, Cambridge, 1997) p. 589, and references therein.