Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Detector Response”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Quantitative Monte Carlo-based holmium-166 SPECT reconstruction

Purpose: Quantitative imaging of the radionuclide distribution is of increasing interest for microsphere radioembolization (RE) of liver malignancies, to aid treatment planning and dosimetry. For this purpose, holmium-166 ({sup 166}Ho) microspheres have been developed, which can be visualized with a gamma camera. The objective of this work is to develop and evaluate a new reconstruction method for quantitative {sup 166}Ho SPECT, including Monte Carlo-based modeling of photon contributions from the full energy spectrum.Methods: A fast Monte Carlo (MC) simulator was developed for simulation of {sup 166}Ho projection images and incorporated in a statistical reconstruction algorithm (SPECT-fMC). Photon scatter and attenuation for all photons sampled from the full {sup 166}Ho energy spectrum were modeled during reconstruction by Monte Carlo simulations. The energy- and distance-dependent collimator-detector response was modeled using precalculated convolution kernels. Phantom experiments were performed to quantitatively evaluate image contrast, image noise, count errors, and activity recovery coefficients (ARCs) of SPECT-fMC in comparison with those of an energy window-based method for correction of down-scattered high-energy photons (SPECT-DSW) and a previously presented hybrid method that combines MC simulation of photopeak scatter with energy window-based estimation of down-scattered high-energy contributions (SPECT-ppMC+DSW). Additionally, the impact of SPECT-fMC on whole-body recovered activities (A{sup est}) and estimated radiation absorbed doses was evaluated using clinical SPECT data of six {sup 166}Ho RE patients.Results: At the same noise level, SPECT-fMC images showed substantially higher contrast than SPECT-DSW and SPECT-ppMC+DSW in spheres ≥17 mm in diameter. The count error was reduced from 29% (SPECT-DSW) and 25% (SPECT-ppMC+DSW) to 12% (SPECT-fMC). ARCs in five spherical volumes of 1.96–106.21 ml were improved from 32%–63% (SPECT-DSW) and 50%–80% (SPECT-ppMC+DSW) to 76%–103% (SPECT-fMC). Furthermore, SPECT-fMC recovered whole-body activities were most accurate (A{sup est}= 1.06 × A − 5.90 MBq, R{sup 2}= 0.97) and SPECT-fMC tumor absorbed doses were significantly higher than with SPECT-DSW (p = 0.031) and SPECT-ppMC+DSW (p = 0.031).Conclusions: The quantitative accuracy of {sup 166}Ho SPECT is improved by Monte Carlo-based modeling of the image degrading factors. Consequently, the proposed reconstruction method enables accurate estimation of the radiation absorbed dose in clinical practice.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2014

TH-C-BRD-01: Analytical Computation of Prompt Gamma Ray Emission and Detection for Proton Range Verification

Purpose: A prompt gamma (PG) slit camera prototype demonstrated that on-line range monitoring within 1–2 mm could be performed by comparing expected and measured PG detection profiles. Monte Carlo (MC) can simulate the expected PG profile but this would result in prohibitive computation time for a complete pencil beam treatment plan. We implemented a much faster method that is based on analytical processing of pre-computed MC data. Methods: The formation of the PG detection signal can be separated into: 1) production of PGs and 2) detection by the camera detectors after PG transport in geometry. For proton energies from 40 to 230 MeV, PG productions in depth were pre-computed by MC (PENH) for 12C, 14N, 16O, 31P and 40Ca. The PG production was then modeled analytically by adding the PG production for each element according to local proton energy and tissue composition.PG transport in the patient/camera geometries and the detector response were modeled by convolving the PG production profile with a transfer function. The latter is interpolated from a database of transfer functions fitted to pre-computed MC data (PENELOPE). The database was generated for a photon source in a cylindrical phantom with various radiuses and a camera placed at various positions.As a benchmark, the analytical model was compared to PENH for a water phantom, a phantom with different slabs (adipose, muscle, lung) and a thoracic CT. Results: Good agreement (within 5%) was observed between the analytical model and PENH for the PG production. Similar accuracy for detecting range shifts was also observed. Speed of around 250 ms per profile was achieved (single CPU) using a non-optimized MatLab implementation. Conclusion: We devised a fast analytical model for generating PG detection profiles. In the test cases considered in this study, similar accuracy than MC was achieved for detecting range shifts. This research is supported by IBA.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Sci—Fri PM: Dosimetry—06: Commissioning of a 3D patient specific QA system for hypofractionated prostate treatments

