54 const G4double tGamma = 87.65*fermi*fermi*fermi*fermi*fermi*
fermi;
60 G4int nBarion = Barions.size();
72 Etot += t1*std::pow(Alpha/
pi, 3/2)*std::exp(-Alpha*r12*r12);
75 Etot += Vo*0.5/r12*std::exp(1/(4*Alpha*GammaY*GammaY))*
76 (std::exp(-r12/GammaY)*(1 -
Erf(0.5/GammaY/std::sqrt(Alpha) - std::sqrt(Alpha)*r12)) -
77 std::exp( r12/GammaY)*(1 -
Erf(0.5/GammaY/std::sqrt(Alpha) + std::sqrt(Alpha)*r12)));
91 Etot = tGamma*std::pow(4*Alpha*Alpha/3/
pi/
pi, 1.5)*std::exp(-Alpha*(r12*r12 + r13*r13));
112 const G4double P2[5] = { 7.3738883, 6.8650185, 3.0317993, 0.56316962, 4.3187787e-5 };
113 const G4double Q2[5] = { 7.3739609, 15.184908, 12.79553, 5.3542168, 1. };
127 H = X*(P10 + P11*Y)/(Q10+Y);
135 for(c1 = 3; c1 >= 0; c1--)
140 H = 1 - std::exp(-V*V)*AP/AQ;
145 H = 1 - std::exp(-V*V)*(C1+Y*(P30 + P31*Y)/(Q30 + Y))/V;
209 for(
size_t cParticle = 0; cParticle < theBarions.size(); cParticle++)
G4double GetPionPotential(G4double radius)
static c2_factory< G4double > c2
static const G4double a_pion
! for pions it has todiffer from kaons
const G4ThreeVector & GetPosition() const
static const G4double P10[nE]
static const G4double P11[nE]
G4ParticleDefinition * GetDefinition() const
void Transport(G4KineticTrackVector &theActive, const G4KineticTrackVector &theSpectators, G4double theTimeStep)
void SetPosition(const G4ThreeVector aPosition)
G4double CalculateTotalEnergy(const G4KineticTrackVector &Barions)
static const G4double a_antiproton
static const G4double coulomb
G4double GetKaonPotential(G4double radius)
static const G4double * P2[nN]
G4double GetProtonPotential(G4double radius)
G4double GetExcitationEnergy(G4int nHitNucleons, const G4KineticTrackVector &theParticles)
void Integrate(const G4KineticTrackVector &theActive, G4double theTimeStep)
G4double GetNeutronPotential(G4double radius)
G4double GetPDGCharge() const
const G4LorentzVector & Get4Momentum() const
static const double fermi
static const G4double a_kaon
G4double GetAntiprotonPotential(G4double radius)