#include <G4CompetitiveFission.hh>

Inheritance diagram for G4CompetitiveFission:

Public Member Functions
virtual G4bool	BreakUpChain (G4FragmentVector theResult, G4Fragment theNucleus)

G4FragmentVector *	BreakUpFragment (G4Fragment *theNucleus)

virtual void	Dump () const

G4Fragment *	EmittedFragment (G4Fragment *theNucleus) override

	G4CompetitiveFission ()

G4double	GetEmissionProbability (G4Fragment *theNucleus) override

G4double	GetFissionBarrier (void) const

G4double	GetLevelDensityParameter (void) const

virtual G4double	GetLifeTime (G4Fragment *theNucleus)

G4double	GetMaximalKineticEnergy (void) const

virtual void	Initialise ()

virtual void	RDMForced (G4bool)

void	SetEmissionStrategy (G4VEmissionProbability *aFissionProb)

void	SetFissionBarrier (G4VFissionBarrier *aBarrier)

virtual void	SetICM (G4bool)

void	SetLevelDensityParameter (G4VLevelDensityParameter *aLevelDensity)

void	SetOPTxs (G4int opt)

void	UseSICB (G4bool use)

	~G4CompetitiveFission () override

Protected Attributes
G4int	OPTxs

G4bool	useSICB

Private Member Functions
G4double	AsymmetricRatio (G4int A, G4double A11) const

G4int	FissionAtomicNumber (G4int A)

G4int	FissionCharge (G4int A, G4int Z, G4double Af)

G4double	FissionKineticEnergy (G4int A, G4int Z, G4int Af1, G4int Zf1, G4int Af2, G4int Zf2, G4double U, G4double Tmax)

	G4CompetitiveFission (const G4CompetitiveFission &right)

G4double	LocalExp (G4double x) const

G4double	MassDistribution (G4double x, G4int A)

G4bool	operator!= (const G4CompetitiveFission &right) const

const G4CompetitiveFission &	operator= (const G4CompetitiveFission &right)

G4bool	operator== (const G4CompetitiveFission &right) const

G4double	Ratio (G4double A, G4double A11, G4double B1, G4double A00) const

G4double	SymmetricRatio (G4int A, G4double A11) const

Detailed Description

Definition at line 46 of file G4CompetitiveFission.hh.

Constructor & Destructor Documentation

◆ G4CompetitiveFission() [1/2]

G4CompetitiveFission::G4CompetitiveFission ( )

explicit

Definition at line 50 of file G4CompetitiveFission.cc.

                                           : G4VEvaporationChannel("fission"), theSecID(-1)
{
  theFissionBarrierPtr = new G4FissionBarrier;
  myOwnFissionBarrier = true;
 
  theFissionProbabilityPtr = new G4FissionProbability;
  myOwnFissionProbability = true;
  
  theLevelDensityPtr = new G4FissionLevelDensityParameter;
  myOwnLevelDensity = true;
 
  maxKineticEnergy = fissionBarrier = fissionProbability = 0.0;
  pairingCorrection = G4NuclearLevelData::GetInstance()->GetPairingCorrection();
 
  theSecID = G4PhysicsModelCatalog::GetModelID("model_G4CompetitiveFission");
}

References fissionBarrier, fissionProbability, G4NuclearLevelData::GetInstance(), G4PhysicsModelCatalog::GetModelID(), G4NuclearLevelData::GetPairingCorrection(), maxKineticEnergy, myOwnFissionBarrier, myOwnFissionProbability, myOwnLevelDensity, pairingCorrection, theFissionBarrierPtr, theFissionProbabilityPtr, theLevelDensityPtr, and theSecID.

◆ ~G4CompetitiveFission()

G4CompetitiveFission::~G4CompetitiveFission ( )

override

Definition at line 67 of file G4CompetitiveFission.cc.

{
  if (myOwnFissionBarrier) delete theFissionBarrierPtr;
  if (myOwnFissionProbability) delete theFissionProbabilityPtr;
  if (myOwnLevelDensity) delete theLevelDensityPtr;
}

References myOwnFissionBarrier, myOwnFissionProbability, myOwnLevelDensity, theFissionBarrierPtr, theFissionProbabilityPtr, and theLevelDensityPtr.

