114 ((AT==1 && ZT==1) || (AP==1 && ZP==1) ||
115 (AT==1 && ZT==0) || (AP==1 && ZP==0) ||
116 (AT==2 && ZT==1) || (AP==2 && ZP==1) ||
117 (AT==3 && ZT==2) || (AP==3 && ZP==2) ||
118 (AT==4 && ZT==2) || (AP==4 && ZP==2))) { result =
true; }
160 if (E_cm <= 0.0) {
return 0.; }
161 if (E_cm <= (0.8 + 0.04*ZT)*xAP && !
lowEnergyCheck) {
return 0.; }
174 G4double Radius = (r_p + r_t)/
fermi + 1.2*(AT13 + AP13)/E_cm13;
185 if ((AT==1 && ZT==1) || (AP==1 && ZP==1)) {
187 D = 1.85 + 0.16/(1+
G4Exp((500.0-E)/200.0));
189 }
else if ((AT==1 && ZT==0) || (AP==1 && ZP==0)) {
191 D = 1.85 + 0.16/(1+
G4Exp((500.0-E)/200.0));
193 }
else if ((AT==2 && ZT==1) || (AP==2 && ZP==1)) {
195 D = 1.65 + 0.1/(1+
G4Exp((500.0-E)/200.0));
197 }
else if ((AT==3 && ZT==2) || (AP==3 && ZP==2)) {
201 }
else if (AP==4 && ZP==2) {
202 if (AT==4 && ZT==2) {T1 = 40.0; G = 300.0;}
203 else if (ZT==4) {T1 = 25.0; G = 300.0;}
204 else if (ZT==7) {T1 = 40.0; G = 500.0;}
205 else if (ZT==13) {T1 = 25.0; G = 300.0;}
206 else if (ZT==26) {T1 = 40.0; G = 300.0;}
207 else {T1 = 40.0; G = 75.0;}
208 D = 2.77 - 8.0E-3*AT + 1.8E-5*AT*AT-0.8/(1.0+
G4Exp((250.0-E)/G));
210 else if (AT==4 && ZT==2) {
211 if (AP==4 && ZP==2) {T1 = 40.0; G = 300.0;}
212 else if (ZP==4) {T1 = 25.0; G = 300.0;}
213 else if (ZP==7) {T1 = 40.0; G = 500.0;}
214 else if (ZP==13) {T1 = 25.0; G = 300.0;}
215 else if (ZP==26) {T1 = 40.0; G = 300.0;}
216 else {T1 = 40.0; G = 75.0;}
217 D = 2.77 - 8.0E-3*AP + 1.8E-5*AP*AP-0.8/(1.0+
G4Exp((250.0-E)/G));
233 deltaE = 1.85*S + 0.16*S/E_cm13 - C_E + 0.91*(AT-2*ZT)*ZP/(xAT*xAP);
234 X1 = 2.83 - 3.1E-2*AT + 1.7E-4*AT*AT;
238 deltaE = 1.85*S + 0.16*S/E_cm13 - C_E + 0.91*(AP-2*ZP)*ZT/(xAT*xAP);
239 X1 = 2.83 - 3.1E-2*AP + 1.7E-4*AP*AP;
253 if (AT==2 && ZT==1) R_c = 13.5;
254 else if (AT==3 && ZT==2) R_c = 21.0;
255 else if (AT==4 && ZT==2) R_c = 27.0;
256 else if (ZT==3) R_c = 2.2;
258 else if (AT==1 && ZT==1)
260 if (AP==2 && ZP==1) R_c = 13.5;
261 else if (AP==3 && ZP==2) R_c = 21.0;
262 else if (AP==4 && ZP==2) R_c = 27.0;
263 else if (ZP==3) R_c = 2.2;
265 else if (AP==2 && ZP==1)
267 if (AT==2 && ZT==1) R_c = 13.5;
268 else if (AT==4 && ZT==2) R_c = 13.5;
269 else if (AT==12 && ZT==6) R_c = 6.0;
271 else if (AT==2 && ZT==1)
273 if (AP==2 && ZP==1) R_c = 13.5;
274 else if (AP==4 && ZP==2) R_c = 13.5;
275 else if (AP==12 && ZP==6) R_c = 6.0;
277 else if ((AP==4 && ZP==2 && (ZT==73 || ZT==79)) ||
278 (AT==4 && ZT==2 && (ZP==73 || ZP==79))) R_c = 0.6;
287 result =
pi * xr * xr * (1.0 - R_c*B/E_cm) * X_m;
302 if (resultp > result) { result = 0.0; }
static G4Pow * GetInstance()
static G4double GetNuclearMass(const G4double A, const G4double Z)
G4double GetWilsonRMSRadius(G4double A)
G4double GetKineticEnergy() const
G4ParticleDefinition * GetDefinition() const
double B(double temperature)
static G4NistManager * Instance()
virtual G4double GetElementCrossSection(const G4DynamicParticle *theProjectile, G4int Z, const G4Material *mat=0)
virtual G4bool IsElementApplicable(const G4DynamicParticle *theProjectile, G4int Z, const G4Material *)
G4double Z13(G4int Z) const
void SetLowEnergyCheck(G4bool)
G4TripathiLightCrossSection()
void SetKineticEnergy(G4double aEnergy)
G4LorentzVector Get4Momentum() const
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
G4double A13(G4double A) const
~G4TripathiLightCrossSection()
G4double GetAtomicMassAmu(const G4String &symb) const
static const double eplus
G4double GetPDGCharge() const
static const double fermi
G4int GetBaryonNumber() const
G4WilsonRadius * theWilsonRadius
CLHEP::HepLorentzVector G4LorentzVector