ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/claudioc/OSNote2010/limit.tex
Revision: 1.26
Committed: Sat Dec 11 01:55:42 2010 UTC (14 years, 5 months ago) by dbarge
Content type: application/x-tex
Branch: MAIN
Changes since 1.25: +2 -2 lines
Log Message:
updated figure 7

File Contents

# Content
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.}
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 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.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 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 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: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 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, They should be the same for the Fall10 samples; this needs to be verified .}
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 Fall10 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