1 |
|
\section{Limit on new physics} |
2 |
|
\label{sec:limit} |
3 |
|
|
4 |
< |
{\bf \color{red} The numbers in this Section need to be double checked.} |
4 |
> |
%{\bf \color{red} The numbers in this Section need to be double checked.} |
5 |
|
|
6 |
|
As discussed in Section~\ref{sec:results}, we see one event |
7 |
|
in the signal region, defined as SumJetPt$>$300 GeV and |
8 |
|
\met/$\sqrt{\rm SumJetPt}>8.5$ GeV$^{\frac{1}{2}}$. |
9 |
|
|
10 |
|
The background prediction from the SM Monte Carlo is |
11 |
< |
1.4 $\pm$ 0.5 events, where the uncertainty comes from |
12 |
< |
the jet energy scale (30\%, see Section~\ref{sec:systematics}), |
13 |
< |
the luminosity (10\%), and the lepton/trigger |
14 |
< |
efficiency (10\%)\footnote{Other uncertainties associated with |
15 |
< |
the modeling of $t\bar{t}$ in MadGraph have not been evaluated. |
16 |
< |
The uncertainty on $pp \to \sigma(t\bar{t})$ is also not included.}. |
11 |
> |
1.3 events. |
12 |
> |
%, where the uncertainty comes from |
13 |
> |
%the jet energy scale (30\%, see Section~\ref{sec:systematics}), |
14 |
> |
%the luminosity (10\%), and the lepton/trigger |
15 |
> |
%efficiency (10\%)\footnote{Other uncertainties associated with |
16 |
> |
%the modeling of $t\bar{t}$ in MadGraph have not been evaluated. |
17 |
> |
%The uncertainty on $pp \to \sigma(t\bar{t})$ is also not included.}. |
18 |
|
The data driven background predictions from the ABCD method |
19 |
< |
and the $P_T(\ell\ell)$ method are 1.5 $\pm$ 0.9 and |
20 |
< |
$2.5 \pm 2.2$ events, respectively. |
19 |
> |
and the $P_T(\ell\ell)$ method are $1.5 \pm 0.9({\rm stat}) \pm 0.3({\rm syst})$ |
20 |
> |
and $4.3 \pm 3.0({\rm stat}) \pm 1.2({\rm syst})$, respectively. |
21 |
|
|
22 |
|
These three predictions are in good agreement with each other |
23 |
|
and with the observation of one event in the signal region. |
24 |
|
We calculate a Bayesian 95\% CL upper limit\cite{ref:bayes.f} |
25 |
|
on the number of non SM events in the signal region to be 4.1. |
26 |
< |
This was calculated using a background prediction of $N_{BG}=1.4 \pm 1.1$ |
26 |
> |
We have also calculated this limit using a profile likelihood method |
27 |
> |
as implemented in the cl95cms software, and we also find 4.1. |
28 |
> |
These limits were calculated using a background prediction of $N_{BG}=1.7 \pm 1.1$ |
29 |
|
events. The upper limit is not very sensitive to the choice of |
30 |
|
$N_{BG}$ and its uncertainty. |
31 |
|
|
32 |
|
To get a feeling for the sensitivity of this search to some |
33 |
|
popular SUSY models, we remind the reader of the number of expected |
34 |
< |
LM0 and LM1 events from Table~\ref{tab:sigcont}: $6.5 \pm 1.3$ |
35 |
< |
events and $2.6 \pm 0.4$ respectively, where the uncertainties |
34 |
> |
LM0 and LM1 events from Table~\ref{tab:sigcont}: $6.3 \pm 1.3$ |
35 |
> |
events and $2.6 \pm 0.4$ |
36 |
> |
respectively, where the uncertainties |
37 |
|
are from energy scale (Section~\ref{sec:systematics}), luminosity, |
38 |
< |
and lepton efficiency. |
38 |
> |
and lepton efficiency. Note that these expected SUSY yields |
39 |
> |
are computed using LO cross-sections, and are therefore underestimated. |
40 |
|
|
41 |
< |
In Figure~\ref{fig:response} we provide the response functions for the |
41 |
> |
Conveying additional useful information about the results of |
42 |
> |
a generic ``signature-based'' search such as the one described |
43 |
> |
in this note is a difficult issue. The next paragraph represent |
44 |
> |
our attempt at doing so. |
45 |
> |
|
46 |
> |
Other models of new physics in the dilepton final state |
47 |
> |
can be confronted in an approximate way by simple |
48 |
> |
generator-level studies that |
49 |
> |
compare the expected number of events in 35 pb$^{-1}$ |
50 |
> |
with our upper limit of 4.1 events. The key ingredients |
51 |
> |
of such studies are the kinematical cuts described |
52 |
> |
in this note, the lepton efficiencies, and the detector |
53 |
> |
responses for SumJetPt and \met/$\sqrt{\rm SumJetPt}$. |
54 |
> |
{LOOKING AT THE 38X MC PLOTS BY EYE, THE FOLLOWING QUANTITIES LOOK ABOUT RIGHT.} |
55 |
> |
The muon identification efficiency is $\approx 95\%$; |
56 |
> |
the electron identification efficiency varies from $\approx$ 63\% at |
57 |
> |
$P_T = 10$ GeV to 91\% for $P_T > 30$ GeV. The isolation |
58 |
> |
efficiency in top events varies from $\approx 83\%$ (muons) |
59 |
> |
and $\approx 89\%$ (electrons) at $P_T=10$ GeV to |
60 |
> |
$\approx 95\%$ for $P_T>60$ GeV. |
61 |
> |
{\bf \color{red} THE FOLLOWING QUANTITIES SHOULD BE RECALCULATED AFTER |
62 |
> |
WE FIX THE BUGS WITH THE MET IN LM SAMPLES} |
63 |
> |
The average detector |
64 |
> |
responses for SumJetPt and $\met/\sqrt{\rm SumJetPt}$ are |
65 |
> |
$1.00 \pm 0.05$ and $0.94 \pm 0.05$ respectively, where |
66 |
> |
the uncertainties are from the jet energy scale uncertainty. |
67 |
> |
The experimental resolutions on these quantities are 10\% and |
68 |
> |
14\% respectively. |
69 |
> |
|
70 |
> |
To justify the statements in the previous paragraph |
71 |
> |
about the detector responses, we plot |
72 |
> |
in Figure~\ref{fig:response} the average response for |
73 |
|
SumJetPt and \met/$\sqrt{\rm SumJetPt}$ in MC, as well as the |
74 |
|
efficiency for the cuts on these quantities used in defining the |
75 |
|
signal region. |
76 |
|
% (SumJetPt $>$ 300 GeV and \met/$\sqrt{\rm SumJetPt} > 8.5$ |
77 |
|
% Gev$^{\frac{1}{2}}$). |
78 |
+ |
{\bf \color{red} THE FOLLOWING QUANTITIES SHOULD BE RECALCULATED AFTER |
79 |
+ |
WE FIX THE BUGS WITH THE MET IN LM SAMPLES} |
80 |
|
We find that the average SumJetPt response |
81 |
|
in the Monte Carlo |
82 |
< |
is very close to one, with an RMS of order 10\%; |
82 |
> |
is very close to one, with an RMS of order 10\% while |
83 |
|
the |
84 |
|
response of \met/$\sqrt{\rm SumJetPt}$ is approximately 0.94 with an |
85 |
|
RMS of 14\%. |
86 |
|
|
87 |
< |
Using this information as well as the kinematical |
88 |
< |
cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies |
89 |
< |
of Figures~\ref{fig:effttbar}, one should be able to confront |
90 |
< |
any existing or future model via a relatively simple generator |
91 |
< |
level study by comparing the expected number of events in 35 pb$^{-1}$ |
92 |
< |
with our upper limit of 4.1 events. |
87 |
> |
%Using this information as well as the kinematical |
88 |
> |
%cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies |
89 |
> |
%of Figures~\ref{fig:effttbar}, one should be able to confront |
90 |
> |
%any existing or future model via a relatively simple generator |
91 |
> |
%level study by comparing the expected number of events in 35 pb$^{-1}$ |
92 |
> |
%with our upper limit of 4.1 events. |
93 |
|
|
94 |
|
\begin{figure}[tbh] |
95 |
|
\begin{center} |
104 |
|
The response is defined as the ratio of the reconstructed quantity |
105 |
|
to the true quantity in MC. These plots are done using the LM0 |
106 |
|
Monte Carlo, but they are not expected to depend strongly on |
107 |
< |
the underlying physics.} |
107 |
> |
the underlying physics. |
108 |
> |
{\bf \color{red} UPDATE AFTER FIXING BUGS WITH LM SAMPLES. } } |
109 |
|
\end{center} |
110 |
|
\end{figure} |