109 if (thePreFragment ==
nullptr)
111 G4cout <<
"G4PreCompoundEmission::PerformEmission : "
112 <<
"I couldn't choose a fragment\n"
113 <<
"while trying to de-excite\n"
123 kinEnergy =
std::max(kinEnergy, 0.0);
130 std::sqrt(kinEnergy*(kinEnergy + 2.0*thePreFragment->
GetNuclearMass()));
150 Rest4Momentum -= Emitted4Momentum;
161 thePreFragment->
GetA());
164 thePreFragment->
GetZ());
196 G4double Eav = 2*p*(p+1)/((p+h)*gg);
204 if (w_num > 0.0 && w_den > 0.0)
206 Eav *= (w_num/w_den);
231 if ( ne > 1 ) {
an /=
static_cast<G4double>(ne); }
234 if (
an > 10. ) {
an = 10.; }
241 if(
an < 0.1) { cost = 1. - 2*random; }
244 cost = 1. +
G4Log(1-random*(1-exp2an))/
an;
245 if(cost > 1.) { cost = 1.; }
246 else if(cost < -1.) {cost = -1.; }
255 G4double sint = std::sqrt((1.0-cost)*(1.0+cost));
269 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);
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double G4Log(G4double x)
G4ThreeVector G4RandomDirection()
static constexpr double pi2
G4GLOB_DLL std::ostream G4cout
void set(double x, double y, double z)
Hep3Vector & rotateUz(const Hep3Vector &)
Hep3Vector boostVector() const
HepLorentzVector & boost(double, double, double)
G4bool UseAngularGen() const
G4double GetFermiEnergy() const
G4int GetNumberOfParticles() const
G4int GetNumberOfHoles() const
void SetNumberOfCharged(G4int value)
G4double GetExcitationEnergy() const
const G4LorentzVector & GetMomentum() const
void SetCreatorModelID(G4int value)
G4int GetNumberOfExcitons() const
void SetMomentum(const G4LorentzVector &value)
void SetNumberOfParticles(G4int value)
G4int GetNumberOfCharged() const
void SetZandA_asInt(G4int Znew, G4int Anew)
G4double GetLevelDensity(G4int Z, G4int A, G4double U)
G4DeexPrecoParameters * GetParameters()
static G4NuclearLevelData * GetInstance()
static G4int GetModelID(const G4int modelIndex)
static G4Pow * GetInstance()
G4double logfactorial(G4int Z) const
G4NuclearLevelData * fNuclData
G4bool fUseAngularGenerator
G4PreCompoundFragmentVector * theFragmentsVector
G4ReactionProduct * PerformEmission(G4Fragment &aFragment)
G4VPreCompoundEmissionFactory * theFragmentsFactory
G4ThreeVector theFinalMomentum
void AngularDistribution(G4VPreCompoundFragment *theFragment, const G4Fragment &aFragment, G4double KineticEnergy)
G4double rho(G4int p, G4int h, G4double gg, G4double E, G4double Ef) const
G4VPreCompoundFragment * ChooseFragment()
void SetVector(pcfvector *avector)
void SetCreatorModelID(const G4int mod)
std::vector< G4VPreCompoundFragment * > * GetFragmentVector()
virtual G4double SampleKineticEnergy(const G4Fragment &aFragment)=0
G4double GetNuclearMass() const
G4double GetBindingEnergy() const
G4ReactionProduct * GetReactionProduct() const
void SetMomentum(const G4LorentzVector &value)
static constexpr double twopi
static constexpr double MeV
T max(const T t1, const T t2)
brief Return the largest of the two arguments