14th Geant4 Users and Collaboration Workshop

15 Oct 2009, 08:30 → 22 Oct 2009, 23:30 Europe/Rome

Laboratori Nazionali del Sud - INFN

This is the timetable of the 2009 Workshop of the Geant4 Collabration. This year the Workshop will be constituted by a first Users Session (15-17 October), where participation of Geant4 users is welcome and by a Collaboration Session (19-22 October) where only Geant4 Members areare admitted. Abstract submission is admitted for the Users Session (15-17 October). Both Geant4 Members and external Users can send their abstracts For Workshop registration, please go to the following page: http://www.shougun.it/CongressiScheda.aspx?IdCongresso=20&lingua=en

more information WorkshopLogisticInfo.pdf

Thursday, 15 October

Participants registration

Opening

slides Workshop_Logistic_Info.pdf Workshop_Logistic_Info.ppt

1 Workshop Logistics

Speaker: Pablo Cirrone (Unknown)

Slides Workshop_Logistic_Info.pdf Workshop_Logistic_Info.ppt

2 Workshop Goals and Opportunities

Slides G4Catania2009-Goals-and-Issues-JApostolakis.pdf G4Catania2009-Goals-and-Issues-JApostolakis.pptx

Plenary Session I - HEP and Hadronics

3 Geant4 applications in high-energy physics

