\section{Introduction}

Recently, F.  A. Mesmer \index{Mesmer}
has reported evidence of a profound influence of
magnets on a variety of aspects of human and animal physiology~\cite{Mesmer}.
Numerous authors have claimed to reproduce certain of the phenomena
reported by Mesmer~\cite{diCenzo,Muller}, but there is little understanding
of the precise mode of action of the magnetic influence.  In this
article, I will present evidence that mesmeric influence can disrupt the
lethal action of sodium cyanide (NaCN) when modulated by a sequence of
high-pitched tones.  After a description of the phenomenon, I will comment
on the implications for the nature of the magnetic coupling.

\section{Observations}

My observations began when I was called on an emergency basis to aid two
Albanian noblemen, Guglielmo B. and Ferrando di N., who had
ingested the poison in an apparent attempt to commit suicide.  (The facts
of the case were rather mysterious, since these young men had been seen in the
harbor earlier in the day, affecting a carefree attitude.  However, this
need not concern us here.)  The standard remedy would have been to induce
vomiting; however, this would have offended the
sensibilities \index{sensibilities}  of two
noble ladies present, the daughters of the duke in whose household this
unfortunate incident took place.  Thus, I thought it wise to begin by
administering magneto-therapy using the devices that I had at hand.

\begin{figure}[htb]
\centering
\includegraphics[height=1.5in]{magnet}
\caption{Plan of the magnet used in the mesmeric studies.}
\label{fig:magnet}
\end{figure}

The source of magnetic influence used in this therapy was the simple
iron pole magnet shown in Figure~\ref{fig:magnet}.  The device, roughly
the size of human palm, was held three fingers' length away from each
of the patients, while the blood cyanide level was monitored.  As this
operator proceeded in silence, only a very small effect on the poison was
observed.  In a moment, however, the sisters began to weep and ululate
in the most painful manner.  Under this blanket of sound, the effects of
the cyanide reaction were observed to be completely reversed.  Indeed,
even before the end of a single session of treatment, the Albanians were
seen to wake from their comas and offer thanks for their rescue.  Some
sample blood cyanide levels are reported in Table~\ref{tab:blood}.

\begin{table}[t]
\begin{center}
\begin{tabular}{l|ccc}
Patient &  Initial level($\mu$g/cc) &  w. Magnet &
w. Magnet and Sound \\ \hline
 Guglielmo B.  &   0.12     &     0.10      &     0.001  \\
 Ferrando di N. &  0.15     &     0.11      &  $< 0.0005$ \\ \hline
\end{tabular}
\caption{Blood cyanide levels for the two patients.}
\label{tab:blood}
\end{center}
\end{table}

\section{Interpretation}

A description of these observations, obviously based on second-hand reports,
has already been reported by L. da Ponte~\cite{daPonte}. \index{da Ponte}
 In his paper,
da Ponte attributes the mode of action of the magnetic force to direct
magnetic attraction.  This would be implausible with the relatively small
magnet used in this study; thus da Ponte was forced to propose an absurdly
large size for the  magnetic device.  His theory does not account for the
acoustic influence, which is known not to be present in magnetic action on
inanimate objects.  A second hypothesis would state that the acoustic
fields perturb the chemistry of the cyanide reactions, allowing the
reagents to couple magnetically.   As with the first hypothesis, we are
skeptical, since again this phenomenon has no inorganic counterpart.  Thus,
we are forced to conclude, with Mesmer, that magnetism couples directly to
the life force, which, as is well known, can be modulated by musical
influences.

