Scintillator Bar Simulation for the Muon Telescope MuTe
Author: Adriana Vásquez Ramírez adrianacvr67@gmail.com
Results published in A. Vásquez-Ramírez et al. 2020, JINST 15 P08004, DOI: 10.1088/1748-0221/15/08/P08004
This repository provides a simulation of the response of the scintillator bars of the Muon Telescope (MuTe) to the incidence of muons at different positions along the bar. The implementation is based on the WLS (WaveLength Shifter) example from GEANT4 version 10.3.2.
To run the simulation, download and extract this repository and proceed as follows:
Compilation requires creating a build directory outside the extracted repository and executing the following commands:
-
mkdir name-build
-
cd name-build
-
cmake -DGeant4_DIR=/home/adriana/GEANT4/geant4.10.3-install/lib/Geant4-10.3.1/ ../Simulacion-barra-MuTe-master/
Note: Here, -DGeant4_DIR must be set to the installation path of Geant4. The path shown corresponds to the local path on my system and is provided only as an example.
- make -j2
This specifies the use of 2 cores; alternatively, adjust the value to match the number of CPU cores available on the system.
Upon successful completion (100% build), the executable bar will be generated.
To launch the Geant4 visualizer:
- ./bar
To execute the simulation using the parameters defined in input.in (recommended for a large number of particles):
- ./bar -m input.in
The configuration file input.in allows specification of the particle type, energy, and initial position, for example:
-
Particle type:
/gun/particle mu- -
Energy:
/gun/energy 1.0 GeV -
Position:
/gun/position 1 0 z cmThe parameter
zmay range from -60 cm to 60 cm to investigate attenuation effects as a function of the distance between the SiPM and the particle impact position. The bar is centered at (0,0,0) cm. -
Number of events:
/run/beamOn 1Note: Try first with run 1 or 10 particles to verify that the code executes correctly. To reproduce the results presented in the paper/thesis, the simulation must be run with 10000 particles.
After configuring the input parameters and executing ./bar -m input.in, the output file scintbar.root will be produced, containing the following histograms:
-
scint_photons: Histogram of the number of scintillation photons generated in the bar. -
sipm_photons: Histogram of the number of photons reaching the SiPM surface. -
photo_elec: Histogram of the number of photoelectrons generated according to the SiPM quantum efficiency.
You can access these histograms using CERN’s ROOT framework.
For further details on these results and the simulation, see the full thesis in the UIS Virtual Library.
Please request the file Document.pdf available at that link.