Interfacing the JQMD to Geant4 by rt3463df


									Interfacing the JQMD and JAM
Nuclear Reaction Codes to
   Stanford Linear Accelerator Center
             Koi, Tatsumi

            Geant4 2003 Workshop @ TRIUMF
Collaboration with

  K. Niita (RIST)
  Y. Nara (University of Arizona)
  M. Asai (SLAC)
  D. H. Wright (SLAC)
  T. Sasaki (SLAC/KEK)
  K. Amako (KEK)

              Geant4 2003 Workshop @ TRIUMF
Introducing joint activity of SLAC, KEK and
other institutes for interfacing Fortran code of
JQMD and JAM to Geant4
- To satisfy the urgent requirements for the beam
test simulation of GLAST experiment.
  What is JQMD
  What is JAM
  Interfacing Fortran reaction codes to Geant4

                Geant4 2003 Workshop @ TRIUMF
  We connected Fortran nuclear reaction codes
      to Geant4 which is written in C++.

 Advantages of this method are
 There are already many well-established reaction
 codes and these codes are often written in Fortran.
 It is more convenient interfacing to Fortran code
 directly than re-writing the code in C++.
 In the process of re-writing, new bugs may
enter into the code. We can avoid this situation.
 Once the interface is established, the Fortran
 code and Geant4 can be updated independently.
 No copyright problems associated with re-writes.
                Geant4 2003 Workshop @ TRIUMF
At the first stage of this activity
we chose two reaction codes, which
satisfy our requirement
(Nucleus-Nucleus Reactions)

Jaeri Quantum Molecular Dynamics
Jet AA Microscopic Transport Model
(JAM)       Geant4 2003 Workshop @ TRIUMF
 QMD (Quantum Molecular Dynamics) is
quantum extension of classical molecular-
 dynamics model.
 QMD model is widely used to analyze
various aspects of heavy ion reactions.
 JQMD ( Jaeri QMD) is a QMD code
 developed by JAERI
( Japan Atomic Energy Research Institute).
 JQMD includes SDM (Statistical Decay Model),
 SDM is used for evaporation and fission decays
 of excited nuclei.

             Geant4 2003 Workshop @ TRIUMF
JQMD (cont.)
 Energy Range of JQMD
 several 10 MeV/N up to about 3 GeV/N
 Projectile particle species
 nuclei (including protons and neutrons) and pions
 Detailed description of JQMD is given in
Niita et al., Physical Review C 52 (1995) 2620
 JQMD is also run within PHITS (Particle and
 Heavy-Ion Transport code System).
Iwase et al., Journal of Nuclear Science and
 Technology 39 (2002) 1142

              Geant4 2003 Workshop @ TRIUMF
JAM (Jet AA Microscopic Transportation
Model) is a hadronic cascade model.
The trajectories of all hadrons as well as
resonances including produced particles are
followed explicitly as a function of space and time.
The inelastic hadron-hadron collisions are
described by resonance formation and decay at
low energies (below ~4GeV). Above 4GeV string
formation and fragmentation into hadrons is
Multiple minijet production is also included at
high energies in the same way as the HIJING
(Heavy Ion Jet INteraction Generator)
               Geant4 2003 Workshop @ TRIUMF
 JAM (cont.)
  Energy Range of JAM
  several 100 MeV/N up to about 100 GeV/N
  and valid for many reactions at LHC and RHIC
  Projectile particle species
  nuclei (including protons and neutrons)
 mesons pions and kaons
  Detailed description of JAM is given in
Nara et al., Physical Review C 561 (1999) 4901

             Geant4 2003 Workshop @ TRIUMF
In order to use JQMD and JAM
 from Geant4,
we had to make interfaces.

Before discussing interfaces,
it is useful to review Geant4
 processes and how they handle
hadronic interactions.
           Geant4 2003 Workshop @ TRIUMF
What Does a Process Do in
 Decides when and where an interaction will
 occur (GetPhysicalInteractionLength)
   Cross section
 Generates the final state (DoIt)
   Reaction model

              Geant4 2003 Workshop @ TRIUMF
What Does a Process Do in Geant4 (cont.)

                                              Cross Section


                    Geant4 2003 Workshop @ TRIUMF
Hadron inelastic models
JAM2G4InelasticModel ::ApplyYourself()

In this function,
Fortran JQMD or JAM routine is called.

             Geant4 2003 Workshop @ TRIUMF
Framework of Hadron Interactions
in Geant4               Fortran JQMD routine
                 PostStepDoIt()                   JQMD2G4InelasticModel

           G4HadronicDiscreteProcess                   G4HadronicInteracton
                  GetCrossSection()                       ApplyYourSelf()

CrossSection                                      JAM2G4InelasticModel
                       Geant4 2003 Workshop @ TRIUMF   Fortran JAM routine
Cross Section for heavy ions
 Triphathi Formula
   NASA Technical Paper 3621 (1997)
 Shen Formula
   Nuclear Physics. A 491 (1989) 130

             Geant4 2003 Workshop @ TRIUMF
Platforms used

 Tested system
   OS: Red Hat Linux 7.2 (6.2)
   Compiler: gcc-2.95.3 (2.91.66 with egcs-1.1.2)
   Geant4: Ver. 5.0.p01 (2003 Feb version.)
 Demonstration System
   Cygwin 1.3.22-1 with gcc 3.2.3
 We did not test this interface on other OS and
 compilers. Perhaps small modifications will be

               Geant4 2003 Workshop @ TRIUMF
   Based on Geant4 novice exampleN03
   Based on Geant4 novice exampleN03
  ICRU Sphere
     International Commission on Radiation Units
   Heavy ion passage in the tissue equivalent material

                 Geant4 2003 Workshop @ TRIUMF
Simulation snapshot 1

       Geant4 2003 Workshop @ TRIUMF
Simulation snapshot 2

       Geant4 2003 Workshop @ TRIUMF
Simulation snapshot 3

       Geant4 2003 Workshop @ TRIUMF
Pion Plus from 17.5 GeV/c Proton on Be

          cos(θ)=0.95                        cos(θ)=0.85

            cos(θ)=0.75                      cos(θ)=0.65

             Geant4 2003 Workshop @ TRIUMF            I. Chemakin et al.,
                                                      Phys. Rev. C 65, 4904 (2002)
Pion from nucleus-nucleus interactions

           Sugaya et al.,
           Nucl. Phys. A634, 115 (1998)

           Geant4 2003 Workshop @ TRIUMF
Pion from nucleus-nucleus interaction


            Geant4 2003 Workshop @ TRIUMF
 We successfully developed the interface
 which connect JQMD and JAM to Geant4.
 Hadronic framework of Geant4 proved its
 flexibility and expandability .
 With this interface, we can simulate
 propagation of heavy ions in complex
 Geant4 geometries.
 Preliminary test results agree well with
 data. Much more validation to be done.
            Geant4 2003 Workshop @ TRIUMF

To top