Measurement of the inclusive and differential t t ¯ $ $ \overline{t} $$ # cross sections in the single-lepton channel and EFT interpretation at s $$ \sqrt{s} $$ = 13 TeV

: The production cross section of a top quark pair in association with a photon is measured in proton-proton collisions at a center-of-mass energy of 13 TeV. The data set, corresponding to an integrated luminosity of 137 fb − 1 , was recorded by the CMS experiment during the 2016–2018 data taking of the LHC. The measurements are performed in a ﬁducial volume deﬁned at the particle level. Events with an isolated, highly energetic lepton, at least three jets from the hadronization of quarks, among which at least one is b tagged, and one isolated photon are selected. The inclusive ﬁducial t¯t γ cross section, for a photon with transverse momentum greater than 20 GeV and pseudorapidity | η | < 1 . 4442 , is measured to be 798 ± 7 (stat) ± 48 (syst) fb, in good agreement with the prediction from the standard model at next-to-leading order in quantum chromodynamics. The diﬀerential cross sections are also measured as a function of several kinematic observables and interpreted in the framework of the standard model eﬀective ﬁeld theory (EFT), leading to the most stringent direct limits to date on anomalous electromagnetic dipole moment interactions of the top quark and the photon.


Introduction
The large amount of proton-proton (pp) collision data at a center-of-mass energy of 13 TeV at the LHC allows for precision measurements of standard model (SM) processes with small production rates. Among these, top quark production provides a testing ground for the SM predictions and for phenomena beyond the SM (BSM). In particular, precise measurements of the inclusive and differential cross sections of top quark pair production in association with a high-energy photon (ttγ) constrain anomalous ttγ electroweak interactions [1][2][3][4].
The CDF Collaboration at the Fermilab Tevatron measured the ttγ production cross section using proton-antiproton collisions at √ s = 1.96 TeV [5], while at the LHC the Collaborations. At 13 TeV, the ATLAS Collaboration measured inclusive and differential ttγ production cross sections in leptonic [9] and in the eµ [10] final states. All of these results are in agreement with the SM.
In this paper, the inclusive and differential ttγ production cross sections are measured in pp collisions at √ s = 13 TeV. The analysis uses a data sample recorded with the CMS detector during Run 2 (2016-2018) of the LHC, which corresponds to an integrated luminosity of 137 fb −1 . The measurement is performed in the single-lepton (electron or muon) final state in a fiducial region defined at particle level. The inclusive fiducial ttγ cross section is measured for a selection on the photon transverse momentum of p T (γ) > 20 GeV and the pseudorapidity of |η(γ)| < 1.4442, corresponding to the barrel region of the CMS electromagnetic calorimeter (ECAL). Differential cross sections are measured in the same fiducial region as a function of p T (γ), |η(γ)|, and the angular separation between the lepton and the photon, ∆R( , γ). The observations are interpreted in the context of the SM effective field theory (SM-EFT) [11], where the c tZ and c I tZ operators, defined in ref. [12], are constrained using the measurement of the distribution of p T (γ). Tabulated results are provided in HEPDATA [13]. Examples of Feynman diagrams at leading order (LO) contributing to the ttγ signal topology are shown in figure 1.
This paper is organized as follows. The CMS detector is briefly introduced in section 2. Details on the simulation of the signal and background processes and their modeling are provided in section 3. The online selection, event reconstruction, and object definitions are described in section 4. The fiducial phase space definition and photon categorization are described in section 5. The event selection and the statistical treatment are discussed in section 6. The procedures to estimate the backgrounds are described in section 7 and the systematic uncertainties are discussed in section 8. The obtained results and the interpretation of the measurements in the context of SM-EFT are presented in section 9. Finally, a summary is provided in section 10. candidates, p T -and η-dependent thresholds are set on the pileup corrected scalar p T sum of photons and neutral and charged hadrons reconstructed by the PF algorithm (I rel (e)) in the range of 5-10%. The chosen "tight" electron working point has a 70% efficiency while rejecting electron candidates originating from jets [58]. A muon candidate is isolated if it satisfies I rel (µ) < 0. 15. The efficiency of the chosen "tight" working point is 90-95%, depending on p T and η of the muon candidate [55]. For photon candidates, the scalar p T sum of the charged particles within a cone of ∆R = 0.3, denoted as the photon chargedhadron isolation, must satisfy I chg (γ) ≤ 1.141 GeV. Depending on the photon candidate p T , there are separate requirements on the photon neutral-hadron and total isolation [58]. The photon reconstruction and selection efficiency for the chosen "medium" working point in simulation is on average 80%. The electron, muon, and photon reconstruction efficiencies are corrected as a function of the p T and η of the reconstructed object to match the efficiency observed in data.
Furthermore, "loose" selection criteria are used to define control regions and to veto events with additional reconstructed leptons and photons. With respect to the tight electron selection, the transverse momentum requirement is relaxed to p T > 15 GeV, the threshold on I rel (e) to the range of 20-25%, depending on p T and η of the electron candidate, and two (three) missed hits in the innermost tracker layers are allowed for electrons in the barrel (endcap) region. The loose muon selection is based on ref.
Jets are reconstructed by clustering PF candidates using the anti-k T algorithm [60, 61] with a distance parameter of 0.4. Selected jets are required to satisfy p T > 30 GeV and |η| < 2.4. Contributions to the clustered energy from pileup interactions are corrected for by requiring charged-hadron candidates to be associated with the PV and an offset correction for the contribution from neutral hadrons falling within the jet area is subtracted from the jet energy. Corrections to the jet energy scale (JES) are applied in simulation and data. The jet energy resolution (JER) is corrected in simulation to match the resolution observed in data [62].
Jets originating from the hadronization of b quarks are identified (b tagged) with a deep neural network algorithm [63] based on tracking and secondary vertex information. A working point is chosen such that the efficiency to identify the b jet is 55-70% for a jet p T of 20-400 GeV. The misidentification rate in this p T range is 1-2% for light-flavor and gluon jets, and up to 12% for charm quark jets. A correction is applied to the simulation to match the b tagging efficiencies observed in data.
The missing transverse momentum vector, p miss T , is defined as the projection onto the plane perpendicular to the beams of the negative vector momentum sum of all PF candidates in an event. The JES and JER corrections are included in the p miss T computation. Its magnitude is referred to as p miss T .

