⚛️ phenomenology

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

Ce papier présente LHEReader, un programme permettant de convertir facilement des fichiers Les Houches Event au format ROOT pour faciliter leur analyse, dont la performance et l'application sont illustrées via des simulations de processus de collision proton-proton et électron-positron.

Auteurs originaux : Aman Desai

Publié 2026-03-03

📖 4 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Aman Desai

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

📂 LHEReader : Le Traducteur Magique pour les Physiciens

Imaginez que vous êtes un physicien qui étudie l'univers en recréant des collisions de particules (comme au CERN). Pour cela, vous utilisez des super-ordinateurs qui génèrent des milliards de "scénarios" de collisions.

1. Le Problème : Des fichiers incompréhensibles

Ces ordinateurs produisent des fichiers appelés LHE (Les Houches Event).

L'analogie : Imaginez que ces fichiers LHE sont comme des recettes de cuisine écrites dans une langue étrangère obscure, avec des codes secrets et des abréviations que seul le chef (le logiciel de simulation) comprend.
Le souci : Vous, le physicien, voulez analyser ces recettes pour voir si le plat final (les particules) est bon. Mais vous ne pouvez pas lire la recette directement. Vous avez besoin d'un outil pour la traduire dans un langage que vous maîtrisez, comme le ROOT (le langage standard des physiciens des particules, un peu comme Excel ou un tableur géant).

2. La Solution : LHEReader

C'est là qu'intervient LHEReader, le programme présenté par Aman Desai.

L'analogie : LHEReader est un traducteur automatique ultra-rapide. Il prend votre fichier LHE (la recette obscure) et le convertit instantanément en un fichier ROOT (un tableau clair et lisible).
Ce qu'il fait : Il ne se contente pas de traduire ; il organise tout. Il prend les informations sur chaque particule (sa masse, son énergie, sa vitesse, sa "famille") et les range dans des colonnes bien ordonnées, prêtes à être analysées.

3. La Vitesse : Un express de haute technologie

L'auteur a testé la vitesse de ce traducteur.

Le résultat : C'est étonnamment rapide. Pour convertir un seul événement (une seule collision simulée), il faut environ 0,2 milliseconde.
L'analogie : C'est comme si vous pouviez traduire un livre entier en une seconde. Même avec des millions de collisions à traiter, le programme ne transpire pas. Il est aussi capable de lire des fichiers compressés (comme des valises fermées à clé .lhe.gz) et de les ouvrir sans effort.

4. L'Exemple Pratique : Chasser le Boson de Higgs

Pour montrer que son outil fonctionne, l'auteur a simulé une collision très célèbre : Proton + Proton → Boson Z + Boson de Higgs.

