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

<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.<strong>数据集：哥伦比亚火山缪子（Muon）层析成像站点（arXiv:1705.09884v1）</strong><br> <strong>MuTe合作组</strong><br> <em>用于生成图6的数据：哥伦比亚火山缪子（Muon）层析成像站点（arXiv:1705.09884v1）</em> 本数据集收录的文件包含由CORSIKA与Magnetocosmic程序开展的广延大气簇射模拟数据，相关模拟针对哥伦比亚马钦火山（Machin Volcano）的某一观测点开展缪子层析成像研究[https://volcano.si.edu/volcano.cfm?vn=351040]。本次研究的目标为统计穿过该火山构造的缪子数目，观测点固定不变。缪子穿过火山岩体的输运过程通过基于CORSIKA-Magnetocosmic输出结果构建的算法进行计算，同时参考粒子数据组（Particle Data Group, PDG）提供的缪子阻止本领表格考量能量损失效应。<br><br>本数据集包含如下内容：<br>6份平坦探测器单次模拟时长为4小时的CORSIKA输出文件（等效于12小时的宇宙线模拟时长）；<br>1份平坦探测器单次模拟时长为12小时的CORSIKA输出文件（等效于6小时的宇宙线模拟时长）；<br>1份平坦探测器单次模拟时长为48小时的CORSIKA输出文件（等效于24小时的宇宙线模拟时长）；<br>2份体探测器单次模拟时长为24小时的CORSIKA输出文件（等效于48小时的宇宙线模拟时长）。<br>所有文件累计模拟时长总计3.75天。<br><br>用于计算穿过火山构造的缪子通量所需的输出文件可通过如下流程获取：从.shw.bz2格式文件中，可提取得到包含粒子x、y、z方向动量信息以及缪子总动量的输出文件，本质为格式化后的(px, py, pz, p)格式文件。可通过终端Shell（Bash代码）执行如下命令构建该文件：<br><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 ",$2, $3, $4, j }}' | sort -n &gt; salida.out</strong><br><br>簇射文件的元数据格式示例如下：<br><pre># # # shw # # 已启用CURVED模式，观测高程为2750米海拔（a.s.l.）<br># # 这是次级粒子文件——CrkTools v3r0<br># # 共12列数据，格式如下：<br># # CORSIKA编号 px py pz x y z 簇射ID 初始粒子ID 初始粒子能量 初始粒子天顶角 初始粒子方位角<br>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<br>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<br>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</pre><br><br>将所有输出文件合并后，可通过https://github.com/AstroparticulasBucaramanga/Propagacion-Muones-en-Roca 提供的两段Python代码计算穿过岩体的缪子通量。该步骤将生成可通过任意绘图软件可视化的最终文件，本次研究同样使用Python完成绘图流程。