◆ G4CompetitiveFission() [2/2]

G4CompetitiveFission::G4CompetitiveFission ( const G4CompetitiveFission & right )

private

Member Function Documentation

◆ AsymmetricRatio()

G4double G4CompetitiveFission::AsymmetricRatio	(	G4int	A,
		G4double	A11
	)		const

inlineprivate

Definition at line 149 of file G4CompetitiveFission.hh.

{
  return Ratio(G4double(A),A11,23.5,134.0);
}

References A, A11, and Ratio().

Referenced by FissionKineticEnergy().

◆ BreakUpChain()

G4bool G4VEvaporationChannel::BreakUpChain	(	G4FragmentVector *	theResult,
		G4Fragment *	theNucleus
	)

virtualinherited

Reimplemented in G4UnstableFragmentBreakUp, and G4PhotonEvaporation.

Definition at line 66 of file G4VEvaporationChannel.cc.

{
  return false;
}

Referenced by G4VEvaporationChannel::BreakUpFragment().

◆ BreakUpFragment()

G4FragmentVector * G4VEvaporationChannel::BreakUpFragment ( G4Fragment * theNucleus )

inlineinherited

Definition at line 106 of file G4VEvaporationChannel.hh.

{
  G4FragmentVector* results = new G4FragmentVector();
  BreakUpChain(results, theNucleus);
  return results;
}

References G4VEvaporationChannel::BreakUpChain().

Referenced by G4NeutronRadCapture::ApplyYourself().

◆ Dump()

void G4VEvaporationChannel::Dump ( ) const

virtualinherited

Reimplemented in G4GEMChannel, and G4GEMChannelVI.

Definition at line 71 of file G4VEvaporationChannel.cc.

72{}

◆ EmittedFragment()

G4Fragment * G4CompetitiveFission::EmittedFragment ( G4Fragment * theNucleus )

overridevirtual

Reimplemented from G4VEvaporationChannel.

Definition at line 95 of file G4CompetitiveFission.cc.

{
  G4Fragment * Fragment1 = nullptr; 
  // Nucleus data
  // Atomic number of nucleus
  G4int A = theNucleus->GetA_asInt();
  // Charge of nucleus
  G4int Z = theNucleus->GetZ_asInt();
  //   Excitation energy (in MeV)
  G4double U = theNucleus->GetExcitationEnergy();
  G4double pcorr = pairingCorrection->GetFissionPairingCorrection(A,Z);
  if (U <= pcorr) { return Fragment1; }
 
  // Atomic Mass of Nucleus (in MeV)
  G4double M = theNucleus->GetGroundStateMass();
 
  // Nucleus Momentum
  G4LorentzVector theNucleusMomentum = theNucleus->GetMomentum();
 
  // Calculate fission parameters
  theParam.DefineParameters(A, Z, U-pcorr, fissionBarrier);
  
  // First fragment
  G4int A1 = 0;
  G4int Z1 = 0;
  G4double M1 = 0.0;
 
  // Second fragment
  G4int A2 = 0;
  G4int Z2 = 0;
  G4double M2 = 0.0;
 
  G4double FragmentsExcitationEnergy = 0.0;
  G4double FragmentsKineticEnergy = 0.0;
 
  G4int Trials = 0;
  do {
 
    // First fragment 
    A1 = FissionAtomicNumber(A);
    Z1 = FissionCharge(A, Z, A1);
    M1 = G4NucleiProperties::GetNuclearMass(A1, Z1);
 
    // Second Fragment
    A2 = A - A1;
    Z2 = Z - Z1;
    if (A2 < 1 || Z2 < 0 || Z2 > A2) {
      FragmentsExcitationEnergy = -1.0;
      continue;
    }
    M2 = G4NucleiProperties::GetNuclearMass(A2, Z2);
    // Maximal Kinetic Energy (available energy for fragments)
    G4double Tmax = M + U - M1 - M2 - pcorr;
 
    // Check that fragment masses are less or equal than total energy
    if (Tmax < 0.0) {
      FragmentsExcitationEnergy = -1.0;
      continue;
    }
 
    FragmentsKineticEnergy = FissionKineticEnergy( A , Z,
                           A1, Z1,
                           A2, Z2,
                           U , Tmax);
    
    // Excitation Energy
    // FragmentsExcitationEnergy = Tmax - FragmentsKineticEnergy;
    // JMQ 04/03/09 BUG FIXED: in order to fulfill energy conservation the
    // fragments carry the fission pairing energy in form of 
    // excitation energy
 
    FragmentsExcitationEnergy = 
      Tmax - FragmentsKineticEnergy + pcorr;
 
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  } while (FragmentsExcitationEnergy < 0.0 && ++Trials < 100);
    
  if (FragmentsExcitationEnergy <= 0.0) { 
    throw G4HadronicException(__FILE__, __LINE__, 
      "G4CompetitiveFission::BreakItUp: Excitation energy for fragments < 0.0!");
  }
 
  // Fragment 1
  M1 += FragmentsExcitationEnergy * A1/static_cast<G4double>(A);
  // Fragment 2
  M2 += FragmentsExcitationEnergy * A2/static_cast<G4double>(A);
  // primary
  M += U;
 
  G4double etot1 = ((M - M2)*(M + M2) + M1*M1)/(2*M);
  G4ParticleMomentum Momentum1 = 
    std::sqrt((etot1 - M1)*(etot1+M1))*G4RandomDirection();
  G4LorentzVector FourMomentum1(Momentum1, etot1);
  FourMomentum1.boost(theNucleusMomentum.boostVector());
    
  // Create Fragments
  Fragment1 = new G4Fragment( A1, Z1, FourMomentum1);
  if (Fragment1 != nullptr) { Fragment1->SetCreatorModelID(theSecID); }
  theNucleusMomentum -= FourMomentum1;
  theNucleus->SetZandA_asInt(Z2, A2);
  theNucleus->SetMomentum(theNucleusMomentum);
  theNucleus->SetCreatorModelID(theSecID);
  return Fragment1;
}

