Beam Energy Measurement using a Bayesian Approach with the Stacked Foil Method

This paper presents a robust Bayesian method for measuring proton beam energy at a medical cyclotron using the stacked foil technique and $^{48}$V activity, demonstrating its accuracy, uncertainty handling, and suitability for diverse experimental conditions without relying on direct charge measurements.

The Gist

The Big Idea: Measuring Speed Without a Speedometer

Imagine you have a high-speed car (a proton beam) zooming out of a factory (a medical cyclotron). You need to know exactly how fast it is going to ensure it hits its target safely and effectively. Usually, you'd use a speedometer or a radar gun. But in this specific lab, the "road" conditions are tricky. Sometimes the car is driving through a vacuum, but other times it's driving through air, or the "speedometer" (electrical current measurement) gets confused by the environment and gives bad readings.

The authors of this paper developed a clever, low-tech way to figure out the car's speed without needing a working speedometer. They call it the "Stacked Foil Method," and they used a mathematical tool called Bayesian inference (think of it as a super-smart detective's logic) to solve the mystery.

The Detective's Toolkit: The "Stacked Foil" Sandwich

Instead of a radar gun, the team built a sandwich made of very thin metal sheets (foils) of Titanium, Copper, and Niobium.

1. The Setup: They place these metal sheets in the path of the proton beam.
2. The Reaction: As the protons smash through the first sheet, they lose a tiny bit of energy (like a runner getting tired). When they hit the second sheet, they are slightly slower. By the time they reach the last sheet, they are much slower.
3. The Clue: When the protons hit the metal, they turn some of the atoms in the metal into a different, radioactive version (like turning a regular apple into a glowing one). This is called "induced activity."
4. The Measurement: After the beam stops, they take the sheets out and measure how "glowing" (radioactive) each one is using a special camera (a gamma spectrometer).

The Analogy: Imagine throwing a ball at a series of five thin walls.

• If the ball is thrown very hard, it punches through all five walls and leaves a big mark on the last one.
• If the ball is thrown gently, it might only punch through the first two walls and leave a tiny mark on the third, with no marks on the rest.
• By looking at which walls have marks and how big those marks are, you can work backward to figure out exactly how hard the ball was thrown, even if you didn't see the throw itself.

The "Magic" Math: Bayesian Inference

The team didn't just guess the speed. They used a method called Bayesian inference.

• The Old Way (Frequentist): Imagine trying to solve a puzzle where you have to guess the speed, calculate what the marks should look like, and then tweak your guess until it matches. If the puzzle is complex (which it is, because the physics is non-linear), this method often gets stuck or underestimates how unsure you really are.
• The New Way (Bayesian): Imagine a detective who starts with a list of possible speeds (e.g., "It's probably between 8 and 19 MeV"). Then, they look at the actual glowing marks on the metal foils. They use a computer to simulate millions of scenarios, asking: "If the speed was X, would we see these marks?"
• The Result: The computer quickly eliminates the impossible speeds and narrows down the list to one very precise answer. It also naturally accounts for "nuisance factors"—things that might mess up the data, like slight variations in the thickness of the metal sheets or small errors in the known physics of the reactions. It treats the total amount of electricity (current) as a "mystery variable" it solves for, rather than needing to measure it perfectly first.

What They Found

The team tested this method in four different scenarios at the Bern Medical Cyclotron:

1. Pristine Beam: Measuring the beam right as it leaves the machine.
2. After Scatterer: Measuring the beam after it passed through a metal screen and some air (which slows it down).
3. Cell Level: Measuring the beam after it passed through a window, an ionization chamber, and the wall of a cell culture flask. This is a "messy" environment where traditional current measurements fail, but their method worked perfectly.
4. Solid Target Station: Measuring the beam at a different exit port.

The Results:

• They successfully measured beam energies ranging from 8 MeV to 19 MeV.
• The method was accurate even when the beam passed through air or other materials that usually confuse standard sensors.
• They found that they didn't need a huge stack of foils; even a smaller stack could give a reliable answer if the math was done right.
• They also checked if their results depended on which "rulebook" (cross-section data) they used for the physics. They found that even if they used slightly different physics data, their speed estimates didn't change much, proving the method is robust.

Why This Matters (According to the Paper)

The paper highlights that this method is calibration-free and simple.

• No Special Gear: You don't need expensive, complex beamline equipment. You just need metal foils and a standard gamma detector.
• Works in "Dirty" Conditions: It works in low-vacuum or air-exposed setups where traditional electrical current measurements are unreliable (because air can mess up the electrical reading).
• Versatile: It can be used in almost any accelerator lab because it relies on standard tools rather than custom-built sensors.

In short, the authors created a "glow-in-the-dark" speed trap for protons that works even when the usual sensors are confused, using a smart mathematical detective to figure out the speed based on how much the metal sheets "glow."

