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.9 by claudioc, Sat Nov 13 17:13:42 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  
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.}.
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 < $2.5 \pm 2.2$  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 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.1$
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: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}: $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 Figure~\ref{fig:response} we provide the response functions for the
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
105 < is very close to one, with an RMS of order 10\%;
106 < the
107 < response of \met/$\sqrt{\rm SumJetPt}$ is approximately 0.94 with an
108 < RMS of 14\%.
109 <
110 < Using this information as well as the kinematical
111 < cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies
112 < of Figures~\ref{fig:effttbar}, one should be able to confront
113 < any existing or future model via a relatively simple generator
53 < level study by comparing the expected number of events in 35 pb$^{-1}$
54 < with our upper limit of 4.1 events.
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]{selectionEff.png}
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
# Line 66 | Line 125 | Right plots: The average response and it
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.}
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