The objective of this work is to commission the 6MV-SRS beam model in COMPASS (v2.1, IBA-Dosimetry) and validate its use for patient specific QA of hypofractionated prostate treatments. The COMPASS system consists of a 2D ion chamber array (MatriXX{sup Evolution}), an independent gantry angle sensor and associated software. The system can either directly calculate or reconstruct (using measured detector responses) a 3D dose distribution on the patient CT dataset for plan verification. Beam models are developed and commissioned in the same manner as a beam model is commissioned in a standard treatment planning system. Model validation was initially performed by comparing both COMPASS calculations and reconstructions to measured open field beam data. Next, 10 hypofractionated prostate RapidArc plans were delivered to both the COMPASS system and a phantom with ion chamber and film inserted. COMPASS dose distributions calculated and reconstructed on the phantom CT dataset were compared to the chamber and film measurements. The mean (± standard deviation) difference between COMPASS reconstructed dose and ion chamber measurement was 1.4 ± 1.0%. The maximum discrepancy was 2.6%. Corresponding values for COMPASS calculation were 0.9 ± 0.9% and 2.6%, respectively. The average gamma agreement index (3%/3mm) for COMPAS reconstruction and film was 96.7% and 95.3% when using 70% and 20% dose thresholds, respectively. The corresponding values for COMPASS calculation were 97.1% and 97.1%, respectively. Based on our results, COMPASS can be used for the patient specific QA of hypofractionated prostate treatments delivered with the 6MV-SRS beam.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Efficient scatter distribution estimation and correction in CBCT using concurrent Monte Carlo fitting

Purpose: X-ray scatter is a significant impediment to image quality improvements in cone-beam CT (CBCT). The authors present and demonstrate a novel scatter correction algorithm using a scatter estimation method that simultaneously combines multiple Monte Carlo (MC) CBCT simulations through the use of a concurrently evaluated fitting function, referred to as concurrent MC fitting (CMCF). Methods: The CMCF method uses concurrently run MC CBCT scatter projection simulations that are a subset of the projection angles used in the projection set, P, to be corrected. The scattered photons reaching the detector in each MC simulation are simultaneously aggregated by an algorithm which computes the scatter detector response, S{sub MC}. S{sub MC} is fit to a function, S{sub F}, and if the fit of S{sub F} is within a specified goodness of fit (GOF), the simulations are terminated. The fit, S{sub F}, is then used to interpolate the scatter distribution over all pixel locations for every projection angle in the set P. The CMCF algorithm was tested using a frequency limited sum of sines and cosines as the fitting function on both simulated and measured data. The simulated data consisted of an anthropomorphic head and a pelvis phantom created from CT data, simulated with and without the use of a compensator. The measured data were a pelvis scan of a phantom and patient taken on an Elekta Synergy platform. The simulated data were used to evaluate various GOF metrics as well as determine a suitable fitness value. The simulated data were also used to quantitatively evaluate the image quality improvements provided by the CMCF method. A qualitative analysis was performed on the measured data by comparing the CMCF scatter corrected reconstruction to the original uncorrected and corrected by a constant scatter correction reconstruction, as well as a reconstruction created using a set of projections taken with a small cone angle. Results: Pearson’s correlation, r, proved to be a suitable GOF metric with strong correlation with the actual error of the scatter fit, S{sub F}. Fitting the scatter distribution to a limited sum of sine and cosine functions using a low-pass filtered fast Fourier transform provided a computationally efficient and accurate fit. The CMCF algorithm reduces the number of photon histories required by over four orders of magnitude. The simulated experiments showed that using a compensator reduced the computational time by a factor between 1.5 and 1.75. The scatter estimates for the simulated and measured data were computed between 35–93 s and 114–122 s, respectively, using 16 Intel Xeon cores (3.0 GHz). The CMCF scatter correction improved the contrast-to-noise ratio by 10%–50% and reduced the reconstruction error to under 3% for the simulated phantoms. Conclusions: The novel CMCF algorithm significantly reduces the computation time required to estimate the scatter distribution by reducing the statistical noise in the MC scatter estimate and limiting the number of projection angles that must be simulated. Using the scatter estimate provided by the CMCF algorithm to correct both simulated and real projection data showed improved reconstruction image quality.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Segmentation-free empirical beam hardening correction for CT

