Hello Geant4 experts:
I am simulating the process of Boron Neutron Capture Therapy (BNCT) in Geant4. Are there any suggestion of the physics list in this work?
In this work I need to simulate the process of nCapture in B10(n,a)Li7 channel, track the dose of 2 heavy fragments, and also track the 478keV gamma from Li7* deexcited to simulate the BNCT-SPECT. Beside B10 nCapture, H(n,gamma)D is also a necessary channel to be simulated, and other light nuclear nCapture as well.
Currently my physics list is :
// Bosons
G4BosonConstructor bs;
bs.ConstructParticle();
// Leptons
G4LeptonConstructor lp;
lp.ConstructParticle();
// Mesons
G4MesonConstructor mc;
mc.ConstructParticle();
// Baryons
G4BaryonConstructor br;
br.ConstructParticle();
// Ions
G4IonConstructor in;
in.ConstructParticle();
// Short Lived
G4ShortLivedConstructor shortLived;
shortLived.ConstructParticle();
// physics process
RegisterPhysics(new G4EmStandardPhysics);
RegisterPhysics(new G4HadronElasticPhysicsHP(verb));
RegisterPhysics(new G4HadronPhysicsQGSP_BIC_HP(verb));
RegisterPhysics( new G4IonElasticPhysics(verb));
RegisterPhysics( new G4IonPhysicsXS(verb));
RegisterPhysics(new G4StoppingPhysics(verb));
RegisterPhysics(new G4RadioactiveDecayPhysics(verb));
RegisterPhysics(new G4IonBinaryCascadePhysics(verb));
RegisterPhysics(new G4DecayPhysics(verb));
// RegisterPhysics(new G4EmExtraPhysics);
// RegisterPhysics(new G4NeutronTrackingCut);
With the above code I simulated the dose distribution of 25 keV monoenergetic neutron induced, we can obviously find that the dose at entrance is higher than that in the target volume (see this figure from top to bottom which is the depth increasing direction, and target is at the middle depth):
For the physics list code above, are there any suggestions, about which one is necessary and which is not, and any other physics is needed to append? Thanks in advance.
_Geant4 Version:_11.0.1
_Operating System:_windows subsystem for linux (WSL2), ubuntu 22.04 LTS
_Compiler/Version:_gcc version 11.4.0
_CMake Version:_3.22.1