<?xml version="1.0" encoding="UTF-8"?>
<article xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id>JHEP</journal-id>
<journal-title-group>
<journal-title>Journal of High Energy Physics</journal-title>
</journal-title-group>
<issn pub-type="epub">1029-8479</issn>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">JHEP05(2015)040</article-id>
<article-id pub-id-type="doi">10.1007/JHEP05(2015)040</article-id>
<title-group><article-title>Improving the sensitivity of stop searches with on-shell constrained invariant mass
variables
</article-title></title-group>
<contrib-group><contrib>
  <string-name>Won Sang Cho</string-name>
  <email>hepkosmos@ibs.re.kr</email>
  <xref ref-type="aff" rid="a1"/>
  <xref ref-type="aff" rid="a2"/>
  </contrib>
<contrib>
  <string-name>James S. Gainer</string-name>
  <email>jgainer137@gmail.com</email>
  <xref ref-type="aff" rid="a1"/>
  </contrib>
<contrib contrib-type="Corresponding author">
  <string-name>Doojin Kim</string-name>
  <email>immworry@ufl.edu</email>
  <xref ref-type="aff" rid="a1"/>
  </contrib>
<contrib>
  <string-name>Konstantin T. Matchev</string-name>
  <email>matchev@ufl.edu</email>
  <xref ref-type="aff" rid="a1"/>
  </contrib>
<contrib>
  <string-name>Filip Moortgat</string-name>
  <email>filip.moortgat@cern.ch</email>
  <xref ref-type="aff" rid="a3"/>
  </contrib>
<contrib>
  <string-name>Luc Pape</string-name>
  <email>luc.pape@cern.ch</email>
  <xref ref-type="aff" rid="a3"/>
  </contrib>
<contrib>
  <string-name>Myeonghun Park</string-name>
  <email>parc@apctp.org</email>
  <xref ref-type="aff" rid="a4"/>
  <xref ref-type="aff" rid="a5"/>
  <xref ref-type="aff" rid="a6"/>
  </contrib>

  <aff id="a1">Physics Department, University of Florida,<break/>Gainesville, FL 32611, U.S.A.</aff>

  <aff id="a2">Center for Theoretical Physics of the Universe, Institute for Basic Science
(IBS),<break/>Daejeon 305-811, Korea</aff>

  <aff id="a3">CERN,<break/>Geneva CH-1211, Switzerland</aff>

  <aff id="a4">Kavli IPMU (WPI), The University of Tokyo,<break/>Kashiwa, Chiba 277-8583, Japan</aff>

  <aff id="a5">Asia Pacific Center for Theoretical Physics,<break/>77 Cheongam-Ro, Nam-Gu, Pohang 790-784, Korea</aff>

  <aff id="a6">Department of Physics, Postech,<break/>Pohang 790-784, Korea</aff>
</contrib-group>
<pub-date><day>07</day><month>05</month><year>2015</year></pub-date>
<volume>2015</volume>
<issue>05</issue>
<fpage>040</fpage>
<history>
  <date date-type="received"><day>12</day><month>11</month><year>2014</year></date>
  <date date-type="corrected"><day>24</day><month>03</month><year>2015</year></date>
  <date date-type="accepted"><day>31</day><month>03</month><year>2015</year></date>
</history>
<permissions><copyright-statement>OPEN ACCESS, © The Authors</copyright-statement>
<copyright-year>2015</copyright-year><license license-type="cc-by" xlink:href="http://creativecommons.org/licenses/by/4.0/">
        <license-p>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.</license-p>
    </license>
</permissions>
<related-article related-article-type="arxiv"><pub-id pub-id-type="arxiv">1411.0664</pub-id></related-article>
<abstract>
<p>The search for light stops is of paramount importance, both in general as a promising
path to the discovery of beyond the standard model physics and more specifically as
a way of evaluating the success of the naturalness paradigm. While the LHC experiments
have ruled out much of the relevant parameter space, there are “stop gaps”, i.e.,
values of sparticle masses for which existing LHC analyses have relatively little
sensitivity to light stops. We point out that techniques involving on-shell constrained
<inline-formula><tex-math><![CDATA[$M_2$]]></tex-math></inline-formula>
variables can do much to enhance sensitivity in this region and hence help close the
stop gaps. We demonstrate the use of these variables for several benchmark
points and describe the effect of realistic complications, such as detector
effects and combinatorial backgrounds, in order to provide a useful toolkit for
light stop searches in particular, and new physics searches at the LHC in
general.
</p>
</abstract>
<kwd-group>
  <kwd>Supersymmetry Phenomenology</kwd>
  <kwd>Hadronic Colliders</kwd>
</kwd-group>
<funding-group>
<open-access>
<p content-type="scoap3">Article funded by SCOAP3</p></open-access></funding-group>
</article-meta>
</front><body>


<sec><title>Introduction</title>
<p><![CDATA[
\label{sec:introduction}




Supersymmetry (SUSY)~\cite{Martin-ml-1997ns} is an important framework
for beyond the Standard Model physics, as, among other features, it
provides an explanation of the relative lightness of the recently
discovered Higgs boson~\cite{Aad-ml-2012tfa, Chatrchyan-ml-2012ufa} and,
potentially, an explanation of dark matter.  The search for SUSY at
the CERN Large Hadron Collider (LHC), is therefore of considerable
importance.  The generic SUSY discovery channel is missing transverse
energy (MET) accompanied by hard jets, which results from the
production of gluinos and/or first generation squarks via the strong
interaction, followed by the subsequent decay of these sparticles to
final states that include an undetected lightest SUSY particle (LSP),
which is a dark matter candidate.
However, the LHC has yet to find evidence of such a
signal~\cite{Aad-ml-2013wta, Aad-ml-2014pda, Aad-ml-2014wea, Aad-ml-2014bva,
Chatrchyan-ml-2013wxa, Chatrchyan-ml-2013fea, Chatrchyan-ml-2014lfa,
Craig-ml-2013cxa}, which strongly motivates looking for SUSY elsewhere.  
Another class of
strong production processes, with somewhat lower cross sections, 
involves the production of third generation squarks: the stop and the sbottom.


]]></p>
<sec><title>Motivation for stop searches</title>
<p><![CDATA[
\label{subsec:motivation_for_stop_searches}

Stop production has long been recognized as a viable SUSY discovery
channel~\cite{Altarelli-ml-1984ve,Bigi-ml-1985aq,Hikasa-ml-1987db,Baer-ml-1991cb,Baer-ml-1994xr,Lopez-ml-1994zw,Chou-ml-1999zb,Demina-ml-1999ty}
 and has been looked for at LEP~\cite{ADLO-ml-2004aa, Heister-ml-2002hp,
   Abdallah-ml-2003xe, Achard-ml-2003ge, Abbiendi-ml-2002mp} and the
Tevatron~\cite{Jaffre-ml-2012gx, Aaltonen-ml-2010uf, Mackin-ml-2010zza}, as
well as at the LHC~\cite{Aad-ml-2014bva, Chatrchyan-ml-2013fea, Aad-ml-2013ija, Aad-ml-2014qaa,
Aad-ml-2014mha, Aad-ml-2014kva, Aad-ml-2014kra, Aad-ml-2014lra, Aad-ml-2014nra,
Chatrchyan-ml-2012paa, Chatrchyan-ml-2013lya, Chatrchyan-ml-2013xsw,
Chatrchyan-ml-2013xna, Chatrchyan-ml-2013iqa, Chatrchyan-ml-2013mya,
Khachatryan-ml-2014doa}.  Searches for the stop are especially
well-motivated theoretically because in many models, the stop is
expected to be the lightest squark for three principal reasons:

\begin{enumerate}
\item The beta function for a squark mass contains a \emph{positive}
  term proportional to the corresponding Yukawa coupling. The effect of 
  such a positive term is to suppress the mass when evolved from a high energy scale.  
  Since the top Yukawa coupling is the largest Yukawa coupling, one generically
  expects the stop soft mass parameter to emerge as the lightest of the
  squark soft masses after RGE evolution from some high energy
  scale~\cite{Feng-ml-1998iq}.
  (At large $\tan\beta$, similar arguments will apply to the sbottom mass as well.)
\item The left-right off-diagonal mixing in the squark mass matrix
  reduces the smaller mass eigenvalue via level repulsion.  The
  smaller of the two eigenvalues is therefore reduced relative to the
  corresponding diagonal element.  The left-right mixing is an
  $\SU(2)$-breaking effect, proportional to the Higgs vacuum
  expectation value and hence to the Yukawa coupling. Thus the stop
  is, again, the squark that would be most affected by a mass-lowering effect. 
\item  The radiative corrections to the tree-level relation among
  $m_Z$ and the Higgs soft mass parameters, which sets the electroweak
  scale, are dominated by stop loops.  Hence a large stop mass would
  destabilize the hierarchy (this has become known as the little
  hierarchy problem).  
  The desire to avoid excessive fine-tuning of the electroweak scale
  has spurred interest in ``natural SUSY'' models
  in which the top squark is among the lightest particles in the
  spectrum~\cite{Dimopoulos-ml-1995mi, Pomarol-ml-1995xc, Cohen-ml-1996vb,
    Feng-ml-1998iq, Perelstein-ml-2007nx}.
\end{enumerate}

Thus in this paper we will consider strategies to discover stops.
Although we have in mind searching for the stop in an arbitrary SUSY
model, everything we will say applies equally to other BSM models with
``top partners".\footnote{Our discussion will also apply to situations where 
there are no new particles and instead the top quark itself 
undergoes a rare decay involving more than one invisible
particle~\cite{Kong-ml-2014jwa}.} 
We assume that all other colored sparticles besides the ``stop'' are
either sufficiently heavy to be ignored or nonexistent.  
We are especially interested in what can be done to 
extend sensitivity into regions of parameter space where existing LHC
searches have not had sufficient sensitivity to discover or rule out
the stop. Our approach is complementary to several recent analyses which
have targeted similarly difficult parameter space regions for stop 
discovery~\cite{Han-ml-2008gy, Carena-ml-2008mj, Perelstein-ml-2008zt,
  Plehn-ml-2010st, Bornhauser-ml-2010mw, Bi-ml-2011ha, He-ml-2011tp, Drees-ml-2012dd,
  Bai-ml-2012gs,Plehn-ml-2012pr, Alves-ml-2012ft, Han-ml-2012fw, Kaplan-ml-2012gd,
  Berger-ml-2012an, Ghosh-ml-2012ud, Chen-ml-2012uw, Kilic-ml-2012kw,
  Graesser-ml-2012qy, Krizka-ml-2012ah, Delgado-ml-2012eu, Dutta-ml-2013sta,
  Buckley-ml-2013lpa, Chakraborty-ml-2013moa, Low-ml-2013aza, Bai-ml-2013ema,
  Belanger-ml-2013oka, Dutta-ml-2013gga, Buckley-ml-2014fqa, Ortiz-ml-2014iza,
  Czakon-ml-2014fka, Ismail-ml-2014cma, Eifert-ml-2014kea}.


]]></p>
</sec>
<sec><title>Stop decays</title>
<p><![CDATA[
\label{subsec:stop_decays}

Having identified and motivated the production mode that we will
consider, we now must decide on the manner in which the particles will
decay.
Unlike other squarks, for which only the gauge couplings are non-negligible
(thereby reducing the number of potentially relevant decay modes), stops have many 
viable decay modes --- there exist two-body decays of stops
to gluinos ($\tilde t \to \tilde g + t$), neutralinos ($\tilde t \to \tilde \chi^0_i + t$), 
charginos ($\tilde t \to \tilde \chi^+_i + b$), gravitinos ($\tilde t
\to \tilde G + t$)~\cite{Chou-ml-1999zb}, 
or even other stops ($\tilde t_2 \to \tilde t_1 +Z$)~\cite{Perelstein-ml-2007nx}
and sbottoms ($\tilde t \to \tilde b + W^+$)~\cite{Datta-ml-2011ef}. 
When the two-body decays are suppressed, there are several 
three-body decays which may dominate, e.g.~$\tilde t \to b W^+ \tilde \chi^0_i$,
$\tilde t \to b W^+ \tilde G$, $\tilde t \to b \ell^+ \tilde \nu$,
$\tilde t \to b \nu \tilde \ell^+$, etc.~\cite{Hikasa-ml-1987db,
  Baer-ml-1991cb, Chou-ml-1999zb, Demina-ml-1999ty}. Finally, 
there can also be loop-induced two-body decays, e.g.~$\tilde t \to c
\tilde\chi^0_i$~\cite{Hikasa-ml-1987db, Baer-ml-1991cb}.

\begin{figure}[t]
\centering
\includegraphics[trim = 0 0 0 0, clip, scale=0.9]{SignalUnderConsider.png}

]]></p>
<fig-group><caption><p><![CDATA[The two signal decay topologies under consideration in this
  paper.  In diagram (a), a stop (antistop) decays to a bottom
  (antibottom) quark and a positively (negatively) charged on-shell
  chargino; the chargino decays into an antilepton (lepton) and
  sneutrino (antisneutrino).  In diagram (b), the stop (antistop)
  decays to a top (antitop) quark and a neutralino.  The top (antitop)
  in turns decays to a bottom (antibottom) quark and a positively
  (negatively) on-shell charged $W^\pm$ boson, which decays, in turn,
  to an antilepton (lepton) and neutrino (antineutrino).  We refer to
  the topology in diagram (a) as ``Topology 1'' and the topology in
  diagram (b) as ``Topology 2''.
]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}




In this paper we shall focus on the most challenging scenario, in
which the stop produces\footnote{As usual, we assume that the 
stop decay proceeds through a ``decay chain'' of sequential 
decays of on-shell intermediate particles.} the same visible 
particles as a top quark decaying leptonically.  In particular, we
shall consider the two signal decay topologies shown in figure~\ref{fig:process}, which commonly occur in realistic 
models.\footnote{In principle, in addition to the two examples
from figure~\ref{fig:process}, there are many other decay topologies
which can mimic a top decay --- for example, there can be additional 
invisible particles emitted in this process~\cite{Agashe-ml-2010gt,
  Agashe-ml-2010tu, Giudice-ml-2011ib, Cho-ml-2012er}, 
or the neutralino in figure~\ref{fig:process}(b) can be emitted in
between the bottom quark and the lepton.}
In the process of figure~\ref{fig:process}(a), which we refer to as
``Topology 1", the stop decay chain is identical to the ``leptonic''
decay of a top quark; the only differences are that here the role of
the $W^\pm$ is played by the chargino $\tilde \chi^\pm_1$, and the
role of the neutrino is played by the sneutrino $\tilde \nu$.
(The sneutrino may further decay invisibly to another DM candidate;
since the sneutrino is on-shell, this does not affect our analysis.)
In the process of figure~\ref{fig:process}(b), which we shall refer to
as ``Topology 2", the stop decays to a top quark and a neutralino
first; the top then decays leptonically.  The resulting \emph{visible}
final state is the same, the difference now is that there are {\em
  two} invisible particles --- a neutralino $\tilde \chi^0_1$ and a
neutrino $\nu$.  When studying Topology 2, we shall assume that the
mass splitting between the stop and the neutralino is large enough
that the top quark produced in this decay is on-shell, as this makes
it \emph{more} difficult to distinguish the signal from top backgrounds.

\begin{figure}[t]
\centering
\includegraphics[trim = 0 0 0 0, clip, scale=0.95]{DecayTopologies.png}

]]></p>
<fig-group><caption><p><![CDATA[The four event topologies considered in this paper.  Diagram
  (a) gives the background event topology, where the top quarks decay
  leptonically, i.e.,~as in figure~\ref{fig:process}(a), but with the
  charginos replaced by $W^\pm$ bosons and the sneutrinos replaced by
  neutrinos. Diagram (b) is the (signal) event topology when both
  stops decay according to decay Topology 1, diagram (c) is the event
  topology when both stops decay following Topology 2, and diagram (d)
  is the mixed event topology when one stop decays according to
  Topology 1 and the other according to Topology
  2. ]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}





Specifically, we will consider how to discriminate between stop
production, where both stops decay according to one of the topologies
in figure~\ref{fig:process}, and the irreducible background from
$t\bar{t}$ dilepton events. Figure~\ref{fig:DecayTopologies}(a)
illustrates the background event topology, while the corresponding
three possible signal event topologies are depicted in figures~\ref{fig:DecayTopologies}(b-d).  In all these processes, the observed
final state consists of two $b$-jets, two opposite sign (OS) leptons,
and MET, which makes it quite challenging to discover stops in this channel.

\sloppy Note in particular that our analysis will allow for the mixed event
topology of \mbox{figure}~\ref{fig:DecayTopologies}(d).  This is because we shall 
\emph{not} make any assumptions about the relative branching
fraction between Topologies 1 and 2 in figure~\ref{fig:process}.  If the
two branching fractions are comparable, there is a sizable fraction of
events of the type depicted in figure~\ref{fig:DecayTopologies}(d);
their number benefits also from the combinatorial factor of 2 relative
to the events in figure~\ref{fig:DecayTopologies}(b) or the events in figure~\ref{fig:DecayTopologies}(c).


\fussy

\begin{figure}[t]
\centering
\includegraphics[scale=1]{DecaySubsystem.png}

]]></p>
<fig-group><caption><p><![CDATA[The background event
topology of figure~\ref{fig:DecayTopologies}(a)
with the corresponding subsystems explicitly delineated. The blue
dotted, green dot-dashed, and black solid lines indicate the
subsystems $(a)$, $(b)$, and $(ab)$, respectively.  See also table~\ref{tab:background-subsystems}.]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}

