41 : particle(part), theCoulombBarrierPtr(aCoulombBarrier),
42 theRestNucleusA(0),theRestNucleusZ(0),theBindingEnergy(0.0),
43 theMaximalKineticEnergy(-
MeV),theRestNucleusMass(0.0),
44 theReducedMass(0.0),theMomentum(0.,0.,0.,0.),
45 theEmissionProbability(0.0),theCoulombBarrier(0.0),
46 OPTxs(3),useSICB(false)
53 theRestNucleusA13 = 0;
72 <<
"PreCompoundModel Emitted Fragment: Z= " << theFragment->
GetZ()
73 <<
" A= " << theFragment->
GetA()
81 theRestNucleusA = aFragment.
GetA_asInt() - theA;
82 theRestNucleusZ = aFragment.
GetZ_asInt() - theZ;
84 if ((theRestNucleusA < theRestNucleusZ) ||
85 (theRestNucleusA < theA) ||
86 (theRestNucleusZ < theZ))
89 theMaximalKineticEnergy = 0.0;
93 theRestNucleusA13 =
g4pow->
Z13(theRestNucleusA);
103 theReducedMass = theRestNucleusMass*theMass/(theRestNucleusMass + theMass);
113 theMaximalKineticEnergy =
114 ((Ecm-theRestNucleusMass)*(Ecm+theRestNucleusMass) + theMass*theMass)/(2.0*Ecm)-theMass;
static G4Pow * GetInstance()
static G4double GetNuclearMass(const G4double A, const G4double Z)
G4VPreCompoundFragment(const G4ParticleDefinition *, G4VCoulombBarrier *aCoulombBarrier)
void Initialize(const G4Fragment &aFragment)
G4double Z13(G4int Z) const
const G4LorentzVector & GetMomentum() const
G4double GetGroundStateMass() const
G4double GetCoulombBarrier() const
virtual ~G4VPreCompoundFragment()
G4double GetNuclearMass() const
G4double theCoulombBarrier
G4double GetPDGMass() const
std::ostream & operator<<(std::ostream &, const BasicVector3D< float > &)
G4PreCompoundParameters * theParameters
G4double GetPDGCharge() const
G4double GetExcitationEnergy() const
G4int GetBaryonNumber() const