[ViewVC] Diff of: cvsroot/UserCode/benhoob/cmsnotes/ZMet2012/bkg.tex

Comparing UserCode/benhoob/cmsnotes/ZMet2012/bkg.tex (file contents):
Revision 1.4 by benhoob, Fri Jun 29 19:41:46 2012 UTC vs.
Revision 1.20 by benhoob, Fri Jan 25 17:46:18 2013 UTC

+\Ht, and Z \pt in the \zjets\ event. The prediction for the \Z sample is simply the sum of all such templates.
+All templates are displayed in App.~\ref{app:templates}.
-<
+While there is in principle a small contribution from backgrounds other than \zjets\ in the preselection regions,
-<
+this contribution is only $\approx$3\% ($\approx$2\%) of the total sample in the inclusive search (targeted search),
-<
+as shown in Table~\ref{table:zyields_2j} (Table~\ref{table:zyields_2j_targeted}), and is therefore negligible compared to the total
-<
+background uncertainty.
->
+After preselection, there is  a small contribution from backgrounds other than \zjets. To correct for this, the \MET\ templates
->
+prediction is scaled such that the total background prediction matches the observed data yield in the \MET\ 0--60 GeV region.
->
+Because the non-\zjets impurity in the low \MET\ region after preselection is very small, this results in
->
+scaling factors of 0.985 (0.995) for the inclusive (targeted) search.
+\subsection{Estimating the Flavor-Symmetric Background with e$\mu$ Events}
+\label{sec:bkg_fs}
+\begin{equation}
+R_{\mu e} = \sqrt{\frac{N_{\mu\mu}^{\rm{offline}}}{N_{ee}^{\rm{offline}}}} = \sqrt{\frac{N_{\mu\mu}^{\rm{trig}}/\epsilon_{\mu\mu}^{\rm{trig}}}{N_{ee}^{\rm{trig}}/\epsilon_{ee}^{\rm{trig}}}}
-<
+= \sqrt{\frac{80367/0.88}{54426/0.95}} = 1.26\pm0.07.
->
+= \sqrt{\frac{234132/0.86}{185555/0.95}} = 1.18\pm0.07.
+\end{equation}
+Here we have used the Z$\to\mu\mu$ and Z$\to$ee yields from Table~\ref{table:zyields_2j} and the trigger efficiencies quoted in Sec.~\ref{sec:datasets}.
+\begin{itemize}
+\item $N_{ee}^{\rm{predicted}}    = \frac {N_{e\mu}^{\rm{trig}}} {\epsilon_{e\mu}^{\rm{trig}}} \frac {\epsilon_{ee}^{\rm{trig}}} {2 R_{\mu e}}
-<
+= \frac{N_{e\mu}^{\rm{trig}}}{0.92}\frac{0.95}{2\times1.26} = (0.41\pm0.05) \times N_{e\mu}^{\rm{trig}}$ ,
->
+= \frac{N_{e\mu}^{\rm{trig}}}{0.93}\frac{0.95}{2\times1.18} = (0.43\pm0.05) \times N_{e\mu}^{\rm{trig}}$ ,
+\item $N_{\mu\mu}^{\rm{predicted}} = \frac {N_{e\mu}^{\rm{trig}}} {\epsilon_{e\mu}^{\rm{trig}}} \frac {\epsilon_{\mu\mu}^{\rm{trig}} R_{\mu e}}  {2}
-<
+= \frac {N_{e\mu}^{\rm{trig}}} {0.95} \frac {0.88 \times 1.26}{2} = (0.58\pm0.07) \times N_{e\mu}^{\rm{trig}}$,
->
+= \frac {N_{e\mu}^{\rm{trig}}} {0.95} \frac {0.86 \times 1.18}{2} = (0.53\pm0.07) \times N_{e\mu}^{\rm{trig}}$,
+\end{itemize}
+and the predicted yield in the combined ee and $\mu\mu$ channel is simply the sum of these two predictions:
+\begin{itemize}
-<
+\item $N_{ee+\mu\mu}^{\rm{predicted}} = (0.99\pm0.06)\times N_{e\mu}^{\rm{trig}}$.
->
+\item $N_{ee+\mu\mu}^{\rm{predicted}} = (0.97\pm0.06)\times N_{e\mu}^{\rm{trig}}$.
+\end{itemize}
+Note that the relative uncertainty in the combined ee and $\mu\mu$ prediction is smaller than those for the individual ee and $\mu\mu$ predictions
+To improve the statistical precision of the FS background estimate, we remove the requirement that the e$\mu$ lepton pair falls in the Z mass window.
+Instead we scale the e$\mu$ yield by $K$, the efficiency for e$\mu$ events to satisfy the Z mass requirement, extracted from simulation. In Fig.~\ref{fig:K_incl}
-<
+we display the value of $K$ in data and simulation, for a variety of \MET\ requirements, for the inclusive analysis. Based on this we chose $K=0.14\pm0.02$
-<
+for all \MET\ regions except for \MET\ $>$ 300 GeV. For this region the statistical precision is reduced, so that we inflate the uncertainty and chose $K=0.14\pm0.08$.
->
+we display the value of $K$ in data and simulation, for a variety of \MET\ requirements, for the inclusive analysis.
->
+Based on this we chose $K=0.14\pm0.02$ for the lower \MET\ regions, $K=0.14\pm0.04$ for the \MET\ $>$ 200 GeV region, and $K=0.14\pm0.08$ for \MET\ $>$ 300 GeV,
->
+where the larger uncertainties reflect the reduced statistical precision at large \MET.
+The corresponding plot for the targeted analysis, including the b-veto, is displayed in Fig.~\ref{fig:K_targeted}.
-<
+Based on this we chose $K=0.13\pm0.02$
-<
+for all \MET\ regions up to  \MET\ $>$ 100 GeV. For higher \MET\ regions (\MET\ $>$ 150 GeV and above) the statistical precision is reduced,
-<
+so that we inflate the uncertainty and chose $K=0.13\pm0.07$.
->
+Based on this we chose $K=0.13\pm0.02$ for all \MET\ regions up to  \MET\ $>$ 100 GeV.
->
+For the \MET\ $>$ 150 GeV region we choose $K=0.13\pm0.03$
->
+and for the \MET\ $>$ 200 GeV region we choose $K=0.13\pm0.05$,
->
+due to the reduced  statistical precision.
+\begin{figure}[!ht]
+\begin{center}
+\begin{tabular}{cc}
-<
+\includegraphics[width=0.4\textwidth]{plots/K_incl.pdf} &
-<
+\includegraphics[width=0.4\textwidth]{plots/K_excl.pdf} \\
->
+\includegraphics[width=0.4\textwidth]{plots/extractK_inclusive_19p5fb.pdf} &
->
+\includegraphics[width=0.4\textwidth]{plots/extractK_exclusive_19p5fb.pdf} \\
+\end{tabular}
-<
+\caption{
->
+\caption{\label{fig:K_incl}
+The efficiency for e$\mu$ events to satisfy the dilepton mass requirement, $K$, in data and simulation for inclusive \MET\ intervals (left) and
-<
+exclusive \MET\ intervals (right) for the inclusive analysis. Based on this we chose $K=0.14\pm0.02$ for all \MET\ regions except \MET\ $>$ 300 GeV,
-<
+where we chose $K=0.14\pm0.08$.
-<
+%{\bf FIXME plots made with 10\% of \zjets\ MC statistics, to be remade with full statistics}
-<
+\label{fig:K_incl}
->
+exclusive \MET\ intervals (right) for the inclusive analysis.
->
+Based on this we chose $K=0.14\pm0.02$ for the lower \MET\ regions, $K=0.14\pm0.04$ for the \MET\ $>$ 200 GeV region, and $K=0.14\pm0.08$ for \MET\ $>$ 300 GeV.
-+
-+
+\begin{comment}
-+
-+
+Using selection : ((((leptype==2)&&(csc==0 && hbhe==1 && hcallaser==1 && ecaltp==1 && trkfail==1 && eebadsc==1 && hbhenew==1))&&(isdata==0 || (run<197556 || run>198913)))&&(njets>=2))&&(lep1.pt()>20 && lep2.pt()>20)
-+
+Using weight    : vtxweight * weight
-+
+OF entries (total)  43808
-+
+OF entries (Z mass) 6021
-+
+K                   0.137441
-+
+Info in <TCanvas::MakeDefCanvas>:  created default TCanvas with name c1
-+
-+
+--------------------------------------------------------------
-+
+pfmet>0
-+
-+
+data  :
-+
+total : 43808
-+
+Z     : 6021
-+
+K     : 0.14 +/- 0.002
-+
-+
+MC    :
-+
+total : 2378.42
-+
+Z     : 344.559
-+
+K     : 0.14 +/- 0.002
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30
-+
-+
+data  :
-+
+total : 36603
-+
+Z     : 5084
-+
+K     : 0.14 +/- 0.002
-+
-+
+MC    :
-+
+total : 2012.6
-+
+Z     : 297.342
-+
+K     : 0.15 +/- 0.002
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60
-+
-+
+data  :
-+
+total : 22692
-+
+Z     : 3209
-+
+K     : 0.14 +/- 0.002
-+
-+
+MC    :
-+
+total : 1285.07
-+
+Z     : 189.292
-+
+K     : 0.15 +/- 0.002
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100
-+
-+
+data  :
-+
+total : 7862
-+
+Z     : 1093
-+
+K     : 0.14 +/- 0.004
-+
-+
+MC    :
-+
+total : 470.932
-+
+Z     : 68.9364
-+
+K     : 0.15 +/- 0.003
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200
-+
-+
+data  :
-+
+total : 424
-+
+Z     : 50
-+
+K     : 0.12 +/- 0.017
-+
-+
+MC    :
-+
+total : 28.2757
-+
+Z     : 2.87288
-+
+K     : 0.10 +/- 0.011
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>300
-+
-+
+data  :
-+
+total : 52
-+
+Z     : 5
-+
+K     : 0.10 +/- 0.043
-+
-+
+MC    :
-+
+total : 3.77378
-+
+Z     : 0.235632
-+
+K     : 0.06 +/- 0.023
-+
+--------------------------------------------------------------
-+
-+
+Using selection : ((((leptype==2)&&(csc==0 && hbhe==1 && hcallaser==1 && ecaltp==1 && trkfail==1 && eebadsc==1 && hbhenew==1))&&(isdata==0 || (run<197556 || run>198913)))&&(njets>=2))&&(lep1.pt()>20 && lep2.pt()>20)
-+
+Using weight    : vtxweight * weight
-+
+OF entries (total)  43808
-+
+OF entries (Z mass) 6021
-+
+K                   0.137441
-+
+Warning in <TROOT::Append>: Replacing existing TH1: htot (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: hZ (Potential memory leak).
-+
-+
+--------------------------------------------------------------
-+
+pfmet>0   && pfmet<30
-+
-+
+data  :
-+
+total : 7205
-+
+Z     : 937
-+
+K     : 0.13 +/- 0.004
-+
-+
+MC    :
-+
+total : 366.332
-+
+Z     : 47.2379
-+
+K     : 0.13 +/- 0.004
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30  && pfmet<60
-+
-+
+data  :
-+
+total : 13911
-+
+Z     : 1875
-+
+K     : 0.13 +/- 0.003
-+
-+
+MC    :
-+
+total : 727.951
-+
+Z     : 108.068
-+
+K     : 0.15 +/- 0.003
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60  && pfmet<100
-+
-+
+data  :
-+
+total : 14830
-+
+Z     : 2116
-+
+K     : 0.14 +/- 0.003
-+
-+
+MC    :
-+
+total : 814.344
-+
+Z     : 120.355
-+
+K     : 0.15 +/- 0.003
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100 && pfmet<200
-+
-+
+data  :
-+
+total : 7438
-+
+Z     : 1043
-+
+K     : 0.14 +/- 0.004
-+
-+
+MC    :
-+
+total : 442.657
-+
+Z     : 66.0631
-+
+K     : 0.15 +/- 0.004
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200 && pfmet<300
-+
-+
+data  :
-+
+total : 372
-+
+Z     : 45
-+
+K     : 0.12 +/- 0.018
-+
-+
+MC    :
-+
+total : 24.502
-+
+Z     : 2.63725
-+
+K     : 0.11 +/- 0.012
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>300
-+
-+
+data  :
-+
+total : 52
-+
+Z     : 5
-+
+K     : 0.10 +/- 0.043
-+
-+
+MC    :
-+
+total : 3.77378
-+
+Z     : 0.235632
-+
+K     : 0.06 +/- 0.023
-+
+--------------------------------------------------------------
-+
-+
-+
+\end{comment}
-+
+\end{center}
+\end{figure}
+\begin{figure}[!hb]
+\begin{center}
+\begin{tabular}{cc}
-<
+\includegraphics[width=0.4\textwidth]{plots/extractK_inclusive_bveto.pdf} &
-<
+\includegraphics[width=0.4\textwidth]{plots/extractK_exclusive_bveto.pdf} \\
->
+\includegraphics[width=0.4\textwidth]{plots/extractK_inclusive_bveto_19p5fb.pdf} &
->
+\includegraphics[width=0.4\textwidth]{plots/extractK_exclusive_bveto_19p5fb.pdf} \\
+\end{tabular}
+\caption{
+The efficiency for e$\mu$ events to satisfy the dilepton mass requirement, $K$, in data and simulation for inclusive \MET\ intervals (left) and
+exclusive \MET\ intervals (right) for the targeted analysis, including the b-veto.
+Based on this we chose $K=0.13\pm0.02$ for the \MET\ regions up to \MET\ $>$ 100 GeV.
-<
+For higher \MET\ regions we chose $K=0.13\pm0.07$.
->
+For \MET\ $>$ 150 we choose $K=0.13\pm0.03$, for \MET\ $>$ 200 GeV we choose $K=0.13\pm0.05$.
+%{\bf FIXME plots made with 10\% of \zjets\ MC statistics, to be remade with full statistics}
+\label{fig:K_targeted}
-+
+\begin{comment}
-+
-+
+Using selection : (((((leptype==2)&&(csc==0 && hbhe==1 && hcallaser==1 && ecaltp==1 && trkfail==1 && eebadsc==1 && hbhenew==1))&&(isdata==0 || (run<197556 || run>198913)))&&(njets>=2))&&(lep1.pt()>20 && lep2.pt()>20))&&(nbcsvm==0)
-+
+Using weight    : vtxweight * weight
-+
+OF entries (total)  11172
-+
+OF entries (Z mass) 1331
-+
+K                   0.119137
-+
+Info in <TCanvas::MakeDefCanvas>:  created default TCanvas with name c1
-+
-+
+--------------------------------------------------------------
-+
+pfmet>0
-+
-+
+data  :
-+
+total : 11172
-+
+Z     : 1331
-+
+K     : 0.12 +/- 0.003
-+
-+
+MC    :
-+
+total : 556.3
-+
+Z     : 72.3357
-+
+K     : 0.13 +/- 0.003
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30
-+
-+
+data  :
-+
+total : 8811
-+
+Z     : 1085
-+
+K     : 0.12 +/- 0.004
-+
-+
+MC    :
-+
+total : 447.641
-+
+Z     : 60.0542
-+
+K     : 0.13 +/- 0.003
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60
-+
-+
+data  :
-+
+total : 5263
-+
+Z     : 677
-+
+K     : 0.13 +/- 0.005
-+
-+
+MC    :
-+
+total : 285.463
-+
+Z     : 39.2608
-+
+K     : 0.14 +/- 0.004
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>80
-+
-+
+data  :
-+
+total : 3325
-+
+Z     : 422
-+
+K     : 0.13 +/- 0.006
-+
-+
+MC    :
-+
+total : 183.689
-+
+Z     : 25.7671
-+
+K     : 0.14 +/- 0.005
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100
-+
-+
+data  :
-+
+total : 1883
-+
+Z     : 234
-+
+K     : 0.12 +/- 0.008
-+
-+
+MC    :
-+
+total : 111.774
-+
+Z     : 14.7812
-+
+K     : 0.13 +/- 0.006
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>150
-+
-+
+data  :
-+
+total : 451
-+
+Z     : 46
-+
+K     : 0.10 +/- 0.015
-+
-+
+MC    :
-+
+total : 29.4551
-+
+Z     : 3.57377
-+
+K     : 0.12 +/- 0.012
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200
-+
-+
+data  :
-+
+total : 138
-+
+Z     : 15
-+
+K     : 0.11 +/- 0.028
-+
-+
+MC    :
-+
+total : 8.60692
-+
+Z     : 0.775129
-+
+K     : 0.09 +/- 0.017
-+
+--------------------------------------------------------------
-+
-+
+Using selection : (((((leptype==2)&&(csc==0 && hbhe==1 && hcallaser==1 && ecaltp==1 && trkfail==1 && eebadsc==1 && hbhenew==1))&&(isdata==0 || (run<197556 || run>198913)))&&(njets>=2))&&(lep1.pt()>20 && lep2.pt()>20))&&(nbcsvm==0)
-+
+Using weight    : vtxweight * weight
-+
+OF entries (total)  11172
-+
+OF entries (Z mass) 1331
-+
+K                   0.119137
-+
+Warning in <TROOT::Append>: Replacing existing TH1: htot (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: hZ (Potential memory leak).
-+
-+
+--------------------------------------------------------------
-+
+pfmet>0   && pfmet<30
-+
-+
+data  :
-+
+total : 2361
-+
+Z     : 246
-+
+K     : 0.10 +/- 0.007
-+
-+
+MC    :
-+
+total : 108.378
-+
+Z     : 12.2795
-+
+K     : 0.11 +/- 0.008
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30  && pfmet<60
-+
-+
+data  :
-+
+total : 3548
-+
+Z     : 408
-+
+K     : 0.11 +/- 0.006
-+
-+
+MC    :
-+
+total : 162.246
-+
+Z     : 20.7935
-+
+K     : 0.13 +/- 0.006
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60  && pfmet<80
-+
-+
+data  :
-+
+total : 1938
-+
+Z     : 255
-+
+K     : 0.13 +/- 0.008
-+
-+
+MC    :
-+
+total : 101.801
-+
+Z     : 13.4941
-+
+K     : 0.13 +/- 0.007
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>80  && pfmet<100
-+
-+
+data  :
-+
+total : 1442
-+
+Z     : 188
-+
+K     : 0.13 +/- 0.010
-+
-+
+MC    :
-+
+total : 71.9073
-+
+Z     : 10.986
-+
+K     : 0.15 +/- 0.009
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100 && pfmet<150
-+
-+
+data  :
-+
+total : 1432
-+
+Z     : 188
-+
+K     : 0.13 +/- 0.010
-+
-+
+MC    :
-+
+total : 82.3186
-+
+Z     : 11.2075
-+
+K     : 0.14 +/- 0.008
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>150 && pfmet<200
-+
-+
+data  :
-+
+total : 313
-+
+Z     : 31
-+
+K     : 0.10 +/- 0.018
-+
-+
+MC    :
-+
+total : 20.8482
-+
+Z     : 2.79864
-+
+K     : 0.13 +/- 0.015
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200
-+
-+
+data  :
-+
+total : 138
-+
+Z     : 15
-+
+K     : 0.11 +/- 0.028
-+
-+
+MC    :
-+
+total : 8.60692
-+
+Z     : 0.775129
-+
+K     : 0.09 +/- 0.017
-+
+--------------------------------------------------------------
-+
-+
-+
-+
-+
+\end{comment}
-+
+\end{center}
+\end{figure}
-+
+\clearpage
+\subsection{Estimating the WZ and ZZ Background with MC}
+Backgrounds from W($\ell\nu$)Z($\ell\ell$) where the W lepton is not identified or is outside acceptance, and Z($\nu\nu$)Z($\ell\ell$),
+are estimated from simulation. The MC modeling of these processes is validated by comparing the MC predictions with data in control samples
+with exactly 3 leptons (WZ control sample) and exactly 4 leptons (ZZ control sample).
-<
+The relevant WZ and ZZ MC samples are:
-<
-<
+\begin{itemize}
-<
+\footnotesize{
-<
+\item \verb=/WZJetsTo3LNu_TuneZ2_8TeV-madgraph-tauola/Summer12-PU_S7_START52_V9-v2/AODSIM= ($\sigma=1.058$ pb),
-<
+\item \verb=/ZZJetsTo4L_TuneZ2star_8TeV-madgraph-tauola/Summer12-PU_S7_START52_V9-v3/AODSIM= ($\sigma=0.093$ pb),
-<
+}
-<
+\end{itemize}
-<
+The WZJetsTo2L2Q, ZZJetsTo2L2Q, and ZZJetsTo2L2Nu samples are also used in this analysis but their contribution to the 3-lepton and 4-lepton
-<
+control samples is negligible.
->
+The critical samples are the WZJetsTo3LNu and ZZJetsTo4L, listed in Table~\ref{tab:mcsamples}
->
+(the WZJetsTo2L2Q, ZZJetsTo2L2Q, and ZZJetsTo2L2Nu samples are also used in this analysis but their contribution to the 3-lepton and 4-lepton
->
+control samples is negligible).
+\subsubsection{WZ Validation Studies}
+\label{sec:bkg_wz}
+\begin{itemize}
+\item Exactly 3 $p_T>20$~GeV leptons passing analysis identication and isolation requirements,
+\item 2 of the 3 leptons must fall in the Z window 81-101 GeV,
-<
+\item \MET $>$ 50 GeV (to suppress DY).
->
+\item \MET $>$ 50 GeV (to suppress DY),
->
+\item veto events with b-tagged jets.
+\end{itemize}
+The data and MC yields passing the above selection are in Table~\ref{tab:wz}.
-<
+The inclusive yields (without any jet requirements) agree within 17\%, which is approximately equal
-<
+to the uncertainty in the measured WZ cross section. A data vs. MC comparison of kinematic
->
+The inclusive yields (without any jet requirements) in the same-flavor final states (660 observed vs. 596 $\pm$ 5.2 MC expected)
->
+agree within 11\%, which is consistent within
->
+the $\approx$15\% uncertainty in the theory prediction for the WZ cross section. A data vs. MC comparison of kinematic
+distributions (jet multiplicity, \MET, Z \pt) is given in Fig.~\ref{fig:wz}. High \MET\
+values in WZ and ZZ events arise from highly boosted W or Z bosons that decay leptonically,
+and we therefore check that the MC does a reasonable job of reproducing the \pt distributions of the
+leptonically decaying \Z. While the inclusive WZ yields are in reasonable agreement, we observe
+an excess in data in events with at least 2 jets, corresponding to the jet multiplicity requirement
-<
+in our preselection. We observe 60 events in data while the MC predicts $34\pm5.2$~(stat), representing an excess of 78\%,
-<
+as indicated in Table~\ref{tab:wz2j}. We note some possible contributions to this discrepancy:
->
+in our preselection. We observe 124 same-flavor events in data while the MC predicts $88\pm1.5$~(stat), representing an excess of 41\%,
->
+as indicated in Table~\ref{tab:wz2j}, and we therefore assess an uncertainty of 50\% on the WZ background.
->
+%We note that the contributions from fake leptons and from \zjets\ with mismeasured \MET\
->
+%is underestimated in the MC.
->
+%This excess will be studied further. For the time being, based on these studies we currently assess an uncertainty of 50\% on the WZ yield.
->
+%A data vs. MC comparison of several kinematic quantities in the sample with 3 leptons and at least 2 jets can be found in App.~\ref{app:WZ}.
-–
+\begin{itemize}
-–
-–
+\item The \zjets\ contribution is under-estimated here, for 2 reasons: first, because the \zjets\
-–
+yield passing a \MET $>$ 50 GeV requirement is under-estimated in MC and second, because the fake
-–
+rate is typically under-estimated in the MC. To get a rough idea for how much the excess depends
-–
+on the \zjets\ yield, if the \zjets\ yield is doubled then the excess is reduced from 78\% to 55\%.
-–
+Also note that we are currently using 10\% of the \zjets\ MC sample and there is 1 event with a weight
-–
+of about 5, so the plots and tables will be remade with full \zjets\ sample.
-–
-–
+\item The \ttbar\ contribution is under-estimated here because the fake
-–
+rate is typically under-estimated in the MC. To get a rough idea for how much the excess depends
-–
+on the \ttbar\ yield, if the \ttbar\ yield is doubled then the excess is reduced from 78\% to 57\%.
-–
-–
+\item Currently no attempt is made to reject jets from pile-up interactions, which may be responsible
-–
+for some of the excess at large \njets. To check this, we increase the jet \pt threhsold to 40 GeV, which
-–
+helps to suppress PU jets, and observe 39 events in data vs. an MC prediction of $25\pm5.2$~(stat),
-–
+decreasing the excess from 78\% to 58\%. In the future this may be improved by explicitly
-–
+requiring the jets to be consistent with originating from the signal primary vertex.
-–
+\end{itemize}
-–
-–
+Based on these studies we currently assess an uncertainty of 80\% on the WZ yield.
+\begin{table}[htb]
+\begin{center}
+\caption{\label{tab:wz} Data and Monte Carlo yields passing the WZ preselection. }
+\begin{tabular}{lccccc}
-+
-+
+%Loading babies at       : ../output/V00-02-13
-+
+%Using selection         : ((((((((leptype==0 && (ee==1 || isdata==0))||(leptype==1 && (mm==1 || isdata==0)))||(leptype==2 && (em==1||me==1||isdata==0)))&&(csc==0 && hbhe==1 && hcallaser==1 && ecaltp==1 && trkfail==1 && eebadsc==1 && hbhenew==1))&&(lep1.pt()>20.0 && lep2.pt()>20.0))&&(nlep==3 && lep3.pt()>20.0))&&(pfmet>50))&&(dilmass>81 && dilmass<101))&&(nbcsvm==0)
-+
+%Using weight            : weight * 19.5 * trgeff * vtxweight
-+
-+
+\hline
+\hline
-+
+         Sample   &             ee   &       $\mu\mu$   &         e$\mu$   &          total  \\
+\hline
-<
+         Sample   &            ee    &        $\mu\mu$   &        e$\mu$   &          total  \\
->
+             WZ   &244.0 $\pm$ 1.6   &301.8 $\pm$ 1.6   & 17.0 $\pm$ 0.4   &562.8 $\pm$ 2.3  \\
->
+             ZZ   & 10.2 $\pm$ 0.1   & 12.8 $\pm$ 0.1   &  0.8 $\pm$ 0.0   & 23.8 $\pm$ 0.1  \\
->
+         \zjets   &  3.3 $\pm$ 2.4   &  5.8 $\pm$ 3.4   &  0.0 $\pm$ 0.0   &  9.0 $\pm$ 4.2  \\
->
+         \ttbar   &  0.6 $\pm$ 0.6   &  4.5 $\pm$ 1.6   &  2.1 $\pm$ 1.1   &  7.3 $\pm$ 2.0  \\
->
+     single top   &  0.0 $\pm$ 0.0   &  0.5 $\pm$ 0.5   &  0.0 $\pm$ 0.0   &  0.5 $\pm$ 0.5  \\
->
+             WW   &  0.0 $\pm$ 0.0   &  0.1 $\pm$ 0.1   &  0.2 $\pm$ 0.1   &  0.3 $\pm$ 0.1  \\
->
+            ttV   &  2.3 $\pm$ 0.2   &  2.6 $\pm$ 0.2   &  0.6 $\pm$ 0.1   &  5.5 $\pm$ 0.3  \\
->
+            VVV   &  3.0 $\pm$ 0.1   &  3.7 $\pm$ 0.1   &  0.6 $\pm$ 0.1   &  7.4 $\pm$ 0.2  \\
+\hline
-<
+             WZ   & 58.9 $\pm$ 0.7   & 82.2 $\pm$ 0.8   &  4.0 $\pm$ 0.2   &145.1 $\pm$ 1.0  \\
-<
+         \ttbar   &  0.6 $\pm$ 0.5   &  4.3 $\pm$ 1.5   &  3.0 $\pm$ 1.2   &  8.0 $\pm$ 2.0  \\
-<
+         \zjets   &  0.4 $\pm$ 0.4   &  4.9 $\pm$ 4.9   &  0.0 $\pm$ 0.0   &  5.3 $\pm$ 4.9  \\
-<
+             ZZ   &  1.4 $\pm$ 0.0   &  2.0 $\pm$ 0.0   &  0.1 $\pm$ 0.0   &  3.5 $\pm$ 0.0  \\
-<
+             WW   &  0.0 $\pm$ 0.0   &  0.2 $\pm$ 0.1   &  0.2 $\pm$ 0.1   &  0.3 $\pm$ 0.1  \\
-<
+     single top   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.1 $\pm$ 0.1  \\
->
+    total SM MC   &263.4 $\pm$ 3.0   &331.8 $\pm$ 4.2   & 21.4 $\pm$ 1.2   &616.6 $\pm$ 5.2  \\
+\hline
-<
+    total SM MC   & 61.3 $\pm$ 0.9   & 93.7 $\pm$ 5.2   &  7.3 $\pm$ 1.3   &162.3 $\pm$ 5.4  \\
-<
+           data   &             68   &            108   &             14   &            190  \\
->
+           data   &            288   &            372   &             36   &            696  \\
+\hline
+\hline
+\begin{table}[htb]
+\begin{center}
-<
+\caption{\label{tab:wz2j} Data and Monte Carlo yields passing the WZ preselection and \njets\ $>$ 2. }
->
+\caption{\label{tab:wz2j} Data and Monte Carlo yields passing the WZ preselection and \njets\ $\geq$ 2. }
+\begin{tabular}{lccccc}
-+
-+
+%Loading babies at       : ../output/V00-02-13
-+
+%Using selection         : (((((((((leptype==0 && (ee==1 || isdata==0))||(leptype==1 && (mm==1 || isdata==0)))||(leptype==2 && (em==1||me==1||isdata==0)))&&(csc==0 && hbhe==1 && hcallaser==1 && ecaltp==1 && trkfail==1 && eebadsc==1 && hbhenew==1))&&(lep1.pt()>20.0 && lep2.pt()>20.0))&&(nlep==3 && lep3.pt()>20.0))&&(pfmet>50))&&(dilmass>81 && dilmass<101))&&(nbcsvm==0))&&(njets>=2)
-+
+%Using weight            : weight * 19.5 * trgeff * vtxweight
-+
+\hline
+\hline
-<
+         Sample   &            ee    &        $\mu\mu$   &        e$\mu$   &          total  \\
->
+         Sample   &             ee   &       $\mu\mu$   &         e$\mu$   &          total  \\
+\hline
-<
+             WZ   &  9.8 $\pm$ 0.3   & 13.3 $\pm$ 0.3   &  0.6 $\pm$ 0.1   & 23.6 $\pm$ 0.4  \\
-<
+         \ttbar   &  0.2 $\pm$ 0.2   &  2.0 $\pm$ 0.9   &  2.2 $\pm$ 1.2   &  4.4 $\pm$ 1.5  \\
-<
+         \zjets   &  0.0 $\pm$ 0.0   &  4.9 $\pm$ 4.9   &  0.0 $\pm$ 0.0   &  4.9 $\pm$ 4.9  \\
-<
+             ZZ   &  0.3 $\pm$ 0.0   &  0.4 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.7 $\pm$ 0.0  \\
-<
+             WW   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.1 $\pm$ 0.0  \\
-<
+     single top   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0  \\
->
+             WZ   & 33.0 $\pm$ 0.6   & 41.1 $\pm$ 0.6   &  2.3 $\pm$ 0.2   & 76.3 $\pm$ 0.9  \\
->
+             ZZ   &  1.7 $\pm$ 0.0   &  2.1 $\pm$ 0.0   &  0.1 $\pm$ 0.0   &  3.9 $\pm$ 0.1  \\
->
+         \zjets   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0  \\
->
+         \ttbar   &  0.6 $\pm$ 0.6   &  1.3 $\pm$ 0.9   &  1.2 $\pm$ 0.9   &  3.1 $\pm$ 1.4  \\
->
+     single top   &  0.0 $\pm$ 0.0   &  0.5 $\pm$ 0.5   &  0.0 $\pm$ 0.0   &  0.5 $\pm$ 0.5  \\
->
+             WW   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0  \\
->
+           \ttV   &  1.8 $\pm$ 0.2   &  2.1 $\pm$ 0.2   &  0.4 $\pm$ 0.1   &  4.3 $\pm$ 0.3  \\
->
+            VVV   &  1.5 $\pm$ 0.1   &  2.0 $\pm$ 0.1   &  0.1 $\pm$ 0.0   &  3.7 $\pm$ 0.1  \\
+\hline
-<
+    total SM MC   & 10.3 $\pm$ 0.3   & 20.8 $\pm$ 5.0   &  2.8 $\pm$ 1.2   & 33.8 $\pm$ 5.2  \\
-<
+           data   &             23   &             32   &              5   &             60  \\
->
+    total SM MC   & 38.7 $\pm$ 0.9   & 49.1 $\pm$ 1.2   &  4.1 $\pm$ 0.9   & 91.8 $\pm$ 1.7  \\
->
+\hline
->
+           data   &             61   &             63   &             11   &            135  \\
+\hline
+\hline
+\begin{figure}[tbh]
+\begin{center}
-<
+\includegraphics[width=1\linewidth]{plots/WZ.pdf}
->
+\includegraphics[width=1\linewidth]{plots/WZ_19p5fb.pdf}
+\caption{\label{fig:wz}\protect
+Data vs. MC comparisons for the WZ selection discussed in the text for \lumi.
+The number of jets, missing transverse energy, and Z boson transverse momentum are displayed.
-+
+\begin{comment}
-+
-+
+Loading babies at       : ../output/V00-02-13
-+
+Using selection         : ((((((((leptype==0 && (ee==1 || isdata==0))||(leptype==1 && (mm==1 || isdata==0)))||(leptype==2 && (em==1||me==1||isdata==0)))&&(csc==0 && hbhe==1 && hcallaser==1 && ecaltp==1 && trkfail==1 && eebadsc==1 && hbhenew==1))&&(lep1.pt()>20.0 && lep2.pt()>20.0))&&(nlep==3 && lep3.pt()>20.0))&&(pfmet>50))&&(dilmass>81 && dilmass<101))&&(nbcsvm==0)
-+
+Using weight            : weight * 19.5 * trgeff * vtxweight
-+
+Plotting var njets flavor sf
-+
+compareDataMC : apply trigeff 1
-+
+MC yield VVV 6.78
-+
+MC yield ttV 4.94
-+
+MC yield WW 0.09
-+
+MC yield single top 0.47
-+
+MC yield ttbar 5.14
-+
+MC yield ZZ 23.01
-+
+MC yield WZ 545.78
-+
+MC yield zjets 9.02
-+
+MC total yield 595.23
-+
+data yield 660
-+
+Plotting var pfmet flavor sf
-+
+compareDataMC : apply trigeff 1
-+
+Warning in <TROOT::Append>: Replacing existing TH1: htemp (Potential memory leak).
-+
+MC yield VVV 6.78
-+
+MC yield ttV 4.94
-+
+MC yield WW 0.09
-+
+MC yield single top 0.47
-+
+MC yield ttbar 5.14
-+
+MC yield ZZ 23.01
-+
+MC yield WZ 545.79
-+
+MC yield zjets 9.02
-+
+MC total yield 595.24
-+
+data yield 660
-+
+Warning in <TROOT::Append>: Replacing existing TH1: datahist (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: mc_0 (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: mc_1 (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: mc_2 (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: mc_3 (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: mc_4 (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: mc_5 (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: mc_6 (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: mc_7 (Potential memory leak).
-+
+Plotting var dileppt flavor sf
-+
+compareDataMC : apply trigeff 1
-+
+Warning in <TROOT::Append>: Replacing existing TH1: htemp (Potential memory leak).
-+
+MC yield VVV 6.78
-+
+MC yield ttV 4.94
-+
+MC yield WW 0.09
-+
+MC yield single top 0.47
-+
+MC yield ttbar 5.14
-+
+MC yield ZZ 23.01
-+
+MC yield WZ 545.78
-+
+MC yield zjets 9.02
-+
+MC total yield 595.23
-+
+data yield 660
-+
-+
+\end{comment}
-+
+\end{center}
+\end{figure}
+\item 2 of the 4 leptons must fall in the $Z$ window 81-101 GeV.
+\end{itemize}
-<
+The data and MC yields passing the above selection are in Table~\ref{tab:zz}. Again we observe an
-<
+excess in data with respect to the MC prediction (29 observed vs. $17.3\pm0.1$~(stat) MC predicted).
-<
+After requiring at least 2 jets, we observe 2 events and the MC predicts $1.5\pm0.1$~(stat).
-<
+However, we have recently discovered that we may be using the wrong (too small) cross section for the ZZ sample,
-<
+and we are in contact with the MC generator group to determine the correct cross section.
-<
+Based on this we currently apply an uncertainty of 80\% to the ZZ background.
->
+The data and MC yields passing the above selection are in Table~\ref{tab:zz}.
->
+In this ZZ-dominated sample we observe good agreement between the data yield and the MC prediction.
->
+After requiring 2 jets (corresponding to the requirement in the analysis selection), we observe 11 events
->
+in data and the MC predicts $11\pm0.2$ events. Due to the limited statistical precision we assign an uncertainty
->
+of 50\% on the ZZ yield.
+\begin{table}[htb]
+\begin{center}
+\caption{\label{tab:zz} Data and Monte Carlo yields for the ZZ preselection. }
+\begin{tabular}{lccccc}
-+
-+
+%Loading babies at       : ../output/V00-02-13
-+
+%Using selection         : ((((((leptype==0 && (ee==1 || isdata==0))||(leptype==1 && (mm==1 || isdata==0)))||(leptype==2 && (em==1||me==1||isdata==0)))&&(csc==0 && hbhe==1 && hcallaser==1 && ecaltp==1 && trkfail==1 && eebadsc==1 && hbhenew==1))&&(lep1.pt()>20.0 && lep2.pt()>20.0))&&(nlep==4 && lep3.pt()>20.0 && lep4.pt()>20.0))&&(dilmass>81 && dilmass<101)
-+
+%Using weight            : weight * 19.5 * trgeff * vtxweight
-+
+\hline
+\hline
+         Sample   &             ee   &       $\mu\mu$   &         e$\mu$   &          total  \\
+\hline
-<
+             ZZ   &  6.6 $\pm$ 0.0   &  9.9 $\pm$ 0.0   &  0.4 $\pm$ 0.0   & 17.0 $\pm$ 0.1  \\
-<
+             WZ   &  0.1 $\pm$ 0.0   &  0.2 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.3 $\pm$ 0.0  \\
->
+             ZZ   & 53.0 $\pm$ 0.2   & 70.7 $\pm$ 0.2   &  3.5 $\pm$ 0.0   &127.2 $\pm$ 0.3  \\
->
+           \ttV   &  1.3 $\pm$ 0.2   &  1.4 $\pm$ 0.2   &  0.3 $\pm$ 0.1   &  3.0 $\pm$ 0.2  \\
->
+            VVV   &  0.7 $\pm$ 0.1   &  0.9 $\pm$ 0.1   &  0.0 $\pm$ 0.0   &  1.7 $\pm$ 0.1  \\
+         \zjets   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0  \\
-+
+             WZ   &  0.1 $\pm$ 0.0   &  0.1 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.3 $\pm$ 0.1  \\
+         \ttbar   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0  \\
-–
+             WW   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0  \\
+     single top   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0  \\
-+
+             WW   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0  \\
-+
+\hline
-+
+    total SM MC   & 55.1 $\pm$ 0.3   & 73.1 $\pm$ 0.3   &  3.9 $\pm$ 0.1   &132.1 $\pm$ 0.4  \\
+\hline
-<
+    total SM MC   &  6.7 $\pm$ 0.0   & 10.1 $\pm$ 0.1   &  0.5 $\pm$ 0.0   & 17.3 $\pm$ 0.1  \\
-<
+           data   &             13   &             16   &              0   &             29  \\
->
+           data   &             57   &             82   &              5   &            144  \\
+\hline
+\hline
-+
+\end{tabular}
+\end{center}
+\end{table}
+\begin{figure}[tbh]
+\begin{center}
-<
+\includegraphics[width=1\linewidth]{plots/ZZ.pdf}
->
+\includegraphics[width=1\linewidth]{plots/ZZ_19p5fb.pdf}
+\caption{\label{fig:zz}\protect
+Data vs. MC comparisons for the ZZ selection discussed in the text for \lumi.
+The number of jets, missing transverse energy, and Z boson transverse momentum are displayed.
+\end{center}
+\end{figure}
-+
+\begin{comment}
-+
-+
+Loading babies at       : ../output/V00-02-13
-+
+Using selection         : ((((((leptype==0 && (ee==1 || isdata==0))||(leptype==1 && (mm==1 || isdata==0)))||(leptype==2 && (em==1||me==1||isdata==0)))&&(csc==0 && hbhe==1 && hcallaser==1 && ecaltp==1 && trkfail==1 && eebadsc==1 && hbhenew==1))&&(lep1.pt()>20.0 && lep2.pt()>20.0))&&(nlep==4 && lep3.pt()>20.0 && lep4.pt()>20.0))&&(dilmass>81 && dilmass<101)
-+
+Using weight            : weight * 19.5 * trgeff * vtxweight
-+
+Plotting var njets flavor sf
-+
+compareDataMC : apply trigeff 1
-+
+MC yield VVV 1.63
-+
+MC yield ttV 2.66
-+
+MC yield WW 0.00
-+
+MC yield single top 0.00
-+
+MC yield ttbar 0.00
-+
+MC yield ZZ 123.68
-+
+MC yield WZ 0.25
-+
+MC yield zjets 0.00
-+
+MC total yield 128.21
-+
+data yield 139
-+
+Plotting var pfmet flavor sf
-+
+compareDataMC : apply trigeff 1
-+
+MC yield VVV 1.63
-+
+MC yield ttV 2.66
-+
+MC yield WW 0.00
-+
+MC yield single top 0.00
-+
+MC yield ttbar 0.00
-+
+MC yield ZZ 123.66
-+
+MC yield WZ 0.25
-+
+MC yield zjets 0.00
-+
+MC total yield 128.19
-+
+data yield 139
-+
+Plotting var dileppt flavor sf
-+
+compareDataMC : apply trigeff 1
-+
+MC yield VVV 1.63
-+
+MC yield ttV 2.66
-+
+MC yield WW 0.00
-+
+MC yield single top 0.00
-+
+MC yield ttbar 0.00
-+
+MC yield ZZ 123.67
-+
+MC yield WZ 0.25
-+
+MC yield zjets 0.00
-+
+MC total yield 128.20
-+
+data yield 139
-+
-+
-+
+\end{comment}

Diff Legend

-–
+Removed lines
-+
+Added lines
-<
+Changed lines
->
+Changed lines

Comparing UserCode/benhoob/cmsnotes/ZMet2012/bkg.tex (file contents): Revision 1.4 by benhoob, Fri Jun 29 19:41:46 2012 UTC vs. Revision 1.20 by benhoob, Fri Jan 25 17:46:18 2013 UTC

Diff Legend

Comparing UserCode/benhoob/cmsnotes/ZMet2012/bkg.tex (file contents):
Revision 1.4 by benhoob, Fri Jun 29 19:41:46 2012 UTC vs.
Revision 1.20 by benhoob, Fri Jan 25 17:46:18 2013 UTC