JHEP12(2021)180
Photon no hadronic origin no hadronic origin ∆R(jet, ) > 0.4 Table 2. Overview of the definition of fiducial regions for various objects at particle level. A photon is isolated, if there are no stable particles (except neutrinos) with p T > 5 GeV within a cone of ∆R = 0.1.

Fiducial phase space definition and photon classification
The fiducial region of the analysis is defined at the particle level by applying an event selection to the stable particles after the event generation, parton showering, and hadronization, but before the detector simulation. Electrons (muons) must have p T > 35 (30) GeV and |η| < 2.4, and must not originate from hadron decays. To account for final-state photon radiation, the four-momenta of photons inside a cone of ∆R = 0.1 are added to the lepton before the lepton selection [64]. Events with leptonically decaying τ leptons in the decay chain of the top quark are considered signal.
Photons are selected if they do not originate from hadron decays, satisfy p T (γ) > 20 GeV and |η(γ)| < 1.4442, and are found outside a cone of ∆R = 0.4 around the leptons. An isolation requirement is applied by removing photons with stable particles (except neutrinos) found within a cone of ∆R = 0.1 that satisfy p T > 5 GeV.
Particle-level jets are clustered using the anti-k T algorithm with a distance parameter of 0.4, using all final-state particles, excluding neutrinos. Jets must satisfy p T > 30 GeV and |η| < 2.4. A ghost matching method [65] is used to determine the flavor of the jets, with those matched to b hadrons tagged as b jets. Finally, the overlap of jets and other candidates is removed by excluding jets with ∆R ≤ 0.4 (0.1) to lepton (photon) candidates. A summary of the object definitions at particle level is provided in table 2.
The fiducial region is constructed by requiring exactly one photon, exactly one lepton, and three or more jets among which at least one must be b tagged. The inclusive fiducial cross section, predicted with MadGraph5_amc@nlo at NLO in QCD, is 773 ± 135 fb. The NLO effects in the decay of the top quarks are not included in this calculation.
To facilitate the estimation of backgrounds with nonprompt and misidentified photons, a photon categorization is based on the matching between the reconstructed photon and simulated particles. Reconstructed photons are matched in ∆R to the corresponding generator-level particle from the primary interaction. The maximum ∆R considered for matching is 0.3 and the p T (γ) is required to be within 50% of the matched particle. Simulated events with a reconstructed photon are subsequently classified into three categories based on the matched generator particle. In the "genuine photon" category, the -7 -JHEP12(2021)180 reconstructed photon is matched to a generated photon that originates from a lepton, a W boson, or a quark. In the "misidentified electron" category, the photon is matched to an electron. The "nonprompt photon" category is comprised of events where the photon is matched to a generated photon that originates from a hadron (71%), or in absence of a match to a generated photon or electron. This category thus includes contributions with misidentified photons and photons that originate from pileup interactions (29%). 6 Analysis strategy

