Abstract: Monte Carlo methods provide detailed and accurate results for radiation transport simulations. Unfortunately, the high computational cost of these methods limits its usage in real-time applications. Moreover, existing computer codes do not provide a methodology for adapting these kinds of simulations to specific problems without advanced knowledge of the corresponding code system, and this restricts their applicability. To help solve these current limitations, we present PenRed, a general-purpose, standalone, extensible and modular framework code based on PENELOPE for parallel Monte Carlo simulations of electron-photon transport through matter. It has been implemented in C++ programming language and takes advantage of modern object-oriented technologies. In addition, PenRed offers the capability to read and process DICOM images as well as to construct and simulate image-based voxelized geometries, so as to facilitate its usage in medical applications. Our framework has been successfully verified against the original PENELOPE Fortran code. Furthermore, the implemented parallelism has been tested showing a significant improvement in the simulation time without any loss in precision of results. Program summary Program title: PenRed: Parallel Engine for Radiation Energy Deposition. CPC Library link to program files: https://doi .org /10 .17632/rkw6tvtngy.1 Licensing provision: GNU Affero General Public License (AGPL). Programming language: C++ standard 2011. Nature of problem: Monte Carlo simulations usually require a huge amount of computation time to achieve low statistical uncertainties. In addition, many applications necessitate particular characteristics or the extraction of specific quantities from the simulation. However, most available Monte Carlo codes do not provide an efficient parallel and truly modular structure which allows users to easily customise their code to suit their needs without an in-depth knowledge of the code system. Solution method: PenRed is a fully parallel, modular and customizable framework for Monte Carlo simulations of the passage of radiation through matter. It is based on the PENELOPE [1] code system, from which inherits its unique physics models and tracking algorithms for charged particles. PenRed has been coded in C++ following an object-oriented programming paradigm restricted to the C++11 standard. Our engine implements parallelism via a double approach: on the one hand, by using standard C++ threads for shared memory, improving the access and usage of the memory, and, on the other hand, via the MPI standard for distributed memory infrastructures. Notice that both kinds of parallelism can be combined together in the same simulation. Moreover, both threads and MPI processes, can be balanced using the builtin load balance system (RUPER-LB [30]) to maximise the performance on heterogeneous infrastructures. In addition, PenRed provides a modular structure with methods designed to easily extend its functionality. Thus, users can create their own independent modules to adapt our engine to their needs without changing the original modules. Furthermore, user extensions will take advantage of the builtin parallelism without any extra effort or knowledge of parallel programming. Additional comments including restrictions and unusual features: PenRed has been compiled in linux systems withg++ of GCC versions 4.8.5, 7.3.1, 8.3.1 and 9; clang version 3.4.2 and intel C++ compiler (icc) version 19.0.5.281. Since it is a C++11-standard compliant code, PenRed should be able to compile with any compiler with C++11 support. In addition, if the code is compiled without MPI support, it does not require any non standard library. To enable MPI capabilities, the user needs to install whatever available MPI implementation, such as openMPI [24] or mpich [25], which can be found in the repositories of any linux distribution. Finally, to provide DICOM processing support, PenRed can be optionally compiled using the dicom toolkit (dcmtk) [32] library. Thus, PenRed has only two optional dependencies, an MPI implementation and the dcmtk library.

Abstract: We have constructed a GEANT4-based detailed software model of photon transport in plastic sontillator blocks and have used it to study the NEMO-3 and SuperNEMO calorimeters employed in experiments designed to search for neutnnoless double beta decay We compare our simulations to measurements using conversion electrons from a calibration source of (BI)-B-207 and show that the agreement is improved if wavelength-dependent properties of the calorimeter are taken into account In this article we briefly describe our modeling approach and results of our studies.

Abstract: We present a search for neutrino induced events containing a single, exclusive photon using data from the NOMAD experiment at the CERN SPS where the average energy of the neutrino flux is similar or equal to 25 GeV. The search is motivated by an excess of electron-like events in the 200-475 MeV energy region as reported by the MiniBooNE experiment. In NOMAD, photons are identified via their conversion to e(+)e(-) in an active target embedded in a magnetic field. The background to the single photon signal is dominated by the asymmetric decay of neutral pions produced either in a coherent neutrino-nucleus interaction, or in a neutrino-nucleon neutral current deep inelastic scattering, or in an interaction occurring outside the fiducial volume. All three backgrounds are determined in situ using control data samples prior to opening the 'signal-box'. In the signal region, we observe 155 events with a predicted background of 129.2 +/- 8.5 +/- 3.3. We interpret this as null evidence for excess of single photon events, and set a limit. Assuming that the hypothetical single photon has a momentum distribution similar to that of a photon from the coherent pi(0) decay, the measurement yields an upper limit on single photon events, < 4.0 x 10(-4) per nu(mu) charged current event. Narrowing the search to events where the photon is approximately collinear with the incident neutrino, we observe 78 events with a predicted background of 76.6 +/- 4.9 +/- 1.9 yielding a more stringent upper limit, < 1.6 x 10(-4) per nu(mu) charged current event.