How to disable full radioactive chain

[Emazzols]

Hello everyone,

I am generating generic ion by the standard way of replacing the geantino particle. Here my lines of code:

if(fParticleGun->GetParticleDefinition() == G4Geantino::Geantino())
{
G4double IonCharge = 0.*eplus;
G4double ExcitEnergy = 0.*keV;

G4ParticleDefinition *ion = G4IonTable::GetIonTable()->GetIon(Z, A, ExcitEnergy);
fParticleGun->SetParticleDefinition(ion);
fParticleGun->SetParticleCharge(IonCharge);
fParticleGun->SetParticleTime(0.0*ns);

}

I use a configuration file to define the Z and A variables.

In particular, I’m interested in Am241. I’d like to know how to stop the generation of the sons of Am241. I’m interested only in the decay of the 241Am (alpha around 5.4 MeV and gamma around 59 keV) and not in the decay of all the sons (from 237Np onward).

Thanks a lot for your help!

_Geant4 Version:_11.1.2
_Operating System:_Ubuntu 22.04.3 LTS on wsl
_Compiler/Version:_gcc 11.4.0
_CMake Version:_3.22.1

---

[Emazzols]

I got it! I put here the solution:

if(fParticleGun->GetParticleDefinition() == G4Geantino::Geantino())
{
G4double IonCharge = 0.eplus;
G4double ExcitEnergy = 0.keV;
G4ParticleDefinition ion = G4IonTable::GetIonTable()->GetIon(Z, A, ExcitEnergy);
fParticleGun->SetParticleDefinition(ion);
fParticleGun->SetParticleCharge(IonCharge);
fParticleGun->SetParticleTime(0.0ns);
if(full_chain==1){
auto radioactiveDecay = static_cast<G4RadioactiveDecay>(G4ProcessTable::GetProcessTable()->FindProcess(“RadioactiveDecay”, “GenericIon”));
if (radioactiveDecay) {
G4NucleusLimits nucleusLimits(A, A, Z, Z);
radioactiveDecay->SetNucleusLimits(nucleusLimits);
G4DeexPrecoParameters deex = G4NuclearLevelData::GetInstance()->GetParameters();
deex->SetCorrelatedGamma(true);
deex->SetStoreAllLevels(true);
// deex->SetInternalConversionFlag(true);
// deex->SetIsomerProduction(true);
} else {
G4cerr << “RadioactiveDecay process not found for GenericIon!” << G4endl;
}
}
}

[maire1]

Why not use the command /process/had/rdm/nucleusLimits ?
Here a macro for example rdecay01:
mazzola.mac.txt (307 Bytes)

[Emazzols]

I didn’t know these commands. I will try, thank you!