ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/claudioc/OSNote2010/limit.tex
Revision: 1.19
Committed: Fri Dec 3 15:17:37 2010 UTC (14 years, 5 months ago) by benhoob
Content type: application/x-tex
Branch: MAIN
Changes since 1.18: +3 -1 lines
Log Message:
Minor updates

File Contents

# User Rev Content
1 claudioc 1.1 \section{Limit on new physics}
2     \label{sec:limit}
3 claudioc 1.2
4 claudioc 1.10 %{\bf \color{red} The numbers in this Section need to be double checked.}
5 claudioc 1.2
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 benhoob 1.15 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 claudioc 1.2 The data driven background predictions from the ABCD method
19 benhoob 1.14 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 claudioc 1.2
22     These three predictions are in good agreement with each other
23     and with the observation of one event in the signal region.
24 benhoob 1.5 We calculate a Bayesian 95\% CL upper limit\cite{ref:bayes.f}
25 benhoob 1.17 on the number of non SM events in the signal region to be 4.1.
26 benhoob 1.18 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 claudioc 1.2 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 benhoob 1.15 LM0 and LM1 events from Table~\ref{tab:sigcont}: $6.3 \pm 1.3$
35 benhoob 1.16 events and $2.6 \pm 0.4$
36 benhoob 1.14 respectively, where the uncertainties
37 claudioc 1.2 are from energy scale (Section~\ref{sec:systematics}), luminosity,
38 claudioc 1.10 and lepton efficiency. Note that these expected SUSY yields
39     are computed using LO cross-sections, and are therefore underestimated.
40 claudioc 1.2
41 claudioc 1.10 Conveying additional useful information about the results of
42     a generic ``signature-based'' search such as the one described
43 claudioc 1.11 in this note is a difficult issue. The next paragraph represent
44 claudioc 1.10 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 benhoob 1.18 responses for SumJetPt and \met/$\sqrt{\rm SumJetPt}$. These
54 benhoob 1.19 quantities have been evaluated with Spring10 MC samples,
55     and we are currently checking if any of them change after
56     switching to Fall10 MC.
57 claudioc 1.10 The muon identification efficiency is $\approx 95\%$;
58     the electron identification efficiency varies from $\approx$ 63\% at
59     $P_T = 10$ GeV to 91\% for $P_T > 30$ GeV. The isolation
60     efficiency in top events varies from $\approx 83\%$ (muons)
61     and $\approx 89\%$ (electrons) at $P_T=10$ GeV to
62     $\approx 95\%$ for $P_T>60$ GeV. The average detector
63     responses for SumJetPt and $\met/\sqrt{\rm SumJetPt}$ are
64     $1.00 \pm 0.05$ and $0.94 \pm 0.05$ respectively, where
65     the uncertainties are from the jet energy scale uncertainty.
66     The experimental resolutions on these quantities are 10\% and
67     14\% respectively.
68    
69     To justify the statements in the previous paragraph
70     about the detector responses, we plot
71     in Figure~\ref{fig:response} the average response for
72 claudioc 1.8 SumJetPt and \met/$\sqrt{\rm SumJetPt}$ in MC, as well as the
73     efficiency for the cuts on these quantities used in defining the
74 claudioc 1.9 signal region.
75     % (SumJetPt $>$ 300 GeV and \met/$\sqrt{\rm SumJetPt} > 8.5$
76     % Gev$^{\frac{1}{2}}$).
77     We find that the average SumJetPt response
78 claudioc 1.8 in the Monte Carlo
79 claudioc 1.10 is very close to one, with an RMS of order 10\% while
80 claudioc 1.9 the
81 claudioc 1.8 response of \met/$\sqrt{\rm SumJetPt}$ is approximately 0.94 with an
82 claudioc 1.9 RMS of 14\%.
83 claudioc 1.8
84 claudioc 1.10 %Using this information as well as the kinematical
85     %cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies
86     %of Figures~\ref{fig:effttbar}, one should be able to confront
87     %any existing or future model via a relatively simple generator
88     %level study by comparing the expected number of events in 35 pb$^{-1}$
89     %with our upper limit of 4.1 events.
90 claudioc 1.8
91     \begin{figure}[tbh]
92     \begin{center}
93     \includegraphics[width=\linewidth]{selectionEff.png}
94     \caption{\label{fig:response} Left plots: the efficiencies
95     as a function of the true quantities for the SumJetPt (top) and
96     tcMET/$\sqrt{\rm SumJetPt}$ (bottom) requirements for the signal
97     region as a function of their true values. The value of the
98     cuts is indicated by the vertical line.
99     Right plots: The average response and its RMS for the SumJetPt
100     (top) and tcMET/$\sqrt{\rm SumJetPt}$ (bottom) measurements.
101     The response is defined as the ratio of the reconstructed quantity
102     to the true quantity in MC. These plots are done using the LM0
103     Monte Carlo, but they are not expected to depend strongly on
104     the underlying physics.}
105     \end{center}
106     \end{figure}