54 EmissionProbability(0.0),
56 theEvaporationProbabilityPtr(aEmissionStrategy),
101 G4double Etot = FragmentMass + ExEnergy;
144 ResidualMomentum -= EvaporatedMomentum;
178 G4double TCN = 1.0/(std::sqrt(aCN/UxCN) - 1.5/UxCN);
185 - 1.25*
G4Log(UxCN/
MeV) + 2.0*std::sqrt(aCN*UxCN));
186 InitialLevelDensity = (
pi/12.0)*
G4Exp((U-E0CN)/TCN)/TCN;
192 InitialLevelDensity = (
pi/12.0)*
G4Exp(2*x1)/(x*std::sqrt(x1));
207 Rb = (1.12*(Aj + Ad) - 0.86*((Aj+Ad)/(Aj*Ad))+2.85)*
fermi;
213 Rb=1.5*(Aj+Ad)*
fermi;
222 G4double ConstantFactor = gg*GeometricalXS*
Alpha*
pi/(InitialLevelDensity*12);
230 for(
G4int i=0; i<100; ++i) {
232 G4double edelta = theEnergy-KineticEnergy-delta0;
233 Probability = ConstantFactor*(KineticEnergy + Beta);
236 G4double T = 1.0/(std::sqrt(a/Ux) - 1.5/Ux);
239 if (theEnergy - KineticEnergy < Ex) {
241 - 1.25*
G4Log(Ux) + 2.0*std::sqrt(a*Ux));
242 Probability *=
G4Exp((theEnergy-KineticEnergy-E0)/T)/T;
251 return KineticEnergy;
static const G4double e2[44]
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double G4Log(G4double x)
G4ThreeVector G4RandomDirection()
static constexpr double pi2
static constexpr double fermi
static constexpr double GeV
static constexpr double MeV
static constexpr double pi
Hep3Vector boostVector() const
HepLorentzVector & boost(double, double, double)
G4double GetGroundStateMass() const
G4double GetExcitationEnergy() const
const G4LorentzVector & GetMomentum() const
void SetCreatorModelID(G4int value)
void SetMomentum(const G4LorentzVector &value)
void SetZandA_asInt(G4int Znew, G4int Anew)
virtual void Dump() const
G4double MaximalKineticEnergy
G4VCoulombBarrier * theCoulombBarrierPtr
G4VLevelDensityParameter * theLevelDensityPtr
G4GEMProbability * theEvaporationProbabilityPtr
G4NuclearLevelData * fNucData
virtual G4Fragment * EmittedFragment(G4Fragment *theNucleus)
G4GEMChannel(G4int theA, G4int theZ, const G4String &aName, G4GEMProbability *aEmissionStrategy)
G4double EmissionProbability
virtual G4double GetEmissionProbability(G4Fragment *theNucleus)
G4double SampleKineticEnergy(const G4Fragment &fragment)
G4double CalcAlphaParam(const G4Fragment &) const
G4double CalcBetaParam(const G4Fragment &) const
G4double GetSpin(void) const
void SetCoulomBarrier(const G4VCoulombBarrier *aCoulombBarrierStrategy)
G4PairingCorrection * GetPairingCorrection()
static G4NuclearLevelData * GetInstance()
static G4double GetNuclearMass(const G4double A, const G4double Z)
static G4int GetModelID(const G4int modelIndex)
static G4Pow * GetInstance()
G4double Z13(G4int Z) const
virtual G4double GetCoulombBarrier(G4int ARes, G4int ZRes, G4double U) const =0
virtual G4double LevelDensityParameter(G4int A, G4int Z, G4double U) const =0