\boldmath

]]></p>
</sec>
<sec><title>On-shell constrained \texorpdfstring{$M_2$}{M(2)} variables</title>
<p><![CDATA[
\label{subsec:on_shell_constrained_M2_variables}
\unboldmath


In this paper we investigate the benefit of the recently proposed
on-shell constrained $M_2$
variables~\cite{Barr-ml-2011xt,Mahbubani-ml-2012kx,Cho-ml-2014naa} 
in discriminating between the signal events of figure~\ref{fig:DecayTopologies}(b-d) and the main background shown
in figure~\ref{fig:DecayTopologies}(a).  The $M_2$ variables are the natural
$3+1$-dimensional generalizations~\cite{Barr-ml-2011xt} of the Cambridge 
$M_{T2}$ variable~\cite{Lester-ml-1999tx,Barr-ml-2003rg}, which is already known
as a useful tool for background
suppression~\cite{Barr-ml-2009wu,Allanach-ml-2011ej,Murayama-ml-2011hj}.
Both $M_2$ and $M_{T2}$ were designed
for events in which particles are pair produced and decay
semi-invisibly.  (I.e., some of the decay products are invisible.)
The variables can then be computed for different subsystems in the event, 
or for the original event as a whole~\cite{Burns-ml-2008va}. 
For example, in the case of the background $t\bar{t}$ events 
from figure~\ref{fig:DecayTopologies}(a), there are three possibilities, which
are shown in figure~\ref{fig:DecaySubsystem}.  
We shall follow the notation of~\cite{Cho-ml-2014naa} and label these
three possibilities as $(ab)$,~$(a)$,~and~$(b)$.  

In the spirit of many other kinematic variables such as
$M_{T2}$~\cite{Lester-ml-1999tx},
$M_{2C}$~\cite{Ross-ml-2007rm}, and $M_{CT2}$~\cite{Cho-ml-2009ve},
the $M_2$ variables are obtained by minimizing some parent invariant mass
with respect to the momenta of the invisible daughter particles.  (The
exact definition and
basic properties of the constrained $M_2$ variables are reviewed in section~\ref{sec:review_of_on_shell_constrained_M2_variables} below.) In
the process of minimization, one may additionally impose certain
kinematic constraints which follow
from the hypothesized event topology. The main advantage of the $M_2$ class
of variables is that, being $3+1$-dimensional, they allow one to incorporate 
all known kinematic
constraints~\cite{Mahbubani-ml-2012kx,Cho-ml-2014naa}. For example, $M_{T2}$
and $M_{CT2}$ are transverse variables, and the only constraint which
can be used in their calculation
is the MET constraint. On the other hand, in calculating an $M_2$ variable,
one is free to impose additional mass shell conditions following either from a
previous measurement (as in the case of $M_{2C}$) or from a
theoretical hypothesis
about the specific nature of the events, e.g.\ that the two decay chains in figure~\ref{fig:DecaySubsystem} are the same and thus $M_{A_1}=M_{A_2}$,
$M_{B_1}=M_{B_2}$ and $M_{C_1}=M_{C_2}$. 

The presence of additional on-shell constraints
will generally increase the calculated value of $M_2$ --- the more constraints 
there are, the larger the value of $M_2$. This simple observation will be at the
center of our discussion below, and we shall use it in several different ways:
\begin{itemize}
\item The imposition of the additional constraints raises the value of $M_2$, 
distorting the shape of the $M_2$ distribution by increasing the 
populations of the higher $M_2$ bins. We can use this effect to increase 
the signal efficiency in those cases where the background $M_2$ distribution is
bounded by an upper kinematic endpoint, while the signal $M_2$ distribution extends 
beyond this kinematic endpoint. This effect will be discussed and illustrated 
in section~\ref{sec:M2_endpoint_study_for_topology_1} for signal events 
with Topology 1 only, and in section~\ref{sec:M2_endpoint_study_for_topology_2} 
for signal events with Topology 2 only. In either case, the two top squarks decay identically, 
giving rise to the ``symmetric"  event topologies of figure~\ref{fig:DecayTopologies}(b) 
and figure~\ref{fig:DecayTopologies}(c), respectively.
\item The kinematic variables $M_{T2}$ and $M_2$ were originally designed 
for symmetric events, and intended to be used for such events only.
But what if the actual event is ``asymmetric" and the two parent particles 
decay in a different manner, e.g.,~as in figure~\ref{fig:DecayTopologies}(d)?
There are two possible approaches. First, one could suitably modify the
definition of $M_{T2}$ in order to adapt it to an asymmetric
case~\cite{Barr-ml-2009jv,Konar-ml-2009qr}.
This would still work, provided that \emph{both} the signal and background have the same 
asymmetric event topology. However, what if the background events are symmetric (as in figure~\ref{fig:DecayTopologies}(a)),
while the signal events are asymmetric (as in figure~\ref{fig:DecayTopologies}(d)) 
or vice versa? This case is the subject of section~\ref{sec:mixed-events}, in which we shall advocate the use of the
conventional ``symmetric" 
on-shell constrained $M_2$ variables for this asymmetric case as well. 
As an illustration, we shall consider a particularly difficult scenario, when
the kinematic endpoints of the symmetric events with Topology 1
(figure~\ref{fig:DecayTopologies}(b))
or Topology 2 (figure~\ref{fig:DecayTopologies}(c)) are too low and are
both ``buried" inside the background distribution. Nevertheless, when
we compute the $M_2$ variables for the mixed event topology of figure~\ref{fig:DecayTopologies}(d), we shall find that the signal
distribution does extend beyond the background endpoint.  The reason
for this apparent ``endpoint violation" is simply the fact that when
calculating $M_2$, we are applying the ``wrong" constraints on the
signal events, but the ``correct" constraints on the background events.
\item The benefit of the on-shell constrained $M_2$ variables is not
  limited only to events in which the background $M_2$ endpoint is
  violated. We can achieve additional separation of signal from
  background by studying the size of the shift caused by the
  application of the on-shell constraints. In section~\ref{sec:detector_resolution_etc} we shall find that this shift
  is generally larger for signal events with the mixed event topology
  of figure~\ref{fig:DecayTopologies}(d) when compared to the
  corresponding shift for background events.  This observation is
  valid even for signal events which do not violate the background
  $M_2$ endpoint.  The study presented in section~\ref{sec:detector_resolution_etc} also includes the effects
  of detector resolution and combinatorial backgrounds.
\end{itemize}

\pagebreak


]]></p>
</sec>
<sec><title>Pr\'{e}cis</title>
<p><![CDATA[
\label{subsec:precis}

In this paper we put forth four ideas for stop discovery.
\begin{enumerate}
\item We propose to use the on-shell constrained $M_2$ variables from
  \emph{each} of the three subsystems ($(a)$, $(b)$, and
  $(ab)$)~\cite{Burns-ml-2008va,Bai-ml-2013ema}.  Previous efforts described
  in the literature have relied mostly on the dileptonic
  $M_{T2}$~\cite{Plehn-ml-2012pr,Kilic-ml-2012kw,Chakraborty-ml-2013moa,Cho-ml-2008cu}.
\item We advocate the use of the on-shell constrained
  $(3+1)$-dimensional $M_2$ variables in place of their transverse
  cousin $M_{T2}$, due to their ability to ``push'' more signal events
  beyond the background endpoints, thus increasing the signal
  efficiency.\footnote{The variable $M_{T2}^W$ introduced in~\cite{Bai-ml-2012gs} is a concrete realization of an on-shell
    constrained $M_2$ variable in the case when one of the leptons in figure~\ref{fig:DecayTopologies} is lost.}
\item We also propose to use the \emph{difference} between the on-shell
  constrained $M_2$ variable and its analogue $M_{T2}$ as an
  additional discriminator against the background.
\item We show that the on-shell constrained $M_2$ variable is able to
  specifically target the mixed signal event topology of figure~\ref{fig:DecayTopologies}(d) and salvage a certain fraction of
  signal events in difficult scenarios when more conventional cuts
  would fail.
\end{enumerate}


]]></p>
</sec>
</sec>
<sec><title>Review of on-shell constrained \texorpdfstring{\boldmath $M_2$}{M(2)} variables</title>
<p><![CDATA[
\label{sec:review_of_on_shell_constrained_M2_variables}

In this section we provide a brief review of on-shell constrained $M_2$
variables. Readers who are familiar with the terminology and notation 
of ref.~\cite{Cho-ml-2014naa} may skip directly to section~\ref{sec:M2_endpoint_study_for_topology_1}.
We begin by reminding the reader that there exists a broad class of kinematic
variables that are useful in the analysis of events with missing energy. 
These variables are defined in two steps~\cite{Barr-ml-2011xt}:
\begin{itemize}
  \item One first assumes a particular event topology consistent with the
    final state particles observed in the event.
  \item One then minimizes an invariant mass quantity in the
    hypothesized topology over the
    unknown invisible momenta, subject to certain kinematic constraints.
\end{itemize}
The best known example of such a variable is the usual transverse mass $M_T$,
~\cite{Smith-ml-1983aa,Barger-ml-1983wf}, which applies to the simple case
of one decay chain with
a single invisible particle. Although a transverse quantity, $M_T$ can be 
thought of as the minimum value of the 3+1 dimensional invariant mass,
consistent with the measured missing transverse
momentum~\cite{Barr-ml-2011xt}.
This example is rather trivial in the sense that
imposing the vector constraint of transverse momentum conservation already 
fixes the (transverse) components of the invisible particle momentum,
and there is only one minimization left to do. A much more interesting
case arises
when there are two decay chains, with one invisible particle in each. Then the 
concept of $M_T$ is generalized to the stransverse mass,
$M_{T2}$~\cite{Lester-ml-1999tx},
in which one finds the minimum, with respect to the
momenta of invisible particles, of the maximum transverse mass of a
given particle on either side of a (symmetric) decay topology, again 
subject to the constraint of total transverse momentum conservation.

The on-shell constrained $M_2$ variables described in~\cite{Cho-ml-2014naa}  
are analogous to $M_{T2}$, with two main differences. First, the quantity 
being minimized is the four-dimensional invariant mass
rather than the transverse mass (see also~\cite{Ross-ml-2007rm,
  Barr-ml-2011xt, Mahbubani-ml-2012kx}). 
Second, in addition to transverse momentum conservation, one is free
to apply additional
on-shell constraints which follow from the assumed event topology.
For concreteness, let us use the event topology of figure~\ref{fig:DecaySubsystem}
to illustrate the procedure of defining the different types of $M_2$
variables, denoted as
\beq
M_{2\sqcup\sqcup}(S;\tilde m).
\label{M2notation}
\eeq
Here $S\in \{ (ab), (a), (b) \}$ denotes the subsystem under consideration, while
$\sqcup$ is an index placeholder to be defined shortly.
\begin{table}[t]
\centering
\renewcommand{\arraystretch}{1.2}
\begin{tabular}{| c | c | c | c |}
\hline
Subsystem & Parents & Daughters & Relatives  \\ 
$S$            &  $P_i$    &  $D_i$       & $R_i$  \\
\hline 
$(ab)$  &   $\{A_i\}=\{t,\bar{t}\,\}$   &    $\{C_i\}=\{\nu,\bar{\nu}\, \}$  &  $\{B_i\}=\{W^+, W^-\}$ \\
$(a)$    &   $\{A_i\}=\{t,\bar{t}\,\}$   &    $\{B_i\}=\{W^+,W^-\}$         &   $\{C_i\}=\{\nu,\bar{\nu}\,\}$ \\
$(b)$    &   $\{B_i\}=\{W^+,W^-\}$  &    $\{C_i\}=\{\nu,\bar{\nu}\,\}$ &     $\{A_i\}=\{t,\bar{t}\,\}$ \\
\hline
\end{tabular}

]]></p>
<table-wrap><caption><p><![CDATA[The roles played by
  different particles in the background process of $t\bar{t}$ production
   (\mbox{figure}~\ref{fig:DecayTopologies}(a)), for each of the three
   subsystems defined in \mbox{figure}~\ref{fig:DecaySubsystem}. ]]></p></caption></table-wrap>
<p><![CDATA[
\end{table}
According to the nomenclature of~\cite{Cho-ml-2014naa}, depending on the 
subsystem being considered,
the intermediate particles $A_i$, $B_i$, and $C_i$ fall into one of the following 
three categories (see table~\ref{tab:background-subsystems}):
\begin{itemize}
\item \emph{Daughters} $D_i$. These are the invisible particles at the end of
the decay chains in the subsystem under consideration. 
Following~\cite{Cho-ml-2014naa}, we shall denote their 3-momenta by
$\vec{q}_1$ and $\vec{q}_2$,
respectively. The value of the $M_2$ variable~(\ref{M2notation}) will
be obtained by
minimizing over all possible values of $\vec{q}_1$ and $\vec{q}_2$,
consistent with the applied kinematic constraints.
As usual, we shall take the daughters' masses to be equal:
\begin{equation}
\label{eq:daughters-equal}
M_{D_1} = M_{D_2} \equiv \tilde m
\end{equation}
and denote them with $\tilde m$, which will be an input parameter for
the $M_2$ calculation.
\item \emph{Parents} $P_i$. These are the two particles at the top of
  the decay chains in the subsystem,
and their masses will be subject to minimization over the invisible momenta in 
order to obtain the variable~(\ref{M2notation}). When performing this
minimization, in addition to the missing transverse momentum constraint
\beq
\vec{q}_{1T}+\vec{q}_{2T} = \mpt, 
\eeq
we can additionally require that the two parent masses 
(when considered as functions of the invisible momenta) are the same:
\begin{equation}
\label{eq:parents-equal}
M^2_{P_1}(\vec{q}_1,\tilde m) = M^2_{P_2}(\vec{q}_2,\tilde m).
\end{equation}
The presence (or absence) of this constraint is indicated by the first
${}_\sqcup$ index in~(\ref{M2notation}),
which takes value $C$ if the constraint is applied, and $X$ otherwise:
\begin{align}
M_{2C\sqcup}(S;\tilde m)  &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\ 
\vec{q}_{1T}+\vec{q}_{2T} = \mpt   \\
M^2_{P_1}(\vec{q}_1,\tilde m) = M^2_{P_2}(\vec{q}_2,\tilde m) } }
\left\{M_{P_1}(\vec{q}_{1},\tilde m) \right\},  
\label{M2Cbla} \\
M_{2X\sqcup}(S;\tilde m)  &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\ 
\vec{q}_{1T}+\vec{q}_{2T} = \mpt} }
\left\{ \max\left[ M_{P_1}(\vec{q}_{1},\tilde m),
    M_{P_2}(\vec{q}_{2},\tilde m) \right] \right\}.
\label{M2Xbla} 
\end{align}
\item \emph{Relatives} $R_i$. As shown in table~\ref{tab:background-subsystems}, the relatives are the
  remaining particles in the event topology --- they are neither
  parents nor daughters. Depending on the subsystem, relatives can
  appear either inside or outside the subsystem.  Their masses can
  also be written as functions of the respective daughters' momenta,
  so by requiring equal masses\footnote{We note that while a parent
    mass squared is always positive, the mass squared of a {\em
      relative} could be negative: keep in mind that the values for
    the invisible momenta $\vec{q}_1$ and $\vec{q}_2$ found in the
    minimization process are not the true momenta and could be
    unphysical, i.e., there is no guarantee that $M^2_{R_i}>0$.  While
    one has the option of adding the further constraint that the
    squared masses of the relative particles be positive, we will not
    do so in this work.} for the relative particles,
\begin{equation}
\label{eq:relatives-equal}
M_{R_1}^2 (\vec{q}_1,\tilde m)= M_{R_2}^2(\vec{q}_2, \tilde m),
\end{equation} 
we are, in effect, imposing an additional constraint on the
minimization over the invisible momenta $\vec{q}_i$. The applicability
of the constraint~(\ref{eq:relatives-equal}) will be indicated by the
second ${}_\sqcup$ index in~(\ref{M2notation}): as before, it will be
equal to $C$ if the constraint~(\ref{eq:relatives-equal}) is applied
and $X$ otherwise.
Altogether, therefore, we have four possible $M_2$ variables:
\begin{align}
M_{2CC}(S;\tilde m)  &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\ 
\vec{q}_{1T}+\vec{q}_{2T} = \mpt   \\
M_{P_1}^2(\vec{q}_1,\tilde m) = M_{P_2}^2(\vec{q}_2,\tilde m) \\
M_{R_1}^2(\vec{q}_1,\tilde m) = M_{R_2}^2(\vec{q}_2,\tilde m)} }
\left\{M_{P_1}(\vec{q}_{1},\tilde m) \right\},  
\label{M2CCdef} \\ M_{2CX}(S;\tilde m)  &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\ 
\vec{q}_{1T}+\vec{q}_{2T} = \mpt \\
M_{P_1}^2(\vec{q}_1,\tilde m) = M_{P_2}^2(\vec{q}_2,\tilde m)} }
\left\{M_{P_1}(\vec{q}_{1},\tilde m) \right\},  
\label{M2CXdef} \\ M_{2XC}(S;\tilde m)  &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\ 
\vec{q}_{1T}+\vec{q}_{2T} = \mpt \\
M_{R_1}^2(\vec{q}_1,\tilde m) = M_{R_2}^2(\vec{q}_2,\tilde m)} }
\left\{ \max\left[ M_{P_1}(\vec{q}_{1},\tilde m), M_{P_2}(\vec{q}_{2},\tilde m) \right] \right\},
\label{M2XCdef} \\% [2mm]
M_{2XX}(S;\tilde m)  &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\ 
\vec{q}_{1T}+\vec{q}_{2T} = \mpt } }
\left\{ \max\left[ M_{P_1}(\vec{q}_{1},\tilde m), M_{P_2}(\vec{q}_{2},\tilde m) \right] \right\}.
\label{M2XXdef} 
\end{align}
\end{itemize}



The definitions~(\ref{M2CCdef})--(\ref{M2XXdef}) should be contrasted to
the analogous definition of the Cambridge $M_{T2}$ variable
\beq
M_{T2}(S;\tilde m)  \equiv
\min_{\substack{\vec{q}_{1T},\vec{q}_{2T}\\ 
\vec{q}_{1T}+\vec{q}_{2T} = \mpt  } }
\left\{ \max\left[ M_{TP_1}(\vec{q}_{1T},\tilde m),
    M_{TP_2}(\vec{q}_{2T},\tilde m) \right] \right\},
\label{MT2def} 
\eeq
where only the transverse components $\vec{q}_{iT}$ are used, and the 
objective function (the function that is minimized) is the larger of
the two \emph{transverse} masses of the parents.




