G4Fissioner Class Reference

#include <G4Fissioner.hh>

Inheritance diagram for G4Fissioner:

Public Member Functions
virtual void	collide (G4InuclParticle *bullet, G4InuclParticle *target, G4CollisionOutput &globalOutput)
virtual void	deExcite (const G4Fragment &target, G4CollisionOutput &output)
	G4Fissioner ()
virtual void	setVerboseLevel (G4int verbose=0)
virtual	~G4Fissioner ()

Protected Member Functions
virtual G4bool	explosion (const G4Fragment &target) const
virtual G4bool	explosion (G4int A, G4int Z, G4double excitation) const
void	getTargetData (const G4Fragment &target)
const G4Fragment &	makeFragment (G4int A, G4int Z, G4double EX=0.)
const G4Fragment &	makeFragment (G4LorentzVector mom, G4int A, G4int Z, G4double EX=0.)
virtual void	setName (const G4String &name)
virtual G4bool	validateOutput (const G4Fragment &target, const std::vector< G4InuclElementaryParticle > &particles)
virtual G4bool	validateOutput (const G4Fragment &target, const std::vector< G4InuclNuclei > &fragments)
virtual G4bool	validateOutput (const G4Fragment &target, G4CollisionOutput &output)

Protected Attributes
G4int	A
G4Fragment	aFragment
G4CascadeCheckBalance *	balance
G4double	EEXS
G4LorentzVector	PEX
G4String	theName
G4int	verboseLevel
G4int	Z

Private Member Functions
	G4Fissioner (const G4Fissioner &)
G4double	getC2 (G4int A1, G4int A2, G4double X3, G4double X4, G4double R12) const
G4double	getZopt (G4int A1, G4int A2, G4int ZT, G4double X3, G4double X4, G4double R12) const
G4Fissioner &	operator= (const G4Fissioner &)
void	potentialMinimization (G4double &VP, G4double(&ED)[2], G4double &VC, G4int AF, G4int AS, G4int ZF, G4int ZS, G4double AL1[2], G4double BET1[2], G4double &R12) const

Private Attributes
G4FissionStore	fissionStore

Detailed Description

Definition at line 49 of file G4Fissioner.hh.

Constructor & Destructor Documentation

G4Fissioner() [1/2]

G4Fissioner::G4Fissioner ( )

inline

Definition at line 51 of file G4Fissioner.hh.

: G4CascadeDeexciteBase("G4Fissioner") {;}

~G4Fissioner()

virtual G4Fissioner::~G4Fissioner ( )

inlinevirtual

Definition at line 52 of file G4Fissioner.hh.

{;}

G4Fissioner() [2/2]

G4Fissioner::G4Fissioner ( const G4Fissioner &  )

private

Member Function Documentation

collide()

void G4VCascadeDeexcitation::collide ( G4InuclParticle *  bullet, G4InuclParticle *  target, G4CollisionOutput &  globalOutput  )

virtualinherited

Implements G4VCascadeCollider.

Definition at line 37 of file G4VCascadeDeexcitation.cc.

                                           {
  if (verboseLevel) {
    G4cout << " >>> G4VCascadeDeexcitation[" << theName << "]::collide "
       << " *** SHOULD NOT BE CALLED ***" << G4endl;
  }
 
  throw G4HadronicException(__FILE__, __LINE__, 
    "G4VCascadeDeexcitation::collide() invalid, must use ::deExcite(G4Fagment*)");
}

References G4cout, G4endl, G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

Referenced by G4InuclEvaporation::BreakItUp().

deExcite()

void G4Fissioner::deExcite ( const G4Fragment &  target, G4CollisionOutput &  output  )

virtual

Implements G4VCascadeDeexcitation.

