Associated production of a top quark pair with a heavy electroweak gauge boson at NLO+NNLL accuracy

We perform threshold resummation of soft gluon corrections to the total cross sections and the invariant mass distributions for production of a top-antitop quark pair associated with a heavy electroweak boson $V = W^+$, $W^-$ or $Z$ in $pp$ collisions at the Large Hadron Collider. The resummation is carried out at next-to-next-to-leading-logarithmic (NNLL) accuracy using the direct QCD Mellin space technique in the three-particle invariant mass kinematics. It is found that for the $t\bar t Z$ process the soft gluon resummation introduces significant corrections to the next-to-leading order (NLO) results. For the central scale equal to the $t\bar t Z$ invariant mass the corrections reach nearly 30%. For this process, the dominant theoretical uncertainty of the cross section due to the scale choice is significantly reduced at the NLO+NNLL level with respect to the NLO results. The effects of resummation are found to be less pronounced in the $t\bar t W^{\pm}$ case. The obtained results are compared to recent measurements performed by CMS and ATLAS collaborations at the LHC.


Introduction
The measurements of associated production of a massive vector boson with a top-antitop quark pair at the LHC [1][2][3][4][5][6][7] provide an important test of the Standard Model (SM). Together with the associated Higgs boson production with a top-quark pair, they belong to a class of processes with the heaviest final states which can be precisely studied at the LHC. Such studies command particular attention as a means to indirectly search for signals of physics Beyond the Standard Model (BSM). Additional, they form dominant background in direct BSM searches, as well as to SM measurements, especially the associated Higgs boson production process. It is therefore necessary to know the theoretical predictions for pp → ttV , V = W + , W − , Z with high accuracy.
Over the years, there has been a great effort to improve the theoretical description of the pp → ttV process. Next-to-leading-order (NLO) QCD corrections were calculated [8][9][10][11][12][13][14][15][16][17] and matched to partons showers [13,18,19]. The electroweak corrections and the combined electroweak corrections with the QCD corrections are also known [20][21][22]. In the light of the full next-to-next-to-leading order (NNLO) QCD calculations being currently out of reach, it is useful to systematically consider at least some part of the higher order corrections to improve the theoretical precision. This can be achieved using resummation techniques for corrections originating from emission of soft gluons. This type of emission happens in the presence of kinematical constraints, where the phase space available for emission of real gluons is restricted. As the kinematical limit is approached, the corrections are dominated by large logarithmic contributions with an argument of the logarithms directly related to the distance to the limit. The observable in question for which the predictions are obtained and the kinematics in which it is considered determines then the exact form of the logarithms.
Two popular approaches to perform soft gluon resummation are either direct calculations in QCD or an application of an effective field theory, in this case soft-collinear effective theory (SCET). Although the physics that is described is obviously the same, and the perturbative accuracy which can be reached is also formally the same, the two approaches differ at the technical level resulting in different treatment of subleading corrections beyond the formal accuracy. However, in practice these corrections can introduce non-negligible effects. Therefore it is valuable to perform calculations using both techniques, firstly as a completely independent check of the calculations and secondly as an indication of the size of subleading effects.
As for any process for which production rates are small, also for the associated toppair production with a heavy boson the first quantity which can be studied with higher precision is the total cross section. The higher order corrections receive then potentially large contributions from soft gluon emission in the threshold limit, i.e. the partonic center of mass energyŝ approaches the energy needed to produce the final state with a given characteristics. For the process pp → ttH soft gluon resummation has been performed both using direct QCD [23][24][25][26] and SCET [28,29] methods. While the next-to-leadinglogarithmic (NLL) calculations [23] were carried out in the absolute threshold limit,ŝ → (2m t + m H ) 2 , the later calculations [24-26, 28, 29] opted for the invariant mass threshold limit, i.e.ŝ → Q 2 with Q 2 = (p t + pt + p H ) 2 . The resummed predictions are now known at the next-to-next-to-leading logarithmic (NNLL) accuracy in both approaches and are matched to the full NLO results to include all available information on the process. In the case of associated top-pair production with a heavy gauge boson, W + , W − or Z, NLO+NNLL predictions obtained within the SCET framework are already available [30][31][32], whereas for calculations in the direct QCD approach only NLO+NLL results have been communicated so far [27]. Here we close this gap and report on soft gluon resummation in this approach at the NLO+NNLL accuracy for the process pp → ttV . Our calculations rely on the techniques described in [25]. We present numerical results for the total cross sections and the invariant mass distributions as well as comment on the comparison between our results and those of [30][31][32].
The paper is structured as follows: in section 2 we review the direct QCD approach applied before in the calculations for the process pp → ttH and now adapted to the pp → ttV case. The numerical results are presented and discussed in section 3. The conclusions and the summary of our work can be found in section 4.
2 NNLL resummation in the triple invariant mass kinematics for 2 → 3 processes with two massive coloured particles in the final state In the following, we use the direct QCD approach to resummation of soft gluon corrections at threshold in Mellin space. In particular, we consider the threshold limit in the three particle invariant mass kinematics. The Mellin transformation of the differential cross section dσ pp→ttV /dQ 2 is then performed w.r.t. the variable ρ = Q 2 /S. Resummation provides a systematic treatment of logarithmic terms of the form α n s [log m (1 − z)/(1 − z)] + , with m ≤ 2n − 1 and z = Q 2 /ŝ which appear at all orders of the perturbative expansion in α s . These logarithms then turn into logarithms of the Mellin moment N in Mellin space, where the threshold limit z → 1 corresponds to the limit N → ∞.
We use the same framework as developed in [25] and in the following we consider a process ij → klV , where i, j are massless coloured partons, k, l two massive quarks and V a massive colour-singlet particle. The collective argument {m 2 } denotes all masses entering the calculations. The resummed partonic cross section up to NNLL accuracy can be written as: where H R ,Ū R , U R andS R are colour matrices and the trace is taken over colour space. ∆ i and ∆ j represent the logarithmic contributions from the (soft-)collinear gluon emission from the initial state partons. They are universal functions, depending only on the emitting parton, and can be found for example in [33,34] up to NLL and in [35] up to NNLL level. The termŪ RSR U R originates from a solution of the renormalization group equation of the soft function, which describes the soft wide angle emission. It consists of the soft function evolution matricesŪ R and U R , as well asS R which plays the role of a boundary condition of the renormalization group equation. In general the evolution matrices are given by path-ordered exponentials of a soft anomalous dimension matrix in the colour space. If the matrix in the argument of the path ordered exponential is diagonal, it reduces to a sum over simple exponentials. All colour matrices in eq. (1) are expressed in the basis in which the one loop soft-anomalous dimension Γ is diagonal. The diagonalization is achieved by a colour basis transformation where R is the diagonalization matrix and λ ij→klV . Correspondingly, all colour matrices in eq. (1) carry a subscript R. Γ ij→klV has to be known up to Γ (1) ij→klV to perform resummation with NLL accuracy and up to Γ (2) ij→klV for NNLL. The one-loop soft anomalous dimension can be found in [23] while the two loop soft anomalous dimension was derived in [36,37].
In practice, we start with a description of the colour structure of the ttW and ttZ processes in the s-channel colour basis, {c q I } and {c g I } given by their basis vectors: Since at leading order ttW state is produced via qq channel we only need the {c q I } basis for its description, whereas both {c q I } and {c g I } basis are needed to describe the ttZ production via the qq and the gg channels.
The functionS R is obtained by transforming the purely eikonal functionS ij→klV , calculated in the s-channel colour basis and NLL accuracy requires knowledge ofS ij→klV while NNLL accuracy requiresS (1) ij→klV . Since the one-loop soft anomalous dimension is in general non diagonal in the triple invariant mass kinematics, in order to calculate the soft function evolution matrices up to NLL we use the diagonalization method of [38]. In this way the path ordered exponentials reduce to a sum over simple exponentials andŪ RSR U R at NLL is given bȳ where λ is defined as and Resummation up to NNLL encounters additional complexity due to the non-commutativity of Γ (1) ij→klV and Γ (2) ij→klV . Therefore we employ the method detailed in [39,40] to recast the soft function evolution matrices into simple exponentials. This results in with and The hard function H R describes the hard scattering contributions and absorbs off-shell effects. It is independent of N and given by a matrix in colour space, which is then also transformed into the R colour space The hard function matrix can be calculated perturbatively: In order to perform resummation up to NNLL knowledge of H ij→klV can be calculated from the LO cross section of the process, the next order includes N -independent non-logarithmic contributions originating from virtual loops, real collinear terms and the evolution matrices U R andŪ R . The virtual contributions are extracted from the PowHel code [11,18,41] and projected on the colour basis. Following the method proposed in [42,43] the real terms are derived from the infrared limit of the real corrections.
The resummed cross sections of different accuracy denoted by "res" in the following are matched with the full NLO cross section according to where "res" = N(N)LL and "matched" = NLO + N(N)LL for the N(N)LL resummed results matched to NLO. The moments of the parton distribution functions f i/h (x, µ 2 F ) are defined in the standard way and dσ (res) ij→klV /dQ 2 | (NLO) represents the perturbative expansion of the resummed cross section truncated at NLO. The inverse Mellin transform (17) is evaluated numerically using a contour C in the complex-N space according to the "Minimal Prescription" method proposed in Ref. [33].
Apart from the NLO+NLL and NLO+NNLL results we also calculate the NLL result improved by including contributions of order O(α s ) terms in H R andS R and matched to NLO which we refer to as NLO+NLL w C . The resummed partonic cross section at this accuracy is given by: where R .