The main goal of this paper is to investigate and contrast the ability of the
variables~(\ref{M2CCdef})--(\ref{MT2def}) to discriminate between the $t\bar{t}$
background of figure~\ref{fig:DecayTopologies}(a) and the three types of
stop signal events of figure~\ref{fig:DecayTopologies}(b-d). We shall
consider the respective variables for all three subsystems of figure~\ref{fig:DecaySubsystem}. Since we know the mass spectrum for the
background event topology, we shall choose the test mass $\tilde m$ to
be equal to the correct, SM value for the respective daughter
particle. In particular,
\begin{itemize}
\item In subsystem $(ab)$, the parent particle is the top quark;
the daughter particle is the neutrino, whose mass is taken to vanish
($\tilde m=0$).  The other, ``relative'', particle is the $W^\pm$ boson. 
Then \emph{for background events}, all 5 variables~(\ref{M2CCdef})--(\ref{MT2def}) 
are bounded from above by the top mass:
\beq
M_{T2}(ab;0), 
M_{2XX}(ab;0),
M_{2CX}(ab;0),
M_{2XC}(ab;0),
M_{2CC}(ab;0)
\le m_t,
\label{mtab}
\eeq
while for signal events, this bound can be violated.
\item In subsystem $(a)$, the parent particle is again the top quark.
The daughter particle is now the $W^\pm$ boson, with mass $\tilde m= m_W$.
The relative particles are the two neutrinos, which in this case appear
downstream outside the subsystem. For background events, the 
variables are again bounded by the top mass
\beq
M_{T2}(a;m_W), 
M_{2XX}(a;m_W),
M_{2CX}(a;m_W),
M_{2XC}(a;m_W),
M_{2CC}(a;m_W)
\le m_t.
\label{mta}
\eeq
\item In subsystem $(b)$ the parent particles are the $W^\pm$ bosons, 
the daughter particles are the two neutrinos with mass $\tilde m=0$,
and the relative particles are the top quarks appearing upstream
outside the subsystem. The background events obey
\beq
M_{T2}(b;0), 
M_{2XX}(b;0),
M_{2CX}(b;0),
M_{2XC}(b;0),
M_{2CC}(b;0)
\le m_W.
\label{mwb}
\eeq
\end{itemize}

In principle, \emph{each} of the bounds~(\ref{mtab})--(\ref{mwb}) allows
us to cut $100\%$ of the
$t\bar{t}$ background events by removing events with values of the
respective subsystem $M_2$ or $M_{T2}$ variable below the appropriate
threshold (``high pass cut''). 
Therefore, as far as just the background is concerned,
we have 15 alternative choices\footnote{Each of the five variables~(\ref{M2CCdef})--(\ref{MT2def}) can be applied for each of the three
  subsystems $(ab)$, $(a)$, and~$(b)$. The explicit definitions of the
  resulting 15 variables can be found in appendix~\ref{sec:appendix}.}
for reducing it, and they should perform comparably well.
The differences between the five variables~(\ref{M2CCdef})--(\ref{MT2def}) begin to emerge when we consider the
effect of such a high pass cut on signal events.  It has been
shown~\cite{Cho-ml-2014naa} that the variables~(\ref{M2CCdef})--(\ref{MT2def}) obey the following hierarchy
\begin{equation}
    \label{eq:inequality}
  M_{T2}(S; \tilde{m}) =  M_{2XX}(S; \tilde{m}) = M_{2CX}(S;
  \tilde{m})  \leq M_{2XC}(S; \tilde{m}) \leq M_{2CC}(S; \tilde{m})
\end{equation} 
for any subsystem $S$ and any value of $\tilde m$.
We should therefore expect the distributions of $M_{2XC}$ and $M_{2CC}$ 
to be more populated at higher values and in particular near their endpoints.
As a result, signal events are more likely to pass if the cut is applied on 
the additionally constrained variables, $M_{2XC}$ and $M_{2CC}$,
as opposed to the less constrained variables $M_{T2}$, $M_{2XX}$, and $M_{2CX}$.
This expectation is borne out by the explicit studies below.  
The property~(\ref{eq:inequality}) offers an opportunity to increase 
the sensitivity of the LHC experiments to the presence of stop signals of the 
type described in figure~\ref{fig:DecayTopologies}(b-d).


]]></p>
</sec>
<sec><title>\texorpdfstring{\boldmath $M_2$}{M(2)} endpoint study for topology 1</title>
<p><![CDATA[
\label{sec:M2_endpoint_study_for_topology_1}

We begin by studying the effectiveness of the on-shell constrained
subsystem $M_2$ variables in the case where both stops in the signal
event decay according to Topology 1. This yields the event topology
pictured in figure~\ref{fig:DecayTopologies}(b). The background, as always in
this paper, consists of dileptonic top production and is shown in figure~\ref{fig:DecayTopologies}(a).
We remind the reader that when applied to \emph{background events},
all of the on-shell constrained subsystem $M_2$ variables
(as well as $M_{T2}$) exhibit very well-defined kinematic endpoints given
by eqs.~(\ref{mtab})--(\ref{mwb}). Therefore, the effectiveness of the $M_2$ 
variables in identifying signal events is determined by how many signal events
violate the bounds~(\ref{mtab})--(\ref{mwb}).

The signal events considered in this section (those of figure~\ref{fig:DecayTopologies}(b)) 
have exactly the same topology as the background. Therefore, when applied to 
\emph{signal events}, the $M_2$ variables will have well-defined upper
kinematic endpoints as well. The precise value of those endpoints will
depend on the underlying signal mass spectrum, i.e., on the true
values of the stop mass $m_{\tilde t}$, 
the chargino mass $m_{\tilde\chi^\pm}$, and the sneutrino mass $m_{\tilde \nu}$.
For any given point $\{ m_{\tilde t}, m_{\tilde\chi^\pm}, m_{\tilde \nu}\}$ 
in the mass parameter space, using the formulas given in ref.~\cite{Burns-ml-2008va}, one can compute the expected $M_2$ kinematic
endpoints for the signal, in each 
of the three subsystems $(ab)$, $(a)$, and $(b)$.\footnote{We remind 
the reader that in this work, the $M_2$ variables are always computed
with test masses corresponding to the background hypothesis; see appendix~\ref{sec:appendix}.} 
Depending on the SUSY mass spectrum, some, all, or none of these
signal endpoints will exceed the corresponding background
endpoints. In the second half of this section, we shall illustrate
each of these three scenarios with specific study points. But first we
shall analyze the relevant mass parameter space and categorize the
different regions, which are defined by the location of the signal endpoints
relative to the background endpoints.  This will be the subject of the
next subsection.



]]></p>
<sec><title>Anatomy of the mass parameter space for Topology 1</title>
<p><![CDATA[
\label{subsec:regions-and-study-points}

In order to divide the stop-chargino-sneutrino mass parameter space
into regions, we start with the analytical expressions for the $M_{T2}$ signal 
endpoints~\cite{Cho-ml-2007qv,Cho-ml-2007dh,Burns-ml-2008va}. 
(The endpoints for the corresponding $M_{2\sqcup\sqcup}$ variables are
given by the exact same expressions~\cite{Cho-ml-2014naa}.)
\begin{align}
M_{T2}^{\rm max}(ab;\tilde m=0) & =   \frac{m^2_{\tilde t}-m^2_{\tilde
    \nu}}{m_{\tilde t}}, 
\label{susymtab}\\
M_{T2}^{\rm max}(a;\tilde m=m_W) & =  \frac{m^2_{\tilde t}-m^2_{\tilde
    \chi^\pm}}{2m_{\tilde t}}
+\sqrt{\left(\frac{m^2_{\tilde t}-m^2_{\tilde \chi^\pm}}{2m_{\tilde
        t}}\right)^2+m_W^2},
\label{susymta} \\
M_{T2}^{\rm max}(b;\tilde m=0) & = 
 \sqrt{ \frac{\left(m^2_{\tilde t}-m^2_{\tilde \nu}\right)\left(m^2_{\tilde
       \chi^\pm}-m^2_{\tilde \nu}\right)}{m^2_{\tilde t}}  }.
\label{susymwb}
\end{align}
Now combining, e.g., eqs.~(\ref{mtab}) and~(\ref{susymtab}), we find that 
for the variables in the $(ab)$ subsystem, the signal endpoints exceed the
background endpoints if
\beq
\frac{m^2_{\tilde t}-m^2_{\tilde \nu}}{m_{\tilde t}} > m_t.
\label{condition_ab}
\eeq
The corresponding region is delineated by the solid black line in figure~\ref{fig:map},
which shows a slice through the 3-dimensional mass parameter space 
for fixed $m_{\tilde \nu}=110$\,GeV\@. For convenience we choose to represent
the remaining two degrees of freedom as the mass differences 
$m_{\tilde\chi^\pm}-m_{\tilde \nu}$ and $m_{\tilde t}-m_{\tilde\chi^\pm}$.
The region satisfying the condition~(\ref{condition_ab}) is above and 
to the right of the solid black line in figure~\ref{fig:map}. If the SUSY mass 
spectrum happens to be in this region, the mass splitting between
the stop and the sneutrino is sufficiently large to cause some number of 
signal events to ``leak" beyond the background endpoint~(\ref{mtab}).
This means that the subsystem $(ab)$ variables are promising
variables to cut on in order to separate signal from background.

\begin{figure}[t]
\centering
\includegraphics[trim = 0 0 0 0, clip, scale=1]{figure4.png}

]]></p>
<fig-group><caption><p><![CDATA[The division of stop-chargino-sneutrino mass parameter space
  into regions, which are defined by which of the background endpoints
  are violated by signal events in which both stops decay according to
  Topology 1.  The sneutrino mass has been set to $110$\,GeV\@.  The four
  study points considered later in this section are also indicated.
]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}

We can apply similar reasoning to the invariant mass variables in subsystems 
$(a)$ and~$(b)$. For example, comparing eqs.~(\ref{mta}) and~(\ref{susymta}),
we find that the bound~(\ref{mta}) applicable to background events will be violated if
the mass spectrum is such that
\beq
\frac{m^2_{\tilde t}-m^2_{\tilde \chi^\pm}}{m_{\tilde t}} >
\frac{m^2_t-m^2_W}{m_t}.
\label{condition_a}
\eeq
The corresponding region extends above the red dashed line in figure~\ref{fig:map}.
Finally, the condition for violating the background endpoints in
subsystem $(b)$ follows from eqs.~(\ref{mwb}) and~(\ref{susymwb}):
\beq
\frac{\left(m^2_{\tilde t}-m^2_{\tilde \nu}\right)\left(m^2_{\tilde
    \chi^\pm}-m^2_{\tilde \nu}\right)}{m^2_{\tilde t}}  > m_W^2.
\label{condition_b}
\eeq
The region where this condition is satisfied is located to the right of the 
blue dot-dashed line in figure~\ref{fig:map}.

For completeness, we shall also consider the variable $m_{bl}$, 
the invariant mass of the lepton and $b$ quark from a given branch of
the decay.  
The endpoint of this quantity for signal events is given by
\beq
m_{b\ell}^{\rm max} = 
\sqrt{\frac{\left(m_{\tilde{t}}^2 -
    m_{\tilde{\chi}^\pm}^2\right)\left(m_{\tilde{\chi}^\pm}^2-
    m_{\tilde{\nu}}^2\right)}{m_{\tilde{\chi}^\pm}^2 }},
\label{susymbl}
\eeq
while background events will obey the bound
\beq
m_{b\ell} \le \sqrt{m_t^2 - M_W^2}.
\label{mblmaxbknd}
\eeq
Therefore the condition for violating the background $m_{b\ell}$ endpoint is
\beq
\frac{\left(m_{\tilde{t}}^2 -
  m_{\tilde{\chi}^\pm}^2\right)\left(m_{\tilde{\chi}^\pm}^2-
  m_{\tilde{\nu}}^2\right)}{m_{\tilde{\chi}^\pm}^2 }
> m_t^2 - M_W^2.
\label{condition_bl}
\eeq
The corresponding region is found above and to the right of the
diagonal magenta thin solid line in figure~\ref{fig:map}.

\begin{table}[t]
\centering
\begin{tabular}{| c | c | c | c | c |}
\hline
~\rule[13pt]{0pt}{0pt} & $m_{b\ell}$  & $M_{T2}(ab)$ & $M_{T2}(a)$ & $M_{T2}(b)$ \\
\textbf{Region} & endpoint  & endpoint & endpoint & endpoint \\
~ & violation  & violation & violation & violation \\ 
\hline 
\textcolor{region1}{i}  \!\!\!  \rule[13pt]{0pt}{0pt}      &  No   &  Yes & Yes & Yes \\
\textcolor{region2}{ii}         &  Yes  &  Yes & Yes & Yes \\
 \textcolor{region3}{iii}       &  Yes  &  Yes & Yes  & No \\
 \textcolor{region4}{iv}       &  No   &  Yes & Yes  & No \\
 \textcolor{region5}{v}        &  No   &  No & Yes  & No \\
 \textcolor{region6}{vi}       &  No   &  No  & No  & No \\
 \textcolor{region7}{vii}      &  No   &  No  & No & Yes \\
 \textcolor{region8}{viii}     &  No   &  Yes  & No & Yes \\
 \textcolor{region9}{ix}       &  Yes  &  Yes  & No & Yes \\
\hline
\end{tabular}

]]></p>
<table-wrap><caption><p><![CDATA[The regions of stop-chargino-sneutrino
  mass space, as depicted in figure~\ref{fig:map}.  The four columns
  describe whether in the given region, the specified background
  endpoint can be violated by signal events where both stops decay
  according to Topology 1.]]></p></caption></table-wrap>
<p><![CDATA[
\end{table}

The conditions implied by eqs.~(\ref{condition_ab})--(\ref{condition_b})
and~(\ref{condition_bl}) divide the mass parameter space of figure~\ref{fig:map}
into nine distinct color-coded regions, which are defined in table~\ref{tab:regions}.
  It is easy to verify analytically that the boundaries of the three
  regions defined by conditions~(\ref{condition_ab})--(\ref{condition_b}), i.e., the red, blue, and black curves in figure~\ref{fig:map}, cross at a single point. For any given sneutrino mass, 
$m_{\tilde \nu}$, the values of the stop mass, $m_{\tilde t}$, and the
chargino mass,
$m_{\tilde{\chi}^\pm}$, corresponding to the triple crossing point are
found from the relations
$$
\frac{m^2_{\tilde t}-m^2_{\tilde \nu}}{m_{\tilde t}} = m_t, \qquad
\frac{m_{\tilde{\chi}^\pm}^2- m_{\tilde{\nu}}^2}{m_{\tilde t}} = \frac{m_W^2}{m_t}. 
$$

The map in figure~\ref{fig:map} serves a dual purpose.  First, it
singles out the regions which might be easier to discover,
as well as the regions which may pose challenges. Second, within each
region, it identifies the variables which might be useful in the
analysis (see table~\ref{tab:regions}). For example, in region
\textcolor{region2}{ii}, the mass spectrum is sufficiently split, and
all four variables exhibit endpoint violations for signal events.
This in turn suggests that separating signal from background should be
relatively easy, since we can use any of the four types of invariant
mass variables in table~\ref{tab:regions} to suppress the background
without much signal loss. To illustrate the expected phenomenology of
region \textcolor{region2}{ii}, in section~\ref{subsec:point-1} we
shall analyze in detail a specific study point from this region. Its
mass spectrum is given in table~\ref{tab:study-points} and its exact
location on the map of figure~\ref{fig:map} is marked with a black cross
($+$).

\begin{table}[t]
\centering
\renewcommand{\arraystretch}{1.1}\begin{tabular}{| c | c | c | c | c |}
\hline
Study Point & Stop Mass & Chargino Mass & Sneutrino Mass & Region \\
\hline 
1  & $1110$\,GeV & $610$\,GeV & $110$\,GeV & \textcolor{region2}{ii} \\
2  & $215$\,GeV & $210$\,GeV & $110$\,GeV & \textcolor{region7}{vii} \\
3  & $630$\,GeV & $130$\,GeV & $110$\,GeV & \textcolor{region3}{iii} \\
4  & $174$\,GeV & $150$\,GeV & $110$\,GeV & \textcolor{region6}{vi} \\
\hline
\end{tabular}

]]></p>
<table-wrap><caption><p><![CDATA[The stop, chargino, and sneutrino masses
  for the four study points considered in this section, as well as the
  region of mass parameter space from figure~\ref{fig:map} that contains
  each point. ]]></p></caption></table-wrap>
