Associated production of a top pair and a Z boson at the LHC to NNLL accuracy

Journal of High Energy Physics, Apr 2017

We study the resummation of soft gluon emission corrections to the production of a top-antitop pair in association with a Z boson at the Large Hadron Collider to next-to-next-to-leading logarithmic accuracy. By means of an in-house parton level Monte Carlo code we evaluate the resummation formula for the total cross section and several differential distributions at a center-of-mass energy of 13 TeV, and we match these calculations to next-to-leading order results.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1007%2FJHEP04%282017%29105.pdf

Associated production of a top pair and a Z boson at the LHC to NNLL accuracy

Received: February Associated production of a top pair and a Z boson at the LHC to NNLL accuracy Alessandro Broggio 0 1 3 6 7 Andrea Ferroglia 0 1 3 4 5 7 Giovanni Ossola 0 1 3 4 5 7 Ben D. Pecjak 0 1 2 3 7 Ray D. Sameshima 0 1 3 4 5 7 Fifth Avenue 0 1 3 7 New York 0 1 3 7 NY 0 1 3 7 U.S.A. 0 1 3 7 Open Access 0 1 3 7 c The Authors. 0 1 3 7 0 Department of Physics, University of Durham , Science Laboratories 1 300 Jay Street, Brooklyn, NY 11201 U.S.A 2 Institute for Particle Physics Phenomenology, Ogden Centre for Fundamental Physics 3 James Franck-Stra e 1 , D-85748 Garching , Germany 4 Physics Department, New York City College of Technology, The City University of New York 5 The Graduate School and University Center, The City University of New York 6 Physik Department T31, Technische Universitat Munchen 7 South Rd, Durham DH1 3LE, U.K We study the resummation of soft gluon emission corrections to the production of a top-antitop pair in association with a Z boson at the Large Hadron Collider to next-tonext-to-leading logarithmic accuracy. By means of an in-house parton level Monte Carlo code we evaluate the resummation formula for the total cross section and several di erential distributions at a center-of-mass energy of 13 TeV, and we match these calculations to nextto-leading order results. NLO Computations; QCD Phenomenology 1 Introduction Outline of the calculation Numerical results Scale choices Total cross section Di erential distributions Introduction The associated production of a top pair and a Z or W boson are the two processes with the heaviest nal states measured to date at the Large Hadron Collider (LHC). The total cross section for these processes was measured during Run I [1, 2], and preliminary measurements at a center-of-mass energy of 13 TeV are also available [3, 4]. The ttZ production process is particularly interesting because it allows one to study the coupling of the Z boson with the top quark. This measurement further tests the Standard Model (SM) of particle physics and probes several Beyond the SM scenarios that predict changes to this coupling with respect to the SM. In addition, these production processes lead to high multiplicity states which are background in the search for new heavy states decaying via long chains, such as dark matter candidates. Given their importance for phenomenological studies, next-to-leading-order (NLO) QCD and electroweak corrections to the associated production of a top pair and a massive vector boson were studied by several groups [5{14]. A full calculation of the QCD corrections to next-to-next-to leading order (NNLO) accuracy would be desirable but it is extremely di cult even with the most up to date techniques for the calculations of higher order corrections. However, the associated production of a top pair and a heavy colorless boson is a multiscale process which is expected to receive potentially large corrections arising from soft gluon emission. The resummation of these e ects to next-to-next-to leading logarithmic (NNLL) accuracy can be carried out by exploiting the factorization properties of the partonic cross section in the soft limit (which can be studied with e ective theory methods1) and by subsequently employing renormalization group improved perturbation theory techniques. In the case of the associated production of a top pair and a Higgs boson the resummation formula in the soft emission limit was discussed in [16], and results for the total cross section and several di erential distributions at NLO+NNLL accuracy 1For an introduction see [15]. were presented in [17]. Studies of the associated production of a top quark pair and a W boson to NLO+NNLL accuracy can be found in [18], where the resummation was carried out in Mellin moment space as in [17], and in [19], where the resummation was instead carried out in momentum space. The results of [17] and [18] were obtained by means of an in-house parton level Monte Carlo code for the numerical evaluation of the resummation formula. The output of this code was then matched to complete NLO calculations obtained by employing MadGraph5_aMC@NLO [20] (which we indicate with MG5 aMC in the following). Building on the results of those two papers, in this work we obtain a resummation formula for the associated production of ttZ nal state, and we evaluate it to NNLL accuracy by means of dedicated parton level Monte Carlo code. We match our results for the total cross section and di erential distributions to NLO calculations in order to obtain predictions at NLO+NNLL accuracy. The paper is organized as follows: in section 2 we introduce some basic notation and we brie y summarize the main steps in our calculations. For a more technical discussion of the methods employed in this paper, we refer the reader to the detailed descriptions provided in [16{18]. In section 3 we present predictions at NLO+NNLL accuracy for the total cross section as well as for several di erential distributions. Finally, we draw our conclusions in section 4. Outline of the calculation The associated production of a top quark pair and a Z boson receives contributions from the partonic process i(p1) + j(p2) ! t(p3) + t(p4) + Z(p5) + X ; where ij 2 fqq; qq; ggg at lowest order in QCD, and X indicates the unobserved partonic nal-state radiation. The two Mandelstam invariants which are relevant for our discuss^ = (p1 + p2)2 = 2p1 p2 ; M 2 = (p3 + p4 + p5)2 : The soft or partonic threshold limit is de ned as the kinematic region in which z 1. In this region, the nal state radiation indicated by X in (2.1) can only be soft. The factorization formula for the QCD cross section in the partonic threshold limit is the same as the one derived in [16] for ttH production, up to the straightforward replacement of the Higgs boson with a Z boson: M 2=s^ ! (s; mt; mZ ) = 1 Z 1 dPSttZ Tr Hij (fpg; ) Sij We indicated with s the square of the hadronic center-of-mass energy and we de ned min = (2mt + mZ )2 =s and hard functions H, soft functions S, and luminosity functions ff, as well as for the nal-state phase-space integration measure, is the same one used in [17, 18] and we refer the reader to these papers for more details. Similarly to LO, the only subprocesses to be considered in the soft limit are those labeled by indices ij 2 fqq; qq; ggg. The hard and soft functions are twoby-two matrices in color space for qq-initiated (quark annihilation) processes, and three-bythree matrices in color space for gg-initiated (gluon fusion) processes. Contributions from other production channels such as qg and qg (collectively referred to as \quark-gluon" or simply \qg" channel in what follows) are subleading in the soft limit. The hard functions satisfy renormalization group equations governed by the channel dependent soft anomalous dimension matrices ij . These anomalous dimension matrices were derived in [21, 22]. In order to carry out the resummation to NNLL accuracy, the hard functions, soft functions, and soft anomalous dimensions must be computed in xed-order perturbation theory up to NLO in s. The NLO soft functions and soft anomalous dimensions are the same ones needed in the calculation of ttH and ttW to NNLL accuracy and can be found in [16{18]. The NLO hard functions are instead process dependent, receive contributions exclusively from the virtual corrections to the tree level amplitudes, and were evaluated by customizing the one-loop provider Openloops [23], which we used in combination with the tensor reduction library Collier [24{27]. The NLO hard function have been cross-checked numerically by means of a customized version of GoSam [28{31], used in combination with the reduction provided by Ninja [32{34]. In this paper we carry out the resummation in Mellin space, starting from the relation (s; mt; mZ ) = 1 Z c+i1 where fefij is the Mellin transform of the luminosity functions, and c is the Mellin transform of the product of the hard and soft function (see [17, 18] for details). Since the soft limit z ! 1 corresponds to the limit N ! 1 in Mellin space, we neglected terms suppressed by The hard and soft functions included in the hard scattering kernels ec in (2.4) can be evaluated in xed order perturbation theory at scales at which they are free from large logarithms. We indicate these scales with h and s, respectively. Subsequently, by solving the renormalization group (RG) equations for the hard and soft functions one can evolve the factor ec to the factorization scale f . Following this procedure one nds s are resummed in the channel-dependent matrix-valued evolution factors Ue . The expression for the evolution factors is formally identical to the one found for ttW and ttH production and can be found for example in equation (3.7) of [18]. s (MZ ) from MMHT 2014 PDFs The l.h.s. of (2.5) is formally independent of h and s. In practice however, one cannot evaluate the hard and soft functions at all orders in perturbation theory; this fact creates a residual dependence on the choice of the scales h and s in any numerical evaluation of ec. The hard and soft functions are free from large logarithms if one chooses h values of N in the hard scattering kernels ec, whose existence is related to the Landau pole . We deal with this issue by choosing the integration path in the complex N plane according to the Minimal Prescription (MP) introduced in [35]. Finally, the parton luminosity functions in Mellin space, which we need in the numerical evaluations, are constructed using techniques described in [36, 37]. Numerical results The main purpose of this section is to present predictions for the associated production of a top pair and a Z boson to NLO+NNLL accuracy. However, we also analyze systematically the relevance of soft emission corrections and their resummation in relation to NLO predictions for the various observables considered in the paper. The NNLL calculations are carried out by means of an in-house parton level Monte Carlo code, while the NLO predictions are obtained by means of MG5 aMC. All of the calculations discussed in this section are carried out with the input parameters listed in table 1. Throughout the paper we employ MMHT 2014 PDFs [38]. In xed order calculations, the order of the PDFs matches the perturbative order of the calculation (i.e. LO calculations are carried out with LO PDFs, NLO calculations employ NLO PDFs, etc.). In matched calculations, we employed NLO PDFs for NLO+NLL accuracy, and NNLO PDFs for NLO+NNLL accuracy. For both the total cross section and several di erential distributions, we consider six di erent types of predictions: i) NLO calculations, obtained with MG5 aMC. ii) Approximate NLO calculations, obtained from the NLO expansion of the NNLL resummation formula. We check that for our choice of scales and input parameters approximate NLO calculations provide a satisfactory approximation to the exact NLO calculation. The approximate NLO formulas obtained by expanding (2.5) account for the single and double powers of ln N as well as N -independent terms but not for terms suppressed by inverse powers of N . N -independent terms depend on the Mandelstam variables, however we refer to them as \constant" terms in what follows. version of (2.5). Approximate NLO calculations are carried out with the in-house parton level Monte Carlo code which was developed speci cally for this project. iii) NLO+NLL calculations, which are obtained by matching NLO results with resummed results at NLL accuracy obtained by means of the in-house Monte Carlo code. The results are matched according to the formula NLO+NLL = NLO + The terms in the square brackets, which contribute to NLO and beyond, depend on the scales s and h in addition to the factorization scale f . Of course the dependence on h is formally of NNLL order; by varying these scales and the factorization scale in (3.1) one can estimate the size of the NNLL corrections. Because of the matching procedure described by equation (3.1), all of the NLO+NLL calculations in this paper include the contribution of the qg channel to NLO. iv) NLO+NNLL predictions are obtained by evaluating the hard scattering kernels in (2.5) to NNLL accuracy with the in-house Monte Carlo code and by matching the results to NLO calculations as follows: NLO+NNLL = NLO + The terms in the squared brackets in (3.2) contribute starting from NNLO and represent the NNLL corrections to be added to the NLO result. We stress that the NLO+NNLL calculations in this paper include the contribution of the qg channel to NLO through the matching procedure in (3.2). v) Approximate NNLO calculations are obtained by the NNLL resummation formula and include all powers of ln N and part of the constant terms from a complete NNLO calculation. The approximate NNLO formulas employed in this paper are constructed as the ones employed in [17, 18] for ttW and ttH production. A detailed description of the constant terms which are included in the approximate NNLO formulas can be found in section 4 of [16]. Approximate NNLO formulas are evaluated with the in-house Monte Carlo code which we developed and they are matched to the NLO calculations as follows nNLO = approx. NNLO where we label the matched result \nNLO" for brevity. By construction nNLO predictions are independent from the hard and soft scales but they do have a residual N3LO dependence on f . vi) NLO+NNLL expanded to NNLO. Finally we consider a second way of expanding the NNLL resummation formula to NNLO. This approach di ers from the approximate NNLO result used above by constant terms, which are formally of N3LL accuracy. This approximation is de ned by the relation NNLO exp. = h NNLL expanded to NNLO while the NLO+NNLL and nNLO curves are obtained using MMHT 2014 NNLO PDFs. The constant pieces in (3.4) contain explicit dependence on h and s, in addition to that on f . This dependence is formally an e ect of N3LL order. By comparing the predictions obtained from (3.4) to the corresponding NLO+NNLL calculations one can see the relative weight of terms of N3LO and higher in the NLO+NNLL calculations. If in the future a complete NNLO calculation for the ttZ production cross section were to become available, it would be possible to match it to the NNLL resummation formula by using precisely this kind of NNLO expansion of the NNLL resummation, as can be seen by replacing N ! NN in all of the superscripts in (3.1). Scale choices Since any numerical evaluation of the resummed expression for the hard scattering kernels must be carried out by evaluating the factors in (2.5) up to a certain order in perturbation theory, the resummed kernels ec will show a residual dependence on the scales s and In order to follow closely the approach adopted in \direct QCD" calculations [35, 39, 40], and s;0 = M=N [17, 18, 41, 42]. In addition, both xed order and resummed calculations depend on the factorization are not large [43]. It is therefore reasonable to choose a dynamical default value for the scale f which is related to M . The dependence of the total ttZ production cross section on f =M at the LHC with p to a di erent perturbative approximation, as indicated in the legend. Figure 1 shows that di er signi cantly for f =M 0:5 and for f =M 0:5. Following this observation, the The uncertainty related to the choice of the factorization scale in xed order results is depend also on the hard and soft scales, consequently, the uncertainty of the resummed results is estimated by varying separately all the three scales around their default values central value of a resummed observable O, which can be the total cross section or the value of the di erential cross section in a given bin, is determined as follows. First one evaluates the quantities in quadrature Oi+ ( Oi ) for i 2 fs; f; hg. Total cross section In this section we analyze the total cross section for the associated production of a top quark pair and a Z boson at the LHC operating at a center-of-mass energy of 13 TeV. The relevant results are collected in table 2. rst compare the approximate NLO cross section, obtained by expanding the resummation formula to NLO (second row of table 2) with the complete NLO cross section (fourth row) and the NLO cross section without the contribution of the quark-gluon channel (third row). The di erence between the approximate NLO result and the NLO result without the qg channel is due to terms in the quark annihilation and gluon fusion channels which are subleading in the partonic threshold limit. We see that the impact of these terms is around 1%. The di erence between these two results is therefore small in spite of the fact that the NLO corrections are large, as can be seen by comparing them with the LO result. However, we see that the approximate NLO result shows a smaller scale uncertainty than the NLO result with the contribution of the qg channel. We conclude that the soft emission corrections provide the bulk of the NLO corrections for this choice of the factorization scale. This motivates us to study the e ect of the resummation of these corrections, keeping in mind that by matching the resummed results to NLO calculations we consider both power corrections and the contribution of the qg channel to that order. The NLO+NLL and NLO+NNLL cross sections, shown in the sixth and seventh line of table 2 are main results of this paper. By looking at the NLO, NLO+NLL, NLO+NNLL results we see that the cross section is progressively increased, but the central value of each prediction falls in the scale uncertainty band of the predictions of lower accuracy. One might also want to evaluate the NLO+NNLL total cross section by employing NLO PDFs; this leads to a total cross section of 10 fb larger than the NLO+NNLL calculation carried out with NNLO PDFs and a scale uncertainty interval which is almost identical to the one obtained by using NNLO PDFs. A comparison of the of the NLO+NNLL total cross section evaluated with NLO PDFs with the NLO total cross section in table 2 allows one to assess directly the numerical impact of the higher-order logarithms which are included in the resummed calculation. One can then look at the NNLO expansions of the NNLL resummation formula, which are shown in the last two lines of table 2. By comparing these results to the NLO+NNLL app. NLO NLO no qg (NLO+NNLL)NNLO exp: in-house MC in-house MC +MG5 aMC in-house MC +MG5 aMC in-house MC +MG5 aMC in-house MC +MG5 aMC 521:4+116156::49 737:7+3684::55 730:4+4614::89 728:3+9930::83 742:0+9300::13 777:8+6615::32 798:7+3263::26 766:2+1570::21 s = 13 TeV and MMHT 2014 PDFs. The default value of the factorization scale is f;0 = M=2, and the uncertainties are estimated through variations of this scale (and of the hard and soft scales s and h when applicable), as explained in the text. cross section, one sees that the e ect of the resummation corrections beyond NNLO are relatively small. As it was observed in the case of the ttH and ttW processes in [16{18], the scale uncertainty a ecting the nNLO result is very small compared to the NLO+NNLL scale uncertainty, and most likely underestimates the residual perturbative uncertainty at NNLO. Experimental collaborations reported measurements of the ttZ total cross section in combination with measurements of the ttW cross section [1{4], where the latter is the sum of the cross sections for ttW + and ttW production. We conclude this section by comparing our predictions for ttW and ttZ with experimental data. The ttW production cross section was evaluated by running the code developed in [18] with the scale choices and input parameters employed in the present work for ttZ production and described above. The results for ttZ and ttW production cross section at 8 and 13 TeV are summarized in gure 2 we follow the structure of gure 12 in [4] in order to compare graphically calculations with the corresponding experimental measurements. The experimental measurements at 8 TeV are taken from [2], while the experimental measurements at 13 TeV are taken from [4]. The green dots and cross-shaped \error bars" correspond to NLO calculations carried out with correspond instead to NLO+NNLL calculations. It is interesting to observe that, while predictions for the ttZ production cross section are in perfect agreement with the measurements at both 8 and 13 TeV, the predictions for the ttW cross section are slightly smaller than measurements for both collider energies. This observation holds for NLO and NLO+NNLL calculations alike. Of course this discrepancy should be taken with a grain of salt, and requires a more detailed discussion with the experimental collaborations. Moreover, we would like to stress that a fully exhaustive comparison between predictions and measurements should also account for the uncertainty associated to the choice of the PDFs and to the value of s. These two sources of uncertainty are not re ected in the error bars of gure 2. 136:7+1155::62 189:8+2244::58 130:7+64::99 203:9+1135::58 356:3+4339::75 182:2+2230::14 728:3+9930::83 341:0+2133::16 177:1+162:9:0 777:8+6615::32 s and pert. order 8 TeV NLO 8 TeV NLO 8 TeV NLO 8 TeV NLO+NNLL 8 TeV NLO+NNLL 8 TeV NLO+NNLL 13 TeV NLO 13 TeV NLO 13 TeV NLO 13 TeV NLO+NNLL 13 TeV NLO+NNLL 13 TeV NLO+NNLL and MMHT 2014 PDFs. The default value of the factorization scale is s = 8 and 13 TeV f;0 = M=2, and the uncertainties are estimated through variations of this scale (and of the resummation scales s and h when applicable). SM theory vs ATLAS data SM theory vs CMS data compared to the ATLAS measurements at 8 TeV [2] (data are represented by the light blue and pink bands). Right panel: total cross section at NLO (green cross) and NLO+NNLL (red cross) compared to the CMS measurements at 13 TeV [4] (light blue and pink bands). Di erential distributions In this section we obtain predictions for four di erential distributions which depend on the momenta of the nal state massive particles. The distributions are i) the distribution di erential with respect to the ttZ invariant mass, M , ii) the distribution di erential with respect to the tt invariant mass, Mtt, iii) the distribution di erential with respect to the 500 600 700 800 900 1000 1100 100 150 200 250 300 100 150 200 250 300 NLO (red band). The default factorization scale is chosen as f;0 = M=2, and the uncertainty bands are generated through scale variations as explained in the text. respect to the transverse momentum of the Z boson, pTZ . transverse momentum of the top quark,2 pt , and iv) the distribution di erential with Figure 3 compares the approximate NLO calculations, carried out with our in-house code, with the complete NLO calculations, carried out with MG5 aMC. We see that the approximate NLO calculations reproduce well the full NLO calculations. The lower part of each panel shows the ratio between the approximate NLO or complete NLO calculations and the central value of the NLO calculation. One can see that the approximate NLO scale uncertainty band is included in the NLO scale uncertainty band. Figure 4 repeats the same 2In this context we refer to top quark in a strict sense, and we evaluate this distribution by considering NLO w/o qg NLO w/o qg NLO w/o qg 500 600 700 800 900 1000 1100 100 150 200 250 300 100 150 200 250 300 tributions without the quark-gluon channel contribution (red band). All settings are as in gure 3. analysis but it compares approximate NLO calculations to NLO calculations without the quark-gluon channel contribution. As expected approximate NLO distributions and NLO distributions without the qg channel have the same shape and scale uncertainty bands of similar size. These two gures show that, for this choice of the factorization scale at least, soft emission corrections provide the bulk of the NLO corrections. Figure 5 provides the main result of this section. This gure compares NLO calculations to the distributions evaluated to NLO+NNLL accuracy. Roughly, one can say that the NLO+NNLL results fall in the upper part of the NLO scale uncertainty interval in each bin. The central value of the NLO+NNLL calculations is slightly larger than the central value of the NLO calculations in all bins shown. The scale uncertainty a ecting NLO 500 600 700 800 900 1000 1100 100 150 200 250 300 100 150 200 250 300 NLO calculation (red band). The uncertainty bands are generated through scale variations of f , s and h as explained in the text. the NLO+NNLL accuracy calculation, which is obtained by varying described above, is smaller than the NLO scale uncertainty band obtained by varying f . Results at NLO+NLL and NLO+NNLL accuracy are compared in gure 6. The main e ect of the NNLL correction with respect to the NLL ones is an increase of the central value of the bins in the tail of the M and Mtt distributions. The scale uncertainty bands turn out to be of similar size at NLO+NLL and NLO+NNLL in almost all bins shown. Figure 7 shows the ratio of distributions at various level of precision to the central value of the NLO+NNLL calculation in each bin. In particular, the blue band refers to NLO+NNLL distributions, the dashed red band to nNLO distributions and the dashed black band to distributions obtained from the NNLO expansion of the NLO+NNLL resummation. The NLO+NNLL expanded distributions di er from the NLO+NNLL distri500 600 700 800 900 1000 1100 100 150 200 250 300 100 150 200 250 300 corresponding NLO+NLL calculation (red band). The uncertainty bands are generated through scale variations. butions by NNLL resummation e ects of order N3LO and higher. These corrections can be as large as 5 to 10 % in all bins shown, and are particularly relevant at higher values of f . The di erence between the nNLO and the NLO+NNLL expanded to NNLO results is due to constant NNLO terms, which are formally of order N3LL. Both the NNLO expansion of the NLO+NNLL calculation and the nNLO calculation underestimate the scale uncertainty which one nds at NLO+NNLL accuracy, a fact which we already observed by looking at the predictions for the total cross section. The envelope of the two NNLO approximations (i.e. the black and red bands) spans almost all of the NLO+NNLL scale uncertainty interval in each bin, with the exception of the tail of the ptT distribution, where this envelope includes the NLO+NNLL scale uncertainty. through scale variations. Conclusions In the present work we carried out the resummation of soft gluon emission corrections to the associated production of a top-antitop quark pair and a Z boson. The resummation was studied in the partonic threshold limit z ! 1 and was implemented to NNLL accuracy. Numerical calculations of the total cross section and di erential distributions to NNLL accuracy were carried out by means of an in-house partonic Monte Carlo code which we developed for this work. The output of this code was matched with NLO calculations obtained from MG5 aMC. The nal outcome of this work is represented by the NLO+NNLL calculations of the total cross section and di erential distributions for the LHC operating at a center-of-mass energy of 13 TeV presented in the previous section. The code can be easily adapted to carry out phenomenological studies which include cuts on the top, antitop and/or Z boson momenta. With the choice of the factorization scale made in this work, we can conclude that the soft emission corrections to ttZ production evaluated to NNLL accuracy lead to a moderate increase of the total cross section and di erential distributions with respect to NLO calculations of the same observables. The residual perturbative uncertainty at NLO+NNLL accuracy, estimated by varying the soft, hard and factorization scales as explained in the text, is smaller than the NLO scale uncertainty, thus making our evaluations of the cross sections and di erential distributions in ttZ production the most precise results currently available in the literature. This work completes a series of papers devoted to the study of the associated production of a top pair and a colorless heavy boson to NLO+NNLL accuracy in the partonic threshold limit. In [18] the associated production of a top pair and a W boson was studied with the methods employed here for ttZ production, while the associated production of a top pair and a Higgs boson at NLO+NNLL accuracy was considered in [17]. In all cases the resummation was carried out in Mellin moment space. The hard and soft scales were chosen in the same way as in the traditional \direct QCD" approach. Codes for the numerical evaluation of the resummation are now available and tested for all of these three processes, and can be further employed in more speci c phenomenological studies, according to the interests of the experimental collaborations. Within such interactions with the experimental community, a detailed study of the uncertainty associated with the choice of the PDFs and to the value of s(MZ ), in the light of a comparison with the new measurements which are expected in the forthcoming months, would be particularly illuminating. At this stage, it would also be interesting to combine the NLO+NNLL calculations of ttW , ttZ and ttH production with the electroweak corrections for these processes [11, 12]. In addition to this, the inclusion of the decays of the heavy particles in the spirit of [44] is also possible. This would allow to put kinematic cuts on the momenta of the detected Acknowledgments The in-house Monte Carlo code which we developed and employed to evaluate the (di erential) cross sections presented in this paper was run on the computer cluster of the Center for Theoretical Physics at the Physics Department of New York City College of Technology. We thank P. Maierhofer and S. Pozzorini for their assistance with the program Openloops. This research was supported by the DFG cluster of excellence \Origin and Structure of the Universe". The work of A.F., G.O., and R.D.S. is supported in part by the National Science Foundation under Grant No. PHY-1417354. Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited. boson in pp collisions at p collisions at p [arXiv:1509.05276] [INSPIRE]. [1] CMS collaboration, Observation of top quark pairs produced in association with a vector [2] ATLAS collaboration, Measurement of the ttW and ttZ production cross sections in pp [3] ATLAS collaboration, Measurement of the ttZ and ttW production cross sections in nal states using 3.2 fb 1 of pp collisions at 13 TeV at the LHC, ATLAS-CONF-2016-003 (2016). boson in pp collisions at 13 TeV, CMS-PAS-TOP-16-017 (2016). [5] A. Lazopoulos, K. Melnikov and F.J. Petriello, NLO QCD corrections to the production of ttZ in gluon fusion, Phys. Rev. D 77 (2008) 034021 [arXiv:0709.4044] [INSPIRE]. [INSPIRE]. Hadroproduction at NLO accuracy in QCD with Parton Shower and Hadronization e ects, JHEP 11 (2012) 056 [arXiv:1208.2665] [INSPIRE]. [8] J.M. Campbell and R.K. Ellis, ttW + production and decay at NLO, JHEP 07 (2012) 052 [arXiv:1204.5678] [INSPIRE]. polarization in ttW production at the LHC, Phys. Lett. B 736 (2014) 252 [arXiv:1406.3262] [INSPIRE]. bosons at NLO in QCD: impact on ttH searches at the LHC, JHEP 02 (2016) 113 [arXiv:1507.05640] [INSPIRE]. hadroproduction in association with a top-quark pair, JHEP 09 (2014) 065 [arXiv:1407.0823] [INSPIRE]. to top-pair hadroproduction in association with heavy bosons, JHEP 06 (2015) 184 [arXiv:1504.03446] [INSPIRE]. [13] R. Rontsch and M. Schulze, Constraining couplings of top quarks to the Z boson in tt + Z production at the LHC, JHEP 07 (2014) 091 [Erratum ibid. 1509 (2015) 132] [arXiv:1404.1005] [INSPIRE]. [14] R. Rontsch and M. Schulze, Probing top-Z dipole moments at the LHC and ILC, JHEP 08 (2015) 044 [arXiv:1501.05939] [INSPIRE]. Notes Phys. 896 (2015) 1 [arXiv:1410.1892] [INSPIRE]. [INSPIRE]. top pair and a Higgs boson beyond NLO, JHEP 03 (2016) 124 [arXiv:1510.01914] production of a top pair and a Higgs boson at the LHC, JHEP 02 (2017) 126 [arXiv:1611.00049] [INSPIRE]. a W boson at next-to-next-to-leading logarithmic accuracy, JHEP 09 (2016) 089 [arXiv:1607.05303] [INSPIRE]. [INSPIRE]. [INSPIRE]. 658 (2003) 175 [hep-ph/0212259] [INSPIRE]. 734 (2006) 62 [hep-ph/0509141] [INSPIRE]. [23] F. Cascioli, P. Maierhofer and S. Pozzorini, Scattering Amplitudes with Open Loops, Phys. Rev. Lett. 108 (2012) 111601 [arXiv:1111.5206] [INSPIRE]. amplitudes with massive partons, Phys. Rev. Lett. 103 (2009) 201601 [arXiv:0907.4791] scattering amplitudes in non-abelian gauge theories, JHEP 11 (2009) 062 [arXiv:0908.3676] [INSPIRE]. 1889 [arXiv:1111.2034] [INSPIRE]. to calculate one-loop tensor integrals with up to six external legs, Comput. Phys. Commun. 180 (2009) 2317 [arXiv:0810.0992] [INSPIRE]. [32] P. Mastrolia, E. Mirabella and T. Peraro, Integrand reduction of one-loop scattering amplitudes through Laurent series expansion, JHEP 06 (2012) 095 [Erratum ibid. 1211 (2012) 128] [arXiv:1203.0291] [INSPIRE]. One-loop Massive Amplitudes from Integrand Reduction via Laurent Expansion, JHEP 03 (2014) 115 [arXiv:1312.6678] [INSPIRE]. [34] T. Peraro, Ninja: Automated Integrand Reduction via Laurent Expansion for One-Loop Amplitudes, Comput. Phys. Commun. 185 (2014) 2771 [arXiv:1403.1229] [INSPIRE]. hadronic collisions, Nucl. Phys. B 478 (1996) 273 [hep-ph/9604351] [INSPIRE]. in SCET vs. direct QCD: Higgs production as a case study, JHEP 01 (2015) 046 [arXiv:1409.0864] [INSPIRE]. [INSPIRE]. [INSPIRE]. [6] A. Lazopoulos , T. McElmurry , K. Melnikov and F. Petriello , Next-to-leading order QCD corrections to ttZ production at the LHC , Phys. Lett . B 666 ( 2008 ) 62 [arXiv:0804.2220] [7] M.V. Garzelli , A. Kardos , C.G. Papadopoulos and Z. Trocsanyi , tt W + [9] F. Maltoni , M.L. Mangano , I. Tsinikos and M. Zaro , Top-quark charge asymmetry and [10] F. Maltoni , D. Pagani and I. Tsinikos , Associated production of a top-quark pair with vector [11] S. Frixione , V. Hirschi , D. Pagani , H.S. Shao and M. Zaro , Weak corrections to Higgs [12] S. Frixione , V. Hirschi , D. Pagani , H.S. Shao and M. Zaro , Electroweak and QCD corrections [15] T. Becher , A. Broggio and A. Ferroglia , Introduction to Soft-Collinear E ective Theory, Lect. [16] A. Broggio , A. Ferroglia , B.D. Pecjak , A. Signer and L.L. Yang , Associated production of a [17] A. Broggio , A. Ferroglia , B.D. Pecjak and L.L. Yang , NNLL resummation for the associated [18] A. Broggio , A. Ferroglia , G. Ossola and B.D. Pecjak , Associated production of a top pair and [24] A. Denner and S. Dittmaier , Reduction of one loop tensor ve point integrals, Nucl . Phys . B [25] A. Denner and S. Dittmaier , Reduction schemes for one-loop tensor integrals, Nucl . Phys . B [26] A. Denner and S. Dittmaier , Scalar one-loop 4-point integrals, Nucl . Phys . B 844 ( 2011 ) 199 [27] A. Denner , S. Dittmaier and L. Hofer , Collier: a fortran-based Complex One-Loop LIbrary in Extended Regularizations, Comput . Phys. Commun . 212 ( 2017 ) 220 [arXiv:1604.06792] [19] H.T. Li , C.S. Li and S.A. Li , Renormalization group improved predictions for ttW production at hadron colliders , Phys. Rev. D 90 ( 2014 ) 094009 [arXiv:1409.1460] [INSPIRE]. [20] J. Alwall et al., The automated computation of tree-level and next-to-leading order di erential cross sections and their matching to parton shower simulations , JHEP 07 ( 2014 ) [21] A. Ferroglia , M. Neubert , B.D. Pecjak and L.L. Yang , Two-loop divergences of scattering [22] A. Ferroglia , M. Neubert , B.D. Pecjak and L.L. Yang , Two-loop divergences of massive [28] G. Cullen et al., Automated One -Loop Calculations with GoSam, Eur. Phys. J. C 72 ( 2012 ) [29] G. Cullen et al., GOSAM-2 . 0 : a tool for automated one-loop calculations within the Standard Model and beyond, Eur. Phys. J. C 74 ( 2014 ) 3001 [arXiv:1404.7096] [INSPIRE]. [30] T. Binoth , J.P. Guillet , G. Heinrich , E. Pilon and T. Reiter , Golem95 : A Numerical program [31] P. Mastrolia , G. Ossola , T. Reiter and F. Tramontano , Scattering AMplitudes from Unitarity-based Reduction Algorithm at the Integrand-level , JHEP 08 ( 2010 ) 080 [33] H. van Deurzen , G. Luisoni , P. Mastrolia , E. Mirabella , G. Ossola and T. Peraro , Multi-leg [35] S. Catani , M.L. Mangano , P. Nason and L. Trentadue , The Resummation of soft gluons in [36] M. Bonvini , Resummation of soft and hard gluon radiation in perturbative QCD , [37] M. Bonvini and S. Marzani , Resummed Higgs cross section at N3LL , JHEP 09 ( 2014 ) 007 [38] L.A. Harland-Lang , A.D. Martin , P. Motylinski and R.S. Thorne , Parton distributions in the LHC era: MMHT 2014 PDFs, Eur . Phys. J. C 75 ( 2015 ) 204 [arXiv:1412.3989] [INSPIRE]. [39] M. Bonvini , S. Forte , M. Ghezzi and G. Ridol , Threshold Resummation in SCET vs. Perturbative QCD: An Analytic Comparison, Nucl . Phys . B 861 ( 2012 ) 337 [40] M. Bonvini , S. Forte , G. Ridol and L. Rottoli , Resummation prescriptions and ambiguities [41] A. Ferroglia , B.D. Pecjak , D.J. Scott and L.L. Yang , QCD resummations for boosted top production, PoS(TOP2015)052 [arXiv:1512 .02535] [INSPIRE]. [42] B.D. Pecjak , D.J. Scott , X. Wang and L.L. Yang , Resummed di erential cross sections for top-quark pairs at the LHC , Phys. Rev. Lett . 116 ( 2016 ) 202001 [arXiv:1601.07020] [43] J.C. Collins , D.E. Soper and G.F. Sterman , Factorization of Hard Processes in QCD, Adv. Ser. Direct. High Energy Phys . 5 ( 1989 ) 1 [hep-ph/0409313] [INSPIRE]. [44] A. Broggio , A.S. Papanastasiou and A. Signer , Renormalization-group improved fully di erential cross sections for top pair production , JHEP 10 ( 2014 ) 98 [arXiv:1407.2532]


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2FJHEP04%282017%29105.pdf

Alessandro Broggio, Andrea Ferroglia, Giovanni Ossola, Ben D. Pecjak, Ray D. Sameshima. Associated production of a top pair and a Z boson at the LHC to NNLL accuracy, Journal of High Energy Physics, 2017, 105, DOI: 10.1007/JHEP04(2017)105