G4HETCNeutron Class Reference

#include <G4HETCNeutron.hh>

Inheritance diagram for G4HETCNeutron:

Public Member Functions
	G4HETCNeutron ()
	~G4HETCNeutron () override=default
G4double	SampleKineticEnergy (const G4Fragment &) override
	G4HETCNeutron (const G4HETCNeutron &right)=delete
const G4HETCNeutron &	operator= (const G4HETCNeutron &right)=delete
G4bool	operator== (const G4HETCNeutron &right) const =delete
G4bool	operator!= (const G4HETCNeutron &right) const =delete
Public Member Functions inherited from G4HETCFragment
	G4HETCFragment (const G4ParticleDefinition *, G4VCoulombBarrier *p=nullptr)
	~G4HETCFragment () override=default
G4double	CalcEmissionProbability (const G4Fragment &aFragment) override
	G4HETCFragment (const G4HETCFragment &right)=delete
const G4HETCFragment &	operator= (const G4HETCFragment &right)=delete
G4bool	operator== (const G4HETCFragment &right) const =delete
G4bool	operator!= (const G4HETCFragment &right) const =delete
Public Member Functions inherited from G4VPreCompoundFragment
	G4VPreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *)
virtual	~G4VPreCompoundFragment ()
void	Initialize (const G4Fragment &aFragment)
G4bool	IsItPossible (const G4Fragment &aFragment) const
G4ReactionProduct *	GetReactionProduct () const
G4int	GetA () const
G4int	GetZ () const
G4int	GetRestA () const
G4int	GetRestZ () const
G4double	GetBindingEnergy () const
G4double	GetEnergyThreshold () const
G4double	GetEmissionProbability () const
G4double	GetNuclearMass () const
G4double	GetRestNuclearMass () const
const G4LorentzVector &	GetMomentum () const
void	SetMomentum (const G4LorentzVector &lv)
void	SetOPTxs (G4int opt)
void	UseSICB (G4bool use)
	G4VPreCompoundFragment (const G4VPreCompoundFragment &right)=delete
const G4VPreCompoundFragment &	operator= (const G4VPreCompoundFragment &right)=delete
G4bool	operator== (const G4VPreCompoundFragment &right) const =delete
G4bool	operator!= (const G4VPreCompoundFragment &right) const =delete

Protected Member Functions
G4double	GetAlpha () const override
G4double	GetBeta () const override
G4double	GetSpinFactor () const override
G4double	K (const G4Fragment &aFragment) const override
Protected Member Functions inherited from G4HETCFragment
G4double	BetaRand (G4int N, G4int L) const
Protected Member Functions inherited from G4VPreCompoundFragment

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundFragment
G4NuclearLevelData *	fNucData
G4DeexPrecoParameters *	theParameters
G4Pow *	g4calc
G4int	theA
G4int	theZ
G4int	theResA {0}
G4int	theResZ {0}
G4int	theFragA {0}
G4int	theFragZ {0}
G4double	theResA13 {0.0}
G4double	theBindingEnergy {0.0}
G4double	theMinKinEnergy {0.0}
G4double	theMaxKinEnergy {0.0}
G4double	theResMass {0.0}
G4double	theReducedMass {0.0}
G4double	theMass
G4double	theEmissionProbability {0.0}
G4double	theCoulombBarrier {0.0}
G4int	OPTxs {3}
G4bool	useSICB {true}

Detailed Description

Definition at line 39 of file G4HETCNeutron.hh.

Constructor & Destructor Documentation

G4HETCNeutron() [1/2]

G4HETCNeutron::G4HETCNeutron

(

)

Definition at line 39 of file G4HETCNeutron.cc.

  : G4HETCFragment(G4Neutron::Neutron())
{}

~G4HETCNeutron()

G4HETCNeutron::~G4HETCNeutron ( )

overridedefault

G4HETCNeutron() [2/2]

G4HETCNeutron::G4HETCNeutron ( const G4HETCNeutron & right )

delete

Member Function Documentation

GetAlpha()

G4double G4HETCNeutron::GetAlpha ( ) const

overrideprotectedvirtual

Implements G4VPreCompoundFragment.

Definition at line 43 of file G4HETCNeutron.cc.

{
  return 0.76+2.2/theResA13;
}

Referenced by GetBeta().

GetBeta()

G4double G4HETCNeutron::GetBeta ( ) const

overrideprotectedvirtual

Reimplemented from G4VPreCompoundFragment.

Definition at line 48 of file G4HETCNeutron.cc.

{
  return (2.12/(theResA13*theResA13)-0.05)/GetAlpha(); // in MeV
}

Referenced by SampleKineticEnergy().

GetSpinFactor()

G4double G4HETCNeutron::GetSpinFactor ( ) const

overrideprotectedvirtual

Implements G4HETCFragment.

Definition at line 53 of file G4HETCNeutron.cc.

{
  // (2s+1)
  return 2.0;
}

K()

G4double G4HETCNeutron::K ( const G4Fragment & aFragment ) const

overrideprotectedvirtual

Implements G4HETCFragment.

Definition at line 59 of file G4HETCNeutron.cc.

{
  // Number of protons in emitted fragment
  G4int Pa = theZ;
  // Number of neutrons in emitted fragment 
  G4int Na = theA - Pa;
 
  G4double r = G4double(theResZ)/G4double(theResA);
  
  G4int P = aFragment.GetNumberOfParticles();
  G4int H = aFragment.GetNumberOfHoles();
  
  G4double result = (P > 0) ? (H + Na/(1.0 - r))/P : 0.0;
  return result;
}

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

SampleKineticEnergy()

G4double G4HETCNeutron::SampleKineticEnergy ( const G4Fragment & frag )

overridevirtual

Implements G4VPreCompoundFragment.

Definition at line 75 of file G4HETCNeutron.cc.

{
  G4int H = frag.GetNumberOfHoles();
  G4int Pb = frag.GetNumberOfParticles();
  G4int Nb = Pb + H;
  G4double U = frag.GetExcitationEnergy();
  G4double g0 = 
    (6.0/CLHEP::pi2)*fNucData->GetLevelDensity(theFragZ,theFragA,U);
  
  G4double Ab = std::max(0.0, (G4double)(Pb*Pb+H*H+Pb-3*H)/(4.0*g0));
  G4double Emax = theMaxKinEnergy - Ab;
  
  G4double cut = GetBeta() / (GetBeta()+Emax/G4double(Nb+1));
  G4double x = (G4UniformRand() <= cut) ? BetaRand(Nb,1) : BetaRand(Nb,2);
 
  return Emax * (1.0 - x);
}

---

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

• G4HETCNeutron.hh
• G4HETCNeutron.cc