<p><![CDATA[
\end{table}
 
In all but one of the remaining regions of figure~\ref{fig:map}, some of
the endpoints are violated while others are not, thus some variables
are expected to perform better than others.  We pick two
representative study points in regions \textcolor{region7}{vii} and
\textcolor{region3}{iii} and study them in sections~\ref{subsec:point-2} and~\ref{subsec:point-3}, respectively.
In figure~\ref{fig:map}, these two study points are marked with the
magenta asterisk $(\ast)$ and the blue circle ($\bullet$). The
corresponding mass spectra are also listed~in~table~\ref{tab:study-points}. 
 
Finally, region \textcolor{region6}{vi} deserves a special mention,
since it represents a particularly challenging scenario.  Here none of
the four types of invariant mass variables exhibits an endpoint
violation, and the signal events are populating the same kinematic
region as the background events.  Our fourth study point 
in table~\ref{tab:study-points} 
(denoted in figure~\ref{fig:map} with the red ($\times$) symbol)
belongs to this challenging region and is considered in section~\ref{subsec:point-4}.


]]></p>
</sec>
<sec><title>Study point 1: split spectrum in region ii</title>
<p><![CDATA[
\label{subsec:point-1}

In this section we shall illustrate the properties of region
\textcolor{region2}{ii} in figure~\ref{fig:map}
with the study point 1 which is marked with the black ($+$) symbol.
As seen in table~\ref{tab:study-points}, this study point 
has a widely split spectrum; both mass differences 
$m_{\tilde t}-m_{\tilde \chi^\pm}$ and $m_{\tilde \chi^\pm} -
m_{\tilde \nu}$ are $500$\,GeV.
\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{M2_Subsystem_ab_Topo2BP3}\quad 
\includegraphics[width=7.0cm]{M2_Subsystem_a_Topo2BP3}\\[3pt]
\includegraphics[width=7.0cm]{M2_Subsystem_b_Topo2BP3}\quad 
\includegraphics[width=7.0cm]{Mbl_Topo2BP3}

]]></p>
<fig-group><caption><p><![CDATA[Signal (thick red lines) and background (thin blue lines)
  distributions for the subsystem variables $M_{T2}(S;\tilde m)$
  (dashed lines) and $M_{2CC}(S;\tilde m)$ (solid lines). We show
  results for each of the three possible subsystems: $S=(ab)$ (upper
  left panel), $S=(a)$ (upper right panel), and $S=(b)$ (lower left
  panel). The lower right panel shows the analogous distributions for
  the invariant mass $m_{b\ell}$. The signal events are for study
  point 1 ($m_{\tilde t}=1110$\,GeV, $m_{\tilde \chi^\pm}=610$\,GeV and
  $m_{\tilde \nu}=110$\,GeV) and have the event topology shown in figure~\ref{fig:DecayTopologies}(b).]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
One can therefore expect that the signal distributions for our
invariant mass variables will extend well beyond the corresponding
background distributions.  This is confirmed in figure~\ref{fig:MT2cuts3}, where we plot distributions for signal events
(red lines) and background events (blue lines) for $M_{T2}(S;\tilde
m)$ (dashed lines), $M_{2CC}(S;\tilde m)$ (solid lines), and
$m_{b\ell}$.   As always in this section, the signal events are
assumed to have the symmetric event topology of figure~\ref{fig:DecayTopologies}(b), i.e., both stops decay according to
Topology 1.  Figure~\ref{fig:MT2cuts3} shows results for all three
subsystems: $S=(ab)$ (upper left panel), $S=(a)$ (upper right panel)
and $S=(b)$ (lower left panel).  In each case, the trial mass $\tilde
m$ for the respective daughter particle has been set to the correct
value for SM events, as explained in more detail in appendix~\ref{sec:appendix}.  In each panel of figure~\ref{fig:MT2cuts3}, the vertical dashed line marks the location
of the kinematic endpoint for background events, as given by eqs.~(\ref{mtab})--(\ref{mwb}) and~(\ref{mblmaxbknd}). As expected, the
blue (background) distributions always obey these kinematic endpoints.
(The figure has been constructed using parton-level events with no
detector modeling or combinatorial backgrounds --- these effects will
be added later on in section~\ref{sec:detector_resolution_etc}.)




On the other hand, the signal events, shown in red, significantly
violate the endpoints.  In fact, for all four variables considered in figure~\ref{fig:MT2cuts3}, the vast majority of signal events violate
the background endpoints. This is also to be expected, since study
point 1 was chosen specifically in region \textcolor{region2}{ii},
where all four endpoints are expected to be violated (see table~\ref{tab:regions}).   This also means that discovery (or
exclusion) for this study point should be relatively straightforward.  

We have noted above that, on an event by event basis,
\begin{equation}
M_{T2}(S; \tilde{m}) \le M_{2CC}(s; \tilde{m}). 
\label{eq:MT2-le-M2CC}
\end{equation}
This is confirmed in figure~\ref{fig:MT2cuts3}, where the (solid line)
$M_{2CC}$ distributions can be seen to be somewhat harder than the
respective $M_{T2}$ dashed line distributions.  As a result, cutting
on $M_{2CC}$ instead of $M_{T2}$ will result in a slightly higher signal
efficiency.  (The effect is most easily seen for subsystems $(a)$ and $(b)$.)

\begin{figure}[t]
\centering
\includegraphics[width=4.7cm]{ET_Topo2BP3}\hfill
\includegraphics[width=4.7cm]{PTl_Topo2BP3}\hfill
\includegraphics[width=4.7cm]{PTb_Topo2BP3}

]]></p>
<fig-group><caption><p><![CDATA[Signal (thick red lines) and background (thin blue lines)
  distributions of some standard variables for study point 1
  ($m_{\tilde t}=1110$\,GeV, $m_{\tilde \chi^\pm}=610$\,GeV and
  $m_{\tilde \nu}=110$\,GeV): the missing transverse energy $MET$ (left
  panel), the lepton transverse momentum $p_{T,\ell}$ (middle panel),
  and the $b$-quark transverse momentum $p_{T,b}$ (right
  panel).]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}



While the variables shown in figure~\ref{fig:MT2cuts3} already
allow study point 1 to be discovered or ruled out rather trivially,
the same can be accomplished with more conventional variables like
the missing transverse energy ($MET$), the lepton transverse momentum
$p_{T,\ell}$, or the $b$-quark transverse momentum $p_{T,b}$,
whose distributions for study point 1 are plotted in figure~\ref{fig:cuts3}.
In all three cases, even though the background distribution does not 
exhibit a strict endpoint, the signal and background distributions 
are very well separated, so the signal can be easily isolated.
Furthermore, as we can measure four endpoints in figure~\ref{fig:MT2cuts3}
and there are only three input mass parameters, full mass reconstruction 
in this case is also possible~\cite{Burns-ml-2008va}.
In conclusion, the analysis of study point 1 demonstrates that region
\textcolor{region2}{ii}
is ``easy" in the sense that the experimenter has a plethora of 
useful tools available for a discovery. It is therefore of interest 
to consider the other, more challenging, regions~of~figure~\ref{fig:map}.

\boldmath

]]></p>
</sec>
<sec><title>Study point 2: soft \texorpdfstring{$b$}{b}-jets in region vii</title>
<p><![CDATA[
\label{subsec:point-2}
\unboldmath


Our second example is in region \textcolor{region7}{vii}, in which
only the variables in subsystem $(b)$ have endpoint violations. The
mass spectrum for study point 2 is given in table~\ref{tab:study-points}.  We notice the relative degeneracy
between the stop and chargino masses, which causes the endpoints of
the signal $M_{T2}(a)$ and $m_{b\ell}$ distributions to be relatively
low.  (See eqs.~(\ref{susymta}) and~(\ref{susymbl}).)  In addition, the
sneutrino mass has been chosen so that the signal endpoint~(\ref{susymtab}) of the $M_{T2}(ab)$ variable is also below the
standard model expectation of~(\ref{mtab}).  This leaves the
$M_{T2}(b)$ variable as the only viable alternative in region
\textcolor{region7}{vii}.

\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{M2_Subsystem_ab_Topo2BP1}\quad 
\includegraphics[width=7.0cm]{M2_Subsystem_a_Topo2BP1}\\[3pt]
\includegraphics[width=7.0cm]{M2_Subsystem_b_Topo2BP1}\quad 
\includegraphics[width=7.0cm]{Mbl_Topo2BP1}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:MT2cuts3}, but with signal events for
study point 2 ($m_{\tilde t}=215$\,GeV,  $m_{\tilde \chi^\pm}=210$\,GeV and $m_{\tilde \nu}=110$\,GeV).]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}



These observations are illustrated in figure~\ref{fig:MT2cuts1} where we
plot the relevant invariant mass distributions for study point 2 using
the same conventions as in figure~\ref{fig:MT2cuts3} above.  Again, we
assume that all signal events have the event topology shown in figure~\ref{fig:DecayTopologies}(b), i.e., both stops decay according to
Topology 1.  Figure~\ref{fig:MT2cuts1} confirms that $M_{T2}(b)$ is a
good variable to cut on: placing a high pass cut with threshold just
above $m_t$ would eliminate all of the background, while leaving
almost half of the signal.  
To determine the optimal value of the threshold and the effectiveness
of the cut requires realistic detector simulation (see section~\ref{sec:detector_resolution_etc}), but it is clear that such
a cut will often be useful. This observation is not new --- the
variable $M_{T2}(b)$ has been discussed in the literature
under various names, e.g., $M_{T2}^{(210)}$~\cite{Burns-ml-2008va} and
dileptonic
$M_{T2}$~\cite{Plehn-ml-2012pr,Kilic-ml-2012kw,Chakraborty-ml-2013moa,Cho-ml-2008cu}.  Here
we would like to contrast $M_{T2}(b)$ to the alternative on-shell
constrained variable $M_{2CC}(b)$. The advantage of the latter is the
slightly higher signal efficiency.  On the other hand, the advantage
of the traditional $M_{T2}$ is its simplicity --- in its calculation,
one does not have to identify the $b$-jets, thus, one avoids
combinatorial ambiguities and the additional penalty due to
$b$-tagging.

Note that the signal and background distributions for the other three
variables in figure~\ref{fig:MT2cuts1}: $M_{T2}(ab)$, $M_{T2}(a)$, and
$m_{b\ell}$, also appear to be quite different, so one might wonder
whether they could be useful if the cut were inverted (i.e., if one
performs a low pass cut). However, we expect other background
processes besides $t\bar{t}$ to contribute events at low values and
swamp the signal~\cite{Barr-ml-2009wu,Allanach-ml-2011ej,Murayama-ml-2011hj}. 
Such backgrounds may not be as well-understood,
which is why in this study we shall only consider \emph{high pass} cuts 
on the invariant mass variables.\footnote{Additionally we remind 
the reader that figure~\ref{fig:MT2cuts1} (and analogous) figures 
throughout the work, depict unit normalized
distributions for signal and background; in reality the background
rates will be far higher than the signal rates in any realistic
model.  This also makes it more challenging to utilize the differences
in shape for a given variable \emph{below} the endpoint; endpoint
violation is the preferred feature for discovery.}



Having identified $M_{T2}(b)$ and $M_{2CC}(b)$ as promising variables,
one might wonder how the more conventional variables would perform in
this case.  In figure~\ref{fig:cuts1}, we show parton-level signal and
background distributions for study point 2 for the three more
traditional variables considered in figure~\ref{fig:cuts3}: $MET$,
$p_{T,\ell}$, and $p_{T,b}$. Since the stop-chargino mass splitting is
rather small, the $b$-jets are quite soft and would often not be
reconstructed.  On the other hand, the $MET$ and $p_{T,\ell}$
distributions show some separation between signal and background,
but the separation is less clear than we observed in the case of $M_{T2}(b)$
(the lower left panel of figure~\ref{fig:MT2cuts3}).  Therefore, placing
cuts on $MET$ and $p_{T,\ell}$ would not be as effective as cutting on
$M_{T2}(b)$.


\begin{figure}[t]
\centering
\includegraphics[width=4.7cm]{ET_Topo2BP1}\hfill
\includegraphics[width=4.7cm]{PTl_Topo2BP1}\hfill
\includegraphics[width=4.7cm]{PTb_Topo2BP1}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:cuts3}, but for signal events for
  study point 2 ($m_{\tilde t}=215$\,GeV, $m_{\tilde \chi^\pm}=210$\,GeV
  and $m_{\tilde \nu}=110$\,GeV).]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}



]]></p>
</sec>
<sec><title>Study point 3: soft leptons in region iii</title>
<p><![CDATA[
\label{subsec:point-3}

Our next example illustrates the complementarity of the subsystem
invariant mass variables.  In the previous subsection~(\ref{subsec:point-2}), we considered a signal study point with soft
$b$-jets and relatively hard leptons, as seen in figure~\ref{fig:cuts1}. Now we shall discuss the opposite situation,
when the leptons are relatively soft, while the jets are hard.  For
this purpose, we focus on study point 3 in region
\textcolor{region3}{iii}, where according to table~\ref{tab:regions}
we expect endpoint violations for $M_{T2}(ab)$, $M_{T2}(a)$, and
$m_{b\ell}$.
\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{M2_Subsystem_ab_Topo2BP2}\quad
\includegraphics[width=7.0cm]{M2_Subsystem_a_Topo2BP2}\\[3pt]
\includegraphics[width=7.0cm]{M2_Subsystem_b_Topo2BP2}\quad
\includegraphics[width=7.0cm]{Mbl_Topo2BP2}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:MT2cuts3}, but with signal events for
study point 3 ($m_{\tilde t}=630$\,GeV,  $m_{\tilde \chi^\pm}=130$\,GeV and $m_{\tilde \nu}=110$\,GeV). ]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
\begin{figure}[t]
\centering
\includegraphics[width=4.7cm]{ET_Topo2BP2}\hfill
\includegraphics[width=4.7cm]{PTl_Topo2BP2}\hfill
\includegraphics[width=4.7cm]{PTb_Topo2BP2}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:cuts3}, but for signal events for
study point 3 ($m_{\tilde t}=630$\,GeV,  $m_{\tilde \chi^\pm}=130$\,GeV and $m_{\tilde \nu}=110$\,GeV).]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
This feature is demonstrated in figure~\ref{fig:MT2cuts2}, where we
again compare the signal and background distributions for the same
four types of variables as in figures~\ref{fig:MT2cuts3} and~\ref{fig:MT2cuts1}.  We see that this time, as expected, the dilepton
$M_{T2}$ variable ($M_{T2}(b)$ in our notation) is suboptimal; due
to the softness of the leptons, the signal $M_{T2}(b)$ and
$M_{2CC}(b)$ distributions lie entirely within the background region.
On the other hand, the other three variables perform very well, unlike
the case in section~\ref{subsec:point-2}. In particular, the subsystem
$(ab)$ variables alone could possibly remove the background with
virtually no loss of signal.  The subsystem $(a)$ variables are also
promising; the use of $M_{2CC}(a)$ seems slightly more effective than
the use of $M_{T2}(a)$. Finally, the usual invariant mass, $m_{b\ell}$,
allows one to separate signal and background, but the signal loss is
more significant for this variable.


The lesson from figures~\ref{fig:MT2cuts1} and~\ref{fig:MT2cuts2} is
that in order to efficiently probe the full mass parameter space of figure~\ref{fig:map}, one would have to design an analysis which
utilizes the full complement of subsystem invariant mass variables,
since different variables are optimal in different regions. 
Of course, one should not overlook the more conventional variables.
In figure~\ref{fig:cuts2} we show the signal and background
distributions of $MET$, $p_{T,\ell}$, and $p_{T,b}$ for study point 3.
As expected, the lepton $p_T$ distribution for the signal is rather
soft, but the signal distributions for both $MET$ and $p_{T,b}$ have
long tails which extend to the right of the bulk of the corresponding
background distribution.  This suggests that $MET$ and $p_{T,b}$ could
also play a useful role in the analysis.


]]></p>
</sec>
<sec><title>Study point 4: a difficult case in region vi</title>
<p><![CDATA[
\label{subsec:point-4}

Our final example for the signal event topology of figure~\ref{fig:DecayTopologies}(b) is a study point in the most
challenging region, \textcolor{region6}{vi}, where no endpoint
violations should occur. The mass spectrum for study point 4 is given
in table~\ref{tab:study-points}, and the resulting signal and
background distributions are displayed in figures~\ref{fig:MT2cuts4} and~\ref{fig:cuts4}.
\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{M2_Subsystem_ab_Topo2BP4}\quad 
\includegraphics[width=7.0cm]{M2_Subsystem_a_Topo2BP4}\\[3pt]
\includegraphics[width=7.0cm]{M2_Subsystem_b_Topo2BP4}\quad
\includegraphics[width=7.0cm]{Mbl_Topo2BP4}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:MT2cuts3}, but with signal events for
study point 4 ($m_{\tilde t}=174$\,GeV,  $m_{\tilde \chi^\pm}=150$\,GeV and $m_{\tilde \nu}=110$\,GeV). ]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
\begin{figure}[t]
\centering
\includegraphics[width=4.7cm]{ET_Topo2BP4}\hfill
\includegraphics[width=4.7cm]{PTl_Topo2BP4}\hfill
\includegraphics[width=4.7cm]{PTb_Topo2BP4}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:cuts3}, but for signal events for
study point 4 ($m_{\tilde t}=174$\,GeV,  $m_{\tilde \chi^\pm}=150$\,GeV and $m_{\tilde \nu}=110$\,GeV).]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
Discovery is clearly challenging in this scenario, as the
$b$-jets will be quite soft, while the lepton $p_T$ and $MET$
distributions for the signal are very similar to those for the
background.  The invariant mass variables in figure~\ref{fig:MT2cuts4}
are not particularly helpful either, since the kinematic endpoints of the
signal distributions are always below those of the background. Whether
stops can be discovered at study point 4 thus remains an open
question, which we shall revisit in section~\ref{sec:mixed-events}.





]]></p>
</sec>
</sec>
<sec><title>\texorpdfstring{\boldmath $M_2$}{M(2)} endpoint study for Topology 2</title>
<p><![CDATA[
\label{sec:M2_endpoint_study_for_topology_2}

In this section we shall focus on the other symmetric signal event
topology in figure~\ref{fig:DecayTopologies}(c), when both stops decay
according to Topology 2 in figure~\ref{fig:process}(b).  The on-shell
constrained invariant mass variables discussed in the previous section
will be useful here as well, since they were constructed with the {\em
  background} topology in mind, which has not changed.  Even though
the signal event topology is now more complicated (there are two
invisible particles in each decay chain), the signal distributions
still exhibit kinematic endpoints.   We find that the $M_{T2}$
endpoints are given by (see, e.g.,~\cite{Mahbubani-ml-2012kx})
\begin{align}
M_{T2}^{\rm max}(ab;\tilde m=0) & =   2\, C_{+}(m_t, m_\nu), 
\label{susymtab2}\\ 
M_{T2}^{\rm max}(a;\tilde m=m_W) & = 
C_{+}(m_{t},m_{W})+\sqrt{C^2_{+}(m_{t},m_{W})+m^2_W},
\label{susymta2} \\ 
M_{T2}^{\rm max}(b;\tilde m=0) & =  2\, \sqrt{C_{+}(m_{W},m_{\nu})\, 
\left[ C_{+}(m_{W},m_{\nu})- C_{-}(m_{t},m_{W}) \right]},
\label{susymwb2}
\end{align}
where
\begin{align}
C_{\pm}(x,y) &\equiv  \frac{x^2-y^2}{4\,m_{\tilde t}\,x^2}\, 
\left\{m^2_{\tilde t}+ m^2_{t} -m^2_{\tilde\chi^0}
\pm \sqrt{\lambda\left(m^2_{\tilde t}, m^2_{t}, m^2_{\tilde\chi^0}\right)}\right\},
\\
\lambda(x,y,z)&\equiv x^2+y^2+z^2 - 2xy - 2yz - 2zx.
\end{align}
Upon careful examination of eqs.~(\ref{susymtab2})--(\ref{susymwb2}), one
can show that these kinematic endpoints are always above the
corresponding background endpoints~(\ref{mtab})--(\ref{mwb}), as long as
the channel $\tilde t \to t \tilde\chi^0$ is open (i.e., the decay is
kinematically allowed).  As we move close
to the threshold for $\tilde t \to t \tilde\chi^0$, the signal
kinematic endpoints~(\ref{susymtab2})--(\ref{susymwb2}) converge to the
corresponding SM values~(\ref{mtab})--(\ref{mwb}), and discovery becomes
very challenging.  In this section, therefore, we shall consider two
study points: one above this threshold and one at threshold. The mass
spectra for those study points  are listed in table~\ref{tab:study-points-topo2}.

\begin{table}[t]
\centering
\renewcommand{\arraystretch}{1.1}
\begin{tabular}{| c | c | c | c | c |}
\hline
Study Point & Stop Mass & Neutralino Mass \\
\hline 
5  & $300$\,GeV & $100$\,GeV  \\
6  & $174$\,GeV & $0$\,GeV  \\
\hline
\end{tabular}

]]></p>
<table-wrap><caption><p><![CDATA[The stop and neutralino masses
  for the two study points considered in this section. ]]></p></caption></table-wrap>