Definition at line 65 of file G4Fissioner.cc.

                                          {
  if (verboseLevel) {
    G4cout << " >>> G4Fissioner::deExcite" << G4endl;
  }
 
  if (verboseLevel > 1) 
    G4cout << " Fissioner input\n" << target << G4endl;
 
  // Initialize buffer for fission possibilities
  fissionStore.setVerboseLevel(verboseLevel);
  fissionStore.clear();
 
  getTargetData(target);
  G4double A13 = G4cbrt(A);
  G4double mass_in = PEX.m();
  G4double e_in = mass_in;      // Mass includes excitation
  G4double PARA = 0.055 * A13*A13 * (G4cbrt(A-Z) + G4cbrt(Z));
  G4double TEM = std::sqrt(EEXS / PARA);
  G4double TETA = 0.494 * A13 * TEM;
  
  TETA = TETA / std::sinh(TETA);
  
  if (A < 246) PARA += (nucleiLevelDensity(A) - PARA) * TETA;
  
  G4int A1 = A/2 + 1;
  G4int Z1;
  G4int A2 = A - A1;
 
  G4double ALMA = -1000.0;
  G4double DM1 = bindingEnergy(A,Z);
  G4double EVV = EEXS - DM1;
  G4double DM2 = bindingEnergyAsymptotic(A, Z);
  G4double DTEM = (A < 220 ? 0.5 : 1.15);
  
  TEM += DTEM;
  
  G4double AL1[2] = { -0.15, -0.15 };
  G4double BET1[2] = { 0.05, 0.05 };
 
  G4double R12 = G4cbrt(A1) + G4cbrt(A2); 
  
  for (G4int i = 0; i < 50 && A1 > 30; i++) {
    A1--;
    A2 = A - A1;
    G4double X3 = 1.0 / G4cbrt(A1);
    G4double X4 = 1.0 / G4cbrt(A2);
    Z1 = G4lrint(getZopt(A1, A2, Z, X3, X4, R12) - 1.);
    G4double EDEF1[2];
    G4int Z2 = Z - Z1;
    G4double VPOT, VCOUL;
    
    potentialMinimization(VPOT, EDEF1, VCOUL, A1, A2, Z1, Z2, AL1, BET1, R12);
    
    G4double DM3 = bindingEnergy(A1,Z1);
    G4double DM4 = bindingEnergyAsymptotic(A1, Z1);
    G4double DM5 = bindingEnergy(A2,Z2);
    G4double DM6 = bindingEnergyAsymptotic(A2, Z2);
    G4double DMT1 = DM4 + DM6 - DM2;
    G4double DMT = DM3 + DM5 - DM1;
    G4double EZL = EEXS + DMT - VPOT;
    
    if(EZL > 0.0) { // generate fluctuations
      //  faster, using randomGauss
      G4double C1 = std::sqrt(getC2(A1, A2, X3, X4, R12) / TEM);
      G4double DZ = randomGauss(C1);
      
      DZ = DZ > 0.0 ? DZ + 0.5 : -std::fabs(DZ - 0.5);
      Z1 += G4int(DZ);
      Z2 -= G4int(DZ);
      
      G4double DEfin = randomGauss(TEM);    
      G4double EZ = (DMT1 + (DMT - DMT1) * TETA - VPOT + DEfin) / TEM;
      
      if (EZ >= ALMA) ALMA = EZ;
      G4double EK = VCOUL + DEfin + 0.5 * TEM;
      G4double EV = EVV + bindingEnergy(A1,Z1) + bindingEnergy(A2,Z2) - EK;
      
      if (EV > 0.0) fissionStore.addConfig(A1, Z1, EZ, EK, EV);
    };
  };
  
  G4int store_size = fissionStore.size();
  if (store_size == 0) return;      // No fission products
 
  G4FissionConfiguration config = 
    fissionStore.generateConfiguration(ALMA, inuclRndm());
  
  A1 = G4int(config.afirst);
  A2 = A - A1;
  Z1 = G4int(config.zfirst);
  
  G4int Z2 = Z - Z1;
  
  G4double mass1 = G4InuclNuclei::getNucleiMass(A1,Z1);
  G4double mass2 = G4InuclNuclei::getNucleiMass(A2,Z2);
  G4double EK = config.ekin;
  G4double pmod = std::sqrt(0.001 * EK * mass1 * mass2 / mass_in);
  
  G4LorentzVector mom1 = generateWithRandomAngles(pmod, mass1);
  G4LorentzVector mom2; mom2.setVectM(-mom1.vect(), mass2);
  
  G4double e_out = mom1.e() + mom2.e();
  G4double EV = 1000.0 * (e_in - e_out) / A;
  if (EV <= 0.0) return;        // No fission energy
 
  G4double EEXS1 = EV*A1;
  G4double EEXS2 = EV*A2;
 
  // Pass only last two nuclear fragments
  output.addRecoilFragment(makeFragment(mom1, A1, Z1, EEXS1));
  output.addRecoilFragment(makeFragment(mom2, A2, Z2, EEXS2));
}

