LHEReader: Simplified Conversion from Les Houches Event Files to ROOT Format

Il documento presenta LHEReader, un programma che semplifica la conversione dei file Les Houches Event in formato ROOT per facilitarne l'analisi, ne valuta le prestazioni attraverso simulazioni di processi come $pp\to jj$ ed $e^+e^- \to ZH$, e ne illustra l'applicazione pratica nell'analisi di eventi a livello di partoni e dopo simulazione rapida.

L'essenziale

🎬 Il Regista che traduce il copione: LHEReader

Immagina di essere un regista cinematografico. Hai appena finito di girare un film (una simulazione di collisioni di particelle) e hai il copione originale scritto in una lingua molto tecnica e complessa, chiamata LHE (Les Houches Event). Questo copione contiene tutto: chi sono gli attori (le particelle), cosa fanno, quanto pesano e come si muovono.

Il problema? Il tuo team di montaggio e gli analisti (i fisici) lavorano con un software diverso, chiamato ROOT, che è come il loro "studio di montaggio" preferito. Ma ROOT non capisce il copione originale LHE. È come se avessi un libro scritto in cinese e volessi analizzarlo usando un software che legge solo l'italiano.

Qui entra in gioco LHEReader, il protagonista di questo paper.

🛠️ Cos'è LHEReader? (Il Traduttore Magico)

LHEReader è un piccolo programma (un "ponte" o un "traduttore") scritto dal dottor Aman Desai. La sua missione è semplice ma potente: prende quel copione complicato (il file .lhe) e lo traduce istantaneamente in un formato che ROOT può leggere (un file .root).

Grazie a questo traduttore, i fisici possono prendere i dati grezzi delle simulazioni e iniziare subito a fare analisi, a creare grafici e a cercare nuove scoperte, senza impazzire a scrivere codice complesso per leggere i file originali.

⏱️ Quanto è veloce? (La gara di corsa)

Gli autori hanno messo alla prova questo traduttore facendogli convertire milioni di "eventi" (immagina di dover tradurre milioni di pagine di copione in un battito di ciglia).

• Hanno simulato due scenari: uno come una collisione di due proiettili (protoni) e uno come un incontro tra elettroni e positroni.
• Il risultato? È velocissimo! Traduce un singolo evento in circa 0,04 millisecondi. È come se potesse leggere e tradurre un'intera enciclopedia in meno tempo di quanto ci voglia a fare un battito di ciglia.

🎯 L'Esempio Pratico: Caccia al "Mostro" (Il Bosone di Higgs)

Per dimostrare quanto sia utile, gli autori hanno usato LHEReader per simulare un esperimento reale al LHC (il grande acceleratore di particelle).
Hanno cercato di ricreare una collisione che produce il Bosone di Higgs (il "Santo Graal" della fisica delle particelle, che dà massa alle cose).

1. La Simulazione: Immagina di far scontrare due treni ad alta velocità. Da questo scontro escono scintille, frammenti e nuove particelle.
2. La Traduzione: LHEReader prende i dati di questo scontro e li trasforma in un file che ROOT può aprire.
3. L'Analisi: Usando uno strumento chiamato EventLoop (una sorta di "robot ispettore" che scorre i dati), i fisici hanno filtrato gli eventi:
   • Hanno cercato coppie di elettroni o muoni (come due gemelli che scappano dalla scena del crimine).
   • Hanno cercato coppie di "quark bottom" (i mattoni fondamentali che formano il Bosone di Higgs).
4. Il Risultato: Hanno ricostruito la massa di queste particelle e hanno visto un picco preciso proprio dove la teoria dice che dovrebbe esserci il Bosone di Higgs. È come se, guardando i detriti di un incidente, fossero riusciti a ricostruire esattamente il modello dell'auto che si era schiantata.

🌟 Perché è importante?

Prima di questo programma, convertire questi file era un lavoro manuale, lento e soggetto a errori.
LHEReader è come un traduttore universale istantaneo:

• Per i ricercatori: Risparmia ore di lavoro, permettendo loro di concentrarsi sulla fisica e non sulla formattazione dei file.
• Per gli studenti: È uno strumento leggero e facile da usare per imparare come funzionano le collisioni di particelle senza dover essere esperti di programmazione.

In sintesi, LHEReader è il piccolo eroe silenzioso che prende il caos dei dati grezzi dell'universo e li organizza in una biblioteca ordinata, pronta per essere esplorata da chiunque voglia scoprire i segreti della materia.

Sommario tecnico

Titolo: LHEReader: Conversione Semplificata da File Les Houches Event a Formato ROOT

1. Il Problema

Nella fisica delle alte energie, i generatori di eventi Monte Carlo (come MadGraph5_aMC@NLO) producono record di eventi a livello di partoni utilizzando un formato standardizzato chiamato Les Houches Event File (LHE). Questi file contengono informazioni critiche come identità delle particelle, codici di stato, quadri-impulsi, spin, informazioni sui genitori, dati di colore, scale di energia e pesi degli eventi.
Tuttavia, per l'analisi successiva, i fisici utilizzano comunemente il framework ROOT (sviluppato dal CERN). Sebbene esistano software per leggere i file LHE, la necessità di convertire efficientemente questi file nel formato binario di ROOT (.root) per sfruttare le funzionalità di analisi, visualizzazione e gestione dati di ROOT rimane un passaggio cruciale. Il problema risiede nella necessità di uno strumento leggero, accessibile e performante che faciliti questa conversione e l'analisi immediata dei dati, sia a livello di partoni che dopo la simulazione rapida (fast simulation).

