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

Ce papier présente une mesure des huit coefficients de polarisation angulaire ($A_0$ à $A_7$) pour la production Drell-Yan $Z \to \mu^+\mu^-$ dans des collisions proton-proton à 13 TeV utilisant 140 fb$^{-1}$ de données CMS, fournissant des résultats doublement différentiels en impulsion transverse et en rapidité qui sont comparés aux prédictions QCD d'ordre suivant-le-suivant-le-plus-élevé.

L'essentiel

La Vue d'Ensemble : Une Machine à Pinceau Cosmique

Imaginez le Grand Collisionneur de Hadrons (LHC) au CERN comme la machine à pinceau la plus puissante du monde. Les scientifiques font entrer en collision des protons à une vitesse proche de celle de la lumière. Habituellement, ces collisions créent un chaos désordonné de particules. Cependant, parfois, la collision génère une particule lourde et instable appelée boson Z, qui se désintègre immédiatement en deux muons (des cousins lourds des électrons).

Ce papier concerne CMS, l'un des détecteurs géants qui surveillent cette machine à pinceau. L'équipe n'a pas seulement compté le nombre de fois où cela s'est produit ; elle a voulu comprendre comment les muons s'envolaient. Sortaient-ils tout droit ? Tournaient-ils sur eux-mêmes ? Privilégiaient-ils une direction par rapport à une autre ?

L'Objectif : Mesurer la « Rotation » de la Collision

Les scientifiques ont mesuré huit nombres différents (étiquetés de $A_0$ à $A_7$). Imaginez ces nombres comme un bulletin de notes détaillé sur la « posture » ou la polarisation du boson Z avant qu'il n'explose.

• L'Analogie : Imaginez un feu d'artifice qui explose dans le ciel. S'il explose parfaitement symétriquement, les étincelles s'éparpillent en une sphère parfaite. S'il est incliné ou en rotation, les étincelles peuvent s'envoler davantage vers la gauche, davantage vers le haut, ou en spirale.
• La Mesure : Les huit coefficients ($A_0$ à $A_7$) nous disent exactement à quoi ressemble cette « forme » de l'explosion. Ils révèlent si le boson Z tournait, vacillait, ou s'il était « étiré » dans une certaine direction.

Comment Ils Ont Fait : Le « Double-Vérification »

L'équipe a examiné 140 billions de collisions (140 fb⁻¹ de données) enregistrées entre 2016 et 2018. Ils n'ont pas seulement regardé le tas de données dans son ensemble ; ils l'ont tranché comme une miche de pain pour voir si la « rotation » changeait en fonction de la force avec laquelle les protons se heurtaient.

1. Vitesse (Impulsion Transverse) : Ils ont observé des muons se déplaçant lentement sur le côté par rapport à ceux se déplaçant très vite.
2. Angle (Rapidité) : Ils ont observé des muons volant tout droit par rapport à ceux volant à un angle aigu.

En mesurant les angles des muons dans ces tranches spécifiques, ils ont pu calculer les huit coefficients avec une extrême précision.

Les Règles du Jeu : La Règle « Lam-Tung »

Le papier discute d'une règle célèbre en physique appelée la relation Lam-Tung.

• L'Analogie : Imaginez une règle qui dit : « Si vous lancez une balle droit vers le haut, elle doit retomber droit. » Dans le monde de la physique des particules, au niveau le plus simple du calcul, deux de ces coefficients ($A_0$ et $A_2$) devraient s'annuler parfaitement ($A_0 - A_2 = 0$).
• La Réalité : Le papier confirme que cette règle tient bien à basse vitesse, mais à mesure que les collisions deviennent plus énergétiques (impulsion plus élevée), la règle commence à se briser. Ce n'est pas un échec ; c'est une caractéristique ! Cela nous indique que les parties « désordonnées » de la collision (comme des particules supplémentaires étant éjectées) commencent à compter.

Les Résultats : Données vs Théorie

Les scientifiques ont comparé leurs mesures aux meilleures simulations informatiques disponibles (les « prédictions théoriques »).

