G4VPreCompoundFragment Class Reference

#include <G4VPreCompoundFragment.hh>

Inheritance diagram for G4VPreCompoundFragment:

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

Protected Attributes
G4NuclearLevelData *	fNucData
G4DeexPrecoParameters *	theParameters
G4Pow *	g4calc
G4InterfaceToXS *	fXSection {nullptr}
const G4Fragment *	pFragment {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 57 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);
  }
  InitialiseIntegrator(0.005, 0.25, 1.10, 0.5*CLHEP::MeV,
               0.2*CLHEP::MeV, 5*CLHEP::MeV);
}

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

~G4VPreCompoundFragment()

G4VPreCompoundFragment::~G4VPreCompoundFragment ( )

override

Definition at line 67 of file G4VPreCompoundFragment.cc.

{
  delete theCoulombBarrierPtr;
  delete fXSection;
}

G4VPreCompoundFragment() [2/2]

G4VPreCompoundFragment::G4VPreCompoundFragment ( const G4VPreCompoundFragment & right )

delete

Member Function Documentation

CalcEmissionProbability()

G4double G4VPreCompoundFragment::CalcEmissionProbability ( const G4Fragment & fr )

virtual

Reimplemented in G4HETCFragment.

Definition at line 139 of file G4VPreCompoundFragment.cc.

{
  G4bool ok = Initialize(fr);
  theEmissionProbability = 0.0;
  if (ok) {
    theEmissionProbability = ComputeIntegral(theMinKinEnergy, theMaxKinEnergy);
  }
  return theEmissionProbability;
}

GetA()

G4int G4VPreCompoundFragment::GetA ( ) const

inline

Definition at line 80 of file G4VPreCompoundFragment.hh.

{ return theA; }

Referenced by operator<<, G4PreCompoundEmission::PerformEmission(), G4PreCompoundEmissionInt::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.

Referenced by G4PreCompoundFragment::CrossSection().

GetBeta()

virtual G4double G4VPreCompoundFragment::GetBeta ( ) const

inlineprotectedvirtual

Reimplemented in G4HETCNeutron, and G4PreCompoundNeutron.

Definition at line 129 of file G4VPreCompoundFragment.hh.

{ return -theCoulombBarrier; }

Referenced by G4PreCompoundFragment::CrossSection().

GetBindingEnergy()

G4double G4VPreCompoundFragment::GetBindingEnergy ( ) const

inline

Definition at line 88 of file G4VPreCompoundFragment.hh.

{ return theBindingEnergy; }

GetEmissionProbability()

G4double G4VPreCompoundFragment::GetEmissionProbability ( ) const

inline

Definition at line 95 of file G4VPreCompoundFragment.hh.

{ return theEmissionProbability; }

GetEnergyThreshold()

G4double G4VPreCompoundFragment::GetEnergyThreshold ( ) const

inline

Definition at line 90 of file G4VPreCompoundFragment.hh.

  {
    return theMaxKinEnergy - theCoulombBarrier;
  }

Referenced by G4HETCFragment::CalcEmissionProbability().

GetMomentum()

const G4LorentzVector & G4VPreCompoundFragment::GetMomentum ( ) const

inline

Definition at line 101 of file G4VPreCompoundFragment.hh.

{ return theMomentum; }

Referenced by GetReactionProduct().

GetNuclearMass()

G4double G4VPreCompoundFragment::GetNuclearMass ( ) const

inline

Definition at line 97 of file G4VPreCompoundFragment.hh.

{ return theMass; }

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

GetReactionProduct()

G4ReactionProduct * G4VPreCompoundFragment::GetReactionProduct ( ) const

inline

Definition at line 168 of file G4VPreCompoundFragment.hh.

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

Referenced by G4PreCompoundEmission::PerformEmission(), and G4PreCompoundEmissionInt::PerformEmission().

GetRestA()

G4int G4VPreCompoundFragment::GetRestA ( ) const

inline

Definition at line 84 of file G4VPreCompoundFragment.hh.

{ return theResA; }

Referenced by G4PreCompoundEmission::PerformEmission().

GetRestNuclearMass()

G4double G4VPreCompoundFragment::GetRestNuclearMass ( ) const

inline

Definition at line 99 of file G4VPreCompoundFragment.hh.

{ return theResMass; }

GetRestZ()

G4int G4VPreCompoundFragment::GetRestZ ( ) const

inline

Definition at line 86 of file G4VPreCompoundFragment.hh.

{ return theResZ; }

Referenced by G4PreCompoundEmission::PerformEmission().

GetZ()

G4int G4VPreCompoundFragment::GetZ ( ) const

inline

Definition at line 82 of file G4VPreCompoundFragment.hh.

{ return theZ; }

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

Initialize()

G4bool G4VPreCompoundFragment::Initialize

(

const G4Fragment &

aFragment

)

