Download
1 / 7
70 likes | 77 Views
Charged Particle HP is a package that uses evaluated nuclear data libraries to simulate charged particle interactions, including total cross sections, inelastic channel cross sections, double differential spectra of outgoing particles, and gamma emission due to nuclear level transitions. This package is an alternative to other models when experimental data cannot be reproduced. The code is compatible with neutron_hp with minor modifications.
E N D
Charged particle HP • Pedro Arce • Emilio Mendoza • (CIEMAT, Madrid) • Patrick Sauvan • (UNED, Madrid)
What is charged_particle_hp? • neutron_hp package uses evaluated nuclear data bases for neutron interactions: • Total cross sections • Inelastic channel cross sections • Double differential spectra of outgoing particles • Gamma emission because of nuclear level transitions • charged_particle_hp: do the same for (inelastic) interactions of other particles (p, d, t, He3, , )
Why charged_particle_hp? Theory models or semi-empirical models sometimes cannot reproduce experimental data (J.M. Quesada agrees): p + 18O n + 18F (barns) IAEA medical DB (exper.) Geant4.9.4.p01 (QGSP_BIC_HP) Energy (MeV)
What nuclear DBs are there? • ENDF-VII: • Uses only experimental data • Only a few isotopes (p:48, d:5, t:3, He3:2 • Only p X reactions (MT=5) • double differential spectra of resulting particles (n,p,d,…), without channel information (n,nn,np,nna,...) • Up to 150 MeV for p (d: 50 MeV, t: 20 MeV, He3: 20 MeV) • ENDF format • TENDL: • Uses some experimental data + TALYS calculations • All isotopes (2400) • All channels (also available a DB with only p X reactions) • Up to 200 MeV • ENDF format IAEA medical database: • Only experimental data • Only a few reaction channels of a few isotopes • Only channel cross sections • Simple text format E. Mendoza translated ENDF format into GEANT4 format
Modifications in the C++ code • No new classes needed • neutron_hp classes modified so that: • - Not assumed projectile particle is neutron • - Not assumed data base directory is given by $G4NEUTRONHPDATA • - Do not instantiate Fission/Capture/Elastic if projectile is not neutron • - Do not use thermal corrections if projectile is not neutron • Add projectile and directory in your physics list • Default is like old neutron_hp: no modifications needed in old physics lists • G4NeutronHPInelastic(G4ParticleDefinition* projectile = G4Neutron::Neutron(), • const G4String& dataDir = "G4NEUTRONHPDATA", const G4String& name = "NeutronHPInelastic" );: • G4NeutronHPInelasticData(G4ParticleDefinition* aPar = G4Neutron::Neutron(),con • st char* envVar = "G4NEUTRONHPDATA");
Tests to do • Auto consistency: • Test that web data are reproduced by running GEANT4 • Check ENDF format Geant4 format • Check Geant4 code • For each isotope: • * Tabular angular distributions • * Tabular residual cross sections • * Tabular spectra • * Tabular total and partial cross sections • Compare with MCNP: • * Cross sections • * Double differential spectra of resulting particles • * Channel by channel comparisons? (extracting data from MCNP is not so easy….) • Collaboration with MCUNED people (UNED, Madrid) people (experts in charged particle interactions with MCNP) Test have been just started
Conclusions • charged_particle_hp package uses evaluated nuclear data libraries for p, d, t, He3, , inelastic interactions (= neutron_hp) • Alternative to other models (some cases others are not OK) • Nuclear data available up to 200 MeV • Reuse neutron_hp code with small modifications • Code is working, detailed tests have started
