Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Error analysis”

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 937 records · Page 52

Subtlety of TiO 2 phase stability: Reliability of the density functional theory predictions and persistence of the self-interaction error

TiO 2 is an important material with broad applications that can exist in different phases with dramatically different properties. Theoretical prediction of their polymorph energetics is therefore critical for the material design and for identifying thermodynamically accessible structures. Determining TiO 2 relative phase stabilities remains challenging for first-principles methods, and density functional theory is the only approach available for studying phase stabilities at finite temperatures with acceptable computational efficiency. Here, we show that density functional theory using the recently developed efficient strongly constrained and appropriately normed (SCAN) [Sun et al., Phys. Rev. Lett. 115, 036402 (2015)] exchange-correlation functional for the first time predicts the phase stability in qualitative agreement with the experimental results at realistic conditions. Further analysis shows that the self-interaction error intrinsic in the density functional persists in the stability prediction. By correcting the self-interaction error through an empirical approach, SCAN predicts the relative stability as well as defect properties in excellent agreement with the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Framework for Error-Bounded Approximate Computing, with an Application to Dot Products

Approximate computing techniques, which trade off the computation accuracy of an algorithm for better performance and energy efficiency, have been successful in reducing computation and power costs in several domains. However, error sensitive applications in high-performance computing are unable to benefit from existing approximate computing strategies that are not developed with guaranteed error bounds. While approximate computing techniques can be developed for individual high-performance computing applications by domain specialists, this often requires additional theoretical analysis and potentially extensive software modification. Hence, the development of low-level error-bounded approximate computing strategies that can be introduced into any high-performance computing application without requiring additional analysis or significant software alterations is desirable. In this paper, we provide a contribution in this direction by proposing a general framework for designing error-bounded approximate computing strategies and apply it to the dot product kernel to develop \bf qdot---an error-bounded approximate dot product kernel. Following the introduction of qdot, here we perform a theoretical analysis that yields a deterministic bound on the relative approximation error introduced by qdot. Empirical tests are performed to illustrate the tightness of the derived error bound and to demonstrate the effectiveness of qdot on a synthetic dataset, as well as two scientific benchmarks---the conjugate gradient (CG) and power methods. In some instances, using qdot for the dot products in CG can result in many components being quantized to half precision without increasing the iteration count required for convergence to the same solution as CG using a double precision dot product.

97 MATHEMATICS AND COMPUTING↗

Emulating Baryonic Effects in Weak Lensing Measurements