2. Metodologia

L'autore presenta LHEReader, un programma scritto in Python (testato con Python 3.10) che automatizza la conversione dei file .lhe (e .lhe.gz) in file .root compatibili con CERN ROOT.

• Architettura e Dipendenze: Il programma richiede l'installazione di ROOT e PyROOT. È disponibile tramite repository GitHub e può essere installato via terminale (git clone).
• Processo di Conversione: Lo script legge la struttura XML del file LHE, estrae i dati per ogni particella e per ogni evento, e li scrive in un oggetto TTree di ROOT.
• Struttura dei Dati in ROOT: Il file di output contiene un TTree chiamato events con diverse "branche" (branch) che mappano i dati originali:
  • Dati dell'evento: numero di particelle, peso dell'evento, scala, $\alpha_{QED}$, $\alpha_{QCD}$.
  • Dati delle particelle (vettori): PID (codice PDG), stato, madri (mother1, mother2), colore, quadri-impulsi ($p_x, p_y, p_z, E$), massa, vita media ($\tau$) e spin.
• Analisi con EventLoop: Il documento fornisce un esempio di codice C++ che utilizza la funzione EventLoop per iterare sugli eventi nel file ROOT. Questo codice dimostra come:
  • Filtrare le particelle in base al codice di stato (es. stato 1 per particelle finali).
  • Identificare particelle specifiche tramite il codice PDG (es. quark bottom con ID 5, leptoni con ID 11/13).
  • Calcolare l'energia mancante ($E_T^{miss}$) sommando i vettori delle particelle invisibili (neutrini).
  • Applicare tagli di selezione (es. $p_T > 25$ GeV, $|\eta| < 2.5$).

3. Contributi Chiave

• Strumento di Conversione Accessibile: Fornisce uno strumento open-source (disponibile su GitHub) che semplifica drasticamente il flusso di lavoro per convertire dati grezzi LHE in un formato analizzabile con ROOT.
• Supporto Multi-formato: Oltre ai file .lhe standard, il programma gestisce file compressi .lhe.gz.
• Esempio Pratico Completo: Il paper non si limita alla conversione, ma fornisce un framework di analisi completo (EventLoop) che mostra come accedere ai dati, filtrare eventi e calcolare osservabili cinematiche.
• Validazione su Casi d'Uso Reali: Il programma è stato testato e utilizzato in progetti di ricerca reali relativi all'LHC e al futuro collisore FCCee.

4. Risultati

L'efficacia del programma è stata valutata simulando due processi: $pp \to jj$ ed $e^+e^- \to ZH$.

• Prestazioni di Conversione:
  • Per il processo $pp \to jj$ (100.000 eventi), il tempo medio di conversione è di 0.20 ms per evento.
  • Per il processo $e^+e^- \to ZH$ (100.000 eventi), il tempo medio scende a 0.04 ms per evento.
  • Il tempo di elaborazione scala linearmente con il numero di eventi, dimostrando alta efficienza anche per grandi campioni (fino a $10^6$ eventi).
• Applicazione Fisica ($pp \to ZH$):
  • È stata simulata la produzione di un bosone Z (che decade in $e^+e^-$ o $\mu^+\mu^-$) e un bosone di Higgs (che decade in $b\bar{b}$) a $\sqrt{s} = 13.6$ TeV.
  • Gli eventi sono stati elaborati attraverso Pythia8 (parton shower e adronizzazione) e MadAnalysis5 (fast simulation e clustering jet con algoritmo anti-$k_T$, $R=0.4$).
  • Risultati dell'analisi: Sono state ricostruite distribuzioni cinematiche chiave:
    • Momento trasverso ($p_T$) di muoni, elettroni e jet $b$.
    • Massa invariante dei leptoni (picco attorno alla massa dello Z).
    • Massa invariante della coppia di jet $b$ (distribuzione attorno alla massa dell'Higgs a livello ricostruito).
  • L'analisi ha confermato la capacità di selezionare eventi con due leptoni isolati e due jet $b$, applicando tagli cinematici standard.

5. Significato

Il lavoro di Desai offre un ponte fondamentale tra la generazione di eventi Monte Carlo e l'analisi dati avanzata.

• Utilità per la Ricerca: Permette ai ricercatori di testare rapidamente distribuzioni cinematiche e strategie di selezione senza dover scrivere parser complessi per i file LHE. È stato già utilizzato in studi per LHC Run 3 e FCCee.
• Utilità Educativa: Si presenta come uno strumento "leggero" e didattico, ideale per studenti e ricercatori che necessitano di convertire file LHE in ROOT per scopi didattici o per prototipazione rapida di analisi.
• Integrazione: Dimostra come integrare flussi di lavoro eterogenei (MadGraph, Pythia, MadAnalysis, ROOT) in un'unica pipeline coerente, facilitando la riproducibilità e l'efficienza nell'analisi dei dati di fisica delle alte energie.

In sintesi, LHEReader risolve un collo di bottiglia tecnico comune nella fisica delle particelle, fornendo una soluzione robusta, veloce e ben documentata per la transizione dai dati grezzi di simulazione all'analisi statistica.