Live pilot
University of Edinburgh
boosted di-boson resonance search
First real-world deployment. Real Geant4 data. Real benchmark.
The analysis
Boosted di-boson resonance search with mass-decorrelated jet tagging
The search targets a heavy resonance X decaying to a pair of electroweak bosons (H → ZV, W′ → WZ, Z′ → ZH) in the boosted regime, where the decay products of each boson merge into a single large-radius fat jet. The analysis reconstructs two fat jets and searches for a bump in the reconstructed di-jet invariant mass m_jj above a smoothly-falling QCD background.
The mass sculpting problem is acute in this topology. A boosted-object tagger trained to separate W/Z/H jets from QCD will learn the jet mass as a discriminating feature: signal jets have a characteristic mass around 80–125 GeV; QCD jets have a steeply falling mass distribution. Any cut on the tagger score preferentially removes low-mass QCD background, carving a bump-shaped enhancement into the m_jj spectrum at the signal hypothesis. A bump-hunt that treats the tagger as mass-decorrelated will attribute this enhancement to a new resonance. The significance is fake.
The existing analysis pipeline uses a neural network tagger trained on eight high-level kinematic variables per fat jet: lep1_pt, lep2_pt, fatjet_pt, fatjet_eta, fatjet_D2, Zll_mass, Zll_pt, MET. Vikshep replaces this tagger with constituent-level wavelet scattering features plus the DisCo mass-decorrelation penalty. The physics reach is measured by two numbers against the same Wilks Δχ² significance pipeline the analysis already uses.
The data
Geant4-simulated ATLAS samples
The pilot runs on Geant4-simulated ATLAS detector response with full detector geometry. Signal samples cover gg → H → ZV at four resonance mass points (700, 1000, 1500, 2000 GeV). Background samples include diboson production (WW, WZ, ZZ), Z+jets (with up to four additional partons), and semi-leptonic top-quark pair production.
Each fat jet is reconstructed with the anti-k_T algorithm at R=1.0 and trimmed with f_cut=0.05. Constituent four-vectors are read directly from the ROOT TTree output of the ATLAS analysis framework. The jet image rasterisation is handled by the Vikshep ingest step at 64×64 bins in (η, φ) centred on the fat jet axis.
| sample | process | mass points | N_events |
|---|---|---|---|
| Signal | gg → H → ZV | 700, 1000, 1500, 2000 GeV | 50k per point |
| Background | Diboson (WW/WZ/ZZ) | — | 200k |
| Background | Z+jets (up to 4j) | — | 500k |
| Background | tt̄ (semi-leptonic) | — | 300k |
The benchmark
Two numbers. Apples-to-apples.
The benchmark is computed against the analysis’s own Wilks Δχ² significance pipeline — the same code the collaboration uses to quote significance in the paper. No new significance metric, no cherry-picked working point. Two numbers are reported:
Δσ — significance gain
TBD
Standard deviations above the NN baseline at the same cut efficiency
ΔJSD — sculpting reduction
TBD
Jensen–Shannon divergence pre→post cut, relative to the NN baseline
Success criterion: Vikshep achieves equal or higher Δσ at strictly lower ΔJSD. A higher-significance result that sculpts more than the NN baseline is not a success — it would trade a physics guarantee for a statistical number. The goal is to demonstrate that the two objectives are not in tension: you can have both.
Benchmark protocol
Status & results
Currently in progress.
Q4 2025
Geant4 sample preparation and ROOT TTree export
Q4 2025
Ingest pipeline and jet image rasterisation validated
Q4 2025
Scattering feature extraction on full background sample
Q1 2026
DisCo classifier training and hyperparameter sweep
Q1 2026
Wilks Δχ² benchmark run against NN baseline
Q1 2026
Results published here and on GitHub
Results expected Q1 2026. When available, the two benchmark numbers (Δσ and ΔJSD) will appear above, and the post-cut m_jj histograms for both taggers will be published as a ROOT file in the pilot directory of the GitHub repo.
Reproduce
The full harness is in the repo.
The pilot analysis harness — data loading, scattering configuration, DisCo training loop, and Wilks significance pipeline — lives in Vikshep/pilot/ in the GitHub repository. To reproduce:
- Clone the repo:
git clone https://github.com/samvardhan03/Vikshep - Install dependencies:
pip install vikshep && cargo install omnipulse-mcp - Obtain the Geant4 samples (see
pilot/README.mdfor access instructions) - Run the harness:
bun run pilot/run.ts --config pilot/uoe-diboson.yaml - Results are written to
pilot/results/as JSON + ROOT files
Acknowledgements
Built with.
Research partner
University of Edinburgh
Nuclear physics researcher (name withheld pending permission). Provided the Geant4 samples, the existing analysis notebook, and the Wilks Δχ² significance pipeline used as the benchmark baseline.
Simulation framework
ATLAS + Geant4
Detector simulation produced with the ATLAS detector Geant4 implementation. Jet reconstruction performed with the FastJet anti-k_T algorithm (R=1.0). ATLAS Open Data release used where applicable.
Theory baseline
DisCo — Kasieczka & Shih (2020)
The DisCo distance-correlation penalty is the closed-form mass-decorrelation guarantee used in the classifier training step. See arXiv:2001.05310.