References A, CLHEP::HepLorentzVector::boost(), CLHEP::HepLorentzVector::boostVector(), G4FissionParameters::DefineParameters(), FissionAtomicNumber(), fissionBarrier, FissionCharge(), FissionKineticEnergy(), G4RandomDirection(), G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4PairingCorrection::GetFissionPairingCorrection(), G4Fragment::GetGroundStateMass(), G4Fragment::GetMomentum(), G4NucleiProperties::GetNuclearMass(), G4Fragment::GetZ_asInt(), M, pairingCorrection, G4Fragment::SetCreatorModelID(), G4Fragment::SetMomentum(), G4Fragment::SetZandA_asInt(), theParam, theSecID, Z, and anonymous_namespace{paraMaker.cc}::Z1.

◆ FissionAtomicNumber()

G4int G4CompetitiveFission::FissionAtomicNumber ( G4int A )

private

Definition at line 201 of file G4CompetitiveFission.cc.

{
 
  // For Simplicity reading code
  G4int A1 = theParam.GetA1();
  G4int A2 = theParam.GetA2();
  G4double As = theParam.GetAs();
  G4double Sigma2 = theParam.GetSigma2();
  G4double SigmaS = theParam.GetSigmaS();
  G4double w = theParam.GetW();
  
  G4double C2A = A2 + 3.72*Sigma2;
  G4double C2S = As + 3.72*SigmaS;
  
  G4double C2 = 0.0;
  if (w > 1000.0 )    { C2 = C2S; }
  else if (w < 0.001) { C2 = C2A; }
  else                { C2 =  std::max(C2A,C2S); }
 
  G4double C1 = A-C2;
  if (C1 < 30.0) {
    C2 = A-30.0;
    C1 = 30.0;
  }
 
  G4double Am1 = (As + A1)*0.5;
  G4double Am2 = (A1 + A2)*0.5;
 
  // Get Mass distributions as sum of symmetric and asymmetric Gasussians
  G4double Mass1 = MassDistribution(As,A); 
  G4double Mass2 = MassDistribution(Am1,A); 
  G4double Mass3 = MassDistribution(G4double(A1),A); 
  G4double Mass4 = MassDistribution(Am2,A); 
  G4double Mass5 = MassDistribution(G4double(A2),A); 
  // get maximal value among Mass1,...,Mass5
  G4double MassMax = Mass1;
  if (Mass2 > MassMax) { MassMax = Mass2; }
  if (Mass3 > MassMax) { MassMax = Mass3; }
  if (Mass4 > MassMax) { MassMax = Mass4; }
  if (Mass5 > MassMax) { MassMax = Mass5; }
 
  // Sample a fragment mass number, which lies between C1 and C2
  G4double xm;
  G4double Pm;
  do {
    xm = C1+G4UniformRand()*(C2-C1);
    Pm = MassDistribution(xm,A); 
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  } while (MassMax*G4UniformRand() > Pm);
  G4int ires = G4lrint(xm);
 
  return ires;
}

