70 delete theFragmentsFactory;
71 delete theFragmentsVector;
76 if (theFragmentsFactory) {
delete theFragmentsFactory; }
78 if (theFragmentsVector) {
88 if (theFragmentsFactory)
delete theFragmentsFactory;
90 if (theFragmentsVector) {
104 if (thePreFragment == 0)
106 G4cout <<
"G4PreCompoundEmission::PerformEmission : "
107 <<
"I couldn't choose a fragment\n"
108 <<
"while trying to de-excite\n"
117 G4double kinEnergyOfEmittedFragment =
124 if(kinEnergyOfEmittedFragment < 0.0) { kinEnergyOfEmittedFragment = 0.0; }
127 AngularDistribution(thePreFragment,aFragment,kinEnergyOfEmittedFragment);
134 EmittedMass + kinEnergyOfEmittedFragment);
146 Rest4Momentum -= Emitted4Momentum;
157 thePreFragment->
GetA());
160 thePreFragment->
GetZ());
176 void G4PreCompoundEmission::AngularDistribution(
195 G4double Eav = 2*p*(p+1)/((p+h)*gg);
201 G4double w_num = rho(p+1, h, gg, Uf, fFermiEnergy);
202 G4double w_den = rho(p, h, gg, Uf, fFermiEnergy);
203 if (w_num > 0.0 && w_den > 0.0)
205 Eav *= (w_num/w_den);
206 Eav += - Uf/(p+h) + fFermiEnergy;
216 G4double Eeff = ekin + Bemission + fFermiEnergy;
227 an = 3*std::sqrt((ProjEnergy+fFermiEnergy)*Eeff)/(zeta*Eav);
230 if ( ne > 1 ) { an /= (
G4double)ne; }
233 if ( an > 10. ) { an = 10.; }
240 if(an < 0.1) { cost = 1. - 2*random; }
243 cost = 1. +
G4Log(1-random*(1-exp2an))/
an;
244 if(cost > 1.) { cost = 1.; }
245 else if(cost < -1.) {cost = -1.; }
254 G4double sint = std::sqrt((1.0-cost)*(1.0+cost));
256 theFinalMomentum.
set(pmag*std::cos(phi)*sint,pmag*std::sin(phi)*sint,
261 theFinalMomentum.
rotateUz(theIncidentDirection);
268 G4double Aph = (p*p + h*h + p - 3.0*h)/(4.0*gg);
271 if ( E - Aph < 0.0) {
return 0.0; }
283 if(logt3 > logmax) { logt3 = logmax; }
289 for(
G4int j=1; j<=h; ++j)
292 if(Eeff < 0.0) {
break; }
295 logt3 = (p+h-1) *
G4Log( Eeff) + logConst;
296 if(logt3 > logmax) { logt3 = logmax; }
297 tot += t1*t2*
G4Exp(logt3);
void set(double x, double y, double z)
static G4Pow * GetInstance()
Hep3Vector boostVector() const
G4ReactionProduct * GetReactionProduct() const
std::vector< G4VPreCompoundFragment * > * GetFragmentVector()
G4int GetNumberOfParticles() const
void SetVector(pcfvector *avector)
G4double GetBindingEnergy() const
G4int GetNumberOfHoles() const
void SetMomentum(const G4LorentzVector &value)
G4GLOB_DLL std::ostream G4cout
const G4LorentzVector & GetMomentum() const
HepLorentzVector & boost(double, double, double)
void SetMomentum(const G4LorentzVector &value)
Hep3Vector & rotateUz(const Hep3Vector &)
G4double GetFermiEnergy()
void SetNumberOfParticles(G4int value)
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)
void SetNumberOfCharged(G4int value)
virtual G4double GetKineticEnergy(const G4Fragment &aFragment)=0
G4int GetNumberOfCharged() const
G4double GetLevelDensity()
G4double GetExcitationEnergy() const