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 for simulate Extensive Atmospheric Showers used for the study of muon flux trought in the Machin Volcano (Colombia) [https://es.wikipedia.org/wiki/Cerro_Machín] for a especific observation point. The objective is count muons that survive after crossed an amount of distance d through volcanic structure for a fixed observation point with the use of a muon detector (Muon Telescope). The individual muons are passed through rock using an algorithm designed for that purpose, which takes into account the energy losses for a given material composition of Particle Data Group (PDG) muon stopping power tables and Corsika simulations outputs are used.

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:
1) 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 construdted. This file can be built by typing in a terminal shell
> 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

2) Concatenate all output files.
3) 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.