References A, C1, C2, G4lrint(), G4UniformRand, G4FissionParameters::GetA1(), G4FissionParameters::GetA2(), G4FissionParameters::GetAs(), G4FissionParameters::GetSigma2(), G4FissionParameters::GetSigmaS(), G4FissionParameters::GetW(), MassDistribution(), G4INCL::Math::max(), and theParam.

Referenced by EmittedFragment().

◆ FissionCharge()

G4int G4CompetitiveFission::FissionCharge	(	G4int	A,
		G4int	Z,
		G4double	Af
	)

private

Definition at line 279 of file G4CompetitiveFission.cc.

{
  static const G4double sigma = 0.6;
  G4double DeltaZ = 0.0;
  if (Af >= 134.0)          { DeltaZ = -0.45; }  
  else if (Af <= (A-134.0)) { DeltaZ = 0.45; }
  else                      { DeltaZ = -0.45*(Af-A*0.5)/(134.0-A*0.5); }
 
  G4double Zmean = (Af/A)*Z + DeltaZ;
 
  G4double theZ;
  do {
    theZ = G4RandGauss::shoot(Zmean,sigma);
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  } while (theZ  < 1.0 || theZ > (Z-1.0) || theZ > Af);
  
  return G4lrint(theZ);
}

References A, G4lrint(), G4INCL::DeJongSpin::shoot(), and Z.

Referenced by EmittedFragment().

◆ FissionKineticEnergy()

G4double G4CompetitiveFission::FissionKineticEnergy	(	G4int	A,
		G4int	Z,
		G4int	Af1,
		G4int	Zf1,
		G4int	Af2,
		G4int	Zf2,
		G4double	U,
		G4double	Tmax
	)

private

Definition at line 300 of file G4CompetitiveFission.cc.

