GEANT4学习 输出粒子信息小计

最新推荐文章于 2022-01-20 15:11:57 发布
最新推荐文章于 2022-01-20 15:11:57 发布
阅读量5.3k 收藏 32
点赞数 7
分类专栏： GEANT4学习
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_42449640/article/details/106318561
版权

GEANT4学习 输出粒子信息小计

写了一个简单的模型，并在step里面输出粒子的信息，通过这些验证G4在粒子输出信息上的规则。

我原本想的是输出
eventID、trackID、stepID、particle name、prepoint process name、postpoint process name、creator name、prepoint Kinetic Energy、postpoint Kinetic Energy、prepoint volume name、postpoint volume name以及energy deposit。

第一次尝试的时候出现了问题，个别的信息输出时会遇到“段错误”，网上解释段错误是栈溢出，但是我感觉好多都是信息的逻辑提取出现了问题，这个我之后讲。

第一次输出信息的代码如下所示：

  G4int eventID = (G4EventManager::GetEventManager())->GetConstCurrentEvent()->GetEventID();
  // 因为在step里面，我使用EventManager这个类去获取EventID
  G4int trackID = step->GetTrack()->GetTrackID();
 
  G4String particleName =  step->GetTrack()->GetDefinition()->GetParticleName();
  
  G4String postProcessName = step->GetPostStepPoint()->GetProcessDefinedStep()->GetProcessName();

  G4String preName = step->GetPreStepPoint()->GetTouchableHandle()->GetVolume()->GetName();

  G4double postKE = step->GetPostStepPoint()->GetKineticEnergy();
  G4double preKE = step->GetPreStepPoint()->GetKineticEnergy();

  G4double edepStep = step->GetTotalEnergyDeposit();

输出语句如下所示：

 outfileInfo<<std::setw(10)<<std::setfill(' ') <<std::left <<particleName                           
                      <<std::setw(4)<<std::setfill(' ') <<std::left <<eventID
                      <<std::setw(4)<<std::setfill(' ') <<std::left <<trackID
                      
                     <<std::setw(8)<<std::setfill(' ') <<std::left <<preName
                    <<std::setw(20)<<std::setfill(' ') <<std::left <<preKE

                   <<std::setw(20)<<std::setfill(' ') <<std::left <<postProcessName
          
                  <<std::setw(20)<<std::setfill(' ') <<std::left <<postKE

                  <<std::setw(20)<<std::setfill(' ') <<std::left <<edepStep
                  <<G4endl;

输出信息的文件如下所示：

particleN eID tID preN    preKE               postProName         postKE              edepStep
neutron   0   1   World   14                  Transportation      14                  0                   
neutron   0   1   A1      14                  Transportation      14                  0                   
neutron   0   1   A2      14                  hadElastic          13.5434             0                   
neutron   0   1   A2      13.5434             hadElastic          9.21031             0                   
neutron   0   1   A2      9.21031             hadElastic          5.6599              0                   
neutron   0   1   A2      5.6599              hadElastic          3.56607             0                   
neutron   0   1   A2      3.56607             hadElastic          2.93008             0                   
neutron   0   1   A2      2.93008             Transportation      2.93008             0                   
neutron   0   1   World   2.93008             Transportation      2.93008             0                   
O16       0   6   A2      0.63616             ionIoni             0                   0.63616             
O16       0   6   A2      0                   RadioactiveDecay    0                   0                   
proton    0   5   A2      2.09497             hIoni               1.38174             0.713232            
proton    0   5   A2      1.38174             hIoni               0.580998            0.800742            
proton    0   5   A2      0.580998            hIoni               0.0178738           0.563124            
proton    0   5   A2      0.0178738           hIoni               0.0114283           0.00644545          
proton    0   5   A2      0.0114283           hIoni               0.00879221          0.00263613          
proton    0   5   A2      0.00879221          hIoni               0                   0.00879221          
proton    0   4   A2      3.55233             hIoni               2.98914             0.563191            
proton    0   4   A2      2.98914             hIoni               2.34926             0.639877            
proton    0   4   A2      2.34926             hIoni               1.67271             0.676547            
proton    0   4   A2      1.67271             hIoni               0.877513            0.7952              
proton    0   4   A2      0.877513            hIoni               0.0726702           0.804843            
proton    0   4   A2      0.0726702           hIoni               0                   0.0726702           
proton    0   3   A2      4.33546             hIoni               3.72805             0.607402            
proton    0   3   A2      3.72805             hIoni               3.15269             0.575363            
proton    0   3   A2      3.15269             hIoni               2.55276             0.599933            
proton    0   3   A2      2.55276             hIoni               1.88505             0.667708            
proton    0   3   A2      1.88505             hIoni               1.15731             0.727738            
proton    0   3   A2      1.15731             hIoni               0.0440143           1.1133              
proton    0   3   A2      0.0440143           hIoni               0                   0.0440143           
proton    0   2   A2      0.456834            hIoni               0.0206327           0.436202            
proton    0   2   A2      0.0206327           hIoni               0.0198335           0.00079924          
proton    0   2   A2      0.0198335           hIoni               0.0121395           0.007694            
proton    0   2   A2      0.0121395           hIoni               0.00741999          0.00471949          
proton    0   2   A2      0.00741999          hIoni               0                   0.00741999

