57 |
|
|
58 |
|
Our choice of ABCD regions is shown in Figure~\ref{fig:abcdMC}. |
59 |
|
The signal region is region D. The expected number of events |
60 |
< |
in the four regions for the SM Monte Carlo, as well as the BG |
61 |
< |
prediction AC/B are given in Table~\ref{tab:abcdMC} for an integrated |
60 |
> |
in the four regions for the SM Monte Carlo, as well as the background |
61 |
> |
prediction A $\times$ C / B are given in Table~\ref{tab:abcdMC} for an integrated |
62 |
|
luminosity of 35 pb$^{-1}$. The ABCD method with chosen boundaries is accurate |
63 |
< |
to about 20\%. As shown in Table~\ref{tab:abcdsyst}, we assess systematic uncertainties |
64 |
< |
by varying the boundaries by an amount consistent with the hadronic energy scale uncertainty, |
65 |
< |
which we take as $\pm$5\% for SumJetPt and $\pm$2.5\% for MET/$\sqrt{\rm SumJetPt}$, since the |
66 |
< |
uncertainty on this quantity partially cancels due to the fact that it is a ratio of correlated |
67 |
< |
quantities. Based on these studies we assess a correction factor $k_{ABCD} = 1.2 \pm 0.2$ to the |
68 |
< |
predicted yield using the ABCD method. |
63 |
> |
to about 20\%, and we assess a corresponding systematic uncertainty. |
64 |
> |
|
65 |
> |
%As shown in Table~\ref{tab:abcdsyst}, we assess systematic uncertainties |
66 |
> |
%by varying the boundaries by an amount consistent with the hadronic energy scale uncertainty, |
67 |
> |
%which we take as $\pm$5\% for SumJetPt and $\pm$2.5\% for MET/$\sqrt{\rm SumJetPt}$, since the |
68 |
> |
%uncertainty on this quantity partially cancels due to the fact that it is a ratio of correlated |
69 |
> |
%quantities. Based on these studies we assess a correction factor $k_{ABCD} = 1.2 \pm 0.2$ to the |
70 |
> |
%predicted yield using the ABCD method. |
71 |
|
|
72 |
|
|
73 |
|
%{\color{red} Avi wants some statement about stability |
197 |
|
\begin{table}[htb] |
198 |
|
\begin{center} |
199 |
|
\caption{\label{tab:victorybad} |
200 |
< |
{\bf \color{red} Need to either update this with 38X MC, or replace it with the systematic studies varying the non-ttdil background yield and jet/met scale. } |
200 |
> |
{\bf \color{red} Need to either update this with 38X MC or remove it } |
201 |
|
Test of the data driven method in Monte Carlo |
202 |
|
under different assumptions. See text for details.} |
203 |
|
\begin{tabular}{|l|c|c|c|c|c|c|c|c|} |
204 |
|
\hline |
205 |
|
& True $t\bar{t}$ dilepton & $t\to W\to\tau$& other SM & GEN or & Lepton $P_T$ & Z veto & \met $>$ 50& obs/pred \\ |
206 |
< |
& included & included & included & RECOSIM & and $\eta$ cuts & & & \\ \hline |
206 |
> |
& included & included & included & RECOSIM & and $\eta$ cuts & & & \\ \hline |
207 |
|
1&Y & N & N & GEN & N & N & N & 1.90 \\ |
208 |
|
2&Y & N & N & GEN & Y & N & N & 1.64 \\ |
209 |
|
3&Y & N & N & GEN & Y & Y & N & 1.59 \\ |
211 |
|
5&Y & N & N & RECOSIM & Y & Y & Y & 1.51 \\ |
212 |
|
6&Y & Y & N & RECOSIM & Y & Y & Y & 1.58 \\ |
213 |
|
7&Y & Y & Y & RECOSIM & Y & Y & Y & 1.38 \\ |
212 |
– |
%%%NOTE: updated value 1.18 -> 1.46 since 2/3 DY events have been removed by updated analysis selections, |
213 |
– |
%%%dpt/pt cut and general lepton veto |
214 |
|
\hline |
215 |
|
\end{tabular} |
216 |
|
\end{center} |
217 |
|
\end{table} |
218 |
|
|
219 |
|
|
220 |
+ |
\begin{table}[htb] |
221 |
+ |
\begin{center} |
222 |
+ |
\caption{\label{tab:victorysyst} |
223 |
+ |
{Summary of uncertainties in $K_C$ due to the MET scale and resolution uncertainty, and to backgrounds other than $t\bar{t} \to$~dilepton. |
224 |
+ |
In the first table, `up' and `down' refer to shifting the hadronic energy scale up and down by 5\%. In the second table, the quoted value |
225 |
+ |
refers to the amount of additional smearing of the MET, as discussed in the text. In the third table, the normalization of all backgrounds |
226 |
+ |
other than $t\bar{t} \to$~dilepton is varied. |
227 |
+ |
{\bf \color{ref} Should I remove `observed' and `predicted' and show only the ratio? }} |
228 |
+ |
|
229 |
+ |
\begin{tabular}{ lcccc } |
230 |
+ |
\hline |
231 |
+ |
MET scale & Predicted & Observed & Obs/pred \\ |
232 |
+ |
\hline |
233 |
+ |
nominal & 0.92 $ \pm $ 0.11 & 1.29 $ \pm $ 0.11 & 1.40 $ \pm $ 0.20 \\ |
234 |
+ |
up & 0.92 $ \pm $ 0.11 & 1.53 $ \pm $ 0.12 & 1.66 $ \pm $ 0.23 \\ |
235 |
+ |
down & 0.81 $ \pm $ 0.07 & 1.08 $ \pm $ 0.11 & 1.32 $ \pm $ 0.17 \\ |
236 |
+ |
\hline |
237 |
+ |
|
238 |
+ |
\hline |
239 |
+ |
MET smearing & Predicted & Observed & Obs/pred \\ |
240 |
+ |
\hline |
241 |
+ |
nominal & 0.92 $ \pm $ 0.11 & 1.29 $ \pm $ 0.11 & 1.40 $ \pm $ 0.20 \\ |
242 |
+ |
10\% & 0.90 $ \pm $ 0.11 & 1.30 $ \pm $ 0.11 & 1.44 $ \pm $ 0.21 \\ |
243 |
+ |
20\% & 0.84 $ \pm $ 0.07 & 1.36 $ \pm $ 0.11 & 1.61 $ \pm $ 0.19 \\ |
244 |
+ |
30\% & 1.05 $ \pm $ 0.18 & 1.32 $ \pm $ 0.11 & 1.27 $ \pm $ 0.24 \\ |
245 |
+ |
40\% & 0.85 $ \pm $ 0.07 & 1.37 $ \pm $ 0.11 & 1.61 $ \pm $ 0.19 \\ |
246 |
+ |
50\% & 1.08 $ \pm $ 0.18 & 1.36 $ \pm $ 0.11 & 1.26 $ \pm $ 0.24 \\ |
247 |
+ |
\hline |
248 |
+ |
|
249 |
+ |
\hline |
250 |
+ |
non-$t\bar{t} \to$~dilepton scale factor & Predicted & Observed & Obs/pred \\ |
251 |
+ |
\hline |
252 |
+ |
ttdil only & 0.77 $ \pm $ 0.07 & 1.05 $ \pm $ 0.06 & 1.36 $ \pm $ 0.14 \\ |
253 |
+ |
nominal & 0.92 $ \pm $ 0.11 & 1.29 $ \pm $ 0.11 & 1.40 $ \pm $ 0.20 \\ |
254 |
+ |
double non-ttdil yield & 1.06 $ \pm $ 0.18 & 1.52 $ \pm $ 0.20 & 1.43 $ \pm $ 0.30 \\ |
255 |
+ |
\hline |
256 |
+ |
\end{tabular} |
257 |
+ |
\end{center} |
258 |
+ |
\end{table} |
259 |
+ |
|
260 |
+ |
|
261 |
+ |
|
262 |
|
The largest discrepancy between prediction and observation occurs on the first |
263 |
|
line of Table~\ref{tab:victorybad}, {\em i.e.}, at the generator level with no |
264 |
|
cuts. We have verified that this effect is due to the polarization of |
281 |
|
We have studied this effect at the generator level using Alpgen. We find |
282 |
|
that the bias is at the few percent level. |
283 |
|
|
242 |
– |
%%%TO BE REPLACED |
243 |
– |
%Based on the results of Table~\ref{tab:victorybad}, we conclude that the |
244 |
– |
%naive data driven background estimate based on $P_T{(\ell\ell)}$ needs to |
245 |
– |
%be corrected by a factor of {\color{red} $ K_{\rm{fudge}} =1.2 \pm 0.3$ |
246 |
– |
%(We still need to settle on thie exact value of this. |
247 |
– |
%For the 11 pb analysis it is taken as =1.)} . The quoted |
248 |
– |
%uncertainty is based on the stability of the Monte Carlo tests under |
249 |
– |
%variations of event selections, choices of \met algorithm, etc. |
250 |
– |
%For example, we find that observed/predicted changes by roughly 0.1 |
251 |
– |
%for each 1.5\% change in the average \met response. |
252 |
– |
|
284 |
|
Based on the results of Table~\ref{tab:victorybad}, we conclude that the |
285 |
< |
naive data driven background estimate based on $P_T{(\ell\ell)}$ needs to |
286 |
< |
be corrected by a factor of $ K_C = X \pm Y$. |
256 |
< |
The value of this correction factor as well as the systematic uncertainty |
257 |
< |
will be assessed using 38X ttbar madgraph MC. In the following we use |
258 |
< |
$K_C = 1$ for simplicity. Based on previous MC studies we foresee a correction |
259 |
< |
factor of $K_C \approx 1.2 - 1.5$, and we will assess an uncertainty |
260 |
< |
based on the stability of the Monte Carlo tests under |
261 |
< |
variations of event selections, choices of \met algorithm, etc. |
262 |
< |
For example, we find that observed/predicted changes by roughly 0.1 |
263 |
< |
for each 1.5\% change in the average \met response. |
285 |
> |
naive data-driven background estimate based on $P_T{(\ell\ell)}$ needs to |
286 |
> |
be corrected by a factor of $ K_C = 1.4 \pm 0.2(stat)$. |
287 |
|
|
288 |
+ |
The 2 dominant sources of systematic uncertainty in $K_C$ are due to non-$t\bar{t} \to$~dilepton backgrounds, |
289 |
+ |
and the MET scale and resolution uncertainties. The impact of non-$t\bar{t}$-dilepton background is assessed |
290 |
+ |
by varying the yield of all backgrounds except for $t\bar{t} \to$~dilepton, as shown in Table~\ref{table_kc}. |
291 |
+ |
The systematic is assessed as the larger of the differences between the nominal $K_C$ value and the values |
292 |
+ |
obtained using only $t\bar{t} \to$~dilepton MC and obtained by doubling the non $t\bar{t} \to$~dilepton component, |
293 |
+ |
giving an uncertainty of $0.04$. |
294 |
+ |
|
295 |
+ |
The uncertainty in $K_C$ due to the MET scale uncertainty is assessed by varying the hadronic energy scale using |
296 |
+ |
the same method as in~\ref{} and checking how much $K_C$ changes, as summarized in Table~\ref{tab:victorysyst}. |
297 |
+ |
This gives an uncertainty of 0.3. We also assess the impact of the MET resolution uncertainty on $K_C$ by applying |
298 |
+ |
a random smearing to the MET. For each event, we determine the expected MET resolution based on the sumJetPt, and |
299 |
+ |
smear the MET to simulate an increase in the resolution of 10\%, 20\%, 30\%, 40\% and 50\%. The results show that |
300 |
+ |
$K_C$ does not depend strongly on the MET resolution and we therefore do not assess any uncertainty. |
301 |
|
|
302 |
+ |
Incorporating all the statistical and systematic uncertainties we find $K_C = 1.4 \pm 0.4$. |
303 |
|
|
304 |
|
\subsection{Signal Contamination} |
305 |
|
\label{sec:sigcont} |