Purpose: The polychromatic nature of the x-ray beams and their effects on the reconstructed image are often disregarded during standard image reconstruction. This leads to cupping and beam hardening artifacts inside the reconstructed volume. To correct for a general cupping, methods like water precorrection exist. They correct the hardening of the spectrum during the penetration of the measured object only for the major tissue class. In contrast, more complex artifacts like streaks between dense objects need other techniques of correction. If using only the information of one single energy scan, there are two types of corrections. The first one is a physical approach. Thereby, artifacts can be reproduced and corrected within the original reconstruction by using assumptions in a polychromatic forward projector. These assumptions could be the used spectrum, the detector response, the physical attenuation and scatter properties of the intersected materials. A second method is an empirical approach, which does not rely on much prior knowledge. This so-called empirical beam hardening correction (EBHC) and the previously mentioned physical-based technique are both relying on a segmentation of the present tissues inside the patient. The difficulty thereby is that beam hardening by itself, scatter, and other effects, which diminish the image quality also disturb the correct tissue classification and thereby reduce the accuracy of the two known classes of correction techniques. The herein proposed method works similar to the empirical beam hardening correction but does not require a tissue segmentation and therefore shows improvements on image data, which are highly degraded by noise and artifacts. Furthermore, the new algorithm is designed in a way that no additional calibration or parameter fitting is needed. Methods: To overcome the segmentation of tissues, the authors propose a histogram deformation of their primary reconstructed CT image. This step is essential for the proposed algorithm to be segmentation-free (sf). This deformation leads to a nonlinear accentuation of higher CT-values. The original volume and the gray value deformed volume are monochromatically forward projected. The two projection sets are then monomially combined and reconstructed to generate sets of basis volumes which are used for correction. This is done by maximization of the image flatness due to adding additionally a weighted sum of these basis images. sfEBHC is evaluated on polychromatic simulations, phantom measurements, and patient data. The raw data sets were acquired by a dual source spiral CT scanner, a digital volume tomograph, and a dual source micro CT. Different phantom and patient data were used to illustrate the performance and wide range of usability of sfEBHC across different scanning scenarios. The artifact correction capabilities are compared to EBHC. Results: All investigated cases show equal or improved image quality compared to the standard EBHC approach. The artifact correction is capable of correcting beam hardening artifacts for different scan parameters and scan scenarios. Conclusions: sfEBHC generates beam hardening-reduced images and is furthermore capable of dealing with images which are affected by high noise and strong artifacts. The algorithm can be used to recover structures which are hardly visible inside the beam hardening-affected regions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

Improved dosimetry for targeted radionuclide therapy using nonrigid registration on sequential SPECT images

Purpose: Voxel-level and patient-specific 3D dosimetry for targeted radionuclide therapy (TRT) typically involves serial nuclear medicine scans. Misalignment of the images can result in reduced dosimetric accuracy. Since the scans are typically performed over a period of several days, there will be patient movement between scans and possible nonrigid organ deformation. This work aims to implement and evaluate the use of nonrigid image registration on a series of quantitative SPECT (QSPECT) images for TRT dosimetry. Methods: A population of 4D extended cardiac torso phantoms, comprised of three In-111 Zevalin biokinetics models and three anatomical variations, was generated based on the patient data. The authors simulated QSPECT acquisitions at five time points. At each time point, individual organ and whole-body deformation between scans were modeled by translating/rotating organs and the body up to 5°/voxels, keeping ≤5% difference in organ volume. An analytical projector was used to generate realistic noisy projections for a medium energy general purpose collimator. Projections were reconstructed using OS-EM algorithm with geometric collimator detector response, attenuation, and scatter corrections. The QSPECT images were registered using organ-based nonrigid image registration method. The cumulative activity in each voxel was obtained by integrating the activity over time. Dose distribution images were obtained by convolving the cumulative activity images with a Y-90 dose kernel. Dose volume histograms (DVHs) for organs-of-interest were analyzed. Results: After nonrigid registration, the mean differences in organ doses compared to the case without misalignment were improved from (−15.50 ± 5.59)% to (−2.12 ± 1.05)% and (−7.28 ± 2.30)% to (−0.23 ± 0.71)% for the spleen and liver, respectively. For all organs, the cumulative DVHs showed improvement after nonrigid registration and the normalized absolute error of differential DVHs ranged from 6.79% to 22.70% for liver and 26.00% to 39.70% for spleen with different segmentation methods. Conclusions: These results demonstrated that nonrigid registration of sequential QSPECT images is feasible for TRT and improves the accuracy of 3D dosimetry.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-C-201-06: Utility of Quantitative 3D SPECT/CT Imaging in Patient Specific Internal Dosimetry of 153-Samarium with GATE Monte Carlo Package