Signal and control region definitions
The ttγ process typically produces events with several jets, up to two b-tagged jets, and an isolated photon with large p T . The measurement is performed in signal regions with exactly one lepton (N = 1), exactly one photon (N γ = 1), and at least three jets (N j ≥ 3), among which at least one is b tagged (N b ≥ 1). Events with additional leptons or photons passing the loose selection are removed. The measurement is performed in the N j = 3 and ≥4 signal selections, denoted by SR3 and SR4p, respectively. Signal events with a jet failing the identification criteria thus enter the SR3 region. The N j ≥ 3 selection is denoted by SR3p. For illustration, figure 2 shows some kinematic distributions in the SR3p region where the simulated event samples are categorized according to the origin of the photon. The backgrounds are normalized according to the methods described in section 7 and the pre-fit systematic uncertainties are shown as a hatched band. In this figure, M 3 denotes the invariant mass of the three-jet combination among all identified jets that maximizes the magnitude of the vector p T sum [66]. This choice preferentially captures the hadronic top quark decay products.
The data-based estimation procedures for the dominant background sources are described in section 7. The simulation predicts a significant background contribution from nonprompt photons (23%), misidentified electrons (19%), and a small contribution from multijet events in the SR3p region. The nonprompt photon contribution is estimated using background-enriched control regions with relaxed criteria on I chg (γ) and σ ηη (γ). The multijet contributions to the signal and control regions are estimated by rescaling suitable normalized distributions (templates) obtained from background-enriched high-I rel ( ) sidebands. The misidentified electron background is estimated in a N b = 0 region where the invariant mass of the electron and photon candidates (m(e, γ)) is consistent with the Z boson hypothesis [67] within 10 GeV, i.e., |m(e, γ) − m Z | ≤ 10 GeV, where m Z is the Z boson mass. The control region is denoted by misDY3 (misDY4p) for N j = 3 (≥4). The Wγ and the Zγ processes contribute events with genuine photons to both the signal regions and the misDY3 and misDY4p control regions. In the electron channel, their contribution is constrained in "low mass" (LM) and "high mass" (HM) regions, defined by m(e, γ) < m Z − 10 GeV and m(e, γ) > m Z + 10 GeV, respectively. In the muon channel, the LM (HM) region is defined by m(µ, γ) < m Z (m(µ, γ) > m Z ), where m(µ, γ) denotes the invariant mass of the muon and the photon. Table 3 provides a summary of the kinematic requirements in the signal and control regions.

Statistical treatment
The signal cross section is extracted from signal and control regions using the statistical procedure detailed in refs. [68,69]. The observed yields, signal and background estimates in each analysis category, and the systematic uncertainties are used to construct a binned likelihood function L(r, θ) as the product of Poisson probabilities of all bins. The nuisances related to the systematic uncertainties in the experiment and in the modeling of signal and background processes are described by log-normal probability density functions. The parameter r is the signal strength modifier, i.e., the ratio between the measured cross section and an arbitrary reference value of 773 fb, chosen as the nominal prediction for the inclusive fiducial cross section. The symbol θ represents the set of nuisance parameters describing the systematic uncertainties. The number of reconstructed ttγ signal events JHEP12(2021)180 generated outside the fiducial phase space is scaled with the same value of r, i.e., no independent production cross section is assumed for this part of the signal.
The used test statistic is the profile likelihood ratio, q(r) = −2 ln L(r,θ r )/L(r,θ), wherê θ r reflects the values of the nuisance parameters that maximize the likelihood function for a signal strength modifier r. The quantitiesr andθ are the values that simultaneously maximize L. A multi-dimensional fit is used to extract the observed cross section of the signal process, the nuisance parameters, and the uncertainties in the nuisance parameters [68, 69].
The LM3, LM4p, HM3, HM4p, misDY3 and misDY4p control regions enter the likelihood fit separately for each data-taking period and lepton flavor. In order to extract the p T (γ) dependence of the background with misidentified electrons, the misDY3 and misDY4p control regions are split into 7 bins separated by the p T (γ) thresholds 20, 35, 50, 65, 80, 120, and 160 GeV. The LM3, LM4p, HM3, and HM4p regions are similarly separated into 3 bins defined by the p T (γ) thresholds 20, 65, and 160 GeV. The binning is chosen to obtain a statistical uncertainty in simulated background yields of less than 15%.
The likelihood fit is performed for the inclusive cross section measurements, and separately for the differential measurement. For the extraction of the inclusive cross section, the SR3 and SR4p signal region are divided in three M 3 bins in the ranges 0-280, 280-420, and >420 GeV. The binning in p T (γ), |η(γ)|, and ∆R( , γ) in the SR3 and SR4p signal regions for the differential measurements is provided in section 9.2. The estimation of contributions from various background processes is performed using control regions binned in p T (γ), which are used in all inclusive and differential measurements.