<p><![CDATA[
\end{table}

\boldmath

]]></p>
<sec><title>Study point 5: a case above the \texorpdfstring{$\tilde t \to t
  \tilde\chi^0$}{tilde(t) --&gt; t tilde(chi)**0} threshold</title>
<p><![CDATA[
\label{subsec:point-5}
\unboldmath


We first discuss study point 5, where the mass splitting is large
enough that the decay $\tilde t \to t \tilde\chi^0$ is open, and the
resulting top quark is on-shell.  The corresponding kinematic
distributions are shown in figures~\ref{fig:MT2cuts5} and~\ref{fig:cuts5}.
\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{M2_Subsystem_ab_Topo1BP1}\quad 
\includegraphics[width=7.0cm]{M2_Subsystem_a_Topo1BP1}\\[3pt]
\includegraphics[width=7.0cm]{M2_Subsystem_b_Topo1BP1}\quad 
\includegraphics[width=7.0cm]{Mbl_Topo1BP1}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:MT2cuts3}, but for signal events with the event topology of figure~\ref{fig:DecayTopologies}(c)
for study point 5 ($m_{\tilde t}=300$\,GeV and $m_{\tilde \chi^0}=100$\,GeV).]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
\begin{figure}[t]
\centering
\includegraphics[width=4.7cm]{ET_Topo1BP1}\hfill
\includegraphics[width=4.7cm]{PTl_Topo1BP1}\hfill
\includegraphics[width=4.7cm]{PTb_Topo1BP1}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:cuts3}, but for signal events with the event topology of figure~\ref{fig:DecayTopologies}(c)
for study point 5 ($m_{\tilde t}=300$\,GeV and $m_{\tilde \chi^0}=100$\,GeV)\@. 
]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
It is clear that this is already a challenging case --- the signal and
background distributions for the ``conventional" variables in figure~\ref{fig:cuts5} are rather similar. The jet and lepton $p_T$
spectra are governed by the known mass differences between the \emph{SM}
  particles $t$, $W^\pm$, and $\nu$, thus, there is very little
distinction between the signal and background $p_T$
distributions. Similarly, the $m_{b\ell}$ distribution in figure~\ref{fig:MT2cuts5} is the same for signal and
background. The $MET$ distribution in figure~\ref{fig:cuts5} is slightly
harder for the signal, due to the presence of two
additional invisible particles. However, the effect is very small and
hence unlikely to be useful in practice.




This motivates the use of the $M_{T2}$ and $M_{2\sqcup\sqcup}$
variables whose distributions are shown in the first three panels of figure~\ref{fig:MT2cuts5}. As anticipated from eqs.~(\ref{susymtab2})--(\ref{susymwb2}), for all three subsystems
$(ab)$, $(a)$, and $(b)$, the signal distributions for the $M_{T2}$
variable have a tail which extends beyond the background
endpoints. This effect is most pronounced for subsystem $(ab)$ and
less so for subsystem $(a)$.

Note how the situation improves if one were to use the on-shell
constrained variable $M_{2CC}$ (solid lines) instead of $M_{T2}$
(dashed lines). For background events, $M_{2CC}$ is computed by
applying the correct kinematic constraints; therefore, the
kinematic endpoints~(\ref{mtab})--(\ref{mwb}) are still obeyed. For
signal events, we get a somewhat different story --- a much larger
fraction of signal events now violate these endpoints, leading
to an improvement in the signal efficiency.  The largest benefit is
observed in the case of subsystem $(a)$, for which previously the
$M_{T2}$ variable was the least helpful.  There are two separate reasons
why $M_{2CC}$ separates signal from background better than $M_{T2}$:
\begin{enumerate}
\item Due to the hierarchy~(\ref{eq:inequality}), the $M_{2CC}$
  distributions are harder than the $M_{T2}$ distributions, thus more
  signal events are expected to migrate above the background
  endpoint.  The shape difference between the $M_{2CC}$ and $M_{T2}$
  distributions is especially noticeable in the case of subsystems
  $(a)$ and $(b)$ in figure~\ref{fig:MT2cuts5}. Notice, in particular,
  the completely different shapes of the $M_{T2}(a)$ and $M_{2CC}(a)$
  distributions, as well as the disappearance of the big spike at
  $M_{T2}(b)=0$. 
\item For signal events, the $M_{2CC}$ kinematic endpoints
  themselves are even higher\footnote{We have not attempted to obtain
    analytical formulas analogous to~(\ref{susymtab2})--(\ref{susymwb2})
    for the $M_{2CC}$ kinematic endpoints, but our numerical studies
    clearly showed that the bounds~(\ref{susymtab2})--(\ref{susymwb2})
    themselves are violated in the case of the $M_{2CC}$ variable. }
  than the $M_{T2}$ kinematic endpoints given in eqs.~(\ref{susymtab2})--(\ref{susymwb2}).
This can be readily observed in figure~\ref{fig:MT2cuts5}, where the
signal $M_{2CC}$ distributions (red solid lines) extend to higher
values than the signal $M_{T2}$ distributions (red dashed lines).
Contrast this situation with the examples considered in section~\ref{sec:M2_endpoint_study_for_topology_1}, when $M_{T2}$ and
$M_{2CC}$ always shared the same kinematic endpoint. There, the
signal event topology (figure~\ref{fig:DecayTopologies}(b)) was the same
as the background event topology
(figure~\ref{fig:DecayTopologies}(a)). As a result, the kinematic
constraints being imposed in the calculation of $M_{2CC}$ did
correspond to the actual physics of the signal events.  Now, in the
case of study point 5, the signal event topology of figure~\ref{fig:DecayTopologies}(c) is completely different --- in a
sense, one is applying ``the wrong" constraints when calculating
$M_{2CC}$. Somewhat paradoxically then, figure~\ref{fig:MT2cuts5}
teaches us that one obtains a beneficial result, despite
applying ``the wrong" constraints. 
\end{enumerate}



\boldmath

]]></p>
</sec>
<sec><title>Study point 6: a case at the \texorpdfstring{$\tilde t \to t \tilde\chi^0$}{tilde(t) --&gt; t tilde(chi)**0} threshold</title>
<p><![CDATA[
\label{subsec:point-6}
\unboldmath


Our last example is a very difficult one: study point 6 in table~\ref{tab:study-points-topo2}.  Here the new physics mass spectrum is
such that the decay $\tilde t \to t \tilde\chi^0$ occurs exactly at
threshold. As a result, the (massless) neutralinos carry away a
negligible amount of momentum, and the signal events look very
top-like. This is illustrated in figures~\ref{fig:MT2cuts6} and~\ref{fig:cuts6}, where we compare our standard set of kinematic
distributions for signal and background. 
\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{M2_Subsystem_ab_Topo1BP2}\quad 
\includegraphics[width=7.0cm]{M2_Subsystem_a_Topo1BP2}\\[3pt]
\includegraphics[width=7.0cm]{M2_Subsystem_b_Topo1BP2}\quad 
\includegraphics[width=7.0cm]{Mbl_Topo1BP2}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:MT2cuts3}, but for signal events with the event topology of figure~\ref{fig:DecayTopologies}(c)
for study point 6 ($m_{\tilde t}=174$\,GeV and $m_{\tilde \chi^0}=0$\,GeV)\@.
]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
\begin{figure}[t]
\centering
\includegraphics[width=4.7cm]{ET_Topo1BP2}\hfill
\includegraphics[width=4.7cm]{PTl_Topo1BP2}\hfill
\includegraphics[width=4.7cm]{PTb_Topo1BP2}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:cuts3}, but for signal events with
  the event topology of figure~\ref{fig:DecayTopologies}(c) for study
  point 6 ($m_{\tilde t}=174$\,GeV and $m_{\tilde \chi^0}=0$\,GeV)\@. 
]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}

Figure~\ref{fig:cuts6} shows that the $p_T$ distributions and the $MET$
distribution are almost identical for signal and background. The
invariant mass distributions from figure~\ref{fig:MT2cuts6} are also
very similar; there are slight differences in the shapes due to the
top quarks in the signal being more likely to be off-shell, but  the
kinematic endpoints are the same.  Thus, barring a shape-based
analysis, there are no obvious cuts which could discriminate signal
from background. Therefore, just like study point 4, this would be a
very difficult, and most likely impossible, scenario for discovery
using these methods. As before, we shall leave this as an open
question to be revisited in section~\ref{sec:mixed-events}.


]]></p>
</sec>
</sec>
<sec><title>\texorpdfstring{\boldmath $M_2$}{M(2)} endpoint study for mixed events </title>
<p><![CDATA[
\label{sec:mixed-events}

In this section, we shall consider signal events with the mixed event topology of figure~\ref{fig:DecayTopologies}(d).  In doing so, we are motivated by two factors:
\begin{itemize}
\item In any realistic model, the stop is likely to have several
  relevant decay modes.  (Here we consider the simplest scenario with
  only the two decay modes from figure~\ref{fig:process}.)  Since the
  stops are pair-produced, the number of signal events in each
  symmetric channel is proportional to the corresponding branching
  ratio squared.  For mixed events, where the two stops decay {\em
    differently}, the number of signal events benefits from an
  additional combinatorial factor of 2.
\item In the course of our study of the symmetric event topologies from figure~\ref{fig:DecayTopologies}(b) (in section~\ref{sec:M2_endpoint_study_for_topology_1}) and figure~\ref{fig:DecayTopologies}(c) (in section~\ref{sec:M2_endpoint_study_for_topology_2}), we determined
that there are ``blind spots" in the mass parameter space, where the
signal resembles the background, and discovery is very
challenging. Study points 4 and 6 are examples of such difficult
cases. In this section, therefore, we shall investigate the question
of whether one can recover some sensitivity by considering mixed
events constructed from precisely those two difficult cases. In other
words, we consider events with the event topology of figure~\ref{fig:DecayTopologies}(d), where the upper (lower) decay chain
corresponds to study point 4 (study point 6). (Study point 4 gives the
stop mass ($174$\,GeV) and the neutralino mass ($0$\,GeV); study point 6
uses the same stop mass, a chargino mass of $150$\,GeV and a sneutrino
mass of $110$\,GeV.)
We shall assume that the two stop decays occur with equal
branching fractions.
\end{itemize}

The idea to use mixed stop events was previously discussed in ref.~\cite{Graesser-ml-2012qy},
which suggested a new variable, ``topness", that quantifies how well an
event can be reconstructed under the top background hypothesis. In
order to calculate the ``topness" of an event, one minimizes the total
$\sqrt{s}$ of the event, making the reasonable ansatz that the
momentum configuration thus obtained provides a good approximation to the
true kinematics of the event~\cite{Graesser-ml-2012qy,Konar-ml-2008ei}. Our
approach is similar to the extent that the
 on-shell constrained invariant mass variables, like $M_{2CC}$, are also
 found by minimization, though not of the total $\sqrt{s}$ but of the parent
 mass in the respective subsystem.
 By imposing the symmetry constraints~(\ref{eq:parents-equal}) and~(\ref{eq:relatives-equal}), we focus on the one key difference
 between the signal and background events: 
 the signal event topology is asymmetric while the background event
 topology is symmetric.
 We can therefore expect that the constraints~(\ref{eq:parents-equal})
 and~(\ref{eq:relatives-equal})
 will affect signal and background events~\emph{differently}.

\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{MT2_Subsystem_ab_No_Combinatorics}\quad 
\includegraphics[width=7.0cm]{MT2_Subsystem_a_No_Combinatorics}\\[3pt]
\includegraphics[width=7.0cm]{MT2_Subsystem_b_No_Combinatorics}\quad
\includegraphics[width=7.0cm]{Mbl_No_Combinatorics}

]]></p>
<fig-group><caption><p><![CDATA[The same as figures~\ref{fig:MT2cuts4} and~\ref{fig:MT2cuts6},
  comparing distributions of the three $M_{T2}$ subsystem variables
  and the invariant mass $m_{b\ell}$ for the four types of events in figure~\ref{fig:DecayTopologies}. The mass spectrum corresponds to
  study points 4 and 6: $m_{\tilde t}=174$\,GeV, $m_{\tilde
    \chi^\pm}=150$\,GeV, $m_{\tilde \nu}=110$\,GeV, and $m_{\tilde
    \chi^0}=0$\,GeV\@. ]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}


\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{M2CC_Subsystem_ab_No_Combinatorics}\quad 
\includegraphics[width=7.0cm]{M2CC_Subsystem_a_No_Combinatorics}\\
\includegraphics[width=7.0cm]{M2CC_Subsystem_b_No_Combinatorics}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:MT2cuts5again}, but for the
  on-shell constrained variable $M_{2CC}$ in the three possible
  subsystems $(ab)$, $(a)$, and $(b)$.]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}

As a point of reference, we begin by showing distributions of
variables for which no improvement can be expected in the case of
mixed events in figure~\ref{fig:MT2cuts5again}. The four variables
depicted in the figure are the three subsystem $M_{T2}$ variables and
the invariant mass $m_{b\ell}$.  All four variables are calculated
here using the assignment of the leptons and b-quarks to the correct
decay chains.
Each panel contains four distributions, one for each event type from figure~\ref{fig:DecayTopologies}: the $t\bar{t}$ background from figure~\ref{fig:DecayTopologies}(a) (blue dot-dashed lines); the
symmetric $\tilde t\to b\tilde\chi^+_1$ signal events from figure~\ref{fig:DecayTopologies}(b) (magenta dotted lines); the
symmetric $\tilde t\to t\tilde\chi^0_1$ signal events from figure~\ref{fig:DecayTopologies}(c) 
(red dashed lines); and the asymmetric signal events from figure~\ref{fig:DecayTopologies}(d) (black solid lines).
The vertical black dashed line in each panel marks the location of the 
upper kinematic endpoint of the $t\bar{t}$ background distribution.
We see that for all three types of signal events (symmetric or asymmetric), 
the respective distributions do not violate the background kinematic
endpoints, thus discovery appears to be just as difficult with mixed
events as it was with the symmetric events considered earlier in sections~\ref{subsec:point-4} and~\ref{subsec:point-6}. 
This conclusion is easy to understand; the $M_{T2}$ variable is a
variable defined on the transverse plane, where it is impossible to
impose a 3+1-dimensional mass constraint like eq.~(\ref{eq:parents-equal})~or~(\ref{eq:relatives-equal}).

\begin{figure}[t]
\centering
\includegraphics[width=.95\textwidth]{figure19.png}

]]></p>
<fig-group><caption><p><![CDATA[The correlation between the two best performing variables
  from figure~\ref{fig:M2CCcuts5}, for the three types of signal events:
  pure Topology 1 from figure~\ref{fig:DecayTopologies}(b) (upper left),
  pure Topology 2 from figure~\ref{fig:DecayTopologies}(c) (upper
  right), mixed topology from figure~\ref{fig:DecayTopologies}(d) (lower
  left), and combined total (lower~right).]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}

The situation is quite different when we consider distributions of
on-shell constrained variables, $M_{2CC}$, for which the constraints 
of eqs.~(\ref{eq:parents-equal}) or~(\ref{eq:relatives-equal}) are imposed. 
As seen in figure~\ref{fig:M2CCcuts5},
the signal distributions for mixed events may now exhibit endpoint violation, 
even when the signal distributions for symmetric events do not.
The effect is most pronounced in the case of the subsystem variable $M_{2CC}(b)$;
for the subsystem variable $M_{2CC}(ab)$ it is less noticeable, while for
$M_{2CC}(a)$ it is absent altogether. Figure~\ref{fig:M2CCcuts5} showcases the main result 
of this section: that with the help of an appropriately chosen 
on-shell constrained variable (in this case $M_{2CC}(b)$), one can obtain 
a relatively good separation of signal from background \emph{for mixed events}. 
It is worth emphasizing that this separation was achieved for a very unfavorable 
choice of mass parameters, as the study points 4 and 6 were not
observable using events where both decay chains had the same topology. 

\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{M2CC_Subsystem_ab_Combinatorics}\quad 
\includegraphics[width=7.0cm]{M2CC_Subsystem_a_Combinatorics}
\includegraphics[width=7.0cm]{M2CC_Subsystem_b_Combinatorics}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:M2CCcuts5}, but including the effects of
combinatorics: for each event, we try both possible assignments of the
lepton-$b$pairs, and plot the smaller of the two resulting $M_{2CC}$
values. ]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}

Given that endpoint violation was observed for both $M_{2CC}(b)$
and $M_{2CC}(ab)$, it is worth investigating the possible correlation
between those two variables.  In figure~\ref{fig:nocomb} we show
two-dimensional plots exhibiting those correlations.  We consider
separately the three types of signal events: pure Topology 1 from figure~\ref{fig:DecayTopologies}(b) (upper left), pure Topology 2 from figure~\ref{fig:DecayTopologies}(c) (upper right), and the mixed
topology from figure~\ref{fig:DecayTopologies}(d) (lower left).
Finally, the lower right panel in figure~\ref{fig:nocomb} shows the
result for the full signal sample, with equal branching fractions for
Topology 1 and Topology 2.   The black dashed lines in figure~\ref{fig:nocomb} mark the locations of the expected upper
kinematic endpoints for background events, following eqs.~(\ref{mwb})
and~(\ref{mtab}).  Any events which appear to the right of the vertical
black dashed lines and/or above the horizontal black dashed lines in figure~\ref{fig:nocomb} are expected to be signal-like. In agreement
with figure~\ref{fig:M2CCcuts5}, we see that for symmetric signal event
topologies (the upper two panels in figure~\ref{fig:nocomb}), the signal
events are contained within the ``background-like'' rectangular region
adjacent to the origin.  On the other hand, for the asymmetric event
topology of figure~\ref{fig:DecayTopologies}(d) (the lower left panel),
many signal events leak out of the background-like box.  The figure
also reveals a linear correlation between $M_{2CC}(b)$ and
$M_{2CC}(ab)$.  Furthermore, the slope is such that if an event
violates the background $M_{2CC}(ab)$ endpoint~(\ref{mtab}), it also
necessarily violates the background $M_{2CC}(b)$ endpoint~(\ref{mwb}),
while the reverse is not true. We therefore conclude that the
$M_{2CC}(b)$ distribution alone is sufficient in separating signal
from background in this scenario with mixed events.

\begin{figure}[t]
\centering
\includegraphics{figure21.png}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:nocomb}, but including the effects of
combinatorics as in figure~\ref{fig:M2CCcuts5comb}.]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}

