# Geometry Factors

Since the intrinsic detection efficiency of the alpha-particle detector in a typical alpha-spectrometry
chamber is nearly 100 %, the overall efficiency of the counting system for a thin alpha-emitting source
is essentially the same as the geometry factor,
defined as the mean solid angle subtended by the detector window at the source, divided by `4π`.
So, if you know the dimensions of the detector and the source, and their relative positions and orientations,
in principle you can predict the efficiency without a calibration.
Or, if you measure the efficiency for a source
in one geometry, you can apply a correction factor to estimate the efficiency for a different geometry.

In the equations below we’ll use spherical coordinates (` ρ, θ, ϕ`), cylindrical coordinates (

`), and Cartesian coordinates (`

*r*,*θ*,*z*`), which are related to each other as follows:`

*x*,*y*,*z*Visualize the ` xy`-plane (

`) as horizontal, with the positive`

*ϕ*= π / 2`-axis, or polar axis (`

*z*`), rising upward.`

*ϕ*= 0## Solid Angle Subtended by a Disk at a Point on its Axis

The solid angle ` Ω` subtended by a surface

`at a point not on the surface is generally given as an integral, which may be a surface integral or a line integral (around the boundary). In spherical coordinates the solid angle subtended by`

*D*`at the origin is given by the surface integral:`

*D*The active window of an alpha-spec detector is typically a flat disk.
Let ` R_{D}` denote the radius of the disk and assume the disk is centered on
the polar axis, parallel to the

`-plane at a fixed height`

*xy*`above it, as shown below`

*h*In this case the surface integral yields the following equation for the solid angle subtended by the disk at the origin:

where ` Φ` is the polar angle

`at the rim of the disk. The cosine of this angle is given by`

*ϕ*So, we get:

Although this expression is theoretically exact, in actual calculations it can produce large round-off
errors if ` h ≫ R_{D}`.
In these situations it is better to apply a few trig identities and compute

`using an equation that is less susceptible to the effects of rounding. For example,`

*Ω*where `tan Φ = R_{D} / h`
and we use the fact that

`0 ≤`. The following version should suffice as long as either

*Φ*< π/2`or`

*h*> 0`and`

*h*= 0`(and neither`

*R*_{D}> 0`nor`

*h*`is huge).`

*R*_{D}## Solid Angle Subtended by an Ellipse

Suppose now that ` D` has an elliptical shape wtih semi-major axis

`, semi-minor axis`

*a*`, eccentricity`

*b*`, and maximum polar angle`

*e*`, but it is still centered on the polar axis at a height`

*Φ*`above the horizontal plane. (Note that`

*h*`tan`.) Now the solid angle subtended by

*Φ*=*a*/*h*`at the origin can be shown to be`

*D*where `Π( n, k)` denotes a complete Legendre elliptic integral of
the third kind, which is defined by the following equation.

**Note**: There are several variants of the definition of the function `Π`. I favor the one used by
Wolfram,
which uses ` k` instead of

`sin`but which in other respects is compatible with the one used by Abramowitz and Stegun. The definition that appears in Numerical Recipes (see below), uses a different sign convention for

*α*`.`

*n*If you can calculate this elliptic integral accurately, you can calculate ` Ω` accurately, at least when the value of

`is not too small. If it is too small, the fact that it is obtained as the difference between`

*Ω*`2π`and another number that is almost equal to

`2π`implies that round-off error can be relatively large.

The value of ` Ω` will be small if either

`(`

*Φ*≈ 0`) or`

*a*/*h*≈ 0`(`

*e*≈ 1`). When`

*b*/*a*≈ 0`, the solid angle is approximated well by the quotient of the area of`

*Φ*≈ 0`and`

*D*`.`

*h*^{2}**Note**: I’m sure there are better approximations.

Suppose instead that ` e ≈ 1`.
As

`from below, the value of the elliptic integral increases without bound. A useful approximation in this case is:`

*e*→ 1where ` E(k, ϕ)` denotes an elliptic integral of the second kind.
If

`is too small, the preceding approximation should be easier to calculate.`

*Φ*## Solid Angle Subtended by a Disk at a Point Not on Its Axis

Next suppose ` D` is again a disk but now consider the solid angle subtended by

`at a point`

*D*`in the`

**P**`-plane at a distance`

*xy*`from the origin.`

*r*The solid angle subtended by the disk ` D` at

`is the same as the solid angle subtended by the image of the disk obtained by a 3-D perspective transformation, where the view plane is oriented so that the image is an ellipse centered on the orthogonal line that passes through the center of projection`

**P**`. The trick then is to find the parameters of that ellipse The following equations provide the necessary values; however, rounding error can be an issue in some circumstances.`

**P**If we define

then we get a substantial simplification, as shown below.

When you use this equation for ` Ω`, there is little reason to worry about round-off error until you
calculate the elliptic integral and perform the final subtraction from

`2π`.

I use Carlson’s method for evaluating the elliptic integral, as presented in Numerical Recipes
[Press et al., 1992, 2007].
For extreme values of the arguments that make ` Ω` very small, you can expect large relative errors
due to rounding in the final subtraction.
The rounding error can be reduced in some cases by the use of a series for the elliptic integral. For example,

where ` n = e^{2} = r^{2} / L^{2}`
and

`and where`

*k*=*r**R*_{D}/*L*^{2}`denotes the Gauss hypergeometric function:`

_{2}*F*_{1}In the equation above for ` Ω`, the value of

`Π(`, multiplied by

*n*,*k*)`4`, is subtracted from

*h*/*L*`2π`. The first term of the series, where

`, can be subtracted from`

*j*= 0`2π`with good accuracy, leaving the higher-order terms to be subtracted. If

`is not large (because`

*e*`is not too far from the axis), those remaining terms are relatively small and do not cause large rounding errors. (If`

**P**`, all the higher-order terms are zero, because`

*e*= 0`.)`

*n*=*k*= 0## The Geometry Factor for an Extended Source

Using this approach you obtain a function ` Ω(r)` for the solid angle defined in terms of the
distance

`of`

*r*`from the origin, assuming fixed values for`

**P**`and`

*R*_{D}`. The geometry factor for an extended source is defined to be the mean solid angle averaged over all points of the source, which is calculated by another integral.`

*h*Given the function ` Ω(r)`, the required integral for a disk-shaped
source of nonzero radius

`(see above) can be calculated using a good technique for numerical integration, such as Gaussian quadrature or Simpson’s Rule.`

*R*_{S}If ` h` is not too small, Gaussian quadrature is likely to give excellent results; but be careful with this technique when

`is very small and`

*h*`or`

*R*_{S}>*R*_{D}`, because the solid angle`

*R*_{S}≈*R*_{D}`is almost a step function in the vicinity of`

*Ω*(*r*)`, dropping steeply from`

*r*=*R*_{D}`to`

*Ω*≈ 2π`as`

*Ω*≈ 0`increases. In this case you can use Simpson’s Rule, which is robust and amenable to brute-force processing; or if`

*r*`is small enough, you can use the fact that:`

*h*

*R*_{S}>*R*_{D}For ` R_{S} < R_{D}`, the limit is just

`2π`.

**To be completed**: (1) There are other equivalent equations for ` Ω(r)`
that can be derived directly from the surface integral; however, my versions of these, shown below, are less simple
than those shown above.

where ` K(k)` denotes the complete elliptic integral of the first kind and where

For ` r ≈ R_{D}` but

`, these alternative equations are probably not great for actual calculations. And when the solid angle is small, you still need to do some work to avoid large rounding errors. (2) It is easier to calculate the solid angle subtended by a polygon than the solid angle subtended by a disk. Approximating the disk by a regular polygon with the same area is a good practical option in the lab. (3) Monte Carlo simulation is another popular approach, although it is inherently inexact. On the other hand it can provide explicit uncertainty estimates. (4) After my coworker described his use of my calculations at a radiochemistry conference in October 2014, I knew I needed to complete my work on evaluating the combined standard uncertainty of the geometry factor. I have done that now but I need to write it up clearly. (5) James Clerk Maxwell discussed the solid angle subtended by an ellipse in A Treatise on Electricity and Magnetism, Chapter XIV. He mentioned a solution using the elliptic integral of the third kind but preferred a different solution in terms of an infinite series of spherical harmonics.`

*r*≠*R*_{D}**Updates**: In 2017 I discovered the work of John T. Conway, who published a nice
paper on this subject in 2006. He provided a closed-form expression for the average geometry factor in terms
of elliptic integrals and also an integral that seems to give highly accurate results without much work.
I have not explored all the cited references, but I recommend the paper anyway:
“Generalizations of Ruby’s formula for the geometric efficiency of a parallel-disk source and detector system,”
Nuclear Instruments & Methods in Physics Research. A 562 (2006) 146–153.