| 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.}. |
| 10 |
> |
The background prediction from the SM Monte Carlo is 1.3 events. |
| 11 |
> |
%, 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.}. |
| 17 |
|
The data driven background predictions from the ABCD method |
| 18 |
< |
and the $P_T(\ell\ell)$ method are 1.5 $\pm$ 0.9 and |
| 19 |
< |
$1.8^{+2.5}_{-1.8}$ events respectively. |
| 18 |
> |
and the $P_T(\ell\ell)$ method are $1.3 \pm 0.8({\rm stat}) \pm 0.3({\rm syst})$ |
| 19 |
> |
and $2.1 \pm 2.1({\rm stat}) \pm 0.6({\rm syst})$, respectively. |
| 20 |
|
|
| 21 |
|
These three predictions are in good agreement with each other |
| 22 |
|
and with the observation of one event in the signal region. |
| 23 |
< |
We calculate a baysean 95\% CL upper limit\cite{ref:bayes.f} |
| 23 |
> |
We calculate a Bayesian 95\% CL upper limit\cite{ref:bayes.f} |
| 24 |
|
on the number of non SM events in the signal region to be 4.1. |
| 25 |
< |
This was calculated using a background prediction of $N_{BG}=1.4 \pm 1.0$ |
| 26 |
< |
events. The upper limit is not very sensitive to the choice of |
| 25 |
> |
We have also calculated this limit using a profile likelihood method |
| 26 |
> |
as implemented in the cl95cms software, and we also find 4.1. |
| 27 |
> |
These limits were calculated using a background prediction of $N_{BG} = 1.4 \pm 0.8$ |
| 28 |
> |
events, the error-weighted average of the ABCD and $P_T(\ell\ell)$ background |
| 29 |
> |
predictions. 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:sigcontABCD}: $5.6 \pm 1.1$ |
| 35 |
< |
events and $2.2 \pm 0.3$ respectively, where the uncertainties |
| 34 |
> |
LM0 and LM1 events from Table~\ref{tab:sigcont}: $8.6 \pm 1.6$ |
| 35 |
> |
events and $3.6 \pm 0.5$ events respectively, where the uncertainties |
| 36 |
|
are from energy scale (Section~\ref{sec:systematics}), luminosity, |
| 37 |
|
and lepton efficiency. |
| 38 |
|
|
| 39 |
< |
In Figures XX and YY we provide the response functions for the |
| 40 |
< |
SumJetPt and \met/$\sqrt{\rm SumJetPt}$ cuts used in our analysis, |
| 41 |
< |
{\em i.e.} the efficiencies of the experimental cuts as a function of |
| 42 |
< |
the true quantities. Using this information as well as the kinematical |
| 43 |
< |
cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies |
| 44 |
< |
of Figures~\ref{fig:effttbar}, one should be able to confront |
| 45 |
< |
any existing or future model via a relatively simple generator |
| 46 |
< |
level study by comparing the expected number of events in 35 pb$^{-1}$ |
| 47 |
< |
with our upper limit of 4.1 events. |
| 39 |
> |
We also performed a scan of the mSUGRA parameter space. We set $\tan\beta=10$, |
| 40 |
> |
sign of $\mu = +$, $A_{0}=0$~GeV, and scan the $m_{0}$ and $m_{1/2}$ parameters |
| 41 |
> |
in steps of 10~GeV. For each scan point, we exclude the point if the expected |
| 42 |
> |
yield in the signal region exceeds 4.7, which is the 95\% CL upper limit |
| 43 |
> |
based on an expected background of $N_{BG}=1.4 \pm 0.8$ and a 20\% acceptance |
| 44 |
> |
uncertainty. The results are shown in Fig.~\ref{fig:msugra}. |
| 45 |
> |
|
| 46 |
> |
\begin{figure}[tbh] |
| 47 |
> |
\begin{center} |
| 48 |
> |
\includegraphics[width=0.6\linewidth]{msugra.png} |
| 49 |
> |
\caption{\label{fig:msugra}\protect Exclusion curve in the mSUGRA parameter space, |
| 50 |
> |
assuming $\tan\beta=10$, sign of $\mu = +$ and $A_{0}=0$~GeVs.} |
| 51 |
> |
\end{center} |
| 52 |
> |
\end{figure} |
| 53 |
> |
|
| 54 |
> |
|
| 55 |
> |
Conveying additional useful information about the results of |
| 56 |
> |
a generic ``signature-based'' search such as the one described |
| 57 |
> |
in this note is a difficult issue. The next paragraph represent |
| 58 |
> |
our attempt at doing so. |
| 59 |
> |
|
| 60 |
> |
Other models of new physics in the dilepton final state |
| 61 |
> |
can be confronted in an approximate way by simple |
| 62 |
> |
generator-level studies that |
| 63 |
> |
compare the expected number of events in 34.0~pb$^{-1}$ |
| 64 |
> |
with our upper limit of 4.1 events. The key ingredients |
| 65 |
> |
of such studies are the kinematical cuts described |
| 66 |
> |
in this note, the lepton efficiencies, and the detector |
| 67 |
> |
responses for SumJetPt and \met/$\sqrt{\rm SumJetPt}$. |
| 68 |
> |
{LOOKING AT THE 38X MC PLOTS BY EYE, THE FOLLOWING QUANTITIES LOOK ABOUT RIGHT.} |
| 69 |
> |
The muon identification efficiency is $\approx 95\%$; |
| 70 |
> |
the electron identification efficiency varies from $\approx$ 63\% at |
| 71 |
> |
$P_T = 10$ GeV to 91\% for $P_T > 30$ GeV. The isolation |
| 72 |
> |
efficiency in top events varies from $\approx 83\%$ (muons) |
| 73 |
> |
and $\approx 89\%$ (electrons) at $P_T=10$ GeV to |
| 74 |
> |
$\approx 95\%$ for $P_T>60$ GeV. {\bf \color{red} The following quantities were calculated |
| 75 |
> |
with Spring10 samples. } The average detector |
| 76 |
> |
responses for SumJetPt and $\met/\sqrt{\rm SumJetPt}$ are |
| 77 |
> |
$1.00 \pm 0.05$ and $0.94 \pm 0.05$ respectively, where |
| 78 |
> |
the uncertainties are from the jet energy scale uncertainty. |
| 79 |
> |
The experimental resolutions on these quantities are 10\% and |
| 80 |
> |
14\% respectively. |
| 81 |
> |
|
| 82 |
> |
To justify the statements in the previous paragraph |
| 83 |
> |
about the detector responses, we plot |
| 84 |
> |
in Figure~\ref{fig:response} the average response for |
| 85 |
> |
SumJetPt and \met/$\sqrt{\rm SumJetPt}$ in MC, as well as the |
| 86 |
> |
efficiency for the cuts on these quantities used in defining the |
| 87 |
> |
signal region. |
| 88 |
> |
% (SumJetPt $>$ 300 GeV and \met/$\sqrt{\rm SumJetPt} > 8.5$ |
| 89 |
> |
% Gev$^{\frac{1}{2}}$). |
| 90 |
> |
{\bf \color{red} The following numbers were derived from Spring10 samples.} |
| 91 |
> |
We find that the average SumJetPt response |
| 92 |
> |
in the Monte Carlo is very close to one, with an RMS of order 10\% while |
| 93 |
> |
the response of \met/$\sqrt{\rm SumJetPt}$ is approximately 0.94 with an |
| 94 |
> |
RMS of 14\%. |
| 95 |
> |
|
| 96 |
> |
%Using this information as well as the kinematical |
| 97 |
> |
%cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies |
| 98 |
> |
%of Figures~\ref{fig:effttbar}, one should be able to confront |
| 99 |
> |
%any existing or future model via a relatively simple generator |
| 100 |
> |
%level study by comparing the expected number of events in 35 pb$^{-1}$ |
| 101 |
> |
%with our upper limit of 4.1 events. |
| 102 |
> |
|
| 103 |
> |
\begin{figure}[tbh] |
| 104 |
> |
\begin{center} |
| 105 |
> |
\includegraphics[width=\linewidth]{selectionEffDec10.png} |
| 106 |
> |
\caption{\label{fig:response} Left plots: the efficiencies |
| 107 |
> |
as a function of the true quantities for the SumJetPt (top) and |
| 108 |
> |
tcMET/$\sqrt{\rm SumJetPt}$ (bottom) requirements for the signal |
| 109 |
> |
region as a function of their true values. The value of the |
| 110 |
> |
cuts is indicated by the vertical line. |
| 111 |
> |
Right plots: The average response and its RMS for the SumJetPt |
| 112 |
> |
(top) and tcMET/$\sqrt{\rm SumJetPt}$ (bottom) measurements. |
| 113 |
> |
The response is defined as the ratio of the reconstructed quantity |
| 114 |
> |
to the true quantity in MC. These plots are done using the LM0 |
| 115 |
> |
Monte Carlo, but they are not expected to depend strongly on |
| 116 |
> |
the underlying physics. |
| 117 |
> |
{\bf \color{red} These plots were made with Spring10 samples. } } |
| 118 |
> |
\end{center} |
| 119 |
> |
\end{figure} |
| 120 |
> |
|
| 121 |
> |
|
| 122 |
> |
|
| 123 |
> |
%%% Nominal |
| 124 |
> |
% ----------------------------------------- |
| 125 |
> |
% observed events 1 |
| 126 |
> |
% relative error on acceptance 0.000 |
| 127 |
> |
% expected background 1.400 |
| 128 |
> |
% absolute error on background 0.770 |
| 129 |
> |
% desired confidence level 0.95 |
| 130 |
> |
% integration upper limit 30.00 |
| 131 |
> |
% integration step size 0.0100 |
| 132 |
> |
% ----------------------------------------- |
| 133 |
> |
% Are the above correct? y |
| 134 |
> |
% 1 16.685 0.29375E-06 |
| 135 |
> |
% |
| 136 |
> |
% limit: less than 4.112 signal events |
| 137 |
> |
|
| 138 |
> |
|
| 139 |
> |
|
| 140 |
> |
%%% Add 20% acceptance uncertainty based on LM0 |
| 141 |
> |
% ----------------------------------------- |
| 142 |
> |
% observed events 1 |
| 143 |
> |
% relative error on acceptance 0.200 |
| 144 |
> |
% expected background 1.400 |
| 145 |
> |
% absolute error on background 0.770 |
| 146 |
> |
% desired confidence level 0.95 |
| 147 |
> |
% integration upper limit 30.00 |
| 148 |
> |
% integration step size 0.0100 |
| 149 |
> |
% ----------------------------------------- |
| 150 |
> |
% Are the above correct? y |
| 151 |
> |
% 1 29.995 0.50457E-06 |
| 152 |
> |
% |
| 153 |
> |
% limit: less than 4.689 signal events |
