Combined NLO EW and QCD corrections to off-shell ttW production at the LHC

The high luminosity that will be accumulated at the LHC will enable precise differential measurements of the hadronic production of a top–antitop-quark pair in association with a W boson. Therefore, an accurate description of this process is needed for realistic final states. In this work we combine for the first time the NLO QCD and electroweak corrections to the full off-shell ttW production at the LHC in the three-charged-lepton channel, including all spin correlations, non-resonant effects, and interferences. To this end, we have computed the NLO electroweak radiative corrections to the leading QCD order as well as the NLO QCD corrections to both the QCD and the electroweak leading orders.


Introduction
The hadronic production of top-antitop pairs in association with a W boson is an interesting process to investigate at the Large Hadron Collider (LHC), as it represents an important probe of the Standard Model (SM) as well as a window to new physics.
This process is one of the heaviest signatures measurable at the LHC. It gives access to the top-quark coupling to weak bosons and to possible deviations from its SM value [1][2][3]. Due to the absence of a neutral initial state at a lower perturbative order than next-to-nextto-leading order (NNLO) in QCD, it is also expected to improve substantially the sensitivity to the tt charge asymmetry [4]. Polarization observables and asymmetries in ttW ± production are capable of enhancing the sensitivity to beyond-the-SM (BSM) interactions featuring a chiral structure different from the one of the SM [4,5]. The hadro-production of ttW ± is in general well suited to directly search for BSM physics, in particular supersymmetry [6,7], supergravity [8], technicolour [9], vectorlike quarks [10], Majorana neutrinos [11] and modified Higgs sectors [12][13][14]. Beyond its own importance in LHC searches, the ttW production is a relevant background to ttH production [15].
The most recent experimental results based on Run 2 show a tension between data and theory predictions in the ttW modelling both in direct measurements [19,20] and in the context of the search for tt associated production with a Higgs boson [24,25]. While the theoretical com-munity has invested a noticeable effort to address this tension, so far no explanation emerged that is capable to fill the gap between the SM predictions and the data.
An improved modelling of the ttW ± process is required to allow for the comparison of SM predictions with future LHC data, particularly those that will be accumulated during the high-luminosity run. The increased statistics will enable not only more precise measurements of ttW ± cross-sections, but also measurements of differential distributions and in different decay channels. This target can only be achieved if the theoretical description of realistic final states embedding the ttW resonance structure is available.
Many theoretical predictions for ttW ± hadro-production are available in the literature. The first next-to-leading order (NLO) QCD calculation for TeV-scale colliders was performed in a spin-correlated narrow-width approximation for the semi-leptonic decay channel [26]. The matching of NLO QCD predictions to a parton shower was first tackled for the same decay channel in Ref. [27]. A number of calculations for ttW ± inclusive production (on-shell top/antitop quarks and W boson) have been carried out, targeting charge asymmetries [4], the impact of ttW ± on the associated production of tt pairs with a Higgs boson [15] at NLO QCD, and the effects of subleading NLO QCD and electroweak (EW) corrections [28][29][30]. Softgluon resummation up to next-to-next-to-leading logarithmic accuracy [31][32][33][34][35] and multi-jet merging [36] have also been investigated. NLO QCD corrections to the EW leading order (LO) have been computed in the narrowwidth approximation (NWA), accounting for complete spin correlations and including parton-shower effects [37]. Very recently a comparison of different fixed-order Monte Carlo generators matched to parton showers has been performed for inclusive ttW ± , with a focus on the two-charged lepton signature (in the NWA) [38]. The first studies that aim at the off-shell modelling of ttW production concern the full NLO QCD corrections in the three-charged-lepton channel [5,39,40].
The calculation of subleading NLO corrections has been performed for inclusive production [28][29][30]37], but is still missing for realistic final states. Our present work targets the complete fixed-order description of the off-shell ttW production, combining the NLO QCD and EW corrections that have a sizeable impact at the LHC@13TeV. This paper is organized as follows. In Section 2.1 we describe the process under investigation, providing details on various NLO corrections that are presented. In Sections 2.2 and 2.3 we provide the SM input parameters and the selection cuts used for numerical simulations, respectively. The integrated results at LO and NLO are presented in Section 3.1, while in Section 3.2 a number of differential distributions are described, focusing on the relative impact of various NLO corrections to the LO predictions. In Section 4 we draw our conclusions.
2 Details of the calculation

