#include <G4PreCompoundNeutron.hh>

Inheritance diagram for G4PreCompoundNeutron:

Public Member Functions
	G4PreCompoundNeutron ()

virtual	~G4PreCompoundNeutron ()

Public Member Functions inherited from G4PreCompoundNucleon
	G4PreCompoundNucleon (const G4ParticleDefinition , G4VCoulombBarrier aCoulombBarrier)

virtual	~G4PreCompoundNucleon ()

Public Member Functions inherited from G4PreCompoundFragment
	G4PreCompoundFragment (const G4ParticleDefinition , G4VCoulombBarrier aCoulombBarrier)

virtual	~G4PreCompoundFragment ()

G4double	CalcEmissionProbability (const G4Fragment &aFragment)

G4double	GetKineticEnergy (const G4Fragment &aFragment)

Public Member Functions inherited from G4VPreCompoundFragment
	G4VPreCompoundFragment (const G4ParticleDefinition , G4VCoulombBarrier aCoulombBarrier)

virtual	~G4VPreCompoundFragment ()

void	Initialize (const G4Fragment &aFragment)

virtual G4double	CalcEmissionProbability (const G4Fragment &aFragment)=0

virtual G4double	GetKineticEnergy (const G4Fragment &aFragment)=0

G4ReactionProduct *	GetReactionProduct () const

G4int	GetA () const

G4int	GetZ () const

G4int	GetRestA () const

G4int	GetRestZ () const

G4double	ResidualA13 () const

G4double	GetCoulombBarrier () const

G4double	GetBindingEnergy () const

G4double	GetMaximalKineticEnergy () const

G4double	GetEnergyThreshold () const

G4double	GetEmissionProbability () const

G4double	GetNuclearMass () const

G4double	GetRestNuclearMass () const

G4double	GetReducedMass () const

const G4LorentzVector &	GetMomentum () const

void	SetMomentum (const G4LorentzVector &value)

const G4String	GetName () const

void	SetOPTxs (G4int)

void	UseSICB (G4bool)

Protected Member Functions
virtual G4double	GetRj (G4int NumberParticles, G4int NumberCharged)

virtual G4double	CrossSection (G4double ekin)

virtual G4double	GetAlpha ()

virtual G4double	GetBeta ()

G4double	GetOpt12 (G4double K)

G4double	GetOpt34 (G4double K)

Protected Member Functions inherited from G4PreCompoundNucleon
virtual G4double	ProbabilityDistributionFunction (G4double eKin, const G4Fragment &aFragment)

virtual G4double	CrossSection (G4double ekin)=0

virtual G4double	GetRj (G4int NumberParticles, G4int NumberCharged)=0

virtual G4double	GetAlpha ()=0

virtual G4double	GetBeta ()=0

G4double	GetOpt0 (G4double ekin)

virtual G4double	ProbabilityDistributionFunction (G4double K, const G4Fragment &aFragment)=0

Protected Member Functions inherited from G4VPreCompoundFragment
G4bool	IsItPossible (const G4Fragment &aFragment) const

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundFragment
G4PreCompoundParameters *	theParameters

G4Pow *	g4pow

G4double	theEmissionProbability

G4double	theCoulombBarrier

G4int	OPTxs

G4bool	useSICB

Detailed Description

Definition at line 42 of file G4PreCompoundNeutron.hh.

Constructor & Destructor Documentation

◆ G4PreCompoundNeutron()

G4PreCompoundNeutron::G4PreCompoundNeutron ( )

Definition at line 48 of file G4PreCompoundNeutron.cc.

  : G4PreCompoundNucleon(G4Neutron::Neutron(), &theNeutronCoulombBarrier)
{
  ResidualA = GetRestA();
  ResidualZ = GetRestZ(); 
  theA = GetA();
  theZ = GetZ();
  ResidualAthrd = ResidualA13();
  FragmentAthrd = ResidualAthrd;
  FragmentA = theA + ResidualA;
}

◆ ~G4PreCompoundNeutron()

G4PreCompoundNeutron::~G4PreCompoundNeutron ( )

virtual

Definition at line 60 of file G4PreCompoundNeutron.cc.

61{}

Member Function Documentation

◆ CrossSection()

G4double G4PreCompoundNeutron::CrossSection ( G4double ekin )

protectedvirtual

Implements G4PreCompoundNucleon.

Definition at line 79 of file G4PreCompoundNeutron.cc.

{
  ResidualA = GetRestA();
  ResidualZ = GetRestZ(); 
  theA = GetA();
  theZ = GetZ();
  ResidualAthrd = ResidualA13();
  FragmentA = theA + ResidualA;
  FragmentAthrd = g4pow->Z13(FragmentA);
 
  if (OPTxs==0) { return GetOpt0( K); }
  else if( OPTxs==1 || OPTxs==2) { return GetOpt12( K); }
  else if (OPTxs==3 || OPTxs==4) { return GetOpt34( K); }
  else{
    std::ostringstream errOs;
    errOs << "BAD NEUTRON CROSS SECTION OPTION !!"  <<G4endl;
    throw G4HadronicException(__FILE__, __LINE__, errOs.str());
    return 0.;
  }
}

◆ GetAlpha()

G4double G4PreCompoundNeutron::GetAlpha ( )

protectedvirtual

Implements G4PreCompoundNucleon.

Definition at line 100 of file G4PreCompoundNeutron.cc.

{
  return 0.76+2.2/ResidualAthrd;
}

Referenced by GetBeta().

◆ GetBeta()

G4double G4PreCompoundNeutron::GetBeta ( )