Numerical results for the pp → ttV processes at NLO+NNLL accuracy
In this section we present our resummed results with different levels of precision e.g. NLL, NLL w C and NNLL matched to NLO. They include distributions differential in Q as well as total cross sections which were calculated by integrating over Q. The resummed results were obtained with two independently developed in-house codes, while the NLO cross sections were calculated with MadGraph5_aMC@NLO [19] for differential distributions and total cross sections, and with PowHel [11,18,41] for NLO total cross sections without the contributions from qg channels. In the calculations we use the PDF4LHC15_30 parton distribution function (pdf) set [44][45][46][47][48][49] and input parameters according to the Higgs Cross Section Working Group (HXSWG) recommendations [50], i.e. m H = 125 GeV, m t = 172.5 GeV, m W = 80.385 GeV, m Z = 91.188 GeV, G F = 1.1663787 · 10 −5 GeV −2 . This is the same choice as the one made in the HXSWG Yellow Report 4 [51], so that we can reproduce the NLO values of the ttV cross sections listed there and compare our new resummed predictions to them. In accordance with the Yellow Report setup, in the calculations of the ttW cross sections the CKM matrix is taken diagonal. NLO pdf sets are used for NLO, NLL matched to NLO and NLL w C matched to NLO results, while NNLO pdf sets for NNLL matched to NLO. The pdf error is only calculated for the NLO cross sections for technical simplicity, since resummation will influence the pdf error only minimally. Two different choices for the central factorisation and renormalisation scales are used for most of the results throughout the section. The first choice is µ 0 = µ F,0 = µ R,0 = Q which is the natural scale for the threshold and kinematics in the resummation, while the second scale is µ 0 = µ F,0 = µ R,0 = M/2, which is often used in fixed order calculations. The scale error is calculated using the seven point method by taking the maximum and minimum of the scale variations (µ F /µ 0 , µ R /µ 0 ) = (0.5, 0.5), (0.5, 1), (1, 0.5), (1, 1), (1, 2), (2, 1), (2, 2). As mentioned in the previous section the hard function includes virtual loop corrections which were extracted from PowHel [11,18,41]. Analytical relations between the basis colour configurations of the colour flow basis used in PowHel and the basis vectors of the s-channel colour basis allow us to extract the full matrix H R,IJ . The colour summed results were then compared to the standalone MadLoop implementation from MadGraph5_aMC@NLO [19].