Description of the process
We consider the process which receives contributions only from quark-induced partonic channels at LO. Gluon-quark and photon-quark channels open up at NLO, while the pure gluonic channel enters only at NNLO in QCD.
Although we consider the final state with three charged leptons with different flavours, the corresponding results for the case of identical positively-charged leptons can be estimated by multiplying our results by a factor 1/2, up to interference contributions, which are expected to be small.
In this work we focus on the production of tt pairs in association with a W + boson, but the calculation of the charge-conjugate process (ttW − ) can be performed with the same techniques and no additional conceptual issues.
At LO, the largest contribution is given by the QCDmediated process of order O(α 2 s α 6 ) (labelled LO QCD ), which always embeds a gluon s-channel propagator if no quark-family mixing is assumed (diagonal quark-mixing matrix with unit entries). The tree-level EW contribution of order O(α 8 ) (labelled LO EW ), despite being characterized by many more diagram topologies, is expected to give a cross-section that is roughly 1% of the LO QCD one owing to the ratio of EW and strong coupling constants. The interference contribution, formally of order O(α s α 7 ), is identically zero due to colour algebra. In Fig. 1 we show sample diagrams for the QCD-mediated and purely-EW process. Note that the diagrams with a resonant topantitop-quark pair are present also in EW tree-level contributions.
At NLO, the ttW process receives contributions from four different perturbative orders, as depicted in Fig. 2.
The corrections that have the largest impact on the NLO cross-section are of order O(α 3 s α 6 ), which are pure QCD corrections to LO QCD . Following the notation of Refs. [29,30], we label this perturbative order as NLO 1 . These corrections have recently been computed for the full off-shell process [39,40]. This perturbative order shows a typical NLO QCD behaviour in the scale dependence, and the NLO relative corrections to LO QCD are at the 10-20% level, depending on the choice of the renormalization and factorization scale [40].
The NLO 2 corrections are known only for on-shell top-antitop quarks and for an on-shell W boson [28][29][30]. They are expected to give a negative contribution of about 4.5% of the inclusive LO cross-section.
In the off-shell calculation, as well as in the on-shell one, the NLO 2 order receives contributions not only from the EW corrections to LO QCD , but also from the QCD corrections to the LO interference, although at Born-level the O(α s α 7 ) contribution vanishes.
Sample contributions to the virtual corrections at order O(α 2 s α 7 ) are shown in Fig. 3. The diagram on the left involves one-loop amplitudes of order O(g 4 s g 6 ) interfered with tree-level EW amplitudes of order O(g 8 ) and is obviously a QCD correction to the LO interference. The diagram on the right involves one-loop amplitudes of order O(g 2 s g 8 ) interfered with tree-level QCD amplitudes of order O(g 2 s g 6 ) and could be naïvely classified as an EW correction to LO QCD . However, it can also be regarded as contributing to the QCD corrections to the LO interference. In fact, the IR singularities of this contribution are partially cancelled by the real photonic corrections to LO QCD and partially by the real gluonic corrections to the LO interference. Diagrams with weakly interacting particles in the loops are more demanding from the computational point of view, as the corresponding one-loop amplitudes include up to 10-point functions, while the latter ones feature at most 7-point functions. A selection of one-loop diagrams which contribute at this perturbative order are shown in Fig. 4.
The real-radiation contributions to NLO 2 corrections are computationally demanding due to the large multiplicity of electrically charged final-state particles. In contrast to the virtual ones, the real NLO 2 corrections can be uniquely classified into two types: the NLO EW corrections to the LO QCD process, which involve a real photon (see Fig. 5 left for an example), and the NLO QCD corrections to the LO interference, which involve a real gluon (see Fig. 5 right for an example). In the first class of contributions, the photon can be either in the final or in the initial state. The processes with a photon in the final state are characterized by many singular regions, as the photon can become soft or collinear to any of the seven charged external particles. This results in a large number of subtraction counterterms that are required to ensure a stable calculation of the NLO cross-section. The real processes with a photon in the initial state possess a smaller number of singular phase-space regions and are suppressed due to the small luminosity of photons in the proton. For on-shell production, the contribution of the photoninduced channels to the leading-order cross-section is at the sub-percent level [28]. The QCD corrections to the Fig. 1. Sample diagrams contributing to LOQCD (left) and to LOEW (right) cross-sections for off-shell ttW + production in the three-charged-lepton channel.
Contributing perturbative orders at LO and NLO for ttW hadro-production in the three-charged-lepton channel. LO interference are non-vanishing only if the radiated gluon is emitted by an initial-state light quark and absorbed by a final-state b quark or top quark (or the other way around). A sample contribution is shown in Fig. 5 right. These corrections, although necessary to account for all O(α 2 s α 7 ) contributions, turn out to be very small, as detailed in Section 3.
To sum up, the full set of real partonic channels that contribute to the NLO 2 corrections is where u and d stand for up-type and down-type quarks, respectively (of the first and second generation). The vanishing LO interference implies that the corresponding EW corrections vanish as well, since additional EW propagators (virtual contributions) and radiated photons (real contributions) do not modify the LO colour structure. Therefore, the only NLO corrections that contribute at order O(α s α 8 ) are genuine QCD corrections to the LO EW cross-section. This order is labelled as NLO 3 in Fig. 2. By simply counting the powers of α s the NLO 3 corrections are expected to give a smaller contribution than the NLO 2 ones. However, at the inclusive level [29] and in the narrow-width approximation [37], they are noticeably larger than the NLO 2 ones. This results from the fact that this perturbative order is dominated by hard real radiation diagrams in the gluon-quark partonic channel that embed the tW scattering process [1]. Sample diagrams are shown in Fig. 6. Thanks to the genuine QCD nature of the NLO 3 corrections, it is possible to match them to a QCD parton shower with no subtleties due to EW corrections, as it has been done in Refs. [37,38].
The last NLO perturbative order, O(α 9 ), is furnished by the EW corrections to the LO EW process. It has been shown at the inclusive level that such contributions are at the sub-permille level [30], as expected by naïve power counting. Even with a substantially larger data set than the one of Run 2 (i.e. 3000 fb −1 at the high-luminosity LHC) these EW effects are out of reach in a realistic Fig. 4. One-loop diagrams of order O(g 2 s g 8 ) contributing to the EW virtual corrections (NLO2) to off-shell ttW production in the three-charged-lepton channel. From left to right: sample diagrams involving 7-, 8-, 9-and 10-point functions.  Fig. 6. Sample diagrams for the partonic channel gu → bb e + νeτ + ντ µ −ν µd that contain tW + scattering as a subprocess and contribute to NLO3 corrections to ttW production in the three-charged-lepton channel.
fiducial region. Therefore, we are not providing results for this order.
In the following we focus on the first three NLO perturbative orders. Tree-level and one-loop SM amplitudes are computed with the Recola matrix-element provider [41,42]. For the tensor reduction and evaluation of loop integrals we use the Collier library [43]. The multichannel Monte Carlo integration is performed with Mo-CaNLO, a generator that has already been used to compute the NLO QCD corrections to ttW [40] and the NLO EW corrections to several LHC processes involving top quarks [44,45]. The subtraction of infrared and collinear singularities is carried out using the dipole formalism of Refs. [46][47][48] both for QCD and for EW corrections. The initial-state collinear singularities are absorbed in the parton distribution functions (PDFs) in the MS factorization scheme.

