Gate icon indicating copy to clipboard operation
Gate copied to clipboard

[ENH] G4Radioactivation process is not implementing decay chain time biasing

Open BishopWolf opened this issue 5 years ago • 0 comments

The old G4RadioactiveDecay process that was implemented in Gate was meant for simple decay chains. However it has some big problems such as not controlling the time biasing of long half lived daughters

The new G4Radioactivation process solves this issue but it must be implemented as in this example. Specially in this file. The relevant lines are:

  // mandatory for G4NuclideTable
  //
  G4NuclideTable::GetInstance()->SetThresholdOfHalfLife(0.1*picosecond);
  G4NuclideTable::GetInstance()->SetLevelTolerance(1.0*eV);

  //read new PhotonEvaporation data set 
  //
  G4DeexPrecoParameters* deex = 
    G4NuclearLevelData::GetInstance()->GetParameters();
  deex->SetCorrelatedGamma(false);
  deex->SetStoreAllLevels(true);
  deex->SetMaxLifeTime(G4NuclideTable::GetInstance()->GetThresholdOfHalfLife()
                /std::log(2.));

As you can see the initial evaporation process defined in G4RadioactiveDecay is not sufficient, it also needs control on the deexitation process to remove very long half lived daughters.

This shall be implemented in Gate to remove some of the emission spectra issues that are mentioned.

BishopWolf avatar Nov 16 '20 17:11 BishopWolf