Skip to content

Photon flux

Three levels of calculation, each building on the previous:

Function What it computes Considers geometry?
calculate_transmission Material attenuation only (exp(-Sigma*t)) No
calculate_flux S / (4piR^2) * exp(-Sigma*t) * B Yes (inverse square law)
calculate_dose Flux * ICRP-116 dose coefficient Yes (inverse square law)

calculate_flux includes source strength (particles/s), inverse-square-law spreading over the total distance, and an optional build-up factor.

Example

import rad_point_kernel as rpk

iron = rpk.Material(composition={"Fe": 1.0}, density=7.874)
layers = [
    rpk.Layer(thickness=1000),               # 10 m void
    rpk.Layer(thickness=10, material=iron),  # 10 cm iron
]

source = rpk.Source(particle="photon", energy=662e3)
result = rpk.calculate_flux(source_strength=1e12, layers=layers, source=source)
print(f"Flux: {result.uncollided_flux} photons/cm2/s")
print(f"Transmission: {result.transmission_fraction}")
print(f"Optical thickness: {result.optical_thickness}")
print(f"Build-up factor: {result.buildup_factor}")
print(f"Distance: {result.total_distance_cm} cm")

Flux: 252.00537358450575 photons/cm2/s
Transmission: 0.0032304454589905237
Optical thickness: 5.735135238256533
Build-up factor: 1.0
Distance: 1010.0 cm

The returned CalcResult object has these properties:

  • uncollided_flux - flux at the outer surface (particles/cm2/s); equals geometry * transmission * B, so it includes the build-up correction when one is supplied
  • transmission_fraction - exp(-Sigma*t)
  • optical_thickness - Sum(Sigma_r,i * t_i)
  • buildup_factor - B (1.0 if no build-up model given)
  • total_distance_cm - total distance from source to detector

With a manual build-up factor

If you already have a build-up factor from tabulated data or a prior run, pass it via BuildupModel.constant. See Calculate build-up with MC for how to compute B yourself.

import rad_point_kernel as rpk

iron = rpk.Material(composition={"Fe": 1.0}, density=7.874)
layers = [
    rpk.Layer(thickness=1000),
    rpk.Layer(thickness=10, material=iron),
]

source = rpk.Source(particle="photon", energy=662e3)

B_flux = 2.0
result = rpk.calculate_flux(
    source_strength=1e12,
    layers=layers,
    source=source,
    buildup=rpk.BuildupModel.constant(B_flux),
)
print(f"Flux (B={B_flux}): {result.uncollided_flux} photons/cm2/s")
print(f"Applied build-up:  {result.buildup_factor}")

Flux (B=2.0): 504.0107471690115 photons/cm2/s
Applied build-up:  2.0