Input parameters
We consider proton-proton collisions at a centre-of-mass energy of 13 TeV. We neglect flavour mixing in the quark sector and use a unit quark-mixing matrix. The three charged leptons that we consider in the final state are massless and characterized by three different flavours. The on-shell masses and widths of weak bosons are set to the following values [49], and then translated into their pole values [50] that enter the Monte Carlo simulations. The Higgs-boson mass and width are fixed, following Ref. [49], to We have computed the LO top-quark width according to Ref. [51], using the pole values for the W-boson mass and width. The NLO top-quark width is obtained applying the NLO QCD and EW correction factors of Ref. [52] to the LO width. The numerical values read m t = 173.0 GeV , The top-quark width is kept fixed when performing variations of the factorization and renormalization scales. The EW coupling is treated in the G µ scheme [53], where G µ = 1.16638 · 10 −5 GeV −2 is the Fermi constant. The masses of weak bosons and of the top quark, and therefore also the EW mixing angle, are treated in the complex-mass scheme [53][54][55][56].
Both for LO and NLO predictions, we employ NNPDF 3.1 luxQED PDFs [57], computed at NLO with α s = 0.118. Using these PDFs, the photon contribution is properly accounted for in the evolution. The evaluation of PDFs and the running of the strong coupling are achieved with the LHAPDF6 interface [58].