Total cross sections
At first we compare the total cross section of the full NLO calculation with our resummed result expanded in α s up to NLO to analyse how well the full NLO cross section can be reproduced. The relatively large significance of the qg channel especially for the scale uncertainties was shown first in [30] for ttW and later in [31,32] for ttW and ttZ. Since the qg channel first appears at NLO, no resummation is performed for this channel. Therefore one has to compare the NLO cross section without the qg channel with the resummed result expanded in α s to judge the quality of the approximation provided by the expansion. The resummed cross sections matched to NLO include the qg channel through the matching to the full NLO calculations. In figure 1 we compare of the full NLO cross section, the NLO cross section without the qg channel and the expansion of the resummed cross section as a function of µ/µ 0 = µ F /µ F,0 = µ R /µ R,0 for ttW at √ S = 13 TeV for the two different scale choices µ 0 = Q and µ 0 = M/2. The corresponding comparison for ttZ is shown in figure 2.
In all cases the NLO cross section without the qg channel is much better approximated NLO (no gq)    as the precision of the calculations increases. Another manifestation of the same effect originating from soft gluon corrections is the decrease in the scale uncertainties calculated for each specific scale choice which is also progressing with increasing precision of the theoretical predictions. This trend is much stronger for ttZ production than for ttW due to the gg channel contributing to the LO and, correspondingly, to the resummed cross section. As the gluon radiate more than quarks, resummation has more relevance for the gg production channel than for qq or qq channels. Correspondingly, we see a decrease in the ttZ scale uncertainty of about 30-40% when increasing the precision from NLO to NLO+NNLL. The ttW cross section scale uncertainty is reduced by 20-30% with the exception of the upwards uncertainty for µ = M/2 which does not receive any significant improvement.
As already noted, the NLO+NNLL predictions with the central scale varied between M/2, Q/2 and Q are closer in value than the corresponding NLO predictions. The NLO+NNLL ttZ results are particularly stable w.r.t. scale variation. Correspondingly, the NNLL K-factors, ranging from 1.04 to 1.29, cf. tables 1 and 2, have to compensate for the scale dependence of the NLO results. Due to limited corrections from resummation for the quark-initiated ttW process the NNLL K-factors are smaller, ranging from 1.01 to 1.07. The scale dependence of the NLO is strongly influenced by the scale dependence of the qg channel (see e.g. [15]) which is formally subleading and not resummed here. Therefore resummation for the qq channel does not fully compensate the scale dependence of the NLO and leads only to moderate improvements at NLO+NNLL. Note that it is also possible to obtain soft gluon approximation of the NNLO corrections by expanding the resummed cross section. We performed such studies for the ttH production at the LHC [25], where we added this approximation to the full NLO result, resulting in the NNLO Approx. predictions. We found them to be fully consistent with the resummed NLO+NNLL cross sections. For the ttV processes, the approximate NNLO predictions were already presented in [31,32]. Since here we are interested in the resummed results, we refer the readers interested in NNLO Approx. predictions to these publications.
The observed improvement in the stability of the predictions w.r.t. scale variation at NLO+NNLL for the ttZ process is akin to the improvement for the ttH process [25]. Similarly, we are encouraged to combine the predictions for our three representative scale choices according to the envelope method proposed by the HXSWG [52]. This way we can obtain theoretical predictions with the most conservative estimate of the scale error. The corresponding result for the ttZ production at 13 TeV is: and at 14 TeV σ ttZ NLO+NNLL = 1045 +8.8%+3.1% −9.9%−3.1% fb. The first uncertainty originates from the scale variation and is calculated using the envelope method, whereas the second one is the pdf+α s uncertainty. These values are in a very good agreement with the NLO results obtained for the scale choice µ 0 = µ F,0 = µ R,0 = M/2, justifying this common choice to obtain theory predictions for this process.
The same treatment can be applied to the ttW + and ttW − production resulting in for √ S = 13 TeV, and at √ S = 14 TeV σ ttW − NLO+NNLL = 224 +26.4%+3.6% −16.4%−3.6% fb, where again the first uncertainty originates from the scale variation and the second is the pdf+α s uncertainty. Due to a worse agreement between cross section predictions for the different central scale choices this treatment leads to a larger uncertainty than the uncertainty for the common choice of µ = M/2.
To further study the scale uncertainty of the total cross sections we show the dependence of the ttW and ttZ cross sections on the choice µ = µ F = µ R in figures 5 and 6. For the associated production of the top quark pair with a W boson, the sum of the ttW + and ttW − production is presented, since the two processes possess a very similar scale dependence. In figure 5, a slight reduction in the scale dependence can be seen with the dominant contribution brought by NLO+NLL w C result, indicating the importance of contributions of hard origin. In addition, a mild increase of the dependence can be seen for the significantly small scales µ 0.3M/2, which can potentially be attributed to the missing quark emission contribution. However, the scale at which this increase happens is not physically motivated and therefore of no relevance in practical studies. Separating the µ F and µ R dependence, i.e. varying the µ F and µ R while keeping µ R and µ F fixed respectively, leads to the conclusion that the ttW scale dependence is almost solely driven by the µ R dependence, cf. figures 7 and 8. For the ttZ production process a more significant reduction in the dependence on µ = µ F = µ R can be seen in figure 6. Similarly to the ttW process, the dominant reduction in the uncertainty can also be attributed to the inclusion of constant contributions in N from the hard and soft functions contained in the difference between the NLO+NLL and NLO+NLL w C results. However, a significant further reduction in the scale dependence originates from the resummation at NLL level and beyond. Additionally, the same increase of the dependence can be seen for the ttZ process at low scales, but again this effect concerns scales which choice is not physically motivated. The figures also illustrate that if we had attributed the uncertainty of the cross section to the scale variation for µ F = µ R , the scale uncertainty would have been drastically reduced, even as low as to approximately 1% for the µ = Q choice. In contrast to the ttW process, the ttZ dependence on µ F = µ R appears to be an effect of cancellations between dependencies on µ F and µ R . Taken separately they show an opposite behaviour, see figures 9 and 10. This behaviour is in fact very similar to the one observed for the process pp → ttH, which also receives significant contributions from the gg channel at LO.

