Hello,
_Geant4 Version:_11.1.2
_Operating System:_Ubuntu 22.04.4 LTS
_Compiler/Version:_11.4.0
_CMake Version:_3.22.1
Hello
I’m trying to use biasing to enhance the spontaneous fission process. According to slide 18 (https://www.slac.stanford.edu/xorg/geant4/SLACTutorial14/EventBiasing.pdf) I need to use G4Radioactivation and G4GenericIon, and construct a process in the physics list, which I have attempted to do:
#include “decayphysics.hh”
MyPhysicsList::MyPhysicsList()
:G4VModularPhysicsList()
{
RegisterPhysics (new G4DecayPhysics());
RegisterPhysics (new G4RadioactiveDecayPhysics());
}
MyPhysicsList::~MyPhysicsList()
{}
void MyPhysicsList::ConstructParticle()
{
G4GenericIon::GenericIon();
}
void MyPhysicsList::ConstructProcess()
{
G4PhysicsListHelper::GetPhysicsListHelper()->
RegisterProcess(new G4Radioactivation(), G4GenericIon::GenericIon());
}
I’m not sure if I’ve done this correctly as the model gets stuck. I’ve had a look at example rdecay02 but I can’t make sense of what’s going in the PhysicsList.cc and BiasedRDPhysics.cc. It seems like there’s some EM physics in here but I don’t need this for my application. I’ve also tried using generic biasing but this doesn’t work either.
Any help would be much appreciated.
Thanks,
Emma