Selection cuts
Coloured partons with |η| < 5 are clustered into jets by means of the anti-k t algorithm [59] with resolution radius R = 0.4. The same algorithm but with R = 0.1 is applied to cluster photons into charged particles.
We choose cuts that mimic the fiducial selections applied by ATLAS in Ref. [20] and that have already been used for the study of NLO QCD corrections to the same final state [40]. We select events with exactly two b jets (assuming perfect b-tagging efficiency) that are required to satisfy Furthermore, we ask for three charged leptons that fulfil standard acceptance and isolation cuts, where the R distance is defined as the sum in quadrature of the azimuthal and rapidity separations, We do not constrain the missing transverse momentum and do not apply any veto to additional light jets.

Results
For the factorization and renormalization scale (µ F = µ R ), we consider three different dynamical choices that have proved to behave better than a fixed scale [39,40]. The first one, introduced in Ref. [39], depends on the transverse-momentum content of the final-state particles, regardless of the top-antitop resonances, 1 1 We use the same notation for the scales as in Ref. [40].
The second and third choices, already used in Ref. [40], are based on the transverse masses of the top and antitop quarks. Due to the ambiguity in choosing the + ν pair that results from the top quark, we pick the pair of leptons that when combined with the bottom quark forms an invariant mass closest to the top-quark mass. We consider two choices that differ by a factor of 1/2, (10) In the following we present the results combining NLO corrections with the so-called additive approach, This approach is exact at the order of truncation of the perturbative expansion. Furthermore, it represents a natural choice for our process, as the combination involves NLO corrections to two different leading orders that do not interfere.

