⚛️ phenomenology

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

Dit paper introduceert LHEReader, een programma dat Les Houches Event-bestanden converteert naar ROOT-formaat voor analyse, met prestatietests en een toepassing op $pp \to ZH$ -simulaties.

Oorspronkelijke auteurs: Aman Desai

Gepubliceerd 2026-03-03

📖 4 min leestijd🧠 Diepgaand

CC BY 4.0

Oorspronkelijke auteurs: Aman Desai

Oorspronkelijk artikel gelicentieerd onder CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Dit is een AI-gegenereerde uitleg van het onderstaande artikel. Het is niet geschreven of goedgekeurd door de auteurs. Raadpleeg het oorspronkelijke artikel voor technische nauwkeurigheid. Lees de volledige disclaimer

LHEReader: De Vertaler voor de Deeltjesfysica

Stel je voor dat je een enorme bibliotheek hebt vol met de meest complexe wetenschappelijke rapporten ter wereld, maar ze zijn allemaal geschreven in een taal die alleen maar door een zeer klein groepje specialisten wordt begrepen. Dat is precies wat er gebeurt in de wereld van de deeltjesfysica.

1. Het Probleem: De "Geheime Code"

Wanneer wetenschappers botsingen simuleren in deeltjesversnellers (zoals de LHC in Genève), gebruiken ze software om te berekenen wat er gebeurt. Deze software schrijft de resultaten op in een bestandstype dat LHE (Les Houches Event) heet.

De Analogie: Stel je voor dat LHE-bestanden zijn als een recept in een vreemde taal, geschreven met een eigen soort symbolen en afkortingen. Het bevat alle informatie: welke deeltjes er zijn, hoe snel ze gaan, wat hun "stamboom" is en hoeveel energie ze hebben. Maar voor de meeste onderzoekers die de data daadwerkelijk willen zien (in grafieken en diagrammen), is dit bestand te onhandig. Het is alsof je een recept hebt, maar je wilt het niet lezen, je wilt het gerecht proeven.

2. De Oplossing: LHEReader

Aman Desai, een onderzoeker van de Universiteit van Adelaide, heeft een programma gemaakt genaamd LHEReader.

De Analogie: LHEReader is als een super-snelle vertaler of een kookmachine. Je gooit het ingewikkelde LHE-recept erin, en het programma spuw het direct uit in een ROOT-bestand.
Wat is ROOT? ROOT is een soort "super-Excel" of een krachtige database die door bijna alle deeltjesfysici wordt gebruikt. Het is de taal waarin je de data kunt visualiseren, kunt meten en kunt analyseren. LHEReader maakt de brug tussen de ruwe, onleesbare data en de bruikbare, visuele data.

3. Hoe snel is het? (De Prestaties)

De auteur heeft getest hoe snel deze vertaler werkt. Hij liet het programma duizenden en miljoenen "virtuele botsingen" verwerken.

Het Resultaat: Het is verrassend snel. Het kost gemiddeld slechts 0,04 tot 0,20 milliseconden om één enkel botsing-gebeurtenis om te zetten.
De Vergelijking: Dat is sneller dan het knipperen van een oog. Je kunt dus een hele berg data (zoals een miljoen botsingen) in een paar minuten omzetten, terwijl het anders uren zou duren.

4. Een Praktijkvoorbeeld: De Higgs-deeltjes Jacht

Om te laten zien hoe goed het werkt, heeft de auteur een specifieke simulatie gedaan: het produceren van een Higgs-boson (het "Goddeeltje") en een Z-boson.

Het Scenario: Ze simuleerden botsingen waarbij deze deeltjes ontstaan en vervolgens vervallen in andere deeltjes (zoals elektronen, muonen en b-quarks).
De Analyse: Met LHEReader en een hulpprogramma genaamd EventLoop (een soort "loop" die door de data loopt), konden ze filteren op de deeltjes die ze zochten. Ze keken bijvoorbeeld: "Hebben we twee elektronen die snel genoeg zijn?" en "Zijn er twee b-quarks die samen de massa van een Higgs-deeltje hebben?"
Het Eindresultaat: Ze kregen prachtige grafieken die lieten zien hoe de energie en massa van deze deeltjes zich verhouden. Dit helpt wetenschappers om te zien of hun theorieën kloppen met wat er in de echte wereld gebeurt.