Speaker: Margar Simonyan (lapp-Laboratoire d'Annecy-le-Vieux de Physique des Particules-In)

Slides main.pdf

4 Geant4 Hadronic Physics Performance: Recent validation and developments

Speaker: Jose Manuel Quesada Molina (Depto. de Fis. Atom., Mole. y Nucl.-Universidad de Sevilla-Unkno)

Slides HadOverview_JMQ.pdf HadOverview_JMQ.pptx

5 Geant4 applications for astroparticle and rare-event physics

Speaker: Mr Markus Horn (Ph.D.student)

Slides Horn_G4workshop_Oct09.pdf

12:45 Short break

Lunch

Parallel Session I - Hadronic Physics Validation

6 Nuclear beams in extended media: issues of nuclear fragmentation

We study the validity of nucleus-nucleus collision models of Geant4 describing fragmentation of light and medium-weight nuclei (from $^{3}$He to $^{58}$Ni) in tissue-like media. Depth-dose distributions and yields of secondary fragments produced by nuclear beams were calculated and compared against experimental data. The Light-Ion Binary Cascade model and the Wilson abrasion model were used to simulate the initial stage of nucleus-nucleus collisions. The both models were successful in simulating the transport of light nuclei, e.g. $^{12}$C, while the cascade model overpredicts the dose from heavy energetic nuclei, possibly due to underpredicting the attenuation of the initial beam due to fragmentation. The calculations were performed with and without the Fermi break-up model and the Statistical Multifragmentation model for multifragment decays of excited light and medium-weight residual nuclei at the final stage of collisions in addition to evaporation model. This allowed us to estimate the role of multifragment decays on the depth-dose, depth-yield and charge distributions of secondary fragments produced in extended media. While the energy deposition is not very sensitive to decay mechanisms of excited nuclei, the fragment yields can not be properly described without accounting for multifragment decays. The validity of the Statistical Multifragmentation model of Geant4 has been confirmed in a set of stand-alone tests, following several fixes introduced to G4StatMF. These corrections improved the stability of the model and consistency of simulation results with the FORTRAN implementation of the model. Comparisons to experimental data on multifragmentation of hot nuclei in central collisions of nuclei at low energies are also presented.

Speaker: Dr Igor Pshenichnov (Frankfurt Institite for Advanced Studies (FIAS) and Institute for Nuclear Research (INR), Moscow)

Slides Pshenichnov_G4Catania.pdf

7 Test of Hadronic Models in GEANT4 using the BESIII Data

The BESIII detector is a conventional solenoidal magnetic spectrometer and operating at the upgraded Beijing Electron-Positron Collider (BEPCII). Since the start of commissioning run in Aug. 2008, about 100M psi(2S) data and 200M J/Psi data have been collected with BESIII detector. Using 10M psi(2S) data taken in 2008, we compare hadronic shower energy, profile and fake photons in the electromagnetic calorimeter with several hadronic models in GEANT4. We find that for pions and protons at low energy (<0.7 GeV) the Bertini model with high precision neutron tracking (QGSP_BERT_HP) is the best among the hadronic models compared, though it takes more simulation time. For anti-proton, the elastic cross section seems to be suppressed and needs to be improved. We also try to replace GHEISHA by CHIPS model in QGSP_BERT_HP for anti-proton inelastic, marked with CHIPS+HP. Although the comparison result of CHIPS+HP is slightly better than other models, it is also not reasonable.

Speakers: Dr Guofu Cao (Institute of High Energy Physics), Prof. Huaimin Liu (Institute of High Energy Physics)

Slides Hadronic_models_caogf.pdf Hadronic_models_caogf.ppt

8 Beam Test and Simulation Results with Highly Granular Calorimeters for the ILC

To evaluate technologies for ILC calorimetry, the CALICE collaboration has constructed a prototypes of highly granular sampling calorimeters. These detectors have been tested extensively in particle beams at DESY, at CERN and at Fermilab. The imaging capabilities of these detector provide three dimensional information of hadronic showers with unprecedented resolution and will thus help to constrain hadronic shower models in simulations. We present results from the analysis of hadronic events including studies of the longitudinal and transverse shower profiles. The results are compared to simulations with a variety of different models. We put particular emphasis on the comparison of our data to the new physics lists proposed by the G4-hadronic team.

Speaker: Mr Benjamin Lutz (Deutsches Elektronen-Synchrotron (DESY))

Slides BLutzCALICE.pdf.pdf

9 slic & lcdd: A Detector Response Simulation Program

Speaker: Norman Graf (SLAC)

Slides

• slic.pdf
• slic.ppt
• slic.ppt

Parallel Session II - EM Physics: Validation and Applications

10 Validation of Geant4 EM physics for gamma rays against the SANDIA, EPDL97 and NIST databases

From version 9.2 of Geant4, the LowEnergy electromagnetic processes (Livermore and Penelope) have been migrated to the design introduced for the Standard EM models. In the new approach there is only one process and multiple models that can be registered to the process. We present a validation of the migrated Geant4 electromagnetic photon models for elements and compounds with respect to several libraries based on experimental data (SANDIA, EPDL97 and NIST). The cross section of all the photon models agrees with NIST within 10%.

Speaker: Zhang Qiwei (INFN-LNS)

Slides 14th_Geant4_Workshop-Zhang.pdf 14th_Geant4_Workshop-Zhang.ppt

11 Evaluation of absorbed fractions for beta-gamma radionuclides in ellipsoidal volumes of soft tissue through Geant4

We developed a simulation in Geant4 to calculate the absorbed fractions for monoenergetic and beta electrons emitted by 199Au, 177Lu, 131I, 153Sm, 186Re and 90Y and for photons between 10 keV and 1 MeV, emitted by sources uniformly distributed in ellipsoidal volumes of soft tissue. Code validation results with respect to reference data for doses, ranges and absorbed fractions in spheres are presented. We discuss the influence Monte Carlo parameters on average energy deposition for the three physics packages available. An analytical relationship between absorbed fraction and a “generalized radius” is introduced, and the dependence of its parameters rho0 and s from average electron or photon energy is discussed. A generalization for the estimation of absorbed fractions for other radionuclides is also proposed. Such results can be useful to improve accuracy and easiness of calculation in dosimetry of beta-gamma emitting radionuclides.

Speaker: Dr Ernesto Amato (University of Messina)

Slides Amato_G4_LNS2009.pdf Amato_G4_LNS2009.ppt

12 Complete Monte Carlo simulation of the optical response of the WArP detector for direct Dark Matter search

WArP is a double phase Dark Matter experiment using Liquid Argon as a target. In order to efficiently collect the 128 nm scintillation light a very careful light collection system, made by a wavelength shifter and a dielectric mirror, is used. A detailed & complete GEANT4 simulation of the detector has been performed, from particle interaction to PMTs photon detection. Optical photons tracking as well as wavelength shifting and boundary optical processes on the light collection system surface have been implemented. A dedicated C++ code, simulating the electronics and DAQ systems, has been also used to obtain WArP like data starting from the PMTs photon hits.

Speaker: Dr Francesco Di Pompeo (Gran Sasso National Laboratory)

Slides dipompeo.pdf

Coffe/Tea break

Parallel Session III - Medical

13 From imaging to dosimetry: Geant4-based study on the application of Medipix2 to neutron personnel dosimetry.

A Geant4 simulation study is performed to optimise a novel neutron personnel dosimeter, developed at the Centre of Medical Radiation Physics (CMRP), University of Wollongong, consisting of a silicon pixelated detector device - Medipix2- covered with a polyethylene converter on top. In this talk we would like to illustrate the Geant4 simulation developed to optimise the novel neutron dosimeter and its validation. First results of the optimization study of Medipix2 polyethylene converter will be shown as well.

Speaker: Dr S. Guatelli (University of Wollongong, Australia)

Slides guatelli_neutron_dosimetry.pdf guatelli_neutron_dosimetry.ppt

14 Skin Dosimetry in MRI-guided Radiotherapy (MRIgRT)

Geant4 has been used to simulate in high-resolution (10 micron thick voxels) the skin doses below 0.5 mm depth occurring in MRIgRT (6 MV photon beam, Varian 2100C). On the entry side lepton contamination removal has been studied to show potential skin dose reductions (at low magnetic field strengths). At high magnetic fields (>2 T) skin dose increases (up to 30% of dmax). In the exit region the Electron Return Effect (ERE) gives rise to massive skin dose increases (>100%). The use of exit bolus is investigated for reducing these increases.

Speaker: Mr Bradley Oborn (Centre for Medical Radiation Physics)

Slides Oborn_Skin_Dosimetry_in_MRIgRT.pdf

15 Geant4-based studies to characterise the tissue-equivalence of SOI and diamond microdosimeter detectors, under development at CMRP.

Silicon On Insulator (SOI) microdosimeters have been under investigation for the past ten years as a possible alternative to tissue equivalent gas counters for microdosimetric measurements in medical physics and radiation protection in earth labs, aviation and space, at Centre of Medical Radiation Physics, University of Wollongong. In the last 5-6 years, extensive Geant4-based simulation studies have been performed to characterise the novel SOI detector. This talk is meant (1) to be a review of the past Geant4-related activities in silicon microdosimetry and (2) to show current Geant4-studies in microdosimetry at CMRP.

Speaker: Dr S. Guatelli (CMRP, University of Wollongong, Australia)

Slides Guatelli_silicon_microdosimetry_G4Workshop.pdf Guatelli_silicon_microdosimetry_G4Workshop.ppt

16 Dosimetric study of photon dose distribution in lungs under different respiratory phases: comparison with GEANT4 simulations

The advances in oncological research have had an important impact on cancer patients survival, but there has been no major or dramatic improvement in long-term survival in the lung cancer despite more innovative treatment techniques and protocols implemented in radiation oncology. The problematic concerning dose calculation in the presence of in-homogeneities has already been afforded by several authors but a new challenge has been rising since radiation therapies based on respiratory gating systems have been introduced in the clinical practice. The MonteCarlo (MC) algorithms have already been proven to be a realistic alternative to analytical algorithms in those regions where charged particle equilibrium does not hold. The aim of this work is to evaluate the capability of GEANT4 respect to some commercial treatment planning systems to provide dose calculation maps with high level of accuracy also when lung densities are changing (i.e. during different respiratory phases, from breathing out to deep inspiration). GEANT4 simulations have been compared to dosimetric data measured by radiochromic films. In this study, a 6 MV beam from a 2100 Clinac Varian was modeled with a GEANT4 MC code. The reported results were obtained with a modified version of the GEANT4 Advanced Example MedLinac. Preliminary measurements were carried out for different square field sizes (20x20 cm2, 7.5x7.5 cm2 and 3x3 cm2 ) with SSD=84 cm. To study the dose distribution a series of three phantoms were constructed and simulated, each characterized by a different lung density according to the following values: 0.030, 0.080 and 0.40 g/cm3. The MC simulation jobs were run on a dedicated Beowulf cluster located in the Technology and Health Department of the Italian National Institute of Health. Two different algorithms released for clinical application from two commercial TPS were included in our study: Varian Eclipse TPS (PB-EqTAR pencil beam with EqTAR algorithm for heterogeneities correction) and Philips Pinnacle TPS (CCC Collapsed Cone Convolution). Results obtained by the MC simulation and calculated data obtained from commercial TPS have been compared with measurements performed with radiochromic films (GafchromicTM EBT). The films were exposed at different depths perpendicularly and parallel to the beam axis at different distances. Images were acquired with a Epson Expression 10000 XL flatbed scanner and analyzed with Picodose X PRO software. Results show a good agreement between GEANT4 simulations and dosimetric data.

Speaker: Dr Barbara Caccia (Technologies and Health Dept. Istituto Superiore di Sanità and INFN, Roma (Italy))

Slides Andenna_Caccia_Dosimetric_study_of_photon_dose.pdf Andenna_Caccia_Dosimetric_study_of_photon_dose.ppt

17 Implementation of a Monte Carlo - GEANT4 Simulation for the dosimetric study of electron beams produced by a mobile accelerator for IORT

In the last few years the use of ionising radiation has obtained a relevant role as a remedy for many types of tumour pathologies. Among the different radiation techniques the Intra Operative radiation therapy (IORT) is a new way to treat some tumours [1,2,3,4], consisting of exposing the surrounding involved tissues to a single high radiation dose during surgery after the obliteration of the tumour [5,6,7]. Many studies have highlighted the success of this approach, especially to treat breast lesions. The employment of conventional linear accelerators (linac) needs complex procedures and offers several limitations making very difficult its adoption in terms of time and money [8]. In the last years, a new generation of linac has been developed specifically for the IORT, characterized by limited energy, limited field size, no bending magnet and high mobility. The design of these linacs is optimized for minimal radiation leakage allowing them to be operated in an unmodified surgery room. Our study is focused on one of these linacs, the NOVAC7 (New Radiant Technology S.p.A.). The aim is to study the dose distribution inside a water phantom changing the dimension and shape of collimators, and the dimension, shape and material of beam stopper generally localized beyond the target in the breast lesion treatments. The first step of our study is to develop a Monte Carlo application using GEANT4 toolkit to simulate the complex geometry of NOVAC7 and the electron beam characteristics, and to compare the dose distribution in the standard dosimetrical configuration with the experimental data, measured using the linac installed at the LATO HSR Giglio (LAboratorio di Tecnologie Oncologiche) at Cefalù (PA). [1] Garton G R, Gunderson L L, Webb M J, Wilton T O, Cha S S and Podratz K C, Intraoperative radiation therapy in gynecologic cancer: update of the experience at a single institution, Int. J. Radiat. Oncol. Biol. Phys. 37 893-43 1997 [2] Krempien et al Long-term results of intraoperative presacral electron boost radiotherapy (IOERT) in combination with total mesorectal excision (TME)and chemoradiation in patient with locally advanced rectal cancer Int. J. Radiat. Oncol. Biol. Phys. 66 1143-51 2006 [3] Roeder F et al. Patterns of failure and local control after intraoperative electron boost radiotherapy to the presacral space in combination with the total mesorectal excision in patients with locally advanced rectal cancer Int. J. Radiat. Oncol. Biol. Phys. 67 1381-8 2007 [4] Tran et al Long-term survivors using intraoperative radiotherapy for recurrent gynecologic malignancies Int. J. Radiat. Oncol. Biol. Phys. 69 504-11 2007 [5] R.Orecchia et al. Intraoperative electron beam radiotherapy (ELIOT)to the breast: a need for quality assurance programme , The Breast 14, 541-546 2005 [6] Merrick III HW et al. Future directions in radiation therapy Surg. Oncol. Clin. N. Am. 12, 1099-105 2003 [7] Willet CG et al. Intraoperative radiation therapy Int. J. Clin. Oncol. 6 209-14 2001 [8] Rapporti ISTISAN 03/1 IT

Speaker: Dr Giorgio Russo (Laboratorio di Tecnologie Oncologiche (LATO) HSR Giglio, Cefalù (PA), Italy)

Slides RussoGiorgioIORT-MC-Geant4.pdf RussoGiorgioIORT-MC-Geant4.ppt

18 The role of Geant4 in the production of a database for an ion therapy Treatment Planning System

In radiotherapy with carbon ions, the beam interaction with tissues has to be carefully estimated, since the spatial changes in the radiation quality due to fragmentation, energy loss and multiple scattering imply variations in the local radiobiological effectiveness. The information of dose deposition pattern alone is not sufficient to predict the biological effect of treatment: a detailed three-dimensional description of beam composition in terms of ion species, kinetic energies and energy releases is needed to estimate the outcome on tumoral and normal cells. In the newly Italian Treatment Planning System (TPS) project, funded by the Istituto Nazionale di Fisica Nucleare (INFN) and the Ion Beam Accelerators (IBA), this knowledge is gathered by simulating the beam propagation by means of both Fluka and Geant4 codes. In this talk will be presented the way this task is accomplished using Geant4, pointing out some apparent irregularities we are finding in the results.

Speaker: Germano Russo (INFN (Istituto Nazionale di Fisica Nucleare), Dipartimento di Fisica Sperimentale - Università di Torino, Italy)

Slides RussoGermanopresentation_G4Workshop.pdf

19 Geant4 based simulation of the Leksell Gamma Knife for treatment planning validation

The Monte Carlo code Geant4 has been used to simulate the Leksell Gamma Knife® and to verify its treatment planning system Leksell GammaPlan® (LGP). By means of this radiosurgical technique, intracranial lesions can be treated in a single session with high precision and critical brain structures can be protected. Radiation from 201 60Co sources comes through a collimation system to the target area, focusing on the isocenter, where the maximum dose amount is released. Simulations were employed to calculate the dose distribution in a spherical water phantom, the same we used for the validation measurements. In order to simplify the code, we simulated just one single source, rotating the phantom at 201 angular positions, the same of the sources in the device. We validated the simulation by measuring the dose distribution with radiochromic films HS, for all the available collimators (4, 8, 14, 18 mm). Good agreement has been found between experimental data and results of Monte Carlo simulations. Then we used the simulation to verify the dose calculated by LGP. Comparisons between Geant4 and LGP outputs in one and two dimensions were performed at different planes (axial, coronal and sagittal) and depths. The outcomes show that LGP calculates correctly the dose distribution in a homogeneous phantom. Different results were found when different density materials are also included, because LGP assumes the phantom only composed by water. Moreover a complete clinical treatment with a multi-shots configuration has also been simulated and also in this case the results were compared with LGP. Some methods for the optimization of the code and the reduction of calculation time are in progress, in order to possibly propose the application as an advanced example of the Geant4 toolkit.

Speaker: Dr Francesco Romano (Istituto Nazionale di Fisica Nucleare, Laboratori Nazionali del Sud, Catania, Italy)

Slides RomanoF_GammaKnife.pdf

Parallel Session IV - Underground

20 Muon Spallation Products in KamLAND

The study of the muon spallation products is essential for rate event detection in neutrino detectors, double beta decay experiments, and dark matter searches. In the KamLAND detector, energetic muons interact mainly with 12C in the liquid scintillator and generate neutrons and isotopes by electromagnetic or hadronic processes. The neutron production yield is evaluated to be (2.8 +- 0.3) x 10^-4 n / (mu (g/cm2)). For other isotopes, the production yields are measured from the observed time correlation related to known isotope lifetimes.

Speaker: Dr Itaru Shimizu (RCNS)

Slides G4Workshop-Oct2009-Shimizu.pdf

21 MonteCarlo Simulations with GEANT4 for the XENON100 Detector

The XENON100 detector is a dual phase (liquid-gas) xenon time-projection chamber for direct dark matter detection. The total amount of liquid xenon is 165 kg, of which 65 kg are in the active target enclosed in a teflon/copper structure, the rest being in the surrounding active veto volume. UV light signals are produced by particle interactions and detected by 242 PMTs (178 in the target and 64 in the veto). In this talk, the importance of MonteCarlo simulations for the complete understanding of the detector performance and the expected sensitivity of the experiment will be discussed. We will report simulation results and comparison with the experimental data, including light simulation and background predictions.

Speaker: Mr Alexander Kish (University of Zürich)

Slides Kish_XENON100.pdf

22 GEANT4 simulation of the Borexino solar neutrino experiment.

The BOREXINO experiment measures the MeV and subMeV Solar neutrinos with the Borexino low background liquid scintillator detector at Gran Sasso National Laboratory. A GEANT4 based simulation was developed and used to model the physical characteristics of the Borexino detector. Object oriented structure of GEANT4 is very suitable to describe the full detector geometry and different physical processes in the Borexino detector . The code is used for debugging and inspection of energy and spatial reconstruction algorithms, simulation of neutrino and background effects , energy scale calibration. The code structure , it's performance and the comparison of the Monte-Carlo simulation results with the experimental data from the operating detector are presented.

Speaker: Igor Machulin (Kurchatov Institute, Moscow)

Slides Geant4-Borex_-PPT97-2003.pdf Geant4-Borex_-PPT97-2003.ppt

Friday, 16 October

Plenary Session II - Space and Physics

23 Applications of Geant4 in space

Speaker: Dr Giovanni Santin (ESA)

Slides Geant409_Catania_Space_ESA_v15.pdf

24 Geant4 EM Physics Performance: Recent developments and validation

Speaker: Dr Ziad Francis

Slides G4W-LowE_VFINAL.pdf G4W-LowE_VFINAL.ppt

25 Physics performance: topical discussion

Speaker: Dr Alberto Ribon (CERN)

Slides phys_perf_catania09.pdf

Cofee/Tea Break

Plenary Session III - Medical and Physics Lists

26 Applications of Geant4 in Proton Radiotherapy at the University of Texas M.D. Anderson Cancer Center

Speaker: Prof. Jeremy C. POLF (UTMD Anderson Cancer Center, Houston, Texas, US)

Slides Polf-GEANT4_Conference_Presentation.pdf Polf-GEANT4_Conference_Presentation.ppt

27 Physics Lists

Lunch

Parallel Session V - HEP

slides g4studiesdualreadout2009_genser.pdf g4studiesdualreadout2009_genser.pptx

28 The Geant4 Simulation of the ATLAS Experiment

Speaker: Makoto Asai (SLAC)

Slides ATLAS.pdf ATLAS.pptx

29 ALICE Experience with Geant3, Fluka and GEANT4

Speaker: Eva Sicking (Institut fuer Kernphysik-Westfaelische Wilhelms-Universitaet Mun)

Slides sicking_ALICE.pdf

30 Geant4 Results from CMS

Speaker: Sunanda Banerjee (Fermilab)

Slides cstf38.pdf

31 Geant4 Studies in the Context of Dual Readout Calorimetry Simulations

Speaker: Krzysztof Genser (Fermilab)

Slides g4studiesdualreadout2009_genser.pdf

Parallel Session VI - Medical

Convener: Joseph Perl (SLAC)

32 The new version of the GATE simulation platform

Monte Carlo simulation is an important tool in medical physics to understand and analyse data taken from medical devices. It is also useful to assist the development of new medical devices and optimise their use. GATE, a Geant4 Application, has been developed since 2001, initially for emission tomography (PET and SPECT) simulations. It encapsulates the GEANT4 libraries in order to achieve a modular, versatile, scripted simulation toolkit adapted to the field of nuclear medicine. In particular, GATE provides the capability for modelling time dependent phenomena such as detector movements or source decay kinetics, thus allowing the simulation of time curves under realistic acquisition conditions. It includes well-validated physics models, geometry modelling tools convenient for complex scanner geometries, models for detector electronic. The version 6 of the GATE simulation platform is built on the previous developments and will be released publicly at the end of 2009. The code has been consolidated and the computational efficiency of GATE has been enhanced. The platform has also been extended to new applications such as modelling of CT imaging or radiotherapy and hadrontherapy treatments. For these applications, new tools have been developed, for example, to easily calculate the deposited dose distribution inside a patient described by a CT image following an intensity modulated radiation therapy (IMRT) photon beam, or a proton/carbon ion beam irradiation. The Geant4 hadronic physic processes with flexible models and cross-section selection are included and tested towards complete simulation of IMRT hadrontherapy. GATE v6 should thus become an invaluable simulation tool for assisting in the development of imaging-assisted radiotherapy and on-line dose control strategies.

Speaker: Dr Thibault Frisson (Centre Léon Bérard)

Slides G4Workshop_GATE.pdf

33 Recent developments and features in Hadrontherapy advanced example

Hadrontherapy is an advanced example of Geant4 toolkit. Hadrontherapy application was originally developed to simulate a well specified proton therapy facility: the passive transport beam line installed at Laboratori Nazionali del Sud (INFN) in Catania, Italy. In the recent past, Hadrontherapy went through different modifies. Nowadays, Hadrontherapy is a multifaceted and more flexible application showing many additional capabilities in the ion radiotherapy field. Our idea is to transform Hadrontherapy into a multi-modules application able to retrieve information of interest for proton-ion radiotherapy. The multi-modules scheme is a power tool that permits to use and integrate different independent components. For example, the source module is dedicated to simulate the characteristics of the initial beam. The geometrical module consents to employ different geometrical setups completely interchangeable switching between different geometrical configurations: passive/active proton beam line and a radiobiological carbon beam line. In the actual version only a 'passive beam line set-up' is available but new geometrical setups will be added in the next future. The detector module allows the simulation of a typical radiotherapy phantom divided in voxel. In the physics module are implemented the electromagnetic and hadronic physical model. And, in the data module all the physical and biological information (like dose, fluence, LET, etc) are collected

Speaker: Dr Giuseppe Antonio Pablo Cirrone (INFN - LNS)

Slides Hadrontherapy.pdf

34 PTSsim Framework for Radiotherapy Monte Carlo

PTSsim Framework for Radiotherapy Monte Carlo

Speaker: Dr Takashi Sasaki (KEK)

Slides Takashi.pdf

35 GAMOS Framework for Medical Monte Carlo

GAMOS Framework for Medical Monte Carlo

Speaker: Pedro Arce (CIEMAT)

Slides g4ws09.gamos.pdf g4ws09.gamos.ppt

Cofee/Tea Break

Parallel Session VII - Nuclear and Space Physics

36 GEANT4 simulation project on the AMS facility, ARTEMIS, at LMC14 in Saclay, France

A project of Ion Beam Simulation on the ARTEMIS Accelerator Mass Spectrometer is carried out at LMC14. The tool GEANT4 is used to take into account all physical processes and predict and control the settings of such a facility.

Speaker: Mrs Anita QUILES (CEA SACLAY LMC14)

Slides GEANT4_simulation_project_on_the_ARTEMIS_AMS_facility.pdf

37 Simulation of MEDEA Response to Neutrons

An application of GEANT4 version 9.2p01 to 4π Medea detector will be presented. MEDEA (1) is a working array aimed at the gamma and light charged particle detection in the range of intermediate energies (up to 100 MeV/A), installed at LNS. MEDEA is made of 180 BaF2 detectors, arranged in a spherical shape of 22 cm radius. Each module is made of 2 separated crystals in optical contact. The interest in neutron detection to extend the knowledge of the reaction dynamics has triggered a study of a possible shortening of a few MEDEA modules with the aim of increasing path length for neutron detection, therefore improving their energy calculation through Time Of Flight determination. The MEDEA geometry has been implemented in the GEANT4 framework and the response of the whole system to sources of monochromatic neutrons up to 18 MeV, impinging on a module has been studied. Results and experienced problems will be shown. 1) E.Migneco, Nuclear Instruments and Methods, A314(1992)31