Invariant mass distributions
Our total cross section predictions are obtained by integrating over invariant mass distributions dσ/dQ 2 . Note that these are the only distributions for which one has got a full control of the resummed contributions while performing threshold resummation in the invariant mass limitŝ → Q 2 . The NLO+NNLL distributions in Q for the two scale choices µ 0 = Q, µ 0 = M/2 for the ttW and ttZ processes are presented in figures 11 and 12, respectively. Apart from the scale choice, the size of the NNLL corrections depends now also on Q. In the ttW case, however, this dependence is moderate and the corrections do not exceed 10%, cf. left plot in figure 11. The corrections to the ttZ invariant mass distribution, on the other hand, show much stronger Q dependence. Figure 12 illustrates that the NNLL corrections can reach up to 40% for the µ 0 = Q scale choice which is a much higher value than the 29% reported for the total cross section in table 1. Similarly as in the case of total cross sections, also for differential distributions inclusion of the NNLL corrections results in a much better agreement between theoretical predictions obtained with various scale choices, and in consequence leads to stabilization of the predictions, see figure 13.

Comparison with other NLO+NNLL predictions in the literature
The NLO+NNLL predictions for the associated ttW and ttZ production calculated in the SCET framework are available [30][31][32]. By comparing our results obtained using the direct QCD approach, we can not only deliver an independent check of the previously published result but also gain insights on the size of the subleading i.e. below formal accuracy effects which are treated differently by the two methods. In order to perform the corresponding comparisons we used the same values of parameters and the same pdf sets as in the above mentioned papers. It has to be noted though, that the scale choices made to obtain results reported in this paper and [30][31][32] are not equivalent. While our resummed expressions depend on µ F and µ R , the formulas obtained in the SCET formalism contain dependence on the hard and soft scales µ h and µ s , as well as µ F . In particular, Ref. [30] uses the choice µ F = µ h = Q and a minimizing procedure to set µ s . Nevertheless, we find a very good agreement with our results calculated using the scale choice µ F = µ R = Q. Specifically, the authors [30] obtain σ NLO+NNLL = 332.99 +5% −4% fb for the ttW + production and σ NLO+NNLL = 169.86 +5% −4% fb for the ttW − production at √ S = 13 TeV, while we have σ NLO+NNLL = 331 +8.9% −8.6% fb and σ NLO+NNLL = 170 +8.8% −8.6% fb, correspondingly. Ref. [31] also reports the NLO+NNLL predictions for the ttW + and ttW − processes at the LHC. Contrary to [30], the µ s scale in [31] is chosen in such a way as to mimic the scale of soft radiation in the Mellin-space framework, i.e. µ s = Q/N . As pointed out in [25], only one choice µ F = µ h = Q and µ s = Q/N in the scale setting procedure of [31] directly corresponds to setting µ = µ F = µ R = Q in our results. With this choice, and using the same pdf and input parameter setup as in [31], we obtain σ NLO+NNLL = 328.6 +29.2 −28.1 fb for the ttW + and σ NLO+NNLL = 171.2 +15.0 −14.6 fb for the ttW − process at √ S=13 TeV, to be compared with σ NLO+NNLL = 333 +14.9 −12.4 fb and σ NLO+NNLL = 173.1 +7.7 −6.6 fb reported in [31]. The differences between the central values obtained by these two different calculations amount thus to 1%. Note that, similarly to [30], scale errors cannot be directly compared due to different methods used for calculating them. In particular, as explained above, our estimates of the scale error are calculated using the seven point method.
The NLO+NNLL predictions for the ttZ production process reported in table 2 of [32] use yet another scale setup with µ F = Q/2, µ h = Q and µ s = Q/N . Therefore these predictions cannot be directly compared with ours. Nevertheless, in figure 14 we present our NLO+NNL, NLO+NLL (w C), NLO+NNLL results as well as the NLO cross section for the ttZ production using the same pdf and input parameter setup as in [32] and choosing µ 0 = µ F,0 = µ R,0 = Q/2. Notably, within the range of the scales considered, we do not see the rising behaviour of the NLO+NNLL cross sections with the growing scale as in figure 1 of [32], instead a relatively flat one, especially at µ ≤ Q/4. Also, as discussed in Section 3.1, compared to NLO+NLL predictions our NLO+NNLL results show decreased dependence on the scale variation. Regarding numerical values for √ S=13 TeV, we obtain σ NLO+NNLL = 819.7 +58. 6 −65.0 fb at the scale µ = µ F = µ R = Q/2 whereas [32] reports σ NLO+NNLL = 777.8 +61. 3 −65.2 fb. For comparison, the difference between the NLO+NNLL total cross section for the ttH process obtained in the SCET formalism with µ F = µ h = Q, µ s = Q/N and the NLO+NNLL total cross section obtained using the direct QCD method with µ F = µ R = Q [25] amounts to 2.5% at √ S=13 TeV when the same pdf and input parameter setup is used. Given this, and the values we can read off from figure 1 of [32], we expect the difference between the two cross sections quoted here to be significantly impacted by the scale setting procedure.

