G4FermiSplitter Class Reference

#include <G4FermiSplitter.hh>

Static Public Member Functions
static G4double	DecayWeight (const G4FermiFragmentVector &split, G4FermiAtomicMass atomicMass, G4double totalEnergy)
static G4double	SplitFactor (const G4FermiFragmentVector &split, G4FermiAtomicMass atomicMass)
static G4double	KineticFactor (const G4FermiFragmentVector &split, G4double totalEnergy)
static void	GenerateSplits (G4FermiNucleiData nucleiData, std::vector< G4FermiFragmentVector > &splits)
static std::vector< G4FermiFragmentVector >	GenerateSplits (G4FermiNucleiData nucleiData)

Detailed Description

Definition at line 38 of file G4FermiSplitter.hh.

Member Function Documentation

DecayWeight()

G4double G4FermiSplitter::DecayWeight ( const G4FermiFragmentVector & split, G4FermiAtomicMass atomicMass, G4double totalEnergy )

static

Definition at line 175 of file G4FermiSplitter.cc.

{
  const auto kineticEnergy = KineticEnergy(split, totalEnergy);  // in MeV
 
  // Check that there is enough energy to produce K fragments
  if (kineticEnergy <= 0.) {
    return 0.;
  }
 
  const auto power = 3.0 * static_cast<G4double>(split.size() - 1) / 2.0 - 1.;
  const auto kineticFactor = std::pow(kineticEnergy, power);
 
  // Spin factor S_n
  const auto spinFactor = SpinFactor(split);
 
  // Calculate MassFactor
  const auto massFactor = MassFactor(split);
 
  // This is the constant (doesn't depend on energy) part
  const auto coef = ConstFactor(atomicMass, split.size());
 
  // Calculation of 1/gamma(3(k-1)/2)
  const auto gamma = GammaFactor(split.size());
 
  // Permutation Factor G_n
  const auto permutationFactor = ConfigurationFactor(split);
 
  return coef * kineticFactor * massFactor * spinFactor / (permutationFactor * gamma);
}

GenerateSplits() [1/2]

std::vector< G4FermiFragmentVector > G4FermiSplitter::GenerateSplits ( G4FermiNucleiData nucleiData )

static

Definition at line 274 of file G4FermiSplitter.cc.

{
  std::vector<G4FermiFragmentVector> splits;
  GenerateSplits(nucleiData, splits);
  return splits;
}

GenerateSplits() [2/2]

void G4FermiSplitter::GenerateSplits ( G4FermiNucleiData nucleiData, std::vector< G4FermiFragmentVector > & splits )

static

Definition at line 281 of file G4FermiSplitter.cc.

{
  ThrowOnInvalidInputs(nucleiData);
 
  splits.reserve(ExpectedSplitSize);
 
  // let's split nucleus into 2, ..., A fragments
  const auto maxFragmentsCount = static_cast<std::uint32_t>(nucleiData.atomicMass);
 
  for (std::uint32_t fragmentCount = 2; fragmentCount <= maxFragmentsCount; ++fragmentCount) {
    // Form all possible partition by combination of A partitions and Z partitions (Z partitions
    // include null parts)
    for (auto& massPartition : G4integerPartition(nucleiData.atomicMass, fragmentCount, 1)) {
      for (auto& chargePartition : G4integerPartition(nucleiData.chargeNumber, fragmentCount, 0)) {
        // Some splits are invalid, some nuclei doesn't exist
        if (auto partitionSplits = PossibleSplits(massPartition, chargePartition);
            !partitionSplits.empty()) {
          splits.insert(splits.end(), std::make_move_iterator(partitionSplits.begin()),
                        std::make_move_iterator(partitionSplits.end()));
        }
      }
    }
  }
}

Referenced by GenerateSplits(), and G4FermiBreakUpAN::Initialise().

KineticFactor()

G4double G4FermiSplitter::KineticFactor ( const G4FermiFragmentVector & split, G4double totalEnergy )

static

SplitFactor()

G4double G4FermiSplitter::SplitFactor ( const G4FermiFragmentVector & split, G4FermiAtomicMass atomicMass )

static

---

The documentation for this class was generated from the following files:

• G4FermiSplitter.hh
• G4FermiSplitter.cc