53 theEvaporationProbabilityPtr(aEmissionStrategy),
54 EmissionProbability(0.0),
98 G4double Etot = FragmentMass + ExEnergy;
137 EvaporatedMomentum.boost(ResidualMomentum.boostVector());
140 ResidualMomentum -= EvaporatedMomentum;
173 G4double TCN = 1.0/(std::sqrt(aCN/UxCN) - 1.5/UxCN);
180 - 1.25*
G4Log(UxCN/MeV) + 2.0*std::sqrt(aCN*UxCN));
181 InitialLevelDensity = (
pi/12.0)*
G4Exp((U-E0CN)/TCN)/TCN;
187 InitialLevelDensity = (
pi/12.0)*
G4Exp(2*x1)/(x*std::sqrt(x1));
203 Rb = (1.12*(Aj + Ad) - 0.86*((Aj+Ad)/(Aj*Ad))+2.85)*
fermi;
209 Rb=1.5*(Aj+Ad)*
fermi;
218 G4double ConstantFactor = gg*GeometricalXS*Alpha*
pi/(InitialLevelDensity*12);
226 for(
G4int i=0; i<100; ++i) {
228 G4double edelta = theEnergy-KineticEnergy-delta0;
229 Probability = ConstantFactor*(KineticEnergy + Beta);
232 G4double T = 1.0/(std::sqrt(a/Ux) - 1.5/Ux);
235 if (theEnergy - KineticEnergy < Ex) {
237 - 1.25*
G4Log(Ux) + 2.0*std::sqrt(a*Ux));
238 Probability *=
G4Exp((theEnergy-KineticEnergy-E0)/T)/T;
242 G4Exp(2*std::sqrt(a*edelta) - 0.25*
G4Log(a*edelta*e2*e2));
247 return KineticEnergy;
255 G4double SinTheta = std::sqrt(1.0 - CosTheta*CosTheta);
258 Magnitude*std::sin(Phi)*SinTheta,
static G4Pow * GetInstance()
G4ThreeVector IsotropicVector(G4double Magnitude=1.0)
static G4double GetNuclearMass(const G4double A, const G4double Z)
CLHEP::Hep3Vector G4ThreeVector
std::vector< ExP01TrackerHit * > a
void SetCoulomBarrier(const G4VCoulombBarrier *aCoulombBarrierStrategy)
G4double CalcAlphaParam(const G4Fragment &) const
G4GEMChannel(G4int theA, G4int theZ, const G4String &aName, G4GEMProbability *aEmissionStrategy)
virtual G4Fragment * EmittedFragment(G4Fragment *theNucleus)
G4double SampleKineticEnergy(const G4Fragment &fragment)
static constexpr double twopi
G4VCoulombBarrier * theCoulombBarrierPtr
G4double Z13(G4int Z) const
double A(double temperature)
G4double GetSpin(void) const
G4PairingCorrection * pairingCorrection
const G4LorentzVector & GetMomentum() const
void SetMomentum(const G4LorentzVector &value)
G4double GetPairingCorrection(G4int A, G4int Z) const
G4double GetGroundStateMass() const
G4double EmissionProbability
G4double CalcBetaParam(const G4Fragment &) const
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
static G4PairingCorrection * GetInstance()
virtual G4double GetEmissionProbability(G4Fragment *theNucleus)
virtual void Dump() const
void SetZandA_asInt(G4int Znew, G4int Anew)
G4VLevelDensityParameter * theLevelDensityPtr
virtual G4double LevelDensityParameter(G4int A, G4int Z, G4double U) const =0
static constexpr double GeV
virtual G4double GetCoulombBarrier(G4int ARes, G4int ZRes, G4double U) const =0
static constexpr double MeV
static constexpr double pi
G4GEMProbability * theEvaporationProbabilityPtr
G4double MaximalKineticEnergy
static constexpr double fermi
G4double GetExcitationEnergy() const
static constexpr double pi2
CLHEP::HepLorentzVector G4LorentzVector