Definition at line 91 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;
  }
  pFragment = &aFragment;
  theResMass = G4NucleiProperties::GetNuclearMass(theResA, theResZ);
  G4double Ecm = aFragment.GetMomentum().m();
  if (Ecm <= theResMass + theMass) { return 0.0; }
 
  theResA13 = g4calc->Z13(theResA);
  G4int nex = aFragment.GetNumberOfExcitons();
 
  G4double elim = 0.0;
  if (0 < theZ) {
    theCoulombBarrier = theCoulombBarrierPtr->
      GetCoulombBarrier(theResA + nex, 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;
  theMinKinEnergy = std::max(theMinKinEnergy, 0.0);
 
  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 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

ProbabilityDensityFunction()

G4double G4VPreCompoundFragment::ProbabilityDensityFunction ( G4double energy )

overridevirtual

Implements G4VSIntegration.

Definition at line 155 of file G4VPreCompoundFragment.cc.

{
  G4double e = std::max(ekin, 0.02);
  return ProbabilityDistributionFunction(e, *pFragment);
}

ProbabilityDistributionFunction()

virtual G4double G4VPreCompoundFragment::ProbabilityDistributionFunction ( G4double , const G4Fragment &  )

inlineprotectedvirtual

Reimplemented in G4PreCompoundIon, and G4PreCompoundNucleon.

Definition at line 124 of file G4VPreCompoundFragment.hh.

  { return 0.0; }; 

Referenced by ProbabilityDensityFunction().

SampleKineticEnergy()

G4double G4VPreCompoundFragment::SampleKineticEnergy ( const G4Fragment & )

virtual

Reimplemented in G4HETCChargedFragment, and G4HETCNeutron.

Definition at line 149 of file G4VPreCompoundFragment.cc.

{
  G4double ekin = SampleValue();
  return ekin;
}

Referenced by G4PreCompoundEmission::PerformEmission(), and G4PreCompoundEmissionInt::PerformEmission().

SetMomentum()

void G4VPreCompoundFragment::SetMomentum ( const G4LorentzVector & lv )

inline

Definition at line 103 of file G4VPreCompoundFragment.hh.

{ theMomentum = lv; }

Referenced by G4PreCompoundEmission::PerformEmission(), and G4PreCompoundEmissionInt::PerformEmission().

SetOPTxs()

void G4VPreCompoundFragment::SetOPTxs ( G4int )

inline

Definition at line 106 of file G4VPreCompoundFragment.hh.

{}

UseSICB()

void G4VPreCompoundFragment::UseSICB ( G4bool use )

inline

Definition at line 108 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 73 of file G4VPreCompoundFragment.cc.

{
  out << &theFragment;
  return out; 
}

operator<< [2/2]

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

friend

Definition at line 80 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 145 of file G4VPreCompoundFragment.hh.

{0};

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

OPTxs

G4int G4VPreCompoundFragment::OPTxs

protected

Definition at line 144 of file G4VPreCompoundFragment.hh.

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

pFragment

const G4Fragment* G4VPreCompoundFragment::pFragment {nullptr}

protected

Definition at line 135 of file G4VPreCompoundFragment.hh.

{nullptr};

Referenced by G4HETCFragment::CalcEmissionProbability(), Initialize(), and ProbabilityDensityFunction().

theA

G4int G4VPreCompoundFragment::theA

protected

Definition at line 137 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 148 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 156 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 155 of file G4VPreCompoundFragment.hh.

{0.0};

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

theFragA

G4int G4VPreCompoundFragment::theFragA {0}

protected

Definition at line 141 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 142 of file G4VPreCompoundFragment.hh.

{0};

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

theMass

G4double G4VPreCompoundFragment::theMass

protected

Definition at line 153 of file G4VPreCompoundFragment.hh.

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

theMaxKinEnergy

G4double G4VPreCompoundFragment::theMaxKinEnergy {0.0}

protected

Definition at line 150 of file G4VPreCompoundFragment.hh.

{0.0};

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

theMinKinEnergy

G4double G4VPreCompoundFragment::theMinKinEnergy {0.0}

protected

Definition at line 149 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by CalcEmissionProbability(), and Initialize().

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 152 of file G4VPreCompoundFragment.hh.

{0.0};

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

theResA

G4int G4VPreCompoundFragment::theResA {0}

protected

Definition at line 139 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 147 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 151 of file G4VPreCompoundFragment.hh.

{0.0};

Referenced by GetRestNuclearMass(), and Initialize().

theResZ

G4int G4VPreCompoundFragment::theResZ {0}

protected

Definition at line 140 of file G4VPreCompoundFragment.hh.

{0};

Referenced by G4PreCompoundFragment::CrossSection(), G4HETCProton::GetAlpha(), G4PreCompoundAlpha::GetAlpha(), G4PreCompoundDeuteron::GetAlpha(), G4PreCompoundHe3::GetAlpha(), G4PreCompoundProton::GetAlpha(), G4PreCompoundTriton::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 138 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 159 of file G4VPreCompoundFragment.hh.

{true};

Referenced by UseSICB().

---

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

• G4VPreCompoundFragment.hh
• G4VPreCompoundFragment.cc