Quasi-real-time range monitoring by in-beam PET: a case for 15O

A fast and reliable range monitoring method is required to take full advantage of the high linear energy transfer provided by therapeutic ion beams like carbon and oxygen while minimizing damage to healthy tissue due to range uncertainties. Quasi-real-time range monitoring using in-beam positron emission tomography (PET) with therapeutic beams of positron-emitters of carbon and oxygen is a promising approach. The number of implanted ions and the time required for an unambiguous range verification are decisive factors for choosing a candidate isotope. An experimental study was performed at the FRS fragment-separator of GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany, to investigate the evolution of positron annihilation activity profiles during the implantation of [Formula: see text]O and [Formula: see text]O ion beams in a PMMA phantom. The positron activity profile was imaged by a dual-panel version of a Siemens Biograph mCT PET scanner. Results from a similar experiment using ion beams of carbon positron-emitters [Formula: see text]C and [Formula: see text]C performed at the same experimental setup were used for comparison. Owing to their shorter half-lives, the number of implanted ions required for a precise positron annihilation activity peak determination is lower for [Formula: see text]C compared to [Formula: see text]C and likewise for [Formula: see text]O compared to [Formula: see text]O, but their lower production cross-sections make it difficult to produce them at therapeutically relevant intensities. With a similar production cross-section and a 10 times shorter half-life than [Formula: see text]C, [Formula: see text]O provides a faster conclusive positron annihilation activity peak position determination for a lower number of implanted ions compared to [Formula: see text]C. A figure of merit formulation was developed for the quantitative comparison of therapy-relevant positron-emitting beams in the context of quasi-real-time beam monitoring. In conclusion, this study demonstrates that among the positron emitters of carbon and oxygen, [Formula: see text]O is the most feasible candidate for quasi-real-time range monitoring by in-beam PET that can be produced at therapeutically relevant intensities. Additionally, this study demonstrated that the in-flight production and separation method can produce beams of therapeutic quality, in terms of purity, energy, and energy spread.

Production and separation of positron emitters for hadron therapy at FRS-Cave M

The FRagment Separator FRS at GSI is a versatile spectrometer and separator for experiments with relativistic in-flight separated short-lived exotic beams. One branch of the FRS is connected to the target hall where the bio-medical cave (Cave M) is located. Recently a joint activity between the experimental groups of the FRS and the biophysics at the GSI and Department of physics at LMU was started to perform biomedical experiments relevant for hadron therapy with positron emitting carbon and oxygen beams. This paper presents the new ion-optical mode and commissioning results of the FRS-Cave M branch where positron emitting ¹⁵O-ions were provided to the medical cave for the first time. An overall conversion efficiency of 2.9±0.2×10^{-4 15}O fragments per primary ¹⁶O ion accelerated in the synchrotron SIS18 was reached.

Precision of the PET activity range during irradiation with10C,11C, and12C beams

Objective. Beams of stable ions have been a well-established tool for radiotherapy for many decades. In the case of ion beam therapy with stable¹²C ions, the positron emitters^10,11C are produced via projectile and target fragmentation, and their decays enable visualization of the beam via positron emission tomography (PET). However, the PET activity peak matches the Bragg peak only roughly and PET counting statistics is low. These issues can be mitigated by using a short-lived positron emitter as a therapeutic beam.Approach.An experiment studying the precision of the measurement of ranges of positron-emitting carbon isotopes by means of PET has been performed at the FRS fragment-separator facility of GSI Helmholtzzentrum für Schwerionenforschung GmbH, Germany. The PET scanner used in the experiment is a dual-panel version of a Siemens Biograph mCT PET scanner.Main results.High-quality in-beam PET images and activity distributions have been measured from the in-flight produced positron emitting isotopes¹¹C and¹⁰C implanted into homogeneous PMMA phantoms. Taking advantage of the high statistics obtained in this experiment, we investigated the time evolution of the uncertainty of the range determined by means of PET during the course of irradiation, and show that the uncertainty improves with the inverse square root of the number of PET counts. The uncertainty is thus fully determined by the PET counting statistics. During the delivery of 1.6 × 10⁷ions in 4 spills for a total duration of 19.2 s, the PET activity range uncertainty for¹⁰C,¹¹C and¹²C is 0.04 mm, 0.7 mm and 1.3 mm, respectively. The gain in precision related to the PET counting statistics is thus much larger when going from¹¹C to¹⁰C than when going from¹²C to¹¹C. The much better precision for¹⁰C is due to its much shorter half-life, which, contrary to the case of¹¹C, also enables to include the in-spill data in the image formation.Significance. Our results can be used to estimate the contribution from PET counting statistics to the precision of range determination in a particular carbon therapy situation, taking into account the irradiation scenario, the required dose and the PET scanner characteristics.