• La Bonne Nouvelle : Pour la plupart des coefficients, les données du monde réel correspondaient très bien aux modèles informatiques. Cela signifie que notre compréhension actuelle de la façon dont ces particules interagissent est solide.
• La Tension Intéressante : Dans la gamme intermédiaire de vitesses, les données pour un coefficient spécifique ($A_0$) étaient légèrement plus élevées que ce que l'ordinateur avait prédit (environ 3 écarts-types d'écart). C'est comme si une prévision météorologique prévoyait 50 % de chances de pluie, mais qu'il a effectivement plu 80 % du temps. Ce n'est pas une catastrophe, mais cela suggère que le modèle informatique pourrait manquer un tout petit détail.
• Les Coefficients « Fantômes » : Trois des coefficients ($A_5, A_6, A_7$) sont censés être nuls ou très proches de zéro. Les données ont montré qu'ils étaient en effet minuscules, cohérents avec zéro, bien que l'un d'eux ($A_6$) ait montré un tout petit, faible indice d'être non nul. C'est comme entendre un chuchotement dans une pièce calme ; il est là, mais vous avez besoin d'oreilles très sensibles pour l'entendre.

Pourquoi Cela Importe

Ce papier est essentiellement un contrôle d'étalonnage de haute précision pour les lois de la physique.

1. Comprendre la « Colle » : Ces mesures nous aident à comprendre la « dynamique partonique » — comment les tout petits blocs de construction à l'intérieur du proton (quarks et gluons) se comportent lorsqu'ils entrent en collision.
2. Tester la Théorie : En comparant la « rotation » du boson Z avec des mathématiques complexes (Chromodynamique Quantique), les scientifiques soumettent notre compréhension de l'univers à un test de stress. Si les mathématiques ne correspondent pas aux données, cela signifie que nous devons inventer une nouvelle physique.
3. La Référence : Ce papier fournit une nouvelle « règle » ultra-précise pour les expériences futures. Toute nouvelle théorie doit pouvoir expliquer ces huit nombres.

Résumé

En bref, l'équipe CMS a pris une immense photo des collisions de particules, mesuré les angles exacts des particules résultantes et calculé huit nombres décrivant la « rotation » de l'événement. Ils ont constaté que, bien que nos théories actuelles soient majoritairement correctes, il existe de minuscules et fascinantes discordances dans la gamme de vitesses intermédiaires qui tiennent les physiciens en éveil, assurant que la recherche d'une compréhension plus profonde de l'univers se poursuit.

Résumé technique

1. Problème et Motivation

Le processus Drell–Yan (DY) ($pp \to \gamma^*/Z \to \ell^+\ell^-$) est une sonde fondamentale pour comprendre la dynamique partonique sous-jacente des collisions proton-proton et les mécanismes de production des bosons vecteurs électrofaibles. La distribution angulaire des leptons de l'état final dans le référentiel au repos du boson encode des informations sur les états de polarisation du photon virtuel ou du boson Z.

Cette distribution est paramétrée par huit coefficients de polarisation angulaire, $A_0$ à $A_7$, définis dans le référentiel de Collins–Soper (CS).

• Objectifs physiques :
  • $A_0$ et $A_2$ : Paramètrent la polarisation longitudinale et l'interférence entre les états transverses. Leur différence ($A_0 - A_2$) teste la relation de Lam–Tung ($A_0 - A_2 = 0$), qui est valable à l'ordre dominant (LO) en QCD mais devrait être violée par le rayonnement QCD d'ordre supérieur et l'impulsion transverse intrinsèque des partons.
  • $A_4$ : Liée à l'asymétrie avant-arrière ($A_{FB}$) et utilisée pour mesurer l'angle de mélange faible effectif ($\sin^2 \theta_{eff}^W$).
  • $A_1, A_3$ : Encodent l'interférence entre les états longitudinaux et transverses et dépendent des couplages axiaux/vectoriels.
  • $A_5, A_6, A_7$ : Définissent les asymétries T-paires. Celles-ci sont attendues nulles aux ordres LO et NLO mais peuvent être non nulles à l'ordre Next-to-Next-to-Leading Order (NNLO) en raison de processus à une boucle.

