ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/claudioc/OSNote2010/limit.tex
(Generate patch)

Comparing UserCode/claudioc/OSNote2010/limit.tex (file contents):
Revision 1.1 by claudioc, Fri Oct 29 02:29:40 2010 UTC vs.
Revision 1.31 by claudioc, Tue Dec 14 05:53:07 2010 UTC

# Line 1 | Line 1
1   \section{Limit on new physics}
2   \label{sec:limit}
3 < Nothing yet.
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.
45 > The results are shown in Fig.~\ref{fig:msugra}.
46 > This figure is still preliminary:
47 > \begin{itemize}
48 > \item The process dependent k-factors from Prospino were not yet available
49 > when the figure was made.  We took a flat k=1.4.
50 > \item The PDF uncertainties were still missing.
51 > \item The limits from other experiments are missing.  Wwe are hoping to
52 > converge on a common format for this plot with other SUSY analyses, so
53 > we have not made any attempt to make the plot look pretty (!).
54 > \item As mentioned above, we took a constant acceptance uncertainty
55 > instead of calculating the uncertainty point by point.
56 > \end{itemize}
57 >
58 > \begin{figure}[tbh]
59 > \begin{center}
60 > \includegraphics[width=0.6\linewidth]{msugra.png}
61 > \caption{\label{fig:msugra}\protect Exclusion curve in the mSUGRA parameter space,
62 > assuming $\tan\beta=10$, sign of $\mu = +$ and $A_{0}=0$~GeVs.}
63 > \end{center}
64 > \end{figure}
65 >
66 >
67 > Conveying additional useful information about the results of
68 > a generic ``signature-based'' search such as the one described
69 > in this note is a difficult issue.  The next paragraph represent
70 > our attempt at doing so.
71 >
72 > Other models of new physics in the dilepton final state
73 > can be confronted in an approximate way by simple
74 > generator-level studies that
75 > compare the expected number of events in 34.0~pb$^{-1}$
76 > with our upper limit of 4.1 events.  The key ingredients
77 > of such studies are the kinematical cuts described
78 > in this note, the lepton efficiencies, and the detector
79 > responses for SumJetPt and \met/$\sqrt{\rm SumJetPt}$.
80 > The muon identification efficiency is $\approx 95\%$;
81 > the electron identification efficiency varies from $\approx$ 63\% at
82 > $P_T = 10$ GeV to 91\% for $P_T > 30$ GeV.  The isolation
83 > efficiency in top events varies from $\approx 83\%$ (muons)
84 > and $\approx 89\%$ (electrons) at $P_T=10$ GeV to
85 > $\approx 95\%$ for $P_T>60$ GeV.
86 > %{\bf \color{red} The following numbers were derived from Fall 10 samples. }
87 > The average detector
88 > responses for SumJetPt and $\met/\sqrt{\rm SumJetPt}$ are
89 > $1.02 \pm 0.05$ and $0.94 \pm 0.05$ respectively, where
90 > the uncertainties are from the jet energy scale uncertainty.
91 > The experimental resolutions on these quantities are 11\% and
92 > 16\% respectively.
93 >
94 > To justify the statements in the previous paragraph
95 > about the detector responses, we plot
96 > in Figure~\ref{fig:response} the average response for
97 > SumJetPt and \met/$\sqrt{\rm SumJetPt}$ in MC, as well as the
98 > efficiency for the cuts on these quantities used in defining the
99 > signal region.
100 > % (SumJetPt $>$ 300 GeV and \met/$\sqrt{\rm SumJetPt} > 8.5$
101 > % Gev$^{\frac{1}{2}}$).  
102 > %{\bf \color{red} The following numbers were derived from Fall10 samples }
103 > We find that the average SumJetPt response
104 > in the Monte Carlo is about 1.02, with an RMS of order 11\% while
105 > the response of \met/$\sqrt{\rm SumJetPt}$ is approximately 0.94 with an
106 > RMS of 16\%.
107 >
108 > %Using this information as well as the kinematical
109 > %cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies
110 > %of Figures~\ref{fig:effttbar}, one should be able to confront
111 > %any existing or future model via a relatively simple generator
112 > %level study by comparing the expected number of events in 35 pb$^{-1}$
113 > %with our upper limit of 4.1 events.
114 >
115 > \begin{figure}[tbh]
116 > \begin{center}
117 > \includegraphics[width=\linewidth]{selectionEffDec10.png}
118 > \caption{\label{fig:response} Left plots: the efficiencies
119 > as a function of the true quantities for the SumJetPt (top) and
120 > tcMET/$\sqrt{\rm SumJetPt}$ (bottom) requirements for the signal
121 > region as a function of their true values.  The value of the
122 > cuts is indicated by the vertical line.
123 > Right plots: The average response and its RMS for the SumJetPt
124 > (top) and tcMET/$\sqrt{\rm SumJetPt}$ (bottom) measurements.
125 > The response is defined as the ratio of the reconstructed quantity
126 > to the true quantity in MC.  These plots are done using the LM0
127 > Monte Carlo, but they are not expected to depend strongly on
128 > the underlying physics.
129 > %{\bf \color{red} These plots were made with Fall10 samples. }
130 > }
131 > \end{center}
132 > \end{figure}
133 >
134 >
135 >
136 > %%%  Nominal
137 > % -----------------------------------------
138 > % observed events                         1
139 > % relative error on acceptance        0.000
140 > % expected background                 1.400
141 > % absolute error on background        0.770
142 > % desired confidence level             0.95
143 > % integration upper limit             30.00
144 > % integration step size              0.0100
145 > % -----------------------------------------
146 > % Are the above correct? y
147 > %    1  16.685     0.29375E-06
148 > %
149 > % limit: less than     4.112 signal events
150 >
151 >
152 >
153 > %%%  Add 20% acceptance uncertainty based on LM0
154 > % -----------------------------------------
155 > % observed events                         1
156 > % relative error on acceptance        0.200
157 > % expected background                 1.400
158 > % absolute error on background        0.770
159 > % desired confidence level             0.95
160 > % integration upper limit             30.00
161 > % integration step size              0.0100
162 > % -----------------------------------------
163 > % Are the above correct? y
164 > %    1  29.995     0.50457E-06
165 > %
166 > % limit: less than     4.689 signal events

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines