74 for (
int i = 0 ; i <
a ; i++ )
104 if ( z == 1 && a == 1 ) {
109 else if ( z == 0 && a == 1 ) {
143 ebin0 = (
ebini - 0.5 ) * 0.001;
144 ebin1 = (
ebini + 0.5 ) * 0.001;
148 ebin0 = (
ebini - 1.5 ) * 0.001;
149 ebin1 = (
ebini + 1.5 ) * 0.001;
164 G4bool areThesePsOK =
false;
191 if ( areThesePsOK ==
false )
continue;
226 rho_l[i] =
cdp * ( rho_a[i] + 1.0 );
228 +
c3p * std::pow ( rho_a[i] ,
gamm )
244 G4bool isThis1stMOK =
false;
256 if ( isThis1stMOK ==
false )
continue;
260 G4bool areTheseMsOK =
true;
272 areTheseMsOK =
false;
283 if ( areTheseMsOK ==
false )
continue;
311 if ( ebinal < ebin0 || ebinal > ebin1 )
345 G4bool isThisEAOK =
false;
353 if ( std::abs ( edif ) <
epse )
377 if ( jfrc == ifrc && std::abs( edif0 ) < std::abs( edif ) )
410 ebinal = ( totalPotentialE + ekinal ) /
double (
GetMassNumber() );
414 if ( isThisEAOK ==
false )
continue;
422 if ( isThisOK ==
false )
424 G4cout <<
"GroundStateNucleus state cannot be created. Try again with another parameters." <<
G4endl;
460 G4cout <<
"GroundStateNucleus state cannot be created. Try again with another parameters." <<
G4endl;
492 rho_l[i] =
cdp * ( rho_a[i] + 1.0 );
494 +
c3p * std::pow ( rho_a[i] ,
gamm )
528 if ( isThisOK ==
true )
break;
536 G4cout <<
"GroundStateNucleus state cannot be created. Try again with another parameters." <<
G4endl;
557 if ( ebinal < ebin0 || ebinal > ebin1 )
590 if ( std::abs ( edif ) <
epse )
612 if ( jfrc == ifrc && std::abs( edif0 ) < std::abs( edif ) )
645 ebinal = ( totalPotentialE + ekinal ) /
double (
GetMassNumber() );
687 rsqr = rx*rx + ry*ry + rz*rz;
689 rrr =
radm * std::sqrt ( rsqr );
690 rwod = 1.0 / ( 1.0 + std::exp ( ( rrr -
rt00 ) /
saa ) );
706 for (
G4int j = 0 ; j < i ; j++ )
730 if ( isThisOK ==
true )
757 pfm = pfm * ( 1.0 + 0.2 * std::sqrt( std::abs( 8.0 +
ebini ) / 8.0 ) );
762 std::vector< G4double > phase;
791 psqr = px*px + py*py + pz*pz;
803 if ( tkdb_i >
maxTrial )
return result;
824 for (
G4int j = 0 ; j < i ; j++ )
840 dist2_p = dist2_p*
cph;
841 expa = expa - dist2_p;
846 phase[j] = std::exp ( expa );
848 if ( phase[j] *
cpc > 0.2 )
870 phase[i] += phase[j];
871 if ( phase[i] *
cpc > 0.3 )
904 if ( isThisOK ==
true )
909 for (
int j = 0 ; j < i ; j++ )
948 rcm_tmp = rcm_tmp/sumMass;
966 for (
G4int i = 0 ; i < 3 ; i++ )
968 for (
G4int j = 0 ; j < 3 ; j++ )
986 rr[0][0] += r.y() * r.y() + r.z() * r.z();
987 rr[1][0] += - r.y() * r.x();
988 rr[2][0] += - r.z() * r.x();
989 rr[0][1] += - r.x() * r.y();
990 rr[1][1] += r.z() * r.z() + r.x() * r.x();
991 rr[2][1] += - r.z() * r.y();
992 rr[2][0] += - r.x() * r.z();
993 rr[2][1] += - r.y() * r.z();
994 rr[2][2] += r.x() * r.x() + r.y() * r.y();
997 for (
G4int i = 0 ; i < 3 ; i++ )
1000 for (
G4int j = 0 ; j < 3 ; j++ )
1002 rr[i][j] = rr[i][j] / x;
1003 ss[i][j] = ss[i][j] / x;
1005 for (
G4int j = 0 ; j < 3 ; j++ )
1010 for (
G4int k = 0 ; k < 3 ; k++ )
1012 rr[j][k] += -y * rr[i][k];
1013 ss[j][k] += -y * ss[i][k];
1026 for (
G4int i = 0 ; i < 3 ; i++ )
1030 for (
G4int j = 0 ; j < 3 ; j++ )
1032 opl[i] += ss[i][j]*rll[j];
1042 p_i += ri.cross(opl_v);
G4double GetRHA(G4int i, G4int j)
std::vector< G4double > rho_l
void killCMMotionAndAngularM()
void SetParticipant(G4QMDParticipant *particle)
CLHEP::Hep3Vector G4ThreeVector
G4ThreeVector GetFFr(G4int i)
void Cal2BodyQuantities()
G4GLOB_DLL std::ostream G4cout
static G4QMDParameters * GetInstance()
G4bool samplingMomentum(G4int)
static G4Proton * Proton()
static G4Neutron * Neutron()
G4double GetTotalPotential()
std::vector< G4QMDParticipant * > participants
G4double GetRHE(G4int i, G4int j)
G4bool samplingPosition(G4int)
static G4double GetBindingEnergy(const G4int A, const G4int Z)
G4double GetRR2(G4int i, G4int j)
void SetSystem(G4QMDSystem *aSystem)
G4ThreeVector GetFFp(G4int i)
std::vector< G4double > d_pot
G4QMDMeanField * meanfield
G4QMDGroundStateNucleus()
std::vector< G4double > phase_g
CLHEP::HepLorentzVector G4LorentzVector