End-to-end simulations of the MUon RAdiography of VESuvius experiment

The MUon RAdiography of VESuvius (MURAVES) project aims at the study of the summital cone of Mt. Vesuvius, an active volcano near Naples (Italy), by measuring its density profile through muon flux attenuation. Its data, combined with those from gravimetric and seismic measurement campaigns, will be used for better defining the volcanic plug at the bottom of the crater. We report on the development of an end-to-end simulation framework, in order to perform accurate investigations of the effects of the experimental constraints and to compare simulations, under various model hypotheses, with the actual observations. The detector simulation setup is developed using GEANT4 and a study of cosmic particle generators has been conducted to identify the most suitable one for our simulation framework. To mimic the real data, GEANT4 raw hits are converted to clusters through a simulated digitization: energy deposits are first summed per scintillator bar, and then converted to number of photoelectrons with a data-driven procedure. This is followed by the same clustering algorithm and same tracking code as in real data. We also report on the study of muon transport through rock using PUMAS and GEANT4. In this paper we elaborate on the rationale for our technical choices, including trade-off between speed and accuracy. The developments reported here are of general interest in muon radiography and can be applied in similar cases.

The design and sensitivity of JUNO’s scintillator radiopurity pre-detector OSIRIS

The OSIRIS detector is a subsystem of the liquid scintillator filling chain of the JUNO reactor neutrino experiment. Its purpose is to validate the radiopurity of the scintillator to assure that all components of the JUNO scintillator system work to specifications and only neutrino-grade scintillator is filled into the JUNO Central Detector. The aspired sensitivity level of $$10^{-16}\hbox { g/g}$$ 10 - 16 g/g of $$^{238}\hbox {U}$$ 238 U and $$^{232}\hbox {Th}$$ 232 Th requires a large ($$\sim 20\,\hbox {m}^3$$ ∼ 20 m 3 ) detection volume and ultralow background levels. The present paper reports on the design and major components of the OSIRIS detector, the detector simulation as well as the measuring strategies foreseen and the sensitivity levels to U/Th that can be reached in this setup.

Full Detector Simulation with Unprecedented Background Occupancy at a Muon Collider

In recent years, a Muon collider has attracted a lot of interest in the high-energy physics community, thanks to its ability of achieving clean interaction signatures at multi-TeV collision energies in the most cost-effective way. Estimation of the physics potential of such an experiment must take into account the impact of beam-induced background on the detector performance, which has to be carefully evaluated using full detector simulation. Tracing of all the background particles entering the detector region in a single bunch crossing is out of reach for any realistic computing facility due to the unprecedented number of such particles. To make it feasible a number of optimisations have been applied to the detector simulation workflow. This contribution presents an overview of the main characteristics of the beam-induced background at a Muon collider, the detector technologies considered for the experiment and how they are taken into account to strongly reduce the number of irrelevant computations performed during the detector simulation. Special attention is dedicated to the optimisation of track reconstruction with the conformal tracking algorithm in this high-occupancy environment, which is the most computationally demanding part of event reconstruction.

The cosmic muon and detector simulation framework of the extreme energy events (EEE) experiment

This paper describes the simulation framework of the extreme energy events (EEE) experiment. EEE is a network of cosmic muon trackers, each made of three multi-gap resistive plate chambers (MRPC), able to precisely measure the absolute muon crossing time and the muon integrated angular flux at the ground level. The response of a single MRPC and the combination of three chambers have been implemented in a GEANT4-based framework (GEMC) to study the telescope response. The detector geometry, as well as details about the surrounding materials and the location of the telescopes have been included in the simulations in order to realistically reproduce the experimental set-up of each telescope. A model based on the latest parametrization of the cosmic muon flux has been used to generate single muon events. After validating the framework by comparing simulations to selected EEE telescope data, it has been used to determine detector parameters not accessible by analysing experimental data only, such as detection efficiency, angular and spatial resolution.

Simplified fast detector simulation in MadAnalysis 5

We introduce a new simplified fast detector simulator in the MadAnalysis 5 platform. The Python-like interpreter of the programme has been augmented by new commands allowing for a detector parametrisation through smearing and efficiency functions. On run time, an associated C++ code is automatically generated and executed to produce reconstructed-level events. In addition, we have extended the MadAnalysis 5 recasting infrastructure to support our detector emulator, and we provide predefined LHC detector configurations. We have compared predictions obtained with our approach to those resulting from the usage of the Delphes 3 software, both for Standard Model processes and a few new physics signals. Results generally agree to a level of about 10% or better, the largest differences in the predictions stemming from the different strategies that are followed to model specific detector effects. Equipped with these new functionalities, MadAnalysis 5 now offers a new user-friendly way to include detector effects when analysing collider events, the simulation of the detector and the analysis being both handled either through a set of intuitive Python commands or directly within the C++ core of the platform.

OASIS: Optimal Analysis-Specific Importance Sampling for event generation

We propose a technique called Optimal Analysis-Specific Importance Sampling (OASIS) to reduce the number of simulated events required for a high-energy experimental analysis to reach a target sensitivity. We provide recipes to obtain the optimal sampling distributions which preferentially focus the event generation on the regions of phase space with high utility to the experimental analyses. OASIS leads to a conservation of resources at all stages of the Monte Carlo pipeline, including full-detector simulation, and is complementary to approaches which seek to speed-up the simulation pipeline.

GeantV

Full detector simulation was among the largest CPU consumers in all CERN experiment software stacks for the first two runs of the Large Hadron Collider. In the early 2010s, it was projected that simulation demands would scale linearly with increasing luminosity, with only partial compensation from increasing computing resources. The extension of fast simulation approaches to cover more use cases that represent a larger fraction of the simulation budget is only part of the solution, because of intrinsic precision limitations. The remainder corresponds to speeding up the simulation software by several factors, which is not achievable by just applying simple optimizations to the current code base. In this context, the GeantV R&D project was launched, aiming to redesign the legacy particle transport code in order to benefit from features of fine-grained parallelism, including vectorization and increased locality of both instruction and data. This paper provides an extensive presentation of the results and achievements of this R&D project, as well as the conclusions and lessons learned from the beta version prototype.

Studies of GEANT4 performance for different ATLAS detector geometries and code compilation methods

Full detector simulation is known to consume a large proportion of computing resources available to the LHC experiments, and reducing time consumed by simulation will allow for more profound physics studies. There are many avenues to exploit, and in this work we investigate those that do not require changes in the GEANT4 simulation suite. In this study, several factors affecting the full GEANT4 simulation execution time are investigated. A broad range of configurations has been tested to ensure consistency of physical results. The effect of a single dynamic library GEANT4 build type has been investigated and the impact of different primary particles at different energies has been evaluated using GDML and GeoModel geometries. Some configurations have an impact on the physics results and are, therefore, excluded from further analysis. Usage of the single dynamic library is shown to increase execution time and does not represent a viable option for optimization. Lastly, the static build type is confirmed as the most effective method to reduce the simulation execution time.