1 / 20

New Hadr02 test for Ion-Ion simulation

This article introduces the new Hadr02 example for ion-ion simulations in Geant4. It discusses the results and comparisons with existing models, as well as the implementation and problems encountered.

dcook
Download Presentation

New Hadr02 test for Ion-Ion simulation

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. New Hadr02 test for Ion-Ion simulation Ivantchenko V., Ivantchenko A.

  2. Outline • Introduction • New example Hadr02 • Results and comparison • Discussions

  3. Introduction

  4. Ion/Ion models in G4 9.4 • There are many requests to have ion-ion simulation within Geant4 • ESA • Recently NA61 • Models existing • BinaryLightIon for all Ions up to 15 GeV/u • LHEP and CHIPS only for light ions • G4QMD for all ions up to 15 GeV/u (very slow)

  5. DPMJET-II.5 interface • High energies (>10GeV/u) Ion/Ion interactions • Software is provided by P. Truscott – it is a version prepared for usage in Geant4 during MARS-REM project sponsored by ESA few years ago from the original DPMJET-II.5 FLUKA code • Our experience in DPMJET within the GRAS tool • Initially installation at SLC4 32 bit gcc3.4.6 was achieved • Installation completely fails at SLC5 default compiler gcc4.1 • Recently installation at SLC5 64 bit gcc4.3.2 was achieved

  6. DPMJET installation in Hadr02

  7. Our Method • For DPMJET interface: • Different from recommendations of MARS-REM • In geant4/examples/extended/hadronic • Extra subdirectory “dpmjet2_5” • Extra subdirectory “GlauberData” • Extended makefiles of the example • Our favorite is the FTF model • Is available since 9.4ref02

  8. Details of implementation • Example is cloned from Hadr01 with some simplification • dpmjet2.5 – separate directory for DPMJET source files • dpmjet2.5/dpm25src – FORTRAN files • dpmjet2.5/include – C++ header files for DPMJET interface • dpmjet2.5/src– source files for DPMJET interface • GlauberData – Glauber cross sections database • $GRAS_USE_DPMJET - new environment variable • install.csh , setup.csh - extra files for installation • Physics Lists details see in Discussions

  9. Results and comparison for 9.4.ref02

  10. First try • Target 10 cm Al • Primary beam 6400 GeV • Ion S32 (Z = 16, A = 32) • CPU for DPMJET and FTFP are similar • Cut 1 km (no delta-electrons) • All secondary particles are killed (no shower) !

  11. Energy log10(GeV) of secondary's d, t, Ions

  12. Z and A of fragments

  13. Summary • FTF and DPMJET are working inside Hadr02 • CPU and physics performance are similar • FTF has no projectile fragmentation

  14. DISCUSSION of problems

  15. Problems of DPMJET • Installation is achieved at SLC5 64 bits gcc 4.3 • Extra libraries are used: • CERNLIB, gmp, mpfr • At pcgeant05 we have a warning : local gmp version is different from one of LCG • GlauberData for DPMJET are limited • He(Z=2) – Ce( Z = 58) • FPE_DEBUG is off, because inside DPMJET there is at least one division by zero

  16. Physics List of Hadr02 • Only reference Physics List is used • Two additional ion physics builders are created • IonDPMJETPhysics.cc for DPMJET • IonFTFPhysics.cc for FTF • Low energy ion-ion models are not included in these builders • These builders are registered additionally to any existing reference Physics List via UI command • This is kind of limitation, because some ion physics is already exist • Results below are obtained with QGSP_BIC

  17. Problems of Physics Lists • We would like to suggest to extended public interface to G4VModularPhysicsList allowing to substitute existing physics builder by a new one provided by users • Likely this will require to add enumeration for Physics List builders

More Related