Purpose: Patient-specific 3-dimensional (3D) internal dosimetry in targeted radionuclide therapy is essential for efficient treatment. Two major steps to achieve reliable results are: 1) generating quantitative 3D images of radionuclide distribution and attenuation coefficients and 2) using a reliable method for dose calculation based on activity and attenuation map. In this research, internal dosimetry for 153-Samarium (153-Sm) was done by SPECT-CT images coupled GATE Monte Carlo package for internal dosimetry. Methods: A 50 years old woman with bone metastases from breast cancer was prescribed 153-Sm treatment (Gamma: 103keV and beta: 0.81MeV). A SPECT/CT scan was performed with the Siemens Simbia-T scanner. SPECT and CT images were registered using default registration software. SPECT quantification was achieved by compensating for all image degrading factors including body attenuation, Compton scattering and collimator-detector response (CDR). Triple energy window method was used to estimate and eliminate the scattered photons. Iterative ordered-subsets expectation maximization (OSEM) with correction for attenuation and distance-dependent CDR was used for image reconstruction. Bilinear energy mapping is used to convert Hounsfield units in CT image to attenuation map. Organ borders were defined by the itk-SNAP toolkit segmentation on CT image. GATE was then used for internal dose calculation. The Specific Absorbed Fractions (SAFs) and S-values were reported as MIRD schema. Results: The results showed that the largest SAFs and S-values are in osseous organs as expected. S-value for lung is the highest after spine that can be important in 153-Sm therapy. Conclusion: We presented the utility of SPECT-CT images and Monte Carlo for patient-specific dosimetry as a reliable and accurate method. It has several advantages over template-based methods or simplified dose estimation methods. With advent of high speed computers, Monte Carlo can be used for treatment planning on a day to day basis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-I-08: Investigation of Deconvolution Methods for Blocker-Based CBCT Scatter Estimation

Purpose: To investigate whether deconvolution methods can improve the scatter estimation under different blurring and noise conditions for blocker-based scatter correction methods for cone-beam X-ray computed tomography (CBCT). Methods: An “ideal” projection image with scatter was first simulated for blocker-based CBCT data acquisition by assuming no blurring effect and no noise. The ideal image was then convolved with long-tail point spread functions (PSF) with different widths to mimic the blurring effect from the finite focal spot and detector response. Different levels of noise were also added. Three deconvolution Methods: 1) inverse filtering; 2) Wiener; and 3) Richardson-Lucy, were used to recover the scatter signal in the blocked region. The root mean square error (RMSE) of estimated scatter serves as a quantitative measure for the performance of different methods under different blurring and noise conditions. Results: Due to the blurring effect, the scatter signal in the blocked region is contaminated by the primary signal in the unblocked region. The direct use of the signal in the blocked region to estimate scatter (“direct method”) leads to large RMSE values, which increase with the increased width of PSF and increased noise. The inverse filtering is very sensitive to noise and practically useless. The Wiener and Richardson-Lucy deconvolution methods significantly improve scatter estimation compared to the direct method. For a typical medium PSF and medium noise condition, both methods (∼20 RMSE) can achieve 4-fold improvement over the direct method (∼80 RMSE). The Wiener method deals better with large noise and Richardson-Lucy works better on wide PSF. Conclusion: We investigated several deconvolution methods to recover the scatter signal in the blocked region for blocker-based scatter correction for CBCT. Our simulation results demonstrate that Wiener and Richardson-Lucy deconvolution can significantly improve the scatter estimation compared to the direct method.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