Technical Summary

Technical Summary: Beam Energy Measurement using a Bayesian Approach with the Stacked Foil Method

Problem Statement
Precise knowledge of proton beam energy is critical for applications involving solid target irradiation, such as the production of medical radioisotopes, cross-section measurements, and radiobiology dosimetry. While established methods like time-of-flight, Rutherford backscattering, and magnetic deflection exist, they often require complex vacuum setups and are unsuitable for low-vacuum or air-exposed environments where direct current or charge measurements are unreliable due to air ionization. At the Bern Medical Cyclotron (BMC), there is a need for a robust, calibration-free method to characterize beam energy across a range of 8–19 MeV, including configurations where the beam passes through degraders and air.

Methodology
The authors propose a method combining the stacked foil technique with Bayesian inference.

• Experimental Setup: A stack of metal foils (Titanium, Copper, and Niobium) is irradiated. Titanium and Copper serve as active monitor foils, while Niobium acts as a beam degrader. The stack is placed in various configurations at the BMC, including the pristine beam line (BTL), after an aluminum scatterer, at the "cell level" (through air and a flask wall), and at a commercial solid target station (STS).
• Physical Model: The method relies on the induced radioactivity of specific monitor reactions: $^{nat}\text{Ti}(p,x)^{48}\text{V}$ and $^{nat}\text{Cu}(p,x)^{65}\text{Zn}$. The expected activity $A_{exp}$ in each foil is calculated analytically based on the integrated current $Q$, the beam energy $E_{in}$, foil thickness, mass stopping power, and reaction cross-sections.
• Bayesian Framework: Instead of standard frequentist fitting, a Bayesian approach is employed to fit the measured activities against the model predictions.
  • Parameters: The initial beam energy $E_0$ is the primary parameter of interest. The integrated current $Q$ is treated as a nuisance parameter (a free variable) rather than a measured input, allowing energy determination without direct current measurement.
  • Uncertainties: The model incorporates uncertainties in cross-sections (using IAEA data with a 10% scaling factor), stopping powers (5% scaling), and foil thicknesses (15% scaling).
  • Validation: Monte Carlo simulations (FLUKA) were used to verify that energy dispersion within the stack has a negligible impact on the yield, justifying the use of a mono-energetic beam model for the analytical calculation.

Key Contributions

1. Current-Independent Measurement: The method successfully infers beam energy without relying on direct current or charge measurements, making it applicable in low-vacuum or air-exposed setups where ionization currents are unreliable.
2. Robust Uncertainty Treatment: The Bayesian framework provides a consistent treatment of nuisance parameters (like $Q$) and systematically handles non-linear dependencies and correlations between parameters (e.g., cross-section scaling and current).
3. Extended Energy Range: By utilizing a combination of Ti and Cu foils, the method covers the 8–19 MeV range. The inclusion of Cu foils allows for reference at lower energies (4–6 MeV) where the Ti cross-section is less sensitive.
4. Cross-Section Sensitivity Analysis: The study evaluates the impact of different cross-section datasets (IAEA recommended vs. custom EXFOR fits) on the final energy determination, demonstrating that the method is robust to these variations within the tested energy ranges.

Results
The method was applied to four distinct experimental configurations at the BMC:

• Pristine BTL Energy ($1^\circ$): Measured at 18.03 ± 0.09 MeV.
• After Scatterer ($2^\circ$): Measured at 15.54 ± 0.12 MeV, consistent with LISE++ simulations of energy degradation.
• Cell Level ($3^\circ$): Measured at 8.14 ± 0.29 MeV, representing the beam after passing through air and a flask wall.
• STS Energy ($4^\circ$): Measured at 17.14 ± 0.13 MeV.

The results show that the measured activities align well with the Bayesian model predictions. The study found that for energies above 13 MeV, a stack of approximately ten 25 µm Titanium foils is sufficient for accurate measurement. For lower energies, the inclusion of Copper foils is necessary. The analysis of reduced datasets suggests that equivalent relative errors can be achieved with fewer foils than the full stack, provided data points are retained in regions of steep cross-section gradients.

Significance and Claims
The paper claims that this approach provides an accurate, practical, and versatile tool for beam energy characterization in medical cyclotron environments. Its primary significance lies in its ability to function in "non-standard" setups (such as air-exposed radiobiology experiments) where traditional current-based methods fail. The authors emphasize that the method requires no special beamline equipment, relying only on standard gamma spectroscopy (HPGe detectors) and simple foil stacks. Furthermore, by treating the integrated current as a free parameter, the method allows for an independent consistency check of current-on-target measurements in configurations where Faraday cups are not feasible. The study concludes that the stacked foil method, enhanced by Bayesian inference, is a reliable alternative for beam characterization in the 8–19 MeV range.