猪头Geant4讲座猪头Geant4讲座.pdfVIP

猪头Geant4 讲座第一讲——预备知识 先来两句套话，本人知识有限，不当之处请大家多多指教:P~~~ 预备知识，不爱聊天兄已经说过了，C++ 这里我重点说一下，重点要学习C++中的类的使用，包括类的调用、类的初始化、类的重构等 其次，基本物理概念和基本Linux 基础也是必要的 当然，Geant4 是要安装好的，G4 手册是要有的。从哪儿找，怎么装，这些我都不再具体叙述了 大家自己翻前面的帖子吧。 猪头Geant4 讲座第二讲——模拟算法 我们要学习Geant4 首先应该学习G4 是如何处理模拟过程的。 在G4 中一个典型的模拟算法是这样的。 首先建立一次模拟，在G4 中称为一次Run Run 建立后，需要对几何结构、物理过程进行初始化 初始化完成后就开始模拟过程了， 首先发射一个粒子，每一步都按照蒙卡方法进行模拟，具体模拟方法请参阅裴鹿成或许淑艳老师的书 这里不具体讲，因为不是重点 在G4 中，发射一个（或一系列）粒子到所有次级粒子死亡的过程成为一次Event。 而每次发射的初始粒子则有粒子发射器进行控制。 而在每一个 event 过程中，粒子与材料反应后会可能生成多个次级粒子，每个粒子都会有一条径迹，称之 为track 而每一个粒子（初始的或次级的）的径迹又是由很多步组成的，称之为step 关于track 和step 的理解请参阅下图 最后总结一下，G4 模拟的基本算法是 A Run Start - 初始化物理模型/几何模型 - An Event Start - 调用粒子发射器发射粒子 - A Track Start - A Step Start - A Step End - Next Step Start - …… - All Step End - A Track End - Next Track Start - …… - All Track End - An Event End - Next Event Strat - …… - All Event End(All Primaries Shot) - A Run End - Next Run Start - …… 猪头Geant4 讲座第三讲——程序结构及基本类 G4 是采用的gcc 编译器，因此其程序结构是和C++一样的。 首先包括有一个主程序main，然后分别包含有子程序src 和头文件include 以及其他调用文件others G4 里面为了与C++相区别，程序后缀都是.cc，头文件后缀都是.hh 其中头文件在.cc 里面写也没问题，但是那样看起来不方便，建议还是按照c++的习惯一一对应比较好。 那么，关键问题是要进行一个模拟我们都需要写哪些src 和include 的？ 下面我们首先看G4 里面的几个基本类，这些基本类基本上是与src 一一对应的。 G4RunManager——对应主程序 这个类在主程序中用以初始化模拟信息的，或者形象地说是用于连接子程序的， 而连接方式是通过Set 函数来完成的 大家可以从$G4INSTALL/source/run/G4RunManager.hh 里面查看各种Set 函数，如 public: // with description inline void SetUserInitialization(G4VUserDetectorConstruction* userInit) { userDetector = userInit; } inline void SetUserInitialization(G4VUserPhysicsList* userInit) { physicsList = userInit; kernel-SetPhysics(userInit); } 可以说G4RunManager 类是贯穿整个程序模拟过程的总线，因此一般说来只能有一个 而开始一次Run 的信号则是通过BeamOn 函数发出的，其格式是 virtual void BeamOn(G4int n_event,const char* macroFile=0,G4int n_select=-1); 可以通过多次调用BeamOn 来实现循环计算。 其余子类包括几何结构类、物理设定类、粒子发射器类（源描述类）、事件处理类、径迹处理类等。 这些