Measurement of the Z $\to$ $\mu^+\mu^-$ angular coefficients in pp collisions at $\sqrt{s}$ = 13 TeV as functions of transverse momentum and rapidity

Diese Arbeit präsentiert eine Messung der acht Winkel-Polarisationskoeffizienten ($A_0$ bis $A_7$) für die Drell-Yan-Produktion $Z \to \mu^+\mu^-$ in Proton-Proton-Kollisionen bei 13 TeV unter Verwendung von 140 fb$^{-1}$ CMS-Daten und liefert doppelt-differenzierte Ergebnisse in transversalem Impuls und Rapidität, die mit Vorhersagen der Quantenchromodynamik im nächsten-zu-nächsten führenden Ordnung verglichen werden.

Das Wesentliche

Das große Ganze: Eine kosmische Flippermaschine

Stellen Sie sich den Large Hadron Collider (LHC) am CERN als die leistungsstärkste Flippermaschine der Welt vor. Wissenschaftler lassen Protonen mit nahezu Lichtgeschwindigkeit aufeinanderprallen. Normalerweise erzeugen diese Kollisionen ein chaotisches Durcheinander von Teilchen. Manchmal entsteht jedoch bei der Kollision ein schweres, instabiles Teilchen, das als Z-Boson bezeichnet wird, welches sich sofort in zwei Myonen (schwere Cousins der Elektronen) aufspaltet.

Dieses Paper handelt von CMS, einem der riesigen Detektoren, die diese Flippermaschine beobachten. Das Team zählte nicht nur, wie oft dies geschah; sie wollten verstehen, wie die Myonen herausflogen. Schossen sie geradeaus? Drehten sie sich? Bevorzugten sie eine Richtung gegenüber einer anderen?

Das Ziel: Messung des „Spins" der Kollision

Die Wissenschaftler maßen acht verschiedene Zahlen (kennzeichnet mit $A_0$ bis $A_7$). Betrachten Sie diese Zahlen als detaillierten Zeugnisbericht über die „Haltung" oder Polarisation des Z-Bosons, bevor es explodierte.

• Die Analogie: Stellen Sie sich ein Feuerwerk vor, das am Himmel explodiert. Wenn es perfekt symmetrisch explodiert, fliegen die Funken in einer perfekten Kugelkugel heraus. Wenn es geneigt ist oder sich dreht, könnten die Funken mehr nach links, mehr nach oben oder spiralförmig herausgeschleudert werden.
• Die Messung: Die acht Koeffizienten ($A_0$ bis $A_7$) verraten uns genau, wie diese „Form" der Explosion aussieht. Sie enthüllen, ob das Z-Boson sich drehte, wackelte oder in eine bestimmte Richtung „gestreckt" war.

Wie sie es taten: Der „Doppel-Check"

Das Team betrachtete 140 Billionen Kollisionen (140 fb⁻¹ an Daten), die zwischen 2016 und 2018 aufgezeichnet wurden. Sie schauten nicht nur auf den gesamten Datenhaufen; sie schnitten ihn wie ein Brotlaib auf, um zu sehen, ob sich der „Spin" änderte, je nachdem, wie hart die Protonen aufeinanderprallten.

1. Geschwindigkeit (Transversaler Impuls): Sie betrachteten Myonen, die sich langsam seitwärts bewegten, im Vergleich zu denen, die sich sehr schnell bewegten.
2. Winkel (Rapidität): Sie betrachteten Myonen, die geradeaus flogen, im Vergleich zu denen, die in einem spitzen Winkel flogen.

Indem sie die Winkel der Myonen in diesen spezifischen Schnitten maßen, konnten sie die acht Koeffizienten mit extremer Präzision berechnen.

Die Regeln des Spiels: Die „Lam-Tung"-Regel

Das Paper diskutiert eine berühmte Regel in der Physik, die Lam-Tung-Beziehung.

• Die Analogie: Denken Sie an eine Regel, die besagt: „Wenn Sie einen Ball gerade nach oben werfen, muss er gerade nach unten kommen." In der Welt der Teilchenphysik sollten sich auf der einfachsten Berechnungsebene zwei dieser Koeffizienten ($A_0$ und $A_2$) perfekt gegenseitig aufheben ($A_0 - A_2 = 0$).
• Die Realität: Das Paper bestätigt, dass diese Regel bei niedrigen Geschwindigkeiten gut hält, aber sobald die Kollisionen energiereicher werden (höherer Impuls), beginnt die Regel zu versagen. Dies ist kein Versagen, sondern ein Merkmal! Es zeigt uns, dass die „chaotischen" Teile der Kollision (wie das Herausschleudern zusätzlicher Teilchen) beginnen, eine Rolle zu spielen.

Die Ergebnisse: Daten vs. Theorie

