113 {
if (rth1 > 0.0 && rth1 <= 1.0)
rth = rth1;}
137 if (
rP-
rT<=
r &&
r<=
rP+
rT) valueP = 0.125*
R*
U*
S - 0.125*(0.5*std::sqrt(
n/
m)*
U-
139 else valueP = (std::sqrt(1.0-m*m)/
n-1.0)*std::sqrt(1.0-
b*
b/
n/
n);
142 if (!(valueP <= 1.0 && valueP>= -1.0))
144 if (valueP > 1.0) valueP = 1.0;
163 if (
rT-
rP<=
r &&
r<=
rT+
rP) valueF = 0.75*
R*
S - 0.125*(3.0*
R-1.0)*
T;
168 if (
rP-
rT<=
r &&
r<=
rP+
rT) valueF = 0.75*
R*
S - 0.125*(3.0*std::sqrt(
n/
m)-
173 if (!(valueF <= 1.0 && valueF>= 0.0))
175 if (valueF > 1.0) valueF = 1.0;
194 if (
AP < 12.0) omega = 1500.0;
195 else if (
AP <= 16.0) omega = 1500.0 - 320.0*(
AP-12.0);
196 Es *= 1.0 + F1*(5.0+omega*F1*F1);
201 else if (Es >
B *
AP)
226 if (
AT < 12.0) omega = 1500.0;
227 else if (
AT <= 16.0) omega = 1500.0 - 320.0*(
AT-12.0);
228 Es *= 1.0 + F1*(5.0+omega*F1*F1);
233 else if (Es >
B *
AT)
236 delete revAbrasionGeometry;
G4double GetExcitationEnergyOfTarget()
static G4Pow * GetInstance()
G4double powA(G4double A, G4double y) const
G4NuclearAbrasionGeometry(G4double AP, G4double AT, G4double r)
G4double powN(G4double x, G4int n) const
static const G4double * P1[nN]
~G4NuclearAbrasionGeometry()
G4double GetPeripheralThreshold()
G4double A23(G4double A) const
G4double GetExcitationEnergyOfProjectile()
void SetPeripheralThreshold(G4double)
static const double fermi
G4double GetWilsonRadius(G4double A)