Abstract The rate for Higgs ( $${\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>H</mml:mi> </mml:math> ) bosons production in association with either one ( $${\mathrm{t}} {\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>t</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> </mml:math> ) or two ( $${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>t</mml:mi> <mml:mover> <mml:mi>t</mml:mi> <mml:mo>¯</mml:mo> </mml:mover> <mml:mi>H</mml:mi> </mml:mrow> </mml:math> ) top quarks is measured in final states containing multiple electrons, muons, or tau leptons decaying to hadrons and a neutrino, using proton–proton collisions recorded at a center-of-mass energy of $$13\,\text {Te}\text {V} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mn>13</mml:mn> <mml:mspace /> <mml:mrow> <mml:mtext>Te</mml:mtext> <mml:mspace /> </mml:mrow> </mml:mrow> </mml:math> by the CMS experiment. The analyzed data correspond to an integrated luminosity of 137 $$\,\text {fb}^{-1}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mspace /> <mml:msup> <mml:mtext>fb</mml:mtext> <mml:mrow> <mml:mo>-</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> . The analysis is aimed at events that contain $${\mathrm{H}} \rightarrow {\mathrm{W}} {\mathrm{W}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>H</mml:mi> <mml:mo>→</mml:mo> <mml:mi>W</mml:mi> <mml:mi>W</mml:mi> </mml:mrow> </mml:math> , $${\mathrm{H}} \rightarrow {\uptau } {\uptau } $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>H</mml:mi> <mml:mo>→</mml:mo> <mml:mi>τ</mml:mi> <mml:mi>τ</mml:mi> </mml:mrow> </mml:math> , or $${\mathrm{H}} \rightarrow {\mathrm{Z}} {\mathrm{Z}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>H</mml:mi> <mml:mo>→</mml:mo> <mml:mi>Z</mml:mi> <mml:mi>Z</mml:mi> </mml:mrow> </mml:math> decays and each of the top quark(s) decays either to lepton+jets or all-jet channels. Sensitivity to signal is maximized by including ten signatures in the analysis, depending on the lepton multiplicity. The separation among $${\mathrm{t}} {\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>t</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> </mml:math> , $${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>t</mml:mi> <mml:mover> <mml:mi>t</mml:mi> <mml:mo>¯</mml:mo> </mml:mover> <mml:mi>H</mml:mi> </mml:mrow> </mml:math> , and the backgrounds is enhanced through machine-learning techniques and matrix-element methods. The measured production rates for the $${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>t</mml:mi> <mml:mover> <mml:mi>t</mml:mi> <mml:mo>¯</mml:mo> </mml:mover> <mml:mi>H</mml:mi> </mml:mrow> </mml:math> and $${\mathrm{t}} {\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>t</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> </mml:math> signals correspond to $$0.92 \pm 0.19\,\text {(stat)} ^{+0.17}_{-0.13}\,\text {(syst)} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mn>0.92</mml:mn> <mml:mo>±</mml:mo> <mml:mn>0.19</mml:mn> <mml:mspace /> <mml:msubsup> <mml:mtext>(stat)</mml:mtext> <mml:mrow> <mml:mo>-</mml:mo> <mml:mn>0.13</mml:mn> </mml:mrow> <mml:mrow> <mml:mo>+</mml:mo> <mml:mn>0.17</mml:mn> </mml:mrow> </mml:msubsup> <mml:mspace /> <mml:mtext>(syst)</mml:mtext> </mml:mrow> </mml:math> and $$5.7 \pm 2.7\,\text {(stat)} \pm 3.0\,\text {(syst)} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mn>5.7</mml:mn> <mml:mo>±</mml:mo> <mml:mn>2.7</mml:mn> <mml:mspace /> <mml:mtext>(stat)</mml:mtext> <mml:mo>±</mml:mo> <mml:mn>3.0</mml:mn> <mml:mspace /> <mml:mtext>(syst)</mml:mtext> </mml:mrow> </mml:math> of their respective standard model (SM) expectations. The corresponding observed (expected) significance amounts to 4.7 (5.2) standard deviations for $${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>t</mml:mi> <mml:mover> <mml:mi>t</mml:mi> <mml:mo>¯</mml:mo> </mml:mover> <mml:mi>H</mml:mi> </mml:mrow> </mml:math> , and to 1.4 (0.3) for $${\mathrm{t}} {\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>t</mml:mi> <mml:mi>H</mml:mi> </mml:mrow> </mml:math> production. Assuming that the Higgs boson coupling to the tau lepton is equal in strength to its expectation in the SM, the coupling $$y_{{\mathrm{t}}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>y</mml:mi> <mml:mi>t</mml:mi> </mml:msub> </mml:math> of the Higgs boson to the top quark divided by its SM expectation, $$\kappa _{{\mathrm{t}}}=y_{{\mathrm{t}}}/y_{{\mathrm{t}}}^{\mathrm {SM}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msub> <mml:mi>κ</mml:mi> <mml:mi>t</mml:mi> </mml:msub> <mml:mo>=</mml:mo> <mml:msub> <mml:mi>y</mml:mi> <mml:mi>t</mml:mi> </mml:msub> <mml:mo>/</mml:mo> <mml:msubsup> <mml:mi>y</mml:mi> <mml:mrow> <mml:mi>t</mml:mi> </mml:mrow> <mml:mi>SM</mml:mi> </mml:msubsup> </mml:mrow> </mml:math> , is constrained to be within $$-0.9&lt; \kappa _{{\mathrm{t}}}&lt; -0.7$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mo>-</mml:mo> <mml:mn>0.9</mml:mn> <mml:mo>&lt;</mml:mo> <mml:msub> <mml:mi>κ</mml:mi> <mml:mi>t</mml:mi> </mml:msub> <mml:mo>&lt;</mml:mo> <mml:mo>-</mml:mo> <mml:mn>0.7</mml:mn> </mml:mrow> </mml:math> or $$0.7&lt; \kappa _{{\mathrm{t}}}&lt; 1.1$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mn>0.7</mml:mn> <mml:mo>&lt;</mml:mo> <mml:msub> <mml:mi>κ</mml:mi> <mml:mi>t</mml:mi> </mml:msub> <mml:mo>&lt;</mml:mo> <mml:mn>1.1</mml:mn> </mml:mrow> </mml:math> , at 95% confidence level. This result is the most sensitive measurement of the $${\mathrm{t}} {{\overline{{{\mathrm{t}}}}}} {\mathrm{H}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>t</mml:mi> <mml:mover> <mml:mi>t</mml:mi> <mml:mo>¯</mml:mo> </mml:mover> <mml:mi>H</mml:mi> </mml:mrow> </mml:math> production rate to date.