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 的模拟中，如果不加限制，是会一直跑下去的，因此会得到整条衰变链的 α 能量。