G4PreCompoundFragment Class Reference

#include <G4PreCompoundFragment.hh>

Inheritance diagram for G4PreCompoundFragment:

Public Member Functions
	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)
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	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 44 of file G4PreCompoundFragment.hh.

Constructor & Destructor Documentation

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

Definition at line 40 of file G4PreCompoundFragment.cc.

  : G4VPreCompoundFragment(part,aCoulombBarrier)
{}

G4PreCompoundFragment::~G4PreCompoundFragment ( )

virtual

Definition at line 45 of file G4PreCompoundFragment.cc.

{}

Member Function Documentation

G4double G4PreCompoundFragment::CalcEmissionProbability ( const G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundFragment.

Definition at line 49 of file G4PreCompoundFragment.cc.

References G4VPreCompoundFragment::GetMaximalKineticEnergy(), G4VPreCompoundFragment::OPTxs, G4VPreCompoundFragment::theCoulombBarrier, G4VPreCompoundFragment::theEmissionProbability, and G4VPreCompoundFragment::useSICB.

{
  //G4cout << theCoulombBarrier << "  " << GetMaximalKineticEnergy() << G4endl;
  // If  theCoulombBarrier effect is included in the emission probabilities
  //if (GetMaximalKineticEnergy() <= 0.0) 
  G4double limit = 0.0;
  if(OPTxs==0 ||  useSICB) { limit = theCoulombBarrier; }
  if (GetMaximalKineticEnergy() <= limit) 
    {
      theEmissionProbability = 0.0;
      return 0.0;
    }    
  // If  theCoulombBarrier effect is included in the emission probabilities
  //  G4double LowerLimit = 0.;
  // Coulomb barrier is the lower limit 
  // of integration over kinetic energy
  G4double LowerLimit = limit;

  // Excitation energy of nucleus after fragment emission is the upper 
  //limit of integration over kinetic energy
  G4double UpperLimit = GetMaximalKineticEnergy();
  
  theEmissionProbability = 
    IntegrateEmissionProbability(LowerLimit,UpperLimit,aFragment);
  /*
  G4cout << "## G4PreCompoundFragment::CalcEmisProb "
         << "Z= " << aFragment.GetZ_asInt() 
     << " A= " << aFragment.GetA_asInt()
     << " Elow= " << LowerLimit/MeV
     << " Eup= " << UpperLimit/MeV
     << " prob= " << theEmissionProbability
     << G4endl;
  */
  return theEmissionProbability;
}

G4double G4PreCompoundFragment::GetKineticEnergy ( const G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundFragment.

Definition at line 128 of file G4PreCompoundFragment.cc.

References G4UniformRand, G4VPreCompoundFragment::GetEmissionProbability(), G4VPreCompoundFragment::GetMaximalKineticEnergy(), G4VPreCompoundFragment::OPTxs, ProbabilityDistributionFunction(), G4VPreCompoundFragment::theCoulombBarrier, and G4VPreCompoundFragment::useSICB.

{
  //let's keep this way for consistency with CalcEmissionProbability method
  G4double V = 0.0;
  if(OPTxs==0 || useSICB) { V = theCoulombBarrier; }

  G4double Tmax = GetMaximalKineticEnergy();  
  if(Tmax < V) { return 0.0; }
  G4double T(0.0);
  G4double Probability(1.0);
  G4double maxProbability = GetEmissionProbability();
  do 
    {
      T = V + G4UniformRand()*(Tmax-V);
      Probability = ProbabilityDistributionFunction(T,aFragment);      
    }   while (maxProbability*G4UniformRand() > Probability);  
  return T;
}

virtual G4double G4PreCompoundFragment::ProbabilityDistributionFunction ( G4double  K, const G4Fragment &  aFragment  )

protectedpure virtual

Implemented in G4PreCompoundIon, and G4PreCompoundNucleon.

Referenced by GetKineticEnergy().

---

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

• G4PreCompoundFragment.hh
• G4PreCompoundFragment.cc