{
  // Find maximal value of A for fragments
  G4int AfMax = std::max(Af1,Af2);
 
  // Weights for symmetric and asymmetric components
  G4double Pas = 0.0;
  if (theParam.GetW() <= 1000) { 
    G4double x1 = (AfMax-theParam.GetA1())/theParam.GetSigma1();
    G4double x2 = (AfMax-theParam.GetA2())/theParam.GetSigma2();
    Pas = 0.5*LocalExp(x1) + LocalExp(x2);
  }
 
  G4double Ps = 0.0;
  if (theParam.GetW() >= 0.001) {
    G4double xs = (AfMax-theParam.GetAs())/theParam.GetSigmaS();
    Ps = theParam.GetW()*LocalExp(xs);
  }
  G4double Psy = (Pas + Ps > 0.0) ? Ps/(Pas+Ps) : 0.5;
 
  // Fission fractions Xsy and Xas formed in symmetric and asymmetric modes
  G4double PPas = theParam.GetSigma1() + 2.0 * theParam.GetSigma2();
  G4double PPsy = theParam.GetW() * theParam.GetSigmaS();
  G4double Xas = (PPas + PPsy > 0.0) ? PPas/(PPas+PPsy) : 0.5;
  G4double Xsy = 1.0 - Xas;
 
  // Average kinetic energy for symmetric and asymmetric components
  G4double Eaverage = (0.1071*(Z*Z)/G4Pow::GetInstance()->Z13(A) + 22.2)*CLHEP::MeV;
 
  // Compute maximal average kinetic energy of fragments and Energy Dispersion 
  G4double TaverageAfMax;
  G4double ESigma = 10*CLHEP::MeV;
  // Select randomly fission mode (symmetric or asymmetric)
  if (G4UniformRand() > Psy) { // Asymmetric Mode
    G4double A11 = theParam.GetA1()-0.7979*theParam.GetSigma1();
    G4double A12 = theParam.GetA1()+0.7979*theParam.GetSigma1();
    G4double A21 = theParam.GetA2()-0.7979*theParam.GetSigma2();
    G4double A22 = theParam.GetA2()+0.7979*theParam.GetSigma2();
    // scale factor
    G4double ScaleFactor = 0.5*theParam.GetSigma1()*
      (AsymmetricRatio(A,A11)+AsymmetricRatio(A,A12))+
      theParam.GetSigma2()*(AsymmetricRatio(A,A21)+AsymmetricRatio(A,A22));
    // Compute average kinetic energy for fragment with AfMax
    TaverageAfMax = (Eaverage + 12.5 * Xsy) * (PPas/ScaleFactor) * 
      AsymmetricRatio(A,G4double(AfMax));
 
  } else { // Symmetric Mode
    G4double As0 = theParam.GetAs() + 0.7979*theParam.GetSigmaS();
    // Compute average kinetic energy for fragment with AfMax
    TaverageAfMax = (Eaverage - 12.5*CLHEP::MeV*Xas)
      *SymmetricRatio(A, G4double(AfMax))/SymmetricRatio(A, As0);
    ESigma = 8.0*CLHEP::MeV;
  }
 
  // Select randomly, in accordance with Gaussian distribution, 
  // fragment kinetic energy
  G4double KineticEnergy;
  G4int i = 0;
  do {
    KineticEnergy = G4RandGauss::shoot(TaverageAfMax, ESigma);
    if (++i > 100) return Eaverage;
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
  } while (KineticEnergy < Eaverage-3.72*ESigma || 
       KineticEnergy > Eaverage+3.72*ESigma ||
       KineticEnergy > Tmax);
  
  return KineticEnergy;
}

References A, A11, A12, A21, A22, AsymmetricRatio(), G4UniformRand, G4FissionParameters::GetA1(), G4FissionParameters::GetA2(), G4FissionParameters::GetAs(), G4Pow::GetInstance(), G4FissionParameters::GetSigma1(), G4FissionParameters::GetSigma2(), G4FissionParameters::GetSigmaS(), G4FissionParameters::GetW(), LocalExp(), G4INCL::Math::max(), CLHEP::MeV, G4INCL::DeJongSpin::shoot(), SymmetricRatio(), theParam, Z, and G4Pow::Z13().

Referenced by EmittedFragment().

◆ GetEmissionProbability()

G4double G4CompetitiveFission::GetEmissionProbability ( G4Fragment * theNucleus )

overridevirtual

Implements G4VEvaporationChannel.

Definition at line 74 of file G4CompetitiveFission.cc.

{
  G4int Z = fragment->GetZ_asInt();
  G4int A = fragment->GetA_asInt();
  fissionProbability = 0.0;
  // Saddle point excitation energy ---> A = 65
  if (A >= 65 && Z > 16) {
    G4double exEnergy = fragment->GetExcitationEnergy() - 
      pairingCorrection->GetFissionPairingCorrection(A, Z);
  
    if (exEnergy > 0.0) {
      fissionBarrier = theFissionBarrierPtr->FissionBarrier(A, Z, exEnergy);
      maxKineticEnergy = exEnergy - fissionBarrier;
      fissionProbability = 
    theFissionProbabilityPtr->EmissionProbability(*fragment,
                              maxKineticEnergy);
    }
  }
  return fissionProbability;
}

References A, G4VEmissionProbability::EmissionProbability(), fissionBarrier, G4VFissionBarrier::FissionBarrier(), fissionProbability, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4PairingCorrection::GetFissionPairingCorrection(), G4Fragment::GetZ_asInt(), maxKineticEnergy, pairingCorrection, theFissionBarrierPtr, theFissionProbabilityPtr, and Z.

◆ GetFissionBarrier()

G4double G4CompetitiveFission::GetFissionBarrier ( void ) const

inline

Definition at line 123 of file G4CompetitiveFission.hh.

{ 
  return fissionBarrier; 
}