References G4CascadeDeexciteBase::A, A13, G4FissionStore::addConfig(), G4CollisionOutput::addRecoilFragment(), G4InuclSpecialFunctions::bindingEnergy(), G4InuclSpecialFunctions::bindingEnergyAsymptotic(), C1, G4FissionStore::clear(), CLHEP::HepLorentzVector::e(), G4CascadeDeexciteBase::EEXS, fissionStore, G4InuclSpecialFunctions::G4cbrt(), G4cout, G4endl, G4lrint(), G4FissionStore::generateConfiguration(), G4InuclSpecialFunctions::generateWithRandomAngles(), getC2(), G4InuclNuclei::getNucleiMass(), G4CascadeDeexciteBase::getTargetData(), getZopt(), G4InuclSpecialFunctions::inuclRndm(), CLHEP::HepLorentzVector::m(), G4CascadeDeexciteBase::makeFragment(), G4InuclSpecialFunctions::nucleiLevelDensity(), G4CascadeDeexciteBase::PEX, potentialMinimization(), G4InuclSpecialFunctions::randomGauss(), CLHEP::HepLorentzVector::setVectM(), G4FissionStore::setVerboseLevel(), G4FissionStore::size(), CLHEP::HepLorentzVector::vect(), G4VCascadeCollider::verboseLevel, G4CascadeDeexciteBase::Z, and anonymous_namespace{paraMaker.cc}::Z1.

Referenced by G4EquilibriumEvaporator::deExcite().

explosion() [1/2]

G4bool G4CascadeDeexciteBase::explosion ( const G4Fragment &  target ) const

protectedvirtualinherited

Reimplemented in G4EquilibriumEvaporator.

Definition at line 109 of file G4CascadeDeexciteBase.cc.

                                                                        {
  return explosion(fragment.GetA_asInt(), fragment.GetZ_asInt(),
           fragment.GetExcitationEnergy());     // in MeV
}

References G4CascadeDeexciteBase::explosion(), G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), and G4Fragment::GetZ_asInt().

Referenced by G4CascadeDeexcitation::deExcite(), G4PreCompoundDeexcitation::deExcite(), G4CascadeDeexciteBase::explosion(), and G4EquilibriumEvaporator::explosion().

explosion() [2/2]

G4bool G4CascadeDeexciteBase::explosion ( G4int  A, G4int  Z, G4double  excitation  ) const

protectedvirtualinherited

Reimplemented in G4EquilibriumEvaporator.

Definition at line 114 of file G4CascadeDeexciteBase.cc.

                                                {
  if (verboseLevel) G4cout << " >>> " << theName << "::explosion ?" << G4endl;
 
  const G4int a_cut = 20;
  const G4double be_cut = 3.0;
 
  // Neutron balls, or small fragments with high excitations can explode
  return ((fragA <= a_cut || fragZ==0) && 
      (excitation >= be_cut * bindingEnergy(fragA,fragZ))
      );
}

References G4InuclSpecialFunctions::bindingEnergy(), G4cout, G4endl, G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

getC2()

G4double G4Fissioner::getC2 ( G4int  A1, G4int  A2, G4double  X3, G4double  X4, G4double  R12  ) const

private

Definition at line 179 of file G4Fissioner.cc.

                                    {
 
  if (verboseLevel > 3) {
    G4cout << " >>> G4Fissioner::getC2" << G4endl;
  }
 
  G4double C2 = 124.57 * (1.0 / A1 + 1.0 / A2) + 0.78 * (X3 + X4) - 176.9 *
    ((X3*X3*X3*X3) + (X4*X4*X4*X4)) + 219.36 * (1.0 / (A1 * A1) + 1.0 / (A2 * A2)) - 1.108 / R12;
 
  return C2;   
}

References C2, G4cout, G4endl, and G4VCascadeCollider::verboseLevel.

Referenced by deExcite(), and getZopt().

getTargetData()

void G4CascadeDeexciteBase::getTargetData ( const G4Fragment &  target )

protectedinherited

Definition at line 66 of file G4CascadeDeexciteBase.cc.

                                                                  {
  A = target.GetA_asInt();
  Z = target.GetZ_asInt();
  PEX = target.GetMomentum()/GeV;   // Convert from G4 to Bertini units
  EEXS = target.GetExcitationEnergy();
}

References G4CascadeDeexciteBase::A, G4CascadeDeexciteBase::EEXS, G4Fragment::GetA_asInt(), G4Fragment::GetExcitationEnergy(), G4Fragment::GetMomentum(), G4Fragment::GetZ_asInt(), GeV, G4CascadeDeexciteBase::PEX, and G4CascadeDeexciteBase::Z.

