Branching ratio biasing for radioactive decay

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


1 Like