ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/benhoob/cmsnotes/StopSearch/overview.tex
Revision: 1.1
Committed: Wed Jun 27 10:50:04 2012 UTC (12 years, 10 months ago) by claudioc
Content type: application/x-tex
Branch: MAIN
Log Message:
beauty

File Contents

# User Rev Content
1 claudioc 1.1 \section{Overview and Analysis Strategy}
2     \label{sec:overview}
3    
4     We are searching for a $t\bar{t}\chi^0\chi^0$ or $W \ell b W \ell \bar{b} \chi^0 \chi^0$ final state
5     (after top decay in the first mode, the final states are actually the same). So to first order
6     this is ``$t\bar{t} +$ extra \met''.
7    
8     We work in the $\ell +$ jets final state, where the main background is $t\bar{t}$. We look for
9     \met inconsistent with $W \to \ell \nu$. We do this by concentrating on the $\ell \nu$ transverse
10     mass ($M_T$), since except for resolution effects, $M_T < M_W$ for $W \to \ell \nu$. Thus, the
11     initial analysis is simply a counting experiment in the tail of the $M_T$ distribution.
12    
13     The event selection is one-and-only-one high $P_T$ isolated lepton, four or more jets, and
14     some moderate \met cut. At least one of the jets has to be btagged to reduce $W+$ jets.
15     The event sample is then dominated by $t\bar{t}$, but there are also contributions from $W+$ jets,
16     single top, dibosons, etc.
17    
18     In order to be sensitive to $\widetilde{t}\widetilde{t}$ production, the background in the $M_T$
19     tail has to be controlled at the level of 10\% or better. So this is (almost) a precision measurement.
20    
21     The $t\bar{t}$ events in the $M_T$ tail can be broken up into two categories:
22     (i) $t\bar{t} \to \ell $+ jets and (ii) $t\bar{t} \to \ell^+ \ell^-$ where one of the two
23     leptons is not found by the second-lepton-veto (here the second lepton can be a hadronically
24     decaying $\tau$).
25     For a reasonable $M_T$ cut, say $M_T >$ 150 GeV, the dilepton background is of order 80\% of
26     the total. This is because in dileptons there are two neutrinos from $W$ decay, thus $M_T$
27     is not bounded by $M_W$. This is a very important point: while it is true that we are looking in
28     the tail of $M_T$, the bulk of the background events end up there not because of some exotic
29     \met reconstruction failure, but because of well understood physics processes. This means that
30     the background estimate can be taken from Monte Carlo (MC), after carefully accounting for possible
31     data/MC differences. Sophisticated fully ``data driven'' techniques are not really needed.
32    
33     Another important point is that in order to minimize systematic uncertainties, the MC background
34     predictions are always normalized to the bulk of the $t\bar{t}$ data, ie, events passing all of the
35     requirements but with $M_T \approx 80$ GeV.
36     This removes uncertainties
37     due to $\sigma(t\bar{t})$, lepton ID, trigger efficiency, luminosity, etc.
38    
39     The $\ell +$ jets background, which is dominated by
40     $t\bar{t} \to \ell $+ jets, but also includes some $W +$ jets as well as single top,
41     is estimated as follows:
42     \begin{enumerate}
43     \item We select a control sample of events passing all cuts, but anti-btagged. This is
44     sample is now dominated by $W +$ jets. The sample is used to understand the
45     $M_T$ tail in $\ell +$ jets processes.
46     \item In MC we measure the ratio of the number of $\ell +$ jets events in the $M_T$ tail to
47     the number of events with $M_T \approx$ 80 GeV. This ratio turns out to be pretty much the
48     same for all sources of $\ell +$ jets.
49     \item In data we measure the same ratio but after correcting for the $t\bar{t} \to$ dilepton
50     contribution, as well as dibosons etc. The dilepton contribution is taken from MC after
51     the correction described below.
52     \item We compare the two ratios, as well as the shapes of the data and MC $M_T$ distributions.
53     If they do not agree, we try to figure out why and fix it. If they agree well enough, we define a
54     data MC scale factor (SF) which is the ratio of the ratios defined in step 2 and 3, keeping track of the
55     uncertainty.
56     \item We next perform the full selection in $t\bar{t} \to \ell +$ jets MC, and measure this ratio
57     again (which should be the same as that in step 2).
58     \item We perform the full selection in data. We count the events with $M_T \approx 80$ GeV, we
59     subtract off the dilepton contribution, we multiply the subtracted event count by the ratio from step 5 (or from
60     step 2), and also by the data/MC SF from step 4. The result is the prediction for the $\ell +$ jets BG in
61     the $M_T$ tail.
62     \end{enumerate}
63    
64     The dilepton background can be broken up into many components depending
65     on the characteristics of the 2nd (undetected) lepton
66     \begin{itemize}
67     \item 3-prong hadronic tau decay
68     \item 1-prong hadronic tau decay
69     \item $e$ or $\mu$ possibly from $\tau$ decay
70     \end{itemize}
71     We have currently no veto against 3-prong taus. For the other two categories, we explicitely
72     veto events with additional electrons and muons above 10 GeV , and
73     we veto events with an isolated track of $P_T > 10$ GeV. This also rejects 1-prong taus
74     (it turns out that the explicit $e$ or $\mu$ veto is redundant with the isolated track veto).
75     Therefore the latter two categories can be broken into
76     \begin{itemize}
77     \item out of acceptance $(|\eta| > 2.50)$
78     \item $P_T < 10$ GeV
79     \item $P_T > 10$ GeV, but survives the additional lepton/track isolation veto
80     \end{itemize}
81     Monte Carlo studies indicate that there is no dominant contribution: it is ``a little bit of this,
82     and a little bit of that''.
83    
84     The high $M_T$ dilepton backgrounds come from MC, but their rate is normalized to the
85     $M_T \approx 80$ GeV peak. In other to perform this normalization in data, the $W +$ jets
86     events in the $M_T$ peak have to be subtracted off. This introduces a systematic uncertainty.
87    
88     There are two types of effects that can influence the MC dilepton prediction: physics effects
89     and instrumental effects. We discuss these next, starting from physics.
90    
91     First of all, many of our $t\bar{t}$ MC samples (eg: MadGraph) have
92     BR$(W \to \ell \nu)=\frac{1}{9} = 0.1111$.
93     PDG says BR$(W \to \ell \nu) = 0.1080 \pm 0.0009$. This difference matter, so the $t\bar{t}$ MC
94     must be corrected to account for this.
95    
96     Second, our selection is $\ell +$ 4 or more jets. A dilepton event passes the selection only if there are
97     two additional jet from ISR, or one jet from ISR and one jet which is reconstructed from the
98     unidentified lepton, {\it e.g.}, a three-prong tau. Therefore, all MC dilepton $t\bar{t}$ samples used
99     in the analysis must have their jet multiplicity corrected (if necessary) to agree with what is
100     seen in $t\bar{t}$ data. We use a data control sample of well identified dilepton events with
101     \met and at least two jets as a template to ``adjust'' the $N_{jet}$ distribution of the $t\bar{t} \to$
102     dileptons MC samples.
103    
104     The final physics effect has to do with the modeling of $t\bar{t}$ production and decay. Different
105     MC models could in principle result in different BG predictions. Therefore we use several different
106     $t\bar{t}$ MC samples using different generators and dfferent parameters, to test the stability
107     of the dilepton BG prediction. All these predictions {\bf after} corrections for branching ratio
108     and $N_{jet}$ dependence, are compared to each other. The spread is a measure of the systematic
109     uncertainty associated with the $t\bar{t}$ generator modeling.
110    
111     The main instrumental effect is associated with the underefficiency of the 2nd lepton veto.
112     We use tag-and-probe to compare the isolated track veto performance in $Z + 4$ jet data and
113     MC, and we extract corrections if necessary. Note that the performance of the isolated track veto
114     is not exactly the same on $e/\mu$ and on one prong hadronic tau decays. This is because
115     the pions from one-prong taus are often accompanied by $\pi^0$'s that can then result in extra
116     tracks due to phton conversions. We let the simulation take care of that. Similarly, at the moment
117     we also let the simulation take care of the possibility of a hadronic tau ``disappearing'' in the
118     detector due to nuclear interaction of the pion.
119    
120     The sample of events failing the last isolated track veto is an important control sample to
121     check that we are doing the right thing.
122    
123     Note that JES uncertainties are effectively ``calibrated away'' by the $N_{jet}$ rescaling described
124     above.
125    
126     Finally, there are possible improvements to this basic analysis strategy that can be added in the future:
127     \begin{itemize}
128     \item Move from counting experiment to shape analysis. But first, we need to get the counting
129     experiment under control.
130     \item Add an explicit three prong tau veto
131     \item Do something to require that three of the jets in the event be consistent with $t \to Wb, W \to q\bar{q}$.
132     This could help rejecting some of the dilepton BG; however, it would also loose efficiency for
133     the $\widetilde{t} \to b \chi^+$ mode
134     \item Consider the $M(\ell b)$ variable, which is not bounded by $M_{top}$ in $\widetilde{t} \to b \chi^+$
135     \end{itemize}