protectedvirtual

Implements G4PreCompoundNucleon.

Definition at line 105 of file G4PreCompoundNeutron.cc.

{
  //   return (2.12/std::pow(GetRestA(),2.0/3.0)-0.05)*MeV/GetAlpha();
  return (2.12/(ResidualAthrd*ResidualAthrd)-0.05)*MeV/GetAlpha();
}

◆ GetOpt12()

G4double G4PreCompoundNeutron::GetOpt12 ( G4double K )

protected

Definition at line 114 of file G4PreCompoundNeutron.cc.

{
  G4double Kc=K;
 
  // Pramana (Bechetti & Greenles) for neutrons is chosen 
 
  // JMQ  xsec is set constat above limit of validity
  if (K > 50*MeV) { Kc = 50*MeV; }
 
  G4double landa, landa0, landa1, mu, mm0, mu1,nu, nu0, nu1, nu2,xs;
 
  landa0 = 18.57;
  landa1 = -22.93;
  mm0 = 381.7;
  mu1 = 24.31;
  nu0 = 0.172;
  nu1 = -15.39;
  nu2 = 804.8;
  landa = landa0/ResidualAthrd + landa1;
  mu = mm0*ResidualAthrd + mu1*ResidualAthrd*ResidualAthrd;
  nu = nu0*ResidualAthrd*ResidualA + nu1*ResidualAthrd*ResidualAthrd + nu2 ;
  xs=landa*Kc + mu + nu/Kc;
  if (xs <= 0.0 ){
    std::ostringstream errOs;
    G4cout<<"WARNING:  NEGATIVE OPT=1 neutron cross section "<<G4endl;     
    errOs << "RESIDUAL: Ar=" << ResidualA << " Zr=" << ResidualZ <<G4endl;
    errOs <<"  xsec("<<Kc<<" MeV) ="<<xs <<G4endl;
    throw G4HadronicException(__FILE__, __LINE__, errOs.str());
              }
  return xs;
}

Referenced by CrossSection().

◆ GetOpt34()

G4double G4PreCompoundNeutron::GetOpt34 ( G4double K )

protected

Definition at line 147 of file G4PreCompoundNeutron.cc.

{
  G4double landa, landa0, landa1, mu, mm0, mu1,nu, nu0, nu1, nu2;
  G4double p, p0;
  G4double flow,ec,ecsq,xnulam,etest(0.),ra(0.),a,signor(1.),sig; 
  G4double b,ecut,cut,ecut2,geom,elab;
 
  flow = 1.e-18;
 
  // PRECO xs for neutrons is choosen
  p0 = -312.;
  landa0 = 12.10;
  landa1=  -11.27;
  mm0 = 234.1;
  mu1 = 38.26;
  nu0 = 1.55;
  nu1 = -106.1;
  nu2 = 1280.8; 
 
  if (ResidualA < 40)  { signor =0.7 + ResidualA*0.0075; }
  if (ResidualA > 210) { signor = 1. + (ResidualA-210)/250.; }
  landa = landa0/ResidualAthrd + landa1;
  mu = mm0*ResidualAthrd + mu1*ResidualAthrd*ResidualAthrd;
  nu = nu0*ResidualAthrd*ResidualA + nu1*ResidualAthrd*ResidualAthrd + nu2;
 
  // JMQ very low energy behaviour corrected (problem  for A (apprx.)>60)
  if (nu < 0.) { nu=-nu; }
 
  ec = 0.5;
  ecsq = 0.25;
  p = p0;
  xnulam = 1.;
  etest = 32.;
  //          ** etest is the energy above which the rxn cross section is
  //          ** compared with the geometrical limit and the max taken.
  //          ** xnulam here is a dummy value to be used later.
 
  a = -2.*p*ec + landa - nu/ecsq;
  b = p*ecsq + mu + 2.*nu/ec;
  ecut = 0.;
  cut = a*a - 4.*p*b;
  if (cut > 0.) { ecut = std::sqrt(cut); }
  ecut = (ecut-a) / (p+p);
  ecut2 = ecut;
  if (cut < 0.) { ecut2 = ecut - 2.; }
  elab = K * FragmentA / G4double(ResidualA);
  sig = 0.;
  if (elab <= ec) { //start for E<Ec 
    if (elab > ecut2) { sig = (p*elab*elab+a*elab+b) * signor; } 
  }              //end for E<Ec
  else {           //start for  E>Ec
    sig = (landa*elab+mu+nu/elab) * signor;
    geom = 0.;
    if (xnulam < flow || elab < etest) { return sig; }
    geom = std::sqrt(theA*K);
    geom = 1.23*ResidualAthrd + ra + 4.573/geom;
    geom = 31.416 * geom * geom;
    sig = std::max(geom,sig);
 
  }
  return sig;
}

Referenced by CrossSection().

◆ GetRj()

G4double G4PreCompoundNeutron::GetRj	(	G4int	NumberParticles,
		G4int	NumberCharged
	)

protectedvirtual

Implements G4PreCompoundNucleon.

Definition at line 63 of file G4PreCompoundNeutron.cc.

{
  G4double rj = 0.0;
  if(nParticles > 0) { 
    rj = static_cast<G4double>(nParticles - nCharged)/
      static_cast<G4double>(nParticles);
  }
  return rj;
}

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

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4PreCompoundNeutron()

◆ ~G4PreCompoundNeutron()

Member Function Documentation

◆ CrossSection()

◆ GetAlpha()

◆ GetBeta()

◆ GetOpt12()

◆ GetOpt34()

◆ GetRj()