ViewVC Help
View File | Revision Log | Show Annotations | Root Listing
root/cvsroot/UserCode/Vuko/Notes/WZCSA07/zjetbackground.tex
Revision: 1.11
Committed: Mon Jun 23 18:17:01 2008 UTC (16 years, 10 months ago) by vuko
Content type: application/x-tex
Branch: MAIN
CVS Tags: pre_release_23Jul08
Changes since 1.10: +30 -28 lines
Log Message:
commenting stuff that is not ready for first pre-release

File Contents

# User Rev Content
1 smorovic 1.1 \section{Signal extraction}
2 beaucero 1.6 \label{sec:SignalExt}
3 ymaravin 1.10 We separate backgrounds into two categories: one with a genuine \Z boson
4     from $\Z+jets$ processes, and the other without a genuine \Z boson from
5     $t\bar{t}$ and $\W+jet$ production. The latter source can be estimated from
6     the invariant mass of the \Z boson candidate, where the background events
7     with no genuine \Z boson should not produce a \Z mass peak and should
8     be relatively smooth.
9    
10     \subsection{Study of the background without a genuine \Z boson}
11     We estimate the background due to events without a genuine \Z boson from
12     fitting a \Z candidate invariant mass to a Gaussian function convoluted
13     with a Breit-Wigner function. The background is parameterized as a straight
14     line. An example of a fit for $3e$ category is given in Fig.~\ref{fig:ZFit}.
15     Both number of signal and background events are calculated for the
16     invariant mass range between 81 and 101 GeV. In Table~\ref{tab:FitVsMC}
17     we summarize the number of background events obtained from the fit
18     and from the Monte Carlo truth information.
19    
20 beaucero 1.8 \begin{figure}[!bp]
21     \begin{center}
22     \scalebox{0.4}{\includegraphics{figs/FitBkg3eTight.eps}}
23 ymaravin 1.10 \caption{The invariant mass distribution of the $Z$ boson candidate that is fit to a signal
24     parameterized as a Gaussian function convoluted with a Breit-Wigner function and
25     a background, parameterized as a straight line.}
26 beaucero 1.8 \label{fig:ZFit}
27     \end{center}
28     \end{figure}
29    
30 vuko 1.11 %%\begin{table}[!tb]
31     %\begin{center}
32     %\begin{tabular}{|l|c|c|c|c|c|c|c|} \hline
33     % & \multicolumn{2}{c|}{Background with genuine \Z} & \multicolumn{4}{c|}{Background without
34     % genuine \Z boson} \\
35     %Channel & $\Z+jets$ & $\Z b\bar{b}$ & $t\bar{t}$ & $\W+jets$ & $t\bar{t}$ + $\W+jets$ & Fit result \\ \hline
36     %$3e$ Loose & 196.5 & 67.4 & 35.7 & 0 & 35.7& 37.8 \\ \hline
37     %$3e$ Tight & 78.9 & 38.5 & 28.1 & 0 & 28.1 & 32.9 \\ \hline
38     %$2\mu 1e$ Loose & 189.6 & 52.6 & 4.7 & 0 & 4.7 & 5.7 \\ \hline
39     %$2\mu 1e$ Tight & 63.1 & 18.2 & 1.3 & 0 & 1.3 & 1.5 \\ \hline
40     %$2e1\mu$ & 10.4 & 9.1 & 30.7 & 2.9 & 33.6 & 29.2 \\ \hline
41     %$3\mu$ & 1.9 & 7.7 & 0.7 & 0 & 0.7 & 0\\ \hline
42     %\end{tabular}%
43 beaucero 1.8
44 vuko 1.11 %\end{center}
45     %\caption{Comparison between Monte Carlo truth information and the results of the fit for the background without genuine \Z boson. Number of events are obtained in the invariant mass range between 81 and 101 GeV.
46 ymaravin 1.10 %I AM NOT SURE I UNDERSTAND WHAT IS WRITTEN HERE
47     %The ``Loose'' and ``Tight'' selection criteria applied for $W\rightarrow e\nu$ final state only. One has to consider that this study as been perform on a smaller sample than the other part of the analysis a 10\% statistics error as to be counted until the study is performed on the whole samples.
48 beaucero 1.8 }
49 vuko 1.11 %\label{tab:FitVsMC}
50     %\end{table}
51 beaucero 1.8
52    
53 ymaravin 1.10 \subsection{Estimation of the background with genuine \Z boson}
54 beaucero 1.8 \label{sec:D0Matrix}
55 ymaravin 1.10 All of the instrumental background events with real \Z boson come from
56     $\Z + jets$ processes where one of the jets is misidentified as a lepton.
57     The probability to misidentify a jet as a muon is very low in CMS, while that
58     for the case of electron can be quite high as jets with large electromagnetic
59     energy fraction can be misidentified as electrons. $\Z+jet$ background
60     is especially high for \WZ\ signal with $\W\to e\nu$. Thus, it is imperative
61     to have a reliable estimation of this background from data to avoid
62     unnecessary systematic uncertainties due to Monte Carlo description of data
63     in startup conditions. Therefore, in the following we describe the data-driven
64     estimation of the $\Z+jets$ background for the $\ell^+\ell^- e$
65     categories. A similar study for the remaining $\ell^+\ell^- \mu$ categories
66     are in progress. However, as the $\Z+jets$ background is sufficiently small, it is
67     possible to use Monte Carlo simulation to estimate $\Z+jet$ background with
68     early data, without incurring significant systematic uncertainty due to data modeling.
69    
70     \subsubsection{$\Z+jets$ background fraction}
71     To estimate the fraction of the $\Z+jets$ events in data
72     we apply a method, commonly referred to as ``matrix'' method.
73     The idea of a method is to apply ``Loose'' identification criteria
74     on the third lepton after \Z boson candidate is identified
75     and count the number of the observed events, $N_{loose}$.
76     These events contain events with real electrons $N_{e}$
77     and events with misidentified jets $N_j$:
78     \begin{equation}
79     \label{eq:matrixEq1}
80     N_{loose} = N_e + N_j.
81     \end{equation}
82 beaucero 1.4
83 ymaravin 1.10 If we are to apply ``Tight'' selection on the third lepton, the number
84     of the observed events $N_{tight}$ would change as following
85     \begin{equation}
86     \label{eq:matrixEq2}
87     N_{tight} = \epsilon_{tight} N_e + p_{fake} N_j,
88     \end{equation}
89     where $\epsilon_{tight}$ and $p_{fake}$ are efficiency of ``Tight''
90     criteria with respect to ``Loose'' requirements for electrons and
91     misidentified jets, respectively. As $N_{loose}$ and $N_{tight}$
92     are directly observable, to extract the number of $Z+jet$ events
93     in the final sample, one needs to measure $\epsilon_{tight}$
94     and $p_{fake}$ in control data samples. Two possible ways
95     to estimate these values are given below.
96 smorovic 1.1
97 beaucero 1.8 \subsubsection{Determination of $\epsilon_{tight}$}
98 smorovic 1.9
99 vuko 1.11 % UNDEERSTAND THIS FIT: THIS CANNOT BE SHOWN AS IT IS!!!
100     %\begin{figure}[bt]
101     % \begin{center}
102     % \scalebox{0.8}{\includegraphics{figs/tag_probe_fit.eps}}
103     % \caption{Invariant mass of the \Z boson candidate for ``Tight-Tight'' (a)
104     % and ``Tight-Loose'' (b) electron selections fitted to a Gaussian with
105     % bifurcated Breit-Wigner functions.}
106     % \label{fig:tagprobe}
107     % \end{center}
108     %\end{figure}
109 smorovic 1.9
110 ymaravin 1.10 To estimate the $\epsilon_{tight}$ we apply ``tag-and-probe'' method
111     using $\Z \to e^+e^-$ from \Z+jets Chowder sample, including \W+jets
112 smorovic 1.9 and $t\bar{t}$ as background. \Z mass distribution is separated for two cases where
113 ymaravin 1.10 electrons from \Z decay either both pass ``Tigh'' selection (``Tight-Tight'' case), or only
114     one passes the ``Tight'' selection, while the other electron passes ``Loose'' but not ``Tight''
115     selection (``Tight-Loose'' case). To estimate signal in the selected \Z candidate invariant mass distribution, we fit it to a Gaussisan with bifurcated Breit-Wigner function as a signal
116     and straight line for a background model. \Z mass distribution and fit are shown in \ref{fig:tagprobe}.
117 smorovic 1.9
118     Equation for determination of signal efficiency is given as
119     \begin{equation}
120     \epsilon_{tight}=\frac{ 2*(N_{TT}-B_{TT}) }{ (N_{TL}-B_{TL})+2*(N_{TT}-B_{TT}) }
121     \end{equation}
122    
123 ymaravin 1.10 where $N_{TT}$,$B_{TT}$,$N_{TL}$ and $B_{TL}$ are, respectively, number of signal+background
124     and background events for ``Tight-Tight'' and ``Loose-Tight'' electron combinations.
125 vuko 1.11 %We estimated an efficiency $\epsilon_{tight}=0.99 \pm 0.01$.
126 smorovic 1.9
127 ymaravin 1.10 \subsubsection{Determination of $p_{fake}$}
128    
129     As the events will be most of the time triggered by the leptons coming
130     from \Z boson, we assume that the third lepton is unbiased toward the
131     trigger requirement. Ideally, we need a sample of pure multi-jet events
132     in order to compute the probability for a jet identified as a loose
133     electron to be also identified as a tight electron. In selecting such
134     a sample in data, one has to avoid any bias from the trigger
135     requirements on the ``Loose'' electron candidate.
136     %Such sample will not
137     %exist in data as they will be bias by the trigger requirement.
138     \begin{figure}[bt]
139     \begin{center}
140     \scalebox{0.6}{\includegraphics{figs/tight_eff_gumbo.eps}}
141     \caption{Fraction of electron candidates passing the ``Tight'' criteria
142 vuko 1.11 in multijet event. No trigger requirement has been applied.{\em NEW PLOT
143     WITH TRIGGER REQUIREMENTS TO COME}}
144 ymaravin 1.10 \label{fig:qcd_efftight_noHLT}
145     \end{center}
146     \end{figure}
147    
148     From a sample of multijet events triggered by an ``OR'' of multi-jet
149     triggers, we select a ``Loose'' electron candidate that are not
150     matched to any of the trigger objects. We also require the
151     electron candidate to be separated from the jet that satisfies
152     the trigger requirement by requiring the candidate to be separated
153     by at least $\Delta R = 0.2$ from the trigger object.
154     This allows us to obtain an unbiased sample of multijet events
155     where an electron candidate is likely to be either a converted
156     photon or a misidentified jet. The $p_{fake}$ function of $p_T$
157     and $\eta$ is simply obtained by dividing the $p_T$ and $\eta$
158     distributions for the electron candidate that satisfied ``Simple Tight''
159     electron identification requirements to that for electron candidates
160     that satisfied ``Simple Loose''. Such distributions are given
161     in Fig.~\ref{fig:qcd_zjet_est}.
162    
163