Le scénario :
1. On fait entrer en collision deux protons (comme deux voitures de course qui s'écrasent).
2. Cela crée un Boson Z et un Boson de Higgs.
3. Le Z se désintègre en deux leptons (électrons ou muons).
4. Le Higgs se désintègre en deux quarks "b" (qui deviennent des jets de particules).
L'analyse : Grâce à LHEReader, le physicien peut maintenant "regarder" les données. Il utilise un petit script (appelé EventLoop) pour dire : "Montre-moi seulement les collisions où il y a deux électrons et deux jets de particules, et calcule leur énergie."
Le résultat : Le programme génère des graphiques (des histogrammes) qui ressemblent à des courbes de température ou de vitesse, montrant si les particules trouvées correspondent bien à ce que la théorie prédit.

5. Pourquoi c'est utile ?

Avant, convertir ces fichiers était une tâche lourde et complexe. Avec LHEReader :

Pour les chercheurs : C'est un outil pour tester rapidement de nouvelles idées sans perdre des heures à coder des convertisseurs.
Pour les étudiants : C'est un outil pédagogique léger pour apprendre comment les données sont structurées sans se perdre dans la complexité du code source.

En résumé

LHEReader est comme un pont entre le monde brut et cryptique des simulations informatiques et le monde clair et visuel de l'analyse scientifique. Il permet aux physiciens de passer du "code source" aux "graphiques intelligibles" en un clin d'œil, facilitant ainsi la découverte de nouveaux secrets de l'univers.

1. Problématique

Dans le domaine de la physique des hautes énergies, les générateurs d'événements Monte Carlo (comme MadGraph5_aMC@NLO) produisent des enregistrements d'événements au niveau des partons au format standard Les Houches Event (LHE). Ces fichiers contiennent des informations cruciales (identités des particules, codes de statut, quadri-impulsions, spins, informations de couleur, poids d'événement, etc.).

Cependant, l'analyse de ces données est souvent réalisée à l'aide du framework ROOT (développé par le CERN), qui est l'outil standard pour le traitement et la visualisation de données en physique. Le défi réside dans la nécessité de convertir les fichiers LHE (souvent volumineux et textuels) en fichiers ROOT binaires (.root) pour pouvoir les exploiter efficacement avec les fonctionnalités de ROOT. Bien que des solutions existent (comme pylhe), l'auteur propose une alternative plus légère et intégrée spécifiquement pour faciliter cette conversion et l'analyse subséquente.

2. Méthodologie

L'auteur présente LHEReader, un programme écrit en Python (compatible avec PyROOT) conçu pour convertir directement les fichiers .lhe (ou .lhe.gz) en fichiers .root.

Architecture et Installation : Le programme est open-source (hébergé sur GitHub). Il nécessite un environnement Linux avec ROOT et PyROOT installés, et a été testé avec Python 3.10. L'installation se fait via git clone et l'exécution se fait par ligne de commande.
Processus de Conversion : LHEReader lit le fichier LHE, extrait les informations pour chaque particule (PID, statut, mère, couleur, quadri-impulsion, etc.) et les stocke dans un objet TTree ROOT nommé events.
Structure des Données : Le fichier ROOT généré contient un TTree avec des branches spécifiques :
- Niveau Événement : Nombre de particules, poids de l'événement, échelle, constantes de couplage ( $\alpha_{QED}$ , $\alpha_{QCD}$ ).
- Niveau Particule : Vecteurs pour les identifiants (PID), codes de statut, indices des particules mères, codes de couleur, et les composantes de l'impulsion ( $p_x, p_y, p_z$ ), l'énergie, la masse, la durée de vie et le spin.
Outil d'Analyse (EventLoop) : Le papier fournit également un exemple de code C++ (EventLoop) démontrant comment parcourir le TTree généré. Ce code illustre comment :
- Filtrer les particules par code de statut (ex: particules finales avec statut = 1).
- Identifier des particules spécifiques via leur code PDG (ex: quarks bottom, leptons).
- Calculer l'énergie manquante ( $E_{miss}^T$ ou MET) en sommant les vecteurs des neutrinos.
- Appliquer des sélections cinématiques (ex: $p_T > 25$ GeV).

3. Contributions Clés

Outil de Conversion Rapide : Développement d'un script Python léger et facile à déployer pour transformer des fichiers LHE en ROOT, éliminant la barrière d'entrée pour les utilisateurs souhaitant utiliser ROOT sans écrire de code de parsing complexe.
Support de Formats Compressés : Mise à jour récente permettant la lecture directe de fichiers .lhe.gz.
Guide Pratique et Éducatif : Fourniture d'un exemple complet de boucle d'événement (EventLoop) en C++ montrant comment accéder aux données, filtrer les événements et reconstruire des observables physiques (comme le MET ou les masses invariantes).
Validation sur Cas Réels : Le programme a été testé et utilisé dans des projets de recherche réels liés au LHC et au FCCee.

4. Résultats

L'évaluation de la performance a été réalisée sur deux processus de simulation Monte Carlo générés par MadGraph5_aMC@NLO :

$pp \to jj$ (collision proton-proton vers jets).
$e^+e^- \to ZH$ (collision électron-positron vers boson Z et Higgs).

Performance de conversion :

Pour $pp \to jj$ , le temps moyen de conversion par événement est d'environ 0,20 ms pour un échantillon de $10^5$ événements.
Pour $e^+e^- \to ZH$ , le temps moyen est encore plus rapide, d'environ 0,04 ms par événement pour le même nombre d'événements.
La performance s'améliore légèrement avec le nombre d'événements traités (passant de 0,27 ms à 0,18 ms pour $pp \to jj$ lorsque le nombre d'événements passe de $10^4$ à $10^6$ ).

Application Physique ( $pp \to ZH$ ) :
L'auteur a appliqué LHEReader à une simulation de $pp \to Z(\to \ell^+\ell^-)H(\to b\bar{b})$ à $\sqrt{s} = 13,6$ TeV (Run 3 du LHC).

Le flux de travail a inclus la génération LHE, le showering et l'hadronisation via Pythia8, et le clustering des jets via MadAnalysis5 (algorithme anti- $k_T$ ).
Une analyse au niveau des partons et au niveau reconstruit a été effectuée.
Les distributions cinématiques (impulsion transverse des leptons et des jets $b$ , masse invariante des paires de leptons et des paires de jets $b$ ) ont été générées avec succès, montrant des pics cohérents avec la masse du boson Z et du boson de Higgs.

5. Signification

Ce travail est significatif pour plusieurs raisons :

Accessibilité : Il fournit un outil « prêt à l'emploi » pour les chercheurs et les étudiants qui souhaitent analyser des données LHE sans maîtriser les bibliothèques complexes de parsing XML ou les interfaces C++ lourdes.
Efficacité : La vitesse de conversion (de l'ordre de la fraction de milliseconde par événement) rend l'outil viable pour traiter de grands échantillons de données Monte Carlo.
Polyvalence : En permettant l'analyse au niveau des partons (via les étiquettes de vérité Monte Carlo) et au niveau reconstruit (via les jets et les PID), il sert de pont entre la génération d'événements et l'analyse physique finale.
Utilité Éducative : Le code fourni sert de modèle pédagogique pour apprendre à manipuler des fichiers ROOT et à écrire des boucles d'événements en physique des hautes énergies.

En résumé, LHEReader simplifie le flux de travail de l'analyse de données en physique des particules en automatisant la conversion de formats de données brutes vers un format d'analyse standard, tout en fournissant les outils nécessaires pour exploiter immédiatement ces données.