Multijet background
The probability for a multijet event to mimic the final state of the signal process is small and subject to large uncertainties. Therefore, the background from multijet events, comprising events with misidentified and nonprompt leptons, is estimated with a data-based procedure in sideband regions with loosened isolation criteria. For each N j requirement, a sideband region is defined by N b = 0 and requiring the lepton to pass the isolation criterion of the loose lepton working point and to fail the tight lepton selection. The N γ = 1 requirement is kept. The resulting selection is dominated by multijet events. After electroweak backgrounds and backgrounds with top quarks (W+jets, t/tt, and Drell-Yan) are subtracted based on the expectation from simulation, templates for the distributions of kinematic observables are extracted.
The template normalization is evaluated from a transfer factor ("TF"), defined as the ratio of the multijet event yield with tightly isolated lepton candidates to the yield with loosely isolated lepton candidates. It is obtained in a selection with N j = 2 and N γ = 0 by fitting the distribution of the transverse mass of the W boson candidate, calculated from the formula where indicates the lepton considered in the event. The distribution is taken from data in a N b = 0 selection with loosely isolated leptons, and electroweak and top quark backgrounds are subtracted. The fit is then performed in the selection with tightly isolated leptons where the total normalization of the electroweak and top quark background is left floating, while its shape is again taken from simulation. For illustration, the fit result for the N j = 2 and N b = 0 region, including the m T (W) multijet distribution from the selection with loosely isolated leptons, is shown in figure 3.
Because the efficiency of the tight lepton selection in multijet events depends on p T and η of the lepton, the estimation procedure, including the TF fit, is performed in a total of 24 bins defined in these observables. Depending on p T and η of the lepton, the TFs vary in the range of 0.9-3.1 (0.1-0.3) for the e channel and 2.0-3.7 (0.6-1.0) for the µ channel, for N b = 0 (≥1). A correction based on simulated multijet events accounts for the TF dependence on N j . Finally, the multijet estimate is obtained by scaling the N b = 0 sideband templates with the corresponding TFs and accumulating the resulting predictions in the 24 bins in lepton p T and η. The total multijet yield is estimated at 12 (8)% in the e (µ) channel in the LM3p, HM3p and misDY3p control regions and below 0.5% in the signal regions.

Nonprompt photon background
The nonprompt photon background component is estimated from data by exploiting the difference between its distribution in the plane defined by the weakly correlated variables σ ηη (γ) with I chg (γ), and the corresponding distribution for genuine photons. In a sideband with a requirement of σ ηη (γ) ≥ 0.011 on the photon candidate, the expected yields with genuine photons, misidentified electrons, and multijet events are subtracted. The sideband is used to obtain the normalization factor r SB , defined as the ratio of the yield passing the I chg (γ) < 1.141 GeV requirement to the event yield failing it. The estimation is obtained by multiplying r SB with the yield in the normalization region, defined by the nominal σ ηη (γ) requirement and the inverted criteria on the photon charged hadron isolation, I chg (γ) > 1.141 GeV. The expected yields with genuine photons, misidentified electrons, and multijet events are subtracted from the observation in the normalization region. The procedure is carried out separately for lepton flavors, N j selections, data-taking period, and for each bin of the differential cross section. The deviation from unity of the double-ratio of r SB to the corresponding ratio in the nominal σ ηη (γ) selection, stemming from the residual correlation between the two variables, is computed from simulation and it amounts to 18%. This value is used to correct the prediction.

Misidentified electron and genuine photon backgrounds
The background from electrons that are misidentified as photons is obtained from control regions defined by the requirements of |m(e, γ) − m Z | ≤ 10 GeV, and exactly three (misDY3), or four or more (misDY4p) jets. In the simulation, these event samples have a combined purity of 58% of Drell-Yan events with Z → ee, where one of the electrons passes the photon selection criteria. The simulated yield of the background component with a misidentified electron is multiplied by the scale factor (SF) defined below, separately for each of the three data-taking periods.
The Wγ (Zγ) process contributes to the LM3p regions for both lepton flavors and has a purity of 41% (21%). In the HM3p regions, the Wγ background is dominant with a purity  of 51%. The SFs for the misidentified electron background and the normalization of the Wγ and Zγ processes are obtained from the likelihood fit as described in section 6.2. The fit includes the data-based multijet estimates. The normalization of the Wγ process is left floating and the normalization of the Zγ process is allowed to vary within its uncertainty. The resulting m( , γ) distributions are shown in figure 4 in the N j ≥ 3 control regions. The background with misidentified electrons is dominant in the misDY3 and misDY4p regions close to the m Z peak. A correction of 15% to the normalization of the Drell-Yan process is measured in a data sample with two well-identified leptons satisfying |m( , ) − m Z | ≤ 10 GeV and N j ≥ 3, and is included in these results.
A summary of the extracted SF values for the misidentified electron background and the normalization of the Zγ and Wγ processes, obtained from a profile likelihood fit excluding the signal regions, is provided in table 4. The observed differences in the SFs for misidentified electrons are a result of the pixel detector replacement in 2017 and its operating conditions in the three data-taking periods. The stability of the procedure to estimate the yields of misidentified electrons and genuine photons is assessed by repeating the fit on individual data-taking periods and separately for the N j = 3 and ≥4 selections. The extracted SFs from these checks agree within the uncertainties. For the measurements of the inclusive and differential cross sections, as well as for setting EFT limits, the SFs are determined in situ by performing the fit simultaneously with the signal regions. Zγ normalization 1.01 ± 0.10 Wγ normalization 1.13 ± 0.08 Table 4. Extracted SFs and the total uncertainty obtained from the likelihood fit for the contribution from misidentified electrons for the three data-taking periods, and for the normalization of the Zγ and Wγ background components.