Die Wissenschaftler verglichen ihre Messungen mit den besten verfügbaren Computersimulationen (den „theoretischen Vorhersagen").

• Die gute Nachricht: Für die meisten Koeffizienten stimmten die Realwelt-Daten sehr gut mit den Computermodellen überein. Dies bedeutet, dass unser aktuelles Verständnis davon, wie diese Teilchen interagieren, solide ist.
• Die interessante Spannung: Im mittleren Geschwindigkeitsbereich war die Daten für einen spezifischen Koeffizienten ($A_0$) leicht höher als vom Computer vorhergesagt (etwa 3 Standardabweichungen abweichend). Es ist, als würde eine Wettervorhersage eine 50-prozentige Regenwahrscheinlichkeit vorhersagen, es aber tatsächlich zu 80 Prozent geregnet hätte. Es ist keine Katastrophe, aber es deutet darauf hin, dass dem Computermodell vielleicht ein winziges Detail fehlt.
• Die „Geister"-Koeffizienten: Drei der Koeffizienten ($A_5, A_6, A_7$) sollten null oder sehr nahe daran sein. Die Daten zeigten, dass sie tatsächlich winzig waren und mit Null konsistent waren, obwohl einer von ihnen ($A_6$) einen winzigen, schwachen Hinweis darauf zeigte, nicht null zu sein. Dies ist wie das Hören eines Flüsterns in einem ruhigen Raum; es ist da, aber man braucht sehr empfindliche Ohren, um es zu hören.

Warum dies wichtig ist

Dieses Paper ist im Wesentlichen eine hochpräzise Kalibrierungsprüfung für die Gesetze der Physik.

1. Verstehen des „Klebstoffs": Diese Messungen helfen uns, die „partonischen Dynamiken" zu verstehen – wie sich die winzigen Bausteine innerhalb des Protons (Quarks und Gluonen) verhalten, wenn sie kollidieren.
2. Testen der Theorie: Indem sie den „Spin" des Z-Bosons mit komplexer Mathematik (Quantenchromodynamik) vergleichen, testen die Wissenschaftler unser Verständnis des Universums auf Belastbarkeit. Wenn die Mathematik nicht mit den Daten übereinstimmt, bedeutet dies, dass wir neue Physik erfinden müssen.
3. Der Maßstab: Dieses Paper liefert einen neuen, ultra-präzisen „Lineal" für zukünftige Experimente. Jede neue Theorie muss in der Lage sein, diese acht Zahlen zu erklären.

Zusammenfassung

Kurz gesagt, machte das CMS-Team einen massiven Schnappschuss von Teilchenkollisionen, maß die genauen Winkel der resultierenden Teilchen und berechnete acht Zahlen, die den „Spin" des Ereignisses beschreiben. Sie fanden heraus, dass unsere aktuellen Theorien zwar größtenteils korrekt sind, aber es winzige, faszinierende Diskrepanzen im mittleren Geschwindigkeitsbereich gibt, die Physiker auf Trab halten und sicherstellen, dass die Suche nach einem tieferen Verständnis des Universums weitergeht.

Technische Zusammenfassung

1. Problem und Motivation

Der Drell-Yan-Prozess (DY) ($pp \to \gamma^*/Z \to \ell^+\ell^-$) ist eine fundamentale Sonde zum Verständnis der zugrunde liegenden partonischen Dynamik von Proton-Proton-Kollisionen und der Produktionsmechanismen elektroschwacher Vektorbosonen. Die Winkelverteilung der Endzustandsleptonen im Ruhesystem des Bosons kodiert Informationen über die Polarisationszustände des virtuellen Photons oder des Z-Bosons.

Diese Verteilung wird durch acht Winkel-Polarisationskoeffizienten, $A_0$ bis $A_7$, im Collins-Soper (CS)-Rahmen parametrisiert.

• Physikalische Ziele:
  • $A_0$ und $A_2$: Parametrisieren die longitudinale Polarisation und die Interferenz zwischen transversalen Zuständen. Ihre Differenz ($A_0 - A_2$) testet die Lam-Tung-Relation ($A_0 - A_2 = 0$), die in der Leading Order (LO) der QCD gilt, aber durch QCD-Strahlung höherer Ordnung und den intrinsischen transversalen Impuls der Partonen verletzt werden sollte.
  • $A_4$: Steht im Zusammenhang mit der Vorwärts-Rückwärts-Asymmetrie ($A_{FB}$) und wird zur Messung des effektiven schwachen Mischungswinkels ($\sin^2 \theta_{eff}^W$) verwendet.
  • $A_1, A_3$: Kodieren die Interferenz zwischen longitudinalen und transversalen Zuständen und hängen von axialen/vektoriellen Kopplungen ab.
  • $A_5, A_6, A_7$: Definieren T-ungerade Asymmetrien. Diese werden in LO und NLO als null erwartet, können jedoch aufgrund von Ein-Schleifen-Prozessen in der Next-to-Next-to-Leading Order (NNLO) von null abweichen.

Während frühere Messungen bei $\sqrt{s}=8$ TeV und im Vorwärtsbereich (LHCb) existierten, war eine umfassende, doppelt differentielle Messung des vollständigen Satzes von Koeffizienten ($A_0$–$A_7$) bei $\sqrt{s}=13$ TeV mit hoher Statistik erforderlich, um theoretische Modelle einzuschränken und NNLO-QCD-Vorhersagen mit größerer Präzision zu testen.

2. Methodik

Daten und Ereignisauswahl

• Datensatz: Proton-Proton-Kollisionsdaten, die vom CMS-Detektor bei $\sqrt{s} = 13$ TeV während 2016–2018 gesammelt wurden, entsprechend einer integrierten Luminosität von 140 fb$^{-1}$.
• Signal: Dilepton-Ereignisse ($Z \to \mu^+\mu^-$) mit einer invariante Masse $81 < m_{\mu\mu} < 101$ GeV (Z-Pol-Bereich).
• Kinematische Schnitte:
  • Transverser Impuls des Dimuonsystems ($p_T^{\mu\mu}$): 0 bis 400 GeV.
  • Rapidität des Dimuonsystems ($|y_{\mu\mu}|$): $< 2.4$.
  • Isolations- und Identifikationskriterien für Myonen angewendet, um eine hohe Reinheit sicherzustellen.
  • Trigger-Anforderungen: Einzelmyon-Triggers mit $p_T > 24$ oder $27$ GeV, je nach Datenerfassungsperiode.

Simulation und Korrekturen

• Signalsimulation: Generiert mit POWHEG v2.0 + MiNNLOPS, gekoppelt mit PYTHIA 8 (für Parton-Showering/Hadronisierung) und PHOTOS++ (für QED-Strahlung im Endzustand). Der NNPDF3.1-PDF-Satz wurde verwendet.
• Schätzung des Untergrunds: Monte-Carlo-Simulationen für Diboson-Prozesse (WW, WZ, ZZ), $t\bar{t}$ und $W$+Jets. Die Beiträge des Untergrunds sind gering (0,03 % bei niedrigem $p_T$, steigend auf ~2,9 % bei hohem $p_T$).
• Korrekturen:
  • Skalierungsfaktoren: Abgeleitet mittels „Tag-and-Probe", um Trigger-, Identifikations- und Isolations-Effizienzen zu korrigieren.
  • Pileup-Neugewichtung: Um die simulierte Pileup-Verteilung an die Daten anzupassen.
  • Kinematische Neugewichtung: Simulierte Ereignisse wurden in $(p_T^{\mu\mu}, y_{\mu\mu})$-Bins neu gewichtet, um die kinematischen Verteilungen der Daten anzupassen und die Modellabhängigkeit zu minimieren.
  • Prefiring-Korrektur: Adressierte Fehlausrichtung von Bunch-Crossings im Level-1-Trigger.

Extraktion der Koeffizienten

Die Extraktion nutzt eine Template-Methode, um Detektorakzeptanz, Auflösung und Migration zwischen $p_T$-Bins zu berücksichtigen.

1. Definition: Der differentielle Wirkungsquerschnitt wird als Summe trigonometrischer Polynome $f_i(\theta^*, \phi^*)$ entwickelt, gewichtet mit Koeffizienten $P_i$ (wobei $A_i = P_i/P_8$).
2. Templates: Neun 2D-Templates ($T_i$) wurden aus rekonstruierten MC-Ereignissen konstruiert. Jedes Template entspricht einem spezifischen Winkelterm $f_i$, gewichtet, um Polarisationsinformationen aus der Ereignisgenerierung zu entfernen, während Detektoreffekte erhalten bleiben.
3. Fitting: Ein zweidimensionales geblocktes Maximum-Likelihood-Fit (unter Verwendung von MINUIT) wurde auf die beobachteten Winkelverteilungen ($\cos \theta^*, \phi^*$) in jedem der 16 kinematischen Bins durchgeführt (8 $p_T$-Bins $\times$ 2 Rapiditäts-Bins). Das Fit extrahiert die Koeffizienten $P_i$ durch Minimierung der Differenz zwischen Daten und der linearen Kombination der Templates.

3. Hauptbeiträge

• Erster vollständiger Satz bei 13 TeV: Dies ist die erste Messung des vollständigen Satzes von Winkelkoeffizienten ($A_0$–$A_7$) im zentralen Rapiditätsbereich ($|y_{\mu\mu}| < 2.4$) bei $\sqrt{s}=13$ TeV.
• Doppelt differentielle Präzision: Die Ergebnisse werden in 8 Bins von $p_T^{\mu\mu}$ und 2 Bins von $|y_{\mu\mu}|$ bereitgestellt und bieten eine mehrdimensionale Sicht auf die Winkelstruktur.
• NNLO-Vergleich: Die Ergebnisse werden mit hochmodernen theoretischen Vorhersagen bei NNLO in QCD (unter Verwendung von POWHEG+MiNNLOPS und MADGRAPH5_aMC@NLO) verglichen.
• Einschränkungen T-ungerader Asymmetrien: Liefert die bisher strengsten Einschränkungen für die T-ungeraden Koeffizienten ($A_5, A_6, A_7$) im zentralen Bereich.

4. Ergebnisse

Allgemeine Übereinstimmung

• Die gemessenen Koeffizienten stimmen im Allgemeinen innerhalb der Unsicherheiten mit NNLO-QCD-Vorhersagen überein.
• Der große Datensatz (140 fb$^{-1}$) reduzierte die statistischen Unsicherheiten im Vergleich zu früheren Run-1-Messungen erheblich.

Spezifische Beobachtungen

• $A_0$ und $A_2$:
  • $A_0$ zeigt eine Abweichung von etwa 3 Standardabweichungen im intermediären $p_T$-Bereich (10–35 GeV) zwischen den Daten und der MADGRAPH5_aMC@NLO-Vorhersage, während POWHEG+MiNNLOPS gut übereinstimmt.
  • $A_2$ wird von MC-Simulationen bei hohem $p_T$ ($>50$ GeV) leicht überschätzt.
• Lam-Tung-Relation ($A_0 - A_2$):
  • Die Relation wird verletzt, wie aufgrund von QCD-Effekten höherer Ordnung erwartet.
  • POWHEG+MiNNLOPS-Vorhersagen stimmen im Bereich 10–35 GeV gut mit den gemessenen $A_0 - A_2$-Werten überein.
  • Beide MC-Vorhersagen unterschätzen die Messung für $p_T > 35$ GeV.
• $A_1$:
  • Wird von POWHEG+MiNNLOPS im Bereich $1 < |y_{\mu\mu}| < 2.4$ für $p_T < 35$ GeV signifikant überschätzt.
  • MADGRAPH5_aMC@NLO überschätzt $A_1$ bei niedrigem $p_T$ in beiden Rapiditätsbereichen ebenfalls.
• T-ungerade Koeffizienten ($A_5, A_6, A_7$):
  • Die Werte liegen in der Größenordnung von $10^{-3}$ und sind in den meisten Bins mit null konsistent.
  • Die signifikanteste Abweichung tritt für $A_6$ im Bin $55 < p_T < 80$ GeV und $1.2 < |y_{\mu\mu}| < 2.4$ auf und zeigt eine 2,8$\sigma$-Abweichung von null. Dies ist konsistent mit NNLO-Erwartungen, aber in der Magnitude kleiner als frühere ATLAS-Beobachtungen bei 8 TeV.
• Systematische Unsicherheiten:
  • Dominant für $A_0$ und $A_1$ bei niedrigem $p_T$ ($<55$ GeV) sowie für $A_2, A_3, A_4$ bis $p_T < 200$ GeV.
  • Statistische Unsicherheiten dominieren für andere Koeffizienten und Bereiche mit höherem $p_T$.

5. Bedeutung

• Theoretische Validierung: Die Ergebnisse liefern einen rigorosen Benchmark für NNLO-QCD-Rechnungen und bestätigen, dass Korrekturen höherer Ordnung entscheidend sind, um die Winkelstruktur der Drell-Yan-Produktion zu beschreiben.
• Partondynamik: Die Messungen bieten neue Einschränkungen für Parton-Verteilungsfunktionen (PDFs) und die Modellierung des intrinsischen transversalen Impulses von Partonen.
• Zukünftige Physik: Die hier etablierte Präzision und mehrdimensionale Binnung dienen als kritische Referenz für zukünftige phänomenologische Studien und zum Testen noch höherer Ordnungen (z. B. N$^3$LO).
• Detektorleistung: Demonstriert die Fähigkeit des CMS-Detektors, subtile Winkelkorrelationen mit hoher Präzision zu messen, selbst in Gegenwart komplexer Detektoreffekte und Pileup.

Zusammenfassend stellt dieses Papier einen Meilenstein in der Präzisionsphysik der elektroschwachen Wechselwirkung dar, bestätigt die Gültigkeit von NNLO-QCD-Beschreibungen des Z-Boson-Produktionsmechanismus und hebt spezifische kinematische Bereiche hervor, in denen theoretische Modelle einer weiteren Verfeinerung bedürfen.