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	SampleKineticEnergy (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	SampleKineticEnergy (const G4Fragment &aFragment)=0
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 &value)
void	SetOPTxs (G4int)
void	UseSICB (G4bool)

Protected Member Functions
virtual G4double	GetAlpha () const =0
virtual G4double	GetBeta () const =0
G4double	CrossSection (G4double ekin) const
virtual G4double	ProbabilityDistributionFunction (G4double K, const G4Fragment &aFragment)=0
virtual G4double	GetAlpha () const =0
virtual G4double	GetBeta () const =0

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

Detailed Description

Definition at line 43 of file G4PreCompoundFragment.hh.

Constructor & Destructor Documentation

G4PreCompoundFragment()

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

Definition at line 42 of file G4PreCompoundFragment.cc.

  : G4VPreCompoundFragment(p, aCoulBarrier)
{
  muu = probmax = 0.0;
  if(0 == theZ)      { index = 0; }
  else if(1 == theZ) { index = theA; }
  else               { index = theA + 1; }
}

~G4PreCompoundFragment()

G4PreCompoundFragment::~G4PreCompoundFragment ( )

virtual

Definition at line 52 of file G4PreCompoundFragment.cc.

{}

Member Function Documentation

CalcEmissionProbability()

G4double G4PreCompoundFragment::CalcEmissionProbability ( const G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundFragment.

Definition at line 55 of file G4PreCompoundFragment.cc.

{
  //G4cout << theCoulombBarrier << "  " << GetMaximalKineticEnergy() << G4endl;
  // If  theCoulombBarrier effect is included in the emission probabilities
  // Coulomb barrier is the lower limit of integration over kinetic energy
 
  theEmissionProbability = 0.0;
 
  if (theMaxKinEnergy <= theMinKinEnergy) { return 0.0; }    
 
  // compute power once
  if(0 < index) { 
    muu = G4KalbachCrossSection::ComputePowerParameter(theResA, index);
  }
  
  theEmissionProbability = 
    IntegrateEmissionProbability(theMinKinEnergy,theMaxKinEnergy,aFragment);
  /*
  G4cout << "## G4PreCompoundFragment::CalcEmisProb "
         << "Z= " << aFragment.GetZ_asInt() 
     << " A= " << aFragment.GetA_asInt()
     << " Elow= " << LowerLimit/MeV
     << " Eup= " << UpperLimit/MeV
     << " prob= " << theEmissionProbability
     << G4endl;
  */
  return theEmissionProbability;
}

CrossSection()

G4double G4PreCompoundFragment::CrossSection ( G4double  ekin ) const

protected

Definition at line 112 of file G4PreCompoundFragment.cc.

{
  G4double res;
  if(OPTxs == 0 || (OPTxs == 4 && theMaxKinEnergy < 10.)) { 
    res = GetOpt0(ekin);
 
  } else if(OPTxs <= 2) { 
    res = G4ChatterjeeCrossSection::ComputeCrossSection(ekin, theCoulombBarrier, 
                            theResA13, muu, 
                            index, theZ, theResA); 
 
  } else { 
    res = G4KalbachCrossSection::ComputeCrossSection(ekin, theCoulombBarrier, 
                             theResA13, muu,
                             index, theZ, theA, theResA);
  }
  return res;
}  

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

GetAlpha()

virtual G4double G4PreCompoundFragment::GetAlpha ( ) const

protectedpure virtual

Implements G4VPreCompoundFragment.

Implemented in G4PreCompoundAlpha, G4PreCompoundDeuteron, G4PreCompoundHe3, G4PreCompoundNeutron, G4PreCompoundProton, and G4PreCompoundTriton.

GetBeta()

virtual G4double G4PreCompoundFragment::GetBeta ( ) const

protectedpure virtual

Implements G4VPreCompoundFragment.

Implemented in G4PreCompoundIon, G4PreCompoundNeutron, and G4PreCompoundProton.

ProbabilityDistributionFunction()

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

protectedpure virtual

Implemented in G4PreCompoundIon, and G4PreCompoundNucleon.

Referenced by SampleKineticEnergy().

SampleKineticEnergy()

G4double G4PreCompoundFragment::SampleKineticEnergy ( const G4Fragment &  aFragment )

virtual

Implements G4VPreCompoundFragment.

Definition at line 140 of file G4PreCompoundFragment.cc.

{
  G4double delta = theMaxKinEnergy - theMinKinEnergy;
  static const G4double toler = 1.25;
  probmax *= toler;
  G4double prob, T(0.0);
  CLHEP::HepRandomEngine* rndm = G4Random::getTheEngine();
  G4int i;
  for(i=0; i<100; ++i) {
    T = theMinKinEnergy + delta*rndm->flat();
    prob = ProbabilityDistributionFunction(T, fragment);
    /*
    if(prob > probmax) { 
      G4cout << "G4PreCompoundFragment WARNING: prob= " << prob 
         << " probmax= " << probmax << G4endl;
      G4cout << "i= " << i << " Z= " << theZ << " A= " << theA 
         << " resZ= " << theResZ << " resA= " << theResA << "\n"
         << " T= " << T << " Tmax= " << theMaxKinEnergy 
         << " Tmin= " << limit
         << G4endl;
      for(G4int i=0; i<N; ++i) { G4cout << " " << probability[i]; }
      G4cout << G4endl; 
    }
    */
    // Loop checking, 05-Aug-2015, Vladimir Ivanchenko
    if(probmax*rndm->flat() <= prob) { break; }
  }
  /*
  G4cout << "G4PreCompoundFragment: i= " << i << " Z= " << theZ << " A= " << theA 
     <<"  T(MeV)= " << T << " Emin(MeV)= " << theMinKinEnergy << " Emax= " 
     << theMaxKinEnergy << G4endl;
  */
  return T;
}

---

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

• G4PreCompoundFragment.hh
• G4PreCompoundFragment.cc