Comparative accuracy and resolution assessment of two prototype proton computed tomography scanners

BACKGROUND

Improving the accuracy of relative stopping power (RSP) in proton therapy may allow reducing range margins. Proton computed tomography (pCT) has been shown to provide state-of-the-art RSP accuracy estimation, and various scanner prototypes have recently been built. The different approaches used in scanner design are expected to impact spatial resolution and RSP accuracy.

PURPOSE

The goal of this study was to perform the first direct comparison, in terms of spatial resolution and RSP accuracy, of two pCT prototype scanners installed at the same facility and by using the same image reconstruction algorithm.

METHODS

A phantom containing cylindrical inserts of known RSP was scanned at the phase-II pCT prototype of the U.S. pCT collaboration and at the commercially oriented ProtonVDA scanner. Following distance-driven binning filtered backprojection reconstruction, the radial edge spread function of high-density inserts was used to estimate the spatial resolution. RSP accuracy was evaluated by the mean absolute percent error (MAPE) over the inserts. No direct imaging dose estimation was possible, which prevented a comparison of the two scanners in terms of RSP noise.

RESULTS

In terms of RSP accuracy, both scanners achieved the same MAPE of 0.72% when excluding the porous sinus insert from the evaluation. The ProtonVDA scanner reached a better overall MAPE when all inserts and the body of the phantom were accounted for (0.81%), compared to the phase-II scanner (1.14%). The spatial resolution with the phase-II scanner was found to be 0.61 lp/mm, while for the ProtonVDA scanner somewhat lower at 0.46 lp/mm.

CONCLUSIONS

The comparison between two prototype pCT scanners operated in the same clinical facility showed that they both fulfill the requirement of an RSP accuracy of about 1%. Their spatial resolution performance reflects the different design choices of either a scanner with full tracking capabilities (phase-II) or of a more compact tracker system which only provides the positions of protons but not their directions (ProtonVDA). This article is protected by copyright. All rights reserved.

Evaluation of the impact of a scanner prototype on proton CT and helium CT image quality and dose efficiency with Monte Carlo simulation

Objective.The use of ion computed tomography (CT) promises to yield improved relative stopping power (RSP) estimation as input to particle therapy treatment planning. Recently, proton CT (pCT) has been shown to yield RSP accuracy on par with state-of-the-art x-ray dual energy CT. There are however concerns that the lower spatial resolution of pCT compared to x-ray CT may limit its potential, which has spurred interest in the use of helium ion CT (HeCT). The goal of this study was to investigate image quality of pCT and HeCT in terms of noise, spatial resolution, RSP accuracy and imaging dose using a detailed Monte Carlo (MC) model of an existing ion CT prototype.Approach.Three phantoms were used in simulated pCT and HeCT scans allowing estimation of noise, spatial resolution and the scoring of dose. An additional phantom was used to evaluate RSP accuracy. The imaging dose required to achieve the same image noise in a water and a head phantom was estimated at both native spatial resolution, and in a scenario where the HeCT spatial resolution was reduced and matched to that of pCT using Hann windowing of the reconstruction filter. A variance reconstruction formalism was adapted to account for Hann windowing.Main results.We confirmed that the scanner prototype would produce higher spatial resolution for HeCT than pCT by a factor 1.8 (0.86 lp mm^-1versus 0.48 lp mm^-1at the center of a 20 cm water phantom). At native resolution, HeCT required a factor 2.9 more dose than pCT to achieve the same noise, while at matched resolution, HeCT required only 38% of the pCT dose. Finally, RSP mean absolute percent error (MAPE) was found to be 0.59% for pCT and 0.67% for HeCT.Significance.This work compared the imaging performance of pCT and HeCT when using an existing scanner prototype, with the spatial resolution advantage of HeCT coming at the cost of increased dose. When matching spatial resolution via Hann windowing, HeCT had a substantial dose advantage. Both modalities provided state-of-the-art RSP MAPE. HeCT might therefore help reduce the dose exposure of patients with comparable image noise to pCT, enhanced spatial resolution and acceptable RSP accuracy at the same time.