References fissionBarrier.

◆ GetLevelDensityParameter()

G4double G4CompetitiveFission::GetLevelDensityParameter ( void ) const

inline

◆ GetLifeTime()

G4double G4VEvaporationChannel::GetLifeTime ( G4Fragment * theNucleus )

virtualinherited

Definition at line 50 of file G4VEvaporationChannel.cc.

{
  return 0.0;
}

◆ GetMaximalKineticEnergy()

G4double G4CompetitiveFission::GetMaximalKineticEnergy ( void ) const

inline

Definition at line 128 of file G4CompetitiveFission.hh.

{ 
  return maxKineticEnergy; 
}

References maxKineticEnergy.

◆ Initialise()

void G4VEvaporationChannel::Initialise ( )

virtualinherited

Reimplemented in G4PhotonEvaporation, and G4EvaporationChannel.

Definition at line 47 of file G4VEvaporationChannel.cc.

48{}

Referenced by G4EvaporationChannel::Initialise(), and G4NeutronRadCapture::InitialiseModel().

◆ LocalExp()

G4double G4CompetitiveFission::LocalExp ( G4double x ) const

inlineprivate

Definition at line 161 of file G4CompetitiveFission.hh.

{
  return (std::abs(x) < 8.) ? G4Exp(-0.5*x*x) : 0.0;
}

References G4Exp().

Referenced by FissionKineticEnergy(), and MassDistribution().

◆ MassDistribution()

G4double G4CompetitiveFission::MassDistribution	(	G4double	x,
		G4int	A
	)

private

Definition at line 257 of file G4CompetitiveFission.cc.

{
  G4double y0 = (x-theParam.GetAs())/theParam.GetSigmaS();
  G4double Xsym = LocalExp(y0);
 
  G4double y1 = (x - theParam.GetA1())/theParam.GetSigma1();
  G4double y2 = (x - theParam.GetA2())/theParam.GetSigma2();
  G4double z1 = (x - A + theParam.GetA1())/theParam.GetSigma1();
  G4double z2 = (x - A + theParam.GetA2())/theParam.GetSigma2();
  G4double Xasym = LocalExp(y1) + LocalExp(y2) 
    + 0.5*(LocalExp(z1) + LocalExp(z2));
 
  G4double res;
  G4double w = theParam.GetW();
  if (w > 1000)       { res = Xsym; }
  else if (w < 0.001) { res = Xasym; }
  else                { res = w*Xsym+Xasym; }
  return res;
}

References A, G4FissionParameters::GetA1(), G4FissionParameters::GetA2(), G4FissionParameters::GetAs(), G4FissionParameters::GetSigma1(), G4FissionParameters::GetSigma2(), G4FissionParameters::GetSigmaS(), G4FissionParameters::GetW(), LocalExp(), and theParam.

Referenced by FissionAtomicNumber().

◆ operator!=()

G4bool G4CompetitiveFission::operator!= ( const G4CompetitiveFission & right ) const

private

◆ operator=()

const G4CompetitiveFission & G4CompetitiveFission::operator= ( const G4CompetitiveFission & right )

private

◆ operator==()

G4bool G4CompetitiveFission::operator== ( const G4CompetitiveFission & right ) const

private

◆ Ratio()

G4double G4CompetitiveFission::Ratio	(	G4double	A,
		G4double	A11,
		G4double	B1,
		G4double	A00
	)		const

inlineprivate

Definition at line 134 of file G4CompetitiveFission.hh.

{
  G4double res;
  if (A11 >= A*0.5 && A11 <= (A00+10.0)) {
    G4double x = (A11-A00)/A;
    res = 1.0 - B1*x*x;
  } else {
    G4double x = 10.0/A;
    res = 1.0 - B1*x*x - 2.0*x*B1*(A11-A00-10.0)/A;
  }
  return res;
}

References A, A00, and A11.

Referenced by AsymmetricRatio(), and SymmetricRatio().

◆ RDMForced()

void G4VEvaporationChannel::RDMForced ( G4bool )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 58 of file G4VEvaporationChannel.cc.

59{}

◆ SetEmissionStrategy()

void G4CompetitiveFission::SetEmissionStrategy ( G4VEmissionProbability * aFissionProb )