Speaker: Dr concettina maiolino (LNS)

Slides Maiolino.pdf Maiolino.ppt

38 The radiation hardness assurance tests at INFN-LNS Catania

This work presents the features of a monitoring system developed for radiation hardness studies on electronics components and systems at Superconducting Cyclotron at INFN-LNS (Istituto Nazionale di Fisica Nucleare, Laboratori Nazionali del Sud, Catania, Italy). Relevant beam parameters as flux, 3-D beam profile and energy are measured with high accuracy with double sided 1.5 mm thick microstrip silicon detector which act as calorimeter. All monitoring setup is fully automatic to perform Single Event Effects studies with heavy ions. In this study gaseous ions Ne20, Ar40, Kr84 and Xe129 with 20 MeV/n energy are used. The operation in air, which means less complication in test setup and faster execution times provides also the advantage of using the air thickness in front of Device Under Test (DUT) as degrader hence to fine tune the energy values available at the DUT surface. A careful evaluation of energy loss in air and of the energy spread at DUT surface is carried out through a full Monte Carlo simulation of test setup and comparing the results with data. The results of careful investigations for beam flux uniformity in space and time as well as effects of fragmentation which may happen in air are also discussed. We have exposed during the same session, two beam calibration systems, the "Reference SEU monitor" developed by ESA/ESTEC and the beam monitoring and dosimetry setup developed by our group. The results are compared and discussed here.