Comparison with the ttV total cross section measurements at the LHC
The currently most precise measurements of ttW ± and ttZ cross sections in pp collisions at √ S = 13 TeV were recently published by CMS [6] and ATLAS collaborations [7] . The data samples correspond to an integrated luminosity of 35.9 fb −1 and 36.1 fb −1 , respectively. While ATLAS measures ttW ± and ttZ production cross sections simultaneously, CMS provides numerical values for individual measurements and a figure with the results of a simultaneous fit. In table 3 we compare the results of these measurements with the central theoretical predictions of this paper at the NLO+NNLL accuracy (19), (20) and (21), to which electroweak corrections reported in [51] are added. For each process, the EW corrections are estimated from the values of the relative EW corrections listed in table 40 of [51] (−0.2% for ttZ, −3.5% for ttW + , −2.6% for ttW − , −3.2% for ttW ), and the corresponding NLO QCD cross sections calculated using the envelope method and the NLO values listed in table 1 in this paper. The QCD uncertainties are applied also to the EW correction effects. In this way we provide theoretical predictions which include   Table 3: Results of experimental measurements by CMS [6] and ATLAS [7] collaborations of total cross sections σ for pp → ttW + /W − /Z at √ S = 13 TeV compared to theory predictions at NLO+NNLL accuracy given in eqs. (19), (20) and (21) with added the electroweak corrections as reported in [51]. The scale and the pdf+α s errors correspond to QCD cross sections.
The comparison shows good agreement between theory and data within errors. The largest relative difference between theory and experiment is found for the CMS measurement of the ttW + production. Even with the most conservative estimate of scale error provided by the envelope method, the overall theory error are smaller than the current experimental errors. It is interesting to note a general trend for the theoretical results to be lower or comparable to the measured experimental values. The same conclusions  Figure 16: NLO+NNLL and NLO predictions for the total ttZ and ttW cross sections, using the envelope method as described in the text, and with added electroweak corrections reported in [51], compared to the CMS [6] (left plot) and ATLAS [7] (right plot) measurements.
In addition, in figures 15 and 16 we show the NLO predictions and our NLO+NNLL predictions for the ttZ and ttW processes, to which we add the electroweak corrections computed in [51]. The NLO+NNLL results are marked by full lines for the central values and darker shaded bands for the errors, while for NLO dashed lines and light shaded bands are used correspondingly. In the figures we also plot the values of the corresponding cross sections measured using a combined fit by the CMS [6] and ATLAS [7] collaborations, together with their confidence level (CL) contours. The theory errors are estimated by adding the the scale errors and pdf+α s errors of the QCD cross sections. In figure 15 the theory values for the central scale µ 0 = M/2 are shown. The NLO QCD+EW results with this scale choice have been reported in the Yellow Report [51], and are taken as the benchmark theory predictions by the experiments. Since, however, the choice of a single fixed scale at the NLO accuracy may lead to underestimation of the theoretical uncertainty, in figure 16 we display also the NLO+NNLL and NLO predictions that include the whole span of the central scales from M/2 to Q. The NLO+NNLL results (19), (20) and (21) with added EW corrections combine predictions for various central scale choices, thus yielding more conservative error estimates. The two-dimensional analyses visualised in figures 15 and 16 confirm good agreement between the theory predictions and the measurements by the ATLAS and the CMS collaborations. Since the NLO+NNLL total cross sections are higher than the NLO ones, the NNLL calculations result in bringing the central values of the theoretical predictions closer to the experimentally measured cross sections. In the case of analysis with more conservative error estimates (figure 16), this distance gets reduced by as much as a factor of two for the ttZ production. The theoretical accuracy for ttZ in this conservative approach is equally well improved due to inclusion of the soft gluon resummation effects in the NNLL approximation, also by around a factor of two w.r.t. the NLO result.