Definition at line 381 of file G4CompetitiveFission.cc.

{
  if (myOwnFissionProbability) delete theFissionProbabilityPtr;
  theFissionProbabilityPtr = aFissionProb;
  myOwnFissionProbability = false;
}

References myOwnFissionProbability, and theFissionProbabilityPtr.

Referenced by G4INCLXXInterface::G4INCLXXInterface().

◆ SetFissionBarrier()

void G4CompetitiveFission::SetFissionBarrier ( G4VFissionBarrier * aBarrier )

Definition at line 373 of file G4CompetitiveFission.cc.

{
  if (myOwnFissionBarrier) delete theFissionBarrierPtr;
  theFissionBarrierPtr = aBarrier;
  myOwnFissionBarrier = false;
}

References myOwnFissionBarrier, and theFissionBarrierPtr.

◆ SetICM()

void G4VEvaporationChannel::SetICM ( G4bool )

virtualinherited

Reimplemented in G4PhotonEvaporation.

Definition at line 55 of file G4VEvaporationChannel.cc.

56{}

Referenced by G4NeutronRadCapture::InitialiseModel().

◆ SetLevelDensityParameter()

void G4CompetitiveFission::SetLevelDensityParameter ( G4VLevelDensityParameter * aLevelDensity )

Definition at line 389 of file G4CompetitiveFission.cc.

{ 
  if (myOwnLevelDensity) delete theLevelDensityPtr;
  theLevelDensityPtr = aLevelDensity;
  myOwnLevelDensity = false;
}

References myOwnLevelDensity, and theLevelDensityPtr.

Referenced by G4INCLXXInterface::G4INCLXXInterface().

◆ SetOPTxs()

void G4VEvaporationChannel::SetOPTxs ( G4int opt )

inlineinherited

Definition at line 113 of file G4VEvaporationChannel.hh.

114{}

◆ SymmetricRatio()

G4double G4CompetitiveFission::SymmetricRatio	(	G4int	A,
		G4double	A11
	)		const

inlineprivate

Definition at line 155 of file G4CompetitiveFission.hh.

{
  G4double A0 = G4double(A);
  return Ratio(A0,A11,5.32,A0*0.5);
}

References A, A11, and Ratio().

Referenced by FissionKineticEnergy().

◆ UseSICB()

void G4VEvaporationChannel::UseSICB ( G4bool use )

inlineinherited

Definition at line 116 of file G4VEvaporationChannel.hh.

117{}

Field Documentation

◆ fissionBarrier

G4double G4CompetitiveFission::fissionBarrier

private

Definition at line 101 of file G4CompetitiveFission.hh.

Referenced by EmittedFragment(), G4CompetitiveFission(), GetEmissionProbability(), and GetFissionBarrier().

◆ fissionProbability

G4double G4CompetitiveFission::fissionProbability

private

Definition at line 102 of file G4CompetitiveFission.hh.

Referenced by G4CompetitiveFission(), and GetEmissionProbability().

◆ maxKineticEnergy

G4double G4CompetitiveFission::maxKineticEnergy

private

Definition at line 100 of file G4CompetitiveFission.hh.

Referenced by G4CompetitiveFission(), GetEmissionProbability(), and GetMaximalKineticEnergy().

◆ myOwnFissionBarrier

G4bool G4CompetitiveFission::myOwnFissionBarrier

private

Definition at line 115 of file G4CompetitiveFission.hh.

Referenced by G4CompetitiveFission(), SetFissionBarrier(), and ~G4CompetitiveFission().

◆ myOwnFissionProbability

G4bool G4CompetitiveFission::myOwnFissionProbability

private

Definition at line 114 of file G4CompetitiveFission.hh.

Referenced by G4CompetitiveFission(), SetEmissionStrategy(), and ~G4CompetitiveFission().

◆ myOwnLevelDensity

G4bool G4CompetitiveFission::myOwnLevelDensity

private

Definition at line 116 of file G4CompetitiveFission.hh.

Referenced by G4CompetitiveFission(), SetLevelDensityParameter(), and ~G4CompetitiveFission().

◆ OPTxs

G4int G4VEvaporationChannel::OPTxs

