Geant4程序的一点技巧总结 2

 

 

4. 关于获得粒子状态

获得粒子的总能量：

G4Track *pTrack = pStep->GetTrack(); const G4DynamicParticle *pParticle = pTrack->GetDynamicParticle(); pParticle->GetTotalEnergy(); 

在Step中获得例子名称：

G4Track *pTrack = pStep->GetTrack(); const G4DynamicParticle *pParticle = pTrack->GetDynamicParticle(); G4String particleName = pParticle->GetDefinition()->GetParticleName(); 

.

其他状态信息（用户手册内容，摘一些常用的）

<1>在一个Step中获得指向前一或后一个Step状态的指针

G4StepPoint* point1 = step->GetPreStepPoint();G4StepPoint* point2 = step->GetPostStepPoint(); 

<2>获得粒子坐标位置

G4ThreeVector pos1 = step->GetPosition(); 

<3>获得粒子所在物理体

G4TouchableHandle touch1 = point1->GetTouchableHandle();G4VPhysicalVolume* volume = touch1->GetVolume(); 

 进一步获得物理体名称和拷贝号

G4String name = volume->GetName(); G4int copyNumber = touch1->GetCopyNumber(); 

 进一步获得物理体对应的逻辑题

G4LogicalVolume* lVolume = volume->GetLogicalVolume(); 

 进一步获得物理体的上级物理体

G4VPhysicalVolume* mother = touch1->GetVolume(depth=1); 

 同理 depth=2 ...可以获得上上级物理体...

<4>判断粒子刚刚进入当前物理体

if (point1->GetStepStatus() == fGeomBoundary) 

 判断粒子即将离开当前物理体

if (point2->GetStepStatus() == fGeomBoundary) 

<5>获得当前过程(Step)沉积的能量

G4double eDeposit = step->GetTotalEnergyDeposit(); 

<6>获得粒子当前动能(粒子能量)

G4double kinEnergy = track->GetKineticEnergy(); 

.

5.在G4例子程序里找“包含指定内容的例子程序”的Linux指令

找到包含 proton 内容的文件清单(用于在example里找内容)：

find -name '*.cc' -exec grep -l 'proton' {} /;【grep的常用参数：】【-i 在字符串比较的时候忽略大小写】【-l 只显示包含匹配模板的行的文件名清单】【-n 在每一行前面打印该行在文件中的行数】 

.

6.关于可视化中粒子的颜色

红色-电子(负电)

蓝色-质子或a粒子(正电)

绿色-光子或中子(不带电)

.

7.关于能量截断

G4和MCNP等蒙卡软件的一个不同点是用户只能指定长度截断，由系统自动转换成每种材料的能量截断。

如何从长度截断知道能量截断：

在运行G4程序后会在控制台输出一些运行参数的信息，其中有一段Table ofregistered couples，这里可以看到指定的长度截断在每个材料中的能量截断，例如：

========= Table of registered couples ==============================Index : 0 used in the geometry : Yes recalculation needed : No Material : Air Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 1 keV Region(s) which use this couple : DefaultRegionForTheWorldIndex : 1 used in the geometry : Yes recalculation needed : No Material : Soil Range cuts : gamma 10 um e- 10 um e+ 10 um proton 10 um Energy thresholds : gamma 6.5864 keV e- 618.791 keV e+ 590.357 keV proton 1 keV Region(s) which use this couple : DefaultRegionForTheWorld==================================================================== 

其中每个材料的Energythresholds项就是能量截断值

.

8.关于单位换算

单位换算的时候，将以个数以MeV单位记录，应该除以MeV(从一个单位到另一个单位)，如：

G4double particleEnergy = (pParticles->GetKineticEnergy())/MeV; 

而赋予变量单位时才是乘(从没有单位到有单位)：

G4double theEnergy = 100*MeV; 

.

9.关于源文件后缀的注意

Geant4的源文件必须是.cc而不能是.cpp格式，否则将不被编译。

但头文件的后缀.hh或者.h都可以（有些IDE对.hh的文件不能提供代码自动识别，就可以改成.h来用），头文件本身不是直接参与编译的。