Efficient cosmic ray generator for particle detector simulations

TL;DR

This paper introduces a more efficient Monte Carlo method for simulating cosmic rays hitting detectors, significantly reducing computation time while maintaining accuracy, especially for non-spherical target geometries.

Contribution

A novel secant method for cosmic ray simulation that improves efficiency by a factor of 37 over traditional approaches, applicable to complex detector geometries.

Findings

Method achieves 37x faster computation

Maintains identical physical results

Effective for non-spherical target volumes

Abstract

Traditional cosmic ray simulations make use of the Montecarlo method in a very naive way to randomise energy and direction for each simulated particle. The flux of cosmic rays is modelled as a rain coming from a plane above the object of interest (detectors in particle physics applications, planes in dosimetry studies, etc.) with an experimental angular and energy distributions. This strategy is very inefficient because many of the particles never touch the detector. Here a refined way of implementing the Montecarlo method is proposed in order to generate a sample of events that hit the target volume whose angular distribution coincides with the one from the naive implementation. It is based on the projection of a sphere containing the target volume onto a plane tangent to it with a fixed angle, we call it the secant method. This configuration allows to compute the probability of a cosmic particle hitting the sphere with this incoming angle as proportional to the area of the corresponding section of a cylinder. The performance of this method is faster in terms of computing time and identical physical results are achieved. It has been implemented in REST-for-Physics framework and it is tested with the geometry of a real detector, the IAXO-D0 Micromegas X-ray detector for the future axion helioscope BabyIAXO. Our method is 37 times more efficient than the traditional Montecarlo schema for the same accuracy, being more useful when the target volume departs from spherical shape