Experimental demonstration of accurate Bragg peak localization with ionoacoustic tandem phase detection (iTPD)

Accurate knowledge of the exact stopping location of ions inside the patient would allow full exploitation of their ballistic properties for patient treatment. The localized energy deposition of a pulsed particle beam induces a rapid temperature increase of the irradiated volume and leads to the emission of ionoacoustic (IA) waves. Detecting the time-of-flight (ToF) of the IA wave allows inferring information on the Bragg peak location and can henceforth be used forin-vivorange verification. A challenge for IA is the poor signal-to-noise ratio at clinically relevant doses and viable machines. We present a frequency-based measurement technique, labeled as ionoacoustic tandem phase detection (iTPD) utilizing lock-in amplifiers. The phase shift of the IA signal to a reference signal is measured to derive theToF. Experimental IA measurements with a 3.5 MHz lead zirconate titanate (PZT) transducer and lock-in amplifiers were performed in water using 22 MeV proton bursts. A digital iTPD was performedin-silicoat clinical dose levels on experimental data obtained from a clinical facility and secondly, on simulations emulating a heterogeneous geometry. For the experimental setup using 22 MeV protons, a localization accuracy and precision obtained through iTPD deviates from a time-based reference analysis by less than 15μm. Several methodological aspects were investigated experimentally in systematic manner. Lastly, iTPD was evaluatedin-silicofor clinical beam energies indicating that iTPD is in reach of sub-mm accuracy for fractionated doses < 5 Gy. iTPD can be used to accurately measure theToFof IA signals online via its phase shift in frequency domain. An application of iTPD to the clinical scenario using a single pulsed beam is feasible but requires further development to reach <1 Gy detection capabilities.

Sub-millimeter precise photon interaction position determination in large monolithic scintillators via convolutional neural network algorithms

In this work, we present the development and application of a convolutional neural network (CNN)-based algorithm to precisely determine the interaction position ofγ-quanta in large monolithic scintillators. Those are used as an absorber component of a Compton camera (CC) system under development for ion beam range verification via prompt-gamma imaging. We examined two scintillation crystals: LaBr₃:Ce and CeBr₃. Each crystal had dimensions of 50.8 mm × 50.8 mm × 30 mm and was coupled to a 64-fold segmented multi-anode photomultiplier tube (PMT) with an 8 × 8 pixel arrangement. We determined the spatial resolution for three photon energies of 662, 1.17 and 1.33 MeV obtained from 2D detector scans with tightly collimated¹³⁷Cs and⁶⁰Co photon sources. With the new algorithm we achieved a spatial resolution for the CeBr3 crystal below 1.11(8) mm and below 0.98(7) mm for the LaBr3:Ce detector for all investigated energies between 662 keV and 1.33 MeV. We thereby improved the performance by more than a factor of 2.5 compared to the previously used categorical average pattern algorithm, which is a variation of the well-established k-nearest neighbor algorithm. The trained CNN has a low memory footprint and enables the reconstruction of up to 10⁴events per second with only one GPU. Those improvements are crucial on the way to future clinicalin vivoapplicability of the CC for ion beam range verification.

Fluence-modulated proton CT optimized with patient-specific dose and variance objectives for proton dose calculation

Particle therapy treatment planning requires accurate volumetric maps of the relative stopping power, which can directly be acquired using proton computed tomography (pCT). With fluence-modulated pCT (FMpCT) imaging fluence is concentrated in a region-of-interest (ROI), which can be the vicinity of the treatment beam path, and imaging dose is reduced elsewhere. In this work we present a novel optimization algorithm for FMpCT which, for the first time, calculates modulated imaging fluences for joint imaging dose and image variance objectives. Thereby, image quality is maintained in the ROI to ensure accurate calculations of the treatment dose, and imaging dose is minimized outside the ROI with stronger minimization penalties given to imaging organs-at-risk. The optimization requires an initial scan at uniform fluence or a previous x-ray CT scan. We simulated and optimized FMpCT images for three pediatric patients with tumors in the head region. We verified that the target image variance inside the ROI was achieved and demonstrated imaging dose reductions outside of the ROI of 74% on average, reducing the imaging dose from 1.2 to 0.3 mGy. Such dose savings are expected to be relevant compared to the therapeutic dose outside of the treatment field. Treatment doses were re-calculated on the FMpCT images and compared to treatment doses re-recalculated on uniform fluence pCT scans using a 1% criterion. Passing rates were above 98.3% for all patients. Passing rates comparing FMpCT treatment doses to the ground truth treatment dose were above 88.5% for all patients. Evaluation of the proton range with a 1 mm criterion resulted in passing rates above 97.5% (FMpCT/pCT) and 95.3% (FMpCT/ground truth). Jointly optimized fluence-modulated pCT images can be used for proton dose calculation maintaining the full dosimetric accuracy of pCT but reducing the required imaging dose considerably by three quarters. This may allow for daily imaging during particle therapy ensuring a safe and accurate delivery of the therapeutic dose and avoiding excess dose from imaging.