Integrated cross-sections
In Table 1, we present the integrated cross-sections in the fiducial region defined in Section 2.2.
The leading corrections to the LO QCD cross-section are expected to come from the corresponding pure QCD radiative corrections (NLO 1 ). Their inclusion has been proved to decrease the theoretical uncertainty due to scale variations and to stabilize the perturbative convergence for this process. Nonetheless, their impact depends on the choice of the scale, as already pointed out both in inclusive [28,29] and off-shell [39,40] computations. In fact, with the resonance-blind dynamical choice µ (c) 0 , the NLO 1 corrections give a 6.6% enhancement to the LO QCD cross-section, with the resonance-aware choices, µ (d) 0 and µ (e) 0 , they give a 18% and a 0.4% correction, respectively. Note that, at variance with Ref. [40] in this paper we compute both LO and NLO predictions with the NLO top-quark width, which gives roughly a 12% enhancement to the LO cross-section.
The NLO 2 corrections are negative and amount to about −5.5% of the LO QCD cross-section for all scale choices. Such supposedly subleading corrections have a sizeable impact on the fiducial NLO cross-section, and this is likely due to large EW Sudakov logarithms enhancing the cross-section in the high-energy regime [60]. This is supported by the fact that the average partonic centre-of-mass energy and H T are quite high, 850 GeV and 520 GeV, respectively. A crude estimate of the Sudakov logarithms gives a result which is of the same order of magnitude as the full NLO 2 corrections we have obtained for this process.
The impact of QCD corrections that can be uniquely attributed to the LO interference is very small, both for real and for virtual corrections, accounting respectively for 5% and for less than 1% of the total NLO 2 result. At O(α 2 s α 7 ), we have also included the contribution from photon-initiated partonic channels, which are positive and account for about 0.1% of the LO QCD crosssection. As already observed at the inclusive level [28], this contribution is very small and its effect will be hardly visible even at the high-luminosity run of LHC (they will yield about 1 event for √ s = 13 TeV and L = 3000 fb −1 ). In inclusive calculations with on-shell top-antitop quarks, the NLO 2 corrections were found to give a −4.5% contribution to the inclusive production cross-section [28,30]. In order to compare our results with those of Ref. [30], we have performed a full off-shell calculation in a very inclusive setup, and divided by the branching ratios for the decays of the top and antitop quarks and of the W boson. The setup is the same as the one of Ref. [30], up to a few unavoidable differences: we use finite top-quark and W widths (same values as those of Section 2.2) and we include Higgs-boson contributions; we apply a minimum invariant-mass cut of 5 GeV to the bb system to protect from infrared singularities and we cluster photons into charged particles with isolation radius R = 0.1; we employ the same dynamical scale as in Ref. [30], but using the kinematics after photon recombination and choosing the top-quark candidate with the same invariant-mass prescription as for the calculation of the scales µ exhibit NLO 2 corrections of −3% of LO QCD , which is not far from the −4.5% of Refs. [28,30]. The remaining discrepancy should be due to both the additional cuts we have applied and the non-resonant effects that are included in our full calculation, while being absent in onshell calculations. The comparison of results in the fiducial and the inclusive setup reveals that the NLO 2 corrections are more sizeable for realistic final states and in the presence of reasonable fiducial cuts.
Coming back to our default fiducial setup (see Section 2.2), the NLO 3 perturbative order is dominated by the real radiation in the ug partonic channel, owing to a PDF enhancement and the tW scattering embedded in this channel. These real corrections account for 85% of the NLO 3 corrections, and are one order of magnitude larger than the corresponding leading order LO EW . The total NLO 3 corrections amount to 13% of the LO QCD crosssection, almost independently of the scale choice. This confirms the 12% effect obtained in the case of on-shell top and antitop quarks [30].
It is worth stressing that the inclusion of NLO 2 and NLO 3 corrections gives a noticeable effect to the ttW cross-section. Therefore such corrections must definitely be accounted for in experimental analyses. Furthermore, their relative contribution to the LO result is rather independent of the scale choice, while the NLO 1 corrections are much more scale dependent.
As a last comment of this section, we point out that the scale uncertainty of the combined LO+NLO crosssection is driven by the NLO 1 corrections which reduce the LO QCD uncertainty roughly from 20% to 5%. Due to their EW nature, the NLO 2 corrections do not diminish the scale uncertainty of the corresponding leading order (LO QCD ). The NLO 3 corrections exceed the corresponding pure EW LO process and, thus, imply a LO-like scale dependence for LO EW +NLO 3 .
So far, we have focused on the relative contributions of various NLO corrections to the fiducial ttW cross-section. However, the interplay among different corrections can be rather different for more exclusive observables. Therefore, it is essential to study differential distributions.