In our discussion so far in this section, we have been ignoring the
combinatorial problem arising when we try to pair up the two $b$-jets
with the two leptons. Since the $b$-quark charge  is not measured, we
have two possible pairings, each resulting in a candidate value for
the kinematic variable. Since we are interested in upper kinematic
endpoints, the simplest solution is to consider both pairings and then
pick the one which gives the \emph{smaller} value for the kinematic
variable. This approach has been followed  in recreating figures~\ref{fig:M2CCcuts5} and~\ref{fig:nocomb} as figures~\ref{fig:M2CCcuts5comb} and~\ref{fig:comb}, respectively.  As
expected, this procedure tends to shift all distributions towards
lower values, thus the number of signal events which violate the
background endpoints is fewer than before; compare, e.g., the
$M_{2CC}(b)$ distributions for mixed events in figures~\ref{fig:M2CCcuts5} and~\ref{fig:M2CCcuts5comb}.  Nevertheless,
the effect is still present, offering hope that difficult cases like
study points 4 and 6 could perhaps best be looked for in such mixed
event topologies instead.


]]></p>
</sec>
<sec><title>Results with realistic detector simulation</title>
<p><![CDATA[
\label{sec:detector_resolution_etc}

In the previous three sections we saw that the $M_{T2}$ and $M_{2CC}$ variables
allow us to identify signal events as tails which extend beyond the
upper kinematic endpoint for background events. However, in a
realistic  experiment, the background distributions themselves may
acquire high tails, for a variety of reasons.  This is why it is
necessary to test our previous observations, which were made at
parton level, with realistic simulation, including the effects of
detector resolution, initial and final state radiation, jet
reconstruction, cuts, etc.  It is clear that our positive conclusions
drawn for fortuitous cases of new physics like study point 1 will
survive all these complications, therefore in this section we shall
only focus on the difficult scenario discussed in section~\ref{sec:mixed-events}, i.e., the mixed events which were a
hybrid between the difficult study points 4 from section~\ref{subsec:point-4} and 6 from section~\ref{subsec:point-6}.


]]></p>
<sec><title>Event simulation details</title>
<p><![CDATA[
\label{sec:simulation}

As before, the parton-level event generation is done 
by \texttt{MadGraph\_aMC@NLO}~\cite{Alwall-ml-2014hca}, where by default
the parton distributions are evaluated by \texttt{NNPDF23}~\cite{Ball-ml-2012cx}.
The relevant output is then piped through \texttt{Pythia 6.4}~\cite{Sjostrand-ml-2006za} 
and \texttt{Delphes3}~\cite{deFavereau-ml-2013fsa}. For both signal and 
$t\bar{t}$ background, the decays of top quarks are handled by \texttt{Pythia 6.4}, 
while Topology 1 is explicitly generated by \texttt{MadGraph\_aMC@NLO} 
without any prior cuts. All simulations are performed at leading order
for a $pp$ collider of $\sqrt{s}=14$\,TeV\@. 

For the signal process, we assume that the branching ratio of
Topology 1 relative to Topology 2 is $1:1$. 
In addition, in Topology 1, the chargino is forced to decay
exclusively into a sneutrino (which may further decay invisibly), and
a lepton (i.e., electron and muon only). In Topology 2, the stop
decays to the lightest neutralino and a top quark, which subsequently
decays with the relevant branching ratios predicted in the SM\@. For our
purposes, we only consider the dilepton final state, in which
both top quarks decay leptonically.  The input top quark mass is set
to $173$\,GeV, while the $W^\pm$ gauge boson mass is $80$\,GeV\@. 
Jets are reconstructed with the anti-$k_t$ algorithm~\cite{Cacciari-ml-2008gp}, 
using a radius parameter $R=0.5$. The $b$-tagging efficiency is taken
to be $70$\%, while light quark jets~are~mis-tagged~at~the~rate~of~$1$\%. 



Given the final state $2b+2\ell+\met$, in principle there are several
sources of SM background that need to be taken into account. In order
to suppress the reducible SM backgrounds, we apply the following
pre-selection cuts:
\begin{itemize}
\item The event must contain exactly two opposite sign leptons with 
$p_T>10$\,GeV and $|\eta|<2.5$ ($|\eta|<2.4$) for electrons (muons).
\item In order to reduce background from low mass resonances, 
in the $ee$ and $\mu\mu$ channels, we demand $m_{ee/\mu\mu}>20$\,GeV\@. Furthermore, to reduce the $Z+$jets
background, events with dilepton masses 
within the $Z$-mass window are vetoed by requiring
$|m_{ee/\mu\mu}-m_Z| > 15$\,GeV\@.
\item To further suppress Drell-Yan, for the $ee$ and $\mu\mu$
  channels, we apply a missing transverse energy cut of $\met>40$\,GeV\@.
 \item The event is required to have $\geq 2$ jets with $p_T > 30$\,GeV and 
$|\eta|<2.4$. It is also required that exactly two of these jets are $b$-tagged. 
\end{itemize}
After these cuts, we are left with $t\bar{t}$ as the dominant (irreducible) 
background, and this will be the only background process we will consider here.
The posterior cuts using $M_2$ variables are imposed for
events already passing the above set of pre-selection cuts. 

\boldmath

]]></p>
</sec>
<sec><title>Results for \texorpdfstring{$M_{T2}$}{M(T2)} and \texorpdfstring{$M_{2CC}$}{M(2CC)}</title>
<p><![CDATA[
\label{sec:MT2}
\unboldmath


We first revisit our results from section~\ref{sec:mixed-events}, this
time including the effects of detector simulation and
combinatorics. As before, we consider both pairings of the two tagged
$b$-jets and the two leptons, and use the \emph{smaller} of the two
resulting values for the kinematic variable under consideration.
However, unlike the plots in section~\ref{sec:mixed-events}, here we
do not separate the three types of signal events (Topology 1,
Topology 2 and mixed), since in the real experiment there is no way to
tell which is which.

\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{MT2_SubSys_bl_Delphes_Log}\quad 
\includegraphics[width=7.0cm]{MT2_SubSys_b_Delphes_Log}
\includegraphics[width=7.0cm]{MT2_SubSys_l_Delphes_Log}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:MT2cuts5again}, 
but including the effects of combinatorics and detector simulation.
The vertical black dashed lines denote the expected $M_{T2}$ endpoints
of the $t\bar{t}$ background in each subsystem. ]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}

Figure~\ref{fig:MT2detector} compares the signal and $t\bar{t}$
background distributions for the three different $M_{T2}$ subsystem variables.
As expected from the parton-level results in section~\ref{sec:mixed-events} (see figure~\ref{fig:MT2cuts5again}),
the discrimination power in the high tail region is relatively poor, since
the signal and the background events obey the same kinematic endpoints.

\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{M2CC_SubSys_bl_Delphes_Log}\quad 
\includegraphics[width=7.0cm]{M2CC_SubSys_b_Delphes_Log}
\includegraphics[width=7.0cm]{M2CC_SubSys_l_Delphes_Log}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:MT2detector}, but for the corresponding
$M_{2CC}$ variables.]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}

Figure~\ref{fig:M2CCdetector} shows the corresponding $M_{2CC}$
distributions for signal and $t\bar{t}$ background events in the three
subsystems.  As anticipated from the parton-level result in figures~\ref{fig:M2CCcuts5} and~\ref{fig:M2CCcuts5comb}, there is a
noticeable improvement in the $(b)$ subsystem (for which the visible
particle is a lepton) as seen in the lower panel
and a slight improvement in the $(ab)$ subsystem as well. 
Therefore, one would expect that a minimum $M_{2CC}$ cut would be beneficial.
The optimal value of the cut would depend on the expected signal
cross-section, and on the assumed systematic uncertainty on the
background normalization in the high tail region.


\pagebreak

A careful comparison of the parton-level results in figures~\ref{fig:MT2cuts5again},~\ref{fig:M2CCcuts5}
and~\ref{fig:M2CCcuts5comb} versus the detector-level results in figures~\ref{fig:MT2detector} and~\ref{fig:M2CCdetector} reveals that at
the detector level the background distributions develop high tails
which, unless properly understood, could be confused with a signal.
We have checked that in the majority of cases, background events
populate the high tail due to imperfect $b$-tagging. A typical event
looks as follows: one of the two $b$-jets is either too soft to pass
the jet ID cuts, or is not tagged as a  $b$-jet.  (Recall that the
$b$-jet tagging inefficiency is $30$\%.)  Instead, a gluon from initial
state radiation (ISR) forms a hard jet which is subsequently mistagged
as a $b$-jet. Thus in computing the $M_{T2}$ and $M_{2CC}$ variables
one is using the wrong $b$-jet object, which leads to the endpoint
violation.
An improvement in the $b$-tagging algorithm, especially one which
lowers the mistag rate for ordinary QCD jets, would help alleviate
this problem.

\boldmath

]]></p>
</sec>
<sec><title>Results for the relative shift from \texorpdfstring{$M_{T2}$}{M(T2)} to \texorpdfstring{$M_{2CC}$}{M(2CC)}</title>
<p><![CDATA[
\label{sec:shift}
\unboldmath


In all of our analysis so far, we have relied on the existence of the background 
kinematic endpoints~(\ref{mtab})--(\ref{mwb}) and focused on the high tails
\emph{above} those endpoints. Ideally, this kinematic region should be
populated only by signal events, even when one accounts for the
two-fold combinatorial ambiguity in pairing the leptons and the
$b$-jets. Unfortunately, as we have already seen, this straightforward
approach has two drawbacks:
\begin{itemize}
\item \emph{Presence of high tails in the background distributions.}
While in theory the background distributions are not supposed to
extend beyond their kinematic endpoints, in practice this is not
always the case. Such tails were readily observed in figures~\ref{fig:MT2detector} and~\ref{fig:M2CCdetector}, which were
obtained using realistic detector simulation. 
\item \emph{Low signal efficiency.} Unless we are dealing with a new
  physics model with a widely split spectrum (see related discussion
  in section~\ref{subsec:regions-and-study-points}), a significant
  fraction of the signal events will also lie below the background
  kinematic endpoints, thus by cutting at or near the endpoint, we
  will be removing a large chunk of signal events as well.
This was very evident in the ``worst case" scenarios like study points
4 and 6, or the mixed event case discussed in sections~\ref{sec:mixed-events} and~\ref{sec:MT2}.
\end{itemize}
These two problems suggest that we should reexamine the region \emph{below} 
the background kinematic endpoints and search for a good
discriminating variable which would be applicable to that region as
well. As in section~\ref{sec:mixed-events}, our goal will be to target
signal events with the mixed event topology of figure~\ref{fig:DecayTopologies}(d).

To begin with, recall the main difference between the background
events described by figure~\ref{fig:DecayTopologies}(a) and the signal
described by figure~\ref{fig:DecayTopologies}(d): the background events
are symmetric
while the signal events are asymmetric. The on-shell constrained
variables $M_{2CX}$, $M_{2XC}$, and $M_{2CC}$ are obtained by applying
the additional constraints of eqs.~(\ref{eq:parents-equal}) and~(\ref{eq:relatives-equal}),
which assume that the events are symmetric. Enforcing these
constraints leads to the hierarchy~(\ref{eq:inequality}) which is
simply due to the fact that a constrained minimum is larger than an
unconstrained minimum. Since the background events are symmetric, the
constraints~(\ref{eq:parents-equal}) and~(\ref{eq:relatives-equal})
will be satisfied for the \emph{true} values of the invisible momenta,
and, as long as the global minimum is not too far away (in momentum
space), one can expect a relatively mild hierarchy~(\ref{eq:inequality}). Conversely, for signal events with the mixed
event topology, the \emph{true} values of the invisible momenta in
general do not satisfy the constraints~(\ref{eq:parents-equal}) and~(\ref{eq:relatives-equal}). Thus one could expect that the effect of
imposing the constraints would be larger, leading to a larger
hierarchy~(\ref{eq:inequality}).

These intuitive considerations suggest that we look at the shift of
the on-shell constrained invariant mass variable which is caused by
the constraint itself. Keeping in mind the identity
$M_{T2}=M_{2XX}$~\cite{Cho-ml-2014naa}, we can take the usual stransverse 
mass $M_{T2}$ as our benchmark variable in the absence of any constraints.
Then, we can ``measure" the effect of the constraints by comparing
$M_{T2}$ to 
$M_{2CC}$, where both~(\ref{eq:parents-equal}) and~(\ref{eq:relatives-equal}) 
have been applied. This motivates the consideration of a new
variable\footnote{Since 
the shift $M_{2CC}- M_{T2}$ is relatively small compared to the individual 
values of $M_{2CC}$ or $M_{T2}$, in eq.~(\ref{eq:deltaM}) we prefer to define
$\Delta M$ in terms of the difference of the \emph{squared} masses.
The square root is then used merely to lower the mass dimension of $\Delta M$ back to GeV.} 
\beq
\Delta M(S;\tilde m) \equiv \sqrt{M^2_{2CC}(S;\tilde m) - M^2_{T2}(S; \tilde m)}.
\label{eq:deltaM}
\eeq 
As indicated in~(\ref{eq:deltaM}), this variable can be computed for
any of the three subsystems.
We have checked that in our example here, the $(ab)$ subsystem shows the
best discrimination between the signal and $t\bar{t}$ background. 
\begin{figure}[t]
\centering
\includegraphics[width=6.0cm]{ttbar_DiffvsM2CC_Subsystem_ab_No_Combinatorics_2}\quad 
\includegraphics[width=6.0cm]{Topo2_DiffvsM2CC_Subsystem_ab_No_Combinatorics_2}\\[3pt]
\includegraphics[width=6.0cm]{Topo1_DiffvsM2CC_Subsystem_ab_No_Combinatorics_2}\quad  
\includegraphics[width=6.0cm]{Mixed_DiffvsM2CC_Subsystem_ab_No_Combinatorics_2}

]]></p>
<fig-group><caption><p><![CDATA[Density plots exhibiting the correlations between the
  variables $\Delta M(ab)$ and $M_{T2}(ab)$ for background events
  (upper left) and for the three types of signal events: with pure
  Topology 1 from figure~\ref{fig:DecayTopologies}(b) (upper right),
  pure Topology 2 from figure~\ref{fig:DecayTopologies}(c) (lower left),
  and the mixed topology from figure~\ref{fig:DecayTopologies}(d) (lower
  right). ]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
Therefore, in figure~\ref{fig:DM}, we contrast the new variable $\Delta M(ab)$ defined
in~(\ref{eq:deltaM}) with the variable $M_{2CC}(ab)$ advocated 
above in section~\ref{sec:mixed-events}. Each panel in figure~\ref{fig:DM} 
shows a specific type of events at parton level:
$t\bar{t}$ background events (upper left panel),
signal events with pure Topology 1 from figure~\ref{fig:DecayTopologies}(b) (upper right panel),
signal events with pure Topology 2 from figure~\ref{fig:DecayTopologies}(c) (lower left panel), and
signal events with mixed topology from figure~\ref{fig:DecayTopologies}(d) (lower right panel). 
We see that, as already observed in figure~\ref{fig:M2CCcuts5},
a certain number of signal events in the mixed channel exceed the
background endpoint for $M_{2CC}$.
More importantly, the figure also shows that there are many more
signal events which do \emph{not}
exceed the background endpoint, yet their value for $\Delta M$ is
significantly larger than 
that for a typical background event. The situation does not change much if we 
account for the two-fold combinatorial ambiguity, as demonstrated by figure~\ref{fig:DMcomb}.
We conclude that $\Delta M$ possesses additional 
discriminating power, and therefore, for an optimal analysis, one should 
use \emph{both} $\Delta M$ and $M_{2CC}$ as discriminating variables.

\begin{figure}[t]
\centering
\includegraphics[width=6.0cm]{ttbar_DiffvsM2CC_Subsystem_ab_Combinatorics_2}\quad 
\includegraphics[width=6.0cm]{Topo2_DiffvsM2CC_Subsystem_ab_Combinatorics_2}\\[3pt]
\includegraphics[width=6.0cm]{Topo1_DiffvsM2CC_Subsystem_ab_Combinatorics_2}\quad 
\includegraphics[width=6.0cm]{Mixed_DiffvsM2CC_Subsystem_ab_Combinatorics_2}

]]></p>
<fig-group><caption><p><![CDATA[The same as figure~\ref{fig:DM}, but including combinatoric
  effects: each variable $M_{T2}$ and $M_{2CC}$ is separately
  calculated with the $b$-$\ell$ pairing which gives the lower value
  for the corresponding variable. ]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}




For completeness, we also present results for the $\Delta M$ variable alone.
Figure~\ref{fig:Rvariable} shows unit-normalized distributions for $\Delta M(ab)$
at the parton level (upper row) and after detector simulation and
selection cuts (lower row).
The upper left panel is done with perfect assignment for the
lepton--$b$-jet pairing,
while the upper right panel accounts for the two-fold combinatorial
ambiguity as before.
The lower right panel shows the observable total signal distribution, 
which is made up of the individual components identified on the lower
left panel.
Clearly, the variable $\Delta M(ab)$ performs quite well for signal events with 
a mixed event topology, and to some extent for signal events with Topology 1.
The effect is diluted, but still visible after detector simulation
(panels in the lower row).

