49 , fpReactionModel(nullptr)
61 double currentStepTime,
62 bool userStepTimeLimit)
71 if (currentStepTime == 0.)
73 userStepTimeLimit =
false;
77 separationDistance, userStepTimeLimit);
96 const G4double D1 = pMoleculeA->GetDiffusionCoefficient();
97 const G4double D2 = pMoleculeB->GetDiffusionCoefficient();
98 const G4double sqrD1 = D1 == 0. ? 0. : std::sqrt(D1);
99 const G4double sqrD2 = D2 == 0. ? 0. : std::sqrt(D2);
100 const G4double inv_numerator = 1./(sqrD1 + sqrD2);
104 for (
G4int j = 0; j < nbProducts; ++j)
106 auto pProduct =
new G4Molecule(pReactionData->GetProduct(j));
109 pProductTrack->SetTrackStatus(
fAlive);