Data: Muon Tomography sites for Colombia volcanoes for generate muon flux trought volcanic structures (arXiv:1705.09884v1)

Data: Muon Tomography sites for Colombia volcanoes (arXiv:1705.09884v1)
The MuTe Collaboration
Data for generating figure 6
Muon Tomography sites for Colombia volcanoes (arXiv:1705.09884v1)

The files in this record contain data from Extensive Atmospheric Shower simulations made by CORSIKA and Magnetocosmic codes, in a muography of Machin Volcano (Colombia) [https://volcano.si.edu/volcano.cfm?vn=351040] for a particular observation point. The objective was to count muons crossing the volcanic structure for a fixed observation point. Muon transport through the volcanic edifice is calculated by using an algorithm taken Corsika-Magnetocosmic output and taking into account the energy losses with the muon stopping power  tables given by Particle Data Group (PDG).

This dataset contains:

• Six (6) Corsika output files of 4 hours of simulation each using a flat detector (equivalent to 12 hours of simulated cosmic rays).
• One (1) Corsika output file of 12 hours of simulation using flat detector (equivalent to 6 hours of simulated cosmic rays).
• One (1) Corsika output file of 48 hours of simulation using flat detector (equivalent to 24 hours of simulated cosmic rays).
• Two (2) Corsika output files of 24 hours of simulation each using volumetric detector (equivalent to 48 hours of simulated cosmic rays).
• Everything makes a total simulated time of 3.75 days.

The output files necessary for the determination of the muon flux through the volcanic structure are obtained through the following process:

• From the .shw.bz2 files it is possible to obtain an output file with the momentum information of the particle in the x, y, and z directions, and also the total momentum of the muons, essentially a formatted file (px, py, pz, p). This file can be built by typing in a terminal shell (bash code):

> bzcat *.shw.bz2 | awk '{if($1==0006 ||$1==0005){j=sqrt(($2*$2)+($3*$3)+($4*$4));printf "%s %s %s %.s\n",$2, $3, $4, j }}' | sort -n > salida.out

• Metadata in the showers file is as this type (for example):

# # # shw

# # CURVED mode is ENABLED and observation level is 2750 m a.s.l.

# # This is the Secondaries file - CrkTools v3r0

# # 12 column format is:

# # CorsikaId px py pz x y z shower_id prm_id prm_energy prm_theta prm_phi

0001 +1.42146e-04 -3.96008e-05 +1.60247e-04 -1.31716e+03 -6.10051e+01 +2.44986e+03 00000001 0703 +1.25065e+02 +43.016 +021.768

0003 +1.80713e-04 +1.89560e-03 +4.49373e-03 -1.32279e+03 -5.62869e+01 +2.44986e+03 00000001 0703 +1.25065e+02 +43.016 +021.768

0003 +9.41713e-03 +2.38845e-03 +1.10020e-02 -1.32098e+03 -5.68323e+01 +2.44986e+03 00000001 0703 +1.25065e+02 +43.016 +021.768

• Concatenate all output files.
• Then, the muon flux trought rock can be calculated from two python codes, available in https://github.com/AstroparticulasBucaramanga/Propagacion-Muones-en-Roca. This step generates the final files to be graphed with any plotter, in our case, also using python.