Simulation capabilities and application results http cern chgeant
Simulation capabilities and application results http: //cern. ch/geant 4. html Maria Grazia Pia INFN Genova on behalf of the Geant 4 Collaboration EPS-HEP 2001 Conference Budapest, 12 -18 July 2001 Maria Grazia Pia, INFN Genova - EPS-HEP 2001 Budker Inst. of Physics IHEP Protvino MEPHI Moscow Pittsburg University
ATLAS Ba. Bar highlights An extensive set of physics processes and models over a wide energy range Courtesy of L 3 Photon attenuation Courtesy of the Italian Nat. Inst. for Cancer Research a particle in a cell 192 UKDM, Boulby Mine Ir E (Me. V) High energy m Low energy photons A rigorous approach to software engineering Maria Grazia Pia, INFN Genova - EPS-HEP 2001 GLAST
A wide domain of applications with a large user community in many fields Geant 4 is a simulation Toolkit designed for a variety of applications It adopts rigorous software engineering methodologies and is based on OO technology It has been developed and is maintained by an international collaboration of > 100 scientists HEP, astrophysics, nuclear physics, space sciences, medical physics, radiation studies etc. It provides a complete set of tools for all the typical domains of simulation - run, event and track management - geometry and materials - tracking - detector response - RD 44 Collaboration (1994 -98) - PDG-compliant particle management - Geant 4 Collaboration - user interface The code is publicly distributed from the WWW, together with ample documentation 1 st production release: end 1998 - 2 new releases/year since then Maria Grazia Pia, INFN Genova - EPS-HEP 2001 - visualisation - persistency - physics processes
Geant 4 architecture Software Engineering plays a fundamental role in Geant 4 Domain decomposition User Requirements • formally collected • systematically updated • PSS-05 standard Software Process hierarchical • spiral iterative approach structure of sub • regular assessments and improvements • monitored following the ISO 15504 model -domains • OOAD Object Oriented methods • use of CASE tools • essential for distributed parallel development Uni-directional flow of • contribute to the transparency of physics dependencies Maria Grazia Pia, INFN Genova - EPS-HEP 2001 • commercial tools Quality Assurance • code inspections • automatic checks of coding guidelines • testing procedures at unit and integration level • dedicated testing team Use of Standards • de jure and de facto
Geometry Role: detailed detector description and efficient navigation Chandra ATLAS Ba. Bar Multiple representations (Same abstract interface) CSG (Constructed Solid Geometries) - simple solids XMMNewton STEP extensions - polyhedra, , spheres, cylinders, cones, toroids, etc. BREPS (Boundary REPresented Solids) - volumes defined by boundary surfaces - include solids defined by NURBS (Non-Uniform Rational B-Splines) CAD exchange: ISO STEP interface Fields: of variable non-uniformity and differentiability Borexino External tool for g 3 tog 4 geometry conversion Maria Grazia Pia, INFN Genova - EPS-HEP 2001 CMS
Guidelines for physics From the Minutes of LCB (LHCC Computing Board) meeting on 21 October, 1997: “It was noted that experiments have requirements for independent, alternative physics models. In Geant 4 these models, differently from the concept of packages, allow the user to understand how the results are produced, and hence improve the physics validation. Geant 4 is developed with a modular architecture and is the ideal framework where existing components are integrated and new models continue to be developed. ” Geant 4 physics keeps evolving Maria Grazia Pia, INFN Genova - EPS-HEP 2001 with attention to UR facilitated by the OO technology
Geant 4 Physics OOD allows to implement or modify any physics process without changing other parts of the software open to extension and evolution Tracking is independent from the physics processes The generation of the final state is independent from the access and use of cross sections Transparent access via virtual functions to - cross sections (formulae, data sets etc. ) - models underlying physics processes An abundant set of electromagnetic and hadronic physics processes ´ a variety of complementary and alternative physics models for most processes Use of public evaluated databases No tracking cuts, only production thresholds - thresholds for producing secondaries are expressed in range, range universal for all media converted into energy for each particle and material The transparency of the physics implementation contributes to the validation of experimental physics results Maria Grazia Pia, INFN Genova - EPS-HEP 2001
Electromagnetic physics energy loss § electrons and positrons n g, X-ray and optical photons n muons n charged hadrons n ions 11. 12. 13. 14. 15. 16. 17. 18. Comparable to Geant 3 already in the 1 st a release (1997) Further extensions (facilitated by the OO technology) î High energy extensions - 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. needed for LHC experiments, cosmic ray experiments… Multiple scattering Bremsstrahlung Ionisation Annihilation Photoelectric effect Compton scattering Rayleigh effect conversion e+e- pair production Synchrotron radiation Transition radiation Cherenkov Refraction Reflection Absorption Scintillation Fluorescence Auger (in progress) î Low energy extensions - fundamental for space and medical applications, experiments, antimatter spectroscopy etc. n î Alternative models for the same process All obeying to the same abstract Process interface transparent to tracking Maria Grazia Pia, INFN Genova - EPS-HEP 2001
Standard e. m. processes 1 ke. V up to O(100 Te. V) Multiple scattering - new model (by L. Urbán) - computes mean free path length and lateral displacement New energy loss algorithm - optimises the generation of d rays near boundaries Variety of models for ionisation and energy loss - including the Photo. Absorption Interaction model Multiple scattering 6. 56 Me. V proton , 92. 6 mm Si Geant 4 Geant 3 data Differential and Integral approach - for ionisation, Bremsstrahlung, positron annihilation, energy loss and multiple scattering J. Vincour and P. Bem Nucl. Instr. Meth. 148. (1978) 399 Maria Grazia Pia, INFN Genova - EPS-HEP 2001
shell effects e, down to 250 e. V Low energy e. m. extensions Fundamental for neutrino/dark matter experiments, space and medical applications, antimatter spectroscopy etc. (EGS 4, ITS to 1 ke. V, Geant 3 to 10 ke. V) Based on EPDL 97, EEDL and EADL evaluated data libraries Photon attenuation Barkas effect (charge dependence) models for negative hadrons protons antiprotons Maria Grazia Pia, INFN Genova - EPS-HEP 2001 Bragg peak Hadron and ion models based on Ziegler and ICRU data and parameterisations
Muons 1 ke. V up to 1000 Pe. V scale simulation of ultra-high energy and cosmic ray physics High energy extensions based on theoretical models Optical photons Production of optical photons in HEP detectors is mainly due to Cherenkov effect and scintillation Processes in Geant 4: - in-flight absorption Rayleigh scattering medium-boundary interactions (reflection, refraction) Photon entering a light concentrator CTF-Borexino Maria Grazia Pia, INFN Genova - EPS-HEP 2001
Parameterised and data-driven hadronic models (1) Based on experimental data Some models originally from GHEISHA - completely reengineered into OO design refined physics parameterisations New parameterisations - pp, elastic differential cross section n. N, total cross section p. N, total cross section np, elastic differential cross section N, total cross section N, coherent elastic scattering Maria Grazia Pia, INFN Genova - EPS-HEP 2001 p elastic scattering on Hydrogen
Parameterised and data-driven hadronic models (2) Other models are completely new, such as: Stopping p stopping particles: - , K- (relevant for m/ PID detectors) absorption Isotope production neutrons nuclear deexcitation Me. V Neutrons Courtesy of CMS Energy All worldwide existing databases used in neutron transport Brond, CENDL, EFF, ENDFB, JEF, JENDL, MENDL etc. Maria Grazia Pia, INFN Genova - EPS-HEP 2001
Theory-driven models Complementary and alternative models Discrete transitions from ENSDF data Geant 4 Theoretical model for continuum Giant Dipole Resonance Evaporation phase Low energy range, pre-equilibrium, O(100 Me. V) Intermediate energy range, O(100 Me. V) to O(5 Ge. V), intra-nuclear transport High energy range, hadronic generator régime Maria Grazia Pia, INFN Genova - EPS-HEP 2001
Other components Materials - elements, isotopes, compounds, chemical formulae Visualisation - Particles - all PDG data and more, for specific Geant 4 use, like ions User Interfaces - Hits & Digi - to describe detector response - Persistency - possibility to run in transient or persistent mode no dependence on any specific persistency model persistency handled through abstract interfaces to ODBMS Maria Grazia Pia, INFN Genova - EPS-HEP 2001 Various drivers Open. GL, Open. Inventor, X 11, Postscript, DAWN, OPACS, VRML Command-line, Tcl/Tk, Tcl/Java, batch+macros, OPACS, GAG, MOMO automatic code generation for geometry and materials Interface to Event Generators - through ASCII file for generators supporting /HEPEVT/ abstract interface to Lund++
Modules for space applications Delayed radioactivity Particle source and spectrum General purpose source particle module INTEGRAL and other science missions Low-energy e. m. extensions Geological surveys of asteroids Sector Shielding Analysis Tool CAD tool front-end Instrument design purposes Maria Grazia Pia, INFN Genova - EPS-HEP 2001 Dose calculations Courtesy of P. Nieminen, ESA
Interface to external tools Through abstract interfaces AIDA No dependence Minimize coupling of components Courtesy of A. Pfeiffer, CERN Example: AIDA & Analysis Tools Similar approach: Lizard Maria Grazia Pia, INFN Genova - EPS-HEP 2001 Java Analysis Studio graphics (G)UI persistency etc.
Ba. Bar Courtesy of D. Wright for the Ba. Bar Collaboration Preliminary Maria Grazia Pia, INFN Genova - EPS-HEP 2001
Example of integrated Fast/Full Simulation application 1. Ba. Bar Object-oriented Geant 4 -based Unified Simulation (BOGUS) 2. Integrated framework for Fast and Full simulation Fast simulation available for public use since February 1999 Integrated in Ba. Bar environment 3. 4. 1. primary generators, reconstruction, OODB persistency 2. parameters for materials and geometry shared with reconstruction applications Exploits Geant 4 parameterisation (new feature) Courtesy of G. Cosmo for the Ba. Bar Collaboration Maria Grazia Pia, INFN Genova - EPS-HEP 2001
ATLAS 300 Ge. V muons 20 Ge. V pions TRT: Energy loss measured in ATLAS test beam compared to Geant 3 and Geant 4 simulations (PAI model) Preliminary Liquid Ar calorimeter Fcal energy resolution Muon detector Courtesy of D. Barberis for ATLAS Collaboration Maria Grazia Pia, INFN Genova - EPS-HEP 2001
HARP with GEANT 4 Courtesy of P. Arce for the HARP Collaboration Maria Grazia Pia, INFN Genova - EPS-HEP 2001
Preliminary Sophisticated geometry Very non-uniform strong magnetic field Primary target as a particle source T 9 beam line Simulation (10 Ge. V/c) Measurement Beam spot width (mm) 3. 27 approx. 4 Beam spot height (mm) 3. 49 approx. 4 Beam spot position (mm) (0. 33: 0. 86) (0. 0: 0. 0) Beam profile and composition at the HARP target Maria Grazia Pia, INFN Genova - EPS-HEP 2001 Crucial to have a precise absolute knowledge of the particle rate incident onto HARP target Impossible to separate experimentally from m in the beam with the accuracy required Courtesy of P. Arce for the HARP Collaboration
GLAST (g-ray telescope) Preliminary GLAST Maria Grazia Pia, INFN Genova - EPS-HEP 2001 Courtesy of F. Longo and R. Giannitrapani, GLAST
Cosmic rays, jovian electrons Other astroparticle applications Solar system explorations Courtesy SOHO EIT Solar X-rays, e, p unique simulation capabilities: low E physics fluorescence radioactivity neutrons space modules etc. . Courtesy of R. Nartallo, ESA X-ray telescope Courtesy of S. Magni, Borexino XMM ZEPLIN III Dark Matter, Boulby mine Maria Grazia Pia, INFN Genova - EPS-HEP 2001 Courtesy of A. Howard, UKDM
anisotropy Technology transfer Medical applications of Geant 4: radiotherapy PET dosimetry etc. Brachytherapy 192 Ir Treatment planning Courtesy National Inst. for Cancer Research, Genova Courtesy LIP & Portuguese Oncological Institute Isodoses data Histogram: Geant 4 Commercial treatment planning system Maria Grazia Pia, INFN Genova - EPS-HEP 2001
Conclusions Geant 4 is a simulation Toolkit, providing advanced tools for all the domains of detector simulation Geant 4 is characterized by a rigorous approach to software engineering Thanks to the OO technology, Geant 4 is open to extension and evolution An abundant set of physics processes is available, often with a variety of complementary and alternative physics models Its areas of application span diverse fields: HEP and nuclear physics, astrophysics and space sciences, medical physics, radiation studies etc. Maria Grazia Pia, INFN Genova - EPS-HEP 2001
- Slides: 26