Speaker: Francesca Renzi (MAPrad s.r.l)

Slides Hardness_assurance_tests_at_INFN-LNS_Catania.pdf Hardness_assurance_tests_at_INFN-LNS_Catania.ppt

Parallel Session VIII - Medical

39 Geant4 simulation of the attenuation properties of plastic shields for beta radionuclides employed in internal radiotherapy

We developed a simulation in Geant4 to compare the attenuation properties and the bremsstrahlung radiation yield of different types of plastic materials employed as shield of beta radioactive sources. Code validation results against Sandia and NIST data are presented. We discuss the influence of cuts in range, step limits and multiple scattering step limitations, for the three physics packages available. For polypropylene, polystyrene, polyamide nylon-6, poly-methyl methacrylate, polycarbonate, polyethylene terephthalate, polyvinyl chloride and polytetrafluoroethylene we evaluated the mean and maximum ranges for electrons originating from 90Sr and 90Y, as well as the number and the spectrum of the bremsstrahlung X-rays produced. Significant differences appear between the various materials, and the choice of the best one depends also on the physical properties requested for each specific application.

Speaker: Dr Domenico Lizio (University of Messina)

Slides geant4_lns_schermi-beta_Definit.pdf geant4_lns_schermi-beta_Definit.ppt

40 USE OF GEANT4 CODE FOR VALIDATION OF RADIOBIOLOGICAL PARAMETERS OBTAINED AFTER PROTON AND CARBON IRRADIATIONS OF MELANOMA CELLS

