BertiniEventTopologyProcess.cxx

#include "SimCore/PhotoNuclearModels/BertiniEventTopologyProcess.h"
namespace simcore {
 
void BertiniEventTopologyProcess::cleanupSecondaries() {
  int secondaries{theParticleChange.GetNumberOfSecondaries()};
  // Geant4 won't clean up this memory for us by default
  for (int i{0}; i < secondaries; ++i) {
    auto secondary{theParticleChange.GetSecondary(i)->GetParticle()};
    delete secondary;
  }
}
 
G4HadFinalState* BertiniEventTopologyProcess::ApplyYourself(
    const G4HadProjectile& projectile, G4Nucleus& targetNucleus) {
  int attempts{1};
  if (!acceptProjectile(projectile) || !acceptTarget(targetNucleus)) {
    // Bertini will handle the particle change on its own here
    return G4CascadeInterface::ApplyYourself(projectile, targetNucleus);
  }
 
  while (true) {
    theParticleChange.Clear();
    theParticleChange.SetStatusChange(stopAndKill);
    G4CascadeInterface::ApplyYourself(projectile, targetNucleus);
    if (acceptEvent()) {
      incrementEventWeight(attempts);
      return &theParticleChange;
    }
    attempts++;
    cleanupSecondaries();
  }
}
 
}  // namespace simcore