Changeset 1347 for trunk/source/processes/hadronic/cross_sections/src/G4TripathiLightCrossSection.cc

Timestamp:

Dec 22, 2010, 3:52:27 PM (14 years ago)

Author:

garnier

Message:

geant4 tag 9.4

File:

• trunk/source/processes/hadronic/cross_sections/src/G4TripathiLightCrossSection.cc (modified) (3 diffs)

trunk/source/processes/hadronic/cross_sections/src/G4TripathiLightCrossSection.cc

@@ -55 +55 @@
 // 15 March 2004, P R Truscott, QinetiQ Ltd, UK
 // Beta release
+//
+// J. M. Quesada 24 November 2010 bug fixed in X_m 
+//(according to eq. 14 in  R.K. Tripathi et al. Nucl. Instr. and Meth. in Phys. Res. B 155 (1999) 349-356)
 //
 // %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -247 +250 @@
   }
   G4double S_L = 1.2 + 1.6*(1.0-std::exp(-E/15.0));
-  G4double X_m = 1.0 - X1*std::exp(-E/X1*S_L);
+//JMQ 241110 bug fixed 
+//  G4double X_m = 1.0 - X1*std::exp(-E/X1*S_L);
+  G4double X_m = 1.0 - X1*std::exp(-E/(X1*S_L));
 //
 //
@@ -341 +346 @@
     G4int ZZ = G4lrint(theTarget->GetZ());
     G4int AA = G4lrint(theTarget->GetN());
-    xsection = GetIsoZACrossSection(theProjectile, ZZ, AA, theTemperature);
+    xsection = GetZandACrossSection(theProjectile, ZZ, AA, theTemperature);
   }