2026-09-28 | | Total: 3
Detector simulation and event reconstruction will be computationally expensive for future $e^+e^-$ collider programs and are currently critical bottlenecks for accurate feasibility studies. To address this challenge, we build a Parnassus model for the CLD detector concept. Parnassus is a framework for automatically tuning a surrogate model, in our case, a conditional flow matching neural network, to emulate a full detector simulation and reconstruction. We train on $e^+e^-\to Z\to q\bar q$ events at $\sqrt{s}=91.2$ GeV processed through a Geant4 simulation of the CLD detector concept and the Pandora particle-flow reconstruction and reproduce single-particle kinematics, particle-IDs, and impact-parameter distributions. We also examine jet- and event-level features, inclusively and split by flavor, and find excellent fidelity, significantly better than the parameterized program Delphes. Furthermore, we train a transformer-based flavor tagger on the reconstructed particle-flow constituents and show that the surrogate preserves the $b/c/s/q$ discrimination of the full CLD reconstruction. The Parnassus CLD model achieves a generation cost of about 1.2 ms (40 ms) per event on a single GPU (CPU), over three (two) orders of magnitude faster than full simulation and reconstruction. Our model is publicly available for feasibility and design studies.
Validating whether the underlying physics of a model has been retained after a full phase-space reweighting poses a unique challenge. Often, validation relies on comparing histograms of 1D observables; however, this can mask correlations and biases in the complete prediction. We present a novel, unbinned approach to comparing the performance of such reweighting schemes based on the "Cross-Section-Mover's Distance", an application of Optimal Transport that quantifies the work required to transform one theoretical prediction into another and enables an interpretation of results in terms of metric spaces. We demonstrate its utility when benchmarking various reweighting schemes that mitigate the effects of negative weights in a Monte Carlo simulation. This approach can be broadly applied in other scenarios where biases in full phase-space reweighting schemes should be studied in an unbinned way.
As the accuracy of experimental results in high energy physics (HEP) increases, so does the demand for precision Monte Carlo (MC) simulation. Higher-accuracy event generation at next-to-leading order (NLO) and beyond brings with it negatively weighted events. These negatively weighted events reduce the statistical power of samples, increasing the number of events that need to be produced and straining already limited computational resources. We present a post-hoc reweighting scheme that employs cell-based resampling using an IRC-safe metric to define the cell radii. An optimal metric embeds kinematically similar events within the same cells, minimizing bias when reweighted. This motivates our exploration of a Optimal Transport (OT) based distance metrics. We compare the performance of the reweighting algorithm with different choices of metric, and explicitly demonstrate the performance on simulated Z+jets events produced at NLO accuracy.