Modeling Production, Interactions and Transport
230 likes | 419 Vues
ILC-CAL. Fermilab. Modeling Production, Interactions and Transport. Nikolai Mokhov, Fermilab. Fermilab November 14, 2005. OUTLINE. General Purpose Event Generators and Transport Codes Neutron Production and Interactions Transport at High and Low Energies.
Modeling Production, Interactions and Transport
E N D
Presentation Transcript
ILC-CAL Fermilab Modeling Production, Interactions and Transport Nikolai Mokhov, Fermilab Fermilab November 14, 2005
OUTLINE • General Purpose Event Generators and Transport Codes • Neutron Production and Interactions • Transport at High and Low Energies n/nbar modeling - N.V. Mokhov
GENERAL-PURPOSE CODES (2005) • Criteria for this express-analysis: • General-purpose widely-used code. • Event generators with consistent soft and hard multiparticle production in hA , ( AA ), g A and n A at MeV to many TeV. • Full particle transport from a fraction of eV to many TeV in complex target/ accelerator/ detector/ shielding systems. • Codes are actively developed, supported and discussed at conferences last decade. n/nbar modeling - N.V. Mokhov
HADRONIC CODES • FLUKA – INFN: DPMJET, PEANUT • GEANT4 – CERN: INC, QMD • MARS15 – FNAL: CEM, LAQGSM, DPMJET, Inclusive, MCNP • MCNPX – LANL: Tables, CEM, no EMS, no magnetic field, work towards high energies in progress • PHITS – JAERI: QMD • SHIELD – INP Troitsk: INC, QGSM, no EMS, no magnetic field, simple geometry n/nbar modeling - N.V. Mokhov
SOME COMMENTS • Current ‘2005 FLUKA and MARS15 are in a very good agreement with data and with each other in all applications • Other codes also perform well for neutronics • Antineutron yields and interactions are reasonably well understood and described by codes. E.g., DPMJET yields for 100-GeV proton on iron: 5.9 protons, 7.4 neutrons, 0.075 pbars, and 0.0735 nbars n/nbar modeling - N.V. Mokhov
NEUTRON YIELD: GEANT4, 256-MeV p+Fe n/nbar modeling - N.V. Mokhov
NEUTRON YIELD: GEANT4, 597-MeV p+Fe n/nbar modeling - N.V. Mokhov
MARS15 EXCLUSIVE EVENT GENERATORS • The improved Cascade-Exciton Model code, CEM03.01, combined with the Fermi break-up model, the coalescence model, and an improved version of the Generalized Evaporation-fission Model (GEM2) is used as a default for hadron-nucleus interactions below 5 GeV. • The Los Alamos Quark-Gluon String Model code, LAQGSM03, was implemented into MARS15 for particle and heavy-ion projectiles at 10 MeV/A to 800 GeV/A. This provides a power of full theoretically consistent modeling of exclusive and inclusive distributions of secondary particles, spallation, fission, and fragmentation products. Further development of this package is underway. • For quite some time, MARS has used the Dual-Parton Model code, DPMJET3,for the very first vertex in a cascade tree. This is used in our numerous studies for the LHC 7x7 TeV collider and its detectors, and at very high energies up to 100 TeV. n/nbar modeling - N.V. Mokhov
NEUTRONS IN CEM03/MARS15 n/nbar modeling - N.V. Mokhov
NEUTRONS IN CEM03/MARS15 n/nbar modeling - N.V. Mokhov
LAQGSM/MARS15: 0.56 and 8 GeV/A REACTIONS n/nbar modeling - N.V. Mokhov
MARS INCLUSIVE EVENT GENERATOR (1) n/nbar modeling - N.V. Mokhov
MARS INCLUSIVE EVENT GENERATOR (2) n/nbar modeling - N.V. Mokhov
INCLUSIVE NUCLEON PRODUCTION IN MARS15 MARS15 vs data n/nbar modeling - N.V. Mokhov
MARS15: NEUTRON ELASTIC SCATTERING n/nbar modeling - N.V. Mokhov
MARS15 NEUTRON AND NBAR X-SECTIONS n/nbar modeling - N.V. Mokhov
NEUTRON SPECTRA FOR 2-GEV ELECTRON BEAM n/nbar modeling - N.V. Mokhov
NEUTRON YIELD FROM LEAD CYLINDER n/nbar modeling - N.V. Mokhov
MARS, FLUKA and other codes for in thick targets and shielding at 0.8, 10 and 50 GeV n/nbar modeling - N.V. Mokhov
Results of PHITS JQMD code inPHITS JQMD-2 Introducing JQMD in PHITS :H. Iwase et.al. J. Nucl. Sci.Technol.39 (2002) 1142 Neutron Spectra from Thick Target n/nbar modeling - N.V. Mokhov
4° 45° 90° 120° MCNPX PHITS Ep=12GeV Cu:1mmt DDX of Secondary Particles Matsuda (1) p(12GeV) + Cu -> neutron neutron calculated by N. Matsuda n/nbar modeling - N.V. Mokhov
Depth Dose Distribution in Cu 10cm radius Ohnishi (2) calculated by H. Ohnishi. in Dr. thesis n/nbar modeling - N.V. Mokhov
Depth Dose Distribution in Cu 10cm radius Ohnishi (3) calculated by H. Ohnishi. in Dr. thesis n/nbar modeling - N.V. Mokhov