Referenced by G4BigBanger::deExcite(), G4EquilibriumEvaporator::deExcite(), deExcite(), and G4NonEquilibriumEvaporator::deExcite().

getZopt()

G4double G4Fissioner::getZopt ( G4int  A1, G4int  A2, G4int  ZT, G4double  X3, G4double  X4, G4double  R12  ) const

private

Definition at line 195 of file G4Fissioner.cc.

                                      {
 
  if (verboseLevel > 3) {
    G4cout << " >>> G4Fissioner::getZopt" << G4endl;
  }
 
  G4double Zopt = (87.7 * (X4 - X3) * (1.0 - 1.25 * (X4 + X3)) +
           ZT * ((124.57 / A2 + 0.78 * X4 - 176.9 * (X4*X4*X4*X4) + 219.36 / (A2 * A2)) - 0.554 / R12)) /
    getC2(A1, A2, X3, X4, R12);
 
  return Zopt;
}        

References G4cout, G4endl, getC2(), and G4VCascadeCollider::verboseLevel.

Referenced by deExcite().

makeFragment() [1/2]

const G4Fragment & G4CascadeDeexciteBase::makeFragment ( G4int  A, G4int  Z, G4double  EX = 0.  )

protectedinherited

Definition at line 81 of file G4CascadeDeexciteBase.cc.

                                                                         {
  return makeFragment(zero, fragA, fragZ, EX);
}

References G4CascadeDeexciteBase::makeFragment(), and anonymous_namespace{G4CascadeDeexciteBase.cc}::zero.

makeFragment() [2/2]

const G4Fragment & G4CascadeDeexciteBase::makeFragment ( G4LorentzVector  mom, G4int  A, G4int  Z, G4double  EX = 0.  )

protectedinherited

Definition at line 86 of file G4CascadeDeexciteBase.cc.

                                              {
  if (verboseLevel>2) {
    G4cout << " >>> " << theName << "::makeFragment " << mom << " " << fragA
       << " " << fragZ << " " << EX << G4endl;
  }
 
  // Adjust four-momentum so that mass is nucleus + excitation
  G4double mass =
    G4InuclNuclei::getNucleiMass(fragA,fragZ) + EX/GeV;
  mom.setVectM(mom.vect(), mass);
 
  // Overwrite previous fragment contents, zeroing out excitons
  aFragment.SetZandA_asInt(fragZ, fragA);
  aFragment.SetMomentum(mom*GeV);           // Bertini uses GeV!
  aFragment.SetNumberOfHoles(0,0);
  aFragment.SetNumberOfExcitedParticle(0,0);
 
  return aFragment;
}

References G4CascadeDeexciteBase::aFragment, G4cout, G4endl, G4InuclNuclei::getNucleiMass(), GeV, G4Fragment::SetMomentum(), G4Fragment::SetNumberOfExcitedParticle(), G4Fragment::SetNumberOfHoles(), CLHEP::HepLorentzVector::setVectM(), G4Fragment::SetZandA_asInt(), G4VCascadeCollider::theName, CLHEP::HepLorentzVector::vect(), and G4VCascadeCollider::verboseLevel.

Referenced by G4EquilibriumEvaporator::deExcite(), deExcite(), G4NonEquilibriumEvaporator::deExcite(), and G4CascadeDeexciteBase::makeFragment().

operator=()

G4Fissioner & G4Fissioner::operator= ( const G4Fissioner &  )

private

potentialMinimization()

void G4Fissioner::potentialMinimization ( G4double &  VP, G4double(&)  ED[2], G4double &  VC, G4int  AF, G4int  AS, G4int  ZF, G4int  ZS, G4double  AL1[2], G4double  BET1[2], G4double &  R12  ) const

private