Bien que des mesures antérieures existent à $\sqrt{s}=8$ TeV et dans la région avant (LHCb), une mesure doublement différentielle complète de l'ensemble des coefficients ($A_0$–$A_7$) à $\sqrt{s}=13$ TeV avec de hautes statistiques était nécessaire pour contraindre les modèles théoriques et tester les prédictions QCD NNLO avec une plus grande précision.

2. Méthodologie

Données et Sélection d'événements

• Jeu de données : Données de collisions proton-proton collectées par le détecteur CMS à $\sqrt{s} = 13$ TeV durant 2016–2018, correspondant à une luminosité intégrée de 140 fb$^{-1}$.
• Signal : Événements dilepton ($Z \to \mu^+\mu^-$) avec une masse invariante $81 < m_{\mu\mu} < 101$ GeV (région du pôle Z).
• Cuts cinématiques :
  • Impulsion transverse du dimuon ($p_T^{\mu\mu}$) : 0 à 400 GeV.
  • Rapidité du dimuon ($|y_{\mu\mu}|$) : $< 2.4$.
  • Critères d'isolement et d'identification des muons appliqués pour assurer une haute pureté.
  • Exigences de déclenchement : Déclenchements à muon unique avec $p_T > 24$ ou $27$ GeV selon la période de prise de données.

Simulation et Corrections

