G4VPreCompoundFragment Class Reference

#include <G4VPreCompoundFragment.hh>

Inheritance diagram for G4VPreCompoundFragment:

Public Member Functions
	G4VPreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *)
virtual	~G4VPreCompoundFragment ()
G4bool	Initialize (const G4Fragment &aFragment)
virtual G4double	CalcEmissionProbability (const G4Fragment &)=0
virtual G4double	SampleKineticEnergy (const G4Fragment &)=0
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)
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
virtual G4double	GetAlpha () const =0
virtual G4double	GetBeta () const

Protected Attributes
G4NuclearLevelData *	fNucData
G4DeexPrecoParameters *	theParameters
G4Pow *	g4calc
G4InterfaceToXS *	fXSection {nullptr}
G4int	theA
G4int	theZ
G4int	theResA {0}
G4int	theResZ {0}
G4int	theFragA {0}
G4int	theFragZ {0}
G4int	OPTxs
G4int	index {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}
G4bool	useSICB {true}

Friends
std::ostream &	operator<< (std::ostream &out, const G4VPreCompoundFragment *theFragment)
std::ostream &	operator<< (std::ostream &out, const G4VPreCompoundFragment &theFragment)

Detailed Description

Definition at line 56 of file G4VPreCompoundFragment.hh.

Constructor & Destructor Documentation

G4VPreCompoundFragment() [1/2]

G4VPreCompoundFragment::G4VPreCompoundFragment ( const G4ParticleDefinition * part, G4VCoulombBarrier * aCoulombBarrier )

explicit

Definition at line 41 of file G4VPreCompoundFragment.cc.

  : theA(part->GetBaryonNumber()),
    theZ(G4lrint(part->GetPDGCharge()/CLHEP::eplus)),
    particle(part),
    theCoulombBarrierPtr(aCoulombBarrier)
{
  theMass = particle->GetPDGMass();
  fNucData = G4NuclearLevelData::GetInstance();
  theParameters = fNucData->GetParameters();
  OPTxs = theParameters->GetDeexModelType();
  g4calc = G4Pow::GetInstance();
 
  if (1 == theZ && 1 == theA) { index = 1; }
  else if (1 == theZ && 2 == theA) { index = 2; }
  else if (1 == theZ && 3 == theA) { index = 3; }
  else if (2 == theZ && 3 == theA) { index = 4; }
  else if (2 == theZ && 4 == theA) { index = 5; }
 
  if (OPTxs == 1) {
    fXSection = new G4InterfaceToXS(particle, index);
  }
}

Referenced by G4HETCFragment::G4HETCFragment(), G4PreCompoundFragment::G4PreCompoundFragment(), G4VPreCompoundFragment(), operator!=(), operator<<, operator<<, operator=(), and operator==().

~G4VPreCompoundFragment()

G4VPreCompoundFragment::~G4VPreCompoundFragment ( )

virtual

Definition at line 65 of file G4VPreCompoundFragment.cc.

{
  delete theCoulombBarrierPtr;
  delete fXSection;
}

G4VPreCompoundFragment() [2/2]

G4VPreCompoundFragment::G4VPreCompoundFragment ( const G4VPreCompoundFragment & right )

delete

Member Function Documentation

CalcEmissionProbability()

virtual G4double G4VPreCompoundFragment::CalcEmissionProbability ( const G4Fragment & )

pure virtual

Implemented in G4HETCFragment, and G4PreCompoundFragment.

GetA()

G4int G4VPreCompoundFragment::GetA ( ) const

inline

Definition at line 89 of file G4VPreCompoundFragment.hh.

{ return theA; }

Referenced by operator<<, G4PreCompoundEmission::PerformEmission(), and G4PreCompoundIon::ProbabilityDistributionFunction().

GetAlpha()

virtual G4double G4VPreCompoundFragment::GetAlpha ( ) const

protectedpure virtual

Implemented in G4HETCAlpha, G4HETCDeuteron, G4HETCHe3, G4HETCNeutron, G4HETCProton, G4HETCTriton, G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, G4PreCompoundNeutron, G4PreCompoundProton, and G4PreCompoundTriton.

GetBeta()

virtual G4double G4VPreCompoundFragment::GetBeta ( ) const

inlineprotectedvirtual

Reimplemented in G4HETCNeutron, and G4PreCompoundNeutron.

Definition at line 129 of file G4VPreCompoundFragment.hh.

{ return -theCoulombBarrier; }

GetBindingEnergy()

G4double G4VPreCompoundFragment::GetBindingEnergy ( ) const

inline

Definition at line 97 of file G4VPreCompoundFragment.hh.

{ return theBindingEnergy; }