protectedinherited

Definition at line 93 of file G4VEvaporationChannel.hh.

Referenced by G4EvaporationChannel::GetEmissionProbability().

◆ pairingCorrection

G4PairingCorrection* G4CompetitiveFission::pairingCorrection

private

Definition at line 112 of file G4CompetitiveFission.hh.

Referenced by EmittedFragment(), G4CompetitiveFission(), and GetEmissionProbability().

◆ theFissionBarrierPtr

G4VFissionBarrier* G4CompetitiveFission::theFissionBarrierPtr

private

Definition at line 105 of file G4CompetitiveFission.hh.

Referenced by G4CompetitiveFission(), GetEmissionProbability(), SetFissionBarrier(), and ~G4CompetitiveFission().

◆ theFissionProbabilityPtr

G4VEmissionProbability* G4CompetitiveFission::theFissionProbabilityPtr

private

Definition at line 108 of file G4CompetitiveFission.hh.

Referenced by G4CompetitiveFission(), GetEmissionProbability(), SetEmissionStrategy(), and ~G4CompetitiveFission().

◆ theLevelDensityPtr

G4VLevelDensityParameter* G4CompetitiveFission::theLevelDensityPtr

private

Definition at line 111 of file G4CompetitiveFission.hh.

Referenced by G4CompetitiveFission(), SetLevelDensityParameter(), and ~G4CompetitiveFission().

◆ theParam

G4FissionParameters G4CompetitiveFission::theParam

private

Definition at line 118 of file G4CompetitiveFission.hh.

Referenced by EmittedFragment(), FissionAtomicNumber(), FissionKineticEnergy(), and MassDistribution().

◆ theSecID

G4int G4CompetitiveFission::theSecID

private

Definition at line 120 of file G4CompetitiveFission.hh.

Referenced by EmittedFragment(), and G4CompetitiveFission().

◆ useSICB

G4bool G4VEvaporationChannel::useSICB

protectedinherited

Definition at line 94 of file G4VEvaporationChannel.hh.

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

source/processes/hadronic/models/de_excitation/fission/include/G4CompetitiveFission.hh
source/processes/hadronic/models/de_excitation/fission/src/G4CompetitiveFission.cc

Private Attributes
G4double	fissionBarrier

G4double	fissionProbability

G4double	maxKineticEnergy

G4bool	myOwnFissionBarrier

G4bool	myOwnFissionProbability

G4bool	myOwnLevelDensity

G4PairingCorrection *	pairingCorrection

G4VFissionBarrier *	theFissionBarrierPtr

G4VEmissionProbability *	theFissionProbabilityPtr

G4VLevelDensityParameter *	theLevelDensityPtr

G4FissionParameters	theParam

G4int	theSecID

Public Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4CompetitiveFission() [1/2]

◆ ~G4CompetitiveFission()

◆ G4CompetitiveFission() [2/2]

Member Function Documentation

◆ AsymmetricRatio()

◆ BreakUpChain()

◆ BreakUpFragment()

◆ Dump()

◆ EmittedFragment()

◆ FissionAtomicNumber()

◆ FissionCharge()

◆ FissionKineticEnergy()

◆ GetEmissionProbability()

◆ GetFissionBarrier()

◆ GetLevelDensityParameter()

◆ GetLifeTime()

◆ GetMaximalKineticEnergy()

◆ Initialise()

◆ LocalExp()

◆ MassDistribution()

◆ operator!=()

◆ operator=()

◆ operator==()

◆ Ratio()

◆ RDMForced()

◆ SetEmissionStrategy()

◆ SetFissionBarrier()

◆ SetICM()

◆ SetLevelDensityParameter()

◆ SetOPTxs()

◆ SymmetricRatio()

◆ UseSICB()

Field Documentation

◆ fissionBarrier

◆ fissionProbability

◆ maxKineticEnergy

◆ myOwnFissionBarrier

◆ myOwnFissionProbability

◆ myOwnLevelDensity

◆ OPTxs

◆ pairingCorrection

◆ theFissionBarrierPtr

◆ theFissionProbabilityPtr

◆ theLevelDensityPtr

◆ theParam

◆ theSecID

◆ useSICB