Report on G4-Med, a Geant4 benchmarking system for medical physics applications developed by the Geant4 Medical Simulation Benchmarking Group.

BACKGROUND: Geant4 is a Monte Carlo code extensively used in medical physics for a wide range of applications, such as dosimetry, micro- and nanodosimetry, imaging, radiation protection, and nuclear medicine. Geant4 is continuously evolving, so it is crucial to have a system that benchmarks this Monte Carlo code for medical physics against reference data and to perform regression testing. AIMS: To respond to these needs, we developed G4-Med, a benchmarking and regression testing system of Geant4 for medical physics. MATERIALS AND METHODS: G4-Med currently includes 18 tests. They range from the benchmarking of fundamental physics quantities to the testing of Monte Carlo simulation setups typical of medical physics applications. Both electromagnetic and hadronic physics processes and models within the prebuilt Geant4 physics lists are tested. The tests included in G4-Med are executed on the CERN computing infrastructure via the use of the geant-val web application, developed at CERN for Geant4 testing. The physical observables can be compared to reference data for benchmarking and to results of previous Geant4 versions for regression testing purposes. RESULTS: This paper describes the tests included in G4-Med and shows the results derived from the benchmarking of Geant4 10.5 against reference data. DISCUSSION: Our results indicate that the Geant4 electromagnetic physics constructor G4EmStandardPhysics_option4 gives a good agreement with the reference data for all the tests. The QGSP_BIC_HP physics list provided an overall adequate description of the physics involved in hadron therapy, including proton and carbon ion therapy. New tests should be included in the next stage of the project to extend the benchmarking to other physical quantities and application scenarios of interest for medical physics. CONCLUSION: The results presented and discussed in this paper will aid users in tailoring physics lists to their particular application.

Dosimetric evaluation of nuclear interaction models in the Geant4 Monte Carlo simulation toolkit for carbon-ion radiotherapy.

We tested the ability of two separate nuclear reaction models, the binary cascade and JQMD (Jaeri version of Quantum Molecular Dynamics), to predict the dose distribution in carbon-ion radiotherapy. This was done by use of a realistic simulation of the experimental irradiation of a water target. Comparison with measurement shows that the binary cascade model does a good job reproducing the spread-out Bragg peak in depth-dose distributions in water irradiated with a 290 MeV/u (per nucleon) beam. However, it significantly overestimates the peak dose for a 400 MeV/u beam. JQMD underestimates the overall dose because of a tendency to break a nucleus into lower-Z fragments than does the binary cascade model. As far as shape of the dose distribution is concerned, JQMD shows fairly good agreement with measurement for both beam energies of 290 and 400 MeV/u, which favors JQMD over the binary cascade model for the calculation of the relative dose distribution in treatment planning.

Solar 8B and hep neutrino measurements from 1258 days of Super-Kamiokande data.

Solar neutrino measurements from 1258 days of data from the Super-Kamiokande detector are presented. The measurements are based on recoil electrons in the energy range 5.0-20.0 MeV. The measured solar neutrino flux is 2.32+/-0.03(stat)+0.08-0.07(syst)x10(6) cm(-2) x s(-1), which is 45.1+/-0.5(stat)+1.6-1.4(syst)% of that predicted by the BP2000 SSM. The day vs night flux asymmetry (Phi(n)-Phi(d))/Phi(average) is 0.033+/-0.022(stat)+0.013-0.012(syst). The recoil electron energy spectrum is consistent with no spectral distortion. For the hep neutrino flux, we set a 90% C.L. upper limit of 40x10(3) cm(-2) x s(-1), which is 4.3 times the BP2000 SSM prediction.

Constraints on neutrino oscillations using 1258 days of Super-Kamiokande solar neutrino data.

We report the result of a search for neutrino oscillations using precise measurements of the recoil electron energy spectrum and zenith angle variations of the solar neutrino flux from 1258 days of neutrino-electron scattering data in Super-Kamiokande. The absence of significant zenith angle variation and spectrum distortion places strong constraints on neutrino mixing and mass difference in a flux-independent way. Using the Super-Kamiokande flux measurement in addition, two allowed regions at large mixing are found.

Tau neutrinos favored over sterile neutrinos in atmospheric muon neutrino oscillations.

The previously published atmospheric neutrino data did not distinguish whether muon neutrinos were oscillating into tau neutrinos or sterile neutrinos, as both hypotheses fit the data. Using data recorded in 1100 live days of the Super-Kamiokande detector, we use three complementary data samples to study the difference in zenith angle distribution due to neutral currents and matter effects. We find no evidence favoring sterile neutrinos, and reject the hypothesis at the 99% confidence level. On the other hand, we find that oscillation between muon and tau neutrinos suffices to explain all the results in hand.