DOI: 10.5281/zenodo.807740
title: Data: Muon Tomography sites for Colombia volcanoes for generate muon flux trought volcanic structures (arXiv:1705.09884v1)

description: <strong>Data: Muon Tomography sites for Colombia volcanoes (arXiv:1705.09884v1)</strong><br> <strong>The MuTe Collaboration</strong><br> <em>Data for generating figure 6<br> Muon Tomography sites for Colombia volcanoes (arXiv:1705.09884v1)</em> 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): <em><strong>&gt; </strong></em><strong>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 &gt; salida.out</strong> 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.
description: The authors acknowledge the financial support of Departamento Administrativo de Ciencia, Tecnolog\'ia e Innovaci\'on of Colombia (ColCiencias) under contracts FP44842-051-2015 and FP44842-661-2015. We are also very thankful to LAGO and to Pierre Auger Collaboration for its continuous support.