\begin{figure}[t]
\centering
\includegraphics[width=7.0cm]{Diff_Subsystem_ab_No_Combinatorics_2}\quad 
\includegraphics[width=7.0cm]{Diff_Subsystem_ab_Combinatorics_2} \\[3pt]
\includegraphics[width=7.0cm]{Diff_SubSys_bl_Delphes_CMS_8TeV_byChannel_2}\quad
\includegraphics[width=7.0cm]{Diff_SubSys_bl_Delphes_CMS_8TeV_2}

]]></p>
<fig-group><caption><p><![CDATA[Unit-normalized distributions of the $\Delta M(ab)$ variable
  for signal and background events, with the same color coding scheme
  used in earlier figures.  The top row shows parton level results
  before any selection cuts, while the bottom row includes detector
  simulation and selection cuts.  The upper right (upper left) panel
  does (does not) account for the two-fold combinatorial ambiguity.
]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}


]]></p>
</sec>
<sec><title>An alternative variable: the ``relative" mass difference</title>
<p><![CDATA[
\label{sec:relative}

In the previous section, we proposed the variable, $\Delta M$,
as a measure of the effect of the constraints~(\ref{eq:parents-equal}) and~(\ref{eq:relatives-equal}). 
The idea was to look at the change in the value of $M_2$ as a result
of enforcing these constraints.
Let us now look at a different way of capturing the same effect.

Recall that as a result of the minimization involved in calculating 
the \emph{unconstrained} $M_{2XX}$ variable, one obtains values 
for the invisible momenta that minimize the maximal parent particle
invariant mass in the specified subsystem.
While these are not necessarily the true momenta of the 
invisible particles in the event, they do provide an useful ansatz 
and can be used to calculate various $3+1$-dimensional
kinematic quantities of interest~\cite{Cho-ml-2014naa}.  (See also the 
MAOS method~\cite{Cho-ml-2008tj,Park-ml-2011uz}.)
In particular, we can compute the masses of the parent particles and the 
relative particles in the event and test whether the constraints~(\ref{eq:parents-equal}) and~(\ref{eq:relatives-equal}) are satisfied or not. 
However, there is one technical complication: the function
which is being minimized in order to compute $M_{2XX}$,
sometimes has a flat direction and does not lead to a unique ansatz for 
the invisible momenta~\cite{Cho-ml-2014naa}. In order to avoid this problem, 
here we prefer to use the variable, $M_{CX}$, where the parent constraint~(\ref{eq:parents-equal}) is already applied.
Thus, we will be comparing the masses of the relative particles instead.
In analogy to~(\ref{eq:deltaM}), we therefore define
\beq
\Delta M_R(S;\tilde m) \equiv \sqrt{|M^2_{R_1}(S;\tilde m) -
  M^2_{R_2}(S; \tilde m)|}.
\label{eq:deltaMR}
\eeq 
\begin{figure}[t]
\centering
\includegraphics[width=6.0cm]{Diff_M2VsMR_tt_SubSys_ab_No_Combinatorics_2}\quad 
\includegraphics[width=6.0cm]{Diff_M2VsMR_1vs2_SubSys_ab_No_Combinatorics_2}

]]></p>
<fig-group><caption><p><![CDATA[The correlation between the parent mass difference $\Delta
  M_R$ defined in~(\ref{eq:deltaM}) and the relative mass difference
  $\Delta M$ defined in~(\ref{eq:deltaMR}) for background events
  (left) and signal events with the mixed event topology of figure~\ref{fig:DecayTopologies}(d) (right).]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
Since they both measure the same effect, namely, the impact of the relative 
constraint~(\ref{eq:relatives-equal}), we expect the two variables $\Delta M$
and $\Delta M_R$ to be correlated. This is illustrated in figure~\ref{fig:MRdiff},
where we compare $\Delta M(ab)$ and $\Delta M_R(ab)$
for background events (left panel) and signal events with the 
mixed event topology of figure~\ref{fig:DecayTopologies}(d) (right panel).
The correlation is very evident and suggests that $\Delta M_R$ can be used 
in place of $\Delta M$. The advantage of using $\Delta M_R$ is convenience:
in order to compute it, one needs to perform a single minimization
(that of the variable $M_{2CX}$), 
while to construct $\Delta M$, one needs to minimize twice: once for
$M_{T2}$ (or, equivalently,
$M_{2XX}$) and then once for $M_{2CC}$. We have also noticed that our
numerical minimization code finds the global minimum of
$M_{2CX}$ more reliably than it finds the minimum of the doubly constrained
variable $M_{2CC}$.


]]></p>
</sec>
</sec>
<sec><title>Conclusions and outlook</title>
<p><![CDATA[
\label{sec:conclusions}

The search for ``top partners," like top squarks in SUSY, will be a
key component of the LHC research program in the next run of the LHC\@. 
This is due to several reasons.
First, particles which behave like top partners are theoretically
well-motivated since they are ubiquitous in models that try to
address the hierarchy problem. 
Second, the experimental limits on third generation partners are
generally weaker, leaving room for improvement in the next run. 
Third, the signatures of top partner production typically resemble
those of SM top production, a process which will continue to be under
close scrutiny because of the intrinsic interest in the top in its own
right.

The main goal of this paper was to tackle certain \emph{difficult}
cases for stop discovery and propose new ideas for improving the
experimental sensitivity in the next LHC run.  We considered stop
signatures which led to an \emph{identical} final state as the main
irreducible top background. Of special interest to us were corners of
parameter space which would evade easy detection by normal means,
either due to small mass splittings, which lead to soft jets and leptons, or 
because the new physics signature involves real SM top quarks. 
Thus we considered the two decay topologies of figure~\ref{fig:process},
which led to the three types of signal events depicted in figure~\ref{fig:DecayTopologies}(b-d).

Given that the signal and background are so similar, discrimination is
only possible if we take full advantage of subtle kinematic
differences. This is why we focused on the recently proposed class of
on-shell constrained variables ($M_{2XX}$, $M_{2CX}$, $M_{2XC}$, and
$M_{2CC}$)~\cite{Mahbubani-ml-2012kx,Cho-ml-2014naa},
which can be suitably defined with the background event topology of \mbox{figure}~\ref{fig:DecaySubsystem} in mind (see appendix~\ref{sec:appendix}). These variables have several useful
properties which can be used for isolating the signal over the
background:
\begin{itemize}
\item \emph{Existence of upper kinematic endpoints.} While the
  background events obey the bounds~(\ref{mtab})--(\ref{mwb}), signal events may violate those bounds, depending on the new physics mass spectrum.
  Thus, by employing suitable high pass cuts on those variables, one can
  remove the majority of the background, leaving some fraction of
  the signal. In section~\ref{subsec:regions-and-study-points}, we
  analyzed the relevant mass parameter space and classified the
  regions where given kinematic endpoints for signal
  events exceed those for the background. The ``easy" regions, where
  the signal endpoints are significantly above the background
  endpoints, should be the first targets in the next LHC runs.
\item \emph{Endpoint violation in the case of the ``wrong" event
    topology.} The on-shell kinematic variables were defined with a
  specific background event topology in mind. If the signal events
  have a different event topology, either because they are asymmetric
  (e.g., the mixed event topology of figure~\ref{fig:DecayTopologies}(d)), or because they contain more
  invisible particles (as in the case of pure Topology 2 in figure~\ref{fig:DecayTopologies}(c)), they may again violate the
  background endpoints, see figures~\ref{fig:MT2cuts5} and~\ref{fig:M2CCcuts5}.
  For concreteness, in this paper we only considered the two specific
  event topologies from figure~\ref{fig:process}, but in realistic
  models, there exist other well-motivated event topologies which
  would lead to the same final state.  (A couple of such examples are
  shown in figure~\ref{fig:process2}.) The multitude of possible stop
  decay modes increases the likelihood that the signal will include
  asymmetric events, which may manifest themselves through violations
  of the expected background endpoints. 
\begin{figure}[t]
\centering
\includegraphics[scale=0.9]{OtherSignal.png}

]]></p>
<fig-group><caption><p><![CDATA[Two other possible decay topologies leading to the same final state as in figure~\ref{fig:process}.]]></p></caption></fig-group>
<p><![CDATA[
\end{figure}
\item \emph{The existence of the hierarchy} (\ref{eq:inequality})
    \emph{between the various $M_{2\sqcup\sqcup}$ variables.}
We showed that the hierarchy is relatively mild in the case of
symmetric events like the $t\bar{t}$ background and gets stronger as
the events become more asymmetric (as in the mixed topology of figure~\ref{fig:DecayTopologies}(d)).  This observation allows us to
also target signal events which are \emph{below} the background
kinematic endpoints. We believe that the related variables $\Delta M$
and $\Delta M_R$ defined in~(\ref{eq:deltaM}) and~(\ref{eq:deltaMR})
respectively, will be a useful addition to the experimenter's arsenal
of tools for new physics searches in missing energy events.
\end{itemize}

In spite of the advances proposed here, these searches remain extremely challenging.
Detector effects, jet combinatorics issues due to ISR and FSR,
and $b$-jet misidentification contribute to the degradation of the 
parton level significance.\footnote{For a more complete discussion, see~\cite{Kim-ml-2014ana}. } Nevertheless, we believe that the techniques presented here 
will prove useful in searches for top partners at the~LHC\@. 

In this paper, we have employed a simplified model approach as shown in figure~\ref{fig:DecayTopologies} in order to best make contact with experimental efforts.
Of course, when interpreting such simplified model experimental limits or discoveries in 
terms of some complete theory, one must compute both signal~\cite{Beenakker-ml-1997ut,Beenakker-ml-2010nq,Beenakker-ml-2011fu}
and background~\cite{Cacciari-ml-2011hy,Czakon-ml-2011xx,Czakon-ml-2012zr,Czakon-ml-2012pz}
production cross-sections and the relevant branching ratios~\cite{Djouadi-ml-2006bz} to a high degree of precision. 

The on-shell constrained variables are suitable generalizations of the
stransverse mass variable, $M_{T2}$, which is being used extensively in
experimental searches, and for which several public codes exist.
In contrast, there is no public code which allows the computation of
the $M_2$ variables.  We are developing such a code for public
release in the near future to facilitate the wider use of $M_2$ variables~\cite{newpaper}.
 

]]></p>
</sec>
</body>
<back>
<app-group>


<app><sec><title>The complete set of \texorpdfstring{\boldmath $M_2$}{M(2)} variables for the \texorpdfstring{\boldmath $t\bar{t}$}{t bar(t)} event topology</title>
<p><![CDATA[
\label{sec:appendix}

In this appendix, we collect the specific definitions of the fifteen
$M_2$ variables used in this paper. We consider the three subsystems
in figure~\ref{fig:DecaySubsystem} as applied to the $t\bar{t}$
background events of figure~\ref{fig:DecayTopologies}(a). We write the
equations in terms of the measured 4-momenta of the $b$ and $\bar{b}$
quarks, $p_b$ and $p_{\bar{b}}$, and the measured 4-momenta of the
lepton and antilepton, $p_{\ell^-}$ and $p_{\ell^+}$.  The invisible
neutrino 4-momenta will be denoted by $q_1$ and $q_2$, where we take
``1'' to refer to the decay chain initiated by a top quark and ``2''
to refer to the decay chain initiated by an anti-top. In each
subsystem, the test mass $\tilde m$ is taken to be the corresponding
true daughter mass.

In the $(ab)$ subsystem the daughter particles are the two neutrinos,
with mass $\tilde m=0$.
The definitions~(\ref{M2CCdef})--(\ref{MT2def}) imply
\begin{align}
M_{2CC}^2(ab;\tilde m = 0) &\equiv 
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\
q_1^2=q_2^2=\tilde m^2=0\\
\vec{q}_{1T}+\vec{q}_{2T} = \mpt  \\
(p_{b}+p_{\ell^+}+q_1)^2 =(p_{\bar b}+p_{\ell^-}+q_2)^2 \\
(p_{\ell^+}+q_1)^2 =(p_{\ell^-}+q_2)^2 }}
\left\{(p_{b}+p_{\ell^+}+q_1)^2 \right\},
\label{eq:m2CCab}
\\ 
M_{2CX}^2(ab;\tilde m = 0) &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\
q_1^2=q_2^2=\tilde m^2=0\\
\vec{q}_{1T}+\vec{q}_{2T} = \mpt  \\
(p_{b}+p_{\ell^+}+q_1)^2 =(p_{\bar b}+p_{\ell^-}+q_2)^2 }}
\left\{(p_{b}+p_{\ell^+}+q_1)^2 \right\},
\label{eq:m2CXab}
\\ 
M_{2XC}^2(ab;\tilde m = 0) &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\
q_1^2=q_2^2=\tilde m^2=0\\
\vec{q}_{1T}+\vec{q}_{2T} = \mpt  \\
(p_{\ell^+}+q_1)^2 =(p_{\ell^-}+q_2)^2 }}
\!\!\!\!\!\!\!\!\!
\left\{
\max\left[(p_{b}+p_{\ell^+}+q_1)^2\;, (p_{\bar b}+p_{\ell^-}+q_2)^2
\right]\right\},
\label{eq:m2XCab}
\\ 
M_{2XX}^2(ab;\tilde m = 0) &\equiv 
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\
q_1^2=q_2^2=\tilde m^2=0\\
\vec{q}_{1T}+\vec{q}_{2T} = \mpt }}
\!\!\!\!\!
\left\{
\max\left[(p_{b}+p_{\ell^+}+q_1)^2\;, (p_{\bar b}+p_{\ell^-}+q_2)^2
\right]\right\},
\label{eq:m2XXab}
\\ 
M_{T2}^2(ab;\tilde m = 0) &\equiv
\min_{\substack{\vec{q}_{1T},\vec{q}_{2T}\\
\vec{q}_{1T}+\vec{q}_{2T} =\mpt }}
\!\!\!\!\!\!\!
\left\{\max\left[M^2_{Tt}(\vec{q}_{1T},\tilde m=0),\;
              M^2_{T\bar{t}}(\vec{q}_{2T},\tilde m=0)\right] \right\}.
\label{eq:mT2ab}
\end{align}
In the last equation $M_{Tt}$ ($M_{T\bar{t}}$) is the transverse mass of the 
hypothesized top quark (anti-top quark):
\begin{align}
M^2_{Tt}(\vec{q}_{1T},\tilde m) &\equiv
 (E_{Tb}+E_{T\ell^+}+E_{T1})^2 -
 (\vec{p}_{Tb}+\vec{p}_{T\ell^+}+\vec{q}_{1T})^2,
\label{MTt} \\ 
M^2_{T\bar{t}}(\vec{q}_{2T},\tilde m) &\equiv
 (E_{T\bar{b}}+E_{T\ell^-}+E_{T2})^2 -
 (\vec{p}_{T\bar{b}}+\vec{p}_{T\ell^-}+\vec{q}_{2T})^2,
 \label{MTt-bar} 
\end{align}
with the transverse energies defined as usual,
e.g.~$E_{Ti}=\sqrt{\tilde m^2 + \vec{q}_{Ti}^{\, 2}}$.

In the $(b)$ subsystem the daughter particles are again the massless
neutrinos, but this time we minimize the masses of the hypothesized
$W^\pm$ particles:
\begin{align}
M_{2CC}^2(b;\tilde m = 0) &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\
q_1^2=q_2^2=\tilde m^2=0\\
\vec{q}_{1T}+\vec{q}_{2T} = \mpt  \\
(p_{b}+p_{\ell^+}+q_1)^2 =(p_{\bar b}+p_{\ell^-}+q_2)^2 \\
(p_{\ell^+}+q_1)^2 =(p_{\ell^-}+q_2)^2 }}
\left\{(p_{\ell^+}+q_1)^2 \right\},
\label{eq:m2CCb}
\\ 
M_{2CX}^2(b;\tilde m = 0) &\equiv 
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\
q_1^2=q_2^2=\tilde m^2=0\\
\vec{q}_{1T}+\vec{q}_{2T} = \mpt  \\
(p_{\ell^+}+q_1)^2 =(p_{\ell^-}+q_2)^2 }}
\left\{(p_{\ell^+}+q_1)^2 \right\},
\label{eq:m2CXb}
\\ 
M_{2XC}^2(b;\tilde m = 0) &\equiv 
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\
q_1^2=q_2^2=\tilde m^2=0\\
\vec{q}_{1T}+\vec{q}_{2T} = \mpt  \\
(p_{b}+p_{\ell^+}+q_1)^2 =(p_{\bar b}+p_{\ell^-}+q_2)^2 }}
\!\!\!\!\!\!\!\!\!
\left\{
\max\left[(p_{\ell^+}+q_1)^2\;, (p_{\ell^-}+q_2)^2 \right]\right\},
\label{eq:m2XCb}
\\ 
M_{2XX}^2(b;\tilde m = 0) &\equiv
\min_{\substack{\vec{q}_{1},\vec{q}_{2}\\
q_1^2=q_2^2=\tilde m^2=0\\
\vec{q}_{1T}+\vec{q}_{2T} = \mpt }}
\left\{
\max\left[(p_{\ell^+}+q_1)^2\;, (p_{\ell^-}+q_2)^2 \right]\right\},
\label{eq:m2XXb}
\\ 
M_{T2}^2(b;\tilde m = 0) &\equiv
\min_{\substack{\vec{q}_{1T},\vec{q}_{2T}\\
\vec{q}_{1T}+\vec{q}_{2T} =\mpt }}
\!\!\!\!
\left\{\max\left[M^2_{TW^+}(\vec{q}_{1T},\tilde m=0),\;
              M^2_{TW^-}(\vec{q}_{2T},\tilde m=0)\right] \right\},
\label{eq:mT2b}
\end{align}
with
\begin{align}
M^2_{TW^+}(\vec{q}_{1T},\tilde m) &\equiv
 (E_{T\ell^+}+E_{T1})^2 - (\vec{p}_{T\ell^+}+\vec{q}_{1T})^2,
 \\ 
M^2_{TW^-}(\vec{q}_{2T},\tilde m) &\equiv
 (E_{T\ell^-}+E_{T2})^2 - (\vec{p}_{T\ell^-}+\vec{q}_{2T})^2.
\end{align}

