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

Comparing UserCode/benhoob/cmsnotes/ZMet2012/bkg.tex (file contents):
Revision 1.5 by benhoob, Wed Sep 12 08:24:01 2012 UTC vs.
Revision 1.6 by benhoob, Fri Sep 14 19:59:12 2012 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{144122/0.88}{110325/0.95}} = 1.19\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.92}\frac{0.95}{2\times1.26} = (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.88 \times 1.26}{2} = (0.55\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:
+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.09$ 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$.
->
+for all \MET\ regions up to  \MET\ $>$ 150 GeV. 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_92fb.pdf} &
->
+\includegraphics[width=0.4\textwidth]{plots/extractK_exclusive_92fb.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.
-+
-+
+\begin{comment}
-+
-+
+----------------------------------------
-+
+  EXCLUSIVE RESULTS
-+
+----------------------------------------
-+
-+
+root [3] extractK(true,false,false)
-+
+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)  21691
-+
+OF entries (Z mass) 2934
-+
+K                   0.135263
-+
+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 : 3650
-+
+Z     : 461
-+
+K     : 0.13 +/- 0.006
-+
-+
+MC    :
-+
+total : 399.019
-+
+Z     : 51.0493
-+
+K     : 0.13 +/- 0.007
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30  && pfmet<60
-+
-+
+data  :
-+
+total : 6951
-+
+Z     : 904
-+
+K     : 0.13 +/- 0.004
-+
-+
+MC    :
-+
+total : 755.309
-+
+Z     : 111.206
-+
+K     : 0.15 +/- 0.003
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60  && pfmet<100
-+
-+
+data  :
-+
+total : 7206
-+
+Z     : 1037
-+
+K     : 0.14 +/- 0.004
-+
-+
+MC    :
-+
+total : 838.418
-+
+Z     : 123.554
-+
+K     : 0.15 +/- 0.003
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100 && pfmet<200
-+
-+
+data  :
-+
+total : 3690
-+
+Z     : 512
-+
+K     : 0.14 +/- 0.006
-+
-+
+MC    :
-+
+total : 451.624
-+
+Z     : 67.7098
-+
+K     : 0.15 +/- 0.004
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200 && pfmet<300
-+
-+
+data  :
-+
+total : 172
-+
+Z     : 17
-+
+K     : 0.10 +/- 0.024
-+
-+
+MC    :
-+
+total : 24.2441
-+
+Z     : 2.67077
-+
+K     : 0.11 +/- 0.013
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>300
-+
-+
+data  :
-+
+total : 22
-+
+Z     : 3
-+
+K     : 0.14 +/- 0.079
-+
-+
+MC    :
-+
+total : 4.53108
-+
+Z     : 0.230071
-+
+K     : 0.05 +/- 0.022
-+
+--------------------------------------------------------------
-+
-+
-+
-+
+----------------------------------------
-+
+  INCLUSIVE RESULTS
-+
+----------------------------------------
-+
-+
+root [4] extractK(false,false,false)
-+
+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)  21691
-+
+OF entries (Z mass) 2934
-+
+K                   0.135263
-+
+Warning in <TROOT::Append>: Replacing existing TH1: htot (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: hZ (Potential memory leak).
-+
-+
+--------------------------------------------------------------
-+
+pfmet>0
-+
-+
+data  :
-+
+total : 21691
-+
+Z     : 2934
-+
+K     : 0.14 +/- 0.002
-+
-+
+MC    :
-+
+total : 2472.89
-+
+Z     : 356.434
-+
+K     : 0.14 +/- 0.002
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30
-+
-+
+data  :
-+
+total : 18041
-+
+Z     : 2473
-+
+K     : 0.14 +/- 0.003
-+
-+
+MC    :
-+
+total : 2074.05
-+
+Z     : 305.382
-+
+K     : 0.15 +/- 0.002
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60
-+
-+
+data  :
-+
+total : 11090
-+
+Z     : 1569
-+
+K     : 0.14 +/- 0.004
-+
-+
+MC    :
-+
+total : 1318.79
-+
+Z     : 194.166
-+
+K     : 0.15 +/- 0.002
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100
-+
-+
+data  :
-+
+total : 3884
-+
+Z     : 532
-+
+K     : 0.14 +/- 0.006
-+
-+
+MC    :
-+
+total : 480.402
-+
+Z     : 70.6107
-+
+K     : 0.15 +/- 0.004
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200
-+
-+
+data  :
-+
+total : 194
-+
+Z     : 20
-+
+K     : 0.10 +/- 0.023
-+
-+
+MC    :
-+
+total : 28.7751
-+
+Z     : 2.90084
-+
+K     : 0.10 +/- 0.012
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>300
-+
-+
+data  :
-+
+total : 22
-+
+Z     : 3
-+
+K     : 0.14 +/- 0.079
-+
-+
+MC    :
-+
+total : 4.53108
-+
+Z     : 0.230071
-+
+K     : 0.05 +/- 0.022
-+
+--------------------------------------------------------------
-+
-+
+\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_bvetoMedium_92fb.pdf} &
->
+\includegraphics[width=0.4\textwidth]{plots/extractK_exclusive_bvetoMedium_92fb.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
+%{\bf FIXME plots made with 10\% of \zjets\ MC statistics, to be remade with full statistics}
+\label{fig:K_targeted}
-+
+\begin{comment}
-+
-+
+root [2] extractK(true,false,true)
-+
+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)  5756
-+
+OF entries (Z mass) 654
-+
+K                   0.113621
-+
+Warning in <TStreamerInfo::Compile>: Counter fNClusterRange should not be skipped from class TTree
-+
+Info in <TCanvas::MakeDefCanvas>:  created default TCanvas with name c1
-+
-+
+--------------------------------------------------------------
-+
+pfmet>0   && pfmet<30
-+
-+
+data  :
-+
+total : 1303
-+
+Z     : 126
-+
+K     : 0.10 +/- 0.009
-+
-+
+MC    :
-+
+total : 131.974
-+
+Z     : 15.1946
-+
+K     : 0.12 +/- 0.020
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30  && pfmet<60
-+
-+
+data  :
-+
+total : 1818
-+
+Z     : 190
-+
+K     : 0.10 +/- 0.008
-+
-+
+MC    :
-+
+total : 172.956
-+
+Z     : 21.9369
-+
+K     : 0.13 +/- 0.007
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60  && pfmet<80
-+
-+
+data  :
-+
+total : 994
-+
+Z     : 122
-+
+K     : 0.12 +/- 0.011
-+
-+
+MC    :
-+
+total : 109.789
-+
+Z     : 13.9792
-+
+K     : 0.13 +/- 0.008
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>80  && pfmet<100
-+
-+
+data  :
-+
+total : 699
-+
+Z     : 99
-+
+K     : 0.14 +/- 0.014
-+
-+
+MC    :
-+
+total : 73.3643
-+
+Z     : 11.5154
-+
+K     : 0.16 +/- 0.010
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100 && pfmet<150
-+
-+
+data  :
-+
+total : 722
-+
+Z     : 93
-+
+K     : 0.13 +/- 0.013
-+
-+
+MC    :
-+
+total : 86.7947
-+
+Z     : 11.742
-+
+K     : 0.14 +/- 0.009
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>150 && pfmet<200
-+
-+
+data  :
-+
+total : 163
-+
+Z     : 18
-+
+K     : 0.11 +/- 0.026
-+
-+
+MC    :
-+
+total : 19.4473
-+
+Z     : 2.97965
-+
+K     : 0.15 +/- 0.017
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200
-+
-+
+data  :
-+
+total : 57
-+
+Z     : 6
-+
+K     : 0.11 +/- 0.043
-+
-+
+MC    :
-+
+total : 8.99801
-+
+Z     : 0.794136
-+
+K     : 0.09 +/- 0.021
-+
+--------------------------------------------------------------
-+
-+
+root [3] Info in <TCanvas::Print>: pdf file /Users/benhoob/tas/ZMet2012/plots/extractK_exclusive_bvetoMedium_92fb.pdf has been created
-+
-+
+root [3]
-+
+root [3] extractK(false,false,true)
-+
+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)  5756
-+
+OF entries (Z mass) 654
-+
+K                   0.113621
-+
+Warning in <TROOT::Append>: Replacing existing TH1: htot (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: hZ (Potential memory leak).
-+
-+
+--------------------------------------------------------------
-+
+pfmet>0
-+
-+
+data  :
-+
+total : 5756
-+
+Z     : 654
-+
+K     : 0.11 +/- 0.004
-+
-+
+MC    :
-+
+total : 603.333
-+
+Z     : 78.1422
-+
+K     : 0.13 +/- 0.005
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30
-+
-+
+data  :
-+
+total : 4453
-+
+Z     : 528
-+
+K     : 0.12 +/- 0.005
-+
-+
+MC    :
-+
+total : 471.396
-+
+Z     : 62.9476
-+
+K     : 0.13 +/- 0.004
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60
-+
-+
+data  :
-+
+total : 2635
-+
+Z     : 338
-+
+K     : 0.13 +/- 0.007
-+
-+
+MC    :
-+
+total : 298.41
-+
+Z     : 41.0107
-+
+K     : 0.14 +/- 0.005
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>80
-+
-+
+data  :
-+
+total : 1641
-+
+Z     : 216
-+
+K     : 0.13 +/- 0.009
-+
-+
+MC    :
-+
+total : 188.602
-+
+Z     : 27.0313
-+
+K     : 0.14 +/- 0.006
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100
-+
-+
+data  :
-+
+total : 942
-+
+Z     : 117
-+
+K     : 0.12 +/- 0.011
-+
-+
+MC    :
-+
+total : 115.24
-+
+Z     : 15.5158
-+
+K     : 0.13 +/- 0.008
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>150
-+
-+
+data  :
-+
+total : 220
-+
+Z     : 24
-+
+K     : 0.11 +/- 0.022
-+
-+
+MC    :
-+
+total : 28.4454
-+
+Z     : 3.77378
-+
+K     : 0.13 +/- 0.014
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200
-+
-+
+data  :
-+
+total : 57
-+
+Z     : 6
-+
+K     : 0.11 +/- 0.043
-+
-+
+MC    :
-+
+total : 8.99801
-+
+Z     : 0.794136
-+
+K     : 0.09 +/- 0.021
-+
+--------------------------------------------------------------
-+
-+
+\end{comment}
-+
-+
+\end{center}
-+
+\end{figure}
-+
-+
-+
+\begin{comment}
-+
-+
+\begin{figure}[!hb]
-+
+\begin{center}
-+
+\begin{tabular}{cc}
-+
+\includegraphics[width=0.4\textwidth]{plots/extractK_inclusive_bvetoLoose_92fb.pdf} &
-+
+\includegraphics[width=0.4\textwidth]{plots/extractK_exclusive_bvetoLoose_92fb.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.
-+
+%{\bf FIXME plots made with 10\% of \zjets\ MC statistics, to be remade with full statistics}
-+
+\label{fig:K_targeted}}
-+
-+
-+
+root [2] extractK(true,false,true)
-+
+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))&&(nbcsvl==0)
-+
+Using weight    : vtxweight * weight
-+
+OF entries (total)  2715
-+
+OF entries (Z mass) 279
-+
+K                   0.102762
-+
+Warning in <TStreamerInfo::Compile>: Counter fNClusterRange should not be skipped from class TTree
-+
+Info in <TCanvas::MakeDefCanvas>:  created default TCanvas with name c1
-+
-+
+--------------------------------------------------------------
-+
+pfmet>0   && pfmet<30
-+
-+
+data  :
-+
+total : 713
-+
+Z     : 59
-+
+K     : 0.08 +/- 0.011
-+
-+
+MC    :
-+
+total : 74.2549
-+
+Z     : 7.09789
-+
+K     : 0.10 +/- 0.025
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30  && pfmet<60
-+
-+
+data  :
-+
+total : 850
-+
+Z     : 79
-+
+K     : 0.09 +/- 0.010
-+
-+
+MC    :
-+
+total : 84.6973
-+
+Z     : 9.55105
-+
+K     : 0.11 +/- 0.009
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60  && pfmet<80
-+
-+
+data  :
-+
+total : 469
-+
+Z     : 61
-+
+K     : 0.13 +/- 0.017
-+
-+
+MC    :
-+
+total : 50.1496
-+
+Z     : 5.92081
-+
+K     : 0.12 +/- 0.012
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>80  && pfmet<100
-+
-+
+data  :
-+
+total : 269
-+
+Z     : 33
-+
+K     : 0.12 +/- 0.021
-+
-+
+MC    :
-+
+total : 30.0547
-+
+Z     : 4.67993
-+
+K     : 0.16 +/- 0.014
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100 && pfmet<150
-+
-+
+data  :
-+
+total : 311
-+
+Z     : 34
-+
+K     : 0.11 +/- 0.019
-+
-+
+MC    :
-+
+total : 39.4475
-+
+Z     : 5.02593
-+
+K     : 0.13 +/- 0.014
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>150 && pfmet<200
-+
-+
+data  :
-+
+total : 79
-+
+Z     : 10
-+
+K     : 0.13 +/- 0.040
-+
-+
+MC    :
-+
+total : 9.96228
-+
+Z     : 1.4975
-+
+K     : 0.15 +/- 0.023
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200
-+
-+
+data  :
-+
+total : 24
-+
+Z     : 3
-+
+K     : 0.12 +/- 0.072
-+
-+
+MC    :
-+
+total : 5.3503
-+
+Z     : 0.425719
-+
+K     : 0.08 +/- 0.027
-+
+--------------------------------------------------------------
-+
-+
+root [3] Info in <TCanvas::Print>: pdf file /Users/benhoob/tas/ZMet2012/plots/extractK_exclusive_bvetoLoose_92fb.pdf has been created
-+
-+
+root [3]
-+
+root [3] extractK(false,false,true)
-+
+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))&&(nbcsvl==0)
-+
+Using weight    : vtxweight * weight
-+
+OF entries (total)  2715
-+
+OF entries (Z mass) 279
-+
+K                   0.102762
-+
+Warning in <TROOT::Append>: Replacing existing TH1: htot (Potential memory leak).
-+
+Warning in <TROOT::Append>: Replacing existing TH1: hZ (Potential memory leak).
-+
-+
+--------------------------------------------------------------
-+
+pfmet>0
-+
-+
+data  :
-+
+total : 2715
-+
+Z     : 279
-+
+K     : 0.10 +/- 0.006
-+
-+
+MC    :
-+
+total : 293.912
-+
+Z     : 34.199
-+
+K     : 0.12 +/- 0.008
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>30
-+
-+
+data  :
-+
+total : 2002
-+
+Z     : 220
-+
+K     : 0.11 +/- 0.007
-+
-+
+MC    :
-+
+total : 219.661
-+
+Z     : 27.101
-+
+K     : 0.12 +/- 0.006
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>60
-+
-+
+data  :
-+
+total : 1152
-+
+Z     : 141
-+
+K     : 0.12 +/- 0.010
-+
-+
+MC    :
-+
+total : 134.962
-+
+Z     : 17.5498
-+
+K     : 0.13 +/- 0.007
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>80
-+
-+
+data  :
-+
+total : 683
-+
+Z     : 80
-+
+K     : 0.12 +/- 0.013
-+
-+
+MC    :
-+
+total : 84.8149
-+
+Z     : 11.629
-+
+K     : 0.14 +/- 0.009
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>100
-+
-+
+data  :
-+
+total : 414
-+
+Z     : 47
-+
+K     : 0.11 +/- 0.017
-+
-+
+MC    :
-+
+total : 54.7604
-+
+Z     : 6.94915
-+
+K     : 0.13 +/- 0.011
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>150
-+
-+
+data  :
-+
+total : 103
-+
+Z     : 13
-+
+K     : 0.13 +/- 0.035
-+
-+
+MC    :
-+
+total : 15.3125
-+
+Z     : 1.92322
-+
+K     : 0.13 +/- 0.019
-+
+--------------------------------------------------------------
-+
-+
-+
+--------------------------------------------------------------
-+
+pfmet>200
-+
-+
+data  :
-+
+total : 24
-+
+Z     : 3
-+
+K     : 0.12 +/- 0.072
-+
-+
+MC    :
-+
+total : 5.3503
-+
+Z     : 0.425719
-+
+K     : 0.08 +/- 0.027
-+
+--------------------------------------------------------------
-+
-+
+\end{center}
+\end{figure}
-+
-+
+\end{comment}
-+
-+
+\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}
+\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) agree within 13\%, which is consistent within
->
+the uncertainty in the CMS measured WZ cross section (17\%). 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\%,
->
+in our preselection. We observe 106 events in data while the MC predicts $62\pm1.5$~(stat), representing an excess of 71\%,
+as indicated in Table~\ref{tab:wz2j}. We note some possible contributions to this discrepancy:
+\begin{itemize}
-+
+\item {\bf The following checks refer to the 5.2 fb$^{-1}$ results and will be updated.}
-+
+\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
+\end{itemize}
-<
+Based on these studies we currently assess an uncertainty of 80\% on the WZ yield.
->
+Based on these studies we currently assess an uncertainty of 70\% on the WZ yield.
+\begin{table}[htb]
+\begin{center}
+\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 4 events
->
+in data and the MC predicts $6.6\pm0.1$ events. Due to the limited statistical precision we assign an uncertainty
->
+fo 50\% on the ZZ yield.
+\begin{table}[htb]
+\begin{center}
+\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  \\
->
->
+%Loading babies at       : ../output/V00-01-04
->
+%Using selection         : (((((isdata==0 || (run<197556 || run>198913))&&(((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 * 9.2 * trgeff * vtxweight
->
+%SCALING ZJETS BY 111/946
->
+%SCALING ZZ BY 1.92
->
->
+             ZZ   & 25.1 $\pm$ 0.1   & 34.9 $\pm$ 0.1   &  1.6 $\pm$ 0.0   & 61.7 $\pm$ 0.1  \\
->
+            ttV   &  0.6 $\pm$ 0.1   &  0.6 $\pm$ 0.1   &  0.2 $\pm$ 0.0   &  1.4 $\pm$ 0.1  \\
->
+            VVV   &  0.3 $\pm$ 0.0   &  0.4 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.7 $\pm$ 0.0  \\
->
+             WZ   &  0.1 $\pm$ 0.0   &  0.1 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.1 $\pm$ 0.0  \\
+         \zjets   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0   &  0.0 $\pm$ 0.0  \\
+         \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
-+
+      tot SM MC   & 26.1 $\pm$ 0.1   & 36.1 $\pm$ 0.1   &  1.8 $\pm$ 0.0   & 63.9 $\pm$ 0.2  \\
+\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   &             24   &             36   &              0   &             60  \\
+\hline
+\hline
+\end{tabular}
+\begin{figure}[tbh]
+\begin{center}
-<
+\includegraphics[width=1\linewidth]{plots/ZZ.pdf}
->
+\includegraphics[width=1\linewidth]{plots/ZZ_92fb.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.

Diff Legend

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

Comparing UserCode/benhoob/cmsnotes/ZMet2012/bkg.tex (file contents): Revision 1.5 by benhoob, Wed Sep 12 08:24:01 2012 UTC vs. Revision 1.6 by benhoob, Fri Sep 14 19:59:12 2012 UTC

Diff Legend

Comparing UserCode/benhoob/cmsnotes/ZMet2012/bkg.tex (file contents):
Revision 1.5 by benhoob, Wed Sep 12 08:24:01 2012 UTC vs.
Revision 1.6 by benhoob, Fri Sep 14 19:59:12 2012 UTC