In conclusion, we note that the unusual mode of action of magnetism is only
one of many which have been reported in the literature.  It is likely
that there are further unusual magnetic phenomena to be discovered.  All
studies to date have involved male patients.  But it is likely that the
prevailing belief that `women are all the same' would be overturned in a
similar investigation.

\bigskip  \bigskip \begin{center} \begin{large
I am grateful to Don Alfonso d'Alba for certain services essential to
this investigation.

\section{Introduction}

All measurements of Higgs boson properties are consistent with the Standard Model (SM) predictions, but the fermion sector was so far less tested.
In the SM, couplings between the Higgs boson and fermions occur through Yukawa interactions, whose strength is proportional to the fermion mass.
The top quark Yukawa coupling is therefore predicted to be by far the largest one, with a value close to unity. Indirect constrains on this coupling can be obtained
via the Higgs boson production by gluon--gluon fusion or its decay to two photons, which proceed through loops. However, New Physics (NP) could also enter in these loops,
hence a direct measurement is needed to disentangle and probe NP. Such a direct measurement can be achieved through the associated production of the Higgs boson and a top
quark pair (\ensuremath{t\bar t H}\xspace), which is a tree-level process and whose cross-section is proportional to the square of the top quark Yukawa coupling. Observing the \ensuremath{t\bar t H}\xspace production would therefore
allow to simultaneously constrain this coupling and possible Beyond the Standard Model (BSM) effects.\\

The \ensuremath{t\bar t H}\xspace production is quite a rare process, with a cross-section at the Large Hadron Collider (LHC) of $\sim 0.5$~pb at a center-of-mass energy of \ensuremath{\sqrt{s}=13~\mathrm{TeV}}.
This is typically 1\% of the inclusive Higgs boson production cross-section and less than 0.1\% of the \ensuremath{t\bar t}\ production one. This rare production has to be convolved with all possible combinations
of \ensuremath{t\bar t}\ and Higgs boson decays, giving rise to a very rich spectrum of possible signatures. A broad spectrum of analyses has therefore been designed to cover these multiple final states,
with their selections designed to avoid any overlap between them.
They will be described in the next sections together with their combination, starting from the most abundant ones (\ensuremath{t\bar t H}\xspace, \Htobb and \ensuremath{t\bar t H}\xspace, \Htomulti), which use 36.1~\mbox{fb$^{-1}$}\
of proton--proton collision data recorded by the ATLAS experiment~\cite{atlas} at a center-of-mass energy of \ensuremath{\sqrt{s}=13~\mathrm{TeV}}, to the more rare ones (\ensuremath{t\bar t H}\xspace, \Htoyy and \ensuremath{t\bar t H}\xspace, \Htofourl ($l=e, \mu$)),
which use 79.8~\mbox{fb$^{-1}$}\ of data.

\section{Measurement in the \Htobb channel}

The search for the \ensuremath{t\bar t H}\xspace production in the \Htobb channel~\cite{bb} benefits from a large branching ratio,
but suffers from a huge combinatoric background at the event reconstruction level and a large irreducible background
(mostly \ensuremath{t\bar t}\ $+$ heavy flavour jets) associated with big theoretical uncertainties. The modelling of this background and the definition of the
associated uncertainty scheme are performed using state-of-the-art Monte Carlo (MC) generators and constraints from data in background-enriched regions.
The nominal \ensuremath{t\bar t}\ MC sample uses an inclusive next-to-leading order (NLO) generator normalised to predictions at next-to-next-to-leading order in perturbative QCD including resummation of
next-to-next-to-leading logarithmic soft gluon terms. The main associated
systematic uncertainties are based on comparisons against alternative MC samples, using different generators, parton shower or hadronisation models.
The normalisations of \ensuremath{t\bar t}+$\ge1b$ and \ensuremath{t\bar t}+$\ge1c$ backgrounds are derived from data.
To define the analysis regions, a selection of semi- and di-leptonic \ensuremath{t\bar t}\ decays is first performed, based on the presence of one or two light leptons,
followed by a further categorisation based on the jet multiplicity and on the value of the $b$-tagging discriminant of each jet (four $b$-tagging working points are used).
This allows to define 9 signal-enriched regions (SR) and 10 background-enriched control regions (CR), with very different signal purities and background fractions
(a ``boosted'' category targeting high $p_\mathrm{T}$ top/Higgs candidates is also defined). The best signal purity in these regions is only $\sim$5\%, so a cascade
of Multivariate Analysis (MVA) classifiers is used to further separate the signal from the backgrounds. A \ensuremath{t\bar t H}\xspace system reconstruction Boosted Decision Tree (BDT)
is designed in most SR to solve the combinatorics to identify the best parton-jet assignments and reconstruct the final state. It is complemented
with a Likelihood discriminant in the single-lepton SR, built from all different jet permutations, and a Matrix Element Method in the single-lepton SR with the best purity.
The outputs of all these classifiers, together with the related Higgs/top candidate properties, event kinematic variables and $b$-tag information, are input into a classification BDT
built in each SR to best separate the signal from the backgrounds. A combined profile likelihood fit is performed across all analysis regions to extract the signal yield.
The significance of the observed (expected) signal is $1.4\sigma$ ($1.6\sigma$). The measured signal strength is in good agreement with the SM prediction within the relative measurement
precision of 75\%, as seen in Figure~\ref{fig_1} (left). This precision is dominated by the systematic uncertainties, in particular those associated with \ensuremath{t\bar t}+$\ge1b$ modelling.

\section{Measurement in multilepton final states}

The search for \ensuremath{t\bar t H}\xspace production in multilepton final states~\cite{ml} targets \HtoWW, \Htott and \HtoZZ decays. The \ensuremath{t\bar t}\ background is suppressed by requiring two same-sign or at least 3 leptons.
Seven orthogonal channels, categorised by the multiplicity of light leptons and hadronically decaying $\tau$ ($\tau_{\mathrm{had}}$) candidates, are defined.
They cover a wide range of yields and S/B purities, with different background composition in the various regions.
The irreducible background, made of SM processes with prompt leptons, is dominated by \ensuremath{t\bar t W}\xspace, \ensuremath{t\bar t Z}\xspace and diboson productions.
It is estimated from NLO MC predictions and validated in data in various dedicated regions.
The reducible background is mainly coming from \ensuremath{t\bar t}\ events, with non-prompt light leptons from $b$-hadron decay, fake $\tau_{\mathrm{had}}$ candidates from
light flavour jet or charge misidentified electrons. These contributions are mostly estimated through data-driven techniques, and validated in various regions.
Dedicated MVA approaches are used in most SR to further enhance the signal to background separation. A combined profile likelihood fit across all SR and CR extracts
the signal with an observed (expected) significance of $4.1\sigma$ ($2.8\sigma$). The measured signal strength is in good agreement with the SM prediction within the relative measurement
precision of 30\%, as seen in Figure~\ref{fig_1} (right). This precision is fairly shared between statistical and systematic errors, dominated by uncertainties
on \ensuremath{t\bar t H}\xspace and \ensuremath{t\bar t V}\xspace ($V=W,Z$) modelling and those associated to non-prompt lepton estimates (with a large contribution from the limited CR statistics).
Similar results are obtained with an alternative fit using a free floating \ensuremath{t\bar t V}\xspace normalisation, albeit with a 15\% degraded sensitivity.

\begin{figure}[htb]
\includegraphics[width=0.52\linewidth]{fig_01.pdf}
\includegraphics[width=0.48\linewidth]{fig_02.pdf}
\vspace*{-1cm}
\caption{Observed best-fit values of the \ensuremath{t\bar t H}\xspace signal strength $\mu$ and their uncertainties by final-state category and combined, for the \ensuremath{t\bar t H}\xspace, \Htobb (left~\cite{bb}) and multileptons (right~\cite{ml}) analyses.
The individual $\mu$ values for the channels are obtained from a simultaneous fit with the signal-strength parameter for each channel floating independently. The SM prediction is $\mu=1$.
\label{fig_1}}
\end{figure}

\section{Measurement in the \Htoyy and \Htofourl ($l=e, \mu$) channels}
The \ensuremath{t\bar t H}\xspace, \Htoyy channel~\cite{obs} is more rare, but offers a clean signature with a narrow peak in the diphoton invariant mass distribution over a smooth continuum background,
extrapolated from sideband regions.  Two high-$E_\mathrm{T}$ isolated photons are required, and events are categorised in leptonic
and hadronic regions depending on the presence of a light lepton. A further signal to background separation power is based on BDT discriminants,
trained with deep neural networks using \ensuremath{t\bar t H}\xspace MC and background events from data control regions. Seven categories (4 hadronic and 3 leptonic) are defined according to the value of the BDT output.
The signal over background ratio is $\sim1.9$ in the purest categories, with the background dominated by the diphoton continuum (including \ensuremath{t\bar t}$+\gamma\gamma$).
A combined unbinned fit to the diphoton invariant mass distributions in the 7 categories is performed (see Figure~\ref{fig_2} (left)), with the Higgs signal modelled with a double-sided Crystal Ball function,
and the continuum background with a functional form whose parameters are free in the fit. The significance of the observed (expected) signal is $4.1\sigma$ ($3.7\sigma$).
The measured signal strength is in good agreement with the SM prediction within the relative measurement
precision of 30\%, dominated by statistical errors. The dominant systematic uncertainties are associated to the signal parton shower model (8\%) and the photon energy resolution (6\%).\\

A similar strategy is performed in the \ensuremath{t\bar t H}\xspace, \Htofourl ($l=e, \mu$) channel~\cite{obs}, with a categorisation in leptonic and hadronic
regions and the use of a BDT discriminant in the hadronic category. The signal over background ratio is $\ge5$ in the purest categories, with the background dominated by \ensuremath{t\bar t V}\xspace
and non-\ensuremath{t\bar t H}\xspace Higgs production modes. The expected significance is $1.2\sigma$ (0.6 \ensuremath{t\bar t H}\xspace events), while 0 event are observed.

\section{Combination and results}

The combination~\cite{obs} of the searches for \ensuremath{t\bar t H}\xspace production described in the previous sections is based on a simultaneous fit to all signal and control regions of the individual analyses.
The contributions from the non-\ensuremath{t\bar t H}\xspace Higgs production modes are fixed to SM predictions.
Using \ensuremath{\sqrt{s}=13~\mathrm{TeV}}\ data only, an observed (expected) signal significance of $5.8\sigma$ ($4.9\sigma$) is achieved.
Combining with \ensuremath{\sqrt{s}=7~\mathrm{TeV}}\ and \ensuremath{\sqrt{s}=8~\mathrm{TeV}}\ data increases the observed (expected) signal significance to $6.3\sigma$ ($5.1\sigma$).
This constitutes an observation of the \ensuremath{t\bar t H}\xspace process, which was a major milestone for LHC Run-2.
The total \ensuremath{t\bar t H}\xspace cross-section at \ensuremath{\sqrt{s}=13~\mathrm{TeV}}\ is measured to be $670\pm90$ (stat) $^{+110} _{-100}$ (syst)~fb, in good agreement with the NLO SM prediction of $507^{+35} _{-50}$~fb.
The measurement precision of 20\% is slightly dominated by systematic uncertainties, in particular from \ensuremath{t\bar t}+heavy flavour modelling (9.9\%) and signal modelling (6.0\%).
The cross-section measurements performed at \ensuremath{\sqrt{s}=13~\mathrm{TeV}}\ and \ensuremath{\sqrt{s}=8~\mathrm{TeV}}\ are shown in Figure~\ref{fig_2} (right).\\

Measurements of Higgs boson couplings combining several ATLAS analyses, including those on \ensuremath{t\bar t H}\xspace reported here, have been performed~\cite{couplings}.
Assuming SM particles only ($gg\to H$ and \Htoyy loops are then resolved with SM particles), the top quark Yukawa coupling, relative to its SM value,
is measured to be $\kappa_{top}=1.03^{+0.12} _{-0.11}$. Allowing new particles to enter in the loops, an effective Higgs-to-gluon coupling is defined,
and a ratio $\kappa_{gluon}/\kappa_{top}=1.09\pm0.14$ is extracted. This is consistent with SM Higgs boson coupling predictions, and allows to put constrains on BSM contributions.

\begin{figure}[htb]
\includegraphics[width=0.55\linewidth]{fig_03.pdf}
\includegraphics[width=0.45\linewidth]{fig_04.pdf}
\vspace*{-1cm}
\caption{Left: Weighted diphoton invariant mass spectrum in the \ensuremath{t\bar t H}\xspace-sensitive BDT bins.
The solid red curve shows the fitted signal-plus-background model.
The non-resonant and total background components of the fit are shown with the dotted blue curve and dashed green curve.
Right: Measured \ensuremath{t\bar t H}\xspace cross-sections in $pp$ collisions at center-of-mass energies of 8 TeV and 13 TeV.
Both the total and statistical-only uncertainties are shown. The measurements are compared with the SM prediction. Ref.~\cite{obs}.
\label{fig_2}}
\end{figure}

\section{Conclusions}

The \ensuremath{t\bar t H}\xspace production process has been observed using $pp$ collision data produced by the LHC and recorded with the ATLAS detector.
The cross-section at \ensuremath{\sqrt{s}=13~\mathrm{TeV}}\ is measured with 20\% precision.
All measurements are in good agreement with SM predictions, allowing to put constrains on BSM contributions.
This constitutes a direct observation of the top quark Yukawa coupling, which, together with the recent observations of \Htott and \Htobb decays
by the ATLAS and CMS collaborations separately, firmly establishes the existence of Yukawa interactions.