我对编码非常陌生,正在为学校编写一个小脚本。它是为HiSparc项目而设计的,用于测量宇宙辐射,因此使用了一个名为sapphire的定制模块,因此我将展示的一些片段可能是新的。我正在努力创建一个图表,显示这些年来,这些检测器检测到的事件量是如何变化的(如果有的话)
这些HiSparc探测器创建了自己的表格,可下载并由Sapphire处理。这里我使用的是quick_download函数,它只下载501站昨天的数据。如果你运行这个脚本,它给出的事件总数约为54000
import sapphire
from sapphire import quick_download
data = quick_download(501)
events_table = data.root.s501.events
events = sum(1 for line in events_table)
print(events)
我想每100000秒左右得到一个数字,HiSparc数据包括GPS时间内每个事件的时间戳,所以这应该是可能的,我不知道如何分割这个图表,然后用这些事件创建一个新文件
100000秒大约是一天,如果我使用这个快速下载功能(24小时)并搜索事件0的时间戳(当然是事件1)和事件50000的时间戳(大约是一整天之后),我得到的时间戳差为:
events_table[0]['timestamp']
Out[7]: 1575936003
events_table[50000]['timestamp']
Out[8]: 1576016547
+-80000,因此,如果我每100000秒运行一次“events=”行,则预期结果大约为60000个事件。这个脚本的想法是找出这些事件的数量是否随着时间的推移而变化
start = datetime.datetime(2000, 1, 1)
end = datetime.datetime(2020, 1, 1)
esd.download_data(data, '/s501/events/timestamp', 501, start, end)
像这样一个脚本,它下载了所有在线年份的数据,可以确定。 如果您需要更多信息,请询问,非常感谢您的帮助
您好, 赛斯
目前没有回答
相关问题 更多 >
编程相关推荐