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

Questo lavoro presenta una misura degli otto coefficienti di polarizzazione angolare ($A_0$–$A_7$) per la produzione Drell-Yan $Z \to \mu^+\mu^-$ in collisioni protone-protone a 13 TeV, utilizzando 140 fb$^{-1}$ di dati CMS, fornendo risultati doppiamente differenziali in impulso trasverso e rapidità che sono confrontati con previsioni QCD al secondo ordine successivo rispetto all'ordine principale.

L'essenziale

Il quadro generale: Una macchina da flipper cosmica

Immaginate il Large Hadron Collider (LHC) al CERN come la macchina da flipper più potente al mondo. Gli scienziati fanno scontrare protoni a velocità prossime a quella della luce. Di solito, queste collisioni creano un caos disordinato di particelle. Tuttavia, a volte, la collisione genera una particella pesante e instabile chiamata bosone Z, che si spezza immediatamente in due muoni (cugini pesanti degli elettroni).

Questo documento riguarda CMS, uno dei giganteschi rivelatori che osservano questa macchina da flipper. Il team non si è limitato a contare quante volte ciò accadeva; voleva comprendere come i muoni venivano espulsi. Sparavano dritti in avanti? Ruotavano? Favorivano una direzione rispetto a un'altra?

L'obiettivo: Misurare lo "spin" della collisione

Gli scienziati hanno misurato otto numeri diversi (etichettati da $A_0$ a $A_7$). Pensate a questi numeri come a una pagella dettagliata sulla "postura" o polarizzazione del bosone Z prima che esplodesse.

• L'analogia: Immaginate un fuoco d'artificio che esplode nel cielo. Se esplode perfettamente simmetricamente, le scintille si disperdono in una sfera perfetta. Se è inclinato o sta ruotando, le scintille potrebbero sparire più verso sinistra, più verso l'alto, o a spirale.
• La misurazione: Gli otto coefficienti ($A_0$ fino a $A_7$) ci dicono esattamente come appare quella "forma" dell'esplosione. Rivelano se il bosone Z stava ruotando, oscillando, o se era "allungato" in una certa direzione.

Come l'hanno fatto: Il "doppio controllo"

Il team ha esaminato 140 trilioni di collisioni (140 fb⁻¹ di dati) registrati tra il 2016 e il 2018. Non hanno guardato semplicemente l'intero mucchio di dati; lo hanno tagliato a fette come un pane per vedere se lo "spin" cambiava a seconda di quanto duramente i protoni si colpivano a vicenda.

1. Velocità (Momento trasverso): Hanno osservato muoni che si muovevano lentamente lateralmente rispetto a quelli che si muovevano molto velocemente.
2. Angolo (Rapidità): Hanno osservato muoni che volavano dritti in avanti rispetto a quelli che volavano con un angolo acuto.

Misurando gli angoli dei muoni in queste specifiche fette, sono riusciti a calcolare gli otto coefficienti con estrema precisione.

Le regole del gioco: La regola "Lam-Tung"

Il documento discute una famosa regola in fisica chiamata relazione Lam-Tung.

• L'analogia: Pensate a una regola che dice: "Se lanci una palla dritta verso l'alto, deve tornare dritta verso il basso". Nel mondo della fisica delle particelle, al livello più semplice di calcolo, due di questi coefficienti ($A_0$ e $A_2$) dovrebbero annullarsi perfettamente a vicenda ($A_0 - A_2 = 0$).
• La realtà: Il documento conferma che questa regola regge bene a basse velocità, ma man mano che le collisioni diventano più energetiche (maggiore momento), la regola inizia a crollare. Questo non è un fallimento; è una caratteristica! Ci dice che le parti "disordinate" della collisione (come particelle extra che vengono scagliate fuori) iniziano a contare.

I risultati: Dati contro teoria

Gli scienziati hanno confrontato le loro misurazioni con le migliori simulazioni al computer disponibili (le "previsioni teoriche").

