# PyMagicTools

**Repository Path**: zhangzj-gitee/PyMagicTools

## Basic Information

- **Project Name**: PyMagicTools
- **Description**: Particle-In-Cell软件MAGIC的相关工具
- **Primary Language**: Unknown
- **License**: Not specified
- **Default Branch**: master
- **Homepage**: None
- **GVP Project**: No

## Statistics

- **Stars**: 0
- **Forks**: 0
- **Created**: 2023-06-28
- **Last Updated**: 2024-10-11

## Categories & Tags

**Categories**: Uncategorized

**Tags**: 后处理, Magic

## README

# MAGIC后处理工具
## 快速入门
不同后缀的文件用对应的parser解析即可。例如，解析par，应采用par_parser.py。

total_parser.py可以对同一.m2d文件生成的几种结果文件进行整合处理。

下面的例子展示了如何解析一个.grd文件（注意，代码中的文件路径、数据标题等字符串请根据实际情况自行替换）：

```python
import grd_parser
grd = grd_parser.GRD(r'path/to/data.grd')

# How to use a GRD object
print(grd.obs.keys()) # 查看此.grd文件中包含的observe数据标题

print(grd.ranges[list(grd.ranges.keys())[0]][0]['data'] # 查看此.grd文件中包含的第0块range数据

```

数据处理大多时候是个走一步看一步的过程。例如，在上面的例子中，一开始我们并不知道grd对象中都有哪些数据标题，因此选择了先将所有数据标题打印出来。然而，读取和解析文件（如上例中创建grd对象的过程）通常耗时较长，要尽可能避免重复执行这些语句。

出于上述原因，推荐使用PyCharm的控制台，交互地执行后处理代码（通过在编辑界面上右键-`Run File in Python Console`）,执行之后效果如下图所示。从下面截图的右侧红框部分可以看出，grd.obs是一个dict对象，其中包括`FIELD_INTEGRAL ...` 等一系列key（即，前文提到的`数据标题`）。

![388d67f0aaddbda7172e2ab6196ed14d.png](.md_attachments/388d67f0aaddbda7172e2ab6196ed14d.png)

## 兼容性

暂时只支持Magic 2005 Single。

double版本及Magic 2017版运行结果中会用1.0+12这种形式来表示浮点数1.0E+12，暂时无法解析。


## 效果展示

### 输出组图到指定文件夹
相关代码见 [此文件](total_parser_2.py)。


![7df2b75a6e930a5a86a283d2ccd00c74.png](.md_attachments/7df2b75a6e930a5a86a283d2ccd00c74.png)
![e5b8f6ee576b1cba60f8b976fea1754e.png](.md_attachments/e5b8f6ee576b1cba60f8b976fea1754e.png)

paper_plot/trivial.py输出精修图：

![d41d8953c378defb8069edf251030904.png](.md_attachments/d41d8953c378defb8069edf251030904.png)

![678f635c7c15cbeb7722e0baea3e4581.png](.md_attachments/678f635c7c15cbeb7722e0baea3e4581.png)
更多数据展示方式：
![a1cb5069d056212f136e1e048782e834.png](.md_attachments/a1cb5069d056212f136e1e048782e834.png)

![238b54c8d91709e4b6038700eaf6b8d4.png](.md_attachments/238b54c8d91709e4b6038700eaf6b8d4.png)
![a072f7196d44eb93e270a42f12203c37.png](.md_attachments/a072f7196d44eb93e270a42f12203c37.png)


### 绘制几何结构
最简代码：

```python
from geom_parser import GEOM
import matplotlib.pyplot as plt
# 文件名后缀无所谓，但必须要有
filename = r"D:\MagicFiles\CherenkovAcc\cascade\min_case_for_gradient_test\test_diffraction-23.grd"
geom = GEOM(filename)
plt.figure()
geom.plot(plt.gca())
```
![24f98c60e811c85ee045b70eef7fd78b.png](.md_attachments/24f98c60e811c85ee045b70eef7fd78b.png)
![2f1e24d84df5afcd85ada993bc97f364.png](.md_attachments/2f1e24d84df5afcd85ada993bc97f364.png)
下面的代码段显示了如何将粒子与几何结构绘制在一张图上：
```python
plt.figure()
geom.plot(plt.gca())
plt.scatter(*par.phasespaces[' ALL PARTICLES @AXES(X1,X2)-#1 $$$PLANE_X1_AND_X2_AT_X0=  0.000'][-5]['data'].values.T,s = 0.0001,c ='r')
```
![aac83820945a4b7f782e547b341c4108.png](.md_attachments/aac83820945a4b7f782e547b341c4108.png)