Definition at line 213 of file G4Fissioner.cc.

                                         {
 
  if (verboseLevel > 3) {
    G4cout << " >>> G4Fissioner::potentialMinimization" << G4endl;
  }
 
  const G4double huge_num = 2.0e35;
  const G4int itry_max = 2000;
  const G4double DSOL1 = 1.0e-6;
  const G4double DS1 = 0.3;
  const G4double DS2 = 1.0 / DS1 / DS1; 
  G4int A1[2] = { AF, AS };
  G4int Z1[2] = { ZF, ZS };
  G4double D = 1.01844 * ZF * ZS;
  G4double D0 = 1.0e-3 * D;
  G4double R[2];
  R12 = 0.0;
  G4double C[2];
  G4double F[2];
  G4double Y1;
  G4double Y2;
  G4int i;
 
  for (i = 0; i < 2; i++) {
    R[i] = G4cbrt(A1[i]);
    Y1 = R[i] * R[i];
    Y2 = Z1[i] * Z1[i] / R[i];
    C[i] = 6.8 * Y1 - 0.142 * Y2;
    F[i] = 12.138 * Y1 - 0.145 * Y2; 
  };
 
  G4double SAL[2];
  G4double SBE[2];
  G4double X[2];
  G4double X1[2];
  G4double X2[2];
  G4double RAL[2];
  G4double RBE[2];
  G4double AA[4][4];
  G4double B[4];
  G4int itry = 0;
 
  while (itry < itry_max) { /* Loop checking 08.06.2015 MHK */
    itry++;
    G4double S = 0.0;
 
    for (i = 0; i < 2; i++) {
      S += R[i] * (1.0 + AL1[i] + BET1[i] - 0.257 * AL1[i] * BET1[i]);
    };
    R12 = 0.0;
    Y1 = 0.0;
    Y2 = 0.0;
 
    for (i = 0; i < 2; i++) {
      SAL[i] = R[i] * (1.0-0.257 * BET1[i]);
      SBE[i] = R[i] * (1.0-0.257 * AL1[i]);
      X[i] = R[i] / S;
      X1[i] = X[i] * X[i];
      X2[i] = X[i] * X1[i];
      Y1 += AL1[i] * X1[i];
      Y2 += BET1[i] * X2[i];
      R12 += R[i] * (1.0 - AL1[i] * (1.0 - 0.6 * X[i]) + BET1[i] * (1.0 - 0.429 * X1[i]));
    }; 
 
    G4double Y3 = -0.6 * Y1 + 0.857 * Y2;
    G4double Y4 = (1.2 * Y1 - 2.571 * Y2) / S;
    G4double R2 = D0 / (R12 * R12);
    G4double R3 = 2.0 * R2 / R12;
 
    for (i = 0; i < 2; i++) {
      RAL[i] = -R[i] * (1.0 - 0.6 * X[i]) + SAL[i] * Y3;
      RBE[i] =  R[i] * (1.0 - 0.429 * X1[i]) + SBE[i] * Y3;
    };
 
    G4double DX1;
    G4double DX2;
  
    for (i = 0; i < 2; i++) {
 
      for (G4int j = 0; j < 2; j++) {
    G4double DEL1 = i == j ? 1.0 : 0.0;
    DX1 = 0.0;
    DX2 = 0.0;
 
    if (std::fabs(AL1[i]) >= DS1) {
      G4double XXX = AL1[i] * AL1[i] * DS2;
      G4double DEX = XXX > 100.0 ? huge_num : G4Exp(XXX);
      DX1 = 2.0 * (1.0 + 2.0 * AL1[i] * AL1[i] * DS2) * DEX * DS2;
    };
 
    if (std::fabs(BET1[i]) >= DS1) {
      G4double XXX = BET1[i] * BET1[i] * DS2;
      G4double DEX = XXX > 100.0 ? huge_num : G4Exp(XXX);
      DX2 = 2.0 * (1.+2.0 * BET1[i] * BET1[i] * DS2) * DEX * DS2;
    };
 
    G4double DEL = 2.0e-3 * DEL1;
    AA[i][j] = R3 * RBE[i] * RBE[j] - 
      R2 * (-0.6 * 
        (X1[i] * SAL[j] + 
         X1[j] * SAL[i]) + SAL[i] * SAL[j] * Y4) + 
      DEL * C[i] + DEL1 * DX1;
    G4int i1 = i + 2;
    G4int j1 = j + 2;
    AA[i1][j1] = R3 * RBE[i] * RBE[j] 
      - R2 * (0.857 * 
          (X2[i] * SBE[j] + 
           X2[j] * SBE[i]) + SBE[i] * SBE[j] * Y4) +
      DEL * F[i] + DEL1 * DX2;
    AA[i][j1] = R3 * RAL[i] * RBE[j] - 
      R2 * (0.857 * 
        (X2[j] * SAL[i] - 
         0.6 * X1[i] * SBE[j]) + SBE[j] * SAL[i] * Y4 - 
        0.257 * R[i] * Y3 * DEL1);
    AA[j1][i] = AA[i][j1];       
      };
    };
 
    for (i = 0; i < 2; i++) {
      DX1 = 0.0;
      DX2 = 0.0;
 
      if (std::fabs(AL1[i]) >= DS1) DX1 = 2.0 * AL1[i] * DS2 * G4Exp(AL1[i] * AL1[i] * DS2);
 
      if (std::fabs(BET1[i]) >= DS1) DX2 = 2.0 * BET1[i] * DS2 * G4Exp(BET1[i] * BET1[i] * DS2);
      B[i] =     R2 * RAL[i] - 2.0e-3 * C[i] * AL1[i] + DX1;
      B[i + 2] = R2 * RBE[i] - 2.0e-3 * F[i] * BET1[i] + DX2;
    };
 
    G4double ST = 0.0;
    G4double ST1 = 0.0;
 
    for (i = 0; i < 4; i++) {
      ST += B[i] * B[i];
 
      for (G4int j = 0; j < 4; j++) ST1 += AA[i][j] * B[i] * B[j];
    };
 
    G4double STEP = ST / ST1;
    G4double DSOL = 0.0;
 
    for (i = 0; i < 2; i++) {
      AL1[i] += B[i] * STEP;
      BET1[i] += B[i + 2] * STEP;
      DSOL += B[i] * B[i] + B[i + 2] * B[i + 2]; 
    };
    DSOL = std::sqrt(DSOL);
 
    if (DSOL < DSOL1) break;
  };
 
  if (verboseLevel > 3) {
  if (itry == itry_max) 
    G4cout << " maximal number of iterations in potentialMinimization " << G4endl
       << " A1 " << AF << " Z1 " << ZF << G4endl; 
 
  };
 
  for (i = 0; i < 2; i++) ED[i] = F[i] * BET1[i] * BET1[i] + C[i] * AL1[i] * AL1[i]; 
 
  VC = D / R12;
  VP = VC + ED[0] + ED[1];
}