• Le buone notizie: Per la maggior parte dei coefficienti, i dati del mondo reale corrispondevano molto bene ai modelli informatici. Questo significa che la nostra attuale comprensione di come queste particelle interagiscono è solida.
• La tensione interessante: Nella gamma intermedia di velocità, i dati per un coefficiente specifico ($A_0$) erano leggermente più alti di quanto previsto dal computer (circa 3 deviazioni standard fuori). È come se una previsione meteorologica prevedesse il 50% di probabilità di pioggia, ma in realtà avesse piovuto l'80% delle volte. Non è un disastro, ma suggerisce che il modello informatico potrebbe perdere un dettaglio minuscolo.
• I coefficienti "fantasma": Tre dei coefficienti ($A_5, A_6, A_7$) dovrebbero essere zero o molto vicini ad esso. I dati hanno mostrato che erano effettivamente minuscoli, coerenti con lo zero, anche se uno di essi ($A_6$) ha mostrato un piccolo, debole accenno di essere non nullo. È come sentire un sussurro in una stanza silenziosa; è lì, ma servono orecchie molto sensibili per sentirlo.

Perché questo è importante

Questo documento è essenzialmente un controllo di calibrazione ad alta precisione per le leggi della fisica.

1. Capire la "colla": Queste misurazioni ci aiutano a comprendere la "dinamica partonica"—come i mattoni fondamentali minuscoli all'interno del protone (quark e gluoni) si comportano quando collidono.
2. Testare la teoria: Confrontando lo "spin" del bosone Z con matematica complessa (Cromodinamica Quantistica), gli scienziati stanno sottoponendo a stress test la nostra comprensione dell'universo. Se la matematica non corrisponde ai dati, significa che dobbiamo inventare nuova fisica.
3. Il punto di riferimento: Questo documento fornisce un nuovo "righello" ultra-preciso per i futuri esperimenti. Qualsiasi nuova teoria deve essere in grado di spiegare questi otto numeri.

Sintesi

In breve, il team CMS ha scattato un'istantanea massiccia di collisioni di particelle, misurato gli angoli esatti delle particelle risultanti e calcolato otto numeri che descrivono lo "spin" dell'evento. Hanno scoperto che, sebbene le nostre teorie attuali siano per lo più corrette, ci sono minuscole e affascinanti discrepanze nella gamma di velocità intermedie che tengono i fisici in allerta, assicurando che la ricerca di una comprensione più profonda dell'universo continui.

Sommario tecnico

1. Problema e Motivazione

Il processo Drell–Yan (DY) ($pp \to \gamma^*/Z \to \ell^+\ell^-$) è una sonda fondamentale per comprendere la dinamica partonica sottostante delle collisioni protone-protone e i meccanismi di produzione dei bosoni vettoriali elettrodeboli. La distribuzione angolare dei leptoni nello stato finale nel sistema di riferimento a riposo del bosone codifica informazioni sugli stati di polarizzazione del fotone virtuale o del bosone Z.

Questa distribuzione è parametrizzata da otto coefficienti di polarizzazione angolare, $A_0$ fino ad $A_7$, definiti nel sistema di riferimento Collins–Soper (CS).

• Obiettivi Fisici:
  • $A_0$ e $A_2$: Parametrizzano la polarizzazione longitudinale e l'interferenza tra stati trasversi. La loro differenza ($A_0 - A_2$) mette alla prova la relazione di Lam–Tung ($A_0 - A_2 = 0$), che vale al Leading Order (LO) in QCD ma si prevede sia violata dalla radiazione QCD di ordine superiore e dalla quantità di moto trasversa intrinseca dei partoni.
  • $A_4$: Correlata all'asimmetria forward-backward ($A_{FB}$) e utilizzata per misurare l'angolo di mixing debole efficace ($\sin^2 \theta_{eff}^W$).
  • $A_1, A_3$: Codificano l'interferenza tra stati longitudinali e trasversi e dipendono dagli accoppiamenti assiali/vettoriali.
  • $A_5, A_6, A_7$: Definiscono asimmetrie T-dispari. Queste sono previste essere nulle a LO e NLO ma possono essere non nulle al Next-to-Next-to-Leading Order (NNLO) a causa di processi a un loop.