Series of numerical simulations with the GEANT4 code were carried out to verify and better interpret radiobiological parameters and biological effects obtained after irradiation of HTB140 human melanoma cells. Irradiations were performed at INFN-LNS with the 62 MeV proton beams at the full energy Bragg peak as well as along the spread out Bragg peak (SOBP) and its distal fall-off part. HTB140 cells were also exposed to the beams of the 62 MeV/u 12C ions, along the Bragg curve including its distal declining edge. The doses ranged from 2 to 24 Gy. The surviving fractions at 2 Gy (SF2) and the corresponding relative biological effectiveness (RBE) obtained with respect to gamma-rays, differed along the proton full energy and the spread out Bragg peak (SOBP), as well as along the carbon Bragg peak. These variations were analysed and compared through the proton and carbon depth dependant dose, fluence, energy and linear energy transfer (LET) distributions produced by numerical simulations. For carbon ions these distributions were also followed for the secondary particles. Depth dependent distributions of fluence and LET revealed different damages in the sense of the ratio of irreparable to reparable breaks that were caused by distinct radiation types. This enabled a more in-depth interpretation of cell proliferation results witch reflect the quality of life of irradiated cells and was particularly useful for the assessment of the events at the distal declining edge of the Bragg curves.

Speaker: Prof. Ivan Petrovic (Vinca Institute of Nuclear Sciences, Belgrade, Serbia)

Slides Presentation_G4_Workshop__Petrovic1.pdf Presentation_G4_Workshop__Petrovic1.ppt

41 Physical and Biological investigations using Geant4 Monte Carlo simulation of the beam delivery line components in particle therapy.

Background: The importance of Monte Carlo (MC) simulations in the field of advanced radiation therapy increases with the necessity to investigate the effect of detectors, monitoring systems and ripple filters at the design phase. In this work we present our investigations using MC simulations of the effects of each of the elements in a beam delivery line of a particle therapy facility. Methods and Materials: A full beam delivery line of the national center of oncologic hadrontherapy (CNAO) is simulated with the Monte Carlo package Geant4 to get the actual distribution of particles and fragments and the corresponding energies, in the treated volume. The evaluation of biological effects was studied using a code based on the Local Effect Model (LEM). The monitoring system on the beam delivery line was fully simulated, as well as the ripple filters and range shifters. The computational effort was performed using the distributed INFN Grid computing resources. Measurements were done within the facilities of INFN Laboratori Nazionali del Sud to compare physical data of the lateral distribution of the beam and the Bragg peak curve (energy loss in depth curve) to MC simulations. The measurements were performed to test the design of the ripple filters to be used at CNAO. We calculated also the corresponding biological effect to evaluate the biological equivalent dose. Results: A good agreement between simulations and experimental measurements for protons was obtained. The test for the design of new filters was successfully realized. Given this good agreement we are confident about other improvements of the nozzle based on the simulation results. We modeled the effect of both ripple filters and monitoring system defining the corresponding transfer functions optimized the design of the monitoring system and especially for the design and test of ripple filters improving some of their features. We provided as well a tool to estimate both the physical and the biological equivalent dose distribution with the definition of corresponding transfer functions. Conclusion: MC simulations using Geant4 helped for the design of the monitoring system and especially for the design and test of a new ripple filters to be installed at CNAO, improving some their features. We provided as well a tool to deduce both the physical and biological equivalent distribution determining their corresponding transfer functions. Physical dose simulated data was largely verified. We are expecting the verification of the biological equivalent dose from future biological measurements using ripple filters.

Speaker: Faiza Bourhaleb (University of Turin, Dpt experimental physics and INFN, Torino, IT)

Slides FaizaBourhalebPhysicalandBioG4Catania.pdf

42 Simulation studies of a therapeutic proton beam delivery system

We have been developing a Monte Carlo simulation software to be applied in clinical usage. To apply a Monte Carlo simulation in clinical support, the reproducibility of real irradiation is prerequisite. From modeling the therapeutic nozzle to mimicking treatment scheme, we have successfully modeled our beam delivery system which is using rotation modulation wheels to generate a spread-out Bragg peaks(SOBP). The initial beam energy for Monte Carlo simulation was estimated based on the Bragg peak comparisons with measurement and used to simulate SOBPs resulting good agreement with measured SOBPs.In this presentation, we will summarize our status of development and introduce user studies applying the developed software.

Speaker: Mr Jungwook Shin (National Cancer Center, Korea)

Slides med88.pdf.pdf

43 New physics processes dedicated to nanometric scale track structure in the Geant4 toolkit