GetEmissionProbability()

G4double G4VPreCompoundFragment::GetEmissionProbability ( ) const

inline

Definition at line 104 of file G4VPreCompoundFragment.hh.

{ return theEmissionProbability; }

GetEnergyThreshold()

G4double G4VPreCompoundFragment::GetEnergyThreshold ( ) const

inline

Definition at line 99 of file G4VPreCompoundFragment.hh.

  {
    return theMaxKinEnergy - theCoulombBarrier;
  }

Referenced by G4HETCFragment::CalcEmissionProbability().

GetMomentum()

const G4LorentzVector & G4VPreCompoundFragment::GetMomentum ( ) const

inline

Definition at line 110 of file G4VPreCompoundFragment.hh.

{ return theMomentum; }

Referenced by GetReactionProduct().

GetNuclearMass()

G4double G4VPreCompoundFragment::GetNuclearMass ( ) const

inline

Definition at line 106 of file G4VPreCompoundFragment.hh.

{ return theMass; }

Referenced by operator<<, and G4PreCompoundEmission::PerformEmission().

GetReactionProduct()

G4ReactionProduct * G4VPreCompoundFragment::GetReactionProduct ( ) const

inline

Definition at line 167 of file G4VPreCompoundFragment.hh.

{
  G4ReactionProduct* theReactionProduct = new G4ReactionProduct(particle);
  theReactionProduct->SetMomentum(GetMomentum().vect());
  theReactionProduct->SetTotalEnergy(GetMomentum().e());
  return theReactionProduct;
}

Referenced by G4PreCompoundEmission::PerformEmission().

GetRestA()

G4int G4VPreCompoundFragment::GetRestA ( ) const

inline

Definition at line 93 of file G4VPreCompoundFragment.hh.

{ return theResA; }

Referenced by G4PreCompoundEmission::PerformEmission().

GetRestNuclearMass()

G4double G4VPreCompoundFragment::GetRestNuclearMass ( ) const

inline

Definition at line 108 of file G4VPreCompoundFragment.hh.

{ return theResMass; }

GetRestZ()

G4int G4VPreCompoundFragment::GetRestZ ( ) const

inline

Definition at line 95 of file G4VPreCompoundFragment.hh.

{ return theResZ; }

Referenced by G4PreCompoundEmission::PerformEmission().

GetZ()

G4int G4VPreCompoundFragment::GetZ ( ) const

inline

Definition at line 91 of file G4VPreCompoundFragment.hh.

{ return theZ; }

Referenced by operator<<, and G4PreCompoundEmission::PerformEmission().

Initialize()

G4bool G4VPreCompoundFragment::Initialize

(

const G4Fragment &

aFragment

)

Definition at line 89 of file G4VPreCompoundFragment.cc.

{
  theFragA = aFragment.GetA_asInt();
  theFragZ = aFragment.GetZ_asInt();
  theResA = theFragA - theA;
  theResZ = theFragZ - theZ;
 
  theMinKinEnergy = theMaxKinEnergy = theCoulombBarrier = 0.0;
  if ((theResA < theResZ) || (theResA < theA) || (theResZ < theZ)
      || (theResA == theA && theResZ < theZ)
      || ((theResA > 1) && (theResA == theResZ || theResZ == 0))) {
    return false;
  }
  theResMass = G4NucleiProperties::GetNuclearMass(theResA, theResZ);
  G4double Ecm = aFragment.GetMomentum().m();
  if (Ecm <= theResMass + theMass) { return 0.0; }
 
  theResA13 = g4calc->Z13(theResA);
 
  G4double elim = 0.0;
  if (0 < theZ) {
    theCoulombBarrier = theCoulombBarrierPtr->
      GetCoulombBarrier(theResA, theResZ, aFragment.GetExcitationEnergy());
    elim = (0 < OPTxs) ? theCoulombBarrier*0.5 : theCoulombBarrier;
  }
      
  // Compute Maximal Kinetic Energy which can be carried by fragments 
  // after separation - the true assimptotic value
  theMaxKinEnergy =
    0.5*((Ecm - theResMass)*(Ecm + theResMass) + theMass*theMass)/Ecm - theMass;
  G4double resM = Ecm - theMass - elim;
  if (resM < theResMass) { return false; }
  theMinKinEnergy =
    0.5*((Ecm - resM)*(Ecm + resM) + theMass*theMass)/Ecm - theMass;
 
  if (theMinKinEnergy >= theMaxKinEnergy) { return false; }
  // Calculate masses
  theReducedMass = theResMass*theMass/(theResMass + theMass);
 
  // Compute Binding Energies for fragments 
  // needed to separate a fragment from the nucleus
  theBindingEnergy = theResMass + theMass - aFragment.GetGroundStateMass();
  return true;
}