我想输出的信息有很多还没办法实现，接下来来一个一个讲：
stepID：同学说没有stepID这个讲法，官网上好像确实没有相关的语句；

preProcessName:

G4String preProcessName = step->GetPreStepPoint()->GetProcessDefinedStep()->GetProcessName();

这段一放上去就会出现段错误，我判断是初始粒子的preProcessName是没有的，因为是初始粒子，由初始粒子生成器生成的。

因此我加了一句来判断是否是初始粒子：

G4String preProcessName;
 if (trackID==1)  preProcessName="primary";
 else preProcessName = step->GetPreStepPoint()->GetProcessDefinedStep()->GetProcessName();

但是还是会出现“段错误”，有可能这个prePoint就没有processname。
所以在输出的时候我又换成下面的语句。对于一个step 前点使用creatorProcessName，后点使用postPointProcessName。

G4String creatorName;
if (trackID==1)  creatorName="primary";
else creatorName = step->GetTrack()->GetCreatorProcess()->GetProcessName();

postPointVolumName：
后点volume的名字出现问题，我觉得是当粒子输运接触world边界时，无法获取后点的volume，因此我也加一个判断：

G4String postName;
  if (step->GetPostStepPoint()->GetStepStatus() == fWorldBoundary) postName="out";
  else postName = step->GetPostStepPoint()->GetTouchableHandle()->GetVolume()->GetName();

这样就可以把想要的信息输出了，摘录部分输出结果：

neutron   6   1   primary             World   14                  Transportation      A1     14                  0                   
neutron   6   1   primary             A1      14                  Transportation      A2     14                  0                   
neutron   6   1   primary             A2      14                  neutronInelastic    A2     0                   0                   
gamma     6   7   neutronInelastic    A2      0.0210016           phot                A2     0                   2.834e-05           
e-        6   9   phot                A2      0.00050894          msc                 A2     0.00050894          0                   
e-        6   9   phot                A2      0.00050894          eIoni               A2     0                   0.00050894          
e-        6   8   phot                A2      0.0204643           msc                 A2     0.000334888         0.0201294           
e-        6   8   phot                A2      0.000334888         eIoni               A2     0                   0.000334888         
gamma     6   6   neutronInelastic    A2      0.0210016           phot                A2     0                   4.338e-05           
e-        6   11  phot                A2      0.0004939           msc                 A2     0.000458443         3.54571e-05         
e-        6   11  phot                A2      0.000458443         eIoni               A2     0                   0.000458443         
e-        6   10  phot                A2      0.0204643           msc                 A2     0.0190617           0.00140257          
e-        6   10  phot                A2      0.0190617           eIoni               A2     0                   0.0190617           
gamma     6   5   neutronInelastic    A2      0.00267715          phot                A2     0                   5.846e-05           
e-        6   13  phot                A2      0.00047882          msc                 A2     0.00047882          0                   
e-        6   13  phot                A2      0.00047882          eIoni               A2     0                   0.00047882          
e-        6   12  phot                A2      0.00213987          msc                 A2     0.00186088          0.000278987         
e-        6   12  phot                A2      0.00186088          eIoni               A2     0                   0.00186088          
gamma     6   4   neutronInelastic    A2      3.85317             Transportation      World  3.85317             0                   
gamma     6   4   neutronInelastic    World   3.85317             Transportation      out    3.85317             0                   
C13       6   3   neutronInelastic    A2      1.38632             ionIoni             A2     0.272673            1.11365             
C13       6   3   neutronInelastic    A2      0.272673            ionIoni             A2     0                   0.272673            
C13       6   3   neutronInelastic    A2      0                   RadioactiveDecay    A2     0                   0                   
alpha     6   2   neutronInelastic    A2      6.5013              ionIoni             A2     3.15038             3.35092             
alpha     6   2   neutronInelastic    A2      3.15038             ionIoni             A2     0                   3.15038
于谦烫头
关注
  • 7
    点赞
  • 踩
  • 32
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 14
    评论