Integration of spatial distortion effects in a 4D computational phantom for simulation studies in extra-cranial MRI-guided radiation therapy: Initial results

PURPOSE

Spatial distortions in magnetic resonance imaging (MRI) are mainly caused by inhomogeneities of the static magnetic field, nonlinearities in the applied gradients, and tissue-specific magnetic susceptibility variations. These factors may significantly alter the geometrical accuracy of the reconstructed MR image, thus questioning the reliability of MRI for guidance in image-guided radiation therapy. In this work, we quantified MRI spatial distortions and created a quantitative model where different sources of distortions can be separated. The generated model was then integrated into a four-dimensional (4D) computational phantom for simulation studies in MRI-guided radiation therapy at extra-cranial sites.

METHODS

A geometrical spatial distortion phantom was designed in four modules embedding laser-cut PMMA grids, providing 3520 landmarks in a field of view of (345 × 260 × 480) mm³ . The construction accuracy of the phantom was verified experimentally. Two fast MRI sequences for extra-cranial imaging at 1.5 T were investigated, considering axial slices acquired with online distortion correction, in order to mimic practical use in MRI-guided radiotherapy. Distortions were separated into their sources by acquisition of images with gradient polarity reversal and dedicated susceptibility calculations. Such a separation yielded a quantitative spatial distortion model to be used for MR imaging simulations. Finally, the obtained spatial distortion model was embedded into an anthropomorphic 4D computational phantom, providing registered virtual CT/MR images where spatial distortions in MRI acquisition can be simulated.

RESULTS

The manufacturing accuracy of the geometrical distortion phantom was quantified to be within 0.2 mm in the grid planes and 0.5 mm in depth, including thickness variations and bending effects of individual grids. Residual spatial distortions after MRI distortion correction were strongly influenced by the applied correction mode, with larger effects in the trans-axial direction. In the axial plane, gradient nonlinearities caused the main distortions, with values up to 3 mm in a 1.5 T magnet, whereas static field and susceptibility effects were below 1 mm. The integration in the 4D anthropomorphic computational phantom highlighted that deformations can be severe in the region of the thoracic diaphragm, especially when using axial imaging with 2D distortion correction. Adaptation of the phantom based on patient-specific measurements was also verified, aiming at increased realism in the simulation.

CONCLUSIONS

The implemented framework provides an integrated approach for MRI spatial distortion modeling, where different sources of distortion can be quantified in time-dependent geometries. The computational phantom represents a valuable platform to study motion management strategies in extra-cranial MRI-guided radiotherapy, where the effects of spatial distortions can be modeled on synthetic images in a virtual environment.

Laser-driven x-ray and proton micro-source and application to simultaneous single-shot bi-modal radiographic imaging

Radiographic imaging with x-rays and protons is an omnipresent tool in basic research and applications in industry, material science and medical diagnostics. The information contained in both modalities can often be valuable in principle, but difficult to access simultaneously. Laser-driven solid-density plasma-sources deliver both kinds of radiation, but mostly single modalities have been explored for applications. Their potential for bi-modal radiographic imaging has never been fully realized, due to problems in generating appropriate sources and separating image modalities. Here, we report on the generation of proton and x-ray micro-sources in laser-plasma interactions of the focused Texas Petawatt laser with solid-density, micrometer-sized tungsten needles. We apply them for bi-modal radiographic imaging of biological and technological objects in a single laser shot. Thereby, advantages of laser-driven sources could be enriched beyond their small footprint by embracing their additional unique properties, including the spectral bandwidth, small source size and multi-mode emission.

