G4PreCompoundFragment Class Reference

#include <G4PreCompoundFragment.hh>

Inheritance diagram for G4PreCompoundFragment:

Public Member Functions
	G4PreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
	~G4PreCompoundFragment () override=default
G4double	CalcEmissionProbability (const G4Fragment &aFragment) override
G4double	SampleKineticEnergy (const G4Fragment &aFragment) override
G4double	CrossSection (G4double ekin)
G4double	RecentXS () const
	G4PreCompoundFragment (const G4PreCompoundFragment &right)=delete
const G4PreCompoundFragment &	operator= (const G4PreCompoundFragment &right)=delete
G4bool	operator== (const G4PreCompoundFragment &right) const =delete
G4bool	operator!= (const G4PreCompoundFragment &right) const =delete
Public Member Functions inherited from G4VPreCompoundFragment
	G4VPreCompoundFragment (const G4ParticleDefinition *, G4VCoulombBarrier *)
virtual	~G4VPreCompoundFragment ()
G4bool	Initialize (const G4Fragment &aFragment)
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 ekin, const G4Fragment &aFragment)=0
Protected Member Functions inherited from G4VPreCompoundFragment
virtual G4double	GetAlpha () const =0
virtual G4double	GetBeta () const

Additional Inherited Members
Protected Attributes inherited from G4VPreCompoundFragment
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}

Detailed Description

Definition at line 43 of file G4PreCompoundFragment.hh.

Constructor & Destructor Documentation

G4PreCompoundFragment() [1/2]

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

Definition at line 44 of file G4PreCompoundFragment.cc.

  : G4VPreCompoundFragment(p, aCoulBarrier)
{}

Referenced by G4PreCompoundFragment(), G4PreCompoundIon::G4PreCompoundIon(), G4PreCompoundNucleon::G4PreCompoundNucleon(), operator!=(), operator=(), and operator==().

~G4PreCompoundFragment()

G4PreCompoundFragment::~G4PreCompoundFragment ( )

overridedefault

G4PreCompoundFragment() [2/2]

G4PreCompoundFragment::G4PreCompoundFragment ( const G4PreCompoundFragment & right )

delete

Member Function Documentation

CalcEmissionProbability()

G4double G4PreCompoundFragment::CalcEmissionProbability ( const G4Fragment & aFragment )

overridevirtual

Implements G4VPreCompoundFragment.

Definition at line 49 of file G4PreCompoundFragment.cc.

{
  theEmissionProbability = (Initialize(fr)) ?
    IntegrateEmissionProbability(theMinKinEnergy, theMaxKinEnergy, fr) : 0.0;
  /*  
  G4cout << "## G4PreCompoundFragment::CalcEmisProb "
         << "Zf= " << fr.GetZ_asInt()
     << " Af= " << fr.GetA_asInt()
     << " Elow= " << theMinKinEnergy
     << " Eup= " << theMaxKinEnergy
     << " prob= " << theEmissionProbability
     << " index=" << index << " Z=" << theZ << " A=" << theA
     << G4endl;
  */
  return theEmissionProbability;
}

CrossSection()

G4double G4PreCompoundFragment::CrossSection

(

G4double

ekin

)

Definition at line 94 of file G4PreCompoundFragment.cc.

{
  /*
  G4cout << "G4PreCompoundFragment::CrossSection OPTxs=" << OPTxs << " E=" << ekin
     << " resZ=" << theResZ << " resA=" << theResA << " index=" << index
     << " fXSection:" << fXSection << G4endl;
  */
  // compute power once
  if (OPTxs > 1 && 0 < index && theResA != lastA) {
    lastA = theResA;
    muu = G4KalbachCrossSection::ComputePowerParameter(lastA, index);
  }
  if (OPTxs == 0) { 
    recentXS = GetOpt0(ekin);
  } else if (OPTxs == 1) {
    G4int Z = std::min(theResZ, ZMAXNUCLEARDATA);
    //G4double e = std::max(ekin, lowEnergyLimitMeV[Z]);
    recentXS = fXSection->GetElementCrossSection(ekin, Z)/CLHEP::millibarn;
 
  } else if (OPTxs == 2) { 
    recentXS = G4ChatterjeeCrossSection::ComputeCrossSection(ekin, 
                                                             theCoulombBarrier, 
                                 theResA13, muu, 
                                 index, theZ, theResA); 
 
  } else { 
    recentXS = G4KalbachCrossSection::ComputeCrossSection(ekin, theCoulombBarrier, 
                                  theResA13, muu, index,
                                  theZ, theA, theResA);
  }
  return recentXS;
}  

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

operator!=()

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

delete

operator=()

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

delete

operator==()

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

delete

ProbabilityDistributionFunction()

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

protectedpure virtual

Implemented in G4PreCompoundIon, and G4PreCompoundNucleon.

Referenced by SampleKineticEnergy().

RecentXS()

G4double G4PreCompoundFragment::RecentXS ( ) const

inline

Definition at line 64 of file G4PreCompoundFragment.hh.

{ return recentXS; };

SampleKineticEnergy()

G4double G4PreCompoundFragment::SampleKineticEnergy ( const G4Fragment & aFragment )

overridevirtual

Implements G4VPreCompoundFragment.

Definition at line 136 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