#include <G4PreCompoundFragment.hh>

Inheritance diagram for G4PreCompoundFragment:

Collaboration 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)

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

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundFragment
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 44 of file G4PreCompoundFragment.hh.

Constructor & Destructor Documentation

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 ( )

virtual

Definition at line 52 of file G4PreCompoundFragment.cc.

53 {}

Member Function Documentation

G4double G4PreCompoundFragment::CalcEmissionProbability ( const G4Fragment & aFragment )

virtual

Implements G4VPreCompoundFragment.

Definition at line 56 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(OPTxs <= 2) { 
     muu = G4ChatterjeeCrossSection::ComputePowerParameter(theResA, index);
   } else {
     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;
 }

Here is the call graph for this function:

G4double G4PreCompoundFragment::CrossSection ( G4double ekin ) const

protected

Definition at line 116 of file G4PreCompoundFragment.cc.

 {
   G4double res;
   if(OPTxs == 0) { 
     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;
 }  

Here is the call graph for this function:

Here is the caller graph for this function:

virtual G4double G4PreCompoundFragment::GetAlpha ( ) const

protectedpure virtual

Implements G4VPreCompoundFragment.

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

virtual G4double G4PreCompoundFragment::GetBeta ( ) const

protectedpure virtual

Implements G4VPreCompoundFragment.

Implemented in G4PreCompoundIon, G4PreCompoundNeutron, and G4PreCompoundProton.

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

protectedpure virtual

Implemented in G4PreCompoundIon, and G4PreCompoundNucleon.

Here is the caller graph for this function:

G4double G4PreCompoundFragment::SampleKineticEnergy ( const G4Fragment & aFragment )

virtual

Implements G4VPreCompoundFragment.

Definition at line 144 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;
 }

Here is the call graph for this function:

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

source/geant4.10.03.p03/source/processes/hadronic/models/pre_equilibrium/exciton_model/include/G4PreCompoundFragment.hh
source/geant4.10.03.p03/source/processes/hadronic/models/pre_equilibrium/exciton_model/src/G4PreCompoundFragment.cc

Public Member Functions

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation