57 const G4double tGamma = 87.65*fermi*fermi*fermi*fermi*fermi*
fermi;
63 G4int nBarion = Barions.size();
78 Etot += Vo*0.5/r12*
G4Exp(1/(4*Alpha*GammaY*GammaY))*
79 (
G4Exp(-r12/GammaY)*(1 -
Erf(0.5/GammaY/std::sqrt(Alpha) - std::sqrt(Alpha)*r12)) -
80 G4Exp( r12/GammaY)*(1 -
Erf(0.5/GammaY/std::sqrt(Alpha) + std::sqrt(Alpha)*r12)));
88 for(
G4int c3 =
c2 + 1; c3 < nBarion; c3++)
115 const G4double P2[5] = { 7.3738883, 6.8650185, 3.0317993, 0.56316962, 4.3187787e-5 };
116 const G4double Q2[5] = { 7.3739609, 15.184908, 12.79553, 5.3542168, 1. };
130 H = X*(P10 + P11*
Y)/(Q10+Y);
138 for(c1 = 3; c1 >= 0; c1--)
143 H = 1 -
G4Exp(-V*V)*AP/AQ;
148 H = 1 -
G4Exp(-V*V)*(C1+Y*(P30 + P31*
Y)/(Q30 + Y))/V;
212 for(
size_t cParticle = 0; cParticle < theBarions.size(); cParticle++)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double GetExcitationEnergy(G4int nHitNucleons, const G4KineticTrackVector &theParticles)
static const G4double a_pion
! for pions it has todiffer from kaons
const G4LorentzVector & Get4Momentum() const
G4double GetKaonPotential(G4double radius)
static constexpr double MeV
G4double CalculateTotalEnergy(const G4KineticTrackVector &Barions)
G4double A23(G4double A) const
static const G4double * P2[nN]
const G4ParticleDefinition * GetDefinition() const
G4double GetPDGCharge() const
G4double G4Log(G4double x)
static constexpr double fermi
static G4Pow * GetInstance()
G4double powA(G4double A, G4double y) const
void Transport(G4KineticTrackVector &theActive, const G4KineticTrackVector &theSpectators, G4double theTimeStep)
static const G4double P11[nE]
static const G4double P10[nE]
G4double GetProtonPotential(G4double radius)
void Integrate(const G4KineticTrackVector &theActive, G4double theTimeStep)
static const G4double coulomb
void SetPosition(const G4ThreeVector aPosition)
Hep3Vector boostVector() const
static const G4double a_antiproton
static constexpr double pi
const G4ThreeVector & GetPosition() const
G4double GetPionPotential(G4double radius)
G4double GetAntiprotonPotential(G4double radius)
typedef void(XMLCALL *XML_ElementDeclHandler)(void *userData
G4double GetNeutronPotential(G4double radius)
static const G4double a_kaon