Referenced by G4PreCompoundFragment::CalcEmissionProbability().

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

SampleKineticEnergy()

virtual G4double G4VPreCompoundFragment::SampleKineticEnergy ( const G4Fragment & )

pure virtual

Implemented in G4HETCChargedFragment, G4HETCNeutron, and G4PreCompoundFragment.

Referenced by G4PreCompoundEmission::PerformEmission().

SetMomentum()

void G4VPreCompoundFragment::SetMomentum ( const G4LorentzVector & lv )

inline

Definition at line 112 of file G4VPreCompoundFragment.hh.

{ theMomentum = lv; }

Referenced by G4PreCompoundEmission::PerformEmission().

SetOPTxs()

void G4VPreCompoundFragment::SetOPTxs ( G4int )

inline

Definition at line 115 of file G4VPreCompoundFragment.hh.

{}

UseSICB()

void G4VPreCompoundFragment::UseSICB ( G4bool use )

inline

Definition at line 117 of file G4VPreCompoundFragment.hh.

{ useSICB = use; } 

Friends And Related Symbol Documentation

operator<< [1/2]

std::ostream & operator<< ( std::ostream & out, const G4VPreCompoundFragment & theFragment )

friend

Definition at line 71 of file G4VPreCompoundFragment.cc.

{
  out << &theFragment;
  return out; 
}

operator<< [2/2]

std::ostream & operator<< ( std::ostream & out, const G4VPreCompoundFragment * theFragment )

friend

Definition at line 78 of file G4VPreCompoundFragment.cc.

{
  out 
    << "PreCompoundModel Emitted Fragment: Z= " << theFragment->GetZ() 
    << " A= " << theFragment->GetA()
    << " Mass(GeV)= " << theFragment->GetNuclearMass()/CLHEP::GeV;
  return out;
}

Member Data Documentation

fNucData

G4NuclearLevelData* G4VPreCompoundFragment::fNucData

protected

Definition at line 131 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment(), G4PreCompoundIon::ProbabilityDistributionFunction(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), and G4HETCNeutron::SampleKineticEnergy().

fXSection

G4InterfaceToXS* G4VPreCompoundFragment::fXSection {nullptr}

protected

Definition at line 134 of file G4VPreCompoundFragment.hh.

{nullptr};

Referenced by G4PreCompoundFragment::CrossSection(), G4VPreCompoundFragment(), and ~G4VPreCompoundFragment().

g4calc

G4Pow* G4VPreCompoundFragment::g4calc

protected

Definition at line 133 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment(), Initialize(), G4PreCompoundIon::ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

index

G4int G4VPreCompoundFragment::index {0}

protected

Definition at line 144 of file G4VPreCompoundFragment.hh.

{0};

Referenced by G4PreCompoundFragment::CrossSection(), and G4VPreCompoundFragment().

OPTxs

G4int G4VPreCompoundFragment::OPTxs

protected

Definition at line 143 of file G4VPreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CrossSection(), G4VPreCompoundFragment(), and Initialize().

theA

G4int G4VPreCompoundFragment::theA

protected

Definition at line 136 of file G4VPreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CrossSection(), G4VPreCompoundFragment(), GetA(), Initialize(), G4HETCAlpha::K(), G4HETCDeuteron::K(), G4HETCHe3::K(), G4HETCNeutron::K(), and G4HETCTriton::K().

theBindingEnergy

G4double G4VPreCompoundFragment::theBindingEnergy {0.0}

protected

Definition at line 147 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by GetBindingEnergy(), Initialize(), G4PreCompoundIon::ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

theCoulombBarrier

G4double G4VPreCompoundFragment::theCoulombBarrier {0.0}

protected

Definition at line 155 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by G4HETCFragment::CalcEmissionProbability(), G4PreCompoundFragment::CrossSection(), GetBeta(), GetEnergyThreshold(), Initialize(), and G4HETCChargedFragment::SampleKineticEnergy().

theEmissionProbability

G4double G4VPreCompoundFragment::theEmissionProbability {0.0}

protected

Definition at line 154 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by G4HETCFragment::CalcEmissionProbability(), G4PreCompoundFragment::CalcEmissionProbability(), and GetEmissionProbability().

theFragA

G4int G4VPreCompoundFragment::theFragA {0}

protected

Definition at line 140 of file G4VPreCompoundFragment.hh.

{0};

Referenced by Initialize(), G4PreCompoundIon::ProbabilityDistributionFunction(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), and G4HETCNeutron::SampleKineticEnergy().

theFragZ