Baryonic effects are one of the most serious systematics to the tomographic analysis of weak lensing statistical measurements. We are working on a novel method for correcting cosmological observations for baryonic effects - a term used to label the fact that we observe galaxies, while theoretical modeling at present is only accurate enough for the dynamics of dark matter. Although modern cosmological simulations of dark matter (the “dark matter only", or DMO simulations) consistently reach the required 1% level of precision, the effect of baryonic physics on cosmological scales is expected to be a source of a much more severe systematic error in weak lensing measurements. In this analysis we apply a general framework to consistently include systematic effects into a likelihood analysis for cosmological surveys, such as the Dark Energy Survey (DES), and future Large Synoptic Survey Telescope survey (LSST).

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of over-magnetization of elemental transition metal solids from the SCAN density functional

Current investigations have revealed that the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation exchange-correlation functional significantly over-magnetizes elemental Fe, Co, and Ni solids. For the paradigmatic case, bcc Fe, the error relative to experiment is ≳ 20 % . A comparative analysis of magnetization results from SCAN and its deorbitalized counterpart, SCAN-L, leads to identification of the source of the discrepancy. It is not from the difference between Kohn-Sham (SCAN-L) and generalized Kohn-Sham (SCAN) procedures. The key is the iso-orbital indicator α (the ratio of the local Pauli and Thomas-Fermi kinetic energy densities). Its deorbitalized counterpart α L has more dispersion in both spin channels with respect to magnetization in an approximate region between 0.6 and 1.2 bohrs around an Fe nucleus. The overall effect is that the SCAN switching function evaluated with α L reduces the energetic disadvantage of the down channel with respect to up compared to the original α , which in turn reduces the magnetization. This identifies the origin of the SCAN magnetization error as insensitivity of the SCAN switching function to α values in the approximate range 0.5 ≲ α ≲ 0.8 and oversensitivity for α ≳ 0.8 .

36 MATERIALS SCIENCE↗

Comparative analysis of image binarization methods for crack identification in concrete structures

Surface cracks in concrete structures are critical indicators of structural damage and durability. Manual visual inspection, the most commonly used method in practice, is inefficient from cost, time, accuracy, and safety perspectives. A promising alternative is computer vision-based methods that can automatically extract crack information from images. Image binarization, developed for text detection, is appropriate for crack identification, as texts and cracks are similar, consisting of distinguishable lines and curves. However, standardizing crack identification using image binarization is challenging, because binarization depends on the method and associated parameters. We investigate image binarization for crack identification, focusing on optimal parameter determination and comparative performance evaluation for five common binarization methods. Crack images are prepared to obtain optimal parameters by minimizing errors in estimated crack widths. Subsequently, comparative analysis is conducted using crack images with different conditions based on three performance evaluation criteria: crack width and length measurement accuracy and computation time.

36 MATERIALS SCIENCE↗

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: BAO and RSD measurements from the anisotropic power spectrum of the quasar sample between redshift 0.8 and 2.2

ABSTRACT We measure the clustering of quasars of the final data release (DR16) of eBOSS. The sample contains $343\, 708$ quasars between redshifts 0.8 ≤ z ≤ 2.2 over $4699\, \mathrm{deg}^2$. We calculate the Legendre multipoles (0,2,4) of the anisotropic power spectrum and perform a BAO and a Full-Shape (FS) analysis at the effective redshift zeff = 1.480. The errors include systematic errors that amount to 1/3 of the statistical error. The systematic errors comprise a modelling part studied using a blind N-body mock challenge and observational effects studied with approximate mocks to account for various types of redshift smearing and fibre collisions. For the BAO analysis, we measure the transverse comoving distance DM(zeff)/rdrag = 30.60 ± 0.90 and the Hubble distance DH(zeff)/rdrag = 13.34 ± 0.60. This agrees with the configuration space analysis, and the consensus yields: DM(zeff)/rdrag = 30.69 ± 0.80 and DH(zeff)/rdrag = 13.26 ± 0.55. In the FS analysis, we fit the power spectrum using a model based on Regularised Perturbation Theory, which includes redshift space distortions and the Alcock–Paczynski effect. The results are DM(zeff)/rdrag = 30.68 ± 0.90 and DH(zeff)/rdrag = 13.52 ± 0.51 and we constrain the linear growth rate of structure f(zeff)σ8(zeff) = 0.476 ± 0.047. Our results agree with the configuration space analysis. The consensus analysis of the eBOSS quasar sample yields: DM(zeff)/rdrag = 30.21 ± 0.79, DH(zeff)/rdrag = 3.23 ± 0.47, and f(zeff)σ8(zeff) = 0.462 ± 0.045 and is consistent with a flat ΛCDM cosmological model using Planck results.

Neveux, Richard↗

Image-based spectral distortion correction for photon-counting x-ray detectors

Purpose: To investigate the feasibility of using an image-based method to correct for distortions induced by various artifacts in the x-ray spectrum recorded with photon-counting detectors for their application in breast computed tomography (CT). Methods: The polyenergetic incident spectrum was simulated with the tungsten anode spectral model using the interpolating polynomials (TASMIP) code and carefully calibrated to match the x-ray tube in this study. Experiments were performed on a Cadmium-Zinc-Telluride (CZT) photon-counting detector with five energy thresholds. Energy bins were adjusted to evenly distribute the recorded counts above the noise floor. BR12 phantoms of various thicknesses were used for calibration. A nonlinear function was selected to fit the count correlation between the simulated and the measured spectra in the calibration process. To evaluate the proposed spectral distortion correction method, an empirical fitting derived from the calibration process was applied on the raw images recorded for polymethyl methacrylate (PMMA) phantoms of 8.7, 48.8, and 100.0 mm. Both the corrected counts and the effective attenuation coefficient were compared to the simulated values for each of the five energy bins. The feasibility of applying the proposed method to quantitative material decomposition was tested using a dual-energy imaging technique with a three-material phantom that consisted of water, lipid, and protein. The performance of the spectral distortion correction method was quantified using the relative root-mean-square (RMS) error with respect to the expected values from simulations or areal analysis of the decomposition phantom. Results: The implementation of the proposed method reduced the relative RMS error of the output counts in the five energy bins with respect to the simulated incident counts from 23.0%, 33.0%, and 54.0% to 1.2%, 1.8%, and 7.7% for 8.7, 48.8, and 100.0 mm PMMA phantoms, respectively. The accuracy of the effective attenuation coefficient of PMMA estimate was also improved with the proposed spectral distortion correction. Finally, the relative RMS error of water, lipid, and protein decompositions in dual-energy imaging was significantly reduced from 53.4% to 6.8% after correction was applied. Conclusions: The study demonstrated that dramatic distortions in the recorded raw image yielded from a photon-counting detector could be expected, which presents great challenges for applying the quantitative material decomposition method in spectral CT. The proposed semi-empirical correction method can effectively reduce these errors caused by various artifacts, including pulse pileup and charge sharing effects. Furthermore, rather than detector-specific simulation packages, the method requires a relatively simple calibration process and knowledge about the incident spectrum. Therefore, it may be used as a generalized procedure for the spectral distortion correction of different photon-counting detectors in clinical breast CT systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Quantum process tomography of a Mølmer-Sørensen gate via a global beam

Here, we present a framework for quantum process tomography of two-ion interactions that leverages modulations of the trapping potential and composite pulses from a global laser beam to achieve individual-ion addressing. Tomographic analysis of identity and delay processes reveals dominant error contributions from laser decoherence and slow qubit frequency drift during the tomography experiment. We use this framework on two co-trapped 40 Ca + ions to analyze both an optimized and an overpowered Mølmer–Sørensen gate and to compare the results of this analysis to a less informative Bell-state tomography measurement and to predictions based on a simplified noise model. These results show that the technique is effective for the characterization of two-ion quantum processes and for the extraction of meaningful information about the errors present in the system. The experimental convenience of this method will allow for more widespread use of process tomography for characterizing entangling gates in trapped-ion systems.

quantum computing↗

Lightweight and Statistical Techniques for Petascale PetaScale Debugging

This project investigated novel techniques for debugging scientific applications on petascale architectures. In particular, we developed lightweight tools that narrow the problem space when bugs are encountered. We also developed techniques that either limit the number of tasks and the code regions to which a developer must apply a traditional debugger or that apply statistical techniques to provide direct suggestions of the location and type of error. We extend previous work on the Stack Trace Analysis Tool (STAT), that has already demonstrated scalability to over one hundred thousand MPI tasks. We also extended statistical techniques developed to isolate programming errors in widely used sequential or threaded applications in the Cooperative Bug Isolation (CBI) project to large scale parallel applications. Overall, our research substantially improved productivity on petascale platforms through a tool set for debugging that complements existing commercial tools.

97 MATHEMATICS AND COMPUTING↗

In vivo dosimetry for gynaecological brachytherapy using a novel position sensitive radiation detector: Feasibility study

Purpose: In gynecological radiotherapy with high dose rate (HDR){sup 192}Ir brachytherapy, the treatment complexity has increased due to improved optimization techniques and dose constraints. As a consequence, it has become more important to verify the dose delivery to the target and also to the organs at risk (e.g., the bladder). In vivo dosimetry, where dosimeters are placed in or on the patient, is one way of verifying the dose but until recently this was hampered by motion of the radiation detectors with respect to the source. The authors present a novel dosimetry method using a position sensitive radiation detector. Methods: The prototype RADPOS system (Best Medical Canada) consists of a metal oxide field effect transistor (MOSFET) dosimeter coupled to a position-sensor, which deduces its 3D position in a magnetic field. To assess the feasibility of in vivo dosimetry based on the RADPOS system, different characteristics of the detector need to be investigated. Using a PMMA phantom, the positioning accuracy of the RADPOS system was quantified by comparing position readouts with the known position of the detector along the x and y-axes. RADPOS dose measurements were performed at various distances from a Nucletron{sup 192}Ir source in a PMMA phantom to evaluate the energy dependence of the MOSFET. A sensitivity analysis was performed by calculating the dose after varying (1) the position of the RADPOS detector to simulate organ motion and (2) the position of the first dwell position to simulate errors in delivery. The authors also performed an uncertainty analysis to determine the action level (AL) that should be used during in vivo dosimetry. Results: Positioning accuracy is found to be within 1 mm in the 1-10 cm range from the origin along the x-axis (away from the transmitter), meeting the requirements for in vivo dosimetry. Similar results are obtained for the other axes. The ALs are chosen to take into account the total uncertainty on the measurements. As a consequence for in vivo dosimetry, it is determined that the RADPOS sensor, if placed, for example, in the bladder Foley balloon, would detect a 2 mm motion of the bladder, at a 5% chance of a false positive, with an AL limit of 9% of the dose delivered. The authors found that source position errors, caused by, e.g., a wrong first dwell position, are more difficult to detect; indeed, with our single RADPOS detector, positioned in the bladder, dwell position errors below 5 mm and resulting in a dose error within 10%, could be detected in the tandem but not in the colpostats. A possible solution to improve error detection is to use multiple MOSFETs to obtain multiple dose values. Conclusions: In this study, the authors proposed a dosimetry procedure, based on the novel RADPOS system, to accurately determine the position of the radiation dosimeter with respect to the applicator. The authors found that it is possible to monitor the delivered dose in a point and compare it to the predetermined dose. This allows in principle the detection of problems such as bladder motion/filling or source mispositioning. Further clinical investigation is warranted.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Eigenstate thermalization hypothesis and approximate quantum error correction

The eigenstate thermalization hypothesis (ETH) is a powerful conjecture for understanding how statistical mechanics emerges in a large class of many-body quantum systems. It has also been interpreted in a CFT context, and, in particular, holographic CFTs are expected to satisfy ETH. Recently, it was observed that the ETH condition corresponds to a necessary and sufficient condition for an approximate quantum error correcting code (AQECC), implying the presence of AQECCs in systems satisfying ETH. In this paper, we explore the properties of ETH as an error correcting code and show that there exists an explicit universal recovery channel for the code. Based on the analysis, we discuss a generalization that all chaotic theories contain error correcting codes. We then specialize to AdS/CFT to demonstrate the possibility of total bulk reconstruction in black holes with a well-defined macroscopic geometry. When combined with the existing AdS/CFT error correction story, this shows that black holes are enormously robust against erasure errors.

97 MATHEMATICS AND COMPUTING↗

TU-F-18C-01: Breast Tissue Decomposition Using Spectral CT After Distortion Correction

Purpose: To investigate the feasibility of accurate breast tissue compositional characterization by using spectral-distortion-corrected dual energy images from a photon-counting spectral CT. Methods: Thirty eight postmortem breasts were imaged with a Cadmium-Zinc-Telluride (CZT)-based photon-counting spectral CT system at beam energy of 100 kVp. The energy-resolved detector sorted photons into low and high energy bins with a splitting energy of 42 keV. The estimated mean glandular dose (MGD) for each breast was approximately 2.0 mGy. Dual energy technique was used to decompose breast tissue into water, lipid, and protein contents. Two image-based methods were investigated to improve the accuracy of tissue compositional characterization. The first method simply limited the recorded spectra up to 90 keV. This reduced the pulse pile-up artifacts but it has some dose penalty. The second method corrected the spectral information of all measured photons by using a spectral distortion correction technique. Breasts were then chemically decomposed into their respective water, lipid, and protein contents, which was used as the reference standard. The accuracy of the tissue compositional measurement with spectral CT was evaluated by the root-mean-square (RMS) errors in percentage composition. Results: The errors in quantitative material decomposition were significantly reduced after the appropriate image processing methods. As compared to the chemical analysis as the reference standard, the averages of the RMS errors were estimated to be 15.5%, 3.3%, and 2.8% for the raw, energy-limited, and spectral-corrected images, respectively. Conclusion: Spectral CT can be used to accurately quantify the water, lipid, and protein contents in breast tissues by implementing a spectral distortion correction algorithm. The tissue compositional information can potentially improve the sensitivity and specificity for breast cancer diagnosis.

60 APPLIED LIFE SCIENCES↗

Harmonic analysis for the characterization and correction of geometric distortion in MRI

Purpose: Magnetic resonance imaging (MRI) is gaining widespread use in radiation therapy planning, patient setup verification, and real-time guidance of radiation delivery. Successful implementation of these technologies relies on the development of simple and efficient methods to characterize and monitor the geometric distortions arising due to system imperfections and gradient nonlinearities. To this end, the authors present the theory and validation of a novel harmonic approach to the quantification of system-related distortions in MRI. Methods: The theory of spatial encoding in MRI is applied to demonstrate that the 3D distortion vector field (DVF) is given by the solution of a second-order boundary value problem (BVP). This BVP is comprised of Laplace’s equation and a limited measurement of the distortion on the boundary of a specified region of interest (ROI). An analytical series expansion solving this BVP within a spherical ROI is obtained, and a statistical uncertainty analysis is performed to determine how random errors in the boundary measurements propagate to the ROI interior. This series expansion is then evaluated to obtain volumetric DVF mappings that are compared to reference data obtained on a 3 T full-body scanner. This validation is performed within two spheres of 20 cm diameter (one centered at the scanner origin and the other offset +3 cm along each of the transverse directions). Initially, a high-order mapping requiring measurements at 5810 boundary points is used. Then, after exploring the impact of the boundary sampling density and the effect of series truncation, a reduced-order mapping requiring measurements at 302 boundary points is evaluated. Results: The volumetric DVF mappings obtained from the harmonic analysis are in good agreement with the reference data. Following distortion correction using the high-order mapping, the authors estimate a reduction in the mean distortion magnitude from 0.86 to 0.42 mm and from 0.93 to 0.39 mm within the central and offset ROIs, respectively. In addition, the fraction of points with a distortion magnitude greater than 1 mm is reduced from 35.6% to 2.8% and from 40.4% to 1.5%, respectively. Similarly, following correction using the reduced-order mapping, the mean distortion magnitude reduces to 0.45–0.42 mm within the central and offset ROIs, and the fraction of points with a distortion magnitude greater than 1 mm is reduced to 2.8% and 1.5%, respectively. Conclusions: A novel harmonic approach to the characterization of system-related distortions in MRI is presented. This method permits a complete and accurate mapping of the DVF within a specified ROI using a limited measurement of the distortion on the ROI boundary. This technique eliminates the requirement to exhaustively sample the DVF at a dense 3D array of points, thereby permitting the design of simple, inexpensive phantoms that may incorporate additional modules for auxiliary quality assurance objectives.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

SU-F-P-07: Applying Failure Modes and Effects Analysis to Treatment Planning System QA

Purpose: A small-scale implementation of Failure Modes and Effects Analysis (FMEA) for treatment planning system QA by utilizing methodology of AAPM TG-100 report. Methods: FMEA requires numerical values for severity (S), occurrence (O) and detectability (D) of each mode of failure. The product of these three values gives a risk priority number (RPN). We have implemented FMEA for the treatment planning system (TPS) QA for two clinics which use Pinnacle and Eclipse TPS. Quantitative monthly QA data dating back to 4 years for Pinnacle and 1 year for Eclipse have been used to determine values for severity (deviations from predetermined doses at points or volumes), and occurrence of such deviations. The TPS QA protocol includes a phantom containing solid water and lung- and bone-equivalent heterogeneities. Photon and electron plans have been evaluated in both systems. The dose values at multiple distinct points of interest (POI) within the solid water, lung, and bone-equivalent slabs, as well as mean doses to several volumes of interest (VOI), have been re-calculated monthly using the available algorithms. Results: The computed doses vary slightly month-over-month. There have been more significant deviations following software upgrades, especially if the upgrade involved re-modeling of the beams. TG-100 guidance and the data presented here suggest an occurrence (O) of 2 depending on the frequency of re-commissioning the beams, severity (S) of 3, and detectability (D) of 2, giving an RPN of 12. Conclusion: Computerized treatment planning systems could pose a risk due to dosimetric errors and suboptimal treatment plans. The FMEA analysis presented here suggests that TPS QA should immediately follow software upgrades, but does not need to be performed every month.

60 APPLIED LIFE SCIENCES↗

The updated ITPA global H-mode confinement database: description and analysis

The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.

Database↗

Easy Web Interfaces to IDL Code for NSTX Data Analysis

Reusing code is a well-known Software Engineering practice to substantially increase the efficiency of code production, as well as to reduce errors and debugging time. A variety of "Web Tools" for the analysis and display of raw and analyzed physics data are in use on NSTX [1], and new ones can be produced quickly from existing IDL [2] code. A Web Tool with only a few inputs, and which calls an IDL routine written in the proper style, can be created in less than an hour; more typical Web Tools with dozens of inputs, and the need for some adaptation of existing IDL code, can be working in a day or so. Efficiency is also increased for users of Web Tools because o f the familiar interface of the web browser, and not needing X-windows, accounts, passwords, etc. Web Tools were adapted for use by PPPL physicists accessing EAST data stored in MDSplus with only a few man-weeks of effort; adapting to additional sites should now be even easier. An overview of Web Tools in use on NSTX, and a list of the most useful features, is also presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-E-T-798: Verification of 3DVH Dose Distribution Before Clinical Implementation for Patient-Specific IMRT QA

Purpose: In recent years patient-specific IMRT QA has transitioned from film and chamber measurements to beam-by-beam 2D array measurements. 3DVH takes this transition a step further by estimating the 3D dose distribution delivered using 2D per beam diode array measurements. In this study, the 3D dose distribution generated by 3DVH is compared to film and chamber measurements. In addition, the accuracy ROI volume and error detection is investigated. Methods: Composite film and ion chamber measurements in a solid water phantom were performed for 9 IMRT PINNACLE patient plans for 4 treatment sites. The film and chamber measurements were compared to the dose distribution predicted by 3DVH using MAPCHECK2 per beam measurements. The absolute point dose measurement (CAX) was extracted from the predicted 3DVH and PINNACLE dose distribution and was compared by taking the ratio of measured to predicted doses. The dose distribution measured with film was compared to the distribution in the corresponding plane (AX, SAG, COR) extracted from predicted dose distribution by 3DVH and PINNACLE using a 2D gamma analysis. Gamma analysis was performed with 2% dose, 2 mm DTA, 20% threshold, and global normalization. In addition, the percent difference between 3DVH and PINNACLE ROI volumes was calculated. Results: The average ratio of the measured point dose vs the 3DVH predicted dose was 1.017 (σ=0.011). The average gamma passing rate for measured vs 3DVH dose distributions was 95.1% (σ=2.53%). The average percent difference of 3DVH vs PINNACLE ROI volume was 2.29% (σ=2.5%). Conclusion: The dose distributions predicted by 3DVH using MAPCHECK2 measurements are the same as the distributions that would have been obtained using film and chamber. The ROI volumes used in 3DVH are not an exact match to those in PINNACLE; the effect requires more investigation. The accuracy of error detection by 3DVH is currently being investigated.

60 APPLIED LIFE SCIENCES↗

SU-E-T-142: Automatic Linac Log File: Analysis and Reporting

Purpose: End to end QA for IMRT/VMAT is time consuming. Automated linac log file analysis and recalculation of daily recorded fluence, and hence dose, distribution bring this closer. Methods: Matlab (R2014b, Mathworks) software was written to read in and analyse IMRT/VMAT trajectory log files (TrueBeam 1.5, Varian Medical Systems) overnight, and are archived on a backed-up network drive (figure). A summary report (PDF) is sent by email to the duty linac physicist. A structured summary report (PDF) for each patient is automatically updated for embedding into the R&V system (Mosaiq 2.5, Elekta AG). The report contains cross-referenced hyperlinks to ease navigation between treatment fractions. Gamma analysis can be performed on planned (DICOM RTPlan) and treated (trajectory log) fluence distributions. Trajectory log files can be converted into RTPlan files for dose distribution calculation (Eclipse, AAA10.0.28, VMS). Results: All leaf positions are within +/−0.10mm: 57% within +/−0.01mm; 89% within 0.05mm. Mean leaf position deviation is 0.02mm. Gantry angle variations lie in the range −0.1 to 0.3 degrees, mean 0.04 degrees. Fluence verification shows excellent agreement between planned and treated fluence. Agreement between planned and treated dose distribution, the derived from log files, is very good. Conclusion: Automated log file analysis is a valuable tool for the busy physicist, enabling potential treated fluence distribution errors to be quickly identified. In the near future we will correlate trajectory log analysis with routine IMRT/VMAT QA analysis. This has the potential to reduce, but not eliminate, the QA workload.

60 APPLIED LIFE SCIENCES↗