Here the authors show a synchronized single-shot bi-modal x-ray and proton source based on laser-generated plasma. This source can be useful for radiographic and tomographic imaging.

A filtering approach for PET and PG predictions in a proton treatment planning system

Positron emission tomography (PET) and prompt gamma (PG) detection are promising proton therapy monitoring modalities. Fast calculation of the expected distributions is desirable for comparison to measurements and to develop/train algorithms for automatic treatment error detection. A filtering formalism was used for positron-emitter predictions and adapted to allow for its use for the beamline of any proton therapy centre. A novel approach based on a filtering formalism was developed for the prediction of energy-resolved PG distributions for arbitrary tissues. The method estimates PG yields and their energy spectra in the entire treatment field. Both approaches were implemented in a research version of the RayStation treatment planning system. The method was validated against PET monitoring data and Monte Carlo simulations for four patients treated with scanned proton beams. Longitudinal shifts between profiles from analytical and Monte Carlo calculations were within -1.7 and 0.9 mm, with maximum standard deviation of 0.9 mm and 1.1 mm, for positron-emitters and PG shifts, respectively. Normalized mean absolute errors were within 1.2 and 5.3%. When comparing measured and predicted PET data, the same more complex case yielded an average shift of 3 mm, while all other cases were below absolute average shifts of 1.1 mm. Normalized mean absolute errors were below 7.2% for all cases. A novel solution to predict positron-emitter and PG distributions in a treatment planning system is proposed, enabling calculation times of only a few seconds to minutes for entire patient cases, which is suitable for integration in daily clinical routine.

Offline imaging of positron emitters induced by therapeutic helium, carbon and oxygen ion beams with a full-ring PET/CT scanner: experiments in reference targets

In vivo verification of light ion therapy based on positron-emission tomography (PET) imaging of irradiation induced patient activation relies on activity predictions from Monte-Carlo (MC) or analytical computational engines for comparison to the measurements. In order to achieve the necessary accuracy, experimental data are indispensable for the validation of the calculation models. For this we irradiated thick reference targets with mono-energetic helium, carbon and oxygen ion beams and measured the resulting material activation offline with a commercial full-ring PET/CT scanner located nearby the treatment room. Acquired PET data were analysed over time to separate the activity contribution of different radionuclides. Determined production yields were compared to published findings obtained from in-beam activation measurements with a limited-angle double-head PET camera. In addition, we investigated the time-dependence of the measured radionuclide-specific contributions and of the distal activity range, as well as the lateral spread of the activity signal as a function of beam penetration depth. We present radionuclide-specific depth-resolved activity distributions and production yields for the radionuclides [Formula: see text], [Formula: see text] and [Formula: see text], dominating irradiation-induced patient activation. We observe systematically lower production yields with a ratio between the dual-head and our full-ring PET measurements of, on average, 1.7 and 1.3 for the oxygen and carbon beam irradiations, and 1.7 (2.1) for the high (low) energy helium beam irradiations. Findings on the temporal development of the activity range confirm the expectation, with the oxygen beam induced signal being the most sensitive scenario. The experimental data reported in this work, acquired with a state-of-the-art full ring PET scanner, provide a comprehensive and consistent basis for the benchmarking of PET signal calculation engines. In particular, they can support a fine-tuning of the underlying physics models used by the respective implementation and therefore improve the accuracy of PET-based therapy verifications at current and future treatment facilities.

Prediction of positron emitter distributions for range monitoring in carbon ion therapy: an analytical approach

Range verification is one of the most relevant tasks in ion beam therapy. In the case of carbon ion therapy, positron emission tomography (PET) is the most widely used method for this purpose, which images the [Formula: see text]-activation following nuclear interactions of the ions with the tissue nuclei. Since the positron emitter activity profile is not directly proportional to the dose distribution, until today only its comparison to a prediction of the PET profile allows for treatment verification. Usually, this prediction is obtained from time-consuming Monte Carlo simulations of high computational effort, which impacts the clinical workflow. To solve this issue in proton therapy, a convolution approach was suggested to predict positron emitter activity profiles from depth dose distributions analytically. In this work, we introduce an approach to predict positron emitter distributions from depth dose profiles in carbon ion therapy. While the distal fall-off position of the positron emitter profile is predicted from a convolution approach similar to the one suggested for protons, additional analytical functions are introduced to describe the characteristics of the positron emitter distribution in tissue. The feasibility of this approach is demonstrated with monoenergetic depth dose profiles and spread out Bragg peaks in homogeneous and heterogeneous phantoms. In all cases, the positron emitter profile is predicted with high precision and the distal fall-off position is reproduced with millimeter accuracy.