G4int G4VPreCompoundFragment::theFragZ {0}

protected

Definition at line 141 of file G4VPreCompoundFragment.hh.

{0};

Referenced by G4HETCAlpha::GetAlpha(), G4HETCDeuteron::GetAlpha(), G4HETCHe3::GetAlpha(), G4HETCTriton::GetAlpha(), G4PreCompoundAlpha::GetAlpha(), G4PreCompoundDeuteron::GetAlpha(), G4PreCompoundHe3::GetAlpha(), G4PreCompoundTriton::GetAlpha(), Initialize(), G4PreCompoundIon::ProbabilityDistributionFunction(), G4PreCompoundNucleon::ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), and G4HETCNeutron::SampleKineticEnergy().

theMass

G4double G4VPreCompoundFragment::theMass

protected

Definition at line 152 of file G4VPreCompoundFragment.hh.

Referenced by G4VPreCompoundFragment(), GetNuclearMass(), and Initialize().

theMaxKinEnergy

G4double G4VPreCompoundFragment::theMaxKinEnergy {0.0}

protected

Definition at line 149 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by G4HETCFragment::CalcEmissionProbability(), G4PreCompoundFragment::CalcEmissionProbability(), GetEnergyThreshold(), Initialize(), G4PreCompoundIon::ProbabilityDistributionFunction(), G4HETCChargedFragment::SampleKineticEnergy(), G4HETCNeutron::SampleKineticEnergy(), and G4PreCompoundFragment::SampleKineticEnergy().

theMinKinEnergy

G4double G4VPreCompoundFragment::theMinKinEnergy {0.0}

protected

Definition at line 148 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by G4PreCompoundFragment::CalcEmissionProbability(), Initialize(), and G4PreCompoundFragment::SampleKineticEnergy().

theParameters

G4DeexPrecoParameters* G4VPreCompoundFragment::theParameters

protected

Definition at line 132 of file G4VPreCompoundFragment.hh.

Referenced by G4HETCFragment::G4HETCFragment(), G4PreCompoundIon::G4PreCompoundIon(), and G4VPreCompoundFragment().

theReducedMass

G4double G4VPreCompoundFragment::theReducedMass {0.0}

protected

Definition at line 151 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by Initialize(), G4PreCompoundIon::ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

theResA

G4int G4VPreCompoundFragment::theResA {0}

protected

Definition at line 138 of file G4VPreCompoundFragment.hh.

{0};

Referenced by G4PreCompoundFragment::CrossSection(), GetRestA(), Initialize(), G4HETCAlpha::K(), G4HETCDeuteron::K(), G4HETCHe3::K(), G4HETCNeutron::K(), G4HETCProton::K(), G4HETCTriton::K(), G4PreCompoundIon::ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

theResA13

G4double G4VPreCompoundFragment::theResA13 {0.0}

protected

Definition at line 146 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by G4PreCompoundFragment::CrossSection(), G4HETCNeutron::GetAlpha(), G4PreCompoundNeutron::GetAlpha(), G4HETCNeutron::GetBeta(), G4PreCompoundNeutron::GetBeta(), and Initialize().

theResMass

G4double G4VPreCompoundFragment::theResMass {0.0}

protected

Definition at line 150 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by GetRestNuclearMass(), and Initialize().

theResZ

G4int G4VPreCompoundFragment::theResZ {0}

protected

Definition at line 139 of file G4VPreCompoundFragment.hh.

{0};

Referenced by G4PreCompoundFragment::CrossSection(), G4HETCProton::GetAlpha(), G4PreCompoundProton::GetAlpha(), GetRestZ(), Initialize(), G4HETCAlpha::K(), G4HETCDeuteron::K(), G4HETCHe3::K(), G4HETCNeutron::K(), G4HETCProton::K(), G4HETCTriton::K(), G4PreCompoundIon::ProbabilityDistributionFunction(), and G4PreCompoundNucleon::ProbabilityDistributionFunction().

theZ

G4int G4VPreCompoundFragment::theZ

protected

Definition at line 137 of file G4VPreCompoundFragment.hh.

Referenced by G4PreCompoundFragment::CrossSection(), G4VPreCompoundFragment(), GetZ(), Initialize(), G4HETCAlpha::K(), G4HETCDeuteron::K(), G4HETCHe3::K(), G4HETCNeutron::K(), G4HETCProton::K(), and G4HETCTriton::K().

useSICB

G4bool G4VPreCompoundFragment::useSICB {true}

protected

Definition at line 158 of file G4VPreCompoundFragment.hh.

{true};

Referenced by UseSICB().

---

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

• G4VPreCompoundFragment.hh
• G4VPreCompoundFragment.cc