_Geant4 Version: 4.11
_Operating System: Windows 10
_Compiler/Version: MSVC/14.29.30133
_CMake Version: 3.22.1.0
Hello Geant4 community,
I’m trying to simulate a point source of Co-60 (2 GBq) that emits two γ-rays:
- 1.17323 MeV with intensity 99.85%
- 1.33249 MeV with intensity 99.9826%
My goal is to sample both lines in a single GPS run using a histogram energy distribution, rather than running two separate monoenergetic simulations and summing the results. I’ve tried several approaches with /gps/source/intensity
and /gps/source/add
, but the deposited dose rate at 200 cm (10 cm probe, 10 M primaries) varies wildly:
Case | Source definition | Result (μGy/h) |
---|---|---|
1 | /gps/source/intensity 0.9985 /gps/particle gamma /gps/pos/type Point /gps/pos/centre 0 0 0 cm /gps/ang/type iso /gps/ene/type Mono /gps/ene/mono 1.17323 MeV /gps/source/add 0.999826 /gps/particle gamma /gps/pos/type Point /gps/pos/centre 0 0 0 cm /gps/ang/type iso /gps/ene/type Mono /gps/ene/mono 1.33249 MeV |
72.46 |
2 | /gps/source/add 0.9985 /gps/particle gamma /gps/pos/type Point /gps/pos/centre 0 0 0 cm /gps/ang/type iso /gps/ene/type Mono /gps/ene/mono 1.17323 MeV /gps/source/add 0.999826 /gps/particle gamma /gps/pos/type Point /gps/pos/centre 0 0 0 cm /gps/ang/type iso /gps/ene/type Mono /gps/ene/mono 1.33249 MeV |
37.19 |
3 | /gps/source/intensity 0.9985 /gps/particle gamma /gps/pos/type Point /gps/pos/centre 0 0 0 cm /gps/ang/type iso /gps/ene/type Mono /gps/ene/mono 1.17323 MeV /gps/source/intensity 0.999826 /gps/particle gamma /gps/pos/type Point /gps/pos/centre 0 0 0 cm /gps/ang/type iso /gps/ene/type Mono /gps/ene/mono 1.33249 MeV |
67.54 |
4 | Run A:/gps/source/intensity 0.9985 … /gps/ene/mono 1.17323 MeV Run B: /gps/source/intensity 0.999826 … /gps/ene/mono 1.33249 MeV |
57.05 + 67.54 = 124.59 |
5 | Run A:/gps/source/add 0.9985 … /gps/ene/mono 1.17323 MeV Run B: /gps/source/add 0.999826 … /gps/ene/mono 1.33249 MeV |
24.56 + 26.14 = 50.71 |
None of these matches my expectation of roughly ~ 134 μGy/h for the correctly weighted double line. I suspect I misunderstand:
- Difference between
/gps/source/intensity
and/gps/source/add
- How to normalize branching ratios so total probability = 1
- Proper use of the
/gps/hist/
commands for a two-point energy spectrum
Minimal macro for context
/control/verbose 1
/geometry/source geometry.tg
/physics_lists/select Shielding
/run/initialize
/control/alias SourcePosition 0.0 0.0 0.0
/control/execute Co-60_Point.dat
/score/create/probe ProbeVol 10.0 cm
/score/probe/material G4_AIR
/score/probe/locate 200.0 0.0 0.0 cm
/score/quantity/doseDeposit dDep
/vis/disable
/tracking/verbose 0
/run/beamOn 10000000
/score/dumpQuantityToFile ProbeVol dDep output.csv
/score/close
My questions:
- What is the correct way to define a two-line Co-60 source in one GPS histogram run?
- How should the histogram points be specified (e.g.
/gps/hist/point
) and normalized? - Can someone share a minimal example using
/gps/hist/type energy
(or/gps/hist/file
) to sample exactly these two energies with the correct branching ratios? - Are there pitfalls in mixing
/gps/source/intensity
vs/gps/source/add
that I should watch out for?