Annotation of ttbar/p20_taujets_note/Strategy_and_Dataset.tex, revision 1.1
1.1 ! uid12904 1:
! 2: \section{Strategy}
! 3:
! 4:
! 5: \subsection{\label{sub:Signal-characteristics}Signal characteristics}
! 6:
! 7:
! 8: \subsubsection{Parton level MC}
! 9:
! 10: First of all we want to examine the properties of our signal. For
! 11: these purposes we used MC simulated samples (\textasciitilde{}10K
! 12: events each) generated with ALPGEN \cite{ALPGEN} interfaced to Pythia
! 13: \cite{PYTHIA} for showering and fragmentation.
! 14:
! 15:
! 16: \paragraph{$t\bar{t}$}
! 17:
! 18: First of all we wanted to look at the properties of the top quark
! 19: itself. Figure \ref{t-jets} shows the $\eta$, $\phi$ and $P_{T}$
! 20: distributions.
! 21:
! 22: %
! 23: \begin{figure}
! 24: \subfigure[$P_{T}$ of b-quarks]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/pt_tjets}}\subfigure[$\eta$ of b-quarks]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/eta_tjets}}
! 25:
! 26: \subfigure[$\phi$ of b-quarks]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/phi_tjets}}
! 27:
! 28:
! 29: \caption{Properties of t-quarks at the parton level in signal MC sample, $P_{T}$>
! 30: 15 GeV}
! 31:
! 32: \begin{centering}\label{t-jets}\par\end{centering}
! 33: \end{figure}
! 34:
! 35:
! 36:
! 37: \paragraph{$\tau$ and $\not\!\! E_{T}$}
! 38:
! 39: The unique property of $\tau$ (compared to other leptons) is that
! 40: it emits neutrino in its decay before it even reaches detector volume,
! 41: contributing to the missing energy of the event. For this reason the
! 42: $\tau$, available for measurement is not the same as physical $\tau$
! 43: produced. Figure \ref{cap:MC taus} demonstrates how a sizable fraction
! 44: of $\tau$ momentum goes missing. The plots of $\tau$ $\eta$ (Figure
! 45: \ref{cap:MC taus}) and transverse mass with $\not\!\! E_{T}$ thus
! 46: are done for the visible part of $\tau$
! 47:
! 48: This is to be compared with the $e+jets$ channel (Figure \ref{cap:MC elecs}).
! 49: As one can observe, the {}``total'' $\tau$ leptons behave very
! 50: simmilar to electrons, as of cause expected. However, after taking
! 51: into account the lost part of the $\tau$ energy situation becomes
! 52: very different.
! 53:
! 54: %
! 55: \begin{figure}
! 56: \subfigure[Parton level $\not E_{T}$ for the signal MC (red is total, green- visible)]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/met_mc}}
! 57: \subfigure[mT of MC $\tau$ and $\not E_{T} $ (red is total, green- visible)]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/mt_tau_met_mc}}
! 58:
! 59: \subfigure[$\eta$ of the MC $\tau$ in signal MC (red is total, green- visible)]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/eta_tau_mc}}
! 60: \subfigure[$P_{T}$ of the MC $\tau$ in signal MC (red is total, green- visible)]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/pt_tau_mc}}
! 61:
! 62: \subfigure[Visible energy fraction of $\tau$ in signal MC]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/visable_fraction}}
! 63:
! 64:
! 65: \caption{Properties of $\not\!\! E_{T}$ and $\tau$ at the parton level in
! 66: signal MC sample}
! 67:
! 68: \begin{centering}\label{cap:MC taus}\par\end{centering}
! 69: \end{figure}
! 70:
! 71:
! 72: %
! 73: \begin{figure}
! 74: \subfigure[Parton level $\not E_{T}$ for the signal MC (red is total, green- visible)]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/met_mc_elec}}
! 75: \subfigure[mT of MC electron and $\not E_{T} $ (red is total, green- visible)]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/mt_tau_met_mc_elec}}
! 76:
! 77: \subfigure[$\eta$ of the MC electron in signal MC (red is total, green- visible)]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/eta_tau_mc_elec}}
! 78: \subfigure[$P_{T}$ of the MC electron in signal MC (red is total, green- visible)]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/pt_tau_mc_elec}}
! 79:
! 80:
! 81: \caption{Properties of $\not\!\! E_{T}$ and electron at the parton level
! 82: in $t\bar{t}\rightarrow e+jets$ MC sample}
! 83:
! 84: \begin{centering}\label{cap:MC elecs}\par\end{centering}
! 85: \end{figure}
! 86:
! 87:
! 88:
! 89: \paragraph{Jets}
! 90:
! 91: b - jets are shown on Figure \ref{b-jets} while the product of W
! 92: decay are shown on Figure \ref{not b-jets}. This however doesn't
! 93: account for all the jets that will be reconstructed. Figure \ref{not b-jets-all}
! 94: demonstrates all the non-b quarks and gluons in a $t\bar{t}$ event.
! 95:
! 96: %
! 97: \begin{figure}
! 98: \subfigure[Number of b-quarks]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/nbquarks}}
! 99: \subfigure[$P_{T}$ of b-quarks]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/pt_bjets}}
! 100:
! 101: \subfigure[$\eta$ of b-quarks]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/eta_bjets}}
! 102: \subfigure[$\phi$ of b-quarks]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/phi_bjets}}
! 103:
! 104:
! 105: \caption{Properties of b-quarks at the parton level in signal MC sample, $P_{T}$>
! 106: 15 GeV}
! 107:
! 108: \begin{centering}\label{b-jets}\par\end{centering}
! 109: \end{figure}
! 110:
! 111:
! 112: %
! 113: \begin{figure}
! 114: \subfigure[Number of light quarks from the W decay]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/nlightquarks}}
! 115: \subfigure[$P_{T}$ of light quarks from the W decay]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/pt_notbjets}}
! 116:
! 117: \subfigure[$\eta$ of light quarks from the W decay]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/eta_notbjets}}
! 118: \subfigure[$\phi$ of light quarks from the W decay]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/phi_notbjets}}
! 119:
! 120:
! 121: \caption{Properties of light quarks from the W decay at the parton level in
! 122: signal MC sample, $P_{T}$> 15 GeV}
! 123:
! 124: \begin{centering}\label{not b-jets}\par\end{centering}
! 125: \end{figure}
! 126:
! 127:
! 128: %
! 129: \begin{figure}
! 130: \subfigure[Number of not b-quarks]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/nlightquarks_all}}
! 131: \subfigure[$P_{T}$ of not b-quarks]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/pt_notbjets_all}}
! 132:
! 133: \subfigure[$\eta$ of not b-quarks ]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/eta_notbjets_all}}
! 134: \subfigure[$\phi$ of not b-quarks ]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/phi_notbjets_all}}
! 135:
! 136:
! 137: \caption{Properties of all the light quarks and gluons in an event at the
! 138: parton level in signal MC sample, $P_{T}$> 15 GeV}
! 139:
! 140: \begin{centering}\label{not b-jets-all}\par\end{centering}
! 141: \end{figure}
! 142:
! 143:
! 144:
! 145: \subsubsection{Detector signature (reconstructed Monte Carlo)}
! 146:
! 147: Now we need to find out how well our experiment observes and reconstructs
! 148: this physical process. The $t\bar{t}\rightarrow\tau+jets$ Monte Carlo
! 149: file was processed through the detailed D0 Detector simulation. Jets
! 150: are reconstructed, using 0.5 radius cone (in $\eta-\phi$). Taus are
! 151: identified using tau ID algorithm and we apply 0.8 cut on the $\tau$
! 152: selection neural net (section \ref{sub:tau--ID}).
! 153:
! 154:
! 155: \paragraph{$\tau$ and $\not\!\! E_{T}$}
! 156:
! 157: We can see that for reconstructed $\tau$'s $mT(of$ $\tau$ $and$
! 158: $\not E_{T})$ (Figure \ref{cap:reco tau}) doesn't look as good as
! 159: for MC taus (Figure \ref{cap:MC taus}). We observe a noticeable \char`\"{}tail\char`\"{}
! 160: above 80 GeV.
! 161:
! 162: %
! 163: \begin{figure}
! 164: \subfigure[Reconstructed $\not E_{T}$ for the signal MC ]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/met_reco}}
! 165: \subfigure[mT of MC $\tau$ and $\not E_{T} $ reconstructed ]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/mt_tau_met_reco}}
! 166:
! 167: \subfigure[$\eta$ of the MC tau in signal MC reconstructed ]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/eta_tau_reco}}
! 168: \subfigure[ $P_{T}$ of the $\tau$ ]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/pt_tau_reco}}
! 169:
! 170:
! 171: \caption{Properties of detector reconstructed $\not\!\! E_{T}$ and $\tau$
! 172: in signal MC sample}
! 173:
! 174: \begin{centering}\label{cap:reco tau}\par\end{centering}
! 175: \end{figure}
! 176:
! 177:
! 178:
! 179: \paragraph{Jets}
! 180:
! 181: Before b-tagging (section \ref{sub:B-tagging}) one can't separate
! 182: b-jets from non-b jets, so we don't make a distinction at this point.
! 183: Most important variables are the number of jets and $\eta$ and $P_{T}$
! 184: distributions (Figure \ref{jets1}).
! 185:
! 186: %
! 187: \begin{figure}
! 188: \subfigure[Number of jets]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/njets}}
! 189: \subfigure[$P_{T}$ of jets]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/jetpt}}
! 190:
! 191: \subfigure[$\eta$ of jets]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/eta_jets}}
! 192: \subfigure[$\phi$ of jets]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/phi_jets}}
! 193:
! 194:
! 195: \caption{Properties of the jets in signal MC sample, $P_{T}$> 15 GeV}
! 196:
! 197: \begin{centering}\label{jets1}\par\end{centering}
! 198: \end{figure}
! 199:
! 200:
! 201: Jets are arranged in the order of their $P_{T}$: leading (highest),
! 202: sub-leading etc. We can see on Figure \ref{jets1} that we typically
! 203: have 4 or 5 jets in an event. It is interesting to compare the leading
! 204: jet to the fourth and fifth jets (Figure \ref{jets2}). We can see
! 205: that jets after the third are very soft, as expected.
! 206:
! 207: %
! 208: \begin{figure}
! 209: \subfigure[1st jet]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/jet0pt}}
! 210: \subfigure[2nd jet]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/jet1pt}}
! 211:
! 212: \subfigure[3rd jet]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/jet2pt}}
! 213: \subfigure[4th jet]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/jet3pt}}
! 214:
! 215: \subfigure[5th jet]{\includegraphics[scale=0.3]{MS_thesis/proposal_plots/jet4pt}}
! 216:
! 217:
! 218: \caption{$P_{T}$ distributions for jets in signal MC sample (including $\tau$)}
! 219:
! 220: \begin{centering}\label{jets2}\par\end{centering}
! 221: \end{figure}
! 222:
! 223:
! 224:
! 225: \subsection{Backgrounds}
! 226:
! 227: Two main distinctive features of the signal limit the spectrum of
! 228: important backgrounds. In order to be relevant the process must have
! 229: high (>3) number of jets as well as sizable (>15 GeV) $\not\!\! E_{T}$.
! 230: All the candidate processes are listed in Table \ref{backgrounds}.
! 231: The cross section listed include the branching fractions into $\tau$
! 232:
! 233: We can conclude that two dominant background sources are QCD ({}``fake
! 234: $\tau$'') and W+4jets. These two sources were taken into account
! 235: in these analysis.
! 236:
! 237: %
! 238: \begin{table}
! 239: \begin{tabular}{|l|l|c|}
! 240: \hline
! 241: \multicolumn{1}{||l|}{Background}&
! 242: \multicolumn{1}{||l|}{Description}&
! 243: \multicolumn{1}{||l||}{Cross Section}\tabularnewline
! 244: \hline
! 245: $W+jjjj\rightarrow\tau\nu jjjj$&
! 246: Has identical signature to the signal&
! 247: $\sim$18 pb \tabularnewline
! 248: \hline
! 249: $Z/\gamma+jjj\rightarrow\tau\tau jjj$ &
! 250: $\tau$ is usually found as a jet&
! 251: $\sim$2.6 pb\tabularnewline
! 252: \hline
! 253: $WZ\rightarrow\tau\nu jj$&
! 254: needs two extra jet (can be gluon emission) &
! 255: $\sim$0.2 pb\tabularnewline
! 256: \hline
! 257: $WW\rightarrow\tau\nu jj$&
! 258: needs two extra jet (can be gluon emission) &
! 259: $\sim$0.5 pb\tabularnewline
! 260: \hline
! 261: single top &
! 262: small cross section, but has b-jets &
! 263: $\sim$0.5 pb\tabularnewline
! 264: \hline
! 265: QCD &
! 266: Any 4-jet event, that doesn't have a real $\tau$ in it&
! 267: $>$100 nb \tabularnewline
! 268: \hline
! 269: \end{tabular}\centering
! 270:
! 271:
! 272: \caption{Background sources, relevant for the $\tau+jets$ analysis. Branching
! 273: into hadronic $\tau$ had been applied \cite{l+jets}}
! 274:
! 275: \label{backgrounds}
! 276: \end{table}
! 277:
! 278:
! 279:
! 280: \section{Dataset}
! 281:
! 282: As will be demonstrated in section \ref{sub:Running-trigsim} the
! 283: optimal (that is most efficient) combination of triggers for this
! 284: analysis is:
! 285:
! 286: \begin{itemize}
! 287: \item The Higgs Missing $H_{T}$ trigger (MHT30\_3CJT5)
! 288: \item The ALLJET trigger (4JT10)
! 289: \end{itemize}
! 290: Together they yield over 85\% signal acceptance and they'd been running
! 291: unprescaled for most of D0 stable operation.
! 292:
! 293: The data skim, utilizing both of these triggers would be optimal for
! 294: this analysis. Until the technical issues involved in its production
! 295: are fully resolved, we are using the ALLJET skim \cite{Luminosity},
! 296: which only contains the data, collected by the 4JT10 (and its subsequent
! 297: versions).
! 298:
! 299: Such skim is only 70\% efficient for the signal, but it's the closest
! 300: available for our needs at the moment. The full PASS2 ALLJET skim
! 301: had been processed through the standard D0 top group data quality
! 302: criteria, discarding bad luminosity blocks, at the same time computing
! 303: the recorded lumi. The results are represented in table \ref{lumi1}.
! 304: Therefore the total luminosity available for the analysis amounts
! 305: to $349\pm23$ $pb^{-1}$ \cite{alljet}
! 306:
! 307: %
! 308: \begin{table}
! 309: \begin{tabular}{|c|c|c|c|}
! 310: \hline
! 311: Stage&
! 312: Luminosity ($pb^{-1})$&
! 313: Relative Size (\%)&
! 314: Absolute Size (\%)\tabularnewline
! 315: \hline
! 316: \hline
! 317: Delivered&
! 318: 482.6&
! 319: 100&
! 320: 100\tabularnewline
! 321: \hline
! 322: Recorded&
! 323: 411.6&
! 324: 85.3&
! 325: 85.3\tabularnewline
! 326: \hline
! 327: Good&
! 328: 352.5&
! 329: 85.6&
! 330: 73.0\tabularnewline
! 331: \hline
! 332: Reconstructed&
! 333: 349.3&
! 334: 99.1&
! 335: 72.4\tabularnewline
! 336: \hline
! 337: \end{tabular}
! 338:
! 339:
! 340: \caption{The results of luminosity calculation for the PASS2 ALLJET top skim}
! 341:
! 342: \label{lumi1}
! 343: \end{table}
! 344:
! 345:
! 346: The table \ref{lumi2} demonstrates the breakdown of this luminosity
! 347: between the different trigger versions.
! 348:
! 349: %
! 350: \begin{table}
! 351: \begin{tabular}{|c|c|c|}
! 352: \hline
! 353: Trigger version&
! 354: Trigger name&
! 355: Luminosity ($pb^{-1})$\tabularnewline
! 356: \hline
! 357: \hline
! 358: 8.0&
! 359: 4JT10&
! 360: 19.4\tabularnewline
! 361: \hline
! 362: 9.0&
! 363: 4JT10&
! 364: 21.2\tabularnewline
! 365: \hline
! 366: 10.0&
! 367: 4JT10&
! 368: 15.1\tabularnewline
! 369: \hline
! 370: 11.0&
! 371: 4JT10&
! 372: 57.3\tabularnewline
! 373: \hline
! 374: 12.0&
! 375: 4JT12&
! 376: 196\tabularnewline
! 377: \hline
! 378: 13.0&
! 379: JT2\_4JT12L\_HT&
! 380: 13.5\tabularnewline
! 381: \hline
! 382: 13.1&
! 383: JT2\_4JT12L\_HT&
! 384: 27.8\tabularnewline
! 385: \hline
! 386: 13.3&
! 387: JT2\_4JT12L\_HT&
! 388: 0\tabularnewline
! 389: \hline
! 390: \end{tabular}
! 391:
! 392:
! 393: \caption{Luminosity of the ALLJET skim for different D0 trigger list versions}
! 394:
! 395: \label{lumi2}
! 396: \end{table}
! 397:
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>