Please fill out the following information to help in answering your question, and also see tips for posting code snippets. If you don’t provide this information it will take more time to help with your problem!
Geant4 Version: 11.01
Operating System:
Compiler/Version:
CMake Version:
I used “TS_Beryllium_Metal”, “TS_H_of_Water”, and “TS_H_of_Polyethylene” in materials.
Additionally, I used G4ParticleHPThermalScattering() in physics ( < 4*eV neutrons).
I conducted simulation for generating thermal neutron with accelerator, beryllium target and HDPE modulator.
Finally, I got significantly unnatural results around 4*eV for neutron energy spectrum.
Is there any known error in TS material or physics?
If I use only G4ParticleHPThermalScattering() without any TS materials, It is meaningful?
It would be so appreciated if you could help me.
Especially, Polyethylene shows the most strange results.
These are the codes in DetectorConstruction.cc and HadronElasticPhysics, respectively.
//DetectorConstruction.cc
G4Element* H_HDPE = new G4Element(“TS_H_of_Polyethylene” ,“H” , 1., 1.0079g/mole);
G4Element C_HDPE = new G4Element(“C_of_HDPE” ,“C” , 6., 12.01g/mole);
HDPE =
new G4Material(“HDPE”, 0.95g/cm3, ncomponents=2,
kStateSolid, 293.15kelvin, 1atmosphere);
HDPE->AddElement(H_HDPE, 0.143711);
HDPE->AddElement(C_HDPE, 0.856289);
//HadronElasticPhysicsHP.cc
void HadronElasticPhysicsHP::ConstructProcess()
{
G4HadronElasticPhysics::ConstructProcess();
GetNeutronModel()->SetMinEnergy(19.5 * MeV);
G4HadronicProcess* process = GetNeutronProcess();
G4ParticleHPElastic* model1 = new G4ParticleHPElastic();
process->RegisterMe(model1);
process->AddDataSet(new G4ParticleHPElasticData());
// if (fThermal) {
if(1){
G4ParticleHPThermalScattering* model2 = new G4ParticleHPThermalScattering();
process->RegisterMe(model2);
process->AddDataSet(new G4ParticleHPThermalScatteringData());
// model1->SetMinEnergy(4. * eV);
}
}
I’m suspecting that no physics is applying to C in HDPE, because it is not TS material.