108 if (thePreFragment ==
nullptr)
110 G4cout <<
"G4PreCompoundEmission::PerformEmission : "
111 <<
"I couldn't choose a fragment\n"
112 <<
"while trying to de-excite\n"
122 kinEnergy =
std::max(kinEnergy, 0.0);
129 std::sqrt(kinEnergy*(kinEnergy + 2.0*thePreFragment->
GetNuclearMass()));
149 Rest4Momentum -= Emitted4Momentum;
160 thePreFragment->
GetA());
163 thePreFragment->
GetZ());
194 G4double Eav = 2*p*(p+1)/((p+h)*gg);
200 G4double w_num =
rho(p+1, h, gg, Uf, fFermiEnergy);
202 if (w_num > 0.0 && w_den > 0.0)
204 Eav *= (w_num/w_den);
205 Eav += - Uf/(p+h) + fFermiEnergy;
226 an = 3*std::sqrt((ProjEnergy+fFermiEnergy)*Eeff)/(zeta*Eav);
229 if ( ne > 1 ) { an /= (
G4double)ne; }
232 if ( an > 10. ) { an = 10.; }
239 if(an < 0.1) { cost = 1. - 2*random; }
242 cost = 1. +
G4Log(1-random*(1-exp2an))/
an;
243 if(cost > 1.) { cost = 1.; }
244 else if(cost < -1.) {cost = -1.; }
253 G4double sint = std::sqrt((1.0-cost)*(1.0+cost));
267 G4double Aph = (p*p + h*h + p - 3.0*h)/(4.0*gg);
270 if ( E - Aph < 0.0) {
return 0.0; }
282 if(logt3 > logmax) { logt3 = logmax; }
288 for(
G4int j=1; j<=h; ++j)
291 if(Eeff < 0.0) {
break; }
294 logt3 = (p+h-1) *
G4Log( Eeff) + logConst;
295 if(logt3 > logmax) { logt3 = logmax; }
296 tot += t1*t2*
G4Exp(logt3);
G4VPreCompoundEmissionFactory * theFragmentsFactory
void set(double x, double y, double z)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double GetFermiEnergy() const
G4ReactionProduct * GetReactionProduct() const
const G4LorentzVector & GetMomentum() const
G4ReactionProduct * PerformEmission(G4Fragment &aFragment)
void SetMomentum(const G4LorentzVector &value)
static constexpr double MeV
G4DeexPrecoParameters * GetParameters()
void SetZandA_asInt(G4int Znew, G4int Anew)
virtual G4double SampleKineticEnergy(const G4Fragment &aFragment)=0
static G4NuclearLevelData * GetInstance()
G4bool UseAngularGen() const
std::vector< G4VPreCompoundFragment * > * GetFragmentVector()
G4int GetNumberOfHoles() const
Hep3Vector & rotateUz(const Hep3Vector &)
void AngularDistribution(G4VPreCompoundFragment *theFragment, const G4Fragment &aFragment, G4double KineticEnergy)
G4ThreeVector G4RandomDirection()
G4int GetNumberOfParticles() const
G4double G4Log(G4double x)
static G4Pow * GetInstance()
G4PreCompoundFragmentVector * theFragmentsVector
G4int GetNumberOfCharged() const
G4double GetBindingEnergy() const
void SetMomentum(const G4LorentzVector &value)
static constexpr double pi2
void SetNumberOfCharged(G4int value)
G4ThreeVector theFinalMomentum
G4double logfactorial(G4int Z) const
G4bool fUseAngularGenerator
static constexpr double twopi
G4double GetLevelDensity() const
G4double GetNuclearMass() const
void SetVector(pcfvector *avector)
void SetNumberOfParticles(G4int value)
G4GLOB_DLL std::ostream G4cout
Hep3Vector boostVector() const
G4int GetNumberOfExcitons() const
G4double rho(G4int p, G4int h, G4double gg, G4double E, G4double Ef) const
G4double GetExcitationEnergy() const
HepLorentzVector & boost(double, double, double)
G4VPreCompoundFragment * ChooseFragment()