Annotation of ttbar/p20_taujets_note/Tools.tex, revision 1.1.1.1
1.1 uid12904 1:
2: \section{Tools}
3:
4:
5: \subsection{Object ID}
6:
7:
8: \subsubsection{\label{sub:tau--ID}$\tau$ ID}
9:
10:
11: \paragraph{Tau decay modes}
12:
13: The $\tau$ lepton have several decay channels, classified by the
14: number of charged particles (tracks) associated with it \cite{PDG}
15: :
16:
17: \begin{itemize}
18: \item electron + muon ($\tau\rightarrow e\nu_{e}\nu_{\tau}$ or $\tau\rightarrow\mu\nu_{\mu}\nu_{\tau})$,
19: BR = 35\%
20: \item charged hadron ($\tau\rightarrow\pi^{-}\nu_{\tau}$), BR =12\%
21: \item charged hadron + $\geq1$ neutral particle (i.e. $\tau\rightarrow\rho^{-}\nu_{\tau}\rightarrow\pi^{0}n+\pi^{-}\nu_{\tau}$)
22: , BR = 38\%
23: \item 3 charged hadrons + $\geq0$ neutral hadrons, BR = 15\% (so-called
24: {}``3-prong'' decays)
25: \end{itemize}
26:
27: \paragraph{Tau ID variables}
28:
29: At D0 $\tau$s are identified in their hadronic modes (contributing
30: to inefficiency of id) as narrow (0.3 cone) jets,isolated and matched
31: to a charged track. The (most important) discriminating variables
32: are \cite{tau ID}:
33:
34: \begin{itemize}
35: \item Profile - $\frac{E_{T}^{1}+E_{T}^{0}}{\sum_{i}E_{T}^{i}}$, where
36: $E_{T}^{i}$ is the $E_{T}$ of the ith highest $E_{T}$ tower in
37: the cluster
38: \item Isolation, defined as $\frac{E(0.5)-E(0.3)}{E(0.3)}$, where $E(R)$
39: is the energy contained in a radius of R around cal cluster centroid
40: \item Track isolation, defined as $\sum p_{T}$ of non-$\tau$ tracks in
41: cone of 0.5 around the calorimeter cluster centroid
42: \end{itemize}
43: Using these and other variables, 2 Neural Networks are trained to
44: identify 3 types of $\tau$ ($\pi$-type, $\rho$-type and 3-prong)
45:
46: The output of these NN provides a set of 3 variables (nnout 1... 3)
47: to be used to select $\tau$ in the event. The high values of NN have
48: to correspond to the physical $\tau$ leptons, while the low ones
49: should indicate fakes. For more details, see \cite{tau ID}.
50:
51:
52: \paragraph{Energy Scale}
53:
54: In \cite{ingo1} the process $Z\rightarrow\tau\tau$ had been studied.
55: In particular the Figure 4 of that note demonstrates an excellent
56: agreement between the data and $Z\rightarrow\tau\tau$ MC in distribution
57: of the invariant mass of the $\tau$ pair. Figure 15 of \cite{ingo2}
58: shows other important properties ($P_{T}$ of $\tau$ and $\not E_{T}$)
59: which also agree very well. Since no energy correction had been applied
60: to the $\tau$ in this work, one can conclude that we can take the
61: energy scale of $\tau$ ID to be 1.
62:
63:
64: \paragraph{Performance }
65:
66: The $\tau$ NN had been trained and optimized for the low jet multiplicity
67: events (i.e. $Z\rightarrow\tau\tau$). We wanted to compare its performance
68: for the high multiplicity signal (top) that we are searching for here.
69:
70: In order to evaluate the ID efficiency reliably one has to match the
71: reconstructed $\tau$ candidate with the true $\tau$ from MC. We
72: start with all the $\tau$ candidates in an event, regardless of the
73: quality. We then want to select those that can be with high confidence
74: correspond to the real $\tau$ leptons. The assumption is that those
75: $\tau$ candidates, whose energy and direction matched to a physical
76: $\tau$ are indeed representing the detector signature of this particle.
77: We can then determine how well does the $\tau$ ID identify this $\tau$
78: lepton. Figure \ref{cap:Matching-of-MC} illustrates the matching
79: procedure - the $\tau$ candidates with $\Delta R$ from a real MC
80: $\tau$ of 0.05 and $\Delta P$ of 10 GeV are deemed to be the {}``real''
81: matched $\tau$
82:
83: For such $\tau$ we plotted the NN for different $\tau$ types (Fig
84: \ref{cap:NN-for-matched}). From these one can determine the efficiency
85: of $\tau$ ID for various cuts on NN (Fig \ref{tauID}).
86:
87: %
88: \begin{figure}
89: \subfigure[$\Delta R$ between reco $\tau$ and MC $\tau$]{\includegraphics[scale=0.3]{plots_for_talk/drmin}}\subfigure[$\Delta R$ between reco $\tau$ and MC $\tau$ (low values)]{\includegraphics[scale=0.3]{plots_for_talk/drmin_zoomed}}
90:
91: \subfigure[Difference in energy between reco and MC $\tau$]{\includegraphics[scale=0.3]{plots_for_talk/dpmin}}\subfigure[Difference in energy between MC and reco $\tau$ that were matched in angle]{\includegraphics[scale=0.3]{plots_for_talk/dpmin005}}
92:
93:
94: \caption{Matching of MC $\tau$ and reco $\tau$. Black is $Z\rightarrow\tau\tau$,
95: red is $t\overline{t}\rightarrow\tau+jets$The histograms are normalized
96: to 1 to enable comparision.}
97:
98: \label{cap:Matching-of-MC}
99: \end{figure}
100:
101:
102: %
103: \begin{figure}
104: \subfigure[NN for ALL types]{\includegraphics[scale=0.3]{plots_for_talk/nnmatched}}\subfigure[NN for type 1]{\includegraphics[scale=0.3]{plots_for_talk/nn1matched}}
105:
106: \subfigure[NN type 2]{\includegraphics[scale=0.3]{plots_for_talk/nn2matched}}\subfigure[NN for type 3]{\includegraphics[scale=0.3]{plots_for_talk/nn3matched}}
107:
108:
109: \caption{NN for matched $\tau$. Black is $Z\rightarrow\tau\tau$, red is
110: $t\overline{t}\rightarrow\tau+jets$. The histograms are normalized
111: to 1 to enable comparision.}
112:
113: \label{cap:NN-for-matched}
114: \end{figure}
115:
116:
117: %
118: \begin{figure}
119: \subfigure[ALL types]{\includegraphics[scale=0.3]{plots_for_talk/eff0}}\subfigure[Type 1]{\includegraphics[scale=0.3]{plots_for_talk/eff1}}
120:
121: \subfigure[Type 2]{\includegraphics[scale=0.3]{plots_for_talk/eff2}}\subfigure[Type 3]{\includegraphics[scale=0.3]{plots_for_talk/eff3}}
122:
123:
124: \caption{$\tau$ ID Efficiencies for different types. Black is $Z\rightarrow\tau\tau$,
125: red is $t\overline{t}\rightarrow\tau+jets$}
126:
127: \label{tauID}
128: \end{figure}
129:
130:
131: In order to choose the best cut on $\tau$ NN one has to also consider
132: the fake rate (the number of fake $\tau$ candidates passing the ID
133: requirements successfully). For this purpose we had examined the $\tau$
134: candidates in the preselected ALLJET data sample (the details of preselection
135: are described in section \ref{sub:Preselection}). Since this dataset
136: is QCD dominated (no more then 0.2\% of electroweak is expected) we
137: can safely assume all $\tau$ in it to be fake (this assumption will
138: be employed again for our QCD background estimation in section \ref{sub:QCD-modeling}).
139: Figure \ref{cap:NN-for-fake} shows the distribution of NN for the
140: $\tau$. From this we can determine the fake rate dependence on NN
141: cut (Fig \ref{tauID_Fake}). We can note that type 3 has noticeably
142: higher fake rate. This is to be expected, since most jets have higher
143: track multiplicities than type 1 and 2 $\tau$ making it harder for
144: them to pass $\tau$ ID requirements.
145:
146: %
147: \begin{figure}
148: \subfigure[NN for ALL types]{\includegraphics[scale=0.3]{plots/NN0fakeNN}}\subfigure[NN for type 1]{\includegraphics[scale=0.3]{plots/NN1akeNN}}
149:
150: \subfigure[NN type 2]{\includegraphics[scale=0.3]{plots/NN2akeNN}}\subfigure[NN for type 3]{\includegraphics[scale=0.3]{plots/NN3akeNN}}
151:
152:
153: \caption{NN for fake $\tau$}
154:
155: \label{cap:NN-for-fake}
156: \end{figure}
157:
158:
159: %
160: \begin{figure}
161: \subfigure[ALL types]{\includegraphics[scale=0.3]{plots/NN0fake}}\subfigure[Type 1]{\includegraphics[scale=0.3]{plots/NN1fake}}
162:
163: \subfigure[Type 2]{\includegraphics[scale=0.3]{plots/NN2fake}}\subfigure[Type 3]{\includegraphics[scale=0.3]{plots/NN3fake}}
164:
165:
166: \caption{$\tau$ ID Fake Rate for different types}
167:
168: \label{tauID_Fake}
169: \end{figure}
170:
171:
172: On Fig \ref{tauID_Fake_Eff} we plot the fake rate vs. efficiency
173: of the $\tau$ ID for our channel. From this we can select the optimal
174: selection cut on $\tau$ NN, based on the $\tau$ID significance,
175: defined as $\frac{Number\, of\, real\, taus}{\sqrt{Number\, of\, real+Number\, fakes}}$
176: (Fig \ref{tauID_Fake_signif}). It is computed on our preselected
177: analysis data set (section \ref{sub:Preselection})
178:
179: We can conclude that D0 $\tau$ ID algorithm has efficiency for $t\bar{t}$
180: comparable with $Z\rightarrow\tau\tau$. The optimal cut on $\tau$
181: NN appears to be 0.95 for all the types.
182:
183: %
184: \begin{figure}
185: \subfigure[ALL types]{\includegraphics[scale=0.3]{plots/NN0fake_eff}}\subfigure[Type 1]{\includegraphics[scale=0.3]{plots/NN1fake_eff}}
186:
187: \subfigure[Type 2]{\includegraphics[scale=0.3]{plots/NN2fake_eff}}\subfigure[Type 3]{\includegraphics[scale=0.3]{plots/NN3fake_eff}}
188:
189:
190: \caption{$\tau$ ID Efficiency vs. the Fake rate. Type 2 is the cleanest,
191: type 3 has highest fake rate, as expected}
192:
193: \label{tauID_Fake_Eff}
194: \end{figure}
195:
196:
197: %
198: \begin{figure}
199: \subfigure[ALL types]{\includegraphics[scale=0.3]{plots/NN0fake_signif}}\subfigure[Type 1]{\includegraphics[scale=0.3]{plots/NN1fake_signif}}
200:
201: \subfigure[Type 2]{\includegraphics[scale=0.3]{plots/NN2fake_signif}}\subfigure[Type 3]{\includegraphics[scale=0.3]{plots/NN3fake_signif}}
202:
203:
204: \caption{$\tau$ ID Significance vs. the NN cut. The 0.95 cut appears to be
205: advantageous for all the types}
206:
207: \label{tauID_Fake_signif}
208: \end{figure}
209:
210:
211:
212: \subsubsection{\label{sub:B-tagging}B-tagging}
213:
214: The chosen b-tagging algorithm is Secondary Vertex Tagger (SVT) \cite{b-ID}.
215: It is characterized by high (compared to other taggers) purity, which
216: is essential for such QCD-dominated channel as ours.
217:
218: The algorithm reconstructs secondary vertices inside a jet, using
219: the jet's associated tracks. The tracks are also required to pass
220: a set of cuts outlined in Table \ref{cap:The-standard-cuts}. Then,
221: the decay length significance is computed. If the jet has this significance
222: grater then 7 (for SVT TIGHT) it is considered b-tagged.
223:
224: As can be seen from Figures \ref{cap:SVT-Signifficance} and \ref{cap:SVT-Signifficance_data},
225: the TIGHT cut is most appropriate for our signal.
226:
227: %
228: \begin{table}
229: \begin{tabular}{|l|l|l|l|}
230: \hline
231: {\large SVT}&
232: &
233: &
234: \tabularnewline
235: \hline
236: &
237: LOOSE &
238: MEDIUM &
239: TIGHT\tabularnewline
240: \hline
241: Number of SMT hits &
242: 2 &
243: 2 &
244: 2\tabularnewline
245: \hline
246: $P_{T}$ of tracks &
247: 1 GeV/c &
248: 1 GeV/c &
249: 1 GeV/c\tabularnewline
250: \hline
251: impact parameter significance of tracks &
252: 3 &
253: 3.5 &
254: 3.5\tabularnewline
255: \hline
256: track $\chi^{2}$&
257: 10 &
258: 10 &
259: 10\tabularnewline
260: \hline
261: max vertex $\chi^{2}$&
262: 100 &
263: 100 &
264: 100\tabularnewline
265: \hline
266: vertex collinearity &
267: 0.9 &
268: 0.9 &
269: 0.9\tabularnewline
270: \hline
271: max vertex decay length &
272: 2.6cm &
273: 2.6cm &
274: 2.6cm\tabularnewline
275: \hline
276: Decay Length Significance Cut &
277: 5 &
278: 6 &
279: 7\tabularnewline
280: \hline
281: \end{tabular}
282:
283:
284: \caption{The standard cuts on SVT \cite{b-ID}}
285:
286: \label{cap:The-standard-cuts}
287: \end{table}
288:
289:
290: %
291: \begin{figure}
292: \includegraphics[scale=0.5]{MS_thesis/proposal_plots/svt}
293:
294:
295: \caption{SVT Decay Length Significance for the b-jets in $t\overline{t}\rightarrow\mu+jets$
296: MC}
297:
298: \label{cap:SVT-Signifficance}
299: \end{figure}
300:
301:
302: %
303: \begin{figure}
304: \subfigure[SVT significance across the entire range]{\includegraphics[scale=0.3]{plots/svt_data}}\subfigure[SVT significance at values near 0]{\includegraphics[scale=0.3]{plots/svt_data_zoomed}}
305:
306:
307: \caption{SVT Decay Length Significance for all the jets in the ALLJET skim
308: data.}
309:
310: \label{cap:SVT-Signifficance_data}
311: \end{figure}
312:
313:
314:
315: \paragraph{Taggability}
316:
317: In order to reconstruct a secondary vertex in a jet, the jet must
318: contain at least 2 tracks. If such tracks are found and their $P_{T}$
319: is greater then 0.5 GeV the jet is called taggable. In MC it is important
320: to distinguish the taggability from the tagging efficiency, since
321: the later depends on the jet's flavor.
322:
323:
324: \paragraph{B-tagging efficiency}
325:
326: It is known that b-tagging applied directly to MC gives an overestimated
327: efficiency. In order to account for this factor SVT had been parameterized
328: on $t\overline{t}\rightarrow\mu+jets$ MC and $\mu+jets$ data to
329: compute the correction factor, which has to be applied to MC. As result
330: we obtain the MC tagging probability and data corrected one (Figure
331: \ref{bID}). It can be noted that the data corrected efficiency is
332: indeed noticeably (>30\%) lower then what we would expect by applying
333: SVT directly to MC.
334:
335:
336: \paragraph{C-tagging efficiency}
337:
338: An assumption is made that the correction factor obtained by dividing
339: the semi-leptonic b-tagging efficiency in data to the one in MC also
340: is correct for c-jets. Hence the MC-obtained inclusive c-taging efficiency
341: is multiplied by this factor (and by it's taggability too) in order
342: to estimate the c-tagging probability
343:
344:
345: \paragraph{Light jet tagging efficiency}
346:
347: The b-tag fake rate from light quarks is computed by measuring the
348: negative tag rate. It is defined by the rate of appearance of secondary
349: vertices with negative decay length significance. It is assumed that
350: the light quarks have equal chances to produce SV with positive and
351: negative decay length significance (due to finite resolution effects)
352: while the heavy flavor jets can only produce SV with positive decay
353: length significance. This however is not quite true and a special
354: scaling factor ($SF_{hf}$) is introduced to correct for the fraction
355: of heavy flavors among the jets with the negative decay length significance.
356: Another correction is for the presence of the long lived particles
357: in light jets ($SF_{ll}$). Both factors are derived from Monte Carlo.
358:
359: %
360: \begin{figure}
361: \subfigure[SVT efficiency]{\includegraphics[scale=0.3]{plots/SVTeff_pt}}\subfigure[SVT efficiency]{\includegraphics[scale=0.3]{plots/SVTeff_eta}}
362:
363: \subfigure[SVT efficiency parametrized on data]{\includegraphics[scale=0.3]{plots/SVTeff}}\subfigure[SVT efficiency parametrized on MC]{\includegraphics[scale=0.3]{plots/SVTeffMC}}
364:
365:
366: \caption{SVT Efficiency for $t\overline{t}\rightarrow\tau+jets$ MC. Red is
367: MC parameterization black is data-corrected. Flavor depancance is
368: taken into account. The lower plots show 2D parametrizations}
369:
370: \label{bID}
371: \end{figure}
372:
373:
374:
375: \paragraph{Event tagging efficiency}
376:
377: The tag rates and the taggability had been combined and used to predict
378: the probability for a jet to be b-tagged (b-tagging weight). The final
379: resulting per-event probability of having at least one such a tag
380: for the $t\overline{t}\rightarrow\tau+jets$ MC is plotted on Figure
381: \ref{cap:The-probability-to}
382:
383: %
384: \begin{figure}
385: \includegraphics[scale=0.5]{plots/ttb_eventprobtag}
386:
387:
388: \caption{The probability to tag at least one jet with SVT for $t\overline{t}\rightarrow\tau+jets$
389: MC}
390:
391: \label{cap:The-probability-to}
392: \end{figure}
393:
394:
395: Finally it has to be noted that we tried to avoid th overlap between
396: $\tau$ ID and b-tagging. That is we remove the jets, matched to a
397: 0.8 $\tau$ candidate within $\Delta R$ = 0.5
398:
399:
400: \subsection{Trigger}
401:
402:
403: \subsubsection{\label{sub:Running-trigsim}Running TRIGSIM}
404:
405: In order to search for a signal one has to collect the data with a
406: chosen set of triggers. We want to find such a combination so to maximize
407: the fraction of signal events written out.
408:
409: For that purpose the trigger simulation program was run on the MC
410: signal. The following efficiencies were obtained (for 12.30 version
411: of the global D0 trigger definition list) (Table \ref{cap:Event-overlaps-between})
412:
413: %
414: \begin{table}
415: \begin{tabular}{|l||l||}
416: \hline
417: {\large Trigger}&
418: {\large Fraction of events passing}\tabularnewline
419: \hline
420: 4JT12 &
421: 0.74$\pm$0.05\tabularnewline
422: 3J15\_2J25\_PVZ &
423: 0.73$\pm$0.05\tabularnewline
424: MHT30\_3CJT5&
425: 0.68$\pm$0.04\tabularnewline
426: MU\_JT20\_L2M0&
427: 0.30$\pm$0.01\tabularnewline
428: \hline
429: \multicolumn{1}{|l|}{MU\_JT20\_L2M0 \&\& MHT30\_3CJT5}&
430: \multicolumn{1}{||l||}{0.2$\pm$0.01}\tabularnewline
431: \hline
432: \multicolumn{1}{||l}{MHT30\_3CJT5 \&\& 4JT12}&
433: \multicolumn{1}{||l||}{0.4$\pm$0.04 }\tabularnewline
434: \hline
435: \multicolumn{1}{||l}{4JT12 \&\& 3J15\_2J25\_PVZ}&
436: \multicolumn{1}{||l||}{0.67$\pm$0.04}\tabularnewline
437: \hline
438: \end{tabular}
439:
440:
441: \caption{Trigger efficiencies and event overlaps between the most efficient
442: (for selecting $t\overline{t}\rightarrow\tau+jets$) unprescaled triggers.
443: As one can see 3J15\_2J25\_PVZ has large overlap with 4JT12, and since
444: 4JT12 is better studied it has been chosen for this analysis.}
445:
446: \label{cap:Event-overlaps-between}
447: \end{table}
448:
449:
450: Taking this into account, we are left with 3 triggers, giving altogether
451: \textasciitilde{}85\% efficiency:
452:
453: $\mathit{MHT30\_3CJT5}$ - $\not\!\! E_{T}$ trigger, requiring at
454: least 30 GeV at level 3, which leads to \textasciitilde{}30\% inefficiency,
455: since our missing $E_{T}$ peaks around 50 GeV
456:
457: $Description$: {\large L1}: At least three Calorimeter trigger towers
458: with $E_{T}$$>$5 GeV. {\large L2}: Require jet $E_{T}$$>$20. {\large L3}:
459: Vector $H_{T}$ sum $>$30 GeV. Also, one in 4000 events is recorded
460: and marked as {}``unbiased''
461:
462: $\mathit{4JT12}$ - Trigger designed for the $t\bar{t}\rightarrow jets$
463: analysis \cite{alljet}. Meets all the jet number requirements, but
464: doesn't often have high enough $\not E_{T}$
465:
466: $Description$: {\large L1}: At least three Calorimeter JET trigger
467: towers having $E_{T}$$>$5 GeV. {\large L2}: Three JET candidates
468: with $E_{T}$$>$8 GeV and HT $>$ 50 GeV. {\large L3}: Four $|$$\eta$$|$$<$3.6
469: jet candidates with $E_{T}$$>$10 GeV found using a simple cone algorithm.
470: Three of those jets must have $E_{T}$$>$15 GeV. Record one in 500
471: events marked as 'unbiased'.
472:
473: $\mathit{MU\_JT20\_L2M0}\mathbb{\,\,}$- Muon trigger with $<$20\%
474: efficiency for our signal, but it's unprescaled and has little overlap
475: with others
476:
477: $Description$: {\large L1}: A single muon trigger based on muon scintillator
478: and also requiring one Calorimeter JET trigger tower with $E_{T}$>3
479: GeV. {\large L2}: At least one muon found meeting MEDIUM quality requirements
480: but no pT or region requirement. Also require at least one jet with
481: $E_{T}$>10 GeV. {\large L3}: At least one jet with $E_{T}$>20 GeV
482: is found using a simple cone algorithm. Additionally, one in 500 of
483: all events is recorded and marked as 'unbiased'.
484:
485:
486: \subsubsection{Triggers in version 13}
487:
488: A new trigger list has recently been used for D0 data-taking. It includes
489: a number of new triggers and modifications to existing ones. Running
490: TRIGSIM on this triglist we discover that the triggers that were best
491: in v12 ($\mathit{4JT12}$ and $\mathit{MHT30\_3CJT5}$) are also most
492: efficient in v13. In fact, OR of just these two triggers gives 90$\pm$5\%
493: efficiency. The names and definitions of these triggers had changed:
494:
495: 4JT12 became JT2\_4JT12L\_HT. An additional $H_{T}$ cut of 120 GeV
496: is being applied:
497:
498: $Description$: {\large L1}: Three calorimeter trigger towers with
499: $E_{T}$>5 GeV. {\large L2}: Pass events with at least three JET candidates
500: with $E_{T}$>6 GeV and $H_{T}$, formed with jets above 6 GeV, greater
501: than 70 GeV. {\large L3}: Requires at least four jets with $E_{T}$
502: > 12 GeV and at at least three jets with $E_{T}$ > 15 GeV. Also require
503: at least two jets to have $E_{T}$ > 25 GeV . Event $H_{T}$ (calculated
504: using all jets with $E_{T}$>9 GeV) > 120 GeV .
505:
506: MHT30\_3CJT5 became JT2\_MHT25\_HT.
507:
508: $Description$:: {\large L1}: Three calorimeter trigger towers with
509: $E_{T}$>4 GeV, |$\eta$|<2.4, and two calorimeter trigger towers
510: with $E_{T}$>5 GeV. {\large L2}: Pass events with at least three
511: JET candidates with $E_{T}$>6 GeV and $H_{T}$, formed with jets
512: above 6 GeV, greater than 70 GeV {\large L3}: Vector $H_{T}$ sum
513: for the event must be above 25 GeV. Also require event scalar $H_{T}$
514: (calculated using all jets with $E_{T}$>9 GeV) > 125 GeV.
515:
516:
517: \subsubsection{Turn-on curves}
518:
519: These features are reflected in the corresponding turn-on curves (integrated
520: efficiencies of signal) . Such a curve for MHT30 and 4JT12 triggers
521: is shown on Figure \ref{trig}.%
522: \begin{figure}
523: \includegraphics[scale=0.8]{MS_thesis/proposal_plots/MET_eff_new__7}
524:
525: \includegraphics[scale=0.6]{MS_thesis/proposal_plots/4JT12jet3pt}
526:
527:
528: \caption{Integrated trigger efficiency for MHT30 and 4JT12 triggers for the
529: $t\overline{t}\rightarrow\tau+jets$ MC, obtained using TRIGSIM.}
530:
531: \begin{centering}\label{trig}\par\end{centering}
532: \end{figure}
533:
534:
535:
536: \subsubsection{Trigger simulation}
537:
538: At the time of the analysis TRIGSIM had not reached the state in which
539: it could reliably reproduce the trigger efficiency on data. Therefore,
540: the accepted practice is to parameterize the trigger turn-ons on data
541: and apply this parameterization to MC files.
542:
543: Such procedure was performed with the top\_trigger package \cite{top_trigger}.
544: On Figure \ref{cap:The-trigger-efficiency} one can see the results
545: of a test to check the validity of such approach. We used the dataset
546: collected by a single muon trigger (MU\_JT20\_L2M0). We can assume
547: that such data has little bias with respect to the 4JT10 trigger.
548: Hence, if we count the number of the 4 jet events that passed 4JT10
549: and compare it with the top\_trigger prediction for the same events
550: we can check how well does top\_trigger perform. As can be noted from
551: Figure \ref{cap:The-trigger-efficiency} the agreement is fairly good,
552: especially in the region which we use in this analysis (we require
553: jets to have $P_{T}>20$ GeV). The efficiency turn-on curve, produced
554: by top\_trigger is shown on Fig. \ref{cap:The-trigger-efficiency_MC}
555: and is in agreement with the TRIGSIM (Fig. \ref{trig}).
556:
557: %
558: \begin{figure}
559: \includegraphics[scale=0.7]{plots/4JT10_MULOOSE}
560:
561:
562: \caption{The integrated trigger efficiency closure plot. Black is the MU\_JT20\_L2M0
563: data, red is top\_trigger prediction for this data.}
564:
565: \label{cap:The-trigger-efficiency}
566: \end{figure}
567:
568:
569: %
570: \begin{figure}
571: \includegraphics[scale=0.7]{plots/4JT12_EFF_toptrig}
572:
573:
574: \caption{The integrated trigger efficiency for the $t\overline{t}\rightarrow\tau+jets$
575: MC obtained using top\_trigger. }
576:
577: \label{cap:The-trigger-efficiency_MC}
578: \end{figure}
579:
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>