Skip to content

adrianacvr/Simulacion-barra-MuTe

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

10 Commits
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

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:

Compiling the code

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.

Running the program "bar"

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

Modifying the input variables

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 cm

    The parameter z may 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 1

    Note: 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.

Data visualization of the results

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.

About

Scintillator bar simulation for the Muon Telescope MuTe using Geant4.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published