Sebbene esistessero misurazioni precedenti a $\sqrt{s}=8$ TeV e nella regione forward (LHCb), era necessaria una misurazione completa, doppio-differenziale dell'intero insieme di coefficienti ($A_0$–$A_7$) a $\sqrt{s}=13$ TeV con alta statistica per vincolare i modelli teorici e testare le previsioni QCD NNLO con maggiore precisione.

2. Metodologia

Dati e Selezione degli Eventi

• Dataset: Dati di collisione protone-protone raccolti dal rivelatore CMS a $\sqrt{s} = 13$ TeV durante il 2016–2018, corrispondenti a una luminosità integrata di 140 fb$^{-1}$.
• Segnale: Eventi dileptonici ($Z \to \mu^+\mu^-$) con massa invariante $81 < m_{\mu\mu} < 101$ GeV (regione del polo Z).
• Tagli Cinematici:
  • Quantità di moto trasversa del dimuone ($p_T^{\mu\mu}$): da 0 a 400 GeV.
  • Rapidità del dimuone ($|y_{\mu\mu}|$): $< 2.4$.
  • Criteri di isolamento e identificazione dei muoni applicati per garantire alta purezza.
  • Requisiti di trigger: trigger a singolo muone con $p_T > 24$ o $27$ GeV a seconda del periodo di presa dati.

Simulazione e Correzioni

• Simulazione del Segnale: Generata utilizzando POWHEG v2.0 + MiNNLOPS interfacciato con PYTHIA 8 (per showering di partoni/adronizzazione) e PHOTOS++ (per radiazione QED dello stato finale). È stato utilizzato il set PDF NNPDF3.1.
• Stima del Fondo: Simulazioni Monte Carlo per processi di dibosone (WW, WZ, ZZ), $t\bar{t}$ e $W$+jets. I contributi di fondo sono piccoli (0,03% a basso $p_T$ che sale a ~2,9% ad alto $p_T$).
• Correzioni:
  • Fattori di Scala: Derivati tramite "tag-and-probe" per correggere l'efficienza di trigger, identificazione e isolamento.
  • Ripesatura del Pileup: Per far corrispondere la distribuzione simulata del pileup ai dati.
  • Ripesatura Cinematica: Gli eventi simulati sono stati ripesati nelle bin $(p_T^{\mu\mu}, y_{\mu\mu})$ per far corrispondere le distribuzioni cinematiche dei dati, minimizzando la dipendenza dal modello.
  • Correzione Prefiring: Ha affrontato il disallineamento degli incroci dei fasci nel trigger di Livello 1.

Estrazione dei Coefficienti

L'estrazione utilizza un metodo a template per tenere conto dell'accettanza del rivelatore, della risoluzione e della migrazione tra le bin di $p_T$.

1. Definizione: La sezione d'urto differenziale è espansa come somma di polinomi trigonometrici $f_i(\theta^*, \phi^*)$ pesati da coefficienti $P_i$ (dove $A_i = P_i/P_8$).
2. Template: Sono stati costruiti nove template 2D ($T_i$) da eventi MC ricostruiti. Ogni template corrisponde a un termine angolare specifico $f_i$, pesato per rimuovere l'informazione di polarizzazione dalla generazione dell'evento preservando gli effetti del rivelatore.
3. Adattamento (Fitting): È stato eseguito un fit di massima verosimiglianza binnato bidimensionale (usando MINUIT) sulle distribuzioni angolari osservate ($\cos \theta^*, \phi^*$) in ciascuna delle 16 bin cinematiche (8 bin di $p_T$ × 2 bin di rapidità). Il fit estrae i coefficienti $P_i$ minimizzando la differenza tra i dati e la combinazione lineare dei template.

3. Contributi Chiave

