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.2 by claudioc, Thu Nov 11 12:34:17 2010 UTC vs.
Revision 1.18 by benhoob, Fri Dec 3 15:01:01 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.}.
11 > 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   The data driven background predictions from the ABCD method
19 < and the $P_T(\ell\ell)$ method are 1.5 $\pm$ 0.9 and
20 < $1.8^{+2.5}_{-1.8}$ events respectively.
19 > 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  
22   These three predictions are in good agreement with each other
23   and with the observation of one event in the signal region.
24 < We calculate a baysean 95\% CL upper limit\cite{ref:bayes.f}
24 > We calculate a Bayesian 95\% CL upper limit\cite{ref:bayes.f}
25   on the number of non SM events in the signal region to be 4.1.
26 < This was calculated using a background prediction of $N_{BG}=1.4 \pm 1.0$
26 > 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   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 < 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}: $6.3 \pm 1.3$
35 > events and $2.6 \pm 0.4$
36 > respectively, where the uncertainties
37   are from energy scale (Section~\ref{sec:systematics}), luminosity,
38 < and lepton efficiency.
38 > and lepton efficiency.  Note that these expected SUSY yields
39 > are computed using LO cross-sections, and are therefore underestimated.
40  
41 < In Figures XX and YY we provide the response functions for the
42 < SumJetPt and \met/$\sqrt{\rm SumJetPt}$ cuts used in our analysis,
43 < {\em i.e.} the efficiencies of the experimental cuts as a function of
44 < the true quantities.  Using this information as well as the kinematical
45 < cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies
46 < of Figures~\ref{fig:effttbar}, one should be able to confront
47 < any existing or future model via a relatively simple generator
48 < level study by comparing the expected number of events in 35 pb$^{-1}$
49 < with our upper limit of 4.1 events.
41 > Conveying additional useful information about the results of
42 > a generic ``signature-based'' search such as the one described
43 > in this note is a difficult issue.  The next paragraph represent
44 > 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 > responses for SumJetPt and \met/$\sqrt{\rm SumJetPt}$. These
54 > quantities have been evaluated with Spring10 MC samples.
55 > The muon identification efficiency is $\approx 95\%$;
56 > the electron identification efficiency varies from $\approx$ 63\% at
57 > $P_T = 10$ GeV to 91\% for $P_T > 30$ GeV.  The isolation
58 > efficiency in top events varies from $\approx 83\%$ (muons)
59 > and $\approx 89\%$ (electrons) at $P_T=10$ GeV to
60 > $\approx 95\%$ for $P_T>60$ GeV.  The average detector
61 > responses for SumJetPt and $\met/\sqrt{\rm SumJetPt}$ are
62 > $1.00 \pm 0.05$ and $0.94 \pm 0.05$ respectively, where
63 > the uncertainties are from the jet energy scale uncertainty.
64 > The experimental resolutions on these quantities are 10\% and
65 > 14\% respectively.
66 >
67 > To justify the statements in the previous paragraph
68 > about the detector responses, we plot
69 > in Figure~\ref{fig:response} the average response for
70 > SumJetPt and \met/$\sqrt{\rm SumJetPt}$ in MC, as well as the
71 > efficiency for the cuts on these quantities used in defining the
72 > signal region.
73 > % (SumJetPt $>$ 300 GeV and \met/$\sqrt{\rm SumJetPt} > 8.5$
74 > % Gev$^{\frac{1}{2}}$).  
75 > We find that the average SumJetPt response
76 > in the Monte Carlo
77 > is very close to one, with an RMS of order 10\% while
78 > the
79 > response of \met/$\sqrt{\rm SumJetPt}$ is approximately 0.94 with an
80 > RMS of 14\%.
81 >
82 > %Using this information as well as the kinematical
83 > %cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies
84 > %of Figures~\ref{fig:effttbar}, one should be able to confront
85 > %any existing or future model via a relatively simple generator
86 > %level study by comparing the expected number of events in 35 pb$^{-1}$
87 > %with our upper limit of 4.1 events.
88 >
89 > \begin{figure}[tbh]
90 > \begin{center}
91 > \includegraphics[width=\linewidth]{selectionEff.png}
92 > \caption{\label{fig:response} Left plots: the efficiencies
93 > as a function of the true quantities for the SumJetPt (top) and
94 > tcMET/$\sqrt{\rm SumJetPt}$ (bottom) requirements for the signal
95 > region as a function of their true values.  The value of the
96 > cuts is indicated by the vertical line.
97 > Right plots: The average response and its RMS for the SumJetPt
98 > (top) and tcMET/$\sqrt{\rm SumJetPt}$ (bottom) measurements.
99 > The response is defined as the ratio of the reconstructed quantity
100 > to the true quantity in MC.  These plots are done using the LM0
101 > Monte Carlo, but they are not expected to depend strongly on
102 > the underlying physics.}
103 > \end{center}
104 > \end{figure}

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines