55 theEvaporationProbabilityPtr(aEmissionStrategy),
56 EmissionProbability(0.0),
57 MaximalKineticEnergy(-CLHEP::
GeV)
101 G4double Etot = FragmentMass + ExEnergy;
143 ResidualMomentum -= EvaporatedMomentum;
176 G4double TCN = 1.0/(std::sqrt(aCN/UxCN) - 1.5/UxCN);
183 - 1.25*
G4Log(UxCN/MeV) + 2.0*std::sqrt(aCN*UxCN));
184 InitialLevelDensity = (
pi/12.0)*
G4Exp((U-E0CN)/TCN)/TCN;
190 InitialLevelDensity = (
pi/12.0)*
G4Exp(2*x1)/(x*std::sqrt(x1));
206 Rb = (1.12*(Aj + Ad) - 0.86*((Aj+Ad)/(Aj*Ad))+2.85)*
fermi;
212 Rb=1.5*(Aj+Ad)*
fermi;
221 G4double ConstantFactor = gg*GeometricalXS*Alpha*
pi/(InitialLevelDensity*12);
229 for(
G4int i=0; i<100; ++i) {
231 G4double edelta = theEnergy-KineticEnergy-delta0;
232 Probability = ConstantFactor*(KineticEnergy + Beta);
235 G4double T = 1.0/(std::sqrt(a/Ux) - 1.5/Ux);
238 if (theEnergy - KineticEnergy < Ex) {
240 - 1.25*
G4Log(Ux) + 2.0*std::sqrt(a*Ux));
241 Probability *=
G4Exp((theEnergy-KineticEnergy-E0)/T)/T;
245 G4Exp(2*std::sqrt(a*edelta) - 0.25*
G4Log(a*edelta*e2*e2));
250 return KineticEnergy;
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
const G4LorentzVector & GetMomentum() const
std::vector< ExP01TrackerHit * > a
G4double Z13(G4int Z) const
G4double MaximalKineticEnergy
static constexpr double MeV
virtual G4double LevelDensityParameter(G4int A, G4int Z, G4double U) const =0
Float_t x1[n_points_granero]
void SetZandA_asInt(G4int Znew, G4int Anew)
static G4NuclearLevelData * GetInstance()
G4PairingCorrection * pairingCorrection
static constexpr double hbarc
G4GEMChannel(G4int theA, G4int theZ, const G4String &aName, G4GEMProbability *aEmissionStrategy)
G4ThreeVector G4RandomDirection()
G4double GetGroundStateMass() const
G4double CalcBetaParam(const G4Fragment &) const
G4double SampleKineticEnergy(const G4Fragment &fragment)
G4double G4Log(G4double x)
G4GEMProbability * theEvaporationProbabilityPtr
static constexpr double fermi
G4double CalcAlphaParam(const G4Fragment &) const
static G4Pow * GetInstance()
virtual void Dump() const
double A(double temperature)
void SetMomentum(const G4LorentzVector &value)
static G4double GetNuclearMass(const G4double A, const G4double Z)
static constexpr double pi2
G4VLevelDensityParameter * theLevelDensityPtr
G4double GetSpin(void) const
virtual G4double GetCoulombBarrier(G4int ARes, G4int ZRes, G4double U) const =0
G4VCoulombBarrier * theCoulombBarrierPtr
G4double GetPairingCorrection(G4int A, G4int Z) const
void SetCoulomBarrier(const G4VCoulombBarrier *aCoulombBarrierStrategy)
Hep3Vector boostVector() const
virtual G4Fragment * EmittedFragment(G4Fragment *theNucleus)
static constexpr double pi
G4PairingCorrection * GetPairingCorrection()
G4double GetExcitationEnergy() const
virtual G4double GetEmissionProbability(G4Fragment *theNucleus)
static constexpr double GeV
G4double EmissionProbability
HepLorentzVector & boost(double, double, double)