5. Waarom is dit belangrijk?

Vroeger was het moeilijk om deze bestanden te lezen; je had vaak complexe, specifieke software nodig die niet altijd werkte.

De Impact: LHEReader is als het openen van een deur. Het maakt het voor studenten, onderzoekers en geïnteresseerden heel makkelijk om de ruwe data van simulaties te nemen en direct te analyseren in een omgeving die ze al kennen (ROOT).
Het is een lichtgewicht tool: het is niet zwaar of ingewikkeld, maar het doet precies wat het moet doen: het zet de "ruwe diamant" (LHE) om in een "geslepen edelsteen" (ROOT) die je kunt bewonderen.

Samenvattend

Dit paper introduceert LHEReader, een handig stukje software dat de taal van deeltjesfysica-simulaties vertaalt naar een formaat dat iedereen kan gebruiken om te analyseren. Het is de sleutel die de deur opent naar het begrijpen van de meest fundamentele bouwstenen van ons universum, zonder dat je eerst jaren hoeft te studeren om de code te ontcijferen.

Probleemstelling

In de hoge-energie fysica worden Monte Carlo (MC) gebeurtenisgeneratoren (zoals MadGraph5_aMC@NLO) gebruikt om deeltjesinteracties op partonniveau te simuleren. Deze generatoren slaan de gegenereerde gebeurtenisrecords op in het standaardformaat Les Houches Event (LHE). Hoewel LHE-bestanden alle noodzakelijke informatie bevatten (deeltjesidentiteiten, statuscodes, viermomenta, spin, ouder-deeltjesinformatie, kleurencodes, schaalfactoren en gewichten), zijn ze niet direct compatibel met de geavanceerde analyse-omgevingen die veelal worden gebruikt in de deeltjesfysica, zoals ROOT (ontwikkeld door CERN).

Fysici moeten LHE-bestanden vaak converteren naar ROOT-bestanden om gebruik te maken van de krachtige analyse-functionaliteiten van ROOT, zoals het maken van histograms, het definiëren van complexe selecties en het uitvoeren van statistische analyses. Bestaande tools voor deze conversie kunnen complex zijn of beperkte functionaliteit bieden. Er is behoefte aan een lichtgewicht, gebruiksvriendelijke en efficiënte oplossing om LHE-bestanden (inclusief gecomprimeerde .lhe.gz bestanden) direct om te zetten naar een ROOT-formaat dat direct leesbaar is voor analyse-scripts.

Methodologie

De auteur, Aman Desai, presenteert LHEReader, een Python-programma dat is ontworpen om LHE-bestanden te parseren en om te zetten naar een ROOT-bestand.

Software Architectuur:
- Het programma is geschreven in Python en maakt gebruik van de PyROOT bibliotheek (de Python-bindingen voor ROOT).
- Het is compatibel met Linux-systemen en vereist Python 3.10 en een geïnstalleerde CERN ROOT-installatie.
- Het ondersteunt zowel .lhe als .lhe.gz (gecomprimeerd) als invoer.
Conversieproces:
- Het script leest het XML-achtige LHE-bestand regel voor regel.
- Het extrahert de relevante gegevens voor elke gebeurtenis en elk deeltje.
- De gegevens worden gestructureerd in een TTree (de standaard datastructuur in ROOT voor grote datasets) genaamd events.
- De TTree bevat diverse "branches" (takken) die corresponderen met de eigenschappen van de gebeurtenis en de deeltjes.
Datastructuur in ROOT:
De gegenereerde ROOT-bestand bevat een TTree met de volgende branches (zie Tabel IV in het paper):
- Gebeurtenis-variabelen: Aantal deeltjes (numParticles), gebeurtenisgewicht (eventweight), schaal (scale), $\alpha_{QED}$ , en $\alpha_{QCD}$ .
- Deeltjes-variabelen (vector per deeltje): PDG-code (pid), statuscode, moeder-deeltjes (mother1, mother2), kleurcodes, viermomentum (px, py, pz, energy), massa, levensduur (tau) en spin.
Analyse Framework (EventLoop):
- Het paper biedt een voorbeeld van een EventLoop in C++ die de gegenereerde ROOT-bestanden kan verwerken.
- Deze loop itereert over de gebeurtenissen en de deeltjes binnen die gebeurtenissen.
- Het demonstreert hoe men statuscodes kan gebruiken om alleen eindtoestandsdeeltjes te selecteren, hoe men PDG-codes gebruikt om specifieke deeltjes (zoals b-quarks of leptonen) te filteren, en hoe men Missing Transverse Energy (MET) berekent op basis van neutrino's.

