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>