Distributions
In the following we present a number of relevant distributions focusing on the impact of the various NLO corrections relative to the LO QCD cross-sections. Since in most of the LHC experimental analyses the theoretical predictions are NLO QCD accurate, we also comment on the distortion of NLO QCD distribution shapes (LO QCD + NLO 1 in our notation) due to the inclusion of NLO 2 and NLO 3 corrections. We choose to present all differential distributions with the µ (e) 0 scale. 2 We start by presenting transverse-momentum distributions in Figs. 7-8.
In Fig. 7(a) we consider the distribution in the transverse momentum of the positron, which is precisely measurable at the LHC. Since we also include EW corrections (NLO 2 ), the positron is understood as dressed (a radiated photon could be clustered into the positron). The distribution peaks around 50 GeV, where the relative impact of QCD and EW corrections follows straight the integrated results. Relative to the LO QCD , all three radiative corrections drop in a monotonic manner. Nonetheless, the decrease of NLO 3 corrections is very mild (14% below 50 GeV, 10% at 380 GeV), while the NLO 1 and NLO 2 corrections decrease steeper: the former become negative around 80 GeV and give −9% at 380 GeV, the latter become lower than −10% already at moderate p T,e + (200 GeV). The behaviour of EW corrections in the tails of this distribution is likely driven by the impact of Sudakov logarithms, which become large at high p T .
The same behaviour of the NLO 2 and NLO 3 corrections characterizes also the distribution in the transverse momentum of the antitop quark, shown in Fig. 7(b). The antitop momentum is computed as the sum of the momenta of the muon, its corresponding antineutrino, and the antibottom quark. This is not observable at the LHC, but its analysis is useful to compare the full off-shell calculation with the on-shell ones. The negative growth of the NLO 2 corrections behaves very similarly in the inclusive calculations, as can be seen for example in Ref. [28] (figure 5 therein). This confirms that for sufficiently inclusive variables the NLO EW effects are dominated by contributions with resonant top and antitop quarks. The NLO 1 corrections increase by roughly 25% in the considered spectrum.
In both transverse-momentum distributions of Fig. 7, the inclusion of subleading NLO corrections (NLO 2 , NLO 3 ) gives a decreasing effect towards large transverse momenta to the NLO QCD cross-section. In fact, the ratio between the combined LO + NLO crosssection and the LO QCD + NLO 1 one ranges between 1.10 and 1.15 for small transverse momenta and drops below 1 already at moderate transverse momenta. We have checked numerically that this conclusion can be drawn 2 In Ref. [40] the scale for the differential distributions is In Fig. 8 we consider the distributions in the transverse momentum of the tt and the bb system. The former is not measurable at the LHC but can be reconstructed from Monte Carlo truth, choosing the positively-charged lepton-neutrino pair that best reconstructs the top-quark mass when combined with the momentum of the bottom quark. The latter variable is directly observable at the LHC.
The NLO 1 corrections to the tt transverse-momentum distribution [ Fig. 8(a)] have already been investigated in Ref. [40]: these QCD corrections grow monotonically and become dramatically large and positive at high p T,tt . Note that at LO this variable coincides with the p T of the recoiling W + boson, and therefore is sensitive to the real QCD radiation which is not clustered into b jets (and that cannot be clustered to the W + -boson decay products). The NLO 1 corrections receive a sizeable contribution by the gq/gq partonic channels, which are enhanced by the gluon PDF. In contrast, the NLO 2 ones feature a typical NLO EW behaviour in the tail of the distribution, giving a negative and monotonically decreasing correction to the LO cross-section. Differently from the QCD corrections, the additional photon can be clustered into any of the external charged particles, thus also into the decay products of the recoiling W boson. Furthermore, the cross-section is not enhanced by the γq/γq partonic channels due to the very small photon luminosity in the proton. The NLO 3 contribution is positive in the whole analyzed spectrum, and increases from 6% (below 50 GeV) to a maximum of 25% around p T,tt = 2m t , then it slowly decreases in the large-p T region. Relative to the LO QCD + NLO 1 result, the combination of all other corrections gives an effect which is about 10% in the soft-p T region and diminishes towards negative values at large p T .
The transverse momentum of the bb system [ Fig. 8(b)] is correlated to the one of tt system. The NLO 2 corrections to this observable behave in the same manner as those for the p T,tt distributions, giving a −20% contribution around 400 GeV. The NLO 3 corrections grow monotonically from +10% (at low transverse momentum) to +30% (around 400 GeV). In the soft part of the spectrum (p T,bb < 150 GeV), the NLO 1 corrections are rather flat, while in the large-p T region they grow positive and become very large, similarly to the p T,tt distribution. The overall NLO corrections are very small below 150 GeV due to mutual cancellations among the three contributions, while at larger transverse momenta the corrections are dominated by the NLO 1 contribution. Furthermore, relatively to the LO QCD + NLO 1 distribution, the combined NLO corrections give a flat and positive effect between 7% and 8% in the whole analyzed spectrum.
In all analyzed transverse-momentum distributions of Figures 7-8, the LO EW contribution increases monotonically (relatively to the LO QCD one) but never exceeds 3% of the LO QCD cross-section.
In Fig. 9(a), we display the distribution in the invariant mass of the antitop quark, which in our setup can be reconstructed from the Monte Carlo truth. The lineshape   is dominated by the Breit-Wigner distribution of the leptonic decay of the antitop quark. The NLO 1 corrections are negative at the peak while below the pole mass they give a very large enhancement to the LO result. Such a radiative tail, coming from unclustered real radiation, is also present, though less sizeable, in the NLO EW corrections (unclustered photons). At values larger than the top-quark mass, the distribution receives an increasingly positive contribution from NLO 1 corrections, while the NLO 2 ones give an almost flat correction of −10% to the LO QCD cross-section. The LO EW contribution shows a slightly wider distribution than the LO QCD one. This could be attributed to the relatively larger contribution of non-resonant background diagrams in the LO EW contribution. Nonetheless, the impact of this difference on the full distribution is almost invisible owing to the very small size of the LO EW contribution. The NLO 3 corrections behave differently from the NLO 1 ones, giving a rather flat enhancement of the fiducial LO cross-section, which is minimal around the peak (+11% at m t ) and mildly increases towards the tails (+20% at 200 GeV, +30% at 150 GeV). This is due to the very large contribution of the u(c)g partonic channel, which has a light d(s) in the final state that cannot come from the radiative decay of the top or of the antitop quark (differently from final-state gluons). As can be seen in the bottom panel of Fig. 9(a), the total LO + NLO result is 12% higher than the LO QCD + NLO 1 one below the top-quark mass. This enhancement is smaller at (4%) and above (7%) the top mass.
Another variable that is often investigated in LHC analyses is H T , whose definition is given in Eq. (9). The LO and NLO distributions in this observable are shown in Fig. 9(b). As in the transverse-momentum distributions studied above, the NLO 2 radiative corrections decrease monotonically towards large values of H T (about −20% for H T ≈ 1.2 TeV). The LO EW contribution grows to 5% of the LO QCD cross-section at 1.5 TeV, where the NLO cross-section is two orders of magnitude lower than its value at the maximum of the distribution. The NLO 3 corrections are rather flat and enhance the LO QCD result between 10% and 15%. The NLO 1 corrections are characterized by a non-flat shape that is increasing for H T < 800 GeV from −10% to +25% and decreasing in the rest of the considered spectrum. We further observe that the combination of the three NLO perturbative orders yields an almost vanishing correction in the soft region of the spectrum, while in the tail of the distribution the overall correction is dominated by the NLO 2 contribution for our scale choice. In a similar fashion as in other   transverse-momentum distributions, the ratio of the combined LO+NLO result over the LO QCD +NLO 1 decreases monotonically from 1.15 to 0.95 in the analyzed range.
In Fig. 10 we study more invariant-mass distributions. The distribution in the invariant mass of the two-b-jet system [ Fig. 10(a)] is characterized by rather flat QCD corrections (NLO 1 and NLO 3 ). The NLO 3 corrections enhance the LO QCD cross-section by 11% to 14% everywhere in the analyzed invariant-mass range. The NLO 2 contribution has a similar behaviour as the one found for the previous variables, growing negative towards the tail of the distribution.
The distribution in the invariant mass of the threecharged-lepton system is considered in Fig. 10(b). The behaviour of the NLO 2 and NLO 3 corrections follows closely the one for the bb system, apart from a less steep decrease of the EW corrections towards large invariant masses. These corrections are at the −10% level for masses larger than 500 GeV. The NLO 1 corrections vary by hardly more than 10% in the studied range.
As shown in the bottom panels of Fig. 10, both for the bb system and for the three-charged-lepton system, the inclusion of NLO 2 and NLO 3 corrections (as well as of LO EW , though hardly visible) gives a non-flat correction to the NLO QCD invariant-mass distributions, decreasing monotonically from +12% to zero in the considered spectra.
After presenting transverse-momentum and invariantmass distributions, we switch to some relevant angular variables. In Fig. 11(a) and Fig. 11(b) we display the distributions in the rapidity of the muon and of the antitop quark, respectively. Since these two variables are correlated (the dominant resonant structure involves the decaȳ t →bµ −ν µ ), the muon rapidity, which is precisely measurable at the LHC, represents a suitable proxy for the rapidity of the corresponding antitop quark (which can only be reconstructed from Monte Carlo truth). Note that the muon rapidity is sharply cut at ±2.5 by fiducial selections, while this is not the case for the antitop quark. However, thanks to the rapidity cut applied to b jets, the cross-section is strongly suppressed for |yt| > 2.5. Both the muon and the antitop quark are produced preferably in the central region. The NLO 2 corrections are rather flat, giving between −4% and −8% decrease to the LO cross-section. The relative NLO 1 corrections to the muonrapidity distribution are characterized by a large variation (about 35%) in the available range. Relative to LO QCD , the differential NLO 3 corrections have a similar shape as the NLO 1 ones, giving in the whole rapidity range a positive correction (8% in the forward regions, 16% in the central region). Almost identical results are found in the rapidity distribution of the antitop quark. Owing to the NLO 3 corrections, the ratio between the complete NLO prediction and the LO QCD +NLO 1 one has   a maximum of 1.11 in the central region and diminishes towards forward regions (close to unity). This holds both for the muon and for the antitop-quark rapidity spectra. In Fig. 12(a) we consider the distribution in the azimuthal separation between the positron and the muon. The two charged leptons tend to be produced in opposite directions both at LO and at NLO, but the inclusion of radiative corrections enhances the fraction of events with small azimuthal separations. The NLO 2 corrections are negative and roughly constant (−5% to −7%) over the full angular range, while the NLO 3 contribution decreases monotonically from +18% to +11% relatively to the LO QCD result. As already observed in Ref. [40], the NLO 1 correction to the LO QCD cross-section decreases with an almost constant negative slope over the full range. The overall NLO corrections to the LO QCD cross-section are positive everywhere except in the vicinity of the peak at ∆φ e + µ − = π. Relative to the LO QCD + NLO 1 prediction, the combination of NLO 2 and NLO 3 corrections gives a pretty flat enhancement (1.11 at ∆φ e + µ − = 0, 1.06 at ∆φ e + µ − = π).
As a last differential result, we present in Fig. 12(b) the distribution in the R distance between the two b jets [see Eq. (8) for its definition]. This distribution is characterized by an absolute maximum around ∆R bb ≈ π. The negative NLO 2 corrections diminish monotonically over the analyzed spectrum. At large distance (∆R bb > 5) they give a contribution of −15%. The positive NLO 3 corrections diminish from +15% at ∆R bb ≈ 0 to +12% at ∆R bb ≈ π and then increase again. The NLO 1 ones show a similar behaviour, however, with larger slopes. The combined NLO 2 and NLO 3 corrections enhance the LO QCD +NLO 1 prediction between 6% and 11%, similarly to the case of the azimuthal distance shown in Fig. 12(a), but with a somewhat different shape.

