%20(2).svg)
项目简介
本项目是基于C++的粒子物理模拟工具，用于对粒子事件进行物理参数的模拟和响应映射。借助蒙特卡洛模拟数据，可精确模拟粒子的位置、能量、动量等参数并生成输出文件，适用于粒子物理实验中的数据分析与模拟计算。
项目的主要特性和功能
- 粒子物理参数模拟：模拟粒子在探测器中的位置、能量、动量等物理参数，生成对应输出文件。
- 响应映射：利用输入参数和响应映射文件，对粒子事件做“模糊”处理，模拟真实探测器响应。
- 灵活的参数设置：用户能自定义探测器的覆盖率、半径等参数，满足不同实验需求。
- 多种输出格式：支持生成ROOT文件和扁平化的ntuple文件，方便后续数据分析。
- 易于扩展：项目结构清晰，可按需添加新功能模块或修改现有模块。
安装使用步骤
假设用户已下载本项目的源码文件，可按以下步骤安装和使用：
1. 链接蒙特卡洛数据目录：将蒙特卡洛数据目录链接到项目的data/montecarlo目录下。
shell
ln -s $(realpath /path/to/truedaq/) data/montecarlo
2. 编译和安装：使用CMake进行编译和安装，可选择安装到系统目录或用户本地目录。
- 安装到系统目录：
shell
cmake . -Bbuild
cmake --build build --target install
- 安装到用户本地目录：
shell
cmake . -Bbuild -DCMAKE_INSTALL_PREFIX=~/.local
cmake --build build --target install
3. 使用核心功能：在C++代码中引入FitterResponseMap.h和PMTGenerator.h头文件，定义探测器参数并生成模拟事件。
```c++
#include
// 定义探测器参数 std::string target = "wbls_3%"; double coverage = 0.20; unsigned radius = 5700;
// 读取探测器响应 FitterResponseMap frp(target, coverage, radius);
// 应用响应映射到事件
for(auto event : event_list) {
ff.GenerateEvent(mcx, mcy, mcz, mcu, mcw, mcv, energy, time);
// 获取拟合后的信息
double x = ff.position_x;
double y = ff.position_y;
double z = ff.position_z;
double u = ff.direction_x;
double v = ff.direction_y;
double w = ff.direction_z;
double energy = ff.energy;
}
4. **运行示例脚本**：项目提供多个示例脚本，可按需运行生成不同输出文件。
- `frprat.cc`：对输入文件应用响应映射，生成`output.root`文件。shell
frprat input.root
- `pathb.cc`：支持从蒙特卡洛数据库生成事件或直接处理输入文件。shell
pathb -t wbls_3pct -c 0.20 -r 6700 -s "reactorSignal" -o signal.root
- `astaire.cc`：生成扁平化的ntuple文件。shell
astaire -t wbls_3pct -c 0.20 -r 6700 -s "reactorSignal" -o signal.ntuple.root
```
注意事项：运行脚本前，建议查看脚本说明和文档，了解参数设置方法和作用。项目代码可能需根据具体环境和需求调整优化。
下载地址
点击下载 【提取码: 4003】【解压密码: www.makuang.net】