References anonymous_namespace{paraMaker.cc}::AA, B(), C(), D(), G4InuclSpecialFunctions::G4cbrt(), G4cout, G4endl, G4Exp(), S(), G4VCascadeCollider::verboseLevel, and anonymous_namespace{paraMaker.cc}::Z1.

Referenced by deExcite().

setName()

virtual void G4VCascadeCollider::setName ( const G4String &  name )

inlineprotectedvirtualinherited

Definition at line 55 of file G4VCascadeCollider.hh.

{ theName = name; }

References G4InuclParticleNames::name(), and G4VCascadeCollider::theName.

Referenced by G4CascadeCheckBalance::setOwner().

setVerboseLevel()

void G4CascadeDeexciteBase::setVerboseLevel ( G4int  verbose = 0 )

virtualinherited

Reimplemented from G4VCascadeCollider.

Reimplemented in G4CascadeDeexcitation, G4EquilibriumEvaporator, and G4PreCompoundDeexcitation.

Definition at line 58 of file G4CascadeDeexciteBase.cc.

                                                         {
  G4VCascadeDeexcitation::setVerboseLevel(verbose);
  if (balance) balance->setVerboseLevel(verbose);
}

References G4CascadeDeexciteBase::balance, and G4VCascadeCollider::setVerboseLevel().

Referenced by G4CascadeDeexcitation::setVerboseLevel(), G4EquilibriumEvaporator::setVerboseLevel(), and G4PreCompoundDeexcitation::setVerboseLevel().

validateOutput() [1/3]

G4bool G4CascadeDeexciteBase::validateOutput ( const G4Fragment &  target, const std::vector< G4InuclElementaryParticle > &  particles  )

protectedvirtualinherited

Definition at line 142 of file G4CascadeDeexciteBase.cc.

                                                                    {
  if (!balance) return true;        // Skip checks unless requested
 
  if (verboseLevel > 1)
    G4cout << " >>> " << theName << "::validateOutput" << G4endl;
 
  balance->setVerboseLevel(verboseLevel);
  balance->collide(target, particles);
  return balance->okay();           // Returns false if violations
}