Experimental validation of the FLUKA Monte Carlo code for dose and [Formula: see text]-emitter predictions of radioactive ion beams

In the context of hadrontherapy, whilst ions are capable of effectively destroying radio resistant, deep seated tumors, their treatment localization must be well assessed to ensure the sparing of surrounding healthy tissue and treatment effectiveness. Thus, range verification techniques, such as online positron-emission-tomography (PET) imaging, hold great potential in clinical practice, providing information on the in vivo beam range and consequent tumor targeting. Furthermore, [Formula: see text] emitting radioactive ions can be an asset in online PET imaging, depending on their half-life, compared to their stable counterparts. It is expected that using these radioactive ions the signal obtained by a PET apparatus during beam delivery will be greatly increased, and exhibit a better correlation to the Bragg Peak. To this end, FLUKA Monte Carlo particle transport and interaction code was used to evaluate, in terms of annihilation events at rest and dose, the figure of merit in using [Formula: see text] emitter, radioactive ion beams (RI [Formula: see text]). For this purpose, the simulation results were compared with experimental data obtained with an openPET prototype in various online PET acquisitions at the Heavy Ion Medical Accelerator in Chiba (HIMAC), in collaboration with colleagues from the National Institute of Radiological Sciences' (NIRS) Imaging Physics Team. The dosimetry performance evaluation with FLUKA benefits from its recent developments in fragmentation production models. The present work estimated that irradiations with RI [Formula: see text], produced via projectile fragmentation and their signal acquisition with state-of-the-art PET scanner, lead to nearly a factor of two more accurate definition of the signals' peak position. In addition to its more advantageous distribution shape, it was observed at least an order magnitude higher signal acquired from ¹¹C and ¹⁵O irradiations, with respect to their stable counterparts.

An overview of recent developments in FLUKA PET tools

The new developments of the FLUKA Positron-Emission-Tomography (PET) tools are detailed. FLUKA is a fully integrated Monte Carlo (MC) particle transport code, used for an extended range of applications, including Medical Physics. Recently, it provided the medical community with dedicated simulation tools for clinical applications, including the PET simulation package. PET is a well-established imaging technique in nuclear medicine, and a promising method for clinical in vivo treatment verification in hadrontherapy. The application of clinically established PET scanners to new irradiation environments such as hadrontherapy requires further experimental and theoretical research to which MC simulations could be applied. The FLUKA PET tools, besides featuring PET scanner models in its library, allow the configuration of new PET prototypes via the FLUKA Graphical User Interface (GUI) Flair. Both the beam time structure and scan time can be specified by the user, reproducing PET acquisitions in time, in a particle therapy scenario. Furthermore, different scoring routines allow the analysis of single and coincident events, and identification of parent isotopes generating annihilation events. Two reconstruction codes are currently supported: the Filtered Back-Projection (FBP) and Maximum-Likelihood Expectation Maximization (MLEM), the latter embedded in the tools. Compatibility with other reconstruction frameworks is also possible. The FLUKA PET tools package has been successfully tested for different detectors and scenarios, including conventional functional PET applications and in beam PET, either using radioactive sources, or simulating hadron beam irradiations. The results obtained so far confirm the FLUKA PET tools suitability to perform PET simulations in R&D environment.

A novel approach to electron data background treatment in an online wide-angle spectrometer for laser-accelerated ion and electron bunches

Laser-based ion acceleration is driven by electrical fields emerging when target electrons absorb laser energy and consecutively leave the target material. A direct correlation between these electrons and the accelerated ions is thus to be expected and predicted by theoretical models. We report on a modified wide-angle spectrometer, allowing the simultaneous characterization of angularly resolved energy distributions of both ions and electrons. Equipped with online pixel detectors, the RadEye1 detectors, the investigation of this correlation gets attainable on a single shot basis. In addition to first insights, we present a novel approach for reliably extracting the primary electron energy distribution from the interfering secondary radiation background. This proves vitally important for quantitative extraction of average electron energies (temperatures) and emitted total charge.