Systematic uncertainties
The systematic uncertainties affecting the signal selection efficiency and background yields are summarized in table 5. The table shows the range of variations in the different bins of the analysis caused by each systematic uncertainty in the signal and background yields, as well as an estimate of the impact of each uncertainty in the measured inclusive cross section. The table also indicates whether the uncertainties are treated as uncorrelated or fully correlated among the data-taking periods. The integrated luminosities of the 2016, 2017, and 2018 data-taking periods are individually known with uncertainties in the 2.3-2.5% range [70-72], while the total Run 2 (2016-2018) integrated luminosity has an uncertainty of 1.8%, the improvement in precision reflecting the (uncorrelated) time evolution of some systematic effects. The uncertainty in the inclusive cross section from these sources is, therefore, 1.8%. Simulated events are reweighted according to the distribution of the number of interactions in each bunch crossing corresponding to a total inelastic pp cross section of 69.2 mb [53]. The uncertainty in the total inelastic pp cross section is 4.6% [73] and affects the pileup estimate. The uncertainty due to the pileup effect is about 2% for the expected yields and less than 0.5% for the inclusive cross section.
The uncertainties in the SFs used to match the simulated trigger selection efficiencies to the ones observed in data are propagated to the results. From the "tag-and-probe" measurement [56, 59], an uncertainty of up to 0.5% is assigned to the yields obtained in simulation. Lepton selection efficiencies are measured in bins of lepton p T and η, and are found to be in the range 50-80 (75-85)% for electrons (muons). These measurements are performed separately in data and simulation and their ratio is used to scale the yields obtained in the simulation. The impact of these uncertainties on the inclusive cross section is 0.5 (0.7)% for the electron (muon) channel.
In the barrel section of the ECAL, an energy resolution of about 1% is achieved for unconverted or late-converting highly energetic photons in the tens of GeV energy range. Total systematic uncertainty 6.0 Statistical uncertainty 0.9 Total 6.0 Table 5. Breakdown of the total uncertainty in its statistical and systematic components in the different signal regions. The first column indicates the source of the uncertainty. The second column shows the correlation between the data-taking periods. The third column shows the typical pre-fit uncertainties in the total simulated yields in the signal region. The last column gives the corresponding systematic uncertainty in the tt γ cross section from the fit to the data. During the 2016 and 2017 data-taking periods, a gradual shift in the timing of the inputs of the ECAL L1 trigger in the forward endcap region (|η| > 2.4) led to a specific inefficiency (labeled "L1 prefiring" in table 5). A correction for this effect was determined using an unbiased data sample and is found to be relevant in events with jets with 2.4 < |η| < 3.0 and p T > 100 GeV. While no reconstructed objects at this η directly enter the measurements, it can affect the p miss T observable. A systematic variation of 20% of this correction for affected objects leads to an uncertainty of 0.3-0.9% in the predicted yields.
To estimate the theoretical uncertainties from missing higher-order corrections in the signal cross section calculation, the choice of µ R and µ F are varied independently up and down by a factor of 2. The acceptance variations are taken as the systematic uncertainty in each bin and are found to be smaller than 4.7%. Two independent nuisance parameters are used for the uncertainty in the choice of µ R and µ F , and their impact on the inclusive cross section measurement in the profile likelihood fit is less than 0.5%. A test with a single nuisance parameter, associated with the envelope of the uncertainties related to the choice of µ R and µ F , leads to negligible differences. The different sets in the NNPDF PDF [49] are used to estimate the corresponding uncertainty in the acceptance for the cross section measurement, which is less than 0.5%. The scale, PDF, and α S uncertainties in the inclusive fiducial cross section of the ttγ process, evaluated with MadGraph5_amc@nlo at NLO in QCD, amount to 17.5%. The limited number of available simulated events is considered by performing the fit using the Barlow-Beeston method [76].
In the parton shower simulation, the uncertainty from the choice of µ F is estimated by varying the scale of initial-and final-state radiation (ISR/FSR) up and down by factors of 2 and √ 2, respectively, as suggested in ref. [45]. The default configuration in pythia includes . The total uncertainty from color reconnection modeling is estimated by taking the maximum deviation from the nominal result and amounts to less than 0.5% in the inclusive cross section. The tWγ background component amounts to at most 3.3% of the total event yield in the SR3 and SR4p signal regions and is predicted by the tW sample, simulated with powheg at NLO precision. To account for uncertainties in the tWγ modeling, we treat the difference between the nominal prediction from the parton shower in the tW sample, normalized to NNLO, and a prediction obtained from MadGraph5_amc@nlo at LO for the 2 → 3 process as an uncertainty. For the SR3 (SR4p) signal regions, the differences of the total tWγ contribution are less than 44% (30%) in the p T (γ) bins, less than 34% (27%) in the |η(γ)| bins, and less than 19% (17%) in the ∆R( , γ) bins and lead to an uncertainty of 1.6% in the inclusive fiducial cross section.
The uncertainty in the normalization of the QCD multijet component is based on the variation of the TF with N j for different N b and amounts to 50%. Independent uncertainties are considered for the contributions to the N b = 0 and ≥1 yields. These have a significant impact only in the LM3, LM4p, HM3, and HM4p control regions, and lead to an uncertainty of 0.9% in the measured inclusive cross section.
The uncertainty in the nonprompt photon prediction is based on the modeling of the I chg (γ) distribution for different requirements on σ ηη (γ) and leads to an uncertainty of 1.8% in the inclusive cross section. The normalization of the Wγ process is left floating in the profile likelihood. To account for the uncertainty in the N j modeling of the Zγ process, we include an uncertainty of 30% in its normalization. In the signal region, the contribution of Wγ and Zγ background events generated with additional b or c quarks is 30%, and we assign an uncertainty of 20% in its normalization. Moreover, 40 (20)% uncertainty is assigned to the normalization of the Zγ (Wγ and misidentified electron) background in the N j ≥ 4 signal and control regions. The corresponding impact of the normalization of the Zγ and Wγ contributions are 0.5 and 2.3%, respectively. The component with misidentified electrons leads to an uncertainty of up to 8% in the predicted background yields with an impact on the inclusive cross section of 1.8%. The 8% uncertainty in the normalization of the Drell-Yan process, the 5% uncertainty in the t/tt normalization, and the uncertainties in the normalization of other small background components lead to additional uncertainties below 1%.

