-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathMasterArbeit.tex
2283 lines (2282 loc) · 255 KB
/
MasterArbeit.tex
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
% --- TO DO ---
%FCNC im LQ Kaptiel braucht sehr wahrscheinlich eine Erklärung!! + schwacher isospin
%
%
%QCD und QED sollte vielleicht mal im SM erwähnt werden
%
% overlap removal bei phys obj sel + trigger erklären! https://arxiv.org/pdf/1611.09661.pdf https://cds.cern.ch/record/2133909/files/ATL-DAQ-PROC-2016-003.pdf (ntuple/settings/*)
%
%-----------------------------------------------------------------------------------------------
% Gliederungsideen
%Experiment-Part: LHC, ATLAS, LQpairprod @pp collisions
%
%Software-Part: Datenprozessierung, MC Simulation, b-tagging, anti-kT(jet reconstruction), tau reconstruction, e+µ detektion? (Unterpunkt von Experimental setup->turning detector signals into physical objects)
%
%Theorie-Part: Standardmodell, beyond SM?(eigentlich klein hakten und eher versteckt von hinten rum in Theorie LQ einführen...), Theorie LQ(begründung, was es erkärt, so bisschen die beyond SM Schiene),
%
%Ausgangspunkt, Forschungsfrage...
%
%Quasi eine Story erzählen: SM sehr erfolgreiche Theorie (->Theorie), aber auch nicht die ganze Geschichte (->beyond), das könnte LQ lösen (->LQ) dafür braucht man aber ein Nachweisgerät (->ATLAS, LHC...) und dass man an die Physik rankommt, braucht man Auswertung (->Datenprocessierung und btagging usw), was ist der aktuelle Stand in der Vorschung, worauf wird sich hier in der Arbeit konzentriert (->ttbar tautau)
%
\input{header.tex}
\begin{document}
\input{titlepage.tex}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter*{Summary}
Grand unified theories has been introduced to solve some of the open questions that are not explained by the standard model of particle physics like the symmetry between quarks and leptons. Within such models the connection between the weak and strong sector are established with new gauge bosons called leptoquarks. The Large Hadron Collider with its capability to collide protons at a center of mass energy of $\SI{13}{\tera\electronvolt}$ has with the {\ATLAS} detector the ideal instrument to search for new exotic particles like leptoquarks. Within the context of pair produced leptoquarks one possible decay within the minimal-Buchm\"{u}ller-R\"{u}ckl-Wyler model involves the $t\tau^{-}+\bar{t}\tau^{+}$ state. Due to this complex final state involving the heaviest quark and heaviest lepton, signal and background studies are performed at the {\ATLAS} experiment in order to determine the tau reconstruction efficiencies and fake rates. \newline
This is achieved with a set of two Monte Carlo samples that model the standard model background ($t\bar{t}$ and $t\bar{t}H$) and two Monte Carlo samples that model a hypothetical leptoquark signal with a mass point of $\SI{500}{\giga\electronvolt}$ and $\SI{1}{\tera\electronvolt}$. Furthermore datasets from the {\ATLAS} data-taking periods in 2015-2018 are used for a comparison between data and simulation. Within an event selection that addresses the final state of the leptoquark pair, fake candidates like jets faking a hadronic tau show a strong jet $p_T$ dependence on the fake rate as determined from the background samples. Tau reconstruction efficiencies for the leptoquark signal samples leads to the conclusion that the efficiency only slightly depends on the taus' reconstructed $p_T$, whereas the efficiency is dependent on the minimum and maximum separation between taus and different objects like jets, b-jets or light leptons.\newline
Kinematic distributions like the reconstructed tau $p_T$ or the sum of all objects' transverse momentum $S_T$ provide insight to the quality of signal extraction from the background and represent good discriminating variables. This fact is quantitatively confirmed in an evaluation of the discovery significance for both leptoquark samples. Within the significance consideration values for the transverse momentum of the reconstructed tau and the kinematic variable $S_T$ are proposed to cut on for a further significance increase. Furthermore the significance is extrapolated depending on the branching fraction of the leptoquarks decaying into a top-pair and tau-pair.
%kinematics
%Discovery significance + cuts
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter*{Zusammenfassung}
Gro{\ss}e vereinheitlichte Theorien stellen eine M\"{o}glichkeit dar, um ungel\"{o}ste Fragen des Standard Modells der Teilchenphyik zu erkl\"{a}ren. Ein Beispiel hierf\"{u}r ist die Symmetrie zwischen Quarks und Leptonen. Die Verbindung zwischen dem schwachen und starken Kraftsektor wird durch neu eingef\"{u}hrte Eichbosonen, den sogenannten Leptoquarks, hergestellt. Am Large Hadron Collider, der Protonen mit einer Schwerpunktsenergie von $\SI{13}{\tera\electronvolt}$ zur Kollision bringt, ist der {\ATLAS} Detektor das ideale Instrument f\"{u}r eine Suche nach neuen, exotischen Teilchen wie den Leptoquarks. Im Rahmen der Paarproduktion von Leptoquarks ist der Zustand $t\tau^{-}+\bar{t}\tau^{+}$ ein m\"{o}glicher Endzustand des Leptoquarkzerfalls, der durch das minimale-Buchm\"{u}ller-R\"{u}ckl-Wyler Modell erlaubt ist. Aufgrund seiner komplexen Zerfallsstruktur, die das massereichste Quark und das massereichste Lepton enth\"{a}lt, wurden Signal- und Untergrundstudien mithilfe des {\ATLAS} Experiments durchgef\"{u}hrt, um die Rekonstruktionseffizienzen und Fehlidentifikationen von Tau Leptonen zu bestimmen.\newline
Dies wurde mithilfe von Monte Carlo Simulationen der Untergrundprozesse ($t\bar{t}$ and $t\bar{t}H$) und von Simulationen, die das hypothetische Leptoquarksignal mit den zwei unterschiedlichen Massepunkten von $\SI{500}{\giga\electronvolt}$ und $\SI{1}{\tera\electronvolt}$ modellieren, erreicht. Au{\ss}erdem wurden Daten des {\ATLAS} Detektors f\"{u}r eine Gegen\"{u}berstellung genutzt, die w\"{a}hrend der Datennahmen aus den Jahren 2015-2018 gemessen wurden. Durch eine Ereignisselektion die dem Endzustand der Leptoquarkpaare angepasst ist, konnte anhand der Simulationen des Untergrunds festgestellt werden, dass die Fehlidentifikation von Jets als Tau Leptonen eine starke Transversalimpulsabh\"{a}ngigkeit aufweist. Die Rekonstruktionseffizienzen von Tau Leptonen f\"{u}r die Leptoquarksimulationen f\"{u}hren zu dem Schluss, dass die Effizienz nur schwach vom Transversalimpuls der rekonstruierten Tau Leptonen abh\"{a}ngen, w\"{a}hrend die Effizienz durchaus eine Abh\"{a}ngigkeit vom minimalem und maximalem Abstand zwischen Taus und verschieden Objekten wie Jets, b-Jets oder Leptonen aufweist.\newline
Kinematische Verteilungen wie der rekonstruierte Transversalimpuls der Tau Leptonen $p_T$ oder die Summe der Transversalimpulse aller physikalischer Objekte $S_T$ geben Einblick in die Qualit\"{a}t der Signalextraktion aus den Untergrundprozessen und repr\"{a}sentieren gute Diskriminanten. Dieser Umstand wird quantitativ in der Bestimmung der Endeckungssignifikanz f\"{u}r beide Leptoquark Simulationen best\"{a}tigt. Innerhalb der Signifikanzbetrachtung wurden au{\ss}erdem Werte f\"{u}r den Transversalimpuls der rekonstruierten Tau Leptonen und der kinematischen Variable $S_T$ bestimmt bei deren Beschnitt die Signifikanz weiter erh\"{o}t werden konnte. Weiterhin wurde die Signifikanz f\"{u}r unterschiedliche Verzweigungsverh\"{a}ltnisse der Leptoquarks beim Zerfall in ein top-Paar und Tau-Paar berechnet.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\tableofcontents
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter{Introduction}
The cornerstones for today's understanding of nature are based on successful theories like quantum mechanics and the theory of general relativity. Over centuries the discovery of new phenomena and the constant development of theoretical concepts helped to expand the frontiers of the knowledge - and this process is still in progress. Theoretical concepts based on empirical evidence are a great opportunity in order to get deeper insight into the nature of matter, forces, space and time with the ambition to deduce an universal description of the laws of nature. Within this evolving gain of knowledge the concept of symmetry was and is an appealing tool and indeed, symmetries are deeply connected with conservation laws formulated in Noether's theorem.\par
Four fundamental forces are known today. That is the electromagnetic, weak and strong force as well as gravity. The standard model of elementary particle physics describes the fundamental components of matter and their interactions induced by the electromagnetic, weak and strong force. Its development profited crucially from symmetry considerations and its formulation is realized with relativistic quantum field theory. Understanding particles as elements of symmetry groups made it possible to formulate the electroweak unification, bringing together the electromagnetic and weak force. Undisputably, the standard model is a powerful theory capable of describing the constituents of matter and their fundamental interactions and explains the phenomena encountered in this context. Nevertheless some open questions remain. Amongst others these are the puzzle of neutrino masses, the hierarchy problem or difficulties in formulating a theory of quantum gravity.\par
Thinking within the picture of symmetry it is very attractive to solve some of the open questions by further unifying the fundamental forces into a grand unified theory, which eventually includes gravity. This would imply the existence of new gauge bosons called leptoquarks. Leptoquarks are colored particles carrying both lepton and baryon number and act as mediators between leptons and quarks. In one of the simpler scenarios possible decay modes involve a top-quark pair and a tau lepton pair and would be a hint to such exotic particles. Leptoquark masses are possibly in reach of today's particle accelerators like the {\LHC} at the research center {\CERN} and could be detected with the {\ATLAS} detector. For this purpose it is crucial to understand tau reconstruction efficiencies for future search strategies in this specific $t\tau^-\bar{t}\tau^+$ channel. This work gives insight into tau reconstruction efficiencies and background contributions to a possible leptoquark signal.\par
The structure of the thesis is as follows: Chapter \ref{theory} gives an overview of the standard model of particle physics. Furthermore it presents issues, which are beyond the scope of the standard model and some approaches to solve these puzzles. One possibility is the introduction of leptoquarks in the context of grand unified theories. Within this chapter the research question of the thesis is discussed. Chapter \ref{experiment} describes the experimental setup for the analysis. That is the hadron collider complex {\LHC} and the {\ATLAS} detector. An introduction into physical object reconstruction can be found in chapter \ref{physObj}. Chapter \ref{analysis} presents the analysis framework with the object and event selection criteria as well as datasets used for the analysis. In chapter \ref{results} the results of the efficiency and background study are presented. Also different kinematic distributions of signal and background events are compared to data. Additionally, insight into the discovery significance of the leptoquark Monte Carlo simulations is given here.
%powerful theories, quantenmechanik und art, four fundamentala forces, tieferes verständnis durch übergeordnete Konzepte, Symmetrie als solch ein Konzept, vereinheitlichung klingt sehr a1111ttraktiv, standardmodel, historischer abriss elektro/magnetische WW->elektromagnetisch+ weak-->electroweak, starke, gravity in art, Quantengravitation? GUT = Weltformel für alles, GUT führt zu LQ, kurz LQ und deren sinn/wirken, CERN+ATLAS in reichweite von third generation lq nach mBRW model im tttautau channel, rolle der tau leptonen, was wird hier in der arbeit gemacht: eff, kin, significance, keine full search, chapter beschreibung, was wo behandelt wird (auch appendix erwähnen?) och neee...
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter{Theoretical background for the search for scalar leptoquarks}\label{theory}
This chapter describes theoretical foundations required for the search for scalar leptoquarks including the successful Standard Model of elementary particle physics evolved from the symbiosis of experimental achievements and theoretical milestones. Besides its success some issues still remain unsolved and could be a hint to physics beyond the Standard Model, giving space to introduce the Leptoquark Model as one possible extension.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{The Standard Model of particle physics}\label{SM}
A remarkable development for understanding nature is the Standard Model (SM) of particle physics, embracing physics at the most fundamental level. This quantum field theory, incorporating the conceptual frameworks of special relativity and quantum mechanics, describes the constituents of matter and the laws governing their interactions. \cite{Mann} Despite its success of being the most promising theory so far capable of explaining the observed results within its domain in agreement with empirical data, it seems not to be the complete story. \cite{Nair}\par
One of the most important concept in physics is that of symmetries, because they are deeply connected with conservation laws, following Noether`s Theorem. A physical property responds to a symmetry transformation in two different ways. It can appear in form of an invariant under symmetry transformation, leaving that property unchanged, or as covariant, changing their property induced by the symmetry transformation. %Continuous symmetries yield additive laws whereas discrete symmetries result in multiplicative laws.
\begin{table}[htbp]
\centering
\begin{tabular*}{\linewidth}{@{\extracolsep{\fill}}cccc}
\hline
\hline
\rule[-6pt]{0pt}{21pt}\textbf{group} & \textbf{defining property} & \textbf{application}
\\
\hline
\rule[-6pt]{0pt}{21pt} $U(n)$ & $n\times n\,$ unitary ($U^\dagger U=1$) & $U(1)$ electromagnetism
\\\\
\rule[-6pt]{0pt}{21pt} \multirow{2}{*}{ $SU(n)$ } & $n\times n\,$ unitary ($U^\dagger U=1$) & $SU(2)$ weak interactions
\\
\rule[-6pt]{0pt}{21pt} & with $\text{det}U=1$ & $SU(3)$ strong interactions
\\\\
\rule[-6pt]{0pt}{21pt} \multirow{2}{*}{$SO(n)$} & $n\times n\,$ orthogonal ($O^\intercal O=1$) & $SO(3)$ rotations
\\
\rule[-6pt]{0pt}{21pt} & with $\text{det}O=1$ & $SO(3,1)$ Lorentz transformations
\\
\hline
\hline
\end{tabular*}
\caption[Lie symmetry groups of the Standard Model.]{Lie symmetry groups for the gauge interactions of the Standard Model \cite{Mann}.}
\label{sym}
\end{table}
Fundamental symmetries of particle physics include space translation symmetry and hence the conservation of momentum, rotational invariance and hence conservation of angular momentum and time translation invariance leading to energy conservation. The formal mathematical description for fundamental symmetries is based on group theory. In case of particle physics almost all groups are Lie groups $\mathcal{G}$, that are a set of objects $\{g_i\}$, which can be combined with a binary operation and has four basic properties: closure, identity, inverse element and associativity. Additionally, the group elements are continuous and differentiable functions of some finite set of parameters $\theta_{a}$:
\begin{align}
g=g(\theta_1, \ldots, \theta_N)=\exp\left[i\theta_a\mathbf{T}^a\right]=\exp\left[i\vec{\theta}\vec{\mathbf{T}}\right]\qquad \text{with}\quad a=1,\ldots,N
\end{align}
Here $\mathbf{T}^a$ are the generators of the group from which all elements of the group can be created. The irreducible representatives of a group can be written as complex matrices\footnote{Irreducible means that not all representing matrices of the group can be decomposed into block-diagonal form simultaneously \cite{Mann}.}, acting on the wave function of the particles and on charges as well as on space-time coordinates. \cite{Mann}
The local symmetry $SU(3)_c\times SU(2)_L\times U(1)_Y$ summarizes the gauge interactions of the SM (see table \ref{sym}). Here $c$ indicates the strong force, $L$ the left handed chirality of the weak regime and $Y=B+s$ the hypercharge calculated from baryon number $B$ and strangeness $s$. \cite{PhysTeV} Besides the continuous symmetries above also important discrete symmetries exist in the SM like parity $P$, referring to the transformation $\vec{x}\rightarrow-\vec{x}$, time reversal $T$, referring to $t\rightarrow-t$, and charge conjugation $C$, corresponding to the exchange of a particle with its anti-particle. The weak force breaks $P$ and $C$, but not the product of $CPT$. \cite{Nair}\par
Matter and its interactions can be described by two basic types of particles and four fundamental forces, i.e. the electromagnetic force, the weak and the strong force and gravity. Fermionic particles, following Fermi-Dirac statistics, make up matter, whereas bosons, following Bose-Einstein statistics, are acting as mediators of the fundamental forces. \cite{Cottingham}\cite{Griffiths}\newline
The fermions can be further categorized into $6$ leptons $l$ characterized through the lepton quantum number $L_l$ and $6$ quarks $q$ characterized through the baryon quantum number $B$ together with their anti-particles ($\bar{l}$ respectively $\bar{q}$). The difference between particle and anti-particle are the contrary charges and contrary lepton and baryon number respectively, whereas mass, mean lifetime and spin stay the same. Leptons occur in three generations with different flavor -- electron (e), muon ($\mu$) and tauon ($\tau$) -- and can carry electrical charge $Q=\pm e$ in units of the elementary charge. The electrically neutral leptons are the neutrinos $\nu_l$. \cite{Griffiths}\newline
Table \ref{SMtable} shows the leptons with selected properties.\newline
\begin{table}[htbp]
\centering
\begin{tabular*}{\linewidth}{@{\extracolsep{\fill}}cccccc}
\hline
\hline
\rule[-6pt]{0pt}{21pt}\textbf{leptons} & & & &
\\
\hline
\rule[-7pt]{0pt}{23pt} $l$ & $L_l$ & $B$ & $Q/e$ & $m$/$\frac{\SI{}{\mega\electronvolt}}{\SI{}{\text{\ensuremath{c}\squared}}}$ & Spin $S_z/\hbar$
\\
\hline
\rule[-6pt]{0pt}{21pt} e$^-$ & \(L_e=1\) & $0$ & $-1$ & $0.511$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt}$\nu_e$& \(L_e=1\) & $0$ & $0$ & $<2\cdot10^{-6}$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt} $\mu^-$ & \(L_{\mu}=1\) & $0$ & $-1$ & $106$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt}$\nu_{\mu}$ & \(L_{\mu}=1\) & $0$ & $0$ & $<2\cdot10^{-6}$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt} $\tau^-$ & \(L_{\tau}=1\) & $0$ & $-1$ & $1.78\cdot10^3$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt}$\nu_{\tau}$& \(L_{\tau}=1\) & $0$ & $0$ & $<2\cdot10^{-6}$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt} & & & &
\\
\rule[-6pt]{0pt}{21pt}\textbf{quarks} & & & &
\\
\hline
\rule[-7pt]{0pt}{23pt} $q$ & \(L_l\) & $B$ & $Q/e$ & $m$/$\frac{\SI{}{\mega\electronvolt}}{\SI{}{\text{\ensuremath{c}\squared}}}$ & Spin $S_z/\hbar$
\\
\hline
\rule[-6pt]{0pt}{21pt} u (up) & $0$ & $\frac{1}{3}$ & $\frac{2}{3}$ & $2.2$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt} d (down) & $0$ & $\frac{1}{3}$ & $-\frac{1}{3}$ & $4.7$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt} c (charm) & $0$ & $\frac{1}{3}$ & $\frac{2}{3}$ & $1.3\cdot10^3$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt} s (strange) & $0$ & $\frac{1}{3}$ & $-\frac{1}{3}$ & $95$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt} t (top) & $0$ & $\frac{1}{3}$ & $\frac{2}{3}$ & $17\cdot10^{4}$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt} b (bottom) & $0$ & $\frac{1}{3}$ & $-\frac{1}{3}$ & $4.2\cdot10^{3}$ & $\frac{1}{2}$
\\
\rule[-6pt]{0pt}{21pt} & & & &
\\
\rule[-6pt]{0pt}{21pt} \textbf{gauge bosons} & & & &
\\
\hline
\rule[-7pt]{0pt}{23pt} boson & \(L_l\) & $B$ & $Q/e$ & $m$/$\frac{\SI{}{\giga\electronvolt}}{\SI{}{\text{\ensuremath{c}\squared}}}$ & Spin $S_z/\hbar$
\\
\hline
\rule[-6pt]{0pt}{21pt} $\gamma$ & $0$ & $0$ & $0$ & $0$ & $1$
\\
\rule[-6pt]{0pt}{21pt} Z$^0$ & $0$ & $0$ & $0$ & $91.2$ & $1$
\\
\rule[-6pt]{0pt}{21pt} W$^-$ & $0$ & $0$ & $-1$ & $80.4$ & $1$
\\
\rule[-6pt]{0pt}{21pt} W$^+$ & $0$ & $0$ & $+1$ & $80.4$ & $1$
\\
\rule[-6pt]{0pt}{21pt} g & $0$ & $0$ & $0$ & $0$ & $1$
\\
\rule[-6pt]{0pt}{21pt} H$^0$ & $0$ & $0$ & $0$ & $125$ & $0$
\\
\hline
\hline
\end{tabular*}
\caption[Overview of elementary particles with some selected properties.]{Overview of leptons $l$, quarks $q$ and gauge bosons as mediators of the fundamental forces with some selected properties and quantum numbers like electrical charge $Q$, mass $m$, lepton number $L_l$ and baryon number $B$. \cite{Nair}\cite{PhysRevD}. The fermion anti-particles are not shown, which would have opposite charges and lepton/baryon number.}
\label{SMtable}
\end{table}
The quarks also occur in three generations and carry electrical charge of either $Q=\mp\frac13 e$ or $Q=\pm\frac23 e$ in units of the elementary charge as well as color charge. The new type of charge can be introduced, because the quarks lie in the fundamental representation of $SU(3)$. Possible color charges are red, green and blue and the additional anti-colors indicate that quarks are interacting with the strong force among the electromagnetic and weak interaction. \cite{Griffiths} Quarks only occur confined in color-neutral compound systems called hadrons. Baryons are three-quark hadronic states with baryon number $1$ and mesons are quark-anti-quark hadronic states, having a baryon number of $0$. \cite{Cottingham} The reason for quark confinement can be found in the potential $V(r)$ between quarks and anti-quarks depending on distance $r$. The potential has the shape $V(r)\propto\frac{-1}{r}+\text{const}\cdot r$. When separating a quark-anti-quark pair, additional potential energy must be supplied, which can exceed the potential $V(R)>2m_q$ for more than two quark masses at distance $R$. Quantum fluctuations result in the origin of a new quark pair in between. The final state now consists of two pairs, again externally color-neutral. \cite{Nair} Table \ref{SMtable} shows the different quark flavors up, down, charm, strange, top and bottom with some characteristic properties. \newline
The bosons shown in table \ref{SMtable} are the quanta of the fundamental forces \cite{Cottingham}:
\begin{itemize}
\item{The photon $\gamma$ is the mediator of the electromagnetic force.}
\item{Three mediators Z$^0$, W$^+$ and W$^-$ for the weak force.}
\item{$8$ colored gluons as mediators for the strong force.}
\item{The Higgs boson H$^0$ as quantum of the Higgs field, providing the masses for the elementary particles.}
\end{itemize}
At the current state of research a quantum field theory of gravity is purely hypothetical. The mediator for such theory would be the graviton (see chapter \ref{beyondSM}). A second aspect is that on elementary particle scales gravity is insignificant compared\footnote{Relative strength of gravity compared to the weak interaction is $10^{-35}$ \cite{Mann}.} to all other forces and is therefore not originally considered in the SM. \cite{Cottingham}\par
%
\begin{figure}[htbp]
\begin{center}
\includegraphics[width=0.5\linewidth]{figures/SMpic.pdf}
\caption[Overview of the Standard Model.]{The Standard Model with its fermions and bosons and the involved interactions. The solid blue line indicates which particles interact with each other. Loops depict self-interaction. \cite{SMpic}}
\label{SMpic}
\end{center}
\end{figure}
%
Figure \ref{SMpic} summarizes the picture of the SM with its fermions and bosons. The lines indicate which particles interact with each other through the mediators, including self-interaction.
%Gedanken zum SM
%strong (n und p of same duplet) and weak isospin (chirality)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Beyond the scope of the Standard Model}\label{beyondSM}
Nevertheless, there are still many puzzles left which are not described by the SM. This circumstance keeps physicists well motivated to gain further progress and to push the frontiers of our understanding. \cite{Nair} In the following some of the issues are briefly described, like neutrino masses, quantum gravity, the hierarchy problem and supersymmetry.
\paragraph{Neutrino masses} are confirmed by various neutrino oscillation experiments \cite{Kajita}\cite{McDonald}, although the SM does not predict neutrino masses. The neutrino flavor states $\nu_\alpha$ with $\alpha=e,\mu,\tau$ are quantum entangled with the mass states $\nu_i$ where $i=1,2,3$ described by an unitary matrix $U_{\alpha i}$ \cite{mixing}. One possible extension of the SM, explaining neutrino masses, is the seesaw mechanism. Because of the absence of the right chirality spinor components $\Psi_{R,j}$, the mass term of the Dirac lagrangian $\bar{\Psi}_L\Psi_R+\bar{\Psi}_R\Psi_L$ cannot be formed in case of neutrinos. One possible solution is to describe the neutrino masses with a Majorana mass term, which introduces very massive right-chiral neutrinos\footnote{The large masses are consistent with the lack of evidence for right-chiral neutrinos at present day energies \cite{Nair}.} besides light weight left-chiral neutrinos. One caveat is that such right-chiral neutrinos as a consequence have to exist, although within the SM there are only left-handed neutrinos known. \cite{Nair}
\paragraph{Quantum gravity} could be the embedding of general relativity into a framework of quantum theory. Quantum theory provides a well confirmed framework for all theories describing particular interactions. Therefore it would be appealing to have a quantum formulation of gravity, following the example of all other fundamental forces and being one step closer to an unified description. From the point of view of cosmology, quantum gravity could be an encompassing theory for a more fundamental understanding where general relativity breaks down, when it comes to the initial conditions of the early universe or conditions of black holes. \cite{quantumgravity}
\paragraph{The hierarchy problem} is a current challenge in particle physics and arises from quadratic corrections to the weak scale. \cite{hierarchy} The hierarchy problem formulates the large differences of scales at which symmetries are broken. Considering a single, unifying symmetry, including the standard model symmetries $SU(3)_c\times SU(2)_L\times U(1)_Y$ (see chapter \ref{SM}), which has to be broken at scale $V$, because it is not manifested at the currently explored energy scales. A lower bound is $V\approx 10^{37}\,\text{GeV}$. At the same time spontaneous symmetry breaking takes place for $SU(2)_L\times U(1)_Y$ at scales of $v\approx\SI{246}{\giga\electronvolt}$\footnote{This is the vacuum expectation value for the Higgs field. \cite{VEV}} in order to allow for massive W$^\pm$ and Z$^0$ bosons as well es massive quarks and leptons. Within the framework of perturbation theory in quantum field theory such great scale differences are difficult to keep due to the fact that corrections of symmetry breaking are controlled by the highest masses involved. One approach for $v$ is $v^2=v_0^2\sum_{i=1}^kc_i\lambda^i V_0^2$ with the initial value $v_0$, the coupling constant $\lambda$ and general coefficients $c_i$. This shows that including the first order of correction, the equation becomes $v^2=v_0^2+c_1\lambda V_0^2$. The consequence is that $v_0^2$ has to be adjusted in such a way that it cancels out several decimal places, for instance 25, of the term $c_1\lambda V_0^2$, but the decimal places beyond that limit should add up to the scale of $\SI{246}{\giga\electronvolt}$. This fine tuning is an unsatisfactory situation, even worsened by adding more orders of quantum corrections. \cite{Nair} A non-perturbative solution to the hierarchy problem is proposed by supersymmetric models. \cite{nonperturbative}
%passt hier wahrscheinlich doch nicht so gut hin, ist ja nicht so wesentlich jetzt für mich... Was ich mir vorstellen könnte ist, dass man das als relativierung zu GUT erklärt, dass auch GUT nicht die eierlegende Wollmilchsau ist--
\paragraph{Supersymmmetry (SUSY)} is a theory, which unifies fermions and bosons. The corresponding transformation can be described by an operator in spinor form $Q_{\alpha i}$ and converts fermion fields into boson fields and vice versa. Here $\alpha$ is the spinor index and $i$ are internal degrees of freedom. The field pairs are called superpartners (sparticles) and belong to the same multiplet. The difference between particles and sparticles is the spin quantum number which distinguishes them by half an unit of $\hbar$. \cite{Vergados} Due to the lack of evidence of sparticles in the same mass range as the already known elementary particles, also supersymmetry has to be broken. \cite{PhysTeV}. {\SUSY} is very attractive, because it can solve the hierarchy problem, achieve gauge unification and can even provide a candidate for dark matter. \cite{Nagashima} There are two main categories of supersymmetric models, namely $R$-parity conserving and $R$-parity violating models. $R$ is a discrete symmetry in the coupling of particles and their superpartners, defined as \cite{Vergados}:
\begin{align}
R=(-1)^{2S+3B+L}= \begin{cases}
+1 & \, \text{for ordinary particles} \\
-1 & \, \text{for {\SUSY} particles} \end{cases}
\text{,}
\label{Rparity}
\end{align}
where $S$ is the spin, $B$ the baryon number and $L$ the lepton number. An exactly conserved $R$-parity will result in a stable lightest supersymmetric particle and the {\SUSY} particles will be produced in pairs. Whereas a $R$-parity violation will result in singly producable sparticles and all of them are intrinsically unstable. \cite{Kuze}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Leptoquarks}\label{LQmodels}
Different shortcomings of the standard model open up chances to extend the current understanding of elementary particle physics with concepts briefly introduced in chapter \ref{beyondSM}. One basic idea develops the symmetry approach further, taking the current symmetries as an example (see chapter \ref{SM}), into a grand unified theory ({\GUT}). This can be made possible by embedding the SM gauge groups $SU(3)_c\times SU(2)_L\times U(1)_Y$ into a higher symmetry $G_\text{{\GUT}}$. This group has to be broken at the so called ${\GUT}$-scale of $\sim 10^{16}\,\SI{}{\giga\electronvolt}$ to such an extent that the three SM interactions occur as separate interactions below this scale. \cite{PhysTeV} The previously mentioned issues in chapter \ref{beyondSM} already gave a clue on this unification concept in different degrees of manifestations and show that {\GUT} can be a solution for different issues, for example the gauge hierarchy problem \cite{hierarchyproblem}\cite{issuesGUT}, the generation problem \cite{issuesGUT}, the parameter problem \cite{issuesGUT} and neutrino masses \cite{neutrinoAndGUT}. {\GUT}s imply new gauge bosons called leptoquarks (LQ), which would explain the strong similarities between leptons and quarks. \cite{Nagashima}
%
%
\subsection{Leptoquarks and grand unified theories}\label{GUTandLQ}
The first {\GUT} model was proposed by Georgi and Glashow in 1974 \cite{GeorgiGlashow}. They proposed that the group $SU(5)$ incorporates all the fermions into one multiplet following the same universal coupling. This enables the possibility to transform leptons and quarks into each other with leptoquarks as mediators. Various global symmetries can embed the SM symmetries, but $SU(5)$ is the simplest one. \cite{Perkins} Although the Georgi-Glashow model had a great impact, this model has been ruled out \cite{GeorgiGlashowruledout}. It predicts a lifetime of the proton of $10^{30}\,\text{years}$, but the current experimental lower limit is $10^{33}\,\text{years}$\footnote{This value is vaild for the decay mode $p\rightarrow e^+K^0$. The predictions of the proton lifetime for other dominant decay modes -- relevant in {\GUT} models -- are higher and can be found in reference \cite{protonlifetime}.} \cite{Griffiths}\cite{protonlifetime}.\par%vllt noch einen Feynmangraph zum Proton Zerfall S. 408, Griffiths
Another unifying group would be $SO(10)$, which is broken to $SO(10)\rightarrow SU(4)_c\times SU(2)_L\times SU(2)_R$, where the leptons are treated as a forth color. \cite{PatiSalamExplanation} This model is known as Pati-Salam {\GUT} model \cite{PatiSalam}. Here also the right-handed leptons are included, acting differently than their left-chiral counterparts as it is expected for solutions of the neutrino masses for example. \cite{PatiSalamExplanation}\par
Other symmetry groups like $E_6$, inspired by superstring models, are also candidates for a {\GUT} theory. \cite{E6} $E_6$ has a rank of $6$ and is a member of the exceptional Lie groups. Furthermore it can be seen as the natural extension of $SU(5)$ and $SU(10)$ and naturally includes a greater set of particles. These additional particles can be, for instance, interpreted as dark matter candidates. \cite{E6explanation}
%
%
\subsection{Effective Leptoquark model}\label{effmodels}
%12 new mediators X (charge +-4/3, 3 colors = 6) und Y (charge +-1/3, 3 colors = 6) they couple to anti quarks and leptons (anti d liegt im selben multiplet wie e-) griffiths
%X(Q,I3,Y)=(4/3,1/2,5/3), >(Q,I3,Y)=(1/3,-1/2,5/3) Nagashima
%gauge bosons both with su3 and su2 quantum numbers. these are associated with generators E_ar^+ a =1,2,3 r=4,5- These will be indicated as (X_mu)^alpha and (Y_mu)^alpha associated with r = 4,5 respectively. They occupy the positions indicated by alpha and r. the gauge bosons associated with E_alphar^- are the adjoined of the above and will be indicated as (X_mu)_alpha and (Y_mu)_alpha or equivalently as (Xbar_mu)^alpha and (Ybar_mu)^alpha. the gauge bosons X und Y are triplets under su3 and doublets under su2 and have hypercharge xxx and xxx (wie bei griffiths) vergados
%The known fermions split into three generations each containing 15 states. for example the first generation comprises u and d each in three colors and helicity states, the e- in two helicity states and the nu_e with one helicity state. by convention su5 multiplets are written down as LH states and RH states are replaced by LH antiparticles (since CP symmetry e-_Land e+_L are equivalent, the RH states of course appear in seperate multiplets) Perkins
\begin{table}[htbp]
\centering
\renewcommand{\arraystretch}{1.2}
\begin{tabular*}{\linewidth}{@{\extracolsep{\fill}}cccc|cccc}
\hline
\hline
\multicolumn{4}{c|}{$|F|=2$ leptoquarks}& \multicolumn{4}{c}{$|F|=0$ leptoquarks}
\\
\hline
LQ & $Q/e$ & $T_3$ & decay & LQ & $Q/e$ & $T_3$ & decay
\\
\hline
$S_{0,L}$ &$-1/3$ &$0$ &$l^{-}_{L}u_L$ or $\nu_Ld_L$ &$V_{0,L}$ &$-2/3$ &$0$ &$l^{-}_{L}\bar{d}_R$ or $\nu_L\bar{u}_R$
\\
$S_{0,R}$ & & &$l^{-}_{R}u_R$ &$V_{0,R}$ & & &$l^{-}_{R}\bar{d}_L$
\\
\hline
$\tilde{S}_{0,R}$ &$-4/3$ &$0$ &$l^{-}_{R}d_R$ &$\tilde{V}_{0,L}$ &$-5/3$ &$0$ &$l^{-}_{R}\bar{u}_L$
\\
\hline
$S_{1,L}$ &$-4/3$ &$-1$ &$l^{-}_{L}d_L$ &$V_{1,L}$ &$-5/3$ &$-1$ &$l^{-}_{L}\bar{u}_R$
\\
&$-1/3$ &$0$ &$l^{-}_{L}u_L$ or $\nu_Ld_L$ & &$-2/3$ &$0$ &$l^{-}_{L}\bar{d}_R$ or $\nu_L\bar{u}_R$
\\
&$2/3$ &$1$ &$\nu_Lu_L$ & &$1/3$ &$1$ &$\nu_L\bar{d}_R$
\\
\hline
$V_{1/2,L}$ &$-4/3$ &$-1/2$ &$l^{-}_{L}d_R$ &$S_{1/2,L}$ &$-5/3$ &$-1/2$ &$l^{-}_{L}\bar{u}_L$
\\
$V_{1/2,R}$ &$-4/3$ & &$l^{-}_{R}d_L$ &$S_{1/2,R}$ &$-5/3$ & &$l^{-}_{R}\bar{u}_R$
\\
&$-1/3$ &$1/2$ &$l^{-}_{R}u_L$ & &$-2/3$ &$1/2$ &$l^{-}_{R}\bar{d}_R$
\\
\hline
$\tilde{V}_{1/2,L}$ &$-1/3$ &$-1/2$ &$l^{-}_{L}u_R$ &$\tilde{S}_{1/2,L}$ &$-2/3$ &$-1/2$ &$l^{-}_{L}\bar{d}_L$
\\
&$2/3$ &$1/2$ &$\nu_Lu_R$ & &$1/3$ &$1/2$ &$\nu_L\bar{d}_L$
\\
\hline
\hline
\end{tabular*}
\caption[Overview of the scalar and vector leptoquarks proposed by the minimal-Buchm\"{u}ller-R\"{u}ckl-Wyler model.]{Overview of the scalar (S) and vector (V) leptoquarks proposed by the minimal-Buchm\"{u}ller-R\"{u}ckl-Wyler model with their third component of the weak isospin $T_3$, electric charge $Q$ and fermion number $F$. The fourth column shows possible decay modes of the leptoquarks. \cite{Kuze}}
\label{LQstable}
\renewcommand{\arraystretch}{1}
\end{table}
%
All the {\GUT} models, which are briefly mentioned in section \ref{GUTandLQ}, share a common implication, namely the introduction of new gauge bosons called leptoquarks. This implication can be formulated as an effective theory of leptoquarks for energy scales below the {\GUT}-scale , describing the underlying physical concepts.\par
The introduction of the effective leptoquark model closely follows \cite{Kuze}.\newline
Leptoquarks are color-triplet scalar $(S)$ or vector $(V)$ bosons having baryon and lepton numbers and carrying fractional electrical charge. \cite{Kuze} A general formulation of an effective Lagrangian for leptoquark interaction with SM fermions was proposed by Buchm\"{u}ller, R\"{u}ckl and Wyler \cite{BRW}. It assumes the following:
\begin{itemize}
\item[(i)] LQs have renormalizable interactions.
\item[(ii)] LQs have interactions invariant under the SM gauge symmetries $SU(3)_c\times SU(2)_L\times U(1)_Y$.
\item[(iii)] LQs couple only to the SM fermions, gauge bosons and the Higgs boson.
\item[(iv)] LQs are required to conserve the lepton number $L$ and the baryon number $B$ separately. They carry the fermion number
\begin{align}
F=3B+L
\label{fermionnumber}
\end{align}
with $|F|=0$ or $|F|=2$.
\item[(v)] Each LQ couples to only one quark-lepton generation, i.e. there are three LQ families.
\item[(vi)] LQs has pure chiral couplings to the SM fermions.
\end{itemize}
Condition (iv) makes sure that the proton instability is avoided. Condition (v) only allows inter-generational interactions and large tree-level flavor changing neutral currents and flavor universalities. Condition (vi) avoids direct contributions to chirally suppressed meson decays like $\pi\rightarrow e\nu$. The LQ model, fulfilling conditions (i)-(vi), is the so called minimal-Buchm\"{u}ller-R\"{u}ckl-Wyler effective model (mBRW) with the Yukawa coupling\footnote{Yukawa coupling describes the interaction of a scalar field with a scalar or pseudoscalar Dirac field \cite{Peskin}.} constant $\lambda$. \cite{Kuze}\newline
Fourteen different leptoquark types are proposed and are listed in table \ref{LQstable}. The same symbol represents LQs of different electric charge within an isospin family, i.e. $S_{1/2,L}$ stands for both, the $S_{1/2}$ state of charge $-\frac53$ and for $-\frac23$. Here $l_X$ are the left-handed lepton doublets in case of $X=L$ and the right handed lepton singlet in case of $X=R$. This differentiation between left-chirality and right-chirality is also valid for the quarks $q_X$ of up-type ($q=u$) or down-type ($q=d$). Seven of the fourteen LQs are scalars ($S_{0,L}$, $S_{0,R}$, $\tilde{S}_{0,R}$, $S_{1,L}$, $S_{1/2,L}$, $S_{1/2,R}$, $\tilde{S}_{1/2,L}$) and seven are vectors ($V_{0,L}$, $V_{0,R}$, $\tilde{V}_{0,R}$, $V_{1,L}$, $V_{1/2,L}$, $V_{1/2,R}$, $\tilde{V}_{1/2,L}$) with their fermion number $|F|=0$ and $|F|=2$. For specific models based on mBRW the branching ratio $\beta(LQ\rightarrow lq)$ is usually fixed to values like $0$, $\frac12$ or $1$. \cite{Kuze}\par
From the experimentalists' point of view it makes sense to expand the phenomenology of leptoquarks and ease some restrictions of the mBRW model, giving rise to more generic models. This means a broader sensitivity to possibly new physics in the face of already determined exclusion limits. One example is to consider the branching ratio $\beta$ as a free parameter. Relaxing condition (iv) or (v) in the lepton sector could open new lepton-flavor violating decays. As a result the search for leptoquarks is more sensitive to various possibilities in the decay modes or to restrictions on the LQ-coupling. \cite{Kuze}
%https://edoc.ub.uni-muenchen.de/9100/1/Krobath_Gernot.pdf
%https://arxiv.org/pdf/hep-ph/9709356.pdf
%https://arxiv.org/pdf/hep-ex/0211048.pdf
%https://ac.els-cdn.com/037026938790637X/1-s2.0-037026938790637X-main.pdf?_tid=7fd73e94-a159-4098-a0c1-5df2ae6f102b&acdnat=1542542854_1685c5ae7ed0ff0824a75e01ecf264ea
%https://journals.aps.org/prd/pdf/10.1103/PhysRevD.8.1240
%https://journals.aps.org/prd/pdf/10.1103/PhysRevD.37.3165
%
%
\subsection{Leptoquark-like couplings in supersymmetry}
An additional generic picture can result from the point of view of supersymmetry. In $R$-parity violating {\SUSY} models squarks can produce leptoquark-like signature due to decay modes involving the Yukawa coupling. The left-chiral $\tilde{u}_L$\footnote{Usually the supersymmetric partner particles are denoted with a tilde $\left(\widetilde{\quad}\right)$.} squark couples to a $e^++d$ pair similar a leptoquark $\tilde{S}_{1/2,L}$ with an electric charge of $|Q|=\frac23 e$ would do. Accordingly the $\tilde{d}_R$ squark couples to a $e^-+u$ pair or $\nu_e+d$ pair and mimics a $S_0$ leptoquark of charge $|Q|=\frac13 e$. From an experimental point of view, the observation of such a decay is not only restricted to LQ models, but also has implications on constraints for the coupling of squarks in {\SUSY} models. \cite{Kuze}
%
%
\subsection{Leptoquark pair production in proton-proton collisions}\label{LQpp}
For the search of leptoquarks, the mBRW model (see chapter \ref{effmodels}) is taken as a basis.\newline
At proton-proton colliders, like the {\LHC}, leptoquarks can be be produced in pairs by gluon-fusion and smaller contributions of quark-fusion. These main processes are shown with feynman diagrams in figure \ref{LQpairs}. Figure \ref{gfa} and \ref{gfb} describes the gluon initiated production in the s-, t- and u-channel and figure \ref{qfc} the quark initiated production in the s- and t-channel. The last feynman diagram (lower right) is proportional to the square of the coupling constant ($\lambda^2$) due to the virtual lepton exchange \cite{hunter}\cite{Hewett}. Cross section calculations for such processes and further details can be found for example in reference \cite{Kramer}.
%
\begin{figure}
\centering
\begin{subfigure}[t]{0.9\textwidth}
\qquad\qquad
\includegraphics[width=0.4\textwidth]{figures/gfa.pdf}
\includegraphics[width=0.4\textwidth]{figures/gfb.pdf}
\subcaption{Leptoquark pair production via gluon-fusion in the t-channel and the u-channel.}
\label{gfa}
\end{subfigure}
%
\begin{subfigure}[t]{0.9\textwidth}
\centering
\includegraphics[width=0.4\textwidth]{figures/gfd.pdf}
\subcaption{Leptoquark pair production via gluon-fusion in the s-channel.}
\label{gfb}
\end{subfigure}
%
\begin{subfigure}[t]{0.9\textwidth}
\qquad\qquad
\includegraphics[width=0.4\textwidth]{figures/qfe.pdf}
\includegraphics[width=0.4\textwidth]{figures/qff.pdf}
\subcaption{Leptoquark pair production via quark-fusion in the s-channel and the t-channel.}
\label{qfc}
\end{subfigure}
%
%
\caption[Feynman diagrams of leptoquark pair production processes at proton-proton colliders like the {\LHC}.]{Feynman diagrams of leptoquark pair production processes dominated by gluon-fusion (a), (b) and smaller contributions by quark-fusion (c) at proton-proton colliders like the {\LHC}. \cite{hunter}}
\label{LQpairs}
\end{figure}
%
%
%
%
%
The detectable final states are governed by the Yukawa coupling $\lambda_{lq}$ of the leptoquark directly to the quark $q$ and lepton $l$ (see figure \ref{YukawaLQ}). The model is defined by two parameters derived from $\lambda_{lq}$: The branching ratio $\beta$ and the coupling constant $\lambda$. They are connected to the Yukawa coupling constant by $\lambda_{lq}=\sqrt{\beta}\lambda$ for a charged lepton and by $\lambda_{lq}=\sqrt{1-\beta}\lambda$ for neutrinos. \cite{currentStatus:13TeVATLAS}\par
%
\begin{figure}[htbp]
\begin{center}
\includegraphics[width=0.4\linewidth]{figures/yukawa.pdf}
\caption[Feynman diagram of a leptoquark decay governed by Yukawa coupling.]{Feynman diagrams of a leptoquark decay into a lepton $l$ and a quark $q$ governed by Yukawa coupling $\lambda_{lq}=\sqrt{\beta}\lambda$. \cite{currentStatus:13TeVATLAS}}
\label{YukawaLQ}
\end{center}
\end{figure}
%
The focus of this thesis lies on the search of pair produced scalar leptoquarks with the final state $LQ+LQ\rightarrow t\tau^{-}+\bar{t}\tau^{+}$. These final states would correspond to leptoquarks $S_{0,L}$ and $S_{1,L}$ (cf. table \ref{LQstable}).
%Futher information on the search of vector leptoquarks can be found for example in \cite{vectorLQ} and \cite{vectorLQQ}.
%https://cds.cern.ch/record/2308307/files/SUgS-18-001-pas.pdf guter Ausgangspunkt auch für LQ Literatur Suche
%
% eigentlich war geplant, dass die folgenden zwei Kapitel zu Kapitel data anaalysis gehören, allerdings findet Mahsana, dass das jetzt gleich kommen sollte
%
%
\subsection{Current status in the search for scalar leptoquarks}\label{currentStatus}
%Literatur dazu: current:<> in der .bib
Leptoquarks are a promising research field for physics beyond the SM. Various efforts are made with the physics program available at today's particle accelerators to search for LQ-like signatures. The following summarizes briefly the current status for the search for scalar leptoquarks at the {\LHC}.\par
LQ considered by the minimal-Buchm\"{u}ller-R\"{u}ckl-Wyler model (cf. section \ref{effmodels}) can result in various final states, which can be potentially detected by the {\ATLAS} or {\CMS} detector.\newline
The data samples recorded by the {\ATLAS} experiment in 2012 at a center of mass energy of $\sqrt{s}=\SI{8}{\tera\electronvolt}$ were searched for signal events of pair-produced LQs decaying to $LQ+LQ\rightarrow e^+q+e^-\bar{q}$, $LQ+LQ\rightarrow \mu^+q+\mu^-\bar{q}$, $LQ+LQ\rightarrow t\nu_\tau+q+\bar{t}\bar{\nu}_\tau$ and $LQ+LQ\rightarrow b\nu_\tau+q+\bar{b}\bar{\nu}_\tau$ at an integrated luminosity of $\SI{20.3}{\per\femto\barn}$. The search for first and second generation LQs was performed in the $e^+qe^-\bar{q}$ and $\mu^+q\mu^-\bar{q}$ channel. The integrated luminosity for the search of third generation LQs ($t\nu_\tau\bar{t}\bar{\nu}_\tau$ and $b\nu_\tau\bar{b}\bar{\nu}_\tau$) is $\SI{20.1}{\per\femto\barn}$. \cite{currentStatus:8TeVATLAS}\footnote{The search for scalar leptoquarks at $\sqrt{s}=\SI{7}{\tera\electronvolt}$ at the {\ATLAS} experiment is addressed by reference \cite{currentStatus:7TeVATLAS}.} For a branching ratio of the LQs of $\beta=1.0$ exclusion limits were found for the first and second generation LQs. The exclusion limit for first generation LQ masses are $m_{LQ1}<\SI{1050}{\giga\electronvolt}$ and for second generation LQ masses are $m_{LQ2}<\SI{1000}{\giga\electronvolt}$ at a confidence limit (CL) of $\SI{95}{\percent}$. For the third generation the mass exclusion limit at CL $\SI{95}{\percent}$ is $m_{LQ3}<\SI{625}{\giga\electronvolt}$ based on the $b\nu_\tau\bar{b}\bar{\nu}_\tau$ decay mode and $\SI{210}{\giga\electronvolt}<m_{LQ3}<\SI{640}{\giga\electronvolt}$ based on the $t\nu_\tau\bar{t}\bar{\nu}_\tau$ decay mode. \cite{currentStatus:8TeVATLAS}\par
A more recent study by the {\ATLAS} collaboration at $\sqrt{s}=\SI{13}{\tera\electronvolt}$ was published in reference \cite{currentStatus:13TeVATLAS} for first and second generation leptoquarks. The observed exclusion limits for the LQ mass at $\SI{95}{\percent}$ CL are $\SI{1100}{\giga\electronvolt}$ and $\SI{1050}{\giga\electronvolt}$ assuming a branching ratio of $\beta=1.0$. The study has increased the sensitivity for the leptoquarks and agrees with the former study \cite{currentStatus:8TeVATLAS}. The exclusion limits are extended \cite{currentStatus:13TeVATLAS} compared to the previous study.\par
The {\CMS} collaboration published a study for the search of third generation scalar leptoquarks decaying to a top quark and $\tau$-lepton at a center of mass energy of $\SI{13}{\tera\electronvolt}$.
%
\begin{figure}[htbp]
\begin{center}
\includegraphics[width=0.6\linewidth]{figures/exclusionLQ.pdf}
\caption[Upper limits of the product of the cross section and $\beta^2$ depending on the leptoquark mass for the $t\tau^-\bar{t}\tau^+$ channel.]{Upper limits at CL $\SI{95}{\percent}$ of the product of the cross section and $\beta^2$ depending on the leptoquark mass for the $t\tau^-\bar{t}\tau^+$ channel (third generation) at $\sqrt{s}=\SI{13}{\tera\electronvolt}$ \cite{currentStatus:13TeVATLAS}.}
\label{exclusionLQ}
\end{center}
\end{figure}
%
The corresponding integrated luminosity is $\SI{35.9}{\per\femto\barn}$. The results agree with the standard model expectation and no evidence for LQ pair production was found. The study excludes LQ with masses below $\SI{900}{\giga\electronvolt}$ at a branching fraction of $\beta=1.0$ and at a CL of $\SI{95}{\percent}$. Figure \ref{exclusionLQ} shows upper limits at CL $\SI{95}{\percent}$ of the product of the cross section and $\beta^2$ depending on the leptoquark mass for the $t\tau^-\bar{t}\tau^+$ channel. These results are the most stringent limits to date for the production of third generation LQs decaying to top quark and $\tau$-lepton. \cite{currentStatus:13TeVCMS}\par
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\subsection{Starting point and research question for the analysis}\label{startingpoint}
This work is enclosed in the general research context briefly described in section \ref{currentStatus} and makes a contribution to the search of pair-produced scalar LQs of type $S_{0,L}$ and $S_{1,L}$ (cf. table \ref{LQstable}) decaying to $LQ+LQ\rightarrow t\tau^{-}+\bar{t}\tau^{+}$ at the {\ATLAS} experiment. The common difficulty in searches for exotic particles is that in the measured datasets signal and background contributions appear as one entity. To reliably distinguish between signal and background events it would be ideally to break the data down into its single constituents. Of course this is not possible, but using Monte Carlo simulations provides up a chance for a better understanding of signal and background contributions in real data (see also section \ref{MCsim}). Although this work does not presents a full search strategy for third generation leptoquarks it provides insight into $\tau$ reconstruction efficiencies and background studies for the search of LQ-signal signatures based on Monte Carlo samples (see section \ref{signalANDbgEff}). Within the context of the $t\tau^{-}+\bar{t}\tau^{+}$ channel, efficiencies reflect the ability to successfully reconstruct a tau lepton in the detector. The thesis makes further contributions to fake rate estimations for objects that incorrectly mimic the signature of a tau lepton. It also presents kinematic variables (see section \ref{Kinematics}) within the event selection (see section \ref{eventselection}) that are sensitive to the decay of third generation leptoquark pairs. By taking into account findings from the kinematic distributions, the efficiency and fake rate contributions, the discovery significances based on signal samples are calculated, which allows first optimization steps (see section \ref{significance}). The used data sets are from the data taking periods $2015$, $2016$, $2017$ and partly $2018$ adding up to an integrated luminosity of $\SI{140}{\per\femto\barn}$.\par
The search for third generation LQ decaying to top quarks and $\tau$-leptons is complex because of the multiple objects resulting from the subsequent decays of the heaviest quark and heaviest lepton.
%
%
\begin{figure}
\centering
\begin{subfigure}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{figures/t_leptonic.pdf}
\subcaption[Leptonic top decay.]{Leptonic single-top decay. Black depicts the top decay, whereas blue depicts the anti-top decay. Leptons of different flavors are depicted as $l$ with the corresponding neutrino flavor $\nu_l$.}
\label{topleptonic}
\end{subfigure}
\quad
\begin{subfigure}[t]{0.45\textwidth}
\includegraphics[width=\textwidth]{figures/t_hadronic.pdf}
\subcaption[Hadronic top decay.]{Hadronic single-top decay. Black depicts the top decay, whereas blue depicts the anti-top decay.}
\label{tophadronic}
\end{subfigure}
\caption[Leptonic and hadronic top decay.]{Leptonic (a) and hadronic (b) single-top decay.}
\label{topdecay}
\end{figure}
%
The top quark decays weakly into $t\rightarrow Wq$ with $q=b,s,d$. The main decay mode involves the bottom quark. The subsequent decay of the weak gauge boson $W$ is either hadronically ($\SI{66.5\pm 1.4}{\percent}$) or leptonically ($\SI{34\pm1.0}{\percent}$) \cite{PhysRevD}. This is shown in figure \ref{topdecay}, where the black diagram paths indicates the top decay and the blue diagram path indicates the the anti-top decay. Leptons of different flavors are depicted as $l$ with their corresponding neutrinos $\nu_l$. This means that top-pairs can produce fully-leptonic decays following both the feynman diagram of figure \ref{topleptonic}, fully-hadronic decays following both the feynman diagram of figure \ref{tophadronic} and semi-leptonic decays, where one top decays leptonically and one hadronically. The overall percentage of the final state evolving from top pair pairs is shown in figure \ref{topcake}.\newline
%
\begin{figure}[htbp]
\begin{center}
\includegraphics[width=0.6\linewidth]{figures/top_pair_branching_frac.eps}
\caption[Branching fractions of top-pairs.]{Branching fractions of top-pairs decaying fully-hadronically (`alljets'), fully-leptonically (`dileptons') and semi-leptonically (`lepton+jets') \cite{D0}. }
\label{topcake}
\end{center}
\end{figure}
%
Tau leptons have various decay modes, which can be categorized by the number of charged particles involved in the decay. Main decay modes for taus include one charged particle (`1-prong') and three charged particles (`3-prong') in the decay. They are shown in table \ref{taudecays} for $\tau^-$. Here $h^{\pm}$ are charged pions $\pi^\pm$ and Kaons $K^\pm$, `neutrals' are photons $\gamma$ or neutral pions $\pi^0$. The decay modes of the $\tau^+$ can be inferred by having corresponding decay products with opposite charges.\par
Signal and background studies with various objects in the final state depend crucially on the purity of the signatures left in the detector (cf. chapter \ref{detectorsignatures}) and how they are traced back to the potential LQ. For example electrons and muons can be reconstructed reliably from their corresponding detector signatures (see chapter \ref{detectorsignatures}). Therefore the dileptonic and semi-leptonic decay modes of the top quark-pair are good candidates to study the signal and background contributions for the search of LQ. A pair of LQ provide, in addition to the hadronic $\tau$ final states, two leptons from the top-pair decay or one lepton and one jet respectively. In the face of high contributions from underlying QCD events, resulting in jet objects, dileptonic and semi-leptonic modes are more accurately to reconstruct compared to the full hadronic decays with jets only in the final state. In order to optimize the LQ signal selection the reconstruction efficiencies and fake rates of hadronically decaying taus is investigated in the $t\tau^{-}+\bar{t}\tau^{+}$ channel.
%
\begin{table}[htbp]
\centering
%\renewcommand{\arraystretch}{1.0}
\begin{tabular*}{\linewidth}{@{\extracolsep{\fill}}clc}
\hline
\hline
&\textbf{mode}& \textbf{fraction }$\Gamma_i/\Gamma$
\\
\hline
\multicolumn{3}{c}{\textbf{1-prong decay}}
\\
$\Gamma_1$ & particle$^-\geq$ 0 neutrals $\geq 0K_L^0\nu_\tau$ & $\SI{84.58\pm0.06}{\percent}$
\\
$\Gamma_2$ & $\qquad \mu^-\bar{\nu}_\mu\nu_\tau$ & $\SI{17.39\pm0.04}{\percent}$
\\
$\Gamma_3$ & $\qquad e^-\bar{\nu}_e\nu_\tau$ & $\SI{17.82\pm0.04}{\percent}$
\\
$\Gamma_4$ & $\qquad h^-\nu_\tau$ & $\SI{11.51\pm0.05}{\percent}$
\\
$\Gamma_5$ & $\qquad h^-\geq 1\pi^0\nu_\tau$ (excluding $K^0$) & $\SI{36.51\pm0.09}{\percent}$
\\
$\Gamma_6$ & $\qquad\qquad h^-\pi^0\nu_\tau$ & $\SI{25.93\pm0.09}{\percent}$
\\
$\Gamma_7$ & $\qquad h^-\geq 2\pi^0\nu_\tau$ & $\SI{10.81\pm0.09}{\percent}$
\\
$\Gamma_8$ & $\qquad\qquad h^-2\pi^0\nu_\tau$ & $\SI{9.48\pm0.10}{\percent}$
\\
\hline
\multicolumn{3}{c}{\textbf{3-prong decay}}
\\
$\Gamma_{9}$ & $h^-h^-h^+ \geq 0\,\text{neutrals}\geq 0K_L^0\nu_\tau$ & $\SI{15.21\pm0.06}{\percent}$
\\
$\Gamma_{10}$ & $\qquad h^-h^-h^+ \geq 0\,\text{neutrals}\,\nu_\tau$ (ex. $K_S^0\rightarrow \pi^+\pi^-$) & $\SI{14.55\pm0.06}{\percent}$
\\
$\Gamma_{11}$ & $\qquad h^-h^-h^+\nu_\tau$ & $\SI{9.80\pm0.05}{\percent}$
\\
$\Gamma_{12}$ & $\qquad \pi^-\pi^+\pi^-\nu_\tau$ & $\SI{9.31\pm0.05}{\percent}$
\\
$\Gamma_{13}$ & $\qquad h^-h^-h^+\geq 1\pi^0\nu_\tau$ (excluding $K^0$) & $\SI{5.09\pm0.05}{\percent}$
\\
$\Gamma_{14}$ & $\qquad\qquad h^-h^-h^+\pi^0\nu_\tau$ & $\SI{4.76\pm0.05}{\percent}$
\\
$\Gamma_{15}$ & $\qquad\qquad \pi^-\pi^+\pi^-\pi^0\nu_\tau$ & $\SI{4.57\pm0.05}{\percent}$
\\
\hline
\hline
\end{tabular*}
\caption[Main decay modes of $\tau^-$.]{Main decay modes of $\tau^-$ including one charged particle (`1-prong') and three charged particles (`3-prong'). $h^{\pm}$ are charged pions $\pi^\pm$ and Kaons $K^\pm$, `neutrals' are photons $\gamma$ or neutral pions $\pi^0$. The decay products of $\tau^+$ have opposite charges correspondingly. \cite{PhysRevD}}
\label{taudecays}
%\renewcommand{\arraystretch}{1}
\end{table}
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter{Experimental setup for the search for scalar leptoquarks}\label{experiment}
For the search for scalar leptoquarks the {\ATLAS} detector at the Large Hadron Collider ({\LHC}) is used as experimental setup, which will be described within this chapter. The general setting of the proton-proton collider located at the {\CERN} research center is the topic of section \ref{LHC}. The particle detection of the resulting collision events, taking place in the {\ATLAS} detector with its different specialized components (section \ref{ATLAS}). Section \ref{LQpp} addresses the possible leptoquark pair production in proton-proton collisions.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{The Large Hadron Collider accelerator complex}\label{LHC}
The {\CERN} (Conseil Europ\'{e}en pour la Recherche Nucl\'{e}aire) research center was founded in $1954$ and is located near Geneva, Switzerland. It became a major European joint venture on elementary and nuclear particle physics. Currently, $22$ member states are participating in that large-scale project. The scientific ambition is manifold and includes to probe the essential constituents of nature and the fundamental forces acting between them. \cite{CERNabout}\par
%
\begin{figure}[htbp]
\begin{center}
\includegraphics[width=0.8\linewidth]{figures/CERNKomplex2.jpg}
\caption[Schematic of the {\CERN} accelerator complex.]{Schematic of the {\CERN} accelerator complex with its different stages and few experiments like {\ATLAS} located at the interaction points of the {\LHC}. \cite{CERNKomplex}}
\label{complex}
\end{center}
\end{figure}
%
In the accelerator complex protons reach energies of $\SI{6.5}{\tera\electronvolt}$ by going through different accelerator stages and are brought to collisions at defined interaction points in time intervals of $\SI{25}{\nano\second}$. Particle detectors, which are positioned at these interaction points, then register signatures of the resulting collision events and the analysis of newly created particles gives insight to the nature of elementary particle physics.\cite{LHCJINST}\newline
The {\LHC} is designed as two-ring proton-proton collider. It is also possible to collide lead ions. These ion collisions are provided for the {\ALICE}\cite{ALICE} experiment. However, the {\ATLAS} and {\CMS} experiment are in principle also capable of measuring lead collisions, but currently they are focused mainly on proton events. Figure \ref{complex} shows the different acceleration stages of the {\LHC} among other {\CERN} experiments.\newline
A hydrogen bottle is the starting point of the proton injection and acceleration. Fully ionized by an electron gun in chambers the protons are guided with electric fields towards the first acceleration stages. Radiofrequency (RF) cavities are responsible for the proton acceleration. Moreover the protons are accelerated in packages, known as bunches. \cite{LHCJINST}\cite{LHCSchule} Starting from the injection, protons will gain a kinetic energy of $\SI{50}{\mega\electronvolt}$ in the linear accelerator {\LINAC}2 and will be further transferred to the Proton Synchrotron Booster ($\SI{1.4}{\giga\electronvolt}$), the Proton Synchrotron ($\SI{25}{\giga\electronvolt}$), the Super Proton Synchrotron ($\SI{450}{\giga\electronvolt}$) and finally to the {\LHC} ring with its $\SI{26.7}{\kilo\meter}$ circumference. \cite{CERNabout} Dipole magnets are responsible for directing the protons through the stages and keep them on closed paths. Quadrupole magnets focus the beam to compensate the repulsive forces between the particles and to ensure small beam diameters at the interaction points for high collision rates. Conditions for a stable proton beam are diverse, including high vacua of $\SI{E-10}{\milli\bar}$ to $\SI{E-11}{\milli\bar}$ to avoid collisions with gas atoms and temperatures of $\SI{1.9}{\kelvin}$ for the superconducting NbTi-magnets of the accelerator. \cite{LHCJINST}\cite{LHCSchule}\par
Different experiments like \ALICE\cite{ALICE}, {{\LHC}}b\cite{LHCb} are located at the {\LHC} due to the variety of research questions. The general-purpose detectors {\ATLAS} and its counterpart {\CMS}\cite{CMS} are specialized to high luminosities\footnote{For the physical definition see eq. (\ref{Lumi}).}, exploiting the large range of physics opportunities made possible by the {\LHC}. Main tasks of the {\ATLAS} experiment are precision measurements of the SM (see chapter \ref{SM}), a better understanding of Quantum Chromo Dynamics (QCD) and the search for supersymmetric models, and new physics, among others. With the {\LHC} collision rate of $10^9$ inelastic events per second, up to $23$ simultaneous events per bunch-crossing at dominating high QCD cross sections require a powerful detector that is capable of recognizing the characteristic signatures. These circumstances are the reasons for the demands on the {\ATLAS} detector, including fast electronic elements, high detector granularity, handling high particles fluxes and reducing overlapping events at a large acceptance and coverage region. \cite{ATLASJINST}
%AußerdeM: Hinweis auf high Lumi LHC als Art mini-Ausblick in diesem Kapitel https://arxiv.org/pdf/1705.08830.pdf, https://cds.cern.ch/record/1711887/files/ATL-COM-UPGRADE-2014-014.pdf ################
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{The ATLAS detector at the LHC}\label{ATLAS}
One of the general purpose detectors for proton-proton collisions is the {\ATLAS} detector. This $\SI{25}{\meter}$ tall detector is located at one interaction point of the {\LHC} where bunches, consisting of approximately $\SI{E11}{}$ protons, collide at a rate of $\SI{40}{\mega\hertz}$ \cite{ATLASJINST}. The number of particles encountered per time is given by \cite{Perkins}
\begin{align}
\dot{N}=\mathcal{L}\sigma
\end{align}
with the cross section $\sigma$ for the present event and the instantaneous luminosity $\mathcal{L}$. The luminosity depends only on the characteristics of the collider and is often used as key parameter in collider physics \cite{LHCJINST}.
\begin{align}
\mathcal{L}=\frac{N_bn_bf_{\text{rev}}\gamma_r}{4\pi\epsilon_n\beta^*}F
\label{Lumi}
\end{align}
Where $N_b$ is the number of particles per bunch, $n_b$ the number of bunches per beam, $f_{\text{rev}}$ the rotational frequency, $\gamma_r$ the Lorentz factor, $\epsilon_n$ the normalized transverse beam emittance, $\beta^*$ the betatron function at the collision point and $F$ respects the geometric luminosity reduction factor due to the crossing angle at the collision point. The luminosity measured at the {\ATLAS} interaction point exceeded the design luminosity of $\mathcal{L}=\SI{2.05E34}{\per\square\centi\meter\per\second}$ by a factor of $2.05$ on the $2^{\text{nd}}$ of November 2017, proving the great success of the {\LHC} machine setup over the years \cite{designLumiExceeded}.\par
%Erklärung für Lumi auch aus Perkins? ##########
To meet the requirements of a general-purpose detector, several considerations for the detector design were taken into account: On the one hand there is the demand to be sensitive to the great variety of particles governed by the fundamental forces (see chapter \ref{SM}). On the other hand the yet unknown physics should be potentially covered as well by the {\ATLAS} detector -- all seen from the angle of performing precision measurements.
%Motivation eher sogar von der Seite, wieso so gebaut, wie er ist --> Teilchen detektieren (Verweis auf SM Kaptiel)-> Schichtaufbau-> etc... #######
The basic structure of {\ATLAS} is shown in figure \ref{ATLASandCoordinate} with its different sub-detector systems together with the convention for the coordinate system. %Koordinatensystem
%
\begin{figure}
\centering
\begin{subfigure}[c]{0.95\textwidth}
\includegraphics[width=\textwidth]{figures/ATLASDesign.png}
\subcaption[The structure of the {\ATLAS} detector and its sub-systems.]{The layered structure of the {\ATLAS} detector at the {\LHC} with its sub-systems inner detector, calorimeter, magnets and muon spectrometer \cite{ATLASJINST}.}
\label{ATLASDesign}
\end{subfigure}
\begin{subfigure}[c]{0.95\textwidth}
\includegraphics[width=1\textwidth]{figures/ATLASDesignCoordinate.png}
\subcaption[Definition of the global {\ATLAS} coordinate system.]{The global {\ATLAS} coordinate system formulated in cylindric coordinates with the $z$-axis parallel to the beam line and the transverse plane defined through azimuthal angle $\phi$ and pseudorapidity $\eta$. Based on \cite{ATLASJINST}.}
\label{Coordinate}
\end{subfigure}
\caption[Structure of the {\ATLAS} detector and the coordinate system.]{Structure of the {\ATLAS} detector and the coordinate system.}
\label{ATLASandCoordinate}
\end{figure}
%
The nominal interaction point acts as origin of the coordinate system, where the $z$-axis follows the beam line counterclockwise. Perpendicular to the $z$-axis lies the transverse $x$-$y$-plane usually described through the azimuthal angle $\phi$. The positive $x$-axis points towards the center of the {\LHC}. The cylindric symmetry of the detector suggests a cylindric coordinate system with the angle $\theta$ starting from the beamline. \cite{ATLASJINST} Since the polar angle is not a Lorentz invariant quantity, it is useful to describe the position in terms of rapidity \cite{LHCJINST} $w=\frac12\ln{\frac{E+p_zc}{E-p_zc}}$ in that highly relativistic regime. In the limit of large momenta, i.e. $|\mathbf{p}|c\approx E$, the rapidity coincides with the pseudorapidity formulated as \cite{ChinaPseudorapidityBook}
\begin{align}
\eta=-\ln{\tan\frac{\theta}{2}}\text{.}
\label{pseudorapidity}
\end{align}
This variable only depends on the polar angle is therefore the adequate quantity in the context of collision experiments, where usually the angle $\theta$ from the beamline is measured. \cite{ChinaPseudorapidityBook}\par
The collision products, originating from the interaction point at the center of the {\ATLAS} detector, will propagate through the different layers of the detector's sub-systems. The inner detector is responsible to track electrically charged particles. The magnetic fields for identifying the particles' electric charge and for momentum measurement is provided by the {\ATLAS} magnet system. The next layer includes the calorimeter. The electromagnetic calorimeter measures mainly the energy deposit of electrons, positrons and photons and allows the distinction between them with the help of the track information given by the inner detector. Hadronic particles will deposit their energy in the hadronic calorimeter. The outermost layer consists of the muon system responsible for a second independent measurement of the muon and anti-muon tracks. \cite{ATLASJINST}\par
\textbf{The magnet system} includes a superconducting solenoid with a field strength of $\SI{2}{\tesla}$ surrounding the inner detector as well as three large superconducting toroid magnets composed in an eight-fold azimuthal symmetry around the calorimeter. The barrel toroid magnet delivers a field strength of $\SI{0.5}{\tesla}$ and in the end-cap a field of $\SI{1}{\tesla}$ is present. \cite{ATLASJINST}\newline%Bending particle path -->identification... ######
\textbf{The inner detector} is responsible for momentum measurements, vertex measurements and pattern recognition for the identification of electrically charged particles. This is achieved with a combination of semiconductor pixel and microstrip trackers (SCT). The Insertable B-Layer (IBL) is the innermost layer of the pixel detectors at a radius of $\SI{3.3}{\centi\meter}$ away from the beam line. Additional straw tube tracking detectors are sensitive to transition radiation (TRT) located in the outer part that are responsible for high vertex and momentum resolution. The $R-\phi$ segmented pixel detectors are of a size of $50\times \SI{400}{\square\micro\meter}$ and the SCTs with their $8$ strip layers cover together a range of $|\eta|<2.5$. Typically $36$ hits per track are provided by the $\SI{4}{\milli\meter}$ straw tubes of the TRTs, which cover the range $|\eta|\leq 2.0$. \cite{IBL}\cite{ATLASJINST}\newline %Unterscheidung teilchen anhand der Übergangsstrahlung ####################
Liquid argon electromagnetic sampling \textbf{calorimeters} with high granularity allow for an excellent energy measurement of electrons and photons. They have a total thickness of more than $22$ radiation lengths $X_0$ in the barrel region ($|\eta|<1.475$) and more than $24X_0$ in the end-cap region ($1.375<|\eta|<3.2$). For hadronic energy measurements a scintillator-tile calorimeter covering $|\eta|<1.7$ is in operation. It is a sampling calorimeter and uses steel as absorber material and scintillating tiles as active material in conjunction with wavelength shifting fibres. Further LAr technology is used for hadronic particles in the outer pseudorapidity range up to $|\eta|=3.2$. Here copper plates provide the absorber material. The forward calorimeters extend the coverage of hadronic and electromagnetic energy measurements to $|\eta|=4.9$ and are $10X_0$ deep. \cite{ATLASJINST}\newline%Energiemessung: provide god res for high energy jets###################
The \textbf{muon system} is located in the outer layer of the {\ATLAS} detector and provides resolution for high energy muon tracks with three layered precision chambers up to the range of $|\eta|<2.7$. This is possible because of the air-cored toroid magnet system including one barrel and two end-cap magnets generating strong bending power in a large volume and delivering a mostly perpendicular magnetic field with respect to the muon trajectories. The bending power $\int{\vec{B}d\vec{l}}$ along the track of the muon $d\vec{l}$ reaches from $\SI{1.5}{\tesla\meter}$ to $\SI{5.5}{\tesla\meter}$ in the range $|\eta|<1.4$ (barrel) and up to $\SI{7.5}{\tesla\meter}$ (end-cap). The precision chambers are Monitored Drift Tubes (MDT) and in the larger pseudorapidity range Cathode Strip Chambers (CSC), which are multiwire proportional chambers. Due to the fact that the overall performance crucially depends on the alignment of the muon detectors with respect to each other and with respect to the inner detector, MDTs are equipped with a optical monitoring system with $12000$ sensors. Resistive Plate Chambers (RPC) and Thin Gap Chambers (TGC) are the constituents of the muon trigger system. \cite{ATLASJINST} \newline %Myonen kommen überall durch-->hier detektieren + Problem der Spurkrümmung für high pt--> eigenes system für die genauere myon-messung ############
Due to technology and resource limitations the data recording rate has to be reduced from $\SI{40}{\mega\hertz}$ to $\SI{1}{\kilo\hertz}$. This poses high demands on an efficient \textbf{trigger system} which is organized in two levels. This is the level $1$ hardware-based trigger and the software-based high level trigger (HLT). Level $1$ uses only a subset of the total detector information making basic decisions to flag so called regions of interest within the detector's coordinate system. Searches include patterns for high transverse momenta of muon tracks, electrons and photons as well as jets or large missing energy balances. The output rate after this first selection accounts for $\SI{100}{\kilo\hertz}$. The HLT is responsible for selecting the level $1$ triggered regions at full granularity and precision, reconstructing the complete events. This event filter is the final stage and achieves data reduction down to the final data-taking rate of $\SI{1}{\kilo\hertz}$, writing events of the size of approximately $\SI{2.4}{\giga\byte\per\second}$ to the disks. The HLT selection criteria are implemented using offline analysis procedures. \cite{ATLASJINST}\cite{newATLASJINST} %efficient triggering with sufficient background rejection ###############
% XXXXXXXXXXXXXXXXXXXXXXXXXXX alte Version XXXXXXXXXXXXXXXXXXXXXXXX
%three levels. Level $1$ uses only a subset of the total detector information making basic decisions to flag so called regions of interest within the coordinate system. Searches include patterns for high transverse momenta of muon tracks, electrons and photons as well as jets or large missing energy balances. The output rate after this first selection accounts for $\SI{75}{\kilo\hertz}$. The high level trigger $2$ and $3$ are responsible for selecting the level $1$ triggered regions at full granularity and precision. The level $3$ event filter is the final stage and achieves data reduction down to the final data-taking rate of $\SI{200}{\hertz}$, writing events of the size of approximately $\SI{1.3}{\mega\byte}$ to the disks. The event filter's selection criteria are implemented using offline analysis procedures. \cite{ATLASJINST}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter{Turning detector signatures into physical objects}\label{detectorsignatures}
All components of the detector system of {\ATLAS} (see chapter \ref{ATLAS}) deliver electronic signals of the proton-proton collisions, which have to be reconstructed to physical objects for the analysis. This chapter describes the reconstruction of the physical objects like electrons, muons, jets and b-flavor jets and taus, which are relevant to the analysis in this thesis. Furthermore, the role of Monte Carlo simulations in the analysis is presented briefly.
\section{Electron reconstruction at ATLAS}\label{eReconstruction}
Electrons and positrons are stable and abundant final states of collision products. Therefore, relying on signatures of such particles in the detector is essential. Electrons and positrons\footnote{In this section, the term `positron' is absorbed in the term `electron'.} give rise to tracks in the inner detector and deposit their total energy in the electromagnetic calorimeter. The tracks and calorimeter signals are used for the electron reconstruction. \cite{ePerformance}\par
For the reconstruction of electrons in the central region ($|\eta|<\SI{2.47}{}$) several steps are necessary. At first, with respect to the granularity of the calorimeter, it is searched for electron cluster seeds in the so called step of \textit{seed-cluster reconstruction}. For more information on clusters see section \ref{jets}. The efficiency of this search is $\SI{95}{\percent}$ to $\SI{99}{\percent}$ for transverse energies up to $\SI{7}{\giga\electronvolt}$ and $>\SI{15}{\giga\electronvolt}$ respectively. The \textit{track reconstruction} is responsible for the pattern recognition in the inner detector taking into account the energy loss due to mainly bremsstrahlung in the detector material of up to $\SI{30}{\percent}$. The track seed, consisting of at least three hits in different layers of the silicon detectors, is extended to a full track of at least seven hits and is extrapolated to the regions of interest at the electromagnetic calorimeter. This is done by the {\ATLAS} Global $\chi^2$ Track Fitter \cite{trackfitter}. Afterwards a loose matching of the tracks to the clusters in the electromagnetic calorimeter is done in the \textit{electron specific track fit}. The final step, \textit{electron candidate reconstruction}, is the selection of the final track candidate out of the different tracks with the best matching to the initial cluster seed. \cite{ePerformance}\newline
Both, the information of the track and the energy cluster, are used for the four-momentum calculation of the electrons. Algorithms for reconstructed electron identification are applied to distinguish signal-like or background-like electron candidates. For that the TRT likelihood method plays an important role, which uses the high-threshold hits of each TRT. Several discriminating variables are evaluated\footnote{For further details see \cite{ePerformance}} for this likelihood method (LH). Based on the outcome, three quality categories for the electron identification are defined: \textit{Loose} relies on information of the hadronic calorimeter and the first two layers of the electromagnetic calorimeter. \textit{Medium} adds information from the TRTs, the transverse impact parameter and the third layer of the electromagnetic calorimeter, whereas the \textit{Tight} operating point additionally considers track-cluster matching variables like the ratio of the cluster energy over the track momentum. \cite{ePerformance}\newline
For a well reconstructed electron additional requirements on the isolation of the cluster and the track need to be met. In this thesis the definition of the operating point \textit{Gradient} is relevant. The efficiency depends linearly on the transverse energy $E_T$. For Gradient the isolation efficiency used in the operating point optimization is $\SI{0.1143}{\percent\per\giga\electronvolt}\cdot E_T+\SI{92.14}{\percent}$ for the calorimeter cluster and the track. \cite{ePerformance}\newline
The high level electron trigger has a transverse energy threshold of at least $\SI{24}{\giga\electronvolt}$ and requires at least a medium quality electron based on the likelihood categories described in section \ref{eReconstruction} (\texttt{e24\textunderscore lhmedium}). This high level trigger (HLT) is seeded from the electromagnetic calorimeter of the level 1 trigger system with a threshold of at least $E_T<\SI{20}{\giga\electronvolt}$ (\texttt{L1EM20VH}). Additionally a energy dependent veto against energy deposits in the hadronic calorimeter behind the electromagnetic cluster is implemented. Thus it acts as hadronic veto (denoted with `\texttt{H}'). Furthermore the transverse energy threshold depends slightly on the pseudorapidity to compensate for energy losses in the material directly in front of the electromagnetic calorimeter (denoted with `\texttt{V}'). The next step of a tighter trigger selection includes a transverse energy threshold of $E_T>\SI{60}{\giga\electronvolt}$ and requires a medium quality electron (\texttt{e60\textunderscore lhmedium}). Adapted for high energy electrons the trigger \texttt{e120\textunderscore lhloose} provides electrons with transverse energies of at least $\SI{120}{\giga\electronvolt}$ by requiring at the same time the loose electron working point. \cite{trigger:e:mu} The above electron trigger menu is valid for data taken in 2015. For the data-taking periods in 2016-2018 the following changes are applied:\newline
HLT thresholds are set to $\SI{26}{\giga\electronvolt}$ with the tight quality requirement (\texttt{e26\textunderscore lhtight}) and a variable cone size isolation (\texttt{ivarloose}), which uses the criteria that the sum of the track momentum must not exceed $\SI{10}{\percent}$ of the electron's transverse energy within a cone of $\Delta R=\sqrt{\Delta\eta^2+\Delta\phi^2}<0.2$ \cite{trigger:ivarloose:erkl}\cite{trigger:ivarloose:form}. Additional to that no transverse impact parameter cut is applied (\texttt{nod0}) \cite{trigger:ivarloose:erkl}. The complementary trigger for the high transverse energy regime is increased to $\SI{140}{\giga\electronvolt}$ \cite{trigger:e:mu}.
%electron identification starts from searches of the highest energy deposit in the second layer of em, immer unter dem Aspekt möglichst schneller processierung
%lowest single e trigger e24_lhmedium threshold of ET>24GeV + medium LH identification of the electron
%seeded from the electromagnetic calorimeter at level 1 with threshold ET>20GeV L1EM20VH
%energy dependent veto against energy deposited in the hadronic calorimeter behind the electronmagnetic cluster = hadronic veto H
%ET threshold slighly eta dependent to compensate for energy loss in the passive material in front of the calorimeter V
%complementary trigger for the high transverse energy regime which requires a loose electron identification with ET>120GeV
%---------------cite trigger:e:mu
%nod0 -> no transverse impact parameter cut is applied
%iloose indic ates a variable cone size isolation requ
%----------cite trigger:ivarloose:erkl
%iloose track momentum sum can not be more than 10% of the electron transverse energy within dR<0,2
%------------------cite trigger:ivarloose:form
%
%
%
%
\section{Muon reconstruction at ATLAS}\label{muReconstruction}
Muons are good tracers of particles with shorter lifetimes and thus should be reconstructed as an important information carrier of the particle collisions. Muons\footnote{In this section, the term `anti-muon' is absorbed in the term `muon'.} are reconstructed independently in the inner detector and the muon spectrometer and the resulting combination of the two gives the full muon track. \cite{muPerformance}\par
The first reconstruction step is again pattern recognition for seeds in the silicon layers from the inside towards the outside of the inner detector. Seeds are formed out of three space points. The space points are provided by the pixel detectors with their local two-dimensional space points and from the combination of a pair of SCTs. The default setting of the seeding algorithm allows three combinations for space points: all space points within the pixel detector, all within the SCTs or two within the pixel detector and one within the SCT. \cite{muInner}\newline
The muon reconstruction in the muon spectrometer is initiated by the search for hit patterns in the muon chambers to form segments. Muon track candidates are created by fitting the segments from different layers, which is done by a combinatorial search. Selection criteria include hit multiplicity and fit quality to find the optimal track. The hits associated to each track candidate, after applying an overlap removal algorithm, are fitted using a global $\chi^2$ fit. \cite{muPerformance}\newline
Different types of reconstructed muons are distinguished, depending on which information of the possible muon candidate is available in the inner detector, the calorimeter and the muon spectrometer. After a first signature combination of the involved sub-detector signatures, the energy loss of the muon traveling through the calorimeters and the detector material is included in a refit. A \textit{combined muon} is reconstructed from the independent inner detector and muon spectrometer measurement in a global refit, if both sub-detector components have registered a muon signature. \textit{Segment-tagged muons} are extrapolated from the inner detector muon track to at least one local segment of the muon spectrometer. If the track of the inner detector can be matched to energy deposits in the calorimeter, the reconstructed muon is called \textit{calorimeter-tagged} muon. Muons reconstructed only on the basis of the muon spectrometer are extrapolated back to the interaction point and therefore are called \textit{extrapolated muons}.\cite{muPerformance}\newline
Quality requirements on the track suppress background contributions and guarantee a robust momentum measurement by including discriminating variables like the $\chi^2$ of the combined track. This categorization finalizes the muon identification process. For the category \textit{Loose muons} all above listed muon types are used. The category \textit{medium muons} minimizes the systematic uncertainties and is based only on combined muons and extrapolated muons. \textit{Tight muons} only uses the combined muon type and is optimized for maximum purity at the cost of some efficiency. The \textit{high-$p_T$ muons} are designed for a maximization of the momentum resolution for tracks above $\SI{100}{\giga\electronvolt}$. \cite{muPerformance}\newline
The relevant isolation requirement used in this thesis for muons is the so called \textit{Gradient} isolation, which includes track-based and calorimeter-based discriminating variables. In this case it is the sum of transverse momenta of the tracks in a cone around the muon with $\Delta R=\sqrt{\Delta\eta^2+\Delta\phi^2}<0.3$, excluding the muon track itself. The calorimeter condition is based on the sum of transverse energies of the topological clusters around the muon in a cone of $\Delta R<0.2$. \cite{muPerformance}\cite{varcone}\newline
Trigger decisions for muons are based on the MS and the quality of the match that combines the track measurement of the MS and the ID. The requirement on the transverse momentum of the combined muon is at least $\SI{20}{\giga\electronvolt}$ with an additional loose isolation (\texttt{mu20\textunderscore iloose}). The loose isolation means that the scalar sum of the track $p_T$ values does not exceed $\SI{12}{\percent}$ of the muon's transverse momentum within a cone of $\Delta R<0.2$. This HLT is based on the decisions of the trigger level 1 seed, which selects muons with a transverse momentum of at least $\SI{15}{\giga\electronvolt}$ (\texttt{L1MU15}). The HLT trigger \texttt{mu50} has a threshold of $p_T<\SI{50}{\giga\electronvolt}$ and complements the trigger's high momentum regime without isolation conditions. The above muon trigger menu is based on the data taking periods in 2015 \cite{trigger:e:mu}. Changes to the data selection of the periods of the year 2016-2018 are as follows: \newline
The HLT transverse momentum requirement is changed to at least $>\SI{26}{\giga\electronvolt}$ and the isolation requirement is changed in such a way that the scalar sum of the track momentum must not exceed $\SI{7}{\percent}$ of the muon's transverse momentum within a cone of $\Delta R<0.3$ (\texttt{e26\textunderscore ivarmedium}) \cite{trigger:mu}.
%MS only und quality of th e matching of muon track measured independently in the ID and MS (physical property = pT)
%pT>20GeV of CB muon candidate + loose isolation: sum of track pT values in cone around dR<0,2 is set to <12% of the muon pT
%seeded by L1MU15 -> pT>15GeV
%complemented by mu50 withoout isolatoin for high pT muons with pT>50GeV
%---------------perf of the atlas trigger system in 2015
%ivarmedium ptrack compared to pt complett <7% in cone dR<0,3
%----------------- cite erkl
%
%
\section{Jet reconstruction at ATLAS}\label{jets}
Note that jets are not physical objects in the sense of elementary particles, but denote collimated streams of colored final states. Due to confinement (cf. chapter \ref{SM}) colored fragments will always result in jet signatures. \cite{PhysTeV}\newline
For the reconstruction of hadronic parts of final states like isolated hadrons, jets (see also section \ref{btagging}) and hadronically decaying taus (see also section \ref{taus}), the {\ATLAS} experiment employs clusters of topologically connected calorimeter cell signals (topo-clusters). The topo-cluster algorithm analyzes the spatial distribution of the cell signals to reconstruct the energy and direction of the incoming particle. \cite{topo}\par
The segmented lateral readout of the high granularity calorimeter allows the distinction between signals from particle showers and background events with the aid of the cell signal significance $\zeta_\text{cell}^\text{EM}$. \cite{topo} Seeds are defined by cells, where the absolute energy $|E|$ exceeds the noise level by a factor of four standard deviations. The topo-clusters are expanded by iteratively adding all neighbor cells with energies two standard deviations above the noise level. \cite{jetPerformance}\newline
For this thesis the relevant clustering algorithm is the anti-$k_t$ algorithm. The algorithm is based on the following definitions \cite{antikt}:
\begin{align}
&d_{ij}=\min(k_{ti}^{-2},k_{tj}^{-2})\frac{\Delta_{ij}^2}{R^2}\,\text{,}\\
&d_{iB}=k_{ti}^{-2}\,\text{,}
\label{antiktalgorithm}
\end{align}
where $d_{ij}$ is the distance between the objects $i$ and $j$ and $d_{iB}$ the distance between $i$ and the beam $B$ in the transverse plane. $k_t$ denotes the transverse momentum and $\Delta_{ij}^2=(\eta_i-\eta_j)^2+(\phi_i-\phi_j)^2$ denotes the cone between $i$ and $j$. The radius parameter $R$ scales the distance $d_{ij}$ with respect to $d_{iB}$ that any final pair of jets $a$ and $b$ is at least separated by $\Delta^2_{ab}=R^2$. The final clustered jet is obtained by the search of the minimum in distances $d_{iB}$ and $d_{ij}$. If the minimum is $d_{ij}$, then the sum of the four-momentum of object $i$ and $j$ is calculated and is treated as a new object. Again all distances between the objects are re-calculated until $d_{iB}$ is a minimum. This means, object $i$ is the final jet. The anti-$k_t$ algorithm is a stable jet reconstruction algorithm and is not prone to fluctuations because of adding low $k_t$-objects. \cite{antikt}
%
%
\section{b-jet reconstruction at ALTAS}\label{btagging}
The third generation quarks, i.e. top (t) and bottom (b), play a crucial role in the SM and its various extension possibilities like the Leptoquark Model due to their large masses \cite{Hansson}. Therefore, it is essential to identify hadrons containing b quarks and separating them from light-flavor quarks at hadron collider detectors like {\ATLAS}. This task is commonly referred to as b-tagging and can be seen as a classification problem with the goal to assign correct jet flavors. To this end the particle tracks in the inner detector and the jet reconstruction of clusters in the electromagnetic and hadronic calorimeter are discriminating objects. \cite{Paganini}\par
%
\begin{figure}[htbp]
\begin{center}
\includegraphics[width=0.55\linewidth]{figures/btagged.pdf}
\caption[Tracks in a b-jet.]{Signature of a b-jet with the primary and secondary vertex created relevant for b-tagging. $d_0$ is the impact parameter. \cite{Hansson}}
\label{btagged}
\end{center}
\end{figure}
%
The long lifetime of B hadrons of $\sim\SI{1.6}{\pico\second}$ allows them to travel a few millimeters in the detector. The subsequent decay of those heavy particles within a secondary vertex produces tracks with a comparably large impact parameter $d_0$ that is the shortest distance of the particle track from the primary vertex (see figure \ref{btagged}). This signature and the deduced impact parameter significance $S(d_0)=\frac{d_o}{\sigma(d_0)}$, where $\sigma(d_0)$ is the uncertainty of the impact parameter, are used by the b-tagging algorithms including five low-level and two high-level taggers. \cite{Hansson} The b-tagging algorithms rely on multivariate combinations of the information and process them to calculate a discriminant value for each jet. Thresholds on these values define the working point to provide efficient identification of b-jets. For better information processing of the combinations of large input parameters neural network models can be used. \cite{Luca} One example for a b-tagging algorithm is the the MV2 tagger which uses 24 input variables of the low-level taggers together with kinematic properties\footnote{For further details on MV2 see \cite{MV2}}. \cite{Paganini} It uses a boosted decision tree (BDT) to discriminate b-jets from light-jets and c-jets. \cite{MV2} The light-flavor and c-jet rejection as a function of b-jet efficiency for MV2c10 (red line) is shown in figure \ref{btageff}. The description `c10' indicates the fraction of c-jet composition of the training background sample, e.g. for MV2c10 the c-jet to light-jet composition is $\SI{10}{\percent}/\SI{90}{\percent}$. Three other tagger configurations are shown in comparison as well as their ratio to MV2c10. Following the MV2c10 efficiency curve different working points can be defined. The relevant working point for the physical object selection in this work is $\SI{77}{\percent}$ (cf. chapter \ref{physObj}). The c-jet rejection rate for this working point is $6$, for light-jet rejection $134$ and for $\tau$-rejection $22$ extracted from $t\bar{t}$ events. \cite{MV2c10}
\begin{figure}
\centering
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{figures/lightjetrejection.pdf}
\subcaption[Light-jet rejection versus b-jet efficiency for the MV2c10 tagger.]{Light-jet rejection versus b-jet efficiency for the MV2c10 tagger (red) in comparison to previous configurations evaluated on $t\bar{t}$ events. The ratio is evaluated with respect to MV2c10. \cite{MV2c10}}
%\label{ATLASDesign}
\end{subfigure}
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=1\textwidth]{figures/cjetrejection.pdf}
\subcaption[c-jet rejection versus b-jet efficiency for the MV2c10 tagger.]{c-jet rejection versus b-jet efficiency for the MV2c10 tagger (red) in comparison to previous configurations evaluated on $t\bar{t}$ events. The ratio is evaluated with respect to MV2c10. \cite{MV2c10}}
%\label{btageff}
\end{subfigure}
\caption[Light-flavor (a) and c-jet rejection (b) versus b-jet efficiency for the MV2c10 tagger.]{Light-flavor and c-jet rejection versus b-jet efficiency for the MV2c10 tagger. \cite{MV2c10}}
\label{btageff}
\end{figure}
%%%%%%%%%%%%%
\section{Tau reconstruction at ATLAS}\label{taus}
Final states with tau leptons\footnote{In this section, the `term anti-tau' is absorbed in the term `tau'.}, decaying hadronically, play an important role for the physics at the {\ATLAS} experiment. Five hadronic decay modes cover over $\SI{90}{\percent}$ of the overall hadronic modes, which results in one (1-prong) or three (3-prong) charged hadrons ($h^\pm$), up to two neutral pions ($\pi^0$) and a tau neutrino $\nu_\tau$. \cite{tauPerformance2}
For the reconstruction of hadronic taus $\tau_\text{had-vis}$ at the {\ATLAS} detector, the anti-$k_t$ algorithm (see eq. (\ref{antiktalgorithm})) for jet formation is used (see section \ref{jets}) with a distance parameter of $R=0.4$. The requirements for the transverse momentum is $p_T>\SI{10}{\giga\electronvolt}$ and for the pseudorapidity is $|\eta|<2.5$. The hadronic identification of the tau relies on topo-clusters in the last layer of the electromagnetic calorimeter and the hadronic calorimeter. The additional track selection requires at least two associated hits in the pixel detector/IBL in the $\tau_\text{had-vis}$ direction and at least seven hits in total in the pixel detector and the SCTs. Also requirements on the distance of closest approach are set to $|d_0|<\SI{1.0}{\milli\meter}$ in the transverse plane and to $|\Delta z_0\sin\theta|<\SI{1.5}{\milli\meter}$.\cite{tauPerformance}\newline
Due to the fact that the reconstruction of tau candidates provides very little rejection against jet background a detailed set of discriminating variables\footnote{They can be found in \cite{tauPerformance}.} are introduced. Tau identification relies on separate boosted decision trees (BDT) for each prong case. The working points for that are labeled \textit{loose}, \textit{medium} and \textit{tight} corresponding to different signal efficiency values. In case of 1-prong they are, respectively, $0.6$, $0.55$ and $0.45$ for loose, medium and tight. The working points for 3-prong decays correspond to $0.5$, $0.4$, $0.3$. \cite{tauPerformance}\newline
To ensure that electrons will be not misinterpreted as 1-prong tau decays an additional BDT is implemented as a tau-electron veto. It has the same classification of \textit{loose}, \textit{medium} and \textit{tight}, but with signal efficiencies of $0.95$, $0.85$ and $0.75$ respectively. \cite{elBDT}
%
%
\section{Monte Carlo simulations}\label{MCsim}
A steady comparison and adaption of theory and experiment is crucial to be able to draw physically reliable conclusions. This fact holds true for high energy particle collision experiments, like the {\LHC}. Therefore predictions for collision events, the creation of new particles, the understanding of decays and fragmentation of these generated particles, modeling the underlying event, predictions to the subsequent hadronisation of quarks, the calculation of interactions with the detector and predictions of final states are necessary for the search of new physics. These scenarios are an extremely challenging endeavor. Typically hundreds of particles of different kinds are involved and their energy and momentum spectra are widely distributed over different orders of magnitude. One important cornerstone to cope with the above mentioned challenges are Monte Carlo (MC) simulations. \cite{bigMC}\par
MC simulations are based on an algorithm relying on sequences of pseudo random numbers. Effectively it is a method to numerically compute multidimensional integrals over a complex region as it is often the case for matrix elements in perturbation theory. \cite{MCbook} The MC chain for the {\ATLAS} detector consists of three main steps. This includes the generation of MC events with event generators and the subsequent simulation of fragmentation and showering processes in the first step. The second step is the simulation of the detector and the propagation of the particles through it with the associated interactions and, as third step, the digitization of energy deposits in the detector's sensitive areas. This digitization step involves the simulation of the resulting signals in the readout system for comparison with real data. The detailed implemented simulation of the detector and the particles interaction with it is done with {\GEANT}4\cite{GEANT4}, which additionally considers realistic detector conditions and misalignments as well as the detector response to be as realistic as possible. This fact enhances the digital copy to furthermore test the effectiveness of proposed search strategies for specific signatures like the ones created from exotic particles. \cite{ATLASSim}
%aus Verenas MA: ATLAS detector and physics performance: Technical design report, 2 1999. CERN-LHCC-99-015
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter{Analysis}\label{analysis}
Before investigating the research question formulated in section \ref{startingpoint}, the framework principles have to be defined. The MC samples and data samples, which are taken into account, are listed in detail in section \ref{MCsamples}. The physical object definition and the event selection for the performed analysis are presented in section \ref{physObj} and \ref{eventselection}.
%das wurde zu LQ verfrachtet, siehe LQ Kapitel: besserer Zusammenhang
%\section{Current status in the search for scalar leptoquarks}
%\section{Starting point and research question for the analysis}
\section{Datasets and Monte Carlo samples}\label{MCsamples}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The motivation is to estimate tau efficiencies for the search for third generation leptoquarks in the $t\tau^{-}+\bar{t}\tau^{+}$ channel. To draw convincing conclusions from the efficiency study, a good estimation of the standard model background events has to be chosen. The choice of background samples involves the following considerations:\newline
%ttbar teil
Dominant background processes, contaminating the potential LQ signal, originate from QCD events producing top-quark pairs and top-pair production in association with Higgs, where Higgs decays into tau leptons.\newline
The differential cross section for top-quark pair generation at $\sqrt{s}=\SI{13}{\tera\electronvolt}$ is \SI{829\pm189}{\pico\barn} \cite{ttbarpairproduction} including the statistical and the systematic uncertainty as well as the uncertainty from the luminosity measurement. Due to the nature of a proton-proton collider mainly gluon fusion ($\SI{90}{\percent}$ \cite{PhysRevD}) is responsible for the top-quark pair production at the {\LHC}. Also quark fusion can lead to $t\bar{t}$ states. The corresponding Feynman diagrams involved in the pair production are shown in figure \ref{ttbarproduction}. Because of its high production cross section and the lack of naturally existing taus, the $t\bar{t}$ sample is the ideal candidate to study jets and especially b-jets that fake as a hadronic tau.
%
\begin{figure}
\centering
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{figures/ttbara.pdf}
\subcaption[Quark fusion process generating a top-quark pair in the s-channel.]{Quark fusion process generating a top-quark pair in the s-channel \cite{Feynmanttbar}.}
%\label{ATLASDesign}
\end{subfigure}
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{figures/ttbarb.pdf}
\subcaption[Gluon fusion process generating a top-quark pair in the t-channel.]{Gluon fusion process generating a top-quark pair in the t-channel \cite{Feynmanttbar}.}
%\label{btageff}
\end{subfigure}
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{figures/ttbarc.pdf}
\subcaption[Gluon fusion process generating a top-quark pair in the u-channel.]{Gluon fusion process generating a top-quark pair in the u-channel \cite{Feynmanttbar}.}
%\label{ATLASDesign}
\end{subfigure}
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{figures/ttbard.pdf}
\subcaption[Gluon fusion process generating a top-quark pair in the s-channel.]{Gluon fusion process generating a top-quark pair in the s-channel \cite{Feynmanttbar}.}
%\label{ATLASDesign}
\end{subfigure}
\caption[Feynman diagrams showing the dominant top-pair production at the {\LHC}.]{Feynman diagrams showing the dominant top-pair production at the {\LHC} including quark fusion (a) and gluon fusion (b)-(d) \cite{Feynmanttbar}.}
\label{ttbarproduction}
\end{figure}
%
%dann noch H teil.
The total cross section measured for Higgs production is $\left( 12^{+25}_{-16} \right)\,\text{pb}$ \cite{HiggsCrossSection}. The dominant Higgs production mechanisms of leading order at hadron colliders are gluon fusion, vector boson fusion, Higgs-strahlung or the associated production with a pair of top (bottom) quarks \cite{PhysRevD}. The relevant background SM process for the search of LQs is the associated Higgs production with a top pair, which can mimic the $t\tau^{-}+\bar{t}\tau^{+}$ state. The $t\bar{t}H$ cross section is $\SI{0.51\pm0.06}{\pico\barn}$. The corresponding Feynman diagram for $t\bar{t}H$ production is depicted in figure \ref{ttHproduction}. The used $t\bar{t}H$ sample (cf. section \ref{MCsamples}) acts as an inclusive sample, which allows all possible Higgs decays. The Higgs coupling strength is proportional to the particle's mass and therefore the direct decay of the Higgs boson into a tau-pair is possible, as it involves the heaviest lepton (cf. \ref{SMtable}). Together with the decay $H\rightarrow b\bar{b}$ the mode $H\rightarrow \tau^-\tau^+$ is the dominant process below the threshold of $\SI{160.8}{\giga\electronvolt}$ and $\SI{182.4}{\giga\electronvolt}$ that is the on-shell $\text{H}^0\rightarrow\text{W}^-\text{W}^+$ and $\text{H}^0\rightarrow\text{Z}^0\text{Z}^0$ production \cite{PhysRevD}. On-shell means that both vector bosons are produced with their likeliest mass rather than producing one vector boson of the pair as a virtual particle. From the point of view of the Higgs boson's coupling to large masses in principle $H\rightarrow t\bar{t}$ is also a possible decay mode. But the di-top mass exceeds the Higgs mass by a factor of roughly three times, which suppresses that decay due to the fact that the top-pair can exist only virtually.
%
\begin{figure}[htbp]
\begin{center}
\includegraphics[width=0.3\linewidth]{figures/ttHproduction.pdf}
\caption[Higgs production associated with a top-quark pair.]{Feynman diagram in leading order for the top-quark pair associated Higgs production \cite{PhysRevD}.}
\label{ttHproduction}
\end{center}
\end{figure}
%
The above mentioned aspects make the $t\bar{t}H$ sample a good candidate to study faked taus from jets with the additional condition that also real taus are naturally present compared to the top-pair production. Furthermore efficiencies for different tau selections can be studied as well in this events as they occur as signal signatures. Findings in the background processes then can be compared to the LQ signal samples (see section \ref{signaleff}). \newline
Together these to background samples cover most of the contamination of the $t\tau^{-}+\bar{t}\tau^{+}$ for the LQ search.\par
The MC samples in the analysis are used to investigate the signal and background composition and to estimate fake rates related to tau reconstruction of LQs in the $t\tau^{-}+\bar{t}\tau^{+}$ channel. The analysis uses two types of MC samples. One category includes the LQ signal following the decays after the mBRW model (cf. section \ref{effmodels}). These signal samples were produced with the event generator {\aMCNLO} (v2.6.0) \cite{aMCNlo} with two different mass points of the LQ of $m_{LQ}=\SI{500}{\giga\electronvolt}$ and $m_{LQ}=\SI{1}{\tera\electronvolt}$ and with a production cross section of $\SI{271}{\femto\barn}$ $\SI{3}{\femto\barn}$ respectively. The top pair is allowed to decay either fully leptonically or semi-leptonically in the LQ sample (cf. section \ref{startingpoint}). The parton showering and fragmentation as well as the underlying event generation were done with {\Pythia} \cite{Pythia8} with the particle distribution function (PDF) sets {\NNPDFd} \cite{NNPDF} and the A14 \cite{A14} tunes at next-to-leading order (NLO) in matrix calculation. \newline
The other category of simulations includes main background contributions of standard model processes as motivated above. That is background from $t\bar{t}$ pairs and inclusive $t\bar{t}H$ decay processes. The $t\bar{t}$ samples are generated with {\POWHEG} \cite{Powheg} and allow the $t\bar{t}$ pair to decay fully leptonically. Showering, fragmentation and the underlying event generation are simulated with the generator {\Pythia} at NLO with the PDF set of {\NNPDFz} and the parameter set of A14. The $h_\text{damp}$ parameter was set to $\SI{258.75}{\giga\electronvolt}$, which is a resummation damping factor and regulates the highest transverse momentum of the first additional parton radiation beyond the Born level \cite{hdamp}\cite{hdamp2}.\newline
%
\begin{table}[htbp]
\centering
\renewcommand{\arraystretch}{1.2}
\begin{tabular*}{\linewidth}{@{\extracolsep{\fill}}cc}
\hline
\hline
\textbf{data-taking year}& \textbf{integrated luminosity}
\\
\hline
$2018$ & $\SI{60}{\per\femto\barn}$
\\
$2017$ & $\SI{44}{\per\femto\barn}$
\\
$2016$ & $\SI{33}{\per\femto\barn}$
\\
$2015$ & $\SI{3}{\per\femto\barn}$
\\
\hline
\textbf{total} & $\SI{140}{\per\femto\barn}$
\\
\hline
\hline
\end{tabular*}
\caption[Integrated luminosity for the data-taking periods 2015-2018.]{Integrated luminosity for the data-taking periods 2015, 2016, 2017 and 2018.}
\label{datalumi}
\renewcommand{\arraystretch}{1}
\end{table}
%
The $t\bar{t}H$ sample is generated with {\POWHEG}. Like for the $t\bar{t}$ samples, the event generator {\Pythia} at NLO with the PDF set of {\NNPDFz} and the parameter set of A14 was used to simulate the showering, fragmentation and underlying event. The mass of the Higgs boson is set to $\SI{125}{\giga\electronvolt}$ and the resummation parameter was set to $\SI{325.5}{\giga\electronvolt}$. \newline
For all samples the {\EvtGen} (v1.6.0) \cite{EvtGen} program was used for the decay of the bottom and charm hadrons produced by the top-quark decay.\newline
Data sets of the data-taking periods 2015, 2016, 2017 and 2018, which were produced at a center of mass energy of $\sqrt{s}=\SI{13}{\tera\electronvolt}$, are studied and compared to the MC predictions. Table \ref{datalumi} shows the integrated luminosity of the different years and the total integrated luminosity of $\SI{140}{\per\femto\barn}$ used in the analysis. Because of the fact that for this study no systematic uncertainties are considered, the luminosity values are rounded to the next whole number. Their uncertainty contributions of usually few percent are small compared to uncertainties from systematics.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%$
\section{Physical object selection}\label{physObj}
The physical object selection for electrons, muons, hadronic taus, jets and b-jets in the analysis is described in this section. For the general reconstruction concepts of physical objects see chapter \ref{detectorsignatures}.\par
The {\ATLAS} software framework for event generation, simulation, reconstruction and derivation production (cf. also section \ref{MCsim}) is Athena \cite{Athena}. %The management and build of the software packages rely on CMake.
The standard package for physical object reconstruction used in this analysis is AnalysisTop. It is derived from the general software chain of AthAnalysisBase specifically adapted for the needs of the top-quark analysis group. The \texttt{C++} framework {\ROOT} is used as data analysis tool. For this work a package based on AnalysisTop of release 21.2.56 (build on November 29, 2018) was used.\par
\subsubsection*{Electrons}
Electrons have to fulfill the tight quality requirement and need to have at least a transverse momentum of $p_T>\SI{27}{\giga\electronvolt}$ and lie within the pseudorapidity range of $|\eta|<2.5$. In order to get a clean signal from a leptonic decaying top quark, the electrons have to meet the $p_T$ dependent isolation requirements. For the electron trigger requirements see the paragraph `Lepton trigger' below.
\subsubsection*{Muons}
Muons have to meet the medium quality requirement and need to have a transverse momentum of $p_T>\SI{27}{\giga\electronvolt}$. The selected muons are reconstructed as a combined muon within a pseudorapidity region of $|\eta|<2.5$. Analogous to the electron selection for a clean lepton signal, the muon has to further fulfill the $p_T$ dependent isolation requirements. For the muon trigger requirements see the paragraph `Lepton trigger' below.
\subsubsection*{Taus}
Taus selected for this analysis need to have a transverse momentum of at least $p_T>\SI{25}{\giga\electronvolt}$ and must lie within the pseudorapidity range of $|\eta|<2.5$. Additional to that the taus have to decay hadronically. Following their basic reconstruction procedure, taus have to meet the conditions of the loose working point (see section \ref{taus}). To avoid misinterpretations of electrons as taus, the corresponding working point is set to loose.
\subsubsection*{Jets}
Jets are reconstructed with the anti-$k_T$ algorithm and need to have a transverse momentum of at least $p_T>\SI{25}{\giga\electronvolt}$.
\subsubsection*{b-jets}\label{physObj:bjet}
For the b-jet selection the b-tagger MV2c10 is used with a working point of $\SI{77}{\percent}$ b-tagging efficiency.
\subsubsection*{Lepton trigger}\label{leptontrigger}
The trigger requirements relevant for the analysis on electrons and muons are described here and are categorized each after the corresponding data-taking period.\par
Electrons are required to to pass the `OR' condition of the single electron triggers. This means that for an electron at least one of the triggers presented in table \ref{etriggersTable} had indicated an electron signature for the specific data-taking period in the detector.\par
%
\begin{table}[htbp]
\centering
\footnotesize
\renewcommand{\arraystretch}{1.2}
\begin{tabular*}{\linewidth}{@{\extracolsep{\fill}}ccc}
\hline
\hline
\textbf{data-taking year}& \textbf{`OR' required trigger} & \textbf{important parameter}
\\
\hline
\multirow{3}{*}{$2018$, $2017$, $2016$} & \texttt{e26\textunderscore lhtight\textunderscore nod0\textunderscore ivarloose} & $p_T>\SI{26}{\giga\electronvolt}$, tight working point
\\
& or \texttt{e60\textunderscore lhmedium\textunderscore nod0} & $p_T>\SI{60}{\giga\electronvolt}$, medium working point
\\
& or \texttt{e140\textunderscore lhloose\textunderscore nod0} & $p_T>\SI{140}{\giga\electronvolt}$, loose working point
\\
\hline
\multirow{3}{*}{$2015$} & \texttt{e24\textunderscore lhmedium\textunderscore L1EM20VH} & $p_T>\SI{24}{\giga\electronvolt}$, medium working point
\\
& or \texttt{e60\textunderscore lhmedium} & $p_T>\SI{60}{\giga\electronvolt}$, medium working point
\\
& or \texttt{e120\textunderscore lhloose} & $p_T>\SI{120}{\giga\electronvolt}$, loose working point
\\
\hline
\hline
\end{tabular*}
\caption[List of single electron triggers.]{List of single electron triggers.}
\label{etriggersTable}
\renewcommand{\arraystretch}{1}
\end{table}
%
The same strategy as it is used for electrons is valid for muons. Muons are required to to pass the `OR' condition of the single muon triggers. This means that for an muon at least one of the triggers presented in table \ref{mutriggersTable} had indicated a muon signature for the specific data-taking period in the detector.\par
%
\begin{table}[htbp]
\centering
\footnotesize
\renewcommand{\arraystretch}{1.2}
\begin{tabular*}{\linewidth}{@{\extracolsep{\fill}}ccc}
\hline
\hline
\textbf{data-taking year}& \textbf{`OR' required trigger} & \textbf{important parameter}
\\
\hline
\multirow{2}{*}{$2018$, $2017$, $2016$} & \texttt{mu26\textunderscore ivarmedium} & $p_T>\SI{26}{\giga\electronvolt}$, medium working point
\\
& or \texttt{mu50} & $p_T>\SI{50}{\giga\electronvolt}$
\\
\hline
$2015$ & \texttt{mu20\textunderscore iloose\textunderscore L1MU15} & $p_T>\SI{26}{\giga\electronvolt}$, loose working point
\\
\hline
\hline
\end{tabular*}
\caption[List of single muon triggers.]{List of single muon triggers.}
\label{mutriggersTable}
\renewcommand{\arraystretch}{1}
\end{table}
%
For more detailed information about the trigger requirements and their explanation see section \ref{eReconstruction} and \ref{muReconstruction}.
\subsubsection*{Overlap removal}\label{overlapremoval}
To avoid misinterpretations an overlap removal between different objects acts like a veto and is applied for the event reconstruction. The overlap is determined by thresholds in the separation variable $\Delta R=\sqrt{\Delta\eta^2+\Delta\phi^2}$. The following overlap removals are applied for objects like electrons $e$, muons $\mu$, tau leptons $\tau$ and jets $j$:
\begin{itemize}
\item electron removal if $\Delta R(e,j)<0.4$
\item muon removal if $\Delta R(\mu,j)<0.4$
\item jet removal if $\Delta R(j,e)<0.2$
\item tau removal if $\Delta R(\tau,l)< 0.2$ with $l=e,\mu$
\end{itemize}
%https://indico.cern.ch/event/784955/contributions/3265310/attachments/1776237/2890543/SR_definitions.pdf
%overlap removal: default between <e,jet>, <mu,jet>, <jet,e> 0.2 nicht tau, 0.4 tau
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Event selection}\label{eventselection}
The signal and background studies for the $t\tau^{-}+\bar{t}\tau^{+}$ channel involves the event selection shown in figure \ref{selection}.\newline
%
\begin{figure}[htbp]
\begin{center}
\includegraphics[width=0.45\linewidth]{figures/selection.pdf}
\caption[Overview of the event selection.]{Overview of the event selection involved for the signal and background studies in the $t\tau^{-}+\bar{t}\tau^{+}$ channel.}
\label{selection}
\end{center}
\end{figure}
%
The selection includes one electron and one muon of opposite electrical charge and a transverse momentum of at least $p_T>\SI{27}{\giga\electronvolt}$. This di-lepton selection covers the leptonic decay of the top-pair mediated by the $W^-$ boson in case of the anti-top quark and mediated by the $W^+$ in case of the top quark. Due to conservation of the electric charge the resulting leptons of the subsequent decay of the $W^\pm$ boson is of opposite sign. In addition to that at least two jets are selected because the weak top-pair decay produces a bottom-pair (cf. figure \ref{topdecay}), which will result in a jet. At least two of the above selected jets have to be b-tagged or at least one of the above jets have to be b-tagged. To cover the full $t\tau^{-}+\bar{t}\tau^{+}$ channel, an additional selection is made to select at least either one or two tau leptons for each cases of b-jet selection. The samples include these categories of one b-tagged jet, two b-tagged jets, one tau or two taus and the corresponding distributions as separate categories.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\chapter{Results}\label{results}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Number of event candidates}\label{eventYield}
The observed and predicted number of event candidates for the event selection described in section \ref{eventselection} is summarized in this chapter.\par
The event yield table \ref{eventYieldTable} includes all the main background contributions from $t\bar{t}$ and $t\bar{t}H$ processes as well as the LQ signal samples generated with two different mass points compared to the number of event candidates from data. The error consists of the statistical uncertainty only. The differentiation of the event number categories follow the event selection described in section \ref{eventselection}. Its main categories include the selection of either at least one b-jet or at least two b-jets. The decay channel of the LQ involves a pair of b-quarks originating from the subsequent top-pair decay mediated by the weak interaction (cf. section \ref{startingpoint}). After the two b-jet selection the sample is dominated by contributions from standard model $t\bar{t}$ productions. Selecting only one b-jet will be useful for a full search strategy for scalar LQs. Due to the leptoquarks' high masses of $\SI{500}{\giga\electronvolt}$ and $\SI{1}{\tera\electronvolt}$ and experimental exclusion limits as well as theoretically considerations, the assumed production cross section at the {\LHC} is small compared to abundant QCD processes of the SM\footnote{For a cross section comparison see \ref{MCsamples}, where the different samples are discussed.}. Thus the looser selection of at least one b-jet and the resulting higher event yield gives room for tighter cuts on kinematic variables to optimize the signal and at the same time looses almost nothing in background rejection. \newline
%
\begin{sidewaystable}[htbp]
\centering
\renewcommand{\arraystretch}{1.5}
\begin{tabular*}{\linewidth}{@{\extracolsep{\fill}}ccccccc}
\hline
\hline
\textbf{channel} & \multicolumn{6}{c}{\textbf{event selection}}
\\
& $\geq 1\,$b-jet & $\geq 1\,$b-jet $+$ $\geq 1\tau$ & $\geq 1\,$b-jet $+$ $\geq 2\tau$ & $\geq 2\,$b-jet & $\geq 2\,$b-jet $+$ $\geq 1\tau$ & $\geq 2\,$b-jet $+$ $\geq 2\tau$
\\
\hline
$t\bar{t}$ & $\SI{555369\pm 303}{}$ & $\SI{2793\pm 21}{}$ & $\SI{7.5\pm 1.1}{}$ & $\SI{270799\pm 210}{}$ & $\SI{765\pm 11}{}$ & $\SI{2.99\pm 0.69}{}$
\\
$t\bar{t}H$ & $\SI{517.4\pm 2.7}{}$ & $\SI{28.02\pm 0.55}{}$ & $\SI{1.554\pm 0.070}{}$ & $\SI{338.1\pm 1.6}{}$ & $\SI{14.24\pm 0.37}{}$ & $\SI{0.735\pm 0.043}{}$
\\
\hline
\textbf{total} & $\SI{555886\pm 303}{}$ & $\SI{2821\pm 21}{}$ & $\SI{9.0\pm 1.1}{}$ & $\SI{271137\pm 210}{}$ & $\SI{779\pm 11}{}$ & $\SI{3.72\pm 0.69}{}$
\\
\hline
$LQ_\text{$500\,$GeV}$ & $\SI{1304\pm 40}{}$ & $\SI{598\pm 27}{}$ & $\SI{96\pm 10}{}$ & $\SI{650\pm 28}{}$ & $\SI{287\pm 19}{}$ & $\SI{45.9\pm 6.7}{}$
\\
$LQ_\text{$1\,$TeV}$ & $\SI{14.69\pm 0.41}{}$ & $\SI{6.65\pm 0.27}{}$ & $\SI{0.745\pm 0.093}{}$ & $\SI{7.44\pm 0.29}{}$ & $\SI{3.38\pm 0.19}{}$ & $\SI{0.371\pm 0.063}{}$
\\
\hline
data & $\SI{554745\pm 745}{}$ & $\SI{3342\pm 58}{}$ & $\SI{12.0\pm 3.5}{}$ & $\SI{265437\pm 515}{}$ & $\SI{1022\pm 32}{}$ & $\SI{3.0\pm 1.7}{}$
\\
\hline
\hline
\end{tabular*}
\caption[Observed and predicted event yield for the event selection.]{The observed and predicted number of event candidates with statistical uncertainty for the event selection described in section \ref{eventselection} including the LQ signal, the background contributions of $t\bar{t}$ and $t\bar{t}H$ processes and data. The integrated luminosity is $\SI{140}{\per\femto\barn}$.}
\label{eventYieldTable}
\renewcommand{\arraystretch}{1}
\end{sidewaystable}
%
The additional selection of at least either one or at least two tau leptons impose an additional reduction on the event yield for both b-jet selections. The difference between the background simulations is that selecting $\geq 1\,\text{b-jet}+\geq 1\,\tau$ additionally results in a strong reduction of the events without any tau lepton selection down to $\SI{0.5\pm 0.75}{\percent}$ in case of the $t\bar{t}$ sample, whereas the amount in case of the $t\bar{t}H$ sample is reduced to $\SI{5.4\pm 2.0}{\percent}$. The expectation is that the events surviving the tau selection contributed from the $t\bar{t}$ sample are mainly due to faked taus. No real $\tau$ lepton is expected here. In contrast, some real $\tau$ contributions from the Higgs sample are expected from the possible decay $H^0\rightarrow \tau^-\tau^+$. This is because the Higgs coupling strength is proportional to the mass and therefore the heaviest lepton (cf. table \ref{SMtable}) is a direct decay product of the Higgs boson. The reduction of the event yield for both signal samples, when at least one tau is selected, is only down to $\SI{45.9\pm 5.5}{\percent}$ (mass point \SI{500}{\giga\electronvolt}) and $\SI{45.3\pm 4.9}{\percent}$ (mass point $\SI{1}{\tera\electronvolt}$). The data event candidates reduce to $\SI{0.6\pm1.7}{\percent}$ in that selection. A more detailed efficiency study is presented in the following section \ref{signalANDbgEff}.
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Signal and background efficiencies}\label{signalANDbgEff}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
Main background contributions are processes from Higgs and top-quark pairs, which contaminate the potential LQ signal. A sub-selection of the event selection (cf. \ref{eventselection}) allows for a more detailed picture of the final states' origin in the $t\tau^{-}+\bar{t}\tau^{+}$ channel for signal and background samples. The MC samples include the so called `truth' information that is the information of the generated particles and the subsequent decay without the digitization and reconstruction step of their signatures in the detector. Taking into account this capability of the MC simulations (cf. \ref{MCsim}) a comparison between the fully reconstructed events and the mere information of the simulated particles is possible. Thus reconstruction efficiencies are the adequate description for that comparison and give insight to the signal and background composition. The efficiency study for signal and background samples are done in two sections \ref{bgeff} and \ref{signaleff} separately.
\subsection{Tau reconstruction efficiencies and fake rates}\label{bgeff}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
The sub-selections for the efficiency study are as follows: The event selection (cf. section \ref{eventselection}) with at least two b-jets and at least one tau lepton is applied additionally to the object selection (cf. section \ref{physObj}), in the following named together as basic selection. Applying the common definition of efficiency, which compares the number of events passing a special selection $Z$ compared to the total event number, different efficiency selections are created. The efficiency definition within the tau reconstruction efficiencies are defined as:
\begin{align}
\epsilon=\frac{\text{\#truth hadronic }\tau\text{ matched to a reconstructed }\tau+Z}{\text{\#truth }\tau\text{ passing the basic selection}+Z}\,\text{.}
\label{effeq}
\end{align}
The applied single selections are summarized in table \ref{ttbarttHEffTable}. The error $\sigma_\epsilon$ is derived from the binomial uncertainty in efficiency following the definition
\begin{align}
\sigma_\epsilon=\frac{\sqrt{n\epsilon\left(1-\epsilon\right)}}{n}\,\text{.}.
\label{binomialerror}
\end{align}
%Tau fake rates
Before discussing the tau reconstruction efficiency a fake rate estimation for the background samples is done at first. The efficiency definition from eq. (\ref{effeq}) can be transferred into the context of fake rates:
\begin{align}
f=\frac{\text{\#truth fake candidate matched to a reconstructed }\tau}{\text{\#truth fake candidate passing the basic selection}}\,\text{.}
\label{fakeeq}
\end{align}
Objects that possibly fake a hadronic tau are given the name `fake candidate'. Four different fake candidates are available, estimating the taus faked by jets in general, by b-jets, by light jets excluding b-jets and electrons and are summarized in table \ref{ttbarttHfakeTable}.
%The total efficiency of reconstructing a hadronic tau matched to its truth particle and passing the basic selection over taking all truth taus fulfilling the basic selections (truth matching for tau) is as low as $\SI{0.635\pm 0.061}{\percent}$ for the top-quark pair events compared to $\SI{29.39\pm 0.26}{\percent}$ efficiency for the $t\bar{t}H$ events. The matching is realized with a geometric matching of the reconstructed and truth particle within a cone of $\Delta R=\sqrt{\Delta\phi^2+\Delta\eta^2}<0.05$. The low efficiency in the $t\bar{t}$ sample reflects the fact that no taus are expected there and some tau leptons are expected from direct Higgs boson decays. The efficiency for Higgs events strongly depends on the $p_T$ of the truth tau up to transverse momenta of $\SI{150}{\giga\electronvolt}$ (see fig. \ref{DividedTau:bg:ttH}). As a conclusion the tau reconstruction efficiency of the $t\bar{t}$ represents a fake tau rate.\par
%
\begin{sidewaystable}[htbp]
\centering
\footnotesize
\renewcommand{\arraystretch}{1.2}
\begin{tabular*}{\linewidth}{@{\extracolsep{\fill}}ccccc}
\hline
\hline
\textbf{sample} & \multicolumn{2}{c}{\textbf{t$\bar{\textbf{t}}$}} & \multicolumn{2}{c}{\textbf{t$\bar{\textbf{t}}$H}}
\\
\hline
\multirow{2}{*}{selection} & reconstructed taus & total number of & reconstructed taus & total number of
\\
& matched to fake object & potential fake objects & matched to fake object & potential fake objects
\\
\hline
%truth matching for tau & $107$ & $16863$ & $8971$ & $30527$
%\\
%efficiency & \multicolumn{2}{c}{$\SI{0.635\pm 0.061}{\percent}$} & \multicolumn{2}{c}{$\SI{29.39\pm 0.26}{\percent}$}
%\\
%\hline
%tau originated from W$^{\pm}$, Z$^0$& $0$ & $0$ & $3251$ & $8843$
%\\
%efficiency & \multicolumn{2}{c}{-} & \multicolumn{2}{c}{$\SI{36.76\pm 0.51}{\percent}$}
%\\
%\hline
%tau originated from W$^{\pm}$, Z$^0$ (1-prong)& $0$ & $0$ & $2488$ & $6195$
%\\
%efficiency & \multicolumn{2}{c}{-} & \multicolumn{2}{c}{$\SI{40.16\pm 0.62}{\percent}$}
%\\
%\hline
%tau originated from W$^{\pm}$, Z$^0$ (3-prong)& $0$ & $0$ & $763$ & $2630$
%\\
%efficiency & \multicolumn{2}{c}{-} & \multicolumn{2}{c}{$\SI{29.01\pm 0.88}{\percent}$}
%\\
%\hline
%tau originated from H$^0$& $0$ & $0$ & $5140$ & $11372$
%\\
%efficiency & \multicolumn{2}{c}{-} & \multicolumn{2}{c}{$\SI{45.20\pm 0.47}{\percent}$}
%\\
%\hline
%tau originated from H$^0$ (1-prong)& $0$ & $0$ & $4132$ & $8390$
%\\
%efficiency & \multicolumn{2}{c}{-} & \multicolumn{2}{c}{$\SI{49.25\pm 0.55}{\percent}$}
%\\
%\hline
%tau originated from H$^0$ (3-prong)& $0$ & $0$ & $1008$ & $2958$
%\\
%efficiency & \multicolumn{2}{c}{-} & \multicolumn{2}{c}{$\SI{34.08\pm 0.87}{\percent}$}
%\\
%\hline
%tau from B-mesons & $107$ & $16860$ & $44$ & $8553$
%\\
%efficiency & \multicolumn{2}{c}{$\SI{0.635\pm 0.061}{\percent}$} & \multicolumn{2}{c}{$\SI{0.514\pm 0.077}{\percent}$}
%\\
%\hline
jet faking as tau & $8232$ & $4728730$ & $21925$ & $2178770$
\\
efficiency & \multicolumn{2}{c}{$\SI{0.1741\pm 0.0019}{\percent}$} & \multicolumn{2}{c}{$\SI{1.0063\pm 0.0068}{\percent}$}
\\
\hline
b-jet faking as tau & $5361$ & $3317260$ & $1993$ & $1332650$
\\
efficiency & \multicolumn{2}{c}{$\SI{0.1616\pm 0.0022}{\percent}$} & \multicolumn{2}{c}{$\SI{0.1496\pm 0.0033}{\percent}$}
\\
\hline
not b-jet faking as tau & $2871$ & $1411470$ & $19932$ & $846121$
\\
efficiency & \multicolumn{2}{c}{$\SI{0.2034\pm 0.0038}{\percent}$} & \multicolumn{2}{c}{$\SI{2.356\pm 0.016}{\percent}$}
\\
\hline
electron faking as tau & $17$ & $1639390$ & $17$ & $488234$
\\
efficiency & \multicolumn{2}{c}{$\SI{0.00104\pm 0.00025}{\percent}$} & \multicolumn{2}{c}{$\SI{0.00348\pm 0.00084}{\percent}$}
\\
\hline
\hline
\end{tabular*}
\caption[Fake rates for the t$\bar{\text{t}}$ and the t$\bar{\text{t}}$H sample.]{Fake rates for different tau lepton categories for the t$\bar{\text{t}}$ and the t$\bar{\text{t}}$H Monte Carlo sample. The fake rate is defined as $f=\frac{\text{\#truth fake candidate matched to a reconstructed }\tau}{\text{\#truth fake candidate passing the basic selection}}$}
\label{ttbarttHfakeTable}
\renewcommand{\arraystretch}{1}
\end{sidewaystable}
%
%
%\begin{figure}
% \centering
% \begin{subfigure}[t]{0.49\textwidth}
% \includegraphics[width=\textwidth]{figures/plots/ttbar/Divided_Tau.pdf}
% \subcaption{Hadronic tau reconstruction efficiency of $\SI{0.635\pm 0.061}{\percent}$ for the top-quark pair events.}
% \label{DividedTau:bg:ttbar}
% \end{subfigure}
% %
% \begin{subfigure}[t]{0.49\textwidth}
% \includegraphics[width=\textwidth]{figures/plots/ttH/Divided_Tau.pdf}
% \subcaption{Hadronic tau reconstruction efficiency of $\SI{29.39\pm 0.26}{\percent}$ for the Higgs events.}
% \label{DividedTau:bg:ttH}
% \end{subfigure}
%\caption[Hadronic tau reconstruction efficiency for the background events.]{Hadronic tau reconstruction efficiency for the background events. The efficiency is defined as number of reconstructed hadronic taus passing the basic selection and is matched to a truth tau over the number of truth taus passing the basic selection (truth matching tau).}
%\label{DividedTau:bg}
%\end{figure}
%
%Jet
The fake rate for jets that is the number of truth jets matched with a reconstructed tau-jet over the number of truth jets satisfying the basic selection is $\SI{0.1741\pm 0.0019}{\percent}$ and $\SI{1.0063\pm 0.0068}{\percent}$ for the top-pair related background and the Higgs related background respectively.\newline
The efficiency depending on the transverse momentum of the truth jet is shown as ratio plot in figure \ref{DividedJet:bg}.
%
\begin{figure}
\centering
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{figures/plots/ttbar/Divided_Jet.pdf}
\subcaption{Estimation of a jet faking a tau lepton of $\SI{0.1741\pm 0.0019}{\percent}$ for the top-quark pair events.}
\label{DividedJet:bg:ttbar}
\end{subfigure}
%
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{figures/plots/ttH/Divided_Jet.pdf}
\subcaption{Estimation of a jet faking a tau lepton of $\SI{1.0063\pm 0.0068}{\percent}$ for the Higgs events.}
\label{DividedJet:bg:ttH}
\end{subfigure}
\caption[Estimation of a jet faking a tau lepton for the background events.]{Estimation of a jet faking a tau lepton for the background events. The efficiency is defined as number of truth jets passing matched to a tau-jet over the number of truth jets passing the basic selection.}
\label{DividedJet:bg}
\end{figure}
%
The qualitative trend is for both background samples the same, showing a peak fake rate at $p_{T,\text{jet}}^{\text{truth}}$ at $\SI{50}{\giga\electronvolt}$. There exists a plateau in the range of the truth jet transverse momentum between $\SI{100}{\giga\electronvolt}$ and $\SI{200}{\giga\electronvolt}$ before the efficiency goes down. As a result the jet fake rate depends directly on the jets' transverse momentum. The main difference is that the fake rate for the $t\bar{t}H$ sample is greater than that of the $t\bar{t}$ sample. These differences in the mean fake rate could be due to the difference in jet composition. For $t\bar{t}$ the jets are gluon-originated (cf. fig. \ref{ttbarproduction}) and for $t\bar{t}H$ quark-originated (cf. fig. \ref{ttHproduction}) jets.\par
%b-jet
The fake rate for b-jets is defined in the same way as for jets above with the difference of selecting only jets of b-flavor. The fake rate values of the two background samples are not in agreement within the uncertainties, Although they are in the same range of $\SI{0.1616\pm 0.0022}{\percent}$ and $\SI{0.1496\pm 0.0033}{\percent}$ corresponding to $t\bar{t}$ and $t\bar{t}H$. The efficiency depending on the truth transverse momentum of the b-jet is shown in figure \ref{DividedBjet:bg}.
%
\begin{figure}
\centering
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{figures/plots/ttbar/Divided_Bjet.pdf}
\subcaption{Estimation of a b-jet faking a tau lepton of $\SI{0.1616\pm 0.0022}{\percent}$ for the top-quark pair events.}
\label{DividedBjet:bg:ttbar}
\end{subfigure}
%
\begin{subfigure}[t]{0.49\textwidth}
\includegraphics[width=\textwidth]{figures/plots/ttH/Divided_Bjet.pdf}
\subcaption{Estimation of a b-jet faking a tau lepton of $\SI{0.1496\pm 0.0033}{\percent}$ for the Higgs events.}
\label{DividedBjet:bg:ttH}