59 std::vector< G4QMDParticipant* >::iterator it;
61 p += (*it)->Get4Momentum();
72 std::vector< G4QMDParticipant* >::iterator it;
92 std::vector< G4QMDParticipant* >::iterator it;
125 - Asym * ( N -
Z )* ( N - Z ) /
A;
151 for (
G4int i= 0; i <
n ; i++ )
155 G4double trans = gamma / ( gamma + 1.0 ) * p_i * beta;
161 pcm0 = pcm0 / double ( n );
165 for (
G4int i= 0; i <
n ; i++ )
177 for (
G4int i= 0; i <
n ; i++ )
180 G4double trans = gamma / ( gamma + 1.0 ) * ri * beta;
186 rcm0 +=
rcm[i]*
es[i];
193 for (
G4int i= 0; i <
n ; i++ )
202 for (
G4int i= 0; i <
n ; i++ )
204 rl +=
rcm[i].cross (
pcm[i] );
207 jj =
int ( std::sqrt ( rl*rl /
hbc ) + 0.5 );
213 for (
G4int i= 0; i <
n ; i++ )
G4double GetNuclearMass()
G4double A13(G4double A) const
static constexpr double MeV
void CalEnergyAndAngularMomentumInCM()
G4int GetTotalNumberOfParticipant()
G4double A23(G4double A) const
static G4Proton * Proton()
G4double GetPDGMass() const
static G4double GetBindingEnergy(const G4int A, const G4int Z)
static G4Pow * GetInstance()
std::vector< G4ThreeVector > pcm
G4double excitationEnergy
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
double A(double temperature)
static G4double GetNuclearMass(const G4double A, const G4double Z)
G4LorentzVector Get4Momentum()
G4ThreeVector GetMomentum()
static G4Neutron * Neutron()
std::vector< G4ThreeVector > rcm
std::vector< G4QMDParticipant * > participants
G4QMDParticipant * GetParticipant(G4int i)
static G4QMDParameters * GetInstance()
G4ThreeVector GetPosition()
G4double bindingEnergy(G4int A, G4int Z)
static constexpr double GeV
std::vector< G4double > es