Annotation of ttbar/p20_taujets_note/Systematics.tex, revision 1.1
1.1 ! uid12904 1: \newpage
! 2: \section{Systematic uncertainties}
! 3:
! 4: Several factors contribute to systematic uncertainties in the measurement. Here we describe such uncertainties.
! 5:
! 6: \subsection{JES}
! 7:
! 8: \noindent The jet energy scale (JES) systematic is determined by shifiting the jet energy scale
! 9: by $\pm 1 \sigma$ in all MC samples.
! 10:
! 11:
! 12: \subsection{TES}
! 13:
! 14: \noindent The tau energy scale (TES) systematic is determined by shifiting the tau energy scale
! 15: by its uncertainty as given in \cite{tes_sys}.
! 16:
! 17:
! 18: \subsection{Jet Energy Resolution and Jet ID}
! 19:
! 20: \noindent The jet energy resolution (JER) systematic is determined by shifiting the jet energy
! 21: resolution by $\pm 1 \sigma$ in all MC samples.
! 22:
! 23: \subsection{Trigger}
! 24: \noindent Each event was corrected by the ratio of the actual and predicted trigger result as a function of $H_{T}$, which was
! 25: used based on the fact that the agreement varies as function of it.
! 26:
! 27:
! 28: \subsection{b-quark fragmentation}
! 29: \noindent This uncertainty is estimated using the standard procedure described in \cite{bfrag} by reweighting $t\bar{t}$ events
! 30: using different fragmentation functions.
! 31:
! 32: \subsection{\boldmath $b$-tagging}
! 33:
! 34: \noindent b-tagging uncertainty effects are taken into account by varying the
! 35: systematic and statistical errors on the MC tagging weights.
! 36:
! 37: These errors (which are computed using standard D\O\ b ID group tools) arise form several independent sources \cite{bID-p20}:
! 38:
! 39: \begin{itemize}
! 40: \item B-jet tagging parameterization.
! 41: \item C-jet tagging parameterization.
! 42: \item Light jet tagging parameterization (negative tag rate).
! 43: %\item Systematic uncertainties on the scale factors $SF_{hf}$ and $SF_{ll}$
! 44: %are derived from the statistical error due to finite MC statistics.
! 45: \item Semi-leptonic b-tagging efficiency parameterization in MC and in data
! 46: (System 8).
! 47: \item Taggability. This includes the statistical error due to finite statistic
! 48: in the samples from which it had been derived and systematic, reflecting
! 49: the (neglected) taggability dependence on the jet multiplicity.
! 50: \end{itemize}
! 51:
! 52:
! 53: \subsection{\boldmath $\tau$ ID systematics}
! 54:
! 55: \noindent Here we include systematics associated to the NN cut (NN $>$ 0.90 for taus types 1 and 2 and NN $>$ 0.95 for taus type 3)
! 56: applied to select hadronic taus. As recommended by the $\tau $-ID group
! 57: these systematics are 9.5\%, 3.5\% and 5.0\% for taus type 1, 2 and 3 respectively. However in this analysis we chose
! 58: treat taus types 1 and 2 together. This led us to combine their uncertainties in the following way
! 59:
! 60: \begin{center}
! 61: \begin{equation}
! 62: sys_{12} = \displaystyle \sqrt{\epsilon_{1}^{2}*f_{1}^{2} + \epsilon_{2}^{2}*f_{2}^{2}}
! 63: \end{equation}
! 64: \end{center}
! 65:
! 66: \noindent where $\epsilon_{1}$ and $\epsilon_{2}$ are the $\tau$ ID efficiencies for taus types 1 and 2 respectively and
! 67: $f_{1}$ and $f_{2}$ are the fractions of taus types 1 (0.16) and 2 (0.84) respectively.
! 68:
! 69:
! 70: \subsection{\label{sub:qcd_syst}QCD systematics}
! 71: As explained in the section \ref{sub:Variables} we use the control (b-veto) data set to
! 72: validate our method of modeling the multijet background. Therefore we have to
! 73: use the same sample to evaluate the associated uncertainties. The way it was
! 74: done is by reweighting topological event NN for QCD template
! 75: (``loose-tight'' $\tau$), so that it matches the one for ``tight'' $\tau$ data exactly
! 76: (electroweak beckgrounds were subtracted).
! 77: %Figure \ref{fig:qcd_reweight} shows that this scaling is close to 1, as it should be, since QCD template
! 78: %models the QCD-dominated data very well.
! 79:
! 80: \subsection{W and Z scale factors}
! 81: We apply a ascale factor of 1.47 to both W + bb and W + cc events with an uncertainty of 15\%. At the same time a scales factors
! 82: of 1.52 and 1.67 are applied to Z + bb and Z + cc events, both with an uncertainty of 20\%.
! 83:
! 84:
! 85: \subsection{Template statistics}
! 86: When we performed the template fit to data (Section \ref{sub:NN}) the QCD template had limited statistics
! 87: (1132 events for taus types 1 and 2 and 4487 events for taus type 3). We have to take the
! 88: statistical uncertainty in this histogram as one of the cross section systematics. It was calculated by
! 89: varying the content of each bin of the QCD template NN distribution within its uncertainty and observing
! 90: how the cross section result changed.
! 91:
! 92: \subsection{$t\bar{t}$ contamination in the loose-tight sample}
! 93: When measuring the cross-section we had to take into account the signal contamination in the loose-tight
! 94: sample we use to model QCD in the high NN region. The systematic uncertainty in this case
! 95: is calculated by varying the final assumed cross section by $\pm 1 \sigma$, re-estimating the signal contamination
! 96: and finally measuring the up and down values of the cross section.
! 97:
! 98: \subsection{PDF}
! 99: Systematics on Parton Distribution Functions (PDF) are estimated by reweighting signal $t\bar{t}$ MC from
! 100: CTEQ6L1 to CTEQ6.1m and its twenty error PDF's. The reweighting of the PDF's is done by using {\tt caf\_pdfreweight}
! 101: package tool on Pythia $t\bar{t}$ MC. We then assigned the relative PDF uncertainty obtained with Pythia on the Alpgen
! 102: $t\bar{t}$ MC.
! 103:
! 104: \subsection{Luminosity}
! 105: Here we take the D\O\ standard measured uncertainty on luminosity of 6.1$\%$ .
! 106:
! 107:
! 108: Tables \ref{cap:Syst1} and \ref{cap:Syst2} summarize all of these uncertainty sources and shows how the resulting cross section shifts.
! 109:
! 110: %\clearpage
! 111: %
! 112:
! 113: %\section{Summary \label{sec:summary}}
! 114:
! 115:
! 116: %
! 117:
! 118: %\section{Summary \label{sec:summary}}
! 119:
! 120: \begin{table}[h]
! 121: \caption{Systematic uncertainties on $\sigma(t\bar{t})$ (in pb) for NNelec $>$ 0.9.}
! 122: %\begin{ruledtabular}
! 123: {\footnotesize }\begin{tabular}{cccc}
! 124: \hline
! 125: Channel&
! 126: {\footnotesize $\tau$+jets types 1 and 2 }&
! 127: {\footnotesize $\tau$+jets type 3 }&
! 128: {\footnotesize Combined }\\
! 129: {\footnotesize Tau Energy Scale }&
! 130: {\footnotesize $_{+0.068, -0.102}$ }&
! 131: {\footnotesize $_{+0.340, -0.306}$ }&
! 132: {\footnotesize $_{+0.136, -0.136}$ }\\
! 133: {\footnotesize Jet Energy Scale }&
! 134: {\footnotesize $_{+0.051, -0.034}$ }&
! 135: {\footnotesize $_{+0.051, -0.085}$ }&
! 136: {\footnotesize $_{+0.051, -0.000}$ }\\
! 137: {\footnotesize Jet Energy Resolution }&
! 138: {\footnotesize $_{+0.102, -0.051}$ }&
! 139: {\footnotesize $_{+0.204, -0.034}$ }&
! 140: {\footnotesize $_{+0.119, -0.052}$ }\\
! 141: {\footnotesize Jet ID }&
! 142: {\footnotesize $_{+0.204, -0.204}$ }&
! 143: {\footnotesize $_{+0.153, -0.153}$ }&
! 144: {\footnotesize $_{+0.204, -0.204}$ }\\
! 145: {\footnotesize b-tag }&
! 146: {\footnotesize $_{+0.562, -0.493}$ }&
! 147: {\footnotesize $_{+0.493, -0.426}$ }&
! 148: {\footnotesize $_{+0.544, -0.477}$ }\\
! 149: {\footnotesize b-fragmentation }&
! 150: {\footnotesize $_{+0.102, -0.102}$ }&
! 151: {\footnotesize $_{+0.068, -0.068}$ }&
! 152: {\footnotesize $_{+0.085, -0.085}$ }\\
! 153: {\footnotesize QCD Modeling }&
! 154: {\footnotesize $_{+0.340, -0.340}$ }&
! 155: {\footnotesize $_{+0.221, -0.221}$ }&
! 156: {\footnotesize $_{+0.324, -0.305}$ }\\
! 157: {\footnotesize $\tau$ ID }&
! 158: {\footnotesize $_{+0.272, -0.272}$ }&
! 159: {\footnotesize $_{+0.306, -0.306}$ }&
! 160: {\footnotesize $_{+0.290, -0.290}$ }\\
! 161: {\footnotesize Trigger }&
! 162: {\footnotesize $_{+0.256, -0.256}$ }&
! 163: {\footnotesize $_{+0.238, -0.238}$ }&
! 164: {\footnotesize $_{+0.256, -0.256}$ }\\
! 165: %{\footnotesize $\tau$ triggering }&
! 166: %{\footnotesize $_{+xxxx, -xxxx}$ }&
! 167: %{\footnotesize $_{+xxxx, -xxxx}$ }&
! 168: %{\footnotesize $_{+xxxx, -xxxx}$ }\\
! 169: {\footnotesize W Scale Factor }&
! 170: {\footnotesize $_{+0.034, -0.034}$ }&
! 171: {\footnotesize $_{+0.034, -0.034}$ }&
! 172: {\footnotesize $_{+0.034, -0.034}$ }\\
! 173: {\footnotesize Z Scale Factor }&
! 174: {\footnotesize $_{+0.072, -0.072}$ }&
! 175: {\footnotesize $_{+0.072, -0.072}$ }&
! 176: {\footnotesize $_{+0.048, -0.048}$ }\\
! 177: {\footnotesize Template statistics }&
! 178: {\footnotesize $_{+0.156, -0.156}$ }&
! 179: {\footnotesize $_{+0.204, -0.204}$ }&
! 180: {\footnotesize $_{+0.168, -0.168}$ }\\
! 181: {\footnotesize Signal contamination}&
! 182: {\footnotesize $_{+0.153, -0.153}$ }&
! 183: {\footnotesize $_{+0.255, -0.272}$ }&
! 184: {\footnotesize $_{+0.188, -0.170}$ }\\
! 185: {\footnotesize PDF }&
! 186: {\footnotesize $_{+0.097, -0.084}$ }&
! 187: {\footnotesize $_{+0.188, -0.198}$ }&
! 188: {\footnotesize $_{+0.092, -0.081}$ }\\
! 189:
! 190: \end{tabular}{\footnotesize \par}
! 191: %\end{ruledtabular}
! 192: \label{cap:Syst1}
! 193: \end{table}
! 194:
! 195:
! 196: \begin{table}[h]
! 197: \caption{Systematic uncertainties on $\sigma(t\bar{t})$ (in pb) when no NNelec cut is applied.}
! 198: %\begin{ruledtabular}
! 199: {\footnotesize }\begin{tabular}{cccc}
! 200: \hline
! 201: Channel&
! 202: {\footnotesize $\tau$+jets types 1 and 2 }&
! 203: {\footnotesize $\tau$+jets type 3 }&
! 204: {\footnotesize Combined }\\
! 205: {\footnotesize Tau Energy Scale }&
! 206: {\footnotesize $_{+0.101, -0.002}$ }&
! 207: {\footnotesize $_{+0.238, -0.255}$ }&
! 208: {\footnotesize $_{+0.102, -0.017}$ }\\
! 209: {\footnotesize Jet Energy Scale }&
! 210: {\footnotesize $_{+0.016, -0.001}$ }&
! 211: {\footnotesize $_{+0.017, -0.000}$ }&
! 212: {\footnotesize $_{+0.016, -0.000}$ }\\
! 213: {\footnotesize Jet Energy Resolution }&
! 214: {\footnotesize $_{+0.084, -0.086}$ }&
! 215: {\footnotesize $_{+0.017, -0.034}$ }&
! 216: {\footnotesize $_{+0.068, -0.085}$ }\\
! 217: {\footnotesize Jet ID }&
! 218: {\footnotesize $_{+0.169, -0.169}$ }&
! 219: {\footnotesize $_{+0.017, -0.017}$ }&
! 220: {\footnotesize $_{+0.153, -0.153}$ }\\
! 221: {\footnotesize b-tag }&
! 222: {\footnotesize $_{+0.424, -0.375}$ }&
! 223: {\footnotesize $_{+0.358, -0.306}$ }&
! 224: {\footnotesize $_{+0.425, -0.375}$ }\\
! 225: {\footnotesize b-fragmentation }&
! 226: {\footnotesize $_{+0.069, -0.069}$ }&
! 227: {\footnotesize $_{+0.102, -0.102}$ }&
! 228: {\footnotesize $_{+0.068, -0.068}$ }\\
! 229: {\footnotesize QCD Modeling }&
! 230: {\footnotesize $_{+0.271, -0.273}$ }&
! 231: {\footnotesize $_{+0.153, -0.136}$ }&
! 232: {\footnotesize $_{+0.225, -0.256}$ }\\
! 233: {\footnotesize $\tau$ ID }&
! 234: {\footnotesize $_{+0.220, -0.220}$ }&
! 235: {\footnotesize $_{+0.204, -0.204}$ }&
! 236: {\footnotesize $_{+0.221, -0.221}$ }\\
! 237: {\footnotesize Trigger }&
! 238: {\footnotesize $_{+0.204, -0.204}$ }&
! 239: {\footnotesize $_{+0.170, -0.170}$ }&
! 240: {\footnotesize $_{+0.204, -0.204}$ }\\
! 241: %{\footnotesize $\tau$ triggering }&
! 242: %{\footnotesize $_{+xxxx, -xxxx}$ }&
! 243: %{\footnotesize $_{+xxxx, -xxxx}$ }&
! 244: %{\footnotesize $_{+xxxx, -xxxx}$ }\\
! 245: {\footnotesize W Scale Factor }&
! 246: {\footnotesize $_{+0.034, -0.034}$ }&
! 247: {\footnotesize $_{+0.034, -0.034}$ }&
! 248: {\footnotesize $_{+0.034, -0.034}$ }\\
! 249: {\footnotesize Z Scale Factor }&
! 250: {\footnotesize $_{+0.072, -0.072}$ }&
! 251: {\footnotesize $_{+0.072, -0.072}$ }&
! 252: {\footnotesize $_{+0.048, -0.048}$ }\\
! 253: {\footnotesize Template statistics }&
! 254: {\footnotesize $_{+0.118, -0.118}$ }&
! 255: {\footnotesize $_{+0.170, -0.170}$ }&
! 256: {\footnotesize $_{+0.102, -0.102}$ }\\
! 257: {\footnotesize Signal contamination}&
! 258: {\footnotesize $_{+0.050, -0.052}$ }&
! 259: {\footnotesize $_{+0.136, -0.187}$ }&
! 260: {\footnotesize $_{+0.051, -0.051}$ }\\
! 261: {\footnotesize PDF }&
! 262: {\footnotesize $_{+0.097, -0.084}$ }&
! 263: {\footnotesize $_{+0.188, -0.198}$ }&
! 264: {\footnotesize $_{+0.092, -0.081}$ }\\
! 265:
! 266: \end{tabular}{\footnotesize \par}
! 267: %\end{ruledtabular}
! 268: \label{cap:Syst2}
! 269: \end{table}
! 270:
! 271:
! 272:
! 273: \clearpage
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>