Conclusions
In this work we have presented the NLO corrections to the off-shell production of ttW + at the LHC in the three-charged-lepton channel. These include the nextto-leading-order (NLO) QCD corrections to the QCD (NLO 1 ) and to the electroweak leading order (NLO 3 ), as well as the NLO electroweak corrections to the QCD leading order (NLO 2 ). It is the first time that the NLO 2 and NLO 3 radiative corrections are computed with full off-shell dependence for a physical final state, accounting for all non-resonant, interference, and spin-correlation effects.
Both integrated and differential cross-sections have been presented and discussed in a realistic fiducial region, keeping in mind the limited statistics of the LHC data and relating the off-shell description of the process to the inclusive predictions that are available in the literature.   The NLO 2 and NLO 3 corrections give a −5.5% and a +13% contribution, respectively, to the LO cross-section, almost independently of the choice of the factorization and renormalization scales. The sizeable impact of NLO 2 and NLO 3 corrections makes it essential to combine them with the NLO 1 ones, in order to arrive at reliable predictions.
The theory uncertainties from 7-point scale variations are driven by the NLO 1 corrections, which are the only corrections which feature a NLO-like scale dependence. Their inclusion reduces the scale uncertainty of the LO cross-section from 20% to 5%. The investigation of differential distributions reveals a more involved interplay among the various perturbative orders compared with the integrated results. The NLO 2 corrections drop by up to −20% in most of the transversemomentum and invariant-mass distributions, showing the typical behaviour of EW corrections with large Sudakov logarithms at high energies. They are rather flat for angular observables. The NLO 3 corrections give a positive enhancement between +10% and +20% (30% in some cases) to the LO cross-section in all analyzed distributions, They are dominated by the ug partonic channel (formally belonging to QCD real corrections to LO EW) that embeds tW scattering. The NLO 1 corrections, which have already been presented in the literature, show quite variable patterns in the various differential K-factors.
We stress that all three NLO contributions usually give non-flat corrections to the LO distributions, also to the angular ones. This indicates that rescaling QCD results (either LO or NLO accurate) by flat K-factors could result in a bad description of some LHC observables.
In the light of an improved experimental description of the ttW process, the inclusion of decay and off-shell effects is mandatory. Although for sufficiently inclusive observables the full computation is well approximated by on-shell calculations, the inclusion of off-shell effects in the modeling of ttW production is definitely needed when studying the tails of transverse-momentum and invariantmass observables.