专栏目录
Geant4 获取粒子信息
whitepu的博客
12-03 4579
Geant4 获取粒子信息 Step文件里 G4StepPoint* preStepPoint = step->GetPreStepPoint();//step的前端点 G4StepPoint* postStepPoint = step->GetPostStepPoint();//step的后端点 G4String particleID = step->GetTrack()->GetParticleDefinition()->GetParticleName();//获取粒子
GEANT4 step、track、event的信息输出
weixin_42449640的博客
05-04 6239
GEANT4 step、track、event的信息输出GEANT4 step、track、event的信息输出提取函数输出反推粒子的反应过程 GEANT4 step、track、event的信息输出 GEANT4输出step、track、event的信息 提取函数 提取信息的函数 // prepoint、postpoint的volume的name G4String postName = step...
文件的读写基本操作
热门推荐
范高伦的博客
09-07 1万+
一、文件是计算机中数据持久化存储的表现形式 读写文件标准操作格式1： 1、打开文件：file1 = open('文件名','读写模式') 2、操作文件 3、关闭文件：file1.close() 文件操作完毕后必须关闭，否则长期保持对文件的连接状态，造成内存溢出的现象发生 读写文件操作格式2： # 1、打开文件 file1 = open('demo.txt','w') # 2、操作文件 file1.write('hello world') # 3、关闭文件 file1.close() 1、打开文件：wi
geant4安装
08-28
geant4是目前功能相对而言最为齐全的物理仿真平台，这个安装操作用于直到windows下如何配置该开源库。
Geant4 能量沉积
12-04
在B1的基础上更改几何形状，并且利用粒子枪产生各向同性的放射性源；且源在圆柱体内均匀分布，粒子动量方向向各个方向发射概率一定，最终探测在自己设定体积内的能量沉积。
G4数据导出保存简明方法 1
08-03
（有兴趣的同学可以考虑如何智能化地解决这个问题，可以考虑在 RunAction 的 Begin 函数中执行文件清空哦，也就是使用清空原有数据的方式打开文本）基本
Geant 4中Run、Event、Step等数据提取（初学笔记）
budong_2017的博客
10-27 8133
Geant 4运行中涉及到Run、Event、Step、Track等过程，其ID、位置、时间、能量、几何体、材料、粒子类型、反应类型等各类数据均可被提取。数据提取方式较多，如SensitiveDetector，Hit和HitsCollection，以及RunAction， EventAction和SteppingAction等等。 相关知识的介绍可以参见官方手册（http://geant4-us...
Geant4 不完全学习指南９（粒子跟踪和物理过程）
Jack_chf的博客
08-06 7315
　学习总结８简单的学习了物理列表的使用情况，９则是学习其中的一些细节。在实际的模拟过程中，物理过程的选择 和粒子的跟踪有着十分重要的作用，模拟所得到的数据，大多和粒子的跟踪有关。通过前面的学习，我们大概知道了Ｇ４的 一个逻辑，可以简单的描述为Ｒun---Event---Track---Step。 　　当运行一个Ｇ４程序时，也就是开始了一次实验的取数，称为一个Run。Ｒｕn这个名字取自于高能物...
Geant4 不完全学习指南10 （关于track 和step的学习）
Jack_chf的博客
08-06 6748
在不完全学习指南１０中，我简单的学习分析了run,G4Run,event,G4Event,stack（栈）之间的基本知识，下面将继续讨论学习 track 和step的基础知识。 2019年８月６日于兰州　下午４：４０ 　　　　　　　　　　　　　什么是step G4Step 存储一个 step 的瞬态信息。包括 step 的两个端点, PreStepPoint 和 PostStepPoin...
Geant4 储存数据到txt文件中
whitepu的博客
12-03 2142
Geant4 储存数据到txt文件中 因为自己至今仍在挣扎学习如何使用root，所以弱弱的使用txt文本储存数据，然后用python或者matlab处理数据。 即：使用C++ 的方法储存数据。 当par1，par2，par3，par4 已经被赋值（已获取到想要的信息），即可在后面添加： std::fstream file("filename.txt", std::ios::app); file <<par1<<" "<<par2<<" "<&lt
geant4中文相关学习手册
10-27
geant4中文相关学习，适合新手。
Geant4UserGuide.zip
04-10
Geant4的用户手册，从官网爬下来的，保存到本地方便查阅 为Geant4 10.7 2020年12月最新版 该手册是全英文的用户手册，各种操作写的非常清晰和详实，建议入门阅读
Shield.zip_Monte Carlo_geant4_shielding_蒙特卡洛_蒙特卡洛粒子
07-15
通过蒙特卡洛方法计算物体对粒子的屏蔽效果
Geant4 physicsListGuide 中文
12-31
Geant4 physicsListGuide 中文
Geant4学习——入门(基本概念、Geant4工具包的结构、强制类、可选类)
M_M-铭的博客
04-26 1万+
Geant4是由欧洲核子研究组织基于C++面向对象技术开发的蒙特卡罗应用软件包，用于模拟粒子在物质中输运的物理过程。由于具有良好的通用性和扩展能力，Geant4在涉及微观粒子与物质相互作用的诸多领域获得了广泛应用。 第一部分 Geant4的基本概念 HEP中的模拟是“虚拟现实”。仿真既可用于在研发阶段帮助设计探测器，又可用于物理研究中了解探测器的响应。为了创建这样的虚拟现实，我们需要对粒子间的相互作用，几何形状和材料进行建模，以便将基本粒子传播到检测器中。我们还需要描述检测器生成原始数据的灵敏度。 Gea
Geant4学习一：写一个简单程序
ccwwNUAA的博客
10-20 4582
1、 基本介绍
关于geant4和root的笔记
zhy18334723216的博客
01-20 1626
一个geant4项目要分为四个部分，把收集到的信息写进输出文件（root）里面输出出来，主要是通过geant4来进行结构的建模添加反应立场，然后通过root输出自己想要的数据的图形以及参数的具体数值。
GEANT4学习 一些小技巧
weixin_42449640的博客
06-29 1109
GEANT4学习 一些小技巧 记录一些使用中的小技巧 1. 更改随机数种子 在主文件中加入以下语句 G4Random::setTheEngine(new CLHEP::MTwistEngine); G4long seed=time(NULL); CLHEP::HepRandom::setTheSeed(seed); 这样每次运行的时候会根据当前的时间来设置随机数种子。 ...
GEANT4学习 第三节 B2中的SD 敏感探测器 2 hit信息输出
weixin_42449640的博客
06-24 1387
GEANT4学习 第三节 B2中的SD 敏感探测器 2 Hit信息输出 之前的博文看了B2a这个例子的设置的类，接下来看输出信息的类。 B2a的输出信息类只有两个： B2RunAction B2EventAction B2RunAction里面没什么可看的，主要还是看一下B2EventAction 1.B2EventAction 在B2EventAction的EndOfEventAction里面，例子实质上已经写了一个输出hit总数的例子，如下： G4VHitsCollection* hc = even
geant4二次粒子信息获取
最新发布
05-24
在 Geant4 中，可以通过实现 G4SteppingVerbose 类来获取二次粒子的信息。这个类是 Geant4 中用于输出每一步模拟信息的基类。你可以继承 G4SteppingVerbose 类并覆盖其成员函数，以输出你需要的信息。 以下是一个简单的例子，演示了如何获取二次粒子的信息： ```cpp #include "G4SteppingVerbose.hh" #include "G4Track.hh" class MySteppingVerbose: public G4SteppingVerbose { public: void TrackingStarted() override { G4SteppingVerbose::TrackingStarted(); G4Track* track = fpSteppingManager->GetTrack(); G4cout << "Primary particle: " << track->GetDefinition()->GetParticleName() << G4endl; } void StepInfo() override { G4SteppingVerbose::StepInfo(); const G4Step* step = fpSteppingManager->GetStep(); G4Track* track = fpSteppingManager->GetTrack(); if (track->GetCurrentStepNumber() == 1) { G4cout << "Secondary particle: " << step->GetSecondary()->at(0)->GetDefinition()->GetParticleName() << G4endl; } } }; int main() { // ... MySteppingVerbose* mySteppingVerbose = new MySteppingVerbose(); runManager->SetUserVerbose(mySteppingVerbose); // ... return 0; } ``` 在上面的例子中，我们继承了 G4SteppingVerbose 类，并覆盖了其中的 TrackingStarted() 和 StepInfo() 函数。TrackingStarted() 函数在开始跟踪一个粒子时调用，我们在这里输出了主粒子的信息。StepInfo() 函数在每一步结束后调用，我们在这里判断当前步骤是否为第一步，如果是，则输出第一个二次粒子的信息。 注意，为了使用自定义的 G4SteppingVerbose 类，我们需要在创建 G4RunManager 对象之后，调用 SetUserVerbose() 函数，将其设置为用户定义的详细模式。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 14
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

于谦烫头

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付：¥1
获取中
扫码支付

您的余额不足，请更换扫码支付或充值

打赏作者

实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值