Inclusive cross section measurement
The observed data, as well as the predicted signal and background yields resulting from the likelihood fit to all signal and control regions, are shown in figures 5 and 6. In these figures, the contributions from the three data-taking periods are summed, accounting for    Table 6. The observed number of events for the SR3 and SR4p signal regions in the e and µ channels, and the predicted yields and total post-fit uncertainties in each background component. the correlation of the systematic uncertainties. The signal cross section is extracted from these categories using the statistical procedure detailed in section 6.2. In the fit, nuisance parameters for the various systematic uncertainties and the normalization of background processes, as described in section 8, are included. The theoretical uncertainty in the inclusive fiducial cross section does not enter the likelihood fit for the inclusive or differential cross section measurements. Using three bins in M 3 reduces the uncertainty in the backgrounds without a hadronically decaying top quark, e.g., the misidentified electron background and the Wγ and Zγ processes, and decreases the total relative uncertainty in the inclusive cross section from 6.7 to 6.0%. The observed number of events for the SR3 and SR4p signal regions in the e and µ channels, and the predicted yields and total uncertainties in each background component are listed in table 6. Figure 7 shows the ranking of the leading systematic uncertainties according to their post-fit impact on the measured inclusive cross section. The normalization of the Wγ background is the largest individual contribution to the uncertainty in the inclusive cross section measurement and amounts to about 2.3%. The post-fit values of the nuisance parameters (pulls) are also shown and are found to lie within the pre-fit uncertainties. The  A comparison of the measured cross sections and the SM prediction is shown in figure 8, providing also the measurements for different choices of N j and the lepton flavor. For the latter results, the likelihood fit is performed separately in the corresponding set of signal regions and the full set of control regions.
Results are obtained simultaneously for the electron and muon channels, the 3 jet and ≥4 jet bins, and for the three data-taking periods. The binning in the SR3 and SR4p selections for the measurement of the differential distributions at the reconstruction level is shown in table 7.