SU-E-T-631: Preliminary Results for Analytical Investigation Into Effects of ArcCHECK Setup Errors

Purpose: As three-dimensional diode arrays increase in popularity for patient-specific quality assurance for intensity-modulated radiation therapy (IMRT), it is important to evaluate an array’s susceptibility to setup errors. The ArcCHECK phantom is set up by manually aligning its outside marks with the linear accelerator’s lasers and light-field. If done correctly, this aligns the ArcCHECK cylinder’s central axis (CAX) with the linear accelerator’s axis of rotation. However, this process is prone to error. This project has developed an analytical expression including a perturbation factor to quantify the effect of shifts. Methods: The ArcCHECK is set up by aligning its machine marks with either the sagittal room lasers or the light-field of the linear accelerator at gantry zero (IEC). ArcCHECK has sixty-six evenly-spaced SunPoint diodes aligned radially in a ring 14.4 cm from CAX. The detector response function (DRF) was measured and combined with inverse-square correction to develop an analytical expression for output. The output was calculated using shifts of 0 (perfect alignment), +/−1, +/−2 and +/−5 mm. The effect on a series of simple inputs was determined: unity, 1-D ramp, steps, and hat-function to represent uniform field, wedge, evenly-spaced modulation, and single sharp modulation, respectively. Results: Geometric expressions were developed with perturbation factor included to represent shifts. DRF was modeled using sixth-degree polynomials with correlation coefficient 0.9997. The output was calculated using simple inputs such as unity, 1-D ramp, steps, and hat-function, with perturbation factors of: 0, +/−1, +/−2 and +/−5 mm. Discrepancies have been observed, but large fluctuations have been somewhat mitigated by aliasing arising from discrete diode placement. Conclusion: An analytical expression with perturbation factors was developed to estimate the impact of setup errors on an ArcCHECK phantom. Presently, this has been applied to simple functions, but is being developed to handle more realistic clinical cases.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2015

SU-E-T-758: To Determine the Source Dwell Positions of HDR Brachytherapy Using 2D 729 Ion Chamber Array

Purpose: Determination of source dwell positions of HDR brachytherapy using 2D 729 ion chamber array Methods: Nucletron microselectron HDR and PTW 2D array were used for the study. Different dwell positions were assigned in the HDR machine. Rigid interstitial needles and vaginal applicator were positioned on the 2D array. The 2D array was exposed for this programmed dwell positions. The positional accuracy of the source was analyzed after the irradiation of the 2D array. This was repeated for different dwell positions. Different test plans were transferred from the Oncentra planning system and irradiated with the same applicator position on the 2D array. The results were analyzed using the in house developed excel program. Results: Assigned dwell positions versus corresponding detector response were analyzed. The results show very good agreement with the film measurements. No significant variation found between the planned and measured dwell positions. Average dose response with 2D array between the planned and nearby dwell positions was found to be 0.0804 Gy for vaginal cylinder applicator and 0.1234 Gy for interstitial rigid needles. Standard deviation between the doses for all the measured dwell positions for interstitial rigid needle for 1 cm spaced positions were found to be 0.33 and 0.37 for 2cm spaced dwell positions. For intracavitory vaginal applicator this was found to be 0.21 for 1 cm spaced dwell positions and 0.06 for 2cm spaced dwell positions. Intracavitory test plans reproduced on the 2D array with the same applicator positions shows the ideal dose distribution with the TPS planned. Conclusion: 2D array is a good tool for determining the dwell position of HDR brachytherapy. With the in-house developed program in excel it is easy and accurate. The traditional way with film analysis can be replaced by this method, as the films will be more costly.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Permanent-magnet energy spectrometer for electron beams from radiotherapy accelerators