Finally, consider the $(a)$ subsystem, in which the the parents are
the top quarks,
while the daughters are the $W^\pm$ bosons, with masses $\tilde m=m_W$. 
Denoting now the 4-momenta of the two hypothesized 
$W^\pm$ bosons with $q_{W^+}$ and $q_{W^-}$, we have 
{\allowdisplaybreaks
\begin{align}
\hskip-20pt M_{2CC}^2(a;\tilde m = m_W) &\equiv 
\min_{\substack{\vec{q}_{W^+},\vec{q}_{W^-}\\
q_{W^+}^2=q_{W^-}^2=\tilde m^2=m_W^2\\
\vec{q}_{TW^+}+\vec{q}_{TW^-} = \mpt +\vec{p}_{T\ell^+} +\vec{p}_{T\ell^-}  \\
(p_{b}+q_{W^+})^2 =(p_{\bar b}+q_{W^-})^2 \\
(q_{W^+}-p_{\ell^+})^2 =(q_{W^-}-p_{\ell^-})^2 }}
\left\{(p_{b}+q_{W^+})^2 \right\},
\label{eq:m2CCa}
\\ 
M_{2CX}^2(a;\tilde m = m_W) &\equiv
\min_{\substack{\vec{q}_{W^+},\vec{q}_{W^-}\\
q_{W^+}^2=q_{W^-}^2=\tilde m^2=m_W^2\\
\vec{q}_{TW^+}+\vec{q}_{TW^-} = \mpt +\vec{p}_{T\ell^+}
+\vec{p}_{T\ell^-}  \\
(p_{b}+q_{W^+})^2 =(p_{\bar b}+q_{W^-})^2 }}
\left\{(p_{b}+q_{W^+})^2 \right\},
\label{eq:m2CXa}
\\ 
M_{2XC}^2(a;\tilde m = m_W) &\equiv\!\!
\min_{\substack{\vec{q}_{W^+},\vec{q}_{W^-}\\
q_{W^+}^2=q_{W^-}^2=\tilde m^2=m_W^2\\
\vec{q}_{TW^+}+\vec{q}_{TW^-} = \mpt +\vec{p}_{T\ell^+}
+\vec{p}_{T\ell^-}  \\
(q_{W^+}-p_{\ell^+})^2 =(q_{W^-}-p_{\ell^-})^2 }}
\!\!\!\!\!\!\!\!\!
\left\{
\max\left[(p_{b}+q_{W^+})^2\;, (p_{\bar b}+q_{W^-})^2 \right]\right\},
\label{eq:m2XCa}
\\ 
M_{2XX}^2(a;\tilde m = m_W) &\equiv\!\!
\min_{\substack{\vec{q}_{W^+},\vec{q}_{W^-}\\
q_{W^+}^2=q_{W^-}^2=\tilde m^2=m_W^2\\
\vec{q}_{TW^+}+\vec{q}_{TW^-} = \mpt +\vec{p}_{T\ell^+} +\vec{p}_{T\ell^-} }}
\!\!\!\!\!\!\!\!\!
\left\{
\max\left[(p_{b}+q_{W^+})^2\;, (p_{\bar b}+q_{W^-})^2 \right]\right\},
\label{eq:m2XXa}
\\ 
M_{T2}^2(a;\tilde m \!=\! m_W) &\equiv\!\!\!\!\!\!\!\!\!
\min_{\substack{\vec{q}_{W^+},\vec{q}_{W^-}\\
\vec{q}_{TW^+}+\vec{q}_{TW^-} = \mpt +\vec{p}_{T\ell^+} +\vec{p}_{T\ell^-} }}
\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!\!
\left\{\max\left[M^2_{Tt}(\vec{q}_{TW^+},m_W),\;
              M^2_{T\bar{t}}(\vec{q}_{TW^-},m_W)\right] \right\},
\label{eq:mT2a}
\end{align}
where now the top quark transverse masses~(\ref{MTt})--(\ref{MTt-bar})
are rewritten in terms of the $W^\pm$ boson momenta, $q_{W^+}$ and
$q_{W^-}$, as
\begin{align}
M^2_{Tt}(\vec{q}_{TW^+},\tilde m=m_W) &\equiv
 (E_{Tb}+E_{TW^+})^2 - (\vec{p}_{Tb}+\vec{q}_{TW^+})^2, \\ 
M^2_{T\bar{t}}(\vec{q}_{TW^-},\tilde m=m_W) &\equiv
 (E_{T\bar{b}}+E_{TW^-})^2 - (\vec{p}_{T\bar{b}}+\vec{q}_{TW^-})^2,
\end{align}}with $E_{TW^\pm}=\sqrt{ m_W^2 + \vec{q}_{TW^\pm}^{\, 2}}$.


]]></p>
</sec>
</app>
<app><sec><ack><title>Acknowledgments</title>
<p><![CDATA[

We thank Kaustubh Agashe, Roberto Franceschini, and Kyoungchul Kong for
useful discussions.  
WC acknowledges support by Project Code (IBS-R018-D1);
WC and DK acknowledge support by LHC-TI
postdoctoral fellowships under grant NSF-PHY-0969510.
WC, JG, DK, and KM acknowledge support from U.S. Department of Energy
Grant ER41990.
JG, KM, FM, and LP thank their CMS colleagues for useful discussions.
MP was supported by World Premier International Research Center
Initiative (WPI), MEXT, Japan. MP also acknowledges the
Max-Planck-Gesellschaft, the Korea Ministry of Education, Science and
Technology, Gyeongsangbuk-Do and Pohang City for the support of the
Independent Junior Research Group at the APCTP.


]]></p>
</ack></sec>
</app>
</app-group>
<ref-list>
 <ref id="Martin-ml-1997ns"> <label>1</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1142/9789814307505_0001</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2012tfa"> <label>2</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.physletb.2012.08.020</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2012ufa"> <label>3</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.physletb.2012.08.021</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2013wta"> <label>4</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP10(2013)130</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2014pda"> <label>5</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP06(2014)035</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2014wea"> <label>6</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP09(2014)176</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2014bva"> <label>7</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP09(2014)015</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2013wxa"> <label>8</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.physletb.2013.06.058</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2013fea"> <label>9</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP01(2014)163</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2014lfa"> <label>10</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP06(2014)055</uri>
  </element-citation>
 </ref>

 <ref id="Craig-ml-2013cxa"> <label>11</label>
  <element-citation>
   <uri>http://arxiv.org/abs/1309.0528</uri>
  </element-citation>
 </ref>

 <ref id="Altarelli-ml-1984ve"> <label>12</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/0370-2693(84)90190-4</uri>
  </element-citation>
 </ref>

 <ref id="Bigi-ml-1985aq"> <label>13</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/0370-2693(85)90560-X</uri>
  </element-citation>
 </ref>

 <ref id="Hikasa-ml-1987db"> <label>14</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.36.724</uri>
  </element-citation>
 </ref>

 <ref id="Baer-ml-1991cb"> <label>15</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.44.725</uri>
  </element-citation>
 </ref>

 <ref id="Baer-ml-1994xr"> <label>16</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.50.4517</uri>
  </element-citation>
 </ref>

 <ref id="Lopez-ml-1994zw"> <label>17</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1142/S0217732395002441</uri>
  </element-citation>
 </ref>

 <ref id="Chou-ml-1999zb"> <label>18</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.61.055004</uri>
  </element-citation>
 </ref>

 <ref id="Demina-ml-1999ty"> <label>19</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.62.035011</uri>
  </element-citation>
 </ref>

 <ref id="ADLO-ml-2004aa"> <label>20</label>
  <mixed-citation>
   LEP2 SUSY Working Group, <italic>Combined LEP stop and sbottom Results
183–208GeV </italic>, Note LEPSUSYWG/04-02.1 (2004).

  </mixed-citation>
 </ref>

 <ref id="Heister-ml-2002hp"> <label>21</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/S0370-2693(02)01827-0</uri>
  </element-citation>
 </ref>

 <ref id="Abdallah-ml-2003xe"> <label>22</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1140/epjc/s2003-01355-5</uri>
  </element-citation>
 </ref>

 <ref id="Achard-ml-2003ge"> <label>23</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.physletb.2003.10.010</uri>
  </element-citation>
 </ref>

 <ref id="Abbiendi-ml-2002mp"> <label>24</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/S0370-2693(02)02593-5</uri>
  </element-citation>
 </ref>

 <ref id="Jaffre-ml-2012gx"> <label>25</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1051/epjconf/20122809006</uri>
  </element-citation>
 </ref>

 <ref id="Aaltonen-ml-2010uf"> <label>26</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.82.092001</uri>
  </element-citation>
 </ref>

 <ref id="Mackin-ml-2010zza"> <label>27</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1063/1.3327567</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2013ija"> <label>28</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP10(2013)189</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2014qaa"> <label>29</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP06(2014)124</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2014mha"> <label>30</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1140/epjc/s10052-014-2883-6</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2014kva"> <label>31</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1140/epjc/s10052-014-3109-7</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2014kra"> <label>32</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP11(2014)118</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2014lra"> <label>33</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP10(2014)024</uri>
  </element-citation>
 </ref>

 <ref id="Aad-ml-2014nra"> <label>34</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.90.052008</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2012paa"> <label>35</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP03(2013)037</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2013lya"> <label>36</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1140/epjc/s10052-013-2568-6</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2013xsw"> <label>37</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevLett.111.221801</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2013xna"> <label>38</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1140/epjc/s10052-013-2677-2</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2013iqa"> <label>39</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.physletb.2014.04.023</uri>
  </element-citation>
 </ref>

 <ref id="Chatrchyan-ml-2013mya"> <label>40</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevLett.112.161802</uri>
  </element-citation>
 </ref>

 <ref id="Khachatryan-ml-2014doa"> <label>41</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.physletb.2014.07.053</uri>
  </element-citation>
 </ref>

 <ref id="Feng-ml-1998iq"> <label>42</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/S0550-3213(99)00026-7</uri>
  </element-citation>
 </ref>

 <ref id="Dimopoulos-ml-1995mi"> <label>43</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/0370-2693(95)00961-J</uri>
  </element-citation>
 </ref>

 <ref id="Pomarol-ml-1995xc"> <label>44</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/0550-3213(96)00074-0</uri>
  </element-citation>
 </ref>

 <ref id="Cohen-ml-1996vb"> <label>45</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/S0370-2693(96)01183-5</uri>
  </element-citation>
 </ref>

 <ref id="Perelstein-ml-2007nx"> <label>46</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2007/04/070</uri>
  </element-citation>
 </ref>

 <ref id="Kong-ml-2014jwa"> <label>47</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.89.074007</uri>
  </element-citation>
 </ref>

 <ref id="Han-ml-2008gy"> <label>48</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2009/05/117</uri>
  </element-citation>
 </ref>

 <ref id="Carena-ml-2008mj"> <label>49</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2008/10/109</uri>
  </element-citation>
 </ref>

 <ref id="Perelstein-ml-2008zt"> <label>50</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2009/03/141</uri>
  </element-citation>
 </ref>

 <ref id="Plehn-ml-2010st"> <label>51</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP10(2010)078</uri>
  </element-citation>
 </ref>

 <ref id="Bornhauser-ml-2010mw"> <label>52</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.83.035008</uri>
  </element-citation>
 </ref>

 <ref id="Bi-ml-2011ha"> <label>53</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.85.035005</uri>
  </element-citation>
 </ref>

 <ref id="He-ml-2011tp"> <label>54</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP05(2012)148</uri>
  </element-citation>
 </ref>

 <ref id="Drees-ml-2012dd"> <label>55</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.86.035024</uri>
  </element-citation>
 </ref>

 <ref id="Bai-ml-2012gs"> <label>56</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP07(2012)110</uri>
  </element-citation>
 </ref>

 <ref id="Plehn-ml-2012pr"> <label>57</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP08(2012)091</uri>
  </element-citation>
 </ref>

 <ref id="Alves-ml-2012ft"> <label>58</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.87.035016</uri>
  </element-citation>
 </ref>

 <ref id="Han-ml-2012fw"> <label>59</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP08(2012)083</uri>
  </element-citation>
 </ref>

 <ref id="Kaplan-ml-2012gd"> <label>60</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP07(2012)119</uri>
  </element-citation>
 </ref>

 <ref id="Berger-ml-2012an"> <label>61</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevLett.109.152004</uri>
  </element-citation>
 </ref>

 <ref id="Ghosh-ml-2012ud"> <label>62</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevLett.110.141801</uri>
  </element-citation>
 </ref>

 <ref id="Chen-ml-2012uw"> <label>63</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP11(2012)124</uri>
  </element-citation>
 </ref>

 <ref id="Kilic-ml-2012kw"> <label>64</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP04(2013)110</uri>
  </element-citation>
 </ref>

 <ref id="Graesser-ml-2012qy"> <label>65</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevLett.111.121802</uri>
  </element-citation>
 </ref>

 <ref id="Krizka-ml-2012ah"> <label>66</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.87.095016</uri>
  </element-citation>
 </ref>

 <ref id="Delgado-ml-2012eu"> <label>67</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1140/epjc/s10052-013-2370-5</uri>
  </element-citation>
 </ref>

 <ref id="Dutta-ml-2013sta"> <label>68</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.87.095007</uri>
  </element-citation>
 </ref>

 <ref id="Buckley-ml-2013lpa"> <label>69</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP08(2013)086</uri>
  </element-citation>
 </ref>

 <ref id="Chakraborty-ml-2013moa"> <label>70</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP10(2013)122</uri>
  </element-citation>
 </ref>

 <ref id="Low-ml-2013aza"> <label>71</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.88.095018</uri>
  </element-citation>
 </ref>

 <ref id="Bai-ml-2013ema"> <label>72</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP08(2013)085</uri>
  </element-citation>
 </ref>

 <ref id="Belanger-ml-2013oka"> <label>73</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.89.015003</uri>
  </element-citation>
 </ref>

 <ref id="Dutta-ml-2013gga"> <label>74</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.90.095022</uri>
  </element-citation>
 </ref>

 <ref id="Buckley-ml-2014fqa"> <label>75</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.90.014046</uri>
  </element-citation>
 </ref>

 <ref id="Ortiz-ml-2014iza"> <label>76</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.90.075009</uri>
  </element-citation>
 </ref>

 <ref id="Czakon-ml-2014fka"> <label>77</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevLett.113.201803</uri>
  </element-citation>
 </ref>

 <ref id="Ismail-ml-2014cma"> <label>78</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.91.074002</uri>
  </element-citation>
 </ref>

 <ref id="Eifert-ml-2014kea"> <label>79</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.physletb.2015.02.039</uri>
  </element-citation>
 </ref>

 <ref id="Datta-ml-2011ef"> <label>80</label>
  <element-citation>
   <uri>http://arxiv.org/abs/1111.0200</uri>
  </element-citation>
 </ref>

 <ref id="Agashe-ml-2010gt"> <label>81</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.82.015007</uri>
  </element-citation>
 </ref>

 <ref id="Agashe-ml-2010tu"> <label>82</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.84.055020</uri>
  </element-citation>
 </ref>

 <ref id="Giudice-ml-2011ib"> <label>83</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.85.075019</uri>
  </element-citation>
 </ref>

 <ref id="Cho-ml-2012er"> <label>84</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevLett.112.211801</uri>
  </element-citation>
 </ref>

 <ref id="Barr-ml-2011xt"> <label>85</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.84.095031</uri>
  </element-citation>
 </ref>

 <ref id="Mahbubani-ml-2012kx"> <label>86</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP03(2013)134</uri>
  </element-citation>
 </ref>

 <ref id="Cho-ml-2014naa"> <label>87</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP08(2014)070</uri>
  </element-citation>
 </ref>

 <ref id="Lester-ml-1999tx"> <label>88</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/S0370-2693(99)00945-4</uri>
  </element-citation>
 </ref>

 <ref id="Barr-ml-2003rg"> <label>89</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/0954-3899/29/10/304</uri>
  </element-citation>
 </ref>

 <ref id="Barr-ml-2009wu"> <label>90</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.80.074007</uri>
  </element-citation>
 </ref>

 <ref id="Allanach-ml-2011ej"> <label>91</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP07(2011)104</uri>
  </element-citation>
 </ref>

 <ref id="Murayama-ml-2011hj"> <label>92</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.84.094015</uri>
  </element-citation>
 </ref>

 <ref id="Burns-ml-2008va"> <label>93</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2009/03/143</uri>
  </element-citation>
 </ref>

 <ref id="Ross-ml-2007rm"> <label>94</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.physletb.2008.06.003</uri>
  </element-citation>
 </ref>

 <ref id="Cho-ml-2009ve"> <label>95</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.81.095010</uri>
  </element-citation>
 </ref>

 <ref id="Barr-ml-2009jv"> <label>96</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2009/11/096</uri>
  </element-citation>
 </ref>

 <ref id="Konar-ml-2009qr"> <label>97</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP04(2010)086</uri>
  </element-citation>
 </ref>

 <ref id="Cho-ml-2008cu"> <label>98</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.78.034019</uri>
  </element-citation>
 </ref>

 <ref id="Smith-ml-1983aa"> <label>99</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevLett.50.1738</uri>
  </element-citation>
 </ref>

 <ref id="Barger-ml-1983wf"> <label>100</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/BF01648783</uri>
  </element-citation>
 </ref>

 <ref id="Cho-ml-2007qv"> <label>101</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevLett.100.171801</uri>
  </element-citation>
 </ref>

 <ref id="Cho-ml-2007dh"> <label>102</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2008/02/035</uri>
  </element-citation>
 </ref>

 <ref id="Konar-ml-2008ei"> <label>103</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2009/03/085</uri>
  </element-citation>
 </ref>

 <ref id="Alwall-ml-2014hca"> <label>104</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP07(2014)079</uri>
  </element-citation>
 </ref>

 <ref id="Ball-ml-2012cx"> <label>105</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.nuclphysb.2012.10.003</uri>
  </element-citation>
 </ref>

 <ref id="Sjostrand-ml-2006za"> <label>106</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2006/05/026</uri>
  </element-citation>
 </ref>

 <ref id="deFavereau-ml-2013fsa"> <label>107</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP02(2014)057</uri>
  </element-citation>
 </ref>

 <ref id="Cacciari-ml-2008gp"> <label>108</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1088/1126-6708/2008/04/063</uri>
  </element-citation>
 </ref>

 <ref id="Cho-ml-2008tj"> <label>109</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.79.031701</uri>
  </element-citation>
 </ref>

 <ref id="Park-ml-2011uz"> <label>110</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1103/PhysRevD.84.096001</uri>
  </element-citation>
 </ref>

 <ref id="Kim-ml-2014ana"> <label>111</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP03(2015)134</uri>
  </element-citation>
 </ref>

 <ref id="Beenakker-ml-1997ut"> <label>112</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/S0550-3213(98)00014-5</uri>
  </element-citation>
 </ref>

 <ref id="Beenakker-ml-2010nq"> <label>113</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP08(2010)098</uri>
  </element-citation>
 </ref>

 <ref id="Beenakker-ml-2011fu"> <label>114</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1142/S0217751X11053560</uri>
  </element-citation>
 </ref>

 <ref id="Cacciari-ml-2011hy"> <label>115</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.physletb.2012.03.013</uri>
  </element-citation>
 </ref>

 <ref id="Czakon-ml-2011xx"> <label>116</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1016/j.cpc.2014.06.021</uri>
  </element-citation>
 </ref>

 <ref id="Czakon-ml-2012zr"> <label>117</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP12(2012)054</uri>
  </element-citation>
 </ref>

 <ref id="Czakon-ml-2012pz"> <label>118</label>
  <element-citation>
   <uri>http://dx.doi.org/10.1007/JHEP01(2013)080</uri>
  </element-citation>
 </ref>

 <ref id="Djouadi-ml-2006bz"> <label>119</label>
  <element-citation>
   <uri>http://arxiv.org/abs/hep-ph/0609292</uri>
  </element-citation>
 </ref>

 <ref id="newpaper"> <label>120</label>
  <mixed-citation>
   W.S. Cho et al., <italic>A computer package for computing constrained invariant mass
variables</italic>, in preparation.

  </mixed-citation>
 </ref>

</ref-list></back>
</article>
