230 likes | 411 Views
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.
E N D
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