1 |
claudioc |
1.1 |
\section{Overview and Analysis Strategy}
|
2 |
|
|
\label{sec:overview}
|
3 |
|
|
|
4 |
|
|
We are searching for a $t\bar{t}\chi^0\chi^0$ or $W \ell b W \ell \bar{b} \chi^0 \chi^0$ final state
|
5 |
|
|
(after top decay in the first mode, the final states are actually the same). So to first order
|
6 |
|
|
this is ``$t\bar{t} +$ extra \met''.
|
7 |
|
|
|
8 |
|
|
We work in the $\ell +$ jets final state, where the main background is $t\bar{t}$. We look for
|
9 |
|
|
\met inconsistent with $W \to \ell \nu$. We do this by concentrating on the $\ell \nu$ transverse
|
10 |
|
|
mass ($M_T$), since except for resolution effects, $M_T < M_W$ for $W \to \ell \nu$. Thus, the
|
11 |
|
|
initial analysis is simply a counting experiment in the tail of the $M_T$ distribution.
|
12 |
|
|
|
13 |
|
|
The event selection is one-and-only-one high $P_T$ isolated lepton, four or more jets, and
|
14 |
|
|
some moderate \met cut. At least one of the jets has to be btagged to reduce $W+$ jets.
|
15 |
|
|
The event sample is then dominated by $t\bar{t}$, but there are also contributions from $W+$ jets,
|
16 |
|
|
single top, dibosons, etc.
|
17 |
|
|
|
18 |
|
|
In order to be sensitive to $\widetilde{t}\widetilde{t}$ production, the background in the $M_T$
|
19 |
|
|
tail has to be controlled at the level of 10\% or better. So this is (almost) a precision measurement.
|
20 |
|
|
|
21 |
|
|
The $t\bar{t}$ events in the $M_T$ tail can be broken up into two categories:
|
22 |
|
|
(i) $t\bar{t} \to \ell $+ jets and (ii) $t\bar{t} \to \ell^+ \ell^-$ where one of the two
|
23 |
|
|
leptons is not found by the second-lepton-veto (here the second lepton can be a hadronically
|
24 |
|
|
decaying $\tau$).
|
25 |
|
|
For a reasonable $M_T$ cut, say $M_T >$ 150 GeV, the dilepton background is of order 80\% of
|
26 |
|
|
the total. This is because in dileptons there are two neutrinos from $W$ decay, thus $M_T$
|
27 |
|
|
is not bounded by $M_W$. This is a very important point: while it is true that we are looking in
|
28 |
|
|
the tail of $M_T$, the bulk of the background events end up there not because of some exotic
|
29 |
|
|
\met reconstruction failure, but because of well understood physics processes. This means that
|
30 |
|
|
the background estimate can be taken from Monte Carlo (MC), after carefully accounting for possible
|
31 |
|
|
data/MC differences. Sophisticated fully ``data driven'' techniques are not really needed.
|
32 |
|
|
|
33 |
|
|
Another important point is that in order to minimize systematic uncertainties, the MC background
|
34 |
|
|
predictions are always normalized to the bulk of the $t\bar{t}$ data, ie, events passing all of the
|
35 |
|
|
requirements but with $M_T \approx 80$ GeV.
|
36 |
|
|
This removes uncertainties
|
37 |
|
|
due to $\sigma(t\bar{t})$, lepton ID, trigger efficiency, luminosity, etc.
|
38 |
|
|
|
39 |
|
|
The $\ell +$ jets background, which is dominated by
|
40 |
|
|
$t\bar{t} \to \ell $+ jets, but also includes some $W +$ jets as well as single top,
|
41 |
|
|
is estimated as follows:
|
42 |
|
|
\begin{enumerate}
|
43 |
|
|
\item We select a control sample of events passing all cuts, but anti-btagged. This is
|
44 |
|
|
sample is now dominated by $W +$ jets. The sample is used to understand the
|
45 |
|
|
$M_T$ tail in $\ell +$ jets processes.
|
46 |
|
|
\item In MC we measure the ratio of the number of $\ell +$ jets events in the $M_T$ tail to
|
47 |
|
|
the number of events with $M_T \approx$ 80 GeV. This ratio turns out to be pretty much the
|
48 |
|
|
same for all sources of $\ell +$ jets.
|
49 |
|
|
\item In data we measure the same ratio but after correcting for the $t\bar{t} \to$ dilepton
|
50 |
|
|
contribution, as well as dibosons etc. The dilepton contribution is taken from MC after
|
51 |
|
|
the correction described below.
|
52 |
|
|
\item We compare the two ratios, as well as the shapes of the data and MC $M_T$ distributions.
|
53 |
|
|
If they do not agree, we try to figure out why and fix it. If they agree well enough, we define a
|
54 |
|
|
data MC scale factor (SF) which is the ratio of the ratios defined in step 2 and 3, keeping track of the
|
55 |
|
|
uncertainty.
|
56 |
|
|
\item We next perform the full selection in $t\bar{t} \to \ell +$ jets MC, and measure this ratio
|
57 |
|
|
again (which should be the same as that in step 2).
|
58 |
|
|
\item We perform the full selection in data. We count the events with $M_T \approx 80$ GeV, we
|
59 |
|
|
subtract off the dilepton contribution, we multiply the subtracted event count by the ratio from step 5 (or from
|
60 |
|
|
step 2), and also by the data/MC SF from step 4. The result is the prediction for the $\ell +$ jets BG in
|
61 |
|
|
the $M_T$ tail.
|
62 |
|
|
\end{enumerate}
|
63 |
|
|
|
64 |
|
|
The dilepton background can be broken up into many components depending
|
65 |
|
|
on the characteristics of the 2nd (undetected) lepton
|
66 |
|
|
\begin{itemize}
|
67 |
|
|
\item 3-prong hadronic tau decay
|
68 |
|
|
\item 1-prong hadronic tau decay
|
69 |
|
|
\item $e$ or $\mu$ possibly from $\tau$ decay
|
70 |
|
|
\end{itemize}
|
71 |
|
|
We have currently no veto against 3-prong taus. For the other two categories, we explicitely
|
72 |
|
|
veto events with additional electrons and muons above 10 GeV , and
|
73 |
|
|
we veto events with an isolated track of $P_T > 10$ GeV. This also rejects 1-prong taus
|
74 |
|
|
(it turns out that the explicit $e$ or $\mu$ veto is redundant with the isolated track veto).
|
75 |
|
|
Therefore the latter two categories can be broken into
|
76 |
|
|
\begin{itemize}
|
77 |
|
|
\item out of acceptance $(|\eta| > 2.50)$
|
78 |
|
|
\item $P_T < 10$ GeV
|
79 |
|
|
\item $P_T > 10$ GeV, but survives the additional lepton/track isolation veto
|
80 |
|
|
\end{itemize}
|
81 |
|
|
Monte Carlo studies indicate that there is no dominant contribution: it is ``a little bit of this,
|
82 |
|
|
and a little bit of that''.
|
83 |
|
|
|
84 |
|
|
The high $M_T$ dilepton backgrounds come from MC, but their rate is normalized to the
|
85 |
|
|
$M_T \approx 80$ GeV peak. In other to perform this normalization in data, the $W +$ jets
|
86 |
|
|
events in the $M_T$ peak have to be subtracted off. This introduces a systematic uncertainty.
|
87 |
|
|
|
88 |
|
|
There are two types of effects that can influence the MC dilepton prediction: physics effects
|
89 |
|
|
and instrumental effects. We discuss these next, starting from physics.
|
90 |
|
|
|
91 |
|
|
First of all, many of our $t\bar{t}$ MC samples (eg: MadGraph) have
|
92 |
|
|
BR$(W \to \ell \nu)=\frac{1}{9} = 0.1111$.
|
93 |
|
|
PDG says BR$(W \to \ell \nu) = 0.1080 \pm 0.0009$. This difference matter, so the $t\bar{t}$ MC
|
94 |
|
|
must be corrected to account for this.
|
95 |
|
|
|
96 |
|
|
Second, our selection is $\ell +$ 4 or more jets. A dilepton event passes the selection only if there are
|
97 |
|
|
two additional jet from ISR, or one jet from ISR and one jet which is reconstructed from the
|
98 |
|
|
unidentified lepton, {\it e.g.}, a three-prong tau. Therefore, all MC dilepton $t\bar{t}$ samples used
|
99 |
|
|
in the analysis must have their jet multiplicity corrected (if necessary) to agree with what is
|
100 |
|
|
seen in $t\bar{t}$ data. We use a data control sample of well identified dilepton events with
|
101 |
|
|
\met and at least two jets as a template to ``adjust'' the $N_{jet}$ distribution of the $t\bar{t} \to$
|
102 |
|
|
dileptons MC samples.
|
103 |
|
|
|
104 |
|
|
The final physics effect has to do with the modeling of $t\bar{t}$ production and decay. Different
|
105 |
|
|
MC models could in principle result in different BG predictions. Therefore we use several different
|
106 |
|
|
$t\bar{t}$ MC samples using different generators and dfferent parameters, to test the stability
|
107 |
|
|
of the dilepton BG prediction. All these predictions {\bf after} corrections for branching ratio
|
108 |
|
|
and $N_{jet}$ dependence, are compared to each other. The spread is a measure of the systematic
|
109 |
|
|
uncertainty associated with the $t\bar{t}$ generator modeling.
|
110 |
|
|
|
111 |
|
|
The main instrumental effect is associated with the underefficiency of the 2nd lepton veto.
|
112 |
|
|
We use tag-and-probe to compare the isolated track veto performance in $Z + 4$ jet data and
|
113 |
|
|
MC, and we extract corrections if necessary. Note that the performance of the isolated track veto
|
114 |
|
|
is not exactly the same on $e/\mu$ and on one prong hadronic tau decays. This is because
|
115 |
|
|
the pions from one-prong taus are often accompanied by $\pi^0$'s that can then result in extra
|
116 |
|
|
tracks due to phton conversions. We let the simulation take care of that. Similarly, at the moment
|
117 |
|
|
we also let the simulation take care of the possibility of a hadronic tau ``disappearing'' in the
|
118 |
|
|
detector due to nuclear interaction of the pion.
|
119 |
|
|
|
120 |
|
|
The sample of events failing the last isolated track veto is an important control sample to
|
121 |
|
|
check that we are doing the right thing.
|
122 |
|
|
|
123 |
|
|
Note that JES uncertainties are effectively ``calibrated away'' by the $N_{jet}$ rescaling described
|
124 |
|
|
above.
|
125 |
|
|
|
126 |
|
|
Finally, there are possible improvements to this basic analysis strategy that can be added in the future:
|
127 |
|
|
\begin{itemize}
|
128 |
|
|
\item Move from counting experiment to shape analysis. But first, we need to get the counting
|
129 |
|
|
experiment under control.
|
130 |
|
|
\item Add an explicit three prong tau veto
|
131 |
|
|
\item Do something to require that three of the jets in the event be consistent with $t \to Wb, W \to q\bar{q}$.
|
132 |
|
|
This could help rejecting some of the dilepton BG; however, it would also loose efficiency for
|
133 |
|
|
the $\widetilde{t} \to b \chi^+$ mode
|
134 |
|
|
\item Consider the $M(\ell b)$ variable, which is not bounded by $M_{top}$ in $\widetilde{t} \to b \chi^+$
|
135 |
|
|
\end{itemize} |