JHEP12(2021)180
As described in section 6.2, the same control regions are used for the inclusive and differential cross section measurements. The signal strength is left floating in the profile likelihood fit separately for each of the differential bins, the N j selection, the lepton flavor, and the data-taking period. The procedure has been tested to reproduce ad-hoc modifications of the simulated signal prediction within the numerical accuracy. The fit is performed separately for each differential distribution.
The distributions of the observables after background subtraction are further unfolded to the fiducial particle level phase space defined in section 5. The unfolded differential cross section is defined in the same phase space as the inclusive cross section reported above, i.e., in the phase space where the top quark pair is produced in association with a photon satisfying p T (γ) > 20 GeV and |η(γ)| < 1.4442. Signal events that are not generated within the fiducial region amount to 5-10% and are subtracted based on simulation. In the simulation, p T (γ) is taken as the transverse momentum after accounting for the effects of QCD and electroweak radiation.
The ttγ MadGraph5_amc@nlo MC sample is used to construct a response matrix that takes into account both detector response and acceptance corrections. The same corrections, SFs, and uncertainties as used in the inclusive cross section are applied. Because of the high momentum-and angular resolutions of photons and leptons, the fraction of events migrating from a specific momentum region at the particle level to another one at the reconstruction level is small for all unfolded distributions. Under such conditions, and with the chosen bin size, no regularization term is required [79]. The TUnfold package [80] is used to obtain the results for the three measured observables using matrix inversion. The binning in the fiducial region is chosen such that two bins at the reconstruction level correspond to one bin in the fiducial region for most cases. This choice provides stability to the unfolding algorithm. The linearity of the unfolding procedure is tested by unfolding suitably reweighted simulated reconstruction-level yields. Differences between the unfolded reweighted distributions and the distributions resulting from the reweighting applied at the fiducial level are found to be negligible.
Uncertainties in the estimated signal yields are propagated through the unfolding procedure, including the effects on the response matrix. Experimental uncertainties from the detector response and efficiency, such as the photon identification, JES, and b tagging uncertainties, are applied as a function of the reconstructed observable. The differential cross sections, obtained by this procedure, are shown in figure 9. It includes a comparison with simulation obtained from MadGraph5_amc@nlo interfaced to her-wig++ [81] v2.7.1 with the EE5C tune [46] and to herwig7 v7.1.4 with the CH3 tune [82] for the parton shower and hadronization. The inclusive fiducial cross section predicted by herwig++ (herwig7) is 8.3% (5.4%) lower than for the nominal simulation.
The bin efficiency, defined as the fraction of generated events that are reconstructed in the corresponding bins at reconstruction level, is in the range of 20-30%. The bin purity, defined as the fraction of reconstructed events that originate from the corresponding bin at the particle level, is in the range of 85-90%. For p T (γ) > 120 GeV, the uncertainties in the JES, the photon identification efficiency, and the color reconnection modeling are the largest sources of systematic uncertainty. The correlation matrices of the systematic uncer- tainties for the unfolded differential measurements are shown in figure 10. The correlations are lower in the tail of p T (γ) due to larger statistical uncertainties in the simulation. The first bin of the ∆R( , γ) measurement is less affected by uncertainties in the normalization of backgrounds, resulting in slightly lower correlations in this case. All correlations from statistical uncertainties originating from the data are below 7%. Including the uncertainty in the fiducial signal cross section, we perform a compatibility test of the unfolded distribution and the nominal prediction. The corresponding χ 2 test statistic evaluates to 12.0 with 9 degrees of freedom (dof) for the p T (γ) distribution, 5.2 with 5 dof for |η(γ)|, and 6.3 with 7 dof for ∆R( , γ).

