104 if (thePreFragment == 0)
106 G4cout <<
"G4PreCompoundEmission::PerformEmission : "
107 <<
"I couldn't choose a fragment\n"
108 <<
"while trying to de-excite\n"
123 if(kinEnergy < 0.0) { kinEnergy = 0.0; }
136 Emitted4Momentum.boost(Rest4Momentum.boostVector());
144 Rest4Momentum -= Emitted4Momentum;
155 thePreFragment->
GetA());
158 thePreFragment->
GetZ());
189 G4double Eav = 2*p*(p+1)/((p+h)*gg);
195 G4double w_num =
rho(p+1, h, gg, Uf, fFermiEnergy);
197 if (w_num > 0.0 && w_den > 0.0)
199 Eav *= (w_num/w_den);
200 Eav += - Uf/(p+h) + fFermiEnergy;
221 an = 3*std::sqrt((ProjEnergy+fFermiEnergy)*Eeff)/(zeta*Eav);
224 if ( ne > 1 ) { an /= (
G4double)ne; }
227 if ( an > 10. ) { an = 10.; }
234 if(an < 0.1) { cost = 1. - 2*random; }
237 cost = 1. +
G4Log(1-random*(1-exp2an))/
an;
238 if(cost > 1.) { cost = 1.; }
239 else if(cost < -1.) {cost = -1.; }
248 G4double sint = std::sqrt((1.0-cost)*(1.0+cost));
262 G4double Aph = (p*p + h*h + p - 3.0*h)/(4.0*gg);
265 if ( E - Aph < 0.0) {
return 0.0; }
277 if(logt3 > logmax) { logt3 = logmax; }
283 for(
G4int j=1; j<=h; ++j)
286 if(Eeff < 0.0) {
break; }
289 logt3 = (p+h-1) *
G4Log( Eeff) + logConst;
290 if(logt3 > logmax) { logt3 = logmax; }
291 tot += t1*t2*
G4Exp(logt3);
static G4Pow * GetInstance()
G4ReactionProduct * GetReactionProduct() const
CLHEP::Hep3Vector G4ThreeVector
std::vector< G4VPreCompoundFragment * > * GetFragmentVector()
void AngularDistribution(G4VPreCompoundFragment *theFragment, const G4Fragment &aFragment, G4double KineticEnergy)
G4int GetNumberOfParticles() const
void SetVector(pcfvector *avector)
G4double rho(G4int p, G4int h, G4double gg, G4double E, G4double Ef) const
G4double GetBindingEnergy() const
G4int GetNumberOfHoles() const
void SetMomentum(const G4LorentzVector &value)
G4GLOB_DLL std::ostream G4cout
const G4LorentzVector & GetMomentum() const
void SetMomentum(const G4LorentzVector &value)
G4ThreeVector theFinalMomentum
static const double twopi
G4double GetFermiEnergy()
G4PreCompoundFragmentVector * theFragmentsVector
void SetNumberOfParticles(G4int value)
G4VPreCompoundEmissionFactory * theFragmentsFactory
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4VPreCompoundFragment * ChooseFragment()
G4double GetNuclearMass() const
G4ReactionProduct * PerformEmission(G4Fragment &aFragment)
G4int GetNumberOfExcitons() const
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double logfactorial(G4int Z) const
void SetZandA_asInt(G4int Znew, G4int Anew)
virtual G4double SampleKineticEnergy(const G4Fragment &aFragment)=0
void SetNumberOfCharged(G4int value)
G4int GetNumberOfCharged() const
G4double GetLevelDensity()
G4double GetExcitationEnergy() const
CLHEP::HepLorentzVector G4LorentzVector