• Simulation du signal : Générée en utilisant POWHEG v2.0 + MiNNLOPS couplé à PYTHIA 8 (pour le showering de partons/hadronisation) et PHOTOS++ (pour le rayonnement QED de l'état final). Le jeu de PDF NNPDF3.1 a été utilisé.
• Estimation du fond : Simulations Monte Carlo pour les processus de dibosons (WW, WZ, ZZ), $t\bar{t}$ et $W$+jets. Les contributions de fond sont faibles (0,03 % à faible $p_T$ augmentant jusqu'à ~2,9 % à haut $p_T$).
• Corrections :
  • Facteurs d'échelle : Déduits via la méthode « tag-and-probe » pour corriger les efficacités de déclenchement, d'identification et d'isolement.
  • Repondération du pileup : Pour faire correspondre la distribution simulée du pileup aux données.
  • Repondération cinématique : Les événements simulés ont été repondérés dans les bins $(p_T^{\mu\mu}, y_{\mu\mu})$ pour correspondre aux distributions cinématiques des données, minimisant ainsi la dépendance au modèle.
  • Correction de pré-déclenchement (Prefiring) : A traité le désalignement des croisements de paquets dans le déclencheur de niveau 1.

Extraction des coefficients

L'extraction utilise une méthode de modèles (template method) pour tenir compte de l'acceptance du détecteur, de la résolution et de la migration entre les bins de $p_T$.

1. Définition : La section efficace différentielle est développée comme une somme de polynômes trigonométriques $f_i(\theta^*, \phi^*)$ pondérés par des coefficients $P_i$ (où $A_i = P_i/P_8$).
2. Modèles (Templates) : Neuf modèles 2D ($T_i$) ont été construits à partir d'événements MC reconstruits. Chaque modèle correspond à un terme angulaire spécifique $f_i$, pondéré pour éliminer l'information de polarisation de la génération d'événements tout en préservant les effets du détecteur.
3. Ajustement : Un ajustement par vraisemblance maximale binné en deux dimensions (utilisant MINUIT) a été effectué sur les distributions angulaires observées ($\cos \theta^*, \phi^*$) dans chacun des 16 bins cinématiques (8 bins $p_T$ × 2 bins de rapidité). L'ajustement extrait les coefficients $P_i$ en minimisant la différence entre les données et la combinaison linéaire des modèles.

3. Contributions clés

• Premier ensemble complet à 13 TeV : Il s'agit de la première mesure de l'ensemble complet des coefficients angulaires ($A_0$–$A_7$) dans la région de rapidité centrale ($|y_{\mu\mu}| < 2.4$) à $\sqrt{s}=13$ TeV.
• Précision doublement différentielle : Les résultats sont fournis dans 8 bins de $p_T^{\mu\mu}$ et 2 bins de $|y_{\mu\mu}|$, offrant une vue multidimensionnelle de la structure angulaire.
• Comparaison NNLO : Les résultats sont comparés aux prédictions théoriques de pointe au NNLO en QCD (utilisant POWHEG+MiNNLOPS et MADGRAPH5_aMC@NLO).
• Contraintes sur les asymétries T-paires : Fournit les contraintes les plus strictes à ce jour sur les coefficients T-paires ($A_5, A_6, A_7$) dans la région centrale.

4. Résultats

Accord général

• Les coefficients mesurés sont généralement en accord avec les prédictions QCD NNLO dans les limites des incertitudes.
• Le grand jeu de données (140 fb$^{-1}$) a considérablement réduit les incertitudes statistiques par rapport aux mesures précédentes du Run 1.

Observations spécifiques

• $A_0$ et $A_2$ :
  • $A_0$ montre un écart d'environ 3 écarts-types dans la gamme intermédiaire de $p_T$ (10–35 GeV) entre les données et la prédiction MADGRAPH5_aMC@NLO, tandis que POWHEG+MiNNLOPS s'accorde bien.
  • $A_2$ est légèrement surestimée par les simulations MC à haut $p_T$ ($>50$ GeV).
• Relation de Lam–Tung ($A_0 - A_2$) :
  • La relation est violée, comme attendu en raison des effets QCD d'ordre supérieur.
  • Les prédictions POWHEG+MiNNLOPS s'accordent bien avec les valeurs mesurées de $A_0 - A_2$ dans la gamme 10–35 GeV.
  • Les deux prédictions MC sous-estiment la mesure pour $p_T > 35$ GeV.
• $A_1$ :
  • Significativement surestimée par POWHEG+MiNNLOPS dans la région $1 < |y_{\mu\mu}| < 2.4$ pour $p_T < 35$ GeV.
  • MADGRAPH5_aMC@NLO surestime également $A_1$ à faible $p_T$ dans les deux régions de rapidité.
• Coefficients T-paires ($A_5, A_6, A_7$) :
  • Les valeurs sont de l'ordre de $10^{-3}$ et compatibles avec zéro dans la plupart des bins.
  • La déviation la plus significative concerne $A_6$ dans le bin $55 < p_T < 80$ GeV et $1.2 < |y_{\mu\mu}| < 2.4$, montrant un écart de 2,8$\sigma$ par rapport à zéro. Ceci est compatible avec les attentes NNLO mais de magnitude plus faible que les observations précédentes d'ATLAS à 8 TeV.
• Incertitudes systématiques :
  • Dominantes pour $A_0$ et $A_1$ à faible $p_T$ ($<55$ GeV) et pour $A_2, A_3, A_4$ jusqu'à $p_T < 200$ GeV.
  • Les incertitudes statistiques dominent pour les autres coefficients et les régions de $p_T$ plus élevées.

5. Signification

• Validation théorique : Les résultats fournissent une référence rigoureuse pour les calculs QCD NNLO, confirmant que les corrections d'ordre supérieur sont essentielles pour décrire la structure angulaire de la production Drell–Yan.
• Dynamique des partons : Les mesures offrent de nouvelles contraintes sur les Fonctions de Distribution de Partons (PDF) et la modélisation de l'impulsion transverse intrinsèque des partons.
• Physique future : La précision et le maillage multidimensionnel établis ici servent de référence critique pour les futures études phénoménologiques et pour tester des calculs d'ordre encore supérieur (par exemple, N$^3$LO).
• Performance du détecteur : Démontre la capacité du détecteur CMS à mesurer des corrélations angulaires subtiles avec une grande précision, même en présence d'effets complexes du détecteur et de pileup.

En résumé, cet article représente une étape importante en physique électrofaible de précision, confirmant la validité des descriptions QCD NNLO du mécanisme de production du boson Z tout en mettant en évidence des régions cinématiques spécifiques où les modèles théoriques pourraient nécessiter un raffinement supplémentaire.