200 likes | 356 Views
New Hadr02 test for Ion-Ion simulation. Ivantchenko V., Ivantchenko A. Outline. Introduction New example Hadr02 Results and comparison Discussions. Introduction. Ion/Ion models in G4 9.4. There are many requests to have ion-ion simulation within Geant4 ESA Recently NA61
E N D
New Hadr02 test for Ion-Ion simulation Ivantchenko V., Ivantchenko A.
Outline • Introduction • New example Hadr02 • Results and comparison • Discussions
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)
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
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
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
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) !
Summary • FTF and DPMJET are working inside Hadr02 • CPU and physics performance are similar • FTF has no projectile fragmentation
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
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
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