Belangrijkste Bijdragen

LHEReader Tool: Een open-source tool (beschikbaar via GitHub) die de conversie van LHE naar ROOT vereenvoudigt. De tool is lichtgewicht en vereist geen complexe afhankelijkheden buiten de standaard HEP-softwarestack.
Ondersteuning voor Gecomprimeerde Bestanden: De tool kan direct .lhe.gz bestanden verwerken, wat essentieel is voor het hanteren van grote datasets.
Voorbeeld Implementatie: Het paper levert niet alleen de conversietool, maar ook een volledig werkend C++-voorbeeld (EventLoop) dat fysici direct kunnen gebruiken om kinematische distributies te analyseren, inclusief het berekenen van MET en het selecteren van specifieke deeltjestypen.
Validatie en Prestaties: De tool is getest op twee verschillende simulatiescenaria en de prestaties zijn kwantitatief geëvalueerd.

Resultaten

De auteurs hebben de prestaties van de conversie geëvalueerd door twee sets Monte Carlo-simulaties te draaien:

$pp \to jj$ (QCD-jet productie).
$e^+e^- \to ZH$ (Higgs-productie in een leptonen-collider omgeving).

Conversiesnelheid:

Voor $pp \to jj$ met $10^5$ gebeurtenissen: gemiddeld 0,20 ms per gebeurtenis.
Voor $e^+e^- \to ZH$ met $10^5$ gebeurtenissen: gemiddeld 0,04 ms per gebeurtenis.
De conversietijd neemt licht af bij grotere aantallen gebeurtenissen (vanwege efficiëntere I/O), wat aangeeft dat de tool schaalbaar is voor grote datasets.

Fysieke Analyse (Toepassing $pp \to ZH$ ):
De tool werd toegepast op een simulatie van $pp \to ZH$ (waarbij $Z \to \ell^+\ell^-$ en $H \to b\bar{b}$ ) bij een bundelenergie van 6,8 TeV (LHC Run 3).

De gegenereerde LHE-bestanden werden verwerkt via Pythia8 (voor parton-shower en hadronisatie) en MadAnalysis5 (voor jet-clustering met de anti- $k_T$ algoritme).
Na conversie naar ROOT en toepassing van de EventLoop werden kinematische distributies gegenereerd:
- Transversale impuls ( $p_T$ ) van muonen, elektronen en b-jets.
- Invariante massa van de leptonparen (herkenning van de Z-boson piek rond 91 GeV).
- Invariante massa van de b-jet paren (herkenning van de Higgs-boson piek rond 125 GeV).
De resultaten tonen aan dat de data-integriteit behouden blijft tijdens de conversie en dat de analyse op zowel partonniveau als reconstrueerd niveau correct verloopt.

Betekenis en Impact

De LHEReader-tool is significant voor de deeltjesfysica-gemeenschap om de volgende redenen:

Toegankelijkheid: Het verlaagt de drempel voor onderzoekers (en studenten) om LHE-bestanden te analyseren zonder complexe scripts te hoeven schrijven. Het fungeert als een "brug" tussen generatie en analyse.
Efficiëntie: Met conversietijden van minder dan 0,2 ms per gebeurtenis is de tool snel genoeg voor het verwerken van grote datasets die nodig zijn voor LHC-analyses.
Educatieve Waarde: Het biedt een lichtgewicht hulpmiddel voor onderwijsdoeleinden, waarmee studenten kunnen leren hoe ze ruwe simulatiegegevens kunnen omzetten naar analyseerbare formaten.
Flexibiliteit: De tool is reeds gebruikt in onderzoeksprojecten gerelateerd aan de LHC en de toekomstige FCCee (Future Circular Collider), wat aantoont dat het een robuuste oplossing is voor zowel huidige als toekomstige experimenten.

Kortom, LHEReader lost een praktisch probleem op in de data-pipeline van de deeltjesfysica door een snelle, betrouwbare en gebruiksvriendelijke conversie van LHE naar ROOT te bieden, ondersteund door een compleet analyse-voorbeeld.