- Timestamp:
- Jun 18, 2010, 11:42:07 AM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/electromagnetic/utils/include/G4EmCorrections.hh
r1228 r1315 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4EmCorrections.hh,v 1.2 4 2008/09/12 14:44:48vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-0 3$26 // $Id: G4EmCorrections.hh,v 1.25 2010/06/04 09:28:46 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-cand-01 $ 28 28 // 29 29 // ------------------------------------------------------------------- … … 56 56 57 57 #include "globals.hh" 58 #include "G4AtomicShells.hh"59 58 #include "G4ionEffectiveCharge.hh" 60 59 #include "G4Material.hh" … … 263 262 G4int nbinCorr; 264 263 265 G4AtomicShells shells;266 264 G4ionEffectiveCharge effCharge; 267 265 … … 356 354 tmax = 2.0*electron_mass_c2*bg2 /(1. + 2.0*gamma*ratio + ratio*ratio); 357 355 charge = p->GetPDGCharge()/eplus; 358 if(charge < 1.5) {q2 = charge*charge;} 359 else { 360 q2 = effCharge.EffectiveChargeSquareRatio(p,mat,kinEnergy); 361 charge = std::sqrt(q2); 362 } 356 //if(charge < 1.5) {q2 = charge*charge;} 357 //else { 358 // q2 = effCharge.EffectiveChargeSquareRatio(p,mat,kinEnergy); 359 // charge = std::sqrt(q2); 360 //} 361 if(charge > 1.5) { charge = effCharge.EffectiveCharge(p,mat,kinEnergy); } 362 q2 = charge*charge; 363 363 } 364 364 if(mat != material) {
Note: See TracChangeset
for help on using the changeset viewer.