109 {
if (rth1 > 0.0 && rth1 <= 1.0) rth = rth1;}
128 if (rT-rP<=r && r<=rT+rP) valueP = 0.125*R*U*S - 0.125*(0.5*R*U+1.0)*T;
133 if (rP-rT<=r && r<=rP+rT) valueP = 0.125*R*U*S - 0.125*(0.5*std::sqrt(n/m)*U-
134 (std::sqrt(1.0-m*m)/n - 1.0)*std::sqrt((2.0-m)/std::pow(m,5.0)))*T;
135 else valueP = (std::sqrt(1.0-m*m)/n-1.0)*std::sqrt(1.0-b*b/n/n);
138 if (!(valueP <= 1.0 && valueP>= -1.0))
140 if (valueP > 1.0) valueP = 1.0;
159 if (rT-rP<=r && r<=rT+rP) valueF = 0.75*R*S - 0.125*(3.0*R-1.0)*T;
164 if (rP-rT<=r && r<=rP+rT) valueF = 0.75*R*S - 0.125*(3.0*std::sqrt(n/m)-
165 (1.0-std::pow(1.0-m*m,3.0/2.0))*std::sqrt(1.0-std::pow(1.0-m,2.0))/std::pow(m,3.0))*T;
166 else valueF = (1.0-std::pow(1.0-m*m,3.0/2.0))*std::sqrt(1.0-b*b/n/n);
169 if (!(valueF <= 1.0 && valueF>= 0.0))
171 if (valueF > 1.0) valueF = 1.0;
185 (1.0+P1-std::pow(1.0-F1,2.0/3.0));
190 if (AP < 12.0) omega = 1500.0;
191 else if (AP <= 16.0) omega = 1500.0 - 320.0*(AP-12.0);
192 Es *= 1.0 + F1*(5.0+omega*F1*F1);
197 else if (Es > B * AP)
217 (1.0+P1-std::pow(1.0-F1,2.0/3.0));
220 if ((r-rT)/rP < rth) {
222 if (AT < 12.0) omega = 1500.0;
223 else if (AT <= 16.0) omega = 1500.0 - 320.0*(AT-12.0);
224 Es *= 1.0 + F1*(5.0+omega*F1*F1);
229 else if (Es > B * AT)
232 delete revAbrasionGeometry;
G4double GetExcitationEnergyOfTarget()
G4NuclearAbrasionGeometry(G4double AP, G4double AT, G4double r)
~G4NuclearAbrasionGeometry()
G4double GetPeripheralThreshold()
G4double GetExcitationEnergyOfProjectile()
void SetPeripheralThreshold(G4double)
G4double GetWilsonRadius(G4double A)