Purpose: The purpose of this work was to adapt a lightweight, permanent magnet electron energy spectrometer for the measurement of energy spectra of therapeutic electron beams. Methods: An irradiation geometry and measurement technique were developed for an approximately 0.54-T, permanent dipole magnet spectrometer to produce suitable latent images on computed radiography (CR) phosphor strips. Dual-pinhole electron collimators created a 0.318-cm diameter, approximately parallel beam incident on the spectrometer and an appropriate dose rate at the image plane (CR strip location). X-ray background in the latent image, reduced by a 7.62-cm thick lead block between the pinhole collimators, was removed using a fitting technique. Theoretical energy-dependent detector response functions (DRFs) were used in an iterative technique to transform CR strip net mean dose profiles into energy spectra on central axis at the entrance to the spectrometer. These spectra were transformed to spectra at 95-cm source to collimator distance (SCD) by correcting for the energy dependence of electron scatter. The spectrometer was calibrated by comparing peak mean positions in the net mean dose profiles, initially to peak mean energies determined from the practical range of central-axis percent depth-dose (%DD) curves, and then to peak mean energies that accounted for how the collimation modified the energy spectra (recalibration). The utility of the spectrometer was demonstrated by measuring the energy spectra for the seven electron beams (7–20 MeV) of an Elekta Infinity radiotherapy accelerator. Results: Plots of DRF illustrated their dependence on energy and position in the imaging plane. Approximately 15 iterations solved for the energy spectra at the spectrometer entrance from the measured net mean dose profiles. Transforming those spectra into ones at 95-cm SCD increased the low energy tail of the spectra, while correspondingly decreasing the peaks and shifting them to slightly lower energies. Energy calibration plots of peak mean energy versus peak mean position of the net mean dose profiles for each of the seven electron beams followed the shape predicted by the Lorentz force law for a uniform z-component of the magnetic field, validating its being modeled as uniform (0.542 ± 0.027 T). Measured Elekta energy spectra and their peak mean energies correlated with the 0.5-cm (7–13 MeV) and the 1.0-cm (13–20 MeV) R{sub 90} spacings of the %DD curves. The full-width-half-maximum of the energy spectra decreased with decreasing peak mean energy with the exception of the 9-MeV beam, which was anomalously wide. Similarly, R{sub 80–20} decreased linearly with peak mean energy with the exception of the 9 MeV beam. Both were attributed to suboptimal tuning of the high power phase shifter for the recycled radiofrequency power reentering the traveling wave accelerator. Conclusions: The apparatus and analysis techniques of the authors demonstrated that an inexpensive, lightweight, permanent magnet electron energy spectrometer can be used for measuring the electron energy distributions of therapeutic electron beams (6–20 MeV). The primary goal of future work is to develop a real-time spectrometer by incorporating a real-time imager, which has potential applications such as beam matching, ongoing beam tune maintenance, and measuring spectra for input into Monte Carlo beam calculations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-J-08: Quantitative SPECT Imaging of Ra-223 in a Phantom

Purpose: Ra-223 therapy of prostate cancer bone metastases is being used to treat patients routinely. However, the absorbed dose distribution at the macroscopic and microscopic scales remains elusive, due to the inability to image the small activities injected. Accurate activity quantification through imaging is essential to calculate the absorbed dose in organs and sub-units in radiopharmaceutical therapy, enabling personalized absorbed dose-based treatment planning methodologies and more effective and optimal treatments. Methods: A 22 cm diameter by 20 cm long cylindrical phantom, containing a 3.52 cm diameter sphere, was used. A total of 2.01 MBq of Ra-223 was placed in the phantom with 177.6 kBq in the sphere. Images were acquired on a dual-head Siemens Symbia T16 gamma camera using three 20% full-width energy windows and centered at 84, 154, and 269 keV (120 projections, 360° rotation, 45 s per view). We have implemented reconstruction of Ra-223 SPECT projections using OS-EM (up to 20 iterations of 10 subsets) with compensation for attenuation using CT-based attenuation maps, collimator-detector response (CDR) (including septal penetration, scatter and Pb x-ray modeling), and scatter in the patient using the effective source scatter estimation (ESSE) method. The CDR functions and scatter kernels required for ESSE were computed using the SIMIND MC simulation code. All Ra-223 photon emissions as well as gamma rays from the daughters Rn-219 and Bi-211 were modeled. Results: The sensitivity of the camera in the three combined windows was 107.3 cps/MBq. The visual quality of the SPECT images was reasonably good and the activity in the sphere was 27% smaller than the true activity. This underestimation is likely due to partial volume effect. Conclusion: Absolute quantitative Ra-223 SPECT imaging is achievable with careful attention to compensate for image degrading factors and system calibration.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-F-T-570: Comparison of Synthetic Diamond, Microionization Chamber, and Radiochromic Film for Absolute Dosimetry of VMAT Radiosurgery

