Hadr03 nuclear reactions QGSP_BIC_HP Vs QGSP_BIC_AllHP

Getting Started
1

Hello,
I am new to Geant4 and I am currently running into some troubles when trying to run the example Hadr03 using two different physic lists, namely QGSP_BIC_HP and QGSP_BIC_AllHP.

QGSP_BIC_HP:
plenty of nuclear reaction are listed along with the produced particles
QGSP_BIC_AllHP
only one reaction listed: proton+Mo98 → N gamma or e- + 2 neutron + Tc97

I installed the TENDL1.4 dataset and my environment variable is pointing to its folder(G4PARTICLEHPDATA, and I also tried with G4PROTONHPDATA).

I need to work with QGSP_BIC_AllHP since I am planning to simulate radioisotope production irradiating a sample of O18-enriched water with a proton beam and as far as I understood from literature QGSP_BIC_AllHP physic list give more accurate result for this specific application.

Geant4 Version: 11.2.2 (same issue with 10.5.1)
Operating System: Ubuntu 20.04
Compiler/Version: gcc version 9.4.0 (Ubuntu 9.4.0-1ubuntu1~20.04.2)
CMake Version: 3.16.3

Thank you for your help

2

@ElyonReed
What is your energy window ?

Your case perfectly fit with HP option, if you are interested in recoils then go for ALlHP case as ALLHP considers 0 threshold value for protons.

VRS

3

Thank you for your answer, I am working with proton energies ranging from 5/6 MeV up to 20 MeV.

I also tried to change the material in the Hadr03 example with 018-enriched water, in this case, using ALLHP, in the nuclear reaction list it appears "proton + O18 → " with nothing on the product side. (whereas using just HP the reaction is the complete one “proton + O18 → n + F18” )

4

Since your goal is the production of radioisotope using proton beams, meaning the activity measurements and emerging particles from isotope production, you may go through a very good example Activation provided in Geant4 examples.

Look at its Physics list and some of functions used to define radioisotope and its emission particles.

VRS

5

Dear @ElyonReed

Welcome to Geant4 forum!

Since you are new to Geant4, maybe it is worth to try first with the tool provided by IAEA for that purpose. They have a web-based calculator,

https://www-nds.iaea.org/relnsd/isotopia/isotopia.html

or you can download it from here

Learning how to use Geant4 for this goal may require more time than the IAEA tool. You can try out both and compare :slight_smile:

Best,
Alvaro

6

Thank you for this advice!

7

Thank you, I didn’t notice this example, I will surely study it.