Simulation about Neutron (E>20MeV) incident bonner sphere with He3 Detector

Physics Lists
1

Hi there, everyone!

Recently, I was working on the topics about the Neutron beam (E>20MeV, up to 100MeV) incident on the bonner sphere with He3 Detector based on the example extended\hadronic\NeutronSource. The original intention of using this example for simulation is that it seems to be relatively complete for the physicslist of neutrons.

However, as we all know, the cross-section data of neutrons with energy higher than 20MeV is unknown. This example seems to have been simulated by extrapolating the database, which might have led to some unsatisfactory results. However, due to the lack of experimental data, I was also unable to conduct data comparisons. The key code of physicslists are as follows:

G4HadronElasticPhysics::ConstructProcess();
GetNeutronModel()->SetMinEnergy(19.5 * MeV);

G4HadronicProcess* process = GetNeutronProcess();
G4ParticleHPElastic* model1 = new G4ParticleHPElastic();
process->RegisterMe(model1);
process->AddDataSet(new G4ParticleHPElasticData());

if (fThermal) {
model1->SetMinEnergy(4 * eV);
G4ParticleHPThermalScattering* model2 = new G4ParticleHPThermalScattering();
process->RegisterMe(model2);
process->AddDataSet(new G4ParticleHPThermalScatteringData());
}

My question is as follows:

  1. I’m not sure if the above statement is correctly understood for this example?

  2. If I don’t want to use the extrapolation method, are there any better theoretical models to choose from?

  3. How should I reorganize my code?

Thanks a lot if you can solve my problem!!

Geant4 Version:11.3
Operating System:Windows
Compiler/Version:
CMake Version:3.31.3

1 Like