Purpose: Patient specific quality assurance of stereotactic radiosurgery (SRS) plans is challenging because of small target sizes and high dose gradients. We compared three detectors for dosimetry of VMAT SRS plans. Methods: The dose at the center of seventeen targets was measured using a synthetic diamond detector (2.2 mm diameter, 1 µm thickness), a 0.007 cm{sup 3} ionization chamber, and radiochromic film. Measurements were made in a PMMA phantom in the clinical geometry – all gantry and table angles were delivered as planned. The diamond and chamber positions were offset by 1 cm from the film plane, so the isocenter was shifted accordingly to place the center of the target at the detector of interest. To ensure accurate detector placement, the phantom was positioned using kV images. To account for the shift-induced difference in geometry and differing prescription doses between plans, the measurements were normalized to the expected dose calculated by the treatment planning system. Results: The target sizes ranged from 2.8 mm to 34.8 mm (median 14.8 mm). The mean measurement-to-plan ratios were 1.054, 1.076, and 1.023 for RCF, diamond, and chamber, respectively. The mean difference between the chamber and film was −3.2% and between diamond and film was 2.2%. For targets larger than 15 mm, the mean difference relative to film was −0.8% and 0.1% for chamber and diamond, respectively, whereas for targets smaller than 15 mm, the difference was −5.3% and 4.2% for chamber and diamond, respectively. The difference was significant (p=0.005) using the two-sample Kolmogorov-Smirnov test. Conclusion: The detectors agree for target sizes larger than 15 mm. Relative to film, for smaller targets the diamond detector over-responds, whereas the ionization chamber under-responds. Further work is needed to characterize detector response in modulated SRS fields.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-JeP1-13: Innovative Tracking Detector for Dose Monitoring in Hadron Therapy: Realization and Monte Carlo Simulations

Purpose: Development of strategies to monitor range uncertainties is necessary to improve treatment planning in Charged Particle Therapy (CPT) and fully exploit the advantages of ion beams. Our group developed (within the framework of the INSIDE project funded by the Italian research ministry) and is currently building a compact detector Dose Profiler (DP) able to backtrack charged secondary particles produced in the patient during the irradiation. Furthermore we are studying monitoring strategy exploiting charged secondary emission profiles to control the range of the ion beam. Methods: This contribution reports on the DP detector design and construction status. The detector consists of a charged secondary tracker composed of scintillating fiber layers and a LYSO calorimeter for particles energy measurement.The detector layout has been optimized using the FLUKA Monte Carlo (MC) simulation software. The simulation of a 220 MeV Carbon beam impinging on a PMMA target has been performed to study the detector response, exploiting previous secondary radiation measurements performed by our group. The emission profile of charged secondary particles was reconstructed backtracking the particles to their generation point to benchmark the DP performances. Results: The DP construction status, including the technological details will be presented. The feasibility of range monitoring with DP will be demonstrated by means of MC studies. The correlation of the charged secondary particles emission shape with the position of the Bragg peak (BP) will be shown, as well as the spatial resolution achievable on the BP position estimation (less than 3 mm) in the clinical like conditions. Conclusion: The simulation studies supported the feasibility of an accurate range monitoring technique exploiting the use of charged secondary fragments emitted during the particle therapy treatment. The DP experimental tests are foreseen in 2016, at CNAO particle therapy center in Pavia.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

SU-G-TeP4-07: Automatic EPID-Based 2D Measurement of MLC Leaf Offset as a Quality Control Tool