Summary
In this paper soft gluon corrections were calculated to ttV production in association with a heavy electroweak boson V , V = W ± or Z in pp collisions. The calculations were performed in the three particle invariant mass kinematics through NNLL accuracy, and the results were matched to existing NLO results. Resummation was achieved using the direct QCD approach in the Mellin space. We computed invariant mass distributions and the total cross sections, obtained by integration of these distributions. In particular, we calculated NLO+NNLL total cross sections for the LHC collisions at √ S = 13 TeV and √ S = 14 TeV. The scale uncertainty of these predictions were estimated by independent variation of the factorisation and renormalisation scales around a central scale µ 0 , using the seven point method. Three different choices of the central scale µ 0 were assumed: µ 0 = M/2, µ 0 = Q/2 and µ 0 = Q, where M = 2m t + m V is the absolute threshold energy, and Q is the invariant mass of the ttV system.
The effect of soft gluon corrections was found to be more important for the ttZ than for the ttW ± production. This was expected, as the LO amplitudes for ttW ± production are driven by quark scattering, while for ttZ the two gluon production channel is dominant, and stronger gluon radiation occurs due to higher colour charges. For the ttZ production we observed a substantial improvement of the theoretical accuracy due to inclusion of the soft gluon corrections. First of all, the results are much more stable w.r.t. the central scale choice: at NLO+NNLL the total cross section increases by only 3% when µ 0 decreases from Q to M/2, while the corresponding increase is 28% at NLO. Moreover, for the fixed central scale, the dominant theory errors from the scale choice uncertainty decrease by 29% -38% by going from NLO to NLO+NNLL. A conservative estimate of the theoretical accuracy obtained as an envelope over results for various scale choices and their errors is improved by up to a factor of two by performing the NNLL soft gluon resummation. As in the case of the Higgs boson production with association of tt quarks, our results are compatible with NLO prediction for the central scale choice µ 0 = µ F,0 = µ R,0 = M/2 justifying that common choice at least for the ttZ process.
The obtained results were compared to the existing predictions for the ttW ± and ttZ cross sections at NLO+NNLL that were calculated in the SCET framework. In order to perform a meaningful comparison, we computed the cross sections employing the same sets of parton distribution functions and the input parameters as in those papers. For equivalent scale choice setups, our NLO+NNLL predictions and the cross sections calculated using the SCET framework agree well.
Finally, the theoretical estimates of ttV total cross sections were compared to the latest ATLAS and CMS measurements at √ S = 13 TeV. A good agreement was found between theory and data. In a two dimensional analysis of ttW and ttZ cross sections, the combined experimental data differ by about one standard deviation from the results of this paper. In comparison with the NLO predictions, the NLO+NNLL calculations result in theoretical predictions with central values closer to the measured experimental cross sections. The errors of the NLO+NNLL predictions are in general smaller than the current experimental errors.