References G4CascadeDeexciteBase::balance, G4CascadeCheckBalance::collide(), G4cout, G4endl, G4CascadeCheckBalance::okay(), G4VCascadeCollider::setVerboseLevel(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

validateOutput() [2/3]

G4bool G4CascadeDeexciteBase::validateOutput ( const G4Fragment &  target, const std::vector< G4InuclNuclei > &  fragments  )

protectedvirtualinherited

Definition at line 154 of file G4CascadeDeexciteBase.cc.

                                                        {
  if (!balance) return true;        // Skip checks unless requested
 
  if (verboseLevel > 1)
    G4cout << " >>> " << theName << "::validateOutput" << G4endl;
 
  balance->setVerboseLevel(verboseLevel);
  balance->collide(target, fragments);
  return balance->okay();           // Returns false if violations
}

References G4CascadeDeexciteBase::balance, G4CascadeCheckBalance::collide(), G4cout, G4endl, G4CascadeCheckBalance::okay(), G4VCascadeCollider::setVerboseLevel(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

validateOutput() [3/3]

G4bool G4CascadeDeexciteBase::validateOutput ( const G4Fragment &  target, G4CollisionOutput &  output  )

protectedvirtualinherited

Definition at line 130 of file G4CascadeDeexciteBase.cc.

                                                    {
  if (!balance) return true;        // Skip checks unless requested
 
  if (verboseLevel > 1)
    G4cout << " >>> " << theName << "::validateOutput" << G4endl;
 
  balance->setVerboseLevel(verboseLevel);
  balance->collide(target, output);
  return balance->okay();           // Returns false if violations
}

References G4CascadeDeexciteBase::balance, G4CascadeCheckBalance::collide(), G4cout, G4endl, G4CascadeCheckBalance::okay(), G4VCascadeCollider::setVerboseLevel(), G4VCascadeCollider::theName, and G4VCascadeCollider::verboseLevel.

Referenced by G4BigBanger::deExcite(), G4EquilibriumEvaporator::deExcite(), and G4NonEquilibriumEvaporator::deExcite().

Field Documentation

A

G4int G4CascadeDeexciteBase::A

protectedinherited

Definition at line 73 of file G4CascadeDeexciteBase.hh.

Referenced by G4BigBanger::deExcite(), G4EquilibriumEvaporator::deExcite(), deExcite(), G4NonEquilibriumEvaporator::deExcite(), G4CascadeDeexciteBase::getTargetData(), and mcscore.MCParticle::printout().

aFragment

G4Fragment G4CascadeDeexciteBase::aFragment

protectedinherited

Definition at line 82 of file G4CascadeDeexciteBase.hh.

Referenced by G4CascadeDeexciteBase::makeFragment().

balance

G4CascadeCheckBalance* G4CascadeDeexciteBase::balance

protectedinherited

Definition at line 57 of file G4CascadeDeexciteBase.hh.

Referenced by G4CascadeDeexciteBase::G4CascadeDeexciteBase(), G4CascadeDeexciteBase::setVerboseLevel(), G4CascadeDeexciteBase::validateOutput(), and G4CascadeDeexciteBase::~G4CascadeDeexciteBase().

EEXS

G4double G4CascadeDeexciteBase::EEXS

protectedinherited

Definition at line 76 of file G4CascadeDeexciteBase.hh.

Referenced by G4BigBanger::deExcite(), G4EquilibriumEvaporator::deExcite(), deExcite(), G4NonEquilibriumEvaporator::deExcite(), and G4CascadeDeexciteBase::getTargetData().

fissionStore

G4FissionStore G4Fissioner::fissionStore

private

Definition at line 57 of file G4Fissioner.hh.

Referenced by deExcite().

PEX

G4LorentzVector G4CascadeDeexciteBase::PEX

protectedinherited

Definition at line 75 of file G4CascadeDeexciteBase.hh.

Referenced by G4BigBanger::deExcite(), G4EquilibriumEvaporator::deExcite(), deExcite(), G4NonEquilibriumEvaporator::deExcite(), and G4CascadeDeexciteBase::getTargetData().

theName

G4String G4VCascadeCollider::theName

protectedinherited

Definition at line 52 of file G4VCascadeCollider.hh.

Referenced by G4CascadeCheckBalance::baryonOkay(), G4CascadeCheckBalance::chargeOkay(), G4VCascadeDeexcitation::collide(), G4CascadeCheckBalance::collide(), G4CascadeCheckBalance::ekinOkay(), G4CascadeCheckBalance::energyOkay(), G4CascadeDeexciteBase::explosion(), G4CascadeRecoilMaker::G4CascadeRecoilMaker(), G4VCascadeCollider::G4VCascadeCollider(), G4CascadeColliderBase::inelasticInteractionPossible(), G4CascadeDeexciteBase::makeFragment(), G4CascadeRecoilMaker::makeRecoilFragment(), G4CascadeRecoilMaker::makeRecoilNuclei(), G4CascadeCheckBalance::momentumOkay(), G4VCascadeCollider::setName(), G4CascadeCheckBalance::strangeOkay(), G4CascadeDeexciteBase::validateOutput(), and G4CascadeColliderBase::validateOutput().

verboseLevel

G4int G4VCascadeCollider::verboseLevel

protectedinherited

Definition at line 53 of file G4VCascadeCollider.hh.

Referenced by G4CascadeCheckBalance::baryonOkay(), G4CascadeCheckBalance::chargeOkay(), G4IntraNucleiCascader::collide(), G4InuclCollider::collide(), G4LightTargetCollider::collide(), G4VCascadeDeexcitation::collide(), G4CascadeCheckBalance::collide(), G4CascadeRecoilMaker::collide(), G4ElementaryParticleCollider::collide(), G4IntraNucleiCascader::copySecondaries(), G4IntraNucleiCascader::copyWoundedNucleus(), G4IntraNucleiCascader::decayTrappedParticle(), G4CascadeDeexcitation::deExcite(), G4InuclCollider::deexcite(), G4PreCompoundDeexcitation::deExcite(), G4BigBanger::deExcite(), G4EquilibriumEvaporator::deExcite(), G4EvaporationInuclCollider::deExcite(), deExcite(), G4NonEquilibriumEvaporator::deExcite(), G4CascadeCheckBalance::ekinOkay(), G4CascadeCheckBalance::energyOkay(), G4EquilibriumEvaporator::explosion(), G4CascadeDeexciteBase::explosion(), G4CascadeRecoilMaker::fillRecoil(), G4IntraNucleiCascader::finalize(), G4IntraNucleiCascader::finishCascade(), G4VCascadeCollider::G4VCascadeCollider(), G4BigBanger::generateBangInSCM(), G4IntraNucleiCascader::generateCascade(), G4BigBanger::generateMomentumModules(), G4ElementaryParticleCollider::generateMultiplicity(), G4ElementaryParticleCollider::generateSCMfinalState(), G4ElementaryParticleCollider::generateSCMmuonAbsorption(), G4ElementaryParticleCollider::generateSCMpionAbsorption(), G4ElementaryParticleCollider::generateSCMpionNAbsorption(), G4BigBanger::generateX(), G4EquilibriumEvaporator::getAF(), getC2(), G4EquilibriumEvaporator::getE0(), G4NonEquilibriumEvaporator::getE0(), G4NonEquilibriumEvaporator::getMatrixElement(), G4NonEquilibriumEvaporator::getParLev(), G4EquilibriumEvaporator::getPARLEVDEN(), G4EquilibriumEvaporator::getQF(), getZopt(), G4CascadeRecoilMaker::goodNucleus(), G4EquilibriumEvaporator::goodRemnant(), G4CascadeColliderBase::inelasticInteractionPossible(), G4IntraNucleiCascader::initialize(), G4CascadeDeexciteBase::makeFragment(), G4CascadeRecoilMaker::makeRecoilFragment(), G4CascadeRecoilMaker::makeRecoilNuclei(), G4BigBanger::maxProbability(), G4CascadeCheckBalance::momentumOkay(), G4IntraNucleiCascader::newCascade(), G4InuclCollider::photonuclearOkay(), G4ElementaryParticleCollider::pionNucleonAbsorption(), potentialMinimization(), G4IntraNucleiCascader::preloadCascade(), G4IntraNucleiCascader::processSecondary(), G4IntraNucleiCascader::processTrappedParticle(), G4IntraNucleiCascader::releaseSecondary(), G4IntraNucleiCascader::rescatter(), G4InuclCollider::rescatter(), G4IntraNucleiCascader::setupCascade(), G4InuclCollider::setVerboseLevel(), G4LightTargetCollider::setVerboseLevel(), G4VCascadeCollider::setVerboseLevel(), G4CascadeCheckBalance::strangeOkay(), G4InuclCollider::useCascadeDeexcitation(), G4InuclCollider::usePreCompoundDeexcitation(), G4CascadeDeexciteBase::validateOutput(), G4CascadeColliderBase::validateOutput(), G4CascadeRecoilMaker::wholeEvent(), and G4BigBanger::xProbability().

Z

G4int G4CascadeDeexciteBase::Z

protectedinherited

Definition at line 74 of file G4CascadeDeexciteBase.hh.

Referenced by G4BigBanger::deExcite(), G4EquilibriumEvaporator::deExcite(), deExcite(), G4NonEquilibriumEvaporator::deExcite(), G4CascadeDeexciteBase::getTargetData(), and mcscore.MCParticle::printout().

---

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

• source/processes/hadronic/models/cascade/cascade/include/G4Fissioner.hh
• source/processes/hadronic/models/cascade/cascade/src/G4Fissioner.cc