Effective field theory interpretation
Many BSM models predict anomalous couplings of the top quark to the electroweak gauge bosons [83][84][85][86][87][88][89]. The differential cross section measurement is interpreted at the reconstruction level in SM-EFT in the Warsaw basis [90], formed by 59 baryon number conserving dimension-six Wilson coefficients. Among them, 15 are relevant for top quark interactions [91]. Anomalous interactions between the top quark and the gluon (chromomagnetic and chromoelectric dipole moment interactions) are tightly constrained by the tt+jets measurements [92,93]. Similarly, the modification of the Wtb vertex is best constrained by measurements of the W helicity fractions in top quark pair production [94] and in t-channel single top quark production [95].
The Wilson coefficients in the Warsaw basis inducing electroweak dipole moments are denoted by C (33) uB and C (33) uW [12]. The SM gauge symmetry provides the ttZ and the ttγ final states with complementary constraining power [1][2][3][4]. The linear relations c tZ = Re − sin θ W C (33) uB + cos θ W C (33) uW , uB + cos θ W C (33) uW , c tγ = Re cos θ W C (33) uB + sin θ W C (33) uW , c I tγ = Im cos θ W C (33) uB + sin θ W C (33) uW , express the modifications of the ttZ interaction vertex, c tZ and c I tZ , and of the ttγ interaction vertex, c tγ and c I tγ , in the Warsaw basis. The constraint C (33) uW = 0 ensures a SM Wtb vertex. Under this assumption, c tZ (c I tZ ) and c tγ (c I tγ ) are dependent and we choose the former to parametrize the BSM hypothesis.
The spectrum of p T (γ) is a sensitive probe to such modifications. Other observables, e.g., |η(γ)| or ∆R( , γ), are found to be largely insensitive. Wilson coefficients that are not considered in this work are kept at their SM values and the SM-EFT expansion parameter is set to a mass scale Λ = 1 TeV. Using the SM-EFT parametrization from ref. [12], simulated samples at the particle level are produced with nonzero values of the Wilson coefficients c tZ and c I tZ . The ttγ signal process and all background processes affected by c tZ or c I tZ at the ME level are included in the simulation. These samples are used to reweight the nominal simulation in the fiducial phase space using the quadratic parametrization detailed in ref. [96]. The reweighting procedure is validated at the reconstruction level with a reduced set of statistically independent samples for nonzero values of c tZ and c I tZ and excellent agreement is found.
The SR3 and SR4p signal regions and the p T (γ) boundaries defining the bins in table 7 are used to construct a binned likelihood function L(θ) as a product of Poisson probabilities from the yields in the signal and control regions. The nuisance parameters are labeled by θ and the profile likelihood ratio q = −2 ln(L(θ, C)/L(θ max )) is the test statistic. Here,θ is the set of nuisance parameters maximizing the likelihood function at a BSM point defined by the Wilson coefficients collectively denoted by C. In the denominator,θ max maximizes the likelihood function in the BSM parameter space. The ttγ signal is normalized according to the SM expectation at NLO in QCD and its uncertainty is included as a nuisance.    Figure 12 displays the one-dimensional scans of the coefficients. In the upper row, one Wilson coefficient is scanned, while the other is profiled. The lower row shows the scans, where the second Wilson coefficient is set to zero. The second local minima in the scans of the log-likelihood as a function of c tZ and c I tZ , visible in figure 12 (lower row), is the result of a mild tension with the SM hypothesis in conjunction with the similarity of the predictions for Wilson coefficients with opposite sign. The corresponding one-dimensional intervals at 68 and 95% confidence interval (CL) are listed in table 8 and are more stringent than previous limits obtained from ttZ final states [97,98]. Models with nonzero electroweak dipole moments predict a harder p T (γ) spectrum that is not observed in data. Figure 13 shows the best fit result in the two-dimensional plane spanned by c tZ and c I tZ and the log-likelihood scan. The SM prediction is within the 95% CL of the best fit value of the c tZ and c I tZ coefficients. In figure 14

Summary
A measurement of the cross section for the top quark pair production in association with a photon using a data sample of proton-proton collisions at √ s = 13 TeV, corresponding to an integrated luminosity of 137 fb −1 , collected with the CMS detector at the LHC has been presented. It is the first result of the CMS Collaboration on measurements in the ttγ final state using 13 TeV data. The analysis has been performed in the single-lepton channel with events with exactly three and four or more jets among which at least one is b tagged. Background components with misidentified electrons, photons originating in the hadronization of jets, the multijet component, and prompt photons from the Wγ and Zγ processes are estimated from data. The measured inclusive cross section in a fiducial region with photon transverse momentum p T (γ) > 20 GeV and jet multiplicity greater than three is measured to be 798 ± 7 (stat) ± 48 (syst) fb, in good agreement with the standard model prediction at next-to-leading order in quantum chromodynamics. The shading quantified by the color scale on the right reflects the negative log-likelihood ratio with respect to the best fit value that is designated by the star. The green and orange lines indicate the 68 and 95% CL contours from the fit, respectively. The allowed areas are those between the two green contours and that inside the orange contour. The dot shows the SM prediction.

JHEP12(2021)180
Differential cross sections for p T (γ) and absolute value of the photon pseudorapidity, as well as for the angular separation of the lepton and the photon, have been measured and unfolded to particle level in the same fiducial volume. The comparison to simulation was performed using different showering algorithms. The measurements are also interpreted in terms of limits on the Wilson coefficients in the context of the standard model effective field theory. The confidence intervals for the Wilson coefficients c tZ and c I tZ are the most stringent to date.

Acknowledgments
We congratulate our colleagues in the CERN accelerator departments for the excellent performance of the LHC and thank the technical and administrative staffs at CERN and at other CMS institutes for their contributions to the success of the CMS effort. In addition, we gratefully acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid and other centers for delivering so effectively the computing infrastructure essential to our analyses. Finally, we acknowledge the enduring support for the construction and operation of the LHC, the CMS detector, and the supporting computing infrastructure provided by the following funding agencies: BMBWF and FWF (Austria); 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.
[6] ATLAS collaboration, Observation of top-quark pair production in association with a photon and measurement of the ttγ production cross section in pp collisions at √ s = 7 TeV using the ATLAS detector, Phys. Rev. D 91 (2015)