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.21 by benhoob, Mon Dec 6 16:37:30 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}$.
54 > {LOOKING AT THE 38X MC PLOTS BY EYE, THE FOLLOWING QUANTITIES LOOK ABOUT RIGHT.}
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.  
61 > {\bf \color{red} THE FOLLOWING QUANTITIES SHOULD BE RECALCULATED AFTER
62 > WE FIX THE BUGS WITH THE MET IN LM SAMPLES}
63 > The average detector
64 > responses for SumJetPt and $\met/\sqrt{\rm SumJetPt}$ are
65 > $1.00 \pm 0.05$ and $0.94 \pm 0.05$ respectively, where
66 > the uncertainties are from the jet energy scale uncertainty.
67 > The experimental resolutions on these quantities are 10\% and
68 > 14\% respectively.
69 >
70 > To justify the statements in the previous paragraph
71 > about the detector responses, we plot
72 > in Figure~\ref{fig:response} the average response for
73 > SumJetPt and \met/$\sqrt{\rm SumJetPt}$ in MC, as well as the
74 > efficiency for the cuts on these quantities used in defining the
75 > signal region.
76 > % (SumJetPt $>$ 300 GeV and \met/$\sqrt{\rm SumJetPt} > 8.5$
77 > % Gev$^{\frac{1}{2}}$).  
78 > {\bf \color{red} THE FOLLOWING QUANTITIES SHOULD BE RECALCULATED AFTER
79 > WE FIX THE BUGS WITH THE MET IN LM SAMPLES}
80 > We find that the average SumJetPt response
81 > in the Monte Carlo
82 > is very close to one, with an RMS of order 10\% while
83 > the
84 > response of \met/$\sqrt{\rm SumJetPt}$ is approximately 0.94 with an
85 > RMS of 14\%.
86 >
87 > %Using this information as well as the kinematical
88 > %cuts described in Section~\ref{sec:eventSel} and the lepton efficiencies
89 > %of Figures~\ref{fig:effttbar}, one should be able to confront
90 > %any existing or future model via a relatively simple generator
91 > %level study by comparing the expected number of events in 35 pb$^{-1}$
92 > %with our upper limit of 4.1 events.
93 >
94 > \begin{figure}[tbh]
95 > \begin{center}
96 > \includegraphics[width=\linewidth]{selectionEff.png}
97 > \caption{\label{fig:response} Left plots: the efficiencies
98 > as a function of the true quantities for the SumJetPt (top) and
99 > tcMET/$\sqrt{\rm SumJetPt}$ (bottom) requirements for the signal
100 > region as a function of their true values.  The value of the
101 > cuts is indicated by the vertical line.
102 > Right plots: The average response and its RMS for the SumJetPt
103 > (top) and tcMET/$\sqrt{\rm SumJetPt}$ (bottom) measurements.
104 > The response is defined as the ratio of the reconstructed quantity
105 > to the true quantity in MC.  These plots are done using the LM0
106 > Monte Carlo, but they are not expected to depend strongly on
107 > the underlying physics.
108 > {\bf \color{red} UPDATE AFTER FIXING BUGS WITH LM SAMPLES. } }
109 > \end{center}
110 > \end{figure}

Diff Legend

Removed lines
+ Added lines
< Changed lines
> Changed lines