• Primo Insieme Completo a 13 TeV: Questa è la prima misurazione dell'insieme completo di coefficienti angolari ($A_0$–$A_7$) nella regione di rapidità centrale ($|y_{\mu\mu}| < 2.4$) a $\sqrt{s}=13$ TeV.
• Precisione Doppio-Differenziale: I risultati sono forniti in 8 bin di $p_T^{\mu\mu}$ e 2 bin di $|y_{\mu\mu}|$, offrendo una visione multidimensionale della struttura angolare.
• Confronto NNLO: I risultati sono confrontati con previsioni teoriche all'avanguardia a NNLO in QCD (usando POWHEG+MiNNLOPS e MADGRAPH5_aMC@NLO).
• Vincoli sulle Asimmetrie T-Dispari: Fornisce i vincoli più stringenti finora sui coefficienti T-dispari ($A_5, A_6, A_7$) nella regione centrale.

4. Risultati

Accordo Generale

• I coefficienti misurati concordano generalmente con le previsioni QCD NNLO entro le incertezze.
• Il grande dataset (140 fb$^{-1}$) ha ridotto significativamente le incertezze statistiche rispetto alle precedenti misurazioni del Run 1.

Osservazioni Specifiche

• $A_0$ e $A_2$:
  • $A_0$ mostra una discrepanza di circa 3 deviazioni standard nell'intervallo intermedio di $p_T$ (10–35 GeV) tra i dati e la previsione MADGRAPH5_aMC@NLO, mentre POWHEG+MiNNLOPS concorda bene.
  • $A_2$ è leggermente sovrastimata dalle simulazioni MC ad alto $p_T$ ($>50$ GeV).
• Relazione di Lam–Tung ($A_0 - A_2$):
  • La relazione è violata, come previsto a causa degli effetti QCD di ordine superiore.
  • Le previsioni POWHEG+MiNNLOPS concordano bene con i valori misurati di $A_0 - A_2$ nell'intervallo 10–35 GeV.
  • Entrambe le previsioni MC sottostimano la misurazione per $p_T > 35$ GeV.
• $A_1$:
  • Significativamente sovrastimata da POWHEG+MiNNLOPS nella regione $1 < |y_{\mu\mu}| < 2.4$ per $p_T < 35$ GeV.
  • MADGRAPH5_aMC@NLO sovrastima anch'essa $A_1$ a basso $p_T$ in entrambe le regioni di rapidità.
• Coefficienti T-Dispari ($A_5, A_6, A_7$):
  • I valori sono dell'ordine di $10^{-3}$ e consistenti con zero nella maggior parte delle bin.
  • La deviazione più significativa è per $A_6$ nella bin $55 < p_T < 80$ GeV e $1.2 < |y_{\mu\mu}| < 2.4$, che mostra una deviazione di 2,8$\sigma$ da zero. Questo è coerente con le aspettative NNLO ma di magnitudine inferiore rispetto alle precedenti osservazioni ATLAS a 8 TeV.
• Incertezze Sistematiche:
  • Dominanti per $A_0$ e $A_1$ a basso $p_T$ ($<55$ GeV) e per $A_2, A_3, A_4$ fino a $p_T < 200$ GeV.
  • Le incertezze statistiche dominano per gli altri coefficienti e nelle regioni di $p_T$ più elevate.

5. Significato

• Validazione Teorica: I risultati forniscono un benchmark rigoroso per i calcoli QCD NNLO, confermando che le correzioni di ordine superiore sono essenziali per descrivere la struttura angolare della produzione Drell–Yan.
• Dinamica dei Partoni: Le misurazioni offrono nuovi vincoli sulle Funzioni di Distribuzione dei Partoni (PDF) e sulla modellazione della quantità di moto trasversa intrinseca dei partoni.
• Fisica Futura: La precisione e la suddivisione in bin multidimensionali stabilite qui servono come riferimento critico per futuri studi fenomenologici e per testare calcoli di ordine ancora superiore (ad es. N$^3$LO).
• Prestazioni del Rivelatore: Dimostra la capacità del rivelatore CMS di misurare sottili correlazioni angolari con alta precisione, anche in presenza di complessi effetti del rivelatore e pileup.

In sintesi, questo documento rappresenta una pietra miliare nella fisica di precisione elettrodebole, confermando la validità delle descrizioni QCD NNLO del meccanismo di produzione del bosone Z, evidenziando al contempo specifiche regioni cinematiche in cui i modelli teorici potrebbero richiedere un ulteriore affinamento.