Geant4 example extend radioactivedecay

如何使用 241Am 放射源？

mac 文件

 1 /control/verbose 2
 2 /run/verbose 1
 3 #
 4 /gun/particle ion
 5 /gun/ion 95 241
 6 #
 7 /tracking/verbose 2
 8 /run/beamOn 1
 9 #
10 /control/cout/ignoreThreadsExcept 0
11 /tracking/verbose 0
12 #
13 /analysis/setFileName 241Am
14 /analysis/h1/set 1  15000  0. 15000 keV    #e+ e-
15 /analysis/h1/set 2  15000  0. 15000 keV    #neutrino
16 /analysis/h1/set 3  15000  0. 15000 keV    #gamma
17 /analysis/h1/set 4  15000  0. 15000 keV    #alpha
18 /analysis/h1/set 6  10000  0. 25000 keV    #EkinTot (Q)
19 /analysis/h1/set 7  15000  0. 15e4 keV    #P balance
20 /analysis/h1/set 8  10000  0. 100. year    #time of life
21 /analysis/h1/set 9  10000  1. 30. MeV      #EvisTot
22 #
23 /run/printProgress 100000
24 #  
25 /run/beamOn 1000000

第 13-21 行，使用 root 进行数据分析是固定的写法。
最终得到的 alpha 数据如图，这里要注意，生成的 root 文件需要用 TBrowser 打开，然后在其中修改 histogram 的名字，使用 "HIST" 参数将直方图的误差棒去掉，画成常见的格式。

问题 1，积分该图中的所有事件，是没有 7999996 个事件的，为什么？
问题 2，为什么在我的测试中，即使是使用 gps，发出 241Am 源，灵敏探测器上获取的能量有十几 MeV 那么大？
在 reachergate 上，有一个类似的问题，https://www.researchgate.net/post/How_can_I_simulate_Co60_source_in_geant4 .

回答：

问题2，事实上，241Am 整个衰变链特别长，十几 MeV 的 α 粒子来源于其衰变的子核，以及子核的子核...而实际情况中，241Am 的放射源不存在测量不到其他的 α 能量，是由于 241Am α 衰变后的子核 237Np，其半衰期有 2.144M y 之久，因此无法看到。而在 G4 的模拟中，如果不加限制，是会一直跑下去的，因此会得到整条衰变链的 α 能量。