Purpose: The MLC dosimetric leaf gap (DLG) and transmission are measured parameters which impact the dosimetric accuracy of IMRT and VMAT plans. This investigation aims to develop an efficient and accurate routine constancy check of the physical DLG in two dimensions. Methods: The manufacturer’s recommended DLG measurement method was modified by using 5 fields instead of 11 and by utilizing the Electronic Portal Imaging Device (EPID). Validations were accomplished using an ion chamber (IC) in solid water and a 2D IC array. EPID data was collected for 6 months on multiple TrueBeam linacs using both Millennium and HD MLCs at 5 different clinics in an international consortium. Matlab code was written to automatically analyze the images and calculate the 2D results. Sensitivity was investigated by introducing deliberate leaf position errors. MLC calibration and initialization history was recorded to allow quantification of their impact. Results were analyzed using statistical process control (SPC). Results: The EPID method took approximately 5 minutes. Due to detector response, the EPID measured DLG and transmission differed from the IC values but were reproducible and consistent with changes measured using the ICs. For the Millennium MLC, the EPID measured DLG and transmission were both consistently lower than IC results. The EPID method was implemented as leaf offset and transmission constancy tests (LOC and TC). Based on 6 months of measurements, the initial leaf-specific action thresholds for changes from baseline were set to 0.1 mm. Upper and lower control limits for variation were developed for each machine. Conclusion: Leaf offset and transmission constancy tests were implemented on Varian HD and Millennium MLCs using an EPID and found to be efficient and accurate. The test is effective for monitoring MLC performance using dynamic delivery and performing process control on the DLG in 2D, thus enhancing dosimetric accuracy. This work was supported by a grant from Varian Medical Systems.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2016

Sci-Sat AM: Radiation Dosimetry and Practical Therapy Solutions - 05: Not all geometries are equivalent for magnetic field Fano cavity tests

The coupling of MRI and radiation treatment systems for the application of magnetic resonance guided radiation therapy necessitates a reliable magnetic field capable Monte Carlo (MC) code. In addition to the influence of the magnetic field on dose distributions, the question of proper calibration has arisen due to the several percent variation of ion chamber and solid state detector responses in magnetic fields when compared to the 0 T case (Reynolds et al., Med Phys, 2013). In the absence of a magnetic field, EGSnrc has been shown to pass the Fano cavity test (a rigorous benchmarking tool of MC codes) at the 0.1 % level (Kawrakow, Med.Phys, 2000), and similar results should be required of magnetic field capable MC algorithms. To properly test such developing MC codes, the Fano cavity theorem has been adapted to function in a magnetic field (Bouchard et al., PMB, 2015). In this work, the Fano cavity test is applied in a slab and ion-chamber-like geometries to test the transport options of an implemented magnetic field algorithm in EGSnrc. Results show that the deviation of the MC dose from the expected Fano cavity theory value is highly sensitive to the choice of geometry, and the ion chamber geometry appears to pass the test more easily than larger slab geometries. As magnetic field MC codes begin to be used for dose simulations and correction factor calculations, care must be taken to apply the most rigorous Fano test geometries to ensure reliability of such algorithms.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2012

Topological cross sections in proton-proton interactions at 50 GeV

Project Thermalization (Experiment SERP-E-190 at IHEP) is aimed to study the pp interactions at 50 GeV with large number of secondary particles. In this report the experimentally measured topological cross sections have been presented taking into account the detector response and procession effciency. These data are in good agreement with gluon dominance model. The comparison with other models is also made and shows no essential discrepancies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2015

Detection of reflected Cherenkov light from extensive air showers in the SPHERE experiment as a method of studying superhigh energy cosmic rays

Although a large number of experiments were carried out during the last few decades, the uncertainty in the spectrum of all nuclei of primary cosmic rays (PCRs) with superhigh energies is still high, and the results of many experiments on nuclear composition of PCRs are contradictory. An overview of the SPHERE experiment on detecting Vavilov-Cherenkov radiation from extensive air shower (EAS) reflected from a ground snow surface is given. A number of experimental studies implementing this method are presented and their results are analyzed. Some other popular methods of studying PCRs with superhigh energies (E{sub 0} > 10{sup 15} eV) and their main advantages and drawbacks are briefly considered. The detecting equipment of the SPHERE-2 experiment and the technique of its calibration are considered. The optical properties of snow, which are important for experiments on reflected Cherenkov light (CL) from EAS, are discussed and the history of observing reflected EAS CL is described. The algorithm of simulating the detector response and calculating the fiducial acceptance of shower detection is described. The procedure of processing the experimental data with a subsequent reconstruction of the spectrum of all PCR nuclei and analysis of the mass composition is shown. The first results of reconstructing the spectrum and separating groups of cosmic-ray nuclei with high energies in the SPHERE-2 experiment are presented. Main sources of systematic errors are considered. The prospects of developing the technique of observation of reflected EAS CL in future experiments are discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