Track structure simulation using the Monte-Carlo method has become increasingly important not only on the macroscopic scale but also on the nanometric and the molecular simulation level. Many studies (mostly done by Nikjoo et al., and Goodhead et al.) showed that damage clustering in cell nucleus is related to the energy deposits distribution inside the cell. However, experimentally it is impossible to guess the exact position of each energy deposit point inside irradiated cells. Monte-Carlo simulations can reproduce the ionizing particle track taking into account all the interactions according to calculated cross sections. Many nano-level codes were developed so far, but generally these codes are dedicated for private use and for a very specific application (Partrac, PITS, OREC…) which makes their access practically difficult and sometimes even impossible due to lack of information, user-support or a simple users manual. The Geant4-DNA project comes as a free open source set of physics processes that can handle nanometric level simulations for electrons (8.22 eV – 1 MeV), protons (100 eV – 100 MeV) and alpha particles (1 keV – 10 MeV). Users can benefit from dedicated documentation, forums and necessary amount of support for new users through workshops and tutorials. This study shows the recent developments in Geant4 concerning the inelastic and the elastic interactions of these particles in liquid water, as water is considered to be the dominant component in a living cell. Models of inelastic scattering cross sections are based on the First Born Approximation (FBA) theory. For low incident energies the FBA becomes inapplicable, correction terms are then used for electrons and semi empirical models (Rudd model for ionisation) are used for Alpha particles and protons. Charge transfer process is also taken into account according to Dingfelder’s model. Relativistic dynamics were recently added for high energy (20 keV – 1 MeV) electrons cross sections. Validation results are still in progress and are not shown in this work. However, a comparison shows a good agreement between microdosimetric calculations (proximity functions and clustering results) obtained using Geant4-DNA and results obtained by other nano-level codes (Nikjoo et al., Chen and Kellerer).

Speaker: Dr Ziad Francis (IRSN FRANCE)

Slides G4DNA_WorkShop_2009_pres1.pdf G4DNA_WorkShop_2009_pres1.ppt

44 INVESTIGATING THE ACCURACY OF HADRONIC MODELS FOR CARBON ION THERAPY

Introduction: An accurate prediction of radiation fields resulting from therapeutic carbon ions in the patient is needed in order to evaluate correctly their biological effectiveness. Materials and Methods: Hadronic models of the Monte Carlo transport codes: Geant4 and Fluka, are benchmarked for regimes relevant for therapeutic carbon ions. The ability of the Monte Carlo codes to reproduce measured fluencies of secondary fragments for a thick-target experiment at differing depths in water is evaluated by simulating the detailed experimental set-up and time-of-flight measurement techniques. Integral and double-differential fluencies are compared. The Geant4 models: BIC, QMD, Fermi Break-up, and Multifragmentation are tested. Conclusions: Agreement within approximate limits of 50% is found for both codes for integral fragment fluencies. Discrepancies between simulations and measurements are specially for forward-directed fragments (0-2 degree). Geant4 tends to underestimate small-angle fluencies whereas larger angles tend to be overestimated.

Speaker: Till Tobias Boehlen (CERN/KI/SU)

Slides BoehlenT_091024_G4UsersWS.pdf BoehlenT_091024_G4UsersWS.ppt

Saturday, 17 October

Plenary Session IV - News from Geant4 and Computing Performance

45 News from Geant4: Recent Developments (non-physics)

Speaker: Ivana Hrivnacova

Slides recent_developments.pdf

46 News from Geant4: Best Practices

Speaker: Joseph Perl (Unknown)

Slides G4BestPractices.pdf G4BestPractices.ppt

47 News from Geant4: Demo of Visualization Capabilities

Speaker: John Allison

Slides New_features_of_Vis_2009.pdf New_features_of_Vis_2009.pptx

48 News from Geant4: Upcoming Geant4 Releases

Speaker: Dr John Apostolakis (CERN)

Slides Geant4-9.3-highlights-v4.pdf Geant4-9.3-highlights-v4.pptx

49 Topic C: Computing Performance - Presentation

Speaker: Daniel Elvira (Fermi National Accelerator Laboratory (FNAL))

Slides G4-Perf-Catania09.pdf G4-Perf-Catania09.ppt

50 Topic C: Computing Performance - Discussion

Speaker: Victor Daniel Elvira (Fermi National Accelerator Laboratory (FNAL))

Cofee/Tea Break

Technical Forum

Lunch

Steering Board Meeting

Free time

Social dinner in Catania

Sunday, 18 October

Organized excursion

Monday, 19 October

Plenary Session I - Goals, Challenges and Performance

• Workshop goals (J.Apostolakis) (30')
  o Action items
  o Delta-review outcomes
  • Emerging requirements from user domains (30' - 5' each)
    o HEP (A.Dotti)
    o Nuclear physics and underground experiments (L.Pandola)
    o Astrophysics and Space application (A.Mantero)
    o Medical application (S.Guatelli)
    o Shielding study and industrial application (T.Koi)
  • Performance improvement (K.Genser) (45')
    o What should do, what should not do

51 Workshop Goals

Speaker: Dr John Apostolakis (CERN)

Slides G4CollabMtg-Wrk2009.pdf G4CollabMtg-Wrk2009.pptx

52 Emerging requirements from user domains - HEP

Speaker: Andrea Dotti (INFN Sezione di Pisa (INFN))

Slides dotti_catania.pdf

53 Emerging requirements from user domains - Nuclear Physics and Underground experiments

Speaker: Luciano Pandola (INFN-INGS)

Slides UR_Pandola.pdf UR_Pandola.ppt

54 Emerging requirements from user domains - Astrophysics and Space Applications

Speaker: Alfonso Mantero (IN2P3/CENBG)

Slides

• SpaceUR09.odp
• SpaceUR09.pdf
• SpaceUR09.ppt

55 Emerging requirements from user domains - Medical applications

Speaker: susanna guatelli (University of Wollongong)

Slides G4MedPhysUR.pdf G4MedPhysUR.ppt

56 Emerging requirements from user domains - Shielding and industrial applications

Speaker: Dr Tatsumi Koi (SLAC)

Slides Emerging_requirements_from_Shielding_applications.pdf Emerging_requirements_from_Shielding_applications.pptx

57 On improving Geant4 performance, robustness and easing code maintenance

Fermilab Geant4 Performance Group Activities and related matters will be presented together with some suggestions and reminders regarding C++ coding techniques and style

Speaker: Krzysztof Genser (Fermilab)

Slides ImprovingG4PerformanceAndRobustness2009_genser.pdf

Coffee/Tea Break

Plenary Session II - SVN and multi-core

First assessment for migration of the Geant4 CVS code repository to SVN.
Status of current thread-safe Geant4 prototype and hints on how to implement thread-safe code in Geant4, based on the work done for the prototype, with the aim of facilitating the automatic translation process.

58 Migration to SVN

Speaker: Gunter Folger (CERN)

Slides MigrationToSubVersion.pdf

59 Multi-core & multi-threading: Report from pilot prototyping project

Speaker: Xin Dong (Northeastern University)

Slides geant4-catania-09.pdf

60 Multi-core & multi-threading: Tips on how to write "thread-safe" code in Geant4

Speaker: Gene Cooperman (Unknown)

Slides geant4-workshop-oct-09-gene.pdf

Lunch

Parallel Session I - Particle tables and processes

Particle tables and processes - coping with multi-threading

61 Context of session - and proposed agenda

** Focus ** Proposed as a working session - focusing on: * understanding the current use cases for Physics tables (in sequential Geant4 applications) * understanding the characteristics (and limitations) when used in a multi-process (copy on write) application; * understanding the characteristics and requirements in the use case of a prototype multi-threaded Geant4 application; * clarifying the requirements for Physics table(s) coming from multi-process and multi-core extensions; * discussing potential ways to overcome these current weaknesses * summarising our considerations, as a starting point for a deeper examination and potential revision of the design of the Physics tables to cope with the new requirements.

Parallel Session II - Model development and comparison with experiments

Status reports on selected models
a. INCL/ABLA
b. precompound/GEM
c. QMD
d. ENDL

Discussion: detailed vs. simplified neutron capture models

62 Developments in INCL/ABLA

Speaker: Pekka Kaitaniemi (Helsinki Institute of Physics HIP)

Slides g4-workshop-09.pdf

63 Progress in GEM and Precompound Models

Speaker: Jose Manuel Quesada Molina (Depto. de Fis. Atom., Mole. y Nucl.-Universidad de Sevilla-Unkno)

Slides Progress_GEM_PreCompound.pdf

64 QMD Update

Speaker: Dr Tatsumi Koi (SLAC)

Slides QMD_Update.pdf QMD_Update.pptx

65 The ENDL Neutron Interface

Speaker: Dr Tatsumi Koi (SLAC)

Slides The_ENDL_interface.pdf The_ENDL_interface.pptx

66 Ion-ion Validation from Catania

Speaker: Dr Giacomo Cuttone (Istituto Nazionale di Fisica Nucleare - Laboratori Nazionali del Sud, Catania, Italy)

Slides ion-ion_validation_catania.pdf

67 Discussion: Detailed vs. Simplified Neutron Capture Models

Slides CHIPScapture.pdf CHIPScapture.ppt

Coffee/Tea Break

Parallel Session III - Performance of Brachytherapy Realistic Application (work session)

• Joseph will demonstrate the a real Brachytherapy application from Univ. Laval, including real (anonymized) patient data. He will show preliminary speed profiling data. We will then discuss the example and workshop ways to speed up performance for specifically this type of application.

Parallel Session IV - Hadronic Code Improvement

Last year this type of session was quite successful. We had a good update of jobs that needed attention, and it was actually in part a working session.

1. Perhaps we can divide it up into two parts
   a. status of improvements
   b. assignment to people and working on problems

2. Topics to be covered
   a. integer Z and A
   b. G4HadronicProcess clean-up
   c. final state rotation
   d. using standard Geant4 features in models
   e. cross sections
   f. code speed up

68 Migration to integer Z and A

Speaker: Gunter Folger (CERN)

Slides A-Z_integer.pdf

69 Final state rotation

Speaker: Dr Tatsumi Koi (SLAC)

Slides Final_State_Rotation.pdf Final_State_Rotation.pptx

70 Using standard Geant4 features, code speed-up

Speaker: Dennis Herbert Wright (Stanford University)

Slides

• SpeedCode.pdf
• SpeedCode.ppt
• UsingStdG4.pdf
• UsingStdG4.ppt

71 Hadronic process clean-up, initialisation of models

Speaker: Prof. Vladimir Ivantchenko (CERN, ESA)

Slides HadronicInfrastructure.pdf

72 Discussion and hands on

Tuesday, 20 October

Plenary Session III - Hadronic modeling and validation

• Hadronic physics validation (to be named) (45')
  o Update and open issues
  • Investigation into model transition regions (to be named) (45')

73 New developments for Bertini cascade

Speaker: Dennis Herbert Wright (Stanford University)

Slides BertNewPi.pdf BertNewPi.ppt

74 Recent developments for CHIPS

Speaker: Dr Mikhail Kosov (CERN, ITEP(MOSCOW))

Slides CHIPSdevel.pdf CHIPSdevel.ppt

75 Recent developments for FTF

Speaker: Vladimir Uzhinskiy (Joint Inst. for Nuclear Research (JINR))

Slides Uzhinsky.pdf Uzhinsky.ppt

76 Validation of Geant4 models in transition energy region

Speaker: Dr Alberto Ribon (CERN)

Slides model_transition_catania09_nf.pdf

77 Discussion

Coffee/Tea Break

Plenary Session IV - Physics Lists and Metrics

• Metrics for comparing Sim/Data (to be named) (30')
  o --- Guided discussion ---
  o Existing metrics
  o New ideas and commitments
  • Physics lists (to be named) (45')
    o Open issues, new developments
    o Documentation/web: status & involvements
    o How do your choose a decent PL? Guidance & short descriptions

78 Physics Lists improvement and documentation

Topics - Status of Physics Lists ( Production HEP, Application Area / Advanced Example) - Improving Physics Lists: development and experimental Physics Lists - Documentation: News, ongoing actions, needs - New Requirements

Slides

• PhysicsLists-Oct20Plenary.pdf
• PhysicsLists-Oct20Plenary-v2.pdf
• PhysicsLists-Oct20Plenary-v2.pptx

79 Metrics for comparing Simulation against data

Lunch

Parallel Session V - Physics lists : documentation

Parallel Session VI - Low-energy EM

80 Future developments and validation of G4LowE extensions: proposed projects by CMRP

Speaker: susanna guatelli (University of Wollongong)

81 Very low energy activities at AIT

Speaker: Marcin Latocha (AIT, Austria)

82 Extended cross sections for microdosimetry

Speaker: Zyad Francis Francis

83 General discussion about workplan and contributions

Coffee/Tea Break

Parallel Session VII - Hadronics validation and testing

1. Status of validation (whatever is not covered in user session above)
   a. FNAL effort
   b. CERN effort
   c. Catania effort

2. Plans and discussion
   a. new validation suites
   b. what is needed going forward

84 Total Cross Section Tests at Process Level

Speaker: Dr Mikhail Kosov (CERN, ITEP(MOSCOW))

Slides CHIPS_XS.pdf CHIPS_XS.ppt

85 Hadronic Generator Tests for Spallation and Low Energies

Speaker: Prof. Vladimir Ivantchenko (CERN, ESA)

Slides Test30IAEA.pdf

86 The Validation Effort at FNAL

Speaker: Sunanda Banerjee (Fermilab)

Slides g421.pdf

87 Geant4 Hadronic Validation Framework

Speaker: Dr Julia Yarba (Fermilab)

Slides G4-HAD-Validation-Fwk.pdf G4-HAD-Validation-Fwk.ppt

88 Plans and Discussion

Parallel Session VIII - Kernel

• Event Biasing (M.Verderi)
  　　* Reverse MC (L.Desorgher)
  • Parallel world navigation (J.Apostolakis)
  • BREP solids by curve revolution and extrusion (G.Camellini)
    o A prototype imprementation

89 Reverse MC - status

Speaker: Laurent Desorgher (SpaceIT)

90 Parallel world navigation - status and work plan

Speaker: Dr John Apostolakis (CERN)

Slides Parallel-Navigation-v4.pdf Parallel-Navigation-v4.pptx

91 BREP solids by curve revolution and extrusion

Speaker: Gabriele Camellini (CERN)

Slides BREP_Solids_.pdf.pdf BREP_Solids.pps.pps

Wednesday, 21 October

Plenary Session V - EM Physics Validation and Development

• EM physics validation and testing
  o Update and open issues

92 New developments and validation for standard EM

Speaker: Prof. Vladimir Ivantchenko (CERN, ESA)

Slides ProgressEM.pdf

93 Migration of low-energy physics to common EM design

Speaker: Luciano Pandola (Unknown)

Slides Migration_Pandola.pdf

94 New high energy models and validation tools

Speaker: Andreas Schaelicke

Slides schaelicke_highenergyEM_v3.pdf

95 Simulating Scintillator Light Collection Using Measured Optical Reflectance

To accurately predict the light collection from a scintillating crystal through Monte Carlo simulations, it is crucial to know the angular distribution from the surface reflectance. Current Monte Carlo codes allow the user to set the optical reflectance to a linear combination of backscatter spike, specular spike, specular lobe, and Lambertian reflections. However, not all light distributions can be expressed in this way. In addition, the user seldom has the detailed knowledge about the surface that is required for accurate modeling. We have previously measured the angular distributions within BGO crystals and now incorporate these data as look-up-tables (LUTs) into modified Geant4 and GATE Monte Carlo codes. The modified codes allow the user to specify the surface treatment (ground, etched, or polished), the attached reflector (Lumirror, Teflon, ESR film, Tyvek, or TiO paint), and the bounding type (air-coupled or glued). Each LUT consists of measured angular distributions with 4 deg by 5 deg resolution in theta and phi, respectively, for incidence angles from 0 deg to 90 deg, in 1 deg steps. We compared the new codes to the original codes by running simulations with a 0.3x1.0x3.0mm**3 BGO crystal coupled to a PMT. The simulations were then compared to measurements. Light output was measured by counting the photons detected by the PMT with the 0.3x1.0, 0.3x3.0, or 1.0x1.0 mm**2 side coupled to the PMT, respectively. Our new code shows better agreement with the measured data than the current Geant4 code. The new code can also simulate reflector materials that are not pure specular or Lambertian reflectors, as was previously required. Our code is also more user friendly, as no detailed knowledge about the surfaces or light distributions is required from the user.

Speaker: Peter Gumplinger

Slides MJanecek.pdf MJanecek.ppt

96 General Discussion

Coffee/Tea Break

Plenary Session VI - Testing, Analysis and CAD/GDML

• New system testing (Test49:Geant4+ROOT) (30')
  • Analysis: Common Solutions for G4 tests, examples (to be named) (60')

slides test49.pdf test49.ppt

97 A stable interface to read and write IAEA phase-space files in Geant4

Speaker: Mr Miguel A. Cortes-Giraldo (University of Sevilla)

Slides G4IAEA_G4WS09.pdf G4IAEA_G4WS09.ppt

98 test49:Geant4+ROOT

Slides test49.pdf test49.ppt

99 Getting CAD Geometries into Geant4 - Discussion

Lunch

Parallel Session IX - Multiple scattering and single scattering

100 Incorporation of Goudsmit-Sounderson electron transport theory into Geant4

Speaker: Dr Omrane Kadri (National Center for Nuclear Sciences and Technologies)

Slides G4GS.pdf

101 Multiple scattering model upgrade

Speaker: Laszlo Urban (Unknown)

Slides Laszlo_msc.pdf

102 MSC benchmarking

Speaker: Joseph Perl (Unknown)

103 Discussion on MSC configuration

Slides MSC.pdf

Parallel Session X - Configure tool / New system testing

104 Replacing Configure/metaconf by modern tool(s)

Speaker: Dr Ben Morgan (University of Warwick)

Slides benmorgan-confreplacement.pdf

105 System testing using nightlies

Slides SystemIntregrationTesting.pdf

Coffee/Tea Break

Parallel Session XI - EM - new models and validation

most presentations will be given remotely

106 Range validation for electrons, protons and alpha particles

Speakers: Christina Zacharatou Jarlskog (Niels Bohr Institute, Copenhagen), Michel Maire (LAPP)

107 ICRU'73 stopping powers

Speaker: Anton Lechner (Atominst. der Oest. Universitaeten - Technische Universitaet Wien)

Slides alechner_ions.pdf

108 Fluctuation model modifications

Speakers: Laszlo Urban (Unknown), Michel Maire (Unknown)

Slides urban_fluct.pdf

109 New test on silicon detectors

Speaker: Frederic Dupertuis

Slides Validation_of_EM_Models.pdf

110 Discussion - Physics List options for 9.3

Speaker: Prof. Vladimir Ivantchenko (CERN, ESA)

Slides PhysListOptions.pdf

Parallel Session XII - [Available slot]

Thursday, 22 October

Plenary Session VII - Summary from the parallel sessions

111 II - Model development and comparison with experiments

Speaker: Dr Mikhail Kosov (CERN, ITEP(MOSCOW))

Slides M.Kossov.pdf M.Kossov.ppt

112 III - Performance of Brachytherapy Realistic Application (work session)

Speaker: Joseph Perl

Slides

• Catania_Ad_Hoc_Parallel_Visualization.pdf
• Catania_Ad_Hoc_Parallel_Visualization.ppt
• Catania_Parallel_III.pdf
• Catania_Parallel_III.ppt

113 IV - Hadronic Code Improvement

Speaker: Prof. Vladimir Ivantchenko (CERN, ESA)

Slides SumHadronicSum.pdf

114 V - Physics lists : documentation and X - Configure tool / New system testing

Speaker: Gunter Folger (CERN)

Slides SummarySessions_v_x.pdf

115 VI - Low-energy EM

Speaker: Sebastien Laurent Incerti

Slides S.Incerti.pdf S.Incerti.ppt

116 I - Physics Tables - and Multi-core

Speaker: Dr John Apostolakis (CERN)

Slides Parallel-Physics-TablesMT.pdf Parallel-Physics-TablesMT.pptx

117 VII - Hadronics validation and testing

Speaker: Prof. Vladimir Ivantchenko (CERN, ESA)

Slides S.Banerjee.pdf S.Banerjee.ppt

118 VIII - Kernel

Speaker: Makoto Asai (SLAC)

Slides Parallel_session_VIII.pdf Parallel_session_VIII.ppt

119 IX - Multiple scattering and single scattering

Speaker: Prof. Vladimir Ivantchenko (CERN, ESA)

Slides SummaryMSC.pdf

120 XI - EM - new models and validation

Speaker: Andreas Schaelicke

Slides SummaryXI.pdf

Coffee/Tea Break

Plenary Session VIII - Future developments and challenges

• Next release (G.Cosmo)
  • Architecture and design review (M.Asai)
  • Outlook (J.Apostolakis)

121 Next release

Speaker: Gabriele Cosmo

Slides Release-schedule-9.3.pdf Release-schedule-9.3.ppt

122 Architecture review

Speaker: Makoto Asai (SLAC)

Slides ArchitectureReview.pdf ArchitectureReview.ppt

123 Outlook

Speaker: Dr John Apostolakis (CERN)